nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

how to create smooth dem and lines from points without changing the vertex number?

AYGÜNALBAYRAK

Hi,

I have points and vertices with z values. I want to create smooth lines and smooth dem. However, my each points represent a value which must not be changed.Therefore, I should not create extra point between each point. How can I give this smooth display?

Find more posts tagged with

Accepted answers

XanderBakker

I get your point, the Interpolate Shape—Help | ArcGIS Desktop tool creates a 3D polyline and inserts vertices using the sample distance which normally corresponds to the cellsize of the DEM. This will create the line that follows the surface and does not respect the rule that only at vertices of the line should be taken into consideration.

The way to resolve this, is using some Python code. The code below will iterate all the vertices of the line and create a 3D polyline featureclass (maintaining all the attributes) and convert the vertices to 3D using the elevation from the DEM/DSM:

<SPAN class="keyword token">import</SPAN> arcpy

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">main</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> os
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

&nbsp;&nbsp;&nbsp; fc_in <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\Polyline3D\data.gdb\my2Dpipeline'</SPAN>
&nbsp;&nbsp;&nbsp; fc_out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\Polyline3D\data.gdb\my3Dresult'</SPAN>
&nbsp;&nbsp;&nbsp; dem <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\Polyline3D\data.gdb\DEM'</SPAN>
&nbsp;&nbsp;&nbsp; z_offset <SPAN class="operator token">=</SPAN> <SPAN class="number token">2</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create output featureclass</SPAN>
&nbsp;&nbsp;&nbsp; sr <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>fc_in<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>spatialReference
&nbsp;&nbsp;&nbsp; ws<SPAN class="punctuation token">,</SPAN> fc_name <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN>fc_out<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CreateFeatureclass_management<SPAN class="punctuation token">(</SPAN>ws<SPAN class="punctuation token">,</SPAN> fc_name<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"POLYLINE"</SPAN><SPAN class="punctuation token">,</SPAN> fc_in<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"DISABLED"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"ENABLED"</SPAN><SPAN class="punctuation token">,</SPAN> sr<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># get fields</SPAN>
&nbsp;&nbsp;&nbsp; flds <SPAN class="operator token">=</SPAN> getFields<SPAN class="punctuation token">(</SPAN>fc_in<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># start insert and search cursors</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>InsertCursor<SPAN class="punctuation token">(</SPAN>fc_out<SPAN class="punctuation token">,</SPAN> flds<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> curs_out<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_in<SPAN class="punctuation token">,</SPAN> flds<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> curs<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> curs<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnts_3d <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polyline_in <SPAN class="operator token">=</SPAN> row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> part <SPAN class="keyword token">in</SPAN> polyline_in<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> pnt <SPAN class="keyword token">in</SPAN> part<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; z <SPAN class="operator token">=</SPAN> getValueFromDEM<SPAN class="punctuation token">(</SPAN>dem<SPAN class="punctuation token">,</SPAN> pnt<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt_3d <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">.</SPAN>X<SPAN class="punctuation token">,</SPAN> pnt<SPAN class="punctuation token">.</SPAN>Y<SPAN class="punctuation token">,</SPAN> z <SPAN class="operator token">+</SPAN> z_offset<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnts_3d<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>pnt_3d<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polyline_out <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN>pnts_3d<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> sr<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">True</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; row_out <SPAN class="operator token">=</SPAN> list<SPAN class="punctuation token">(</SPAN>row<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; row_out<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">=</SPAN> polyline_out
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; row_out <SPAN class="operator token">=</SPAN> tuple<SPAN class="punctuation token">(</SPAN>row_out<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; curs_out<SPAN class="punctuation token">.</SPAN>insertRow<SPAN class="punctuation token">(</SPAN>row_out<SPAN class="punctuation token">)</SPAN>


<SPAN class="keyword token">def</SPAN> <SPAN class="token function">getValueFromDEM</SPAN><SPAN class="punctuation token">(</SPAN>dem<SPAN class="punctuation token">,</SPAN> pnt<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; pnt_txt <SPAN class="operator token">=</SPAN> <SPAN class="string token">"{0} {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">.</SPAN>X<SPAN class="punctuation token">,</SPAN> pnt<SPAN class="punctuation token">.</SPAN>Y<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; result <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetCellValue_management<SPAN class="punctuation token">(</SPAN>dem<SPAN class="punctuation token">,</SPAN> pnt_txt<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; elevation <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>result<SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>replace<SPAN class="punctuation token">(</SPAN><SPAN class="string token">','</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'.'</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> elevation

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">getFields</SPAN><SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; flds <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>fld<SPAN class="punctuation token">.</SPAN>name <SPAN class="keyword token">for</SPAN> fld <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFields<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; fld_shp <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>ShapeFieldname
&nbsp;&nbsp;&nbsp; flds<SPAN class="punctuation token">.</SPAN>remove<SPAN class="punctuation token">(</SPAN>fld_shp<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; flds<SPAN class="punctuation token">.</SPAN>insert<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'SHAPE@'</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> flds

<SPAN class="keyword token">if</SPAN> __name__ <SPAN class="operator token">==</SPAN> <SPAN class="string token">'__main__'</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; main<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Below a 3D view of the result:

the red line is the result of the interpolate shape tool
the blue line is the result of the script
the black points are the vertices of the input line.

All comments

DanPatterson_Retired

spline interpolation is one option. interpolation doesn't produce other points, it produces a raster.

Your interpolation's smoothness is only going to be as smooth as the surface you are trying to interpolate and how well you have gathered points so that the interpolation reflects it.

If you sample poorly, then don't expect great results.

AYGÜNALBAYRAK

Hi Dan,

Thank you for your reply. To explain the problem in detailed, I am working with natural gas pipe lines. I take the z values from the surface DEM. I must display my lines smooth and their own vertices. Therefore I need smooth DEM I think. This method that you suggest is suitable for me?

DanPatterson_Retired

You won't know until you run an interpolation. Smoothing interpolatants like spline, will not necessarily return the actual value at the observation point (overshoot and undershoot). You have to decide which is more important and there are dozens of methods of interpolation applicable to their own situations. There is no one, perfect interpolator. Experiment once you understand your data

XanderBakker

In case you have a gas line, the vertices should represent the parts of the construction where for instance the line changes in direction (XY) of elevation (Z). When you extract the Z values at the vertices intermediate (virtual) points should not have an elevation to make the line go straight between the defined vertices. It should not be necessary to smooth the DEM in case the DEM represents correctly the surface.

In your resulting line, do you have additional intermediate points with elevation? These may create an incorrect representation of the line. Is it a design of the gas line or the as built?

AYGÜNALBAYRAK

Hi Xander,

I have the datas without their z values and DEM. Firstly, I gave z values by interpolating DEM surface and pull down the gas lines under the surface according to their default or spesific depth. This method gave me lines which are faulted. I illustrated these with images. The yellow lines represent gas pipe lines. The lines must go on smoothly but I should not create extra vertices because of other infrastructure elements.

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.