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

Ik heb een Polyline Feature class die ik wil behouden in 3D Editing

MathieuBoonen

<\/HEAD>

Hallo allemaal,<\/P>

Ik ben op zoek naar begeleiding of een wonder..<\/P>

<\/P>

Ik heb een Polyline feature class die Z-bewust is wanneer deze voor het eerst wordt gemaakt, maar om deze aan te passen en nieuwe features toe te voegen. Het lijkt erop dat ik alleen een andere Feature class kan maken met de bijgewerkte Z-bewuste Polyline features. en over het algemeen moet ik data in de database verwijderen en vervangen om de statische dataset te behouden.<\/P>

Ik heb geprobeerd ArcScene te gebruiken om in 3d te bewerken met zeer beperkte succes, d.w.z. het duurt ongeveer 10 minuten om een lijn van 15 vertices van 250 meter lang te digitaliseren (weergaveprestaties zijn traag). En dan vult het vaak niet de Z-attribuut van de vertices langs de lijn in, ik ben slechts één keer succesvol geweest. (Ik heb verschillende Elevation Sources geprobeerd van Raster tot TIN en een paar andere TIF enz.).<\/P>

Als ik dezelfde sessie opnieuw start en het ArcScene-document de volgende dag opnieuw laad en een nieuwe feature toevoeg, duurt het weer 10 minuten en is de z-bewustheid verloren,<\/P>

<\/P>

De Polyline Layer weet dat de hoogtes zijn afgeleid van het 3D Surface maar het vult de Z-waarde niet in.<\/P>

Ik zou eigenlijk liever de Z-bewuste features bewerken in ArcMap omdat het gekoppeld is aan een andere relationele database, dus als ik een polyline doormidden splits, moet ook het ID worden bijgewerkt in de relationele database met passende nieuwe records toegevoegd.<\/P>

De 3D-lijn zelf hoeft eigenlijk alleen de Startpunt- en Eindpunt-Z-waarden te kennen omdat ik probeer de Gradient van de lijn over het oppervlak te berekenen.<\/P>

<\/P>

Voor mij klinkt dit als basis 3D-bewerking, maar ik ben verrast dat ik het niet gemakkelijk of snel kan doen. want alles wat ik wil is dat de Z-waarde geïnterpoleerd van het oppervlak wordt opgeslagen bij de Polyline zelf.<\/P>

<\/P>

Waarom Polyline features wijzigen vraag je? Nou, het gaat allemaal om besluitvormingsprocessen voor het bepalen van Gradient en het binnen fysieke beperkingen houden zoals 12% helling tussen eindpunten. Als je de polyline omhoog of omlaag verplaatst op het oppervlak in XY, moet ook de Z-waarde veranderen, dus ik moet dit zo dicht mogelijk bij realtime tijdens het bewerkingsproces weerspiegelen, terwijl ik de koppelingen naar de relationele database behoud.<\/P>

Ik weet gewoon niet hoe ik dit gemakkelijk kan bereiken. De wijzigingen kunnen 50 -100 aanpassingen omvatten in één sessie, en één gradientverandering beïnvloedt de volgende polyline gradient in het netwerk van lijnfeatures.<\/P>

<\/P>

Bedankt<\/P>

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

Vind meer berichten getagd met

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.

Doe mee aan het gesprek

Aangemelde leden kunnen berichten plaatsen, updates volgen en meer. Nieuw hier? Registreer een gratis account.

Snelle Links

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