[Grass-de] gps daten in pd83

[Markus Neteler]

2009/8/18 Markus Neteler <neteler at osgeo.org>:
> 2009/8/18 Torsten Eckart <torecki at arcor.de>:
>> Hallo,
>> ich habe einige Daten aus Thüringen im amtlichen KBS PD83 und möchte einige
>> Strecken und Punkte mit gps überprüfen. Leider liegen meine gps-Strecken immer
>> um  ca 150m südlich und 60m westlich von dem wo sie eigentlich liegen sollten.
>> Der Fehler rührt wahrscheinlich von einem falsch eingestellten Datum oder der
>> Transformation aus dem WGS84 in PD83. In arcgis löse ich das Problem indem ich
>> die Transformation von wgs84 auf pd83 einstelle wie mache ich das in grass?
>
> Im Prinzip genauso...
> http://spatialreference.org/ref/epsg/3396/
>
> Du muesstest eine "location" mit EPSG 3396 erzeugen, dort
> sollte dann das Datum erfragt werden.
>
> Ahhh -  ich sehe (hab's gerade mal durchgefuehrt):
>
> g.proj -w
> PROJCS["Transverse Mercator",
>    GEOGCS["bessel",
>        DATUM["unknown",       <<<<<-------!
>            SPHEROID["Bessel_1841",6377397.155,299.1528128]],
>        PRIMEM["Greenwich",0],
>        UNIT["degree",0.0174532925199433]],
>    PROJECTION["Transverse_Mercator"],
>    PARAMETER["latitude_of_origin",0],
>    PARAMETER["central_meridian",9],
>    PARAMETER["scale_factor",1],
>    PARAMETER["false_easting",3500000],
>    PARAMETER["false_northing",0],
>    UNIT["Meter",1]]
>
>
> Nicht gut.
>
> http://www.epsg-registry.org/  berichtet:
> "PD/83 is the realisation of DHDN in Thuringen. It is the resultant of
> applying a transformation derived at 13 points on the border between
> East and West Germany to Pulkovo 1942/83 points in Thuringen."
>
> Aber die Parameter sehe ich nicht...

Ich hab sie wohl gefunden:

http://www.crs-geo.eu
-> CRS-Description
  -> national CRS
    -> Germany
      -> DE_PD/83 / GK_3
        -> Descr. of Transf.
            (http://crs.bkg.bund.de/crseu/crs/eu-description.php?crs_id=dERFX1BELzgzIC8gR0tfMw==)
               -> 7 parameters

Wenn den Code hat, kann man auch hierhin gehen:

http://www.epsg-registry.org/
-> retrieve by code
   ->   15867
      -> Coordinate Transformation [PD/83 to ETRS89 (1)]

Wenn man dann "grep" in der Datumstabelle macht von GRASS,
sieht man, dass es haha schon da war:
datumtransform.table:potsdam
"towgs84=599.4,72.4,419.2,-0.063,-0.022,-2.723,6.46" "Germany
(Thüringen)" "Accuracy <1m"

:-)

Also zurueck auf Start!
Kannst Du mal folgenden Befehl in der GRASS -Location machen?

g.proj -w

Was kommt dabei raus?

Gruesse
Markus