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

Paul Kelly paul-grass at stjohnspoint.co.uk
Thu Sep 25 10:03:27 EDT 2003


Hello Markus

On Thu, 25 Sep 2003, Markus Neteler wrote:

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

Seems like maybe the location is not set up properly? Did you do it with
g.setproj or r.in.gdal. Could you try editing the DEFAULT_WIND and WIND
files manually to change the zone to -4? Does that help?





More information about the grass-dev mailing list