[Fossgis-talk] BeTA 2007 in GRASS

Markus Neteler neteler at osgeo.org
Mo Okt 18 23:06:28 CEST 2010


Hallo Tom,

viele Experten sind hier unterwegs:

http://trac.osgeo.org/proj/wiki/WikiStart#MailingList

Am besten waere die Frage wohl dorthin zu senden.

Schoene Gruesse
Markus

2010/10/13 Tom Gottfried <tom.gottfried at wzw.tum.de>:
> 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
> _______________________________________________
> Fossgis-talk-liste mailing list
> Fossgis-talk-liste at fossgis.de
> https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste
>