[GRASS5] Re: [bug #1492] (grass) r.in.gdal/GRASS/GDAL: coordinates

Mon Dec 23 04:09:56 EST 2002

On Sun, Dec 22, 2002 at 06:38:15PM +0100, Request Tracker wrote:
> Hi
> 
> This is not unique to GDAL, other software such as ARC/INFO and ERDAS Imagine also define
> UTM for the Southern Hemisphere differently. One should just be aware of this and be able
> to edit the coordinates of a raster file or imagery group, later on.

Hi Wolfgang, hi all,

perhaps we should modify r.in.gdal to take care of this modification
automagically (in case GDAL send a negative UTM zone)?

But I just see that *GDAL* gets it wrong:

gdalinfo p174r73_5t900416_nn1.tif
Driver: GTiff/GeoTIFF
Size is 8191, 7172
Coordinate System is:
PROJCS["WGS 84 / UTM zone 34N",
    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"]],
    AUTHORITY["EPSG","32634"],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",21],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (701128.500000,-1977985.500000)
Pixel Size = (28.500000,-28.500000)
Corner Coordinates:
Upper Left  (  701128.500,-1977985.500) ( 22d53'54.03"E, 17d52'49.72"S)
Lower Left  (  701128.500,-2182387.500) ( 22d55'9.01"E, 19d43'36.50"S)
Upper Right (  934572.000,-1977985.500) ( 25d 5'55.19"E, 17d50'47.82"S)
Lower Right (  934572.000,-2182387.500) ( 25d 8'36.64"E, 19d41'21.06"S)
Center      (  817850.250,-2080186.500) ( 24d 0'54.24"E, 18d47'20.73"S)
Band 1 Block=8191x1 Type=Byte, ColorInterp=Gray
  Overviews: 1024x897

LANDSAT 173/73 is in Botswana, so GDAL should report UTM zone 34S.

I'll post this to the GDAL bugtracker.

 Markus

> > this bug's URL: http://intevation.de/rt/webrt?serial_num=1492
> > -------------------------------------------------------------------------
> > 
> > Subject: r.in.gdal/GRASS/GDAL: coordinates confusion on southern hemisphere (UTM)
> > 
> > grass binary for platform: Compiled from Sources
> > 
> > Hi,
> > 
> > when importing a LANDSAT-TM5 into GRASS
> > Botswana, UTM Zone 35
> > with r.in.gdal,I get a coordinates confusion.
> > A sample dataset is here:
> >
> ftp://ftp.glcf.umiacs.umd.edu/glcf/Landsat/WRS2/p173/r072/p173r72_5t890625.TM-EarthSat-Orthorectified/
> 
> > 
> > (get the 1st channel, Geotiff).
> > 
> > After importing the TM1 channel, the cell header contains:
> > proj:       1
> > zone:       35
> > north:      -1818499.5
> > south:      -2018512.5
> > east:       499120.5
> > west:       269040
> > cols:       8073
> > rows:       7018
> > e-w resol:  28.5
> > n-s resol:  28.5
> > format:     0
> > compressed: 1
> > 
> > But the region definition in "GRASS speak" is for a southern hemisphere
> > area:
> > g.region -pd
> > projection: 1 (UTM)
> > zone:       -35
> > datum:      ** unknown (default: WGS84) **
> > ellipsoid:  wgs84
> > north:      8230695
> > south:      7629900
> > west:       3600
> > east:       525060
> > nsres:      15
> > ewres:      15
> > rows:       40053
> > cols:       34764
> > 
> > Note (the boundary coordinates are not fully overlapping):
> >  - negative zone (GRASS) vs positive zone (GDAL)
> >  - positive N/S coordinates (GRASS) vs negative N/S coordinates (GDAL)
> > 
> > Do we have to fix that in GDAL or in r.in.gdal or is the convention
> > in GRASS for southern hemisphere bad? If the latter, how to solve the
> > problem (I have several data sets...).
> > 
> > Thanks,
> > 
> >  Markus
> > 
> > 
> > -------------------------------------------- Managed by Request Tracker
> > 
> > _______________________________________________
> > grass5 mailing list
> > grass5 at grass.itc.it
> > http://grass.itc.it/mailman/listinfo/grass5
> 
> 
> ---------------------------------------------
> This message was sent using M-Web Airmail.
> JUST LIKE THAT
> Are you ready for 10-digit dialling?
> To find out how this will affect your Internet connection go to www.mweb.co.za/ten
> http://airmail.mweb.co.za/
> 
> 
> --- Headers Follow ---
> 
> >From wolfluck at mweb.co.za  Sun Dec 22 18:38:12 2002
> Return-Path: <wolfluck at mweb.co.za>
> Delivered-To: grass-bugs at lists.intevation.de
> Received: from mail.intevation.de (aktaia [212.95.126.10])
> 	by lists.intevation.de (Postfix) with ESMTP id D2C8C13942
> 	for <grass-bugs at lists.intevation.de>; Sun, 22 Dec 2002 18:38:12 +0100 (CET)
> Received: from rammstein.mweb.co.za (rammstein.mweb.co.za [196.2.53.175])
> 	by mail.intevation.de (Postfix) with ESMTP id E71B51C7F5
> 	for <grass-bugs at intevation.de>; Sun, 22 Dec 2002 18:38:09 +0100 (CET)
> Received: from air.mweb.co.za ([196.2.53.154] helo=Debug)
> 	by rammstein.mweb.co.za with smtp (Exim 3.33 #1)
> 	id 18Q9rL-0003IT-00
> 	for grass-bugs at intevation.de; Sun, 22 Dec 2002 19:25:51 +0200
> To: Request Tracker <grass-bugs at intevation.de>
> From: wolfluck at mweb.co.za
> Subject: Re: [GRASS5] [bug #1492] (grass) r.in.gdal/GRASS/GDAL: coordinates
> Date: Sun, 22 Dec 2002 17:38:05 GMT
> X-Posting-IP: 196.44.133.64
> X-Mailer: Endymion MailMan  
> Message-Id: <E18Q9rL-0003IT-00 at rammstein.mweb.co.za>
> X-Spam-Status: No, hits=1.6 required=5.0
> 	tests=NO_REAL_NAME,QUOTED_EMAIL_TEXT,SPAM_PHRASE_03_05
> 	version=2.43
> X-Spam-Level: *
> 
> -------------------------------------------- Managed by Request Tracker

-- 
Markus Neteler

ITC-irst, Istituto per la Ricerca Scientifica e Tecnologica
      Project on Predictive Models for the Environment    
Via Sommarive, 18        -       38050 Povo (Trento), Italy
tel +39 0461 314 520 (fax -591)           http://mpa.itc.it