[gdal-dev] gdalwrap and crossing 180W/180E

[Ghislain Picard]

  Hi all,
I'm using gdal 1.7.2 and still have the problem of reprojection with 
data crossing 180W/180E.

The data are in lat-long from:
http://www.naturalearthdata.com/downloads/50m-raster-data/50m-shaded-relief/

The image is then cropped as follow:

gdal_merge.py -ul_lr
-180.0 -50.0 180.0 -90.0 -of GTiff -o SR_50_crop.tif SR_50M.tif


and projected in polar stereographic as follow:

gdalwarp -s_srs EPSG:4326
-t_srs EPSG:3031 -of GTiff SR_50M_crop.tif SR_50_stereo.tif


In the reprojected image, a sector of data is missing around the 
180W/180E meridan. This sector is several degree wide although the 
original data was with a resolution of 1/3 degree... surprising.

Is it a bug ? or something wrong with my gdal commands ?
Cheers,
Ghislain

-- 
Ghislain Picard

Laboratoire de Glaciologie et Géophysique de l'Environnement
54 rue Molière - Domaine Universitaire - BP 96
38402 St Martin d'Hères Cedex, FRANCE
tel : +33 476 82 42 42
fax : +33 476 82 42 01

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20101004/20ab8307/attachment.html