[gdal-dev] Reprojection raster 0 360 longitude
laura0
lauramail at iol.it
Fri Jun 7 08:32:43 PDT 2013
Hi all,
I have a geotiff with the following corner coordinates:
Upper Left ( -0.2500000, 90.2500000) ( 0d15' 0.00"W, 90d15' 0.00"N)
Lower Left ( -0.2500000, -90.2500000) ( 0d15' 0.00"W, 90d15' 0.00"S)
Upper Right ( 359.750, 90.250) (359d45' 0.00"E, 90d15' 0.00"N)
Lower Right ( 359.750, -90.250) (359d45' 0.00"E, 90d15' 0.00"S)
Center ( 179.7500000, 0.0000000)
I would like to reproject it in UTM (then I will crop the image) for
instance with the command
gdalwarp -t_srs "+proj=utm +zone=32 +datum=WGS84" temp.tif proj.tif
As you can see above, the range of longitude is 0 360 and not -180 180 so I
have tried:
gdalwarp -t_srs "+proj=utm +zone=32 +datum=WGS84" -te -180 -90 180 90
temp.tif proj.tif
but I get the error message
ERROR 1: Reprojection failed, err = -14, further errors will be supressed
on the transform object.
Creating output file that is 0P x 0L.
I have also tried the option -wo SOURCE_EXTRA=1000 --config CENTER_LONG 0 I
have found somewhere, but I am not sure this is the right way.
Can someone suggest me the way?
Thanks
Laura
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Reprojection-raster-0-360-longitude-tp5058766.html
Sent from the GDAL - Dev mailing list archive at Nabble.com.
More information about the gdal-dev
mailing list