[GRASS-user] Strange reply or error from g.proj
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
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