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; } } }

Mám třídu prvků Polyline, kterou chci udržovat v 3D úpravách

MathieuBoonen

<\/HEAD>

Ahoj všichni,<\/P>

Hledám radu nebo zázrak..<\/P>

<\/P>

Mám třídu prvků Polyline, která je při prvotním vytvoření Z aware, ale pro úpravy a přidávání nových prvků do ní. Zdá se, že mohu vytvořit pouze další třídu prvků s aktualizovanými Z aware Polyline prvky a obecně musím mazat a nahrazovat data v databázi, abych udržel statický dataset.<\/P>

Zkoušel jsem ArcScene pro úpravy v 3D s velmi omezeným úspěchem, tj. trvá asi 10 minut nakreslit čáru o 15 vrcholech dlouhou 250 metrů (výkon zobrazení je pomalý). A často nezaplní atribut Z vrcholů podél čáry, podařilo se mi to jen jednou. (Zkoušel jsem různé Elevation Sources od Rasteru po TIN a několik dalších TIF atd.).<\/P>

Pokud restartuji stejnou relaci a načtu dokument ArcScene další den a přidám nový prvek, opět to trvá 10 minut a ztratí se povědomí o Z,<\/P>

<\/P>

Vrstva Polyline ví, že výšky jsou odvozeny z 3D povrchu, ale nezaplní hodnotu Z.<\/P>

Ve skutečnosti bych raději upravoval Z aware prvky v ArcMap, protože je propojen s jinou relační databází, takže pokud rozdělím polyline na polovinu, jeho ID musí být také aktualizováno v relační databázi s přidáním odpovídajících nových záznamů.<\/P>

Sama 3D čára potřebuje znát pouze hodnoty Z počátečního a koncového bodu, protože se snažím vypočítat Gradient čáry přes povrch.<\/P>

<\/P>

Pro mě to zní jako základní 3D úpravy, ale jsem překvapený, že to nemohu udělat snadno nebo rychle. Vše, co chci, je uložit interpolovanou hodnotu Z z povrchu přímo s Polyline.<\/P>

<\/P>

Proč měnit Polyline prvky se ptáte? Je to celé o rozhodovacích procesech určování Gradientu a udržení gradientu v rámci fyzických omezení jako je 12% sklon mezi koncovými body. Pokud posunete polyline nahoru nebo dolů na povrchu v XY, hodnota Z se musí také změnit, takže potřebuji toto reflektovat co nejvíce v reálném čase během procesu úprav při zachování vazeb na relační databázi.<\/P>

Jsem prostě bezradný, jak toho snadno dosáhnout. Úpravy by mohly zahrnovat 50 -100 změn během jedné relace a jedna změna gradientu ovlivňuje další gradient polyline v síti liniových prvků.<\/P>

<\/P>

Díky<\/P>

Mat Boonen , Hikurangi . NZ<\/P><\/BODY><\/HTML>

Najděte více příspěvků označených jako

Accepted answers

XanderBakker

In case you want to use the FieldCalculator you could do this:

Select Python as parser
Activate the option "Show Codeblock"
Paste the code below in the codeblock
change the path to your DEM (see line 3 of the code below)
Specify "getSlope( !SHAPE! )" to assign to your output field (in my case the field " test")

Code block:

import arcpy
# change the path to your DEM
dem = "D:/Xander/GeoNet/Pasture/ned10m34101g3.tiff"

def createExtentForPoint(pnt, cellsize):
    return arcpy.Extent(pnt.X - cellsize, pnt.Y - cellsize, pnt.X + cellsize, pnt.Y + cellsize)

def getRowColExtentFromXY(ext, pnt, cellsize):
    col = int(((pnt.X - ext.XMin) - ((pnt.X - ext.XMin) % cellsize)) / cellsize)
    row = int(((ext.YMax - pnt.Y) - ((ext.YMax - pnt.Y) % cellsize)) / cellsize)
    return row, col

def createNumpyArrayUsingExtent(raster, ext, cellsize):
    lowerLeft = arcpy.Point(ext.XMin, ext.YMin)
    ncols = int(ext.width / cellsize)
    nrows = int(ext.height / cellsize)
    return arcpy.RasterToNumPyArray(raster, lowerLeft, ncols, nrows)

def adaptExtentUsingRaster(ext, raster, cellsize):
    ras_ext = arcpy.Describe(raster).extent
    xmin = ext.XMin - ((ext.XMin - ras_ext.XMin) % cellsize)
    ymin = ext.YMin - ((ext.YMin - ras_ext.YMin) % cellsize)
    xmax = ext.XMax + ((ras_ext.XMax - ext.XMax) % cellsize)
    ymax = ext.YMax + ((ras_ext.YMax - ext.YMax) % cellsize)
    return arcpy.Extent(xmin, ymin, xmax, ymax)

def getRasterCellSize(raster):
    desc = arcpy.Describe(raster)
    return (desc.meanCellHeight + desc.meanCellWidth) / 2

def getSlope(polyline):
    cellsize = getRasterCellSize(dem)
    pnts = [polyline.firstPoint, polyline.lastPoint]
    elevs = []
    for pnt in pnts:
        ext1 = createExtentForPoint(pnt, cellsize)
        ext2 = adaptExtentUsingRaster(ext1, dem, cellsize)
        np = createNumpyArrayUsingExtent(dem, ext2, cellsize)
        row, col = getRowColExtentFromXY(ext2, pnt, cellsize)
        z = np.item(row,col)
        elevs.append(z)
    zdif = elevs[1] - elevs[0]
    slope = float(zdif / polyline.length)
    return slope

This was the result on my test set:

You can use the numeric format of the field to represent the calculated fraction as percentage.

All comments

MathieuBoonen

I know what you mean, the truth is engineers place the line at strategic locations to minimize the possibility of creating the variance that you describe, and generally these sections represent the start and finish height of each section of road. to create an overall gradient between the two locations, the peice in the middle is just the effort amount to get to the finished requirement, but we first need to know that the gradient is not too steep overall, the balancing (or lack thereof)of cuts and fills is really the nuts and bolts of getting the grade from one point to the next right.

I have Arc Pro, but need to come up to speed with that too As 3D editing is not so obvious as yet.

I have used the Guidelines for editing features with base height properties in 3D—Help | ArcGIS for Desktop

And yet still am up against it.

XanderBakker

Personally I would be more interested in knowing the elevation difference in the entire trajectory then only at the start and end positions, especially in cases like the one below:

However, I you want to extract start and end height, determine the difference between them and divide it by the length to obtain the slope, this could be scripted with some Python code. This is much faster than editing in ArcScene in a 3D environment.

You may also want to consider editing in ArcGIS Pro (if you have access to it) since it has enhanced performance for 3D editing workflows.

Guidelines for editing features with base height properties in 3D—Help | ArcGIS for Desktop

http://desktop.arcgis.com/en/desktop/latest/guide-books/extensions/3d-analyst/about-setting-feature-heights-z-s-for-3d-e…

MathieuBoonen

Thanks Xander,

But Sadly that only finds Max and Min plus Mean Z it will not interpolate the Z at the beginning and end of the Polyline which is what i need. generally the Max and Min are not always at the start or finish of each line.

I have got a little further with using ArcScene in that if I have a very small surface model (250 hectares in size) it seems to be a lot more stable in that environment. Strictly speaking I dont even wish to render the surface model in 3d anyway . a 2D view to digitize is fine, I just want the extract Z heights as above.

I could break my surface model into hundred of tiles but that also seems a bit of a compromise. especially when trying to associate the Polyline with the correct surface model under it.

XanderBakker

Do you know the "Add Surface Information (3D Analyst)" tool? It will add this information to a polyline based on a surface.

For a polyline it will add:

3D distance of the line along the surface.
Minimum, maximum, and mean elevation and slope of the line along the surface.

Připojte se k diskusi

Přihlášení členové mohou přispívat, sledovat aktualizace a další. Jste tu noví? Zaregistrujte si bezplatný účet.

Rychlé odkazy

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