[Grass-de] GRASS auf Englisch

[Markus Neteler]

Hallo Tom,

2009/3/13 Tom Gottfried <tom.gottfried at wzw.tum.de>:
> Markus Neteler schrieb:
...
>>> "WARNUNG: Datum <unknown> not recognised by GRASS and no parameters found
>>
>> (welches Datum ist es denn? Das koennen wir vielleicht noch reparieren
>> fuer GRASS 6.4.0).
>
> in der tfw-Datei (oder muss ich woanders schauen?)

ja :)

> steht folgendes:
>
> "0.100000000000000010
> 0.0
> 0.0
> -0.100000000000000010
> 4487306.9000000004
> 5361007.0000000000"
>
> Also wohl kein Datum spezifiziert?

Das wird nicht in der tfw-Datei spezifiziert sondern gglfs in einer
.prj Datei.

>>> Die Projektionsinformationen des Eingabedatensatzes und der aktuellen
>>> Location scheinen übereinzustimmen.
>>> WARNUNG: Fehler beim Schreiben der Null-Reihe 8.
>>> [...diese Zeile zigmal...]
>>
>> ich sehe diese Moeglichkeiten:
>> - Platte voll
>
> genau das war's! Hätt ich auch drauf kommen können, da ich meine
> GISDBASE auf einem SMB-Share liegen hatte :(

Bingo... GRASS kann nur begrenzt Fehler auf fremden Systemen
erkennen, aber es sollte alles melden, was geht. Der Rest waere
eher einer Kristallkugel zu entnehmen...

...
>> Was sagen
>>
>> g.region -p
>
> "projection: 99 (Transverse Mercator)
> zone:       0
> datum:      towgs84=597.1,71.4,412.1,0.894,0.068,-1.563,7.58

Ah ok:
grep 597 datumtransform.table
datumtransform.table:potsdam
"towgs84=597.1,71.4,412.1,0.894,0.068,-1.563,7.58" "Germany (West-
South - 47d00N to 50d20'N)" "Accuracy <1m"

Das Potsdam Datum also.

...
>> g.proj -w
>
> "PROJCS["Transverse Mercator",
>    GEOGCS["bessel",
>        DATUM["unknown",
>            SPHEROID["Bessel_1841",6377397.155,299.1528128],
>            TOWGS84[597.1,71.4,412.1,0.894,0.068,-1.563,7.58]],
>        PRIMEM["Greenwich",0],
>        UNIT["degree",0.0174532925199433]],
>    PROJECTION["Transverse_Mercator"],
>    PARAMETER["latitude_of_origin",0],
>    PARAMETER["central_meridian",12],
>    PARAMETER["scale_factor",1],
>    PARAMETER["false_easting",4500000],
>    PARAMETER["false_northing",0],
>    UNIT["meter",1]]"
>
> Habe das irgendwann mal mit `g.proj -c proj4="+proj=tmerc +a=6377397.155
> +rf=299.1528128000008 +towgs84=597.1,71.4,412.1,0.894,0.068,-1.563,7.58
> +no_defs +lon_0=12 +x_0=4500000"' festgelegt; Datumsinformationen
> müssten Potsdam entsprechen

Welche GRASS Version war das? Typischerweise wird man gefragt, welches
Datum man will (am besten Projektion ueber EPSG Code definieren, dann
kommt der Datumsdialog).
[ GRASS ist da um Klassen besser als QGIS, da es eigene Datumslisten hat
  und v.a. sich nicht bei einem +towgs84 String erschreckt. ]

Mit g.setproj koenntest Du das auf Potsdam setzen, aber der +towgs84 ist
ja drin, daher sollte die Definition auch so in Ordnung sein. Ist also nur
Kosmetik in diesem Fall.

Ciao
Markus

-- 
Markus Neteler
Foundation Edmund Mach (FEM) - Research and Innovation Centre
Environment and Natural Resources Area
GIS and Remote Sensing Unit, Trento, Italy
Web:  http://gis.fem-environment.eu/