[Mapserver-DE] Problem mit Projektionen (?)

[Dirk Jesko]

Hallo Guido,

Zunächst gehe ich davon aus, dass im PostgreSQL mit 
"transform(the_geom,'4326')" und nicht mit 4226 transformiert wurde.
Ansonsten habe ich festgestellt, dass die Transformationsparameter im 
proj.4 für WGS84->Bessel große Abweichungen erzeugen (zumindest in 
Deutschland).

Ich habe in einem Projekt Daten in 
Gauß-Krüger-Projektion/Bessel-Ellipsoid und andere in 
latlon-Koordinaten/WGS84. Mit den einfachen Einstellungen für die 
Projektion bekomme ich ebenfalls Lagefehler über 100m. Gebe ich jedoch 
die 7 Parameter für den Datumsübergang zwischen Bessel und WGS84 
explizit an, ist die Lage perfekt. Dazu verwende ich folgende Definition 
für die Projektion der WGS84-Layer:

   PROJECTION
     "init=epsg:4326"
     "towgs84=-598.1,-73.7,-418.2,-0.202,-0.045,2.455,-6.70"
   END

Die Parameter stammen von der Webseite des BKG 
(http://crs.bkg.bund.de/crseu/crs/eu-countrysel.php?country=DE). Ich 
habe die Einstellungen in verschiedenen Gegenden Deutschlands getestet 
und konnte bisher keine nennenswerten Abweichung feststellen.


Mit freundlichen Grüßen,
Dirk

P.S.: Falls Interesse besteht, kann ich eine GPS-Aufzeichnung als 
ESRI-Shape zur Verfügung stellen.


Geografische Informationssysteme wrote:

>Hallo UMN.Benutzer,
>
>Ich hab da ein Problem und/oder Unverständnis mit der Darstellung von Punkten und brauche eure Hilfe. 
>
>Also, zukünftig möchte ich Punkte (GPS-Koordinaten) in einer Karte darstellen. 
>Da mir derzeit noch keine Test-Koordinaten vorliegen, habe ich einfach eine bestehende Postgres-Tabelle kopiert und die Inhalte der geom Spalte mit "transform(the_geom,'4226')" umgewandelt.
>
>Dann hab ich nen Layer erstellt, der diese Tabelle anspricht und die Projection-Sektion folgendermaßen gesetzt:
>1. Variante: 
>PROJECTION
>	"init=EPSG:4326"
>END
>2. Variante:
>PROJECTION
>                  "proj=latlong"
>                  "datum=WGS84"
>                  "ellps=WGS84"
>END
>
>Der "globale" PROJECTION Eintrag im Mapfile sieht so aus:
>PROJECTION
>	"init=EPSG:31466" #PROJECTION # DHDN / Germany zone 2
>END
>
>Bei der Darstellung der Punkte (in beiden Fällen) muss ich dann feststellen dass die geclonten Punkte nach unten links rutschen (Entfernung immer konstant ca 206 Meter). Nach meinem Verständnis müssten doch die Punkte deckungsgleich sein, isn't it?
>
>Hier mal die Ergebnisse aus den select statements:
>select the_geom, transform(the_geom,'4326') from pipunkte ergibt:
>"SRID=31466;POINT(2657882 5538488)"	"SRID=4326;POINT(8.19989675269949 49.9618389782115)"
>
>die Gegenprüfung sieht dann so aus:
>select the_geom, transform(the_geom,'31466') from pipunktegps ergibt:
>"SRID=4326;POINT(8.19989675269949 49.9618389782115)"	"SRID=31466;POINT(2657882.00000032 5538487.99999881)"
>
>Also ergibt hin und zurück den gleichen Punkt. Daran kannst also nicht liegen. Es muss m.was mit der Mapserver-Projektion zu tun haben? Oder? 
>
>MfG
>
>Guido
>
>
>_______________________________________________
>Mapserver-DE mailing list
>Mapserver-DE at freegis.org
>https://freegis.org/mailman/listinfo/mapserver-de
>
>  
>