[Fossgis-talk] BeTA 2007 in GRASS

Tom Gottfried tom.gottfried at wzw.tum.de
Di Okt 19 15:02:36 CEST 2010 

Hallo Markus,

danke für den Hinweis! Während ich einigen Hinweisen von der PROJ.4-Liste nachgegangen bin, bin ich 
eher zufällig auf g.setproj gestoßen und habe festgestellt, dass es damit geht. Und zwar folgendermaßen:

_ die Datei BETA2007.gsb nach $GISBASE/etc/nad/ kopieren
_ in $GISBASE/etc/datumtransform.table die folgende Zeile einfügen:
potsdam "nadgrids=BETA2007.gsb" "Germany (Whole Country)" "Accuracy <1m"

_ g.setproj ausführen und die entsprechenden Parameter auswählen

Anschließend erhalte ich mit v.proj die richtigen Koordinaten für die Testpunkte.
Mit g.proj geht es auch nach dem kopieren der Datei und der Änderung im datumtransform.table nicht.

Der kleine Haken an der Sache:
mit `g.proj -jf' erhalte ich jetzt folgendes:
"+proj=tmerc +k_0=1.0000000000 +lat_0=0.0000000000 +lon_0=12.0000000000 +x_0=4500000.0000000000 
+y_0=0.0000000000 +a=6377397.155 +rf=299.1528128 +no_defs +nadgrids=/opt/grass/etc/nad/BETA2007.gsb 
+to_meter=1.0"

Dieser PROJ.4-string scheint ungültig zu sein:

`g.proj -c proj4="proj=tmerc +k_0=1.0000000000 +lat_0=0.0000000000 +lon_0=12.0000000000 
+x_0=4500000.0000000000 +y_0=0.0000000000 +a=6377397.155 +rf=299.1528128 +no_defs 
+nadgrids=/opt/grass/etc/nad/BETA2007.gsb +to_meter=1.0"'

gibt das zurück
"FEHLER: Kann den PROJ.4 Parameterstring nicht verstehen.",

und

`echo "10.3321172833 48.8323271886
11.8576755 48.3788954" | cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +to proj=tmerc 
+k_0=1.0000000000 +lat_0=0.0000000000 +lon_0=12.0000000000 +x_0=4500000.0000000000 +y_0=0.0000000000 
+a=6377397.155 +rf=299.1528128 +no_defs +nadgrids=/opt/grass/etc/nad/BETA2007.gsb +to_meter=1.0'

gibt das zurück
"Using to definition: k_0=1.0000000000 lat_0=0.0000000000 lon_0=12.0000000000 x_0=4500000.0000000000 
y_0=0.0000000000 a=6377397.155 rf=299.1528128 no_defs nadgrids=/opt/grass/etc/nad/BETA2007.gsb 
to_meter=1.0
Rel. 4.7.1, 23 September 2009
<cs2cs>:
projection initialization failure
cause: projection not named
program abnormally terminated
bash: echo: write error: Datenübergabe unterbrochen (broken pipe)"

Aber hiermit bin ich wohl auch besser auf der PROJ.4-Liste aufgehoben. Wollte hier nur noch 
mitteilen, dass und wie es möglich ist, BETA2007 in GRASS zu verwenden.

Gruß,
Tom

Am 18.10.2010 23:06, schrieb Markus Neteler:
> 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
>>

-- 
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