Jag har en del rasterkartor i RT90-koordinater (scannade papperskartor och kartor från Lantmäteriets onlinebutik och Digitala kartbiblioteket), och jag använder QGIS för att plotta GPS-data med rasterkartorna som bakgrund. Eftersom GPS-datan är i grader i WGS84 så måste den konverteras till RT90 för att passa med rasterkartorna. QGIS kan göra det automatiskt och har ett koordinatsystem färdigdefinierat för RT90 2.5 gon V (som kommer från PostGIS' spatial_ref_sys-tabell) med följande definition (i proj4-format):
+proj=tmerc +lat_0=0 +lon_0=15.80827777777778 +k=1.000000 +x_0=1500000 +y_0=0 +ellps=bessel +units=m +no_defs
+proj är projektionen som används (Transverse Mercator), +ellps är ellipsoiden, +lon_0 är centralmeridianen för TM-projektionen, +k är skalfaktor. Dessa parametrar verkar stämma med de som Lantmäteriet har på sin websida om RT90-systemet - men när jag konverterar från WGS84 till ovanstående koordinatsystem så hamnar de konverterade punkterna alltid ca 200 meter sydväst om den korrekta positionen (i alla fall i Stockholmstrakten).
Om jag däremot använder den här definitionen, som jag hittade någonstans där det stod att det var en approximation av RT90 med ett fel på någon eller några meter, hamnar punkterna precis där de ska:
+ellps=GRS80 +proj=tmerc +lon_0=15.80628 +x_0=1500064.274 +y_0=-667.711 +k=1.00000561024
Är det någon kunnig som ser något uppenbart fel i den första definitionen? Eller har jag hittat en bugg i proj4-biblioteket (kanske något litet numeriskt fel vid konvertering från WGS84 till Bessel) ?
+proj=tmerc +lat_0=0 +lon_0=15.80827777777778 +k=1.000000 +x_0=1500000 +y_0=0 +ellps=bessel +units=m +no_defs
+proj är projektionen som används (Transverse Mercator), +ellps är ellipsoiden, +lon_0 är centralmeridianen för TM-projektionen, +k är skalfaktor. Dessa parametrar verkar stämma med de som Lantmäteriet har på sin websida om RT90-systemet - men när jag konverterar från WGS84 till ovanstående koordinatsystem så hamnar de konverterade punkterna alltid ca 200 meter sydväst om den korrekta positionen (i alla fall i Stockholmstrakten).
Om jag däremot använder den här definitionen, som jag hittade någonstans där det stod att det var en approximation av RT90 med ett fel på någon eller några meter, hamnar punkterna precis där de ska:
+ellps=GRS80 +proj=tmerc +lon_0=15.80628 +x_0=1500064.274 +y_0=-667.711 +k=1.00000561024
Är det någon kunnig som ser något uppenbart fel i den första definitionen? Eller har jag hittat en bugg i proj4-biblioteket (kanske något litet numeriskt fel vid konvertering från WGS84 till Bessel) ?
