[GRASS5] GRASS 5.3/5.7, GDAL and southern hemisphere

Markus Neteler neteler at itc.it
Thu Sep 25 09:37:09 EDT 2003


Today we imported a LANDSAT scene into a southern
hemisphere location (Argentina), UTM. There is a slight
confusion with 'south: defined' and negative zone numbers:

gdalinfo tells us:

GRASS 5.7.0-cvs:~ > gdalinfo tuc6.tif
Driver: GTiff/GeoTIFF
Size is 1024, 692
Coordinate System is:
PROJCS["WGS 84 / UTM zone 4S",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.2572235629972,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",-159],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",10000000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","32704"]]
Origin = (-65.356600,9999973.275500)
Pixel Size = (30.000000,-30.000000)
Corner Coordinates:
Upper Left  (     -65.357, 9999973.275)
Lower Left  (     -65.357, 9979213.275)
Upper Right (   30654.643, 9999973.275)
Lower Right (   30654.643, 9979213.275)
Center      (   15294.643, 9989593.275)
Band 1 Block=1024x1 Type=Byte, ColorInterp=Gray

Importing it:

GRASS 5.7.0-cvs:~ > r.in.gdal in=tuc6.tif out=tuc6
 100%
CREATING SUPPORT FILES FOR tuc6
SETTING GREY COLOR TABLE FOR tuc6 (8bit, full range)

Error on display:

GRASS 5.7.0-cvs:~ > d.rast tuc6
WARNING: [tuc6] in mapset [tucuman] - in different zone [-4] than current
         region [4]
WARNING: unable to open raster map [tuc6 in tucuman]
ERROR: Not able to open cellfile for [tuc6]

How the location was defined:

GRASS 5.7.0-cvs:~ > g.projinfo  -p

-----------------------------------------------------------
PROJ_INFO file:
name:    UTM
datum:   wgs84
towgs84:         0.000,0.000,0.000
proj:    utm
ellps:   wgs84
a:       6378137.0000000000
es:      0.0066943800
f:       298.2572235630
zone:    4
south:   defined
-----------------------------------------------------------
PROJ_UNITS file:
unit:    meter
units:   meters
meters:  1.0
-----------------------------------------------------------

I am not sure where it should be fixed.
Thanks for hints,

 Markus




More information about the grass-dev mailing list