[Talk-de] Mit Google Layers Plugin OSM-Rasterbilder in QuantumGIS laden

[Andre Joost]

M∡rtin Koppenhoefer schrieb:
> Am 24. November 2010 08:45 schrieb André Joost <andre+joost at nurfuerspam.de>:
>> 3785:
>> +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137
>> +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
> 
> 
> das war mal ein Entwurf für eine Spec der EPSG die aber nie offiziell wurde.
> 
> 
>> 3857:
>> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0
>> +k=1.0 +units=m +nadgrids=@null +wktext +no_defs
> 
> 
> die 3857 ist gem. EPSG "valid", allerdings nicht mit den Einstellungen
> (die EPSG hat einen Ellipsoid spezifiziert).
> 
> 
>> 900913:
>> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0
>> +k=1.0 +units=m +nadgrids=@null +wktext +over +no_defs
> 
> 
> die gibts "offiziell" (also bei der EPSG) nicht, ist aber identisch
> mit der 3857 oben (die wie erwähnt nicht die der EPSG ist). Das +over
> ist lediglich ein switch fürs wrapping.
> 

Ich habs mal mit cs2cs aus den FWtools nachrechnen lassen:
Ausgangswert EPSG:4326 7 Grad Ost und 51 Grad Nord
x-Koordinate ist immer 779236
y-Koordiante:
EPSG 3395: 6588067
EPSG 3785: 6621294
EPSG 3857: 6621294
EPSG 900913: 6621294

das gleiche mit obigen proj-Definitionen aus QuantumGis Tethys:
EPSG 3395: 6588067
EPSG 3785: 6588044
EPSG 3857: 6621294
EPSG 900913: 6621294

Da hat sich also bei der Definition von EPSG 3785 ein Fehler
eingeschlichen, der einen Versatz von 33km in unseren Breiten erzeugt,
und zufälligerweise (?) fast der Projektion 3395 entspricht. Ausserdem
ist die Höhenkoordinate -12872, während sie bei allen anderen 0 ist.

-- 
Gruß,
André Joost