<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 14 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
code
{mso-style-priority:99;
font-family:"Courier New";}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Danilo,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Resampling only occurs when an image is being resized or warped. If the dimensions are not changing then resampling is not really appropriate. If you do need to resample, the utilities called gdal_translate.exe and gdalwarp.exe may be used depending on your needs. Other utilities also provide resample options but they may not fit your purposes.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>As a general rule you should NOT use nearest neighbor to resample elevation (floating point) data because it can drastically alter the elevation values from the original values and creates a staircase like effect in the output imagery. Use something like cubic or bilinear. That said, if the input data is bad then nothing you do will fix it. Maybe look on the web for the USGS 3Dep elevation data if you are in the USA or SRTM 1 arc second data for continuous worldwide coverage. Maybe the source data you are using is just bad or has been altered.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>You should definitely NOT resample the slope output. I think your source data may just be bad but it’s hard to know. Download other sources and compare the output. If the output looks similar then review your code or command line options.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Best regards,<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'>Martin Chapman<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1F497D'><o:p> </o:p></span></p><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma","sans-serif"'> gdal-dev [mailto:gdal-dev-bounces@lists.osgeo.org] <b>On Behalf Of </b>Danilo da Rosa<br><b>Sent:</b> Tuesday, March 17, 2020 2:25 PM<br><b>To:</b> gdal-dev@lists.osgeo.org<br><b>Subject:</b> Re: [gdal-dev] Slope obtained with GDAL has weird lines<o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><div><div><div><p class=MsoNormal>Thanks for your answers Joaquim, Alex and Martin!<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>I think the problem is height is measured by 1m intervals, thatś why I see 1 meter contour lines in the slope result. The dem grid was not interpolated from contours, at least not by me, the data comes from Mapzen terrain tiles, I don't know exactly how they are producing it.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Do you think it would be a good idea to do some kind of interpolation to smooth the DEM file or the slope file? Do you have any recommendations on how to do that using gdal? The idea is to use the gdaldem color-relief command to generate a coloured and easy to read map. The problem is that this lines makes the map more difficult to understand, which is a priority in this case.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Thanks again!<o:p></o:p></p></div><div><p class=MsoNormal>Danilo<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div></div><p class=MsoNormal><o:p> </o:p></p><div><div><p class=MsoNormal>El mar., 17 mar. 2020 a las 4:50, Martin Chapman (<<a href="mailto:mchapman@hyperacktive.com" target="_blank">mchapman@hyperacktive.com</a>>) escribió:<o:p></o:p></p></div><blockquote style='border:none;border-left:solid #CCCCCC 1.0pt;padding:0in 0in 0in 6.0pt;margin-left:4.8pt;margin-right:0in'><div><div><p class=MsoNormal>A slope map indicates the maximum rate of change between a pixel elevation value and its neighboring 8 pixels. In other words, the output value for each pixel represents the steepest slope, or decent to the neighboring elevation values. The output values are given as either percent slope or degrees from 0 (flat) to 90 (vertical) depending on the command line switches you choose. <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Therefore, the output raster values no longer represent elevation values from the input but rather slope (steepness) values and should be interpreted differently. A common way to visually interpret them would be to color code them with a color ramp where each color in the ramp represents a range of degrees somewhere between 0 and 90. That way a user can look at the image and quickly identify steep slopes versus not so steep slopes and flat areas. To create a color ramp that matches changes in steepness you will need to use a statistical method on the slope values such as standard deviations or jenks natural breaks to group similar slope values together and where they change. In other words, you will want your color ramp to change colors where flat areas start to get steeper. I hope that makes sense.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Another way to use the results would be to write an algorithm to identify areas that are within a certain range of steepness so you could identify a good place to build a road or perhaps lay train tracks that must not exceed a certain steepness. <o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>I hope that helps.<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>Best regards,<o:p></o:p></p></div><div><p class=MsoNormal>Martin Chapman<o:p></o:p></p></div><div><p class=MsoNormal><br><br><o:p></o:p></p><p class=MsoNormal style='margin-bottom:12.0pt'>On Mar 16, 2020, at 5:09 PM, Danilo da Rosa <<a href="mailto:daniloide@gmail.com" target="_blank">daniloide@gmail.com</a>> wrote:<o:p></o:p></p></div><blockquote style='margin-top:5.0pt;margin-bottom:5.0pt'><div><p class=MsoNormal><span style='font-family:"Tahoma","sans-serif"'></span><o:p></o:p></p><div><div><p class=MsoNormal>Hi!<o:p></o:p></p></div><div><p class=MsoNormal><o:p> </o:p></p></div><div><p class=MsoNormal>This may be a very basic question, but I can't find the answer anyware.<o:p></o:p></p></div><div><p>I am working with this DEM obtained using <a href="https://pypi.org/project/terraincache/" target="_blank">terraincache</a>, it uses mapzen data. The SRC is WGS 84.<o:p></o:p></p><p><a href="https://i.stack.imgur.com/Sp4Sg.jpg" target="_blank"><span style='text-decoration:none'><img border=0 id="_x0000_i1025" src="https://i.stack.imgur.com/Sp4Sg.jpg" alt="DEM File in QGIS"></span></a><o:p></o:p></p><p>I ran this GDAL command to get the slope:<o:p></o:p></p><p><code><span style='font-size:10.0pt'>gdaldem slope -s 111120 dem.tif slope.tif</span></code> (I had to do scaling because of the SRC)<o:p></o:p></p><p>The resulting image looks fine but it's not smooth, it has weird lines over it, like contour lines.<o:p></o:p></p><p><a href="https://i.stack.imgur.com/POcHR.jpg" target="_blank"><span style='text-decoration:none'><img border=0 id="_x0000_i1026" src="https://i.stack.imgur.com/POcHR.jpg" alt="Slope file on QGIS"></span></a><o:p></o:p></p><p>You can download the DEM file here: <a href="https://gofile.io/?c=Du7WnM" target="_blank">https://gofile.io/?c=Du7WnM</a><o:p></o:p></p><p>Does anyone have an idea of how to solve it? <o:p></o:p></p><p>Thanks in advance for your help,<o:p></o:p></p><p>Danilo<o:p></o:p></p><p>P.S. I made this question here too: <a href="https://gis.stackexchange.com/questions/353983/slope-obtained-with-gdal-has-weird-lines" target="_blank">https://gis.stackexchange.com/questions/353983/slope-obtained-with-gdal-has-weird-lines</a><o:p></o:p></p></div></div><p class=MsoNormal>_______________________________________________<br>gdal-dev mailing list<br><a href="mailto:gdal-dev@lists.osgeo.org" target="_blank">gdal-dev@lists.osgeo.org</a><br><a href="https://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">https://lists.osgeo.org/mailman/listinfo/gdal-dev</a><o:p></o:p></p></div></blockquote></div></blockquote></div></div></div></body></html>