Rengifo,<br><br>Can you compare the two contours with the pixel values? Which contours are shifted?<br>FYI, you can read the algorithm rules in the docs of GDALContourGenerate() function[1].<br><br>[1]: <a href="http://www.gdal.org/gdal__alg_8h.html#ceaf98ad40f159cbfb626988c054c085">http://www.gdal.org/gdal__alg_8h.html#ceaf98ad40f159cbfb626988c054c085</a><br>
<br><div class="gmail_quote">On Wed, Aug 31, 2011 at 6:42 PM, Rengifo Ortega <span dir="ltr"><<a href="mailto:rengifoo@yahoo.de">rengifoo@yahoo.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div><div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: tahoma,new york,times,serif; font-size: 12pt;"><div>
</div><div class="MsoNormal"><span lang="EN-US">Dear gdal
users</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US">I hope you
can help me with a question that is puzzling me, since a couple of days. So let
me explain my problem.</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span><span lang="EN-US">2 years ago,
I downloaded a SRTM 90m .asc file into GRASS GIS. </span></div>
<div class="MsoNormal"><span lang="EN-US">A week ago, I
decided to generate some contours lines using gdal_contour command (version FWTools2.4.7)
and compare the contours lines with contours already produced in GRASS GIS. I was quite
disturbed when in fact, I found a horizontal shift of the contours lines to the
south of about -/+35m or to north (unknown distance), with respect to the grass gis- generated contours<br></span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US">Technical aspect of my hardware<br></span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US">I am<span> </span>in Grass GIS <span> </span>MS Windows environment (native using MinGW)
combined with QGIS for visualisation purposes</span></div>
<div class="MsoNormal"><span lang="EN-US">Using FWtools
2.4.7</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US">So, here are
the steps I followed</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US">Part I</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US"><span> </span>-Produced 90m contours lines in GRASS
GIS out of the SRTM DEM 90 by using the r.contour command</span></div>
<div class="MsoNormal"><span lang="EN-US"><span></span>-Visualised the contours lines in QGIS
and saved them as shp format</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US">Part II</span></div>
<div><span lang="EN-US"><span>1)<span style="font: 7pt "Times New Roman";"> </span></span></span><span lang="EN-US">-</span><span lang="EN-US">Export SRTM 90 rasterfile as <span> </span>an ESRI ASCII GRID fomat and converted it to .tif
format using FWTools2.4.7</span></div>
<div><span lang="EN-US"><span>2)<span style="font: 7pt "Times New Roman";"> </span></span></span><span lang="EN-US"><span> </span>Using gdal_contours utility I generated 90
meters contour lines</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US"><span> </span>Results 1<br></span></div>
<div class="MsoNormal"><span lang="EN-US">Comparison
between the GRASS-generated and the GDAL- generated contours lines exhibited a horizontal
shift of -/+ 35m to the south in relation to GRASS GIS generated contours</span></div>
<div class="MsoNormal"><span lang="EN-US"> </span></div>
<div class="MsoNormal"><span lang="EN-US">Part III</span></div>
<div><span lang="EN-US"><span>1<span style="font: 7pt "Times New Roman";"></span></span></span><span lang="EN-US">--Export the raster file SRTM90m from
GRASS GIS; as a Geotiff using gdal</span><span lang="EN-US"><br></span></div><div><span lang="EN-US">-Generate 90m contours using gdal</span></div><div><span lang="EN-US"><br></span></div>
<div><span lang="EN-US"><span><span style="font: 7pt "Times New Roman";"></span></span></span><span lang="EN-US"></span></div>
<div><span lang="EN-US"> </span></div>
<div><span lang="EN-US">Result 2 </span></div>
<div><span lang="EN-US"> </span></div>
<div><span lang="EN-US">Comparing the GRASS GIS contours with the geotiff- generated contours a
displacement to the north was observed, no distance measured.</span></div><div><br><span lang="EN-US"></span></div><div>I hope you can give some hints..</div><div>best regards <br></div><font color="#888888"><div>Rengifo<br>
<span lang="EN-US"></span></div><div><span lang="EN-US"><br></span></div>
</font></div></div><br>_______________________________________________<br>
FWTools mailing list<br>
<a href="mailto:FWTools@lists.maptools.org">FWTools@lists.maptools.org</a><br>
<a href="http://lists.maptools.org/mailman/listinfo/fwtools" target="_blank">http://lists.maptools.org/mailman/listinfo/fwtools</a><br>
<a href="http://fwtools.maptools.org/" target="_blank">http://fwtools.maptools.org/</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Best regards,<br>Chaitanya kumar CH.<br><br>+91-9494447584<br>17.2416N 80.1426E<br>