Greetings
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 =
src_file).strip()
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
or
- 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?
Thanks
Antonio
__________ Information from ESET NOD32 Antivirus, version of virus signature
database 5827 (20110128) __________
The message was checked by ESET NOD32 Antivirus.
http://www.eset.com
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
