<div dir="ltr">George,<div style> That sounds like a resampling issue (either from up/down sampling or in your case from the reprojection stage). During projection, I <b>don't</b> recommend using a nearest neighbor method but a bilinear (or "higher") for DEMs. Nearest neighbor is generally the default for applications because it is fast. Fortunately, GDAL and ArcMap both have options for bilinear or cubic during reprojection. Running a nearest neighbor can generate a stair-step for the DEM (and thus odd slopes).</div>
<div style><br></div><div style>That should solve it...,</div><div style>Trent</div><div style><br></div><div style> </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Feb 11, 2013 at 11:17 PM, George <span dir="ltr"><<a href="mailto:grub80@gmail.com" target="_blank">grub80@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi there,<div><br></div><div>I am using gdaldem to get a slope raster from a DEM (source: ASTER), but it seems like I can't get a correct result. For comparison purposes I did the same steps using Arcgis 10.0. In particular, I followed two procedures:</div>
<div><br></div><div>Procedure 1</div><div>1. Mosaiced different DEM rasters together </div><div>2. Projected the mosaic coordinates to linear </div><div>3. Calculated slope (for gdal: gdaldem slope with no extra options)</div>
<div><br></div><div>The result I get is ok, but there are some linear features that are clearly artificial. ArcGIS gave the same result.</div><div><br></div><div>Procedure 2</div><div>1. Mosaiced different DEM rasters together </div>
<div>2. Calculated slope (for gdal: gdaldem slope -s <span style="font-size:12px;font-family:'Lucida Grande',Verdana,Geneva,Arial,sans-serif">111120</span>)</div><div>3. Projected the slope raster in linear coord.</div>
<div> </div><div>The result is perfect in ArcGIS, while gdaldem returns a rsater that looks all white with some black dots (in greyscale, using the same scale of the ArcGIS one).</div><div><br></div><div>
What am I doing wrong? I am new to the mailing list (and gdal), so I am not sure if it's ok to attach examples. If you need them, I'll be happy to post the rasters.</div><span class="HOEnZb"><font color="#888888"><div>
<br></div><div><br clear="all"><div>George</div>
</div></font></span></div>
<br>_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br></blockquote></div><br></div>