[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