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

_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to