[GRASS-user] Strange reply or error from g.proj

António Rocha antonio.rocha at deimos.com.pt
Fri Jan 28 10:50:30 EST 2011


I have used g.proj command to retrieve current LOCATION proejction
proj_location = grass.read_command('g.proj', flags = 'jf').strip()   
and I obtain:
+proj=utm +south +no_defs +zone=36 +a=6378137 +rf=298.257223563 
+towgs84=0.000,0.000,0.000 +to_meter=1

then I used this expression in gdalwarp to reproject the image. After 
that I was able to import image with r.in.gdal.
But, this temporary file that I created with gdalwarp, if I do
proj_image = grass.read_command('g.proj', flags = 'jf', georef = 
I get :
+proj=utm +south +no_defs +zone=36 +a=6378137 +rf=298.257223563 +to_meter=1
It seems that:
- gdalwarp is not reprojecting it right because it's not using 
+towgs84=0.000,0.000,0.000 but g.proj is not using this parameter
- g.proj is not reading the geoprojection-info in the right way because 
he is considering that projected_image and location projections are the 
same (the difference is +towgs84=0.000,0.000,0.000)

is this an issue or +towgs84=0.000,0.000,0.000 expression is irrelevant?


__________ Information from ESET NOD32 Antivirus, version of virus signature database 5827 (20110128) __________

The message was checked by ESET NOD32 Antivirus.


More information about the grass-user mailing list