[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