[Grass-de] GRASS auf Englisch

Tom Gottfried tom.gottfried at wzw.tum.de
Fr Mär 13 16:40:13 CET 2009


Hallo Markus,

danke für die Local-Variablen!

Markus Neteler schrieb:
>> Trotzdem schonmal mein Problem, falls es hier jemand auf Anhieb lösen
>> kann: Wenn ich folgendes tue
>>
>> `r.in.gdal orthophoto.tif out=orthophoto_r band=1'
>>
>> bekomme ich folgende Meldungen:
>> "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?) steht folgendes:

"0.100000000000000010
0.0
0.0
-0.100000000000000010
4487306.9000000004
5361007.0000000000"

Also wohl kein Datum spezifiziert?

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

> - Platte kaputt
> - Karte ueber 2.x GB gross.
> 
>> wobei ab dem Moment, wo die Warnung auf der Konsole erscheint, der Rest
>> des Bildes einfach schwarz gezeichnet wird.
> 
> Ja, da die Datei(en) der Karte nicht vollstaendig sind.
> 
>> Die gleiche Datei, die ich
>> mit `r.in.gdal' importiere, kann ich in QGIS problemlos öffnen.
> 
> Das ist erstaunlich, vielleicht wird der Fehler nicht angezeigt?

Also direkt die Datei kann ich in QGIS öffnen, ich meinte nicht den
GRASS-Layer.

>> Woran könnte es liegen?
> 
> 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
ellipsoid:  bessel
north:      5360976
south:      5359285
west:       4487356
east:       4489047
nsres:      0.5
ewres:      0.5
rows:       3382
cols:       3382
cells:      11437924"

> 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

Danke!
Gruß,
Tom