<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">
</head>
<body text="#000000" bgcolor="#ffffff">
Hi all,<br>
I'm using gdal 1.7.2 and still have the problem of reprojection with
data crossing 180W/180E.<br>
<br>
The data are in lat-long from:<br>
<a class="moz-txt-link-freetext"
href="http://www.naturalearthdata.com/downloads/50m-raster-data/50m-shaded-relief/">http://www.naturalearthdata.com/downloads/50m-raster-data/50m-shaded-relief/</a><br>
<br>
The image is then cropped as follow:<br>
<br>
<div class="moz-text-html" lang="x-western">
<style type="text/css">p, li { white-space: pre-wrap; }</style>
<p style="margin: 0px; text-indent: 0px;">gdal_merge.py -ul_lr<br>
-180.0 -50.0 180.0 -90.0 -of GTiff -o SR_50_crop.tif SR_50M.tif</p>
<br>
and projected in polar stereographic as follow:<br>
<br>
<style type="text/css">p, li { white-space: pre-wrap; }</style>
<p style="margin: 0px; text-indent: 0px;">gdalwarp -s_srs
EPSG:4326<br>
-t_srs EPSG:3031 -of GTiff SR_50M_crop.tif SR_50_stereo.tif</p>
<br>
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.<br>
<br>
Is it a bug ? or something wrong with my gdal commands ?<br>
Cheers,<br>
Ghislain</div>
<pre class="moz-signature" cols="72">--
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
</pre>
</body>
</html>