JFYI http://trac.osgeo.org/grass/ticket/1905
Seems that GRASS fails to notice difference between "north" and "south" at all. Maris. 2013/3/14 Shaun Langley <[email protected]>: > Here is the projection for my input file: > > Driver: GTiff/GeoTIFF > Files: MOD112_2004_001.LST_Day_1km.tif > Size is 2495, 2504 > Coordinate System is: > PROJCS["WGS 84 / UTM zone 36N", > GEOGCS["WGS 84", > DATUM["WGS_1984", > SPHEROID["WGS 84",6378137,298.257223563, > 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",33], > PARAMETER["scale_factor",0.9996], > PARAMETER["false_easting",500000], > PARAMETER["false_northing",0], > UNIT["metre",1, > AUTHORITY["EPSG","9001"]], > AUTHORITY["EPSG","32636"]] > > Here is the projection on which I initialized my GRASS location: > > Driver: GTiff/GeoTIFF > Files: Mode_LULC_T1_500m.tif > Size is 1792, 2082 > Coordinate System is: > PROJCS["WGS_1984_UTM_Zone_36S", > GEOGCS["WGS 84", > DATUM["WGS_1984", > SPHEROID["WGS 84",6378137,298.257223563, > 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",33], > PARAMETER["scale_factor",0.9996], > PARAMETER["false_easting",500000], > PARAMETER["false_northing",10000000], > UNIT["metre",1, > AUTHORITY["EPSG","9001"]], > AUTHORITY["EPSG","32736"]] > > Here is the projection info from GRASS for my location: > > +proj=utm > +south > +no_defs > +zone=36 > +a=6378137 > +rf=298.257223563 > +towgs84=0.000,0.000,0.000 > +to_meter=1 > > Without reprojecting the MODIS image, I can import it into GRASS even though > the projection systems don't match with: > > GRASS 6.4.2 (kenya):~/PullDir/0302568977ukjJXB/001 > r.in.gdal > input=MOD112_2004_001.LST_Night_1km.tif output=foo.testing > Projection of input dataset and current location appear to match > 100% > r.in.gdal complete. Raster map <foo.testing> created. > > version: GRASS 6.4.2 > OS: MacOS X ML (10.8.2) > > Regards, > Shaun > > (sorry for any duplicates, I sent it to the wrong users list) > > -------------------------- > Shaun Langley > Doctoral Candidate > Department of Geography > Michigan State University > (562) 281-7847 > > > On Thu, Mar 14, 2013 at 3:09 AM, Markus Neteler <[email protected]> wrote: >> >> On Thu, Mar 14, 2013 at 2:57 AM, Shaun Langley <[email protected]> >> wrote: >> > I'm wondering if I've stumbled on a bug… >> > >> > I'm working with a bunch of MODIS data that I've converted to UTM (Zone >> > 36, WGS84) using the Modis Reprojection Tool. Using gdal I can see that >> > these maps have projections of EPSG:32636; however my GRASS location is >> > defined with EPSG:32736. When importing into grass with r.in.gdal, GRASS >> > identifies the two projections as the same and allows me to proceed to >> > import without error. When displaying the maps though, it's obvious that >> > the projections were wrong. It too me a long time to figure this problem >> > out and I'm just wondering if maybe it's a bug? >> >> Please post whch GRASS version on which OS you use and >> ideally also the command lines + messages to give better help. >> >> (We do similar things and did not meet such problems) >> >> Markus > > > > _______________________________________________ > grass-user mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/grass-user > _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
