[Fossgis-talk] BeTA 2007 in GRASS

[Tom Gottfried]

Hallo Listler,

folgendes Problem:
ich habe 2 Punkte in lonlat (ersterer Punkt 2 aus dem Testdatensatz zu BeTA 2007, siehe 
http://crs.bkg.bund.de/crseu/crs/descrtrans/BeTA/de_dhdn2etrs_beta.php, der zweite ein 
Lagefestpunkt, den ich mit RTK-Float eingemessen habe). Wenn ich diese mit cs2cs unter der 
Verwendung von BeTA 2007 in Gauß-Krüger projiziere sieht das folgendermaßen aus:

echo "10.3321172833 48.8323271886
11.8576755 48.3788954" | cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +to +proj=tmerc 
+lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel +units=m +no_defs 
+nadgrids='/home/gottfried/Documents/geobasisdaten/ntv2/BETA2007.gsb' -f "%.3f"
4377657.795     5411879.840 0.000
4489564.124     5360125.555 0.000

Die "amtlichen" Gauß-Krüger-Koordinaten sind folgende:
4377657.795	5411879.840
4489564.391	5360125.452
Die Projektion mit cs2cs passt also perfekt.

Jetzt zu GRASS:
Ich habe die gleichen Punkte in einer WGS84-Location:

g.proj -p
-PROJ_INFO-------------------------------------------------
name       : Lat/Lon
proj       : ll
datum      : wgs84
ellps      : wgs84
no_defs    : defined
-PROJ_UNITS------------------------------------------------
unit       : degree
units      : degrees
meters     : 1.0


Meiner Gauß-Krüger-Location möchte ich BeTA 2007 beibringen:

g.proj -c proj4="+proj=tmerc +lat_0=0 +lon_0=12 +k=1 +x_0=4500000 +y_0=0 +ellps=bessel +units=m 
+no_defs +nadgrids='/home/gottfried/Documents/geobasisdaten/ntv2/BETA2007.gsb'"
WARNUNG: Datum <unknown> nicht erkannt von GRASS und keine Parameter
          gefunden.
Projektionsinformationen aktualisiert!

Die Warnung ist schonmal unsympathisch und die Projektion der Punkte aus der WGS84-Location leider auch:

v.proj testdaten_proj loc=test_world maps=test
v.out.ascii testdaten_proj fs=" "
4377566.08131916 5411769.92497926 1
4489458.2716424 5360021.17394219 2

Ich habe auch schon BETA2007.gsb in $GISBASE/etc/nad/ kopiert und dann "+nadgrids=BETA2007.gsb" und 
"nadgrids=BETA2007" in g.proj verwendet. Das Ergebnis ist beide Male das gleiche.

Weiß jemand, was ich hier falsch mache?

Vielen Dank und viele Grüße,
Tom

-- 
Technische Universität München
Department für Pflanzenwissenschaften
Lehrstuhl für Grünlandlehre
Alte Akademie 12
85350 Freising / Germany
Phone: ++49 (0)8161 715324
Fax:   ++49 (0)8161 713243
email: tom.gottfried at wzw.tum.de
http://www.wzw.tum.de/gruenland