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

J'ai une classe d'entités Polyline que je souhaite maintenir en édition 3D

MathieuBoonen

<\/HEAD>

Bonjour à tous,<\/P>

Je cherche des conseils ou un miracle..<\/P>

<\/P>

J'ai une classe d'entités Polyline qui est consciente de Z lors de sa création initiale mais pour ajuster et ajouter de nouvelles entités. Il semble que je ne puisse créer qu'une autre classe d'entités avec les entités Polyline mises à jour et conscientes de Z, et généralement devoir supprimer et remplacer les données dans la base de données pour maintenir l'ensemble de données statique.<\/P>

J'ai essayé ArcScene pour éditer en 3D avec un succès très limité, c'est-à-dire qu'il faut environ 10 minutes pour numériser une ligne de 15 sommets longue de 250 mètres (la performance d'affichage est lente). Et souvent, il ne remplit pas l'attribut Z des sommets le long de la ligne, je n'ai réussi qu'une seule fois. (J'ai essayé différentes sources d'élévation allant du Raster au TIN et quelques autres TIF, etc.).<\/P>

Si je redémarre la même session et recharge le document ArcScene le lendemain et ajoute une nouvelle entité, cela prend encore 10 minutes et la conscience Z est perdue,<\/P>

<\/P>

La couche Polyline sait que les élévations sont dérivées de la surface 3D mais elle ne remplit pas la valeur Z.<\/P>

Je préférerais en fait éditer les entités conscientes de Z dans ArcMap car elle est liée à une autre base de données relationnelle, donc si je divise une polyline en deux, son ID doit également être mis à jour dans la base relationnelle avec les nouveaux enregistrements appropriés ajoutés.<\/P>

La ligne 3D elle-même n'a vraiment besoin que des valeurs Z du point de départ et du point d'arrivée car j'essaie de calculer la pente de la ligne sur la surface.<\/P>

<\/P>

Pour moi, cela ressemble à une édition 3D basique mais je suis surpris de ne pas pouvoir le faire facilement ou rapidement. Tout ce que je veux, c'est que la valeur Z interpolée depuis la surface soit stockée avec la Polyline elle-même.<\/P>

<\/P>

Pourquoi modifier les entités Polyline, demandez-vous ? Eh bien, tout est question de processus décisionnels pour déterminer la pente et maintenir la pente dans des contraintes physiques telles qu'une inclinaison de 12 % entre les points finaux. Si vous déplacez la polyline vers le haut ou vers le bas sur la surface en XY, la valeur Z doit aussi changer, donc je dois refléter cela aussi près que possible du temps réel pendant le processus d'édition, tout en maintenant les liens avec la base relationnelle.<\/P>

Je suis simplement perdu quant à comment réaliser cela facilement. Les modifications pourraient comprendre 50 à 100 altérations dans une seule session, et un changement de pente affecte la pente suivante dans le réseau d'entités linéaires.<\/P>

<\/P>

Merci<\/P>

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

Trouver plus de publications taguées avec

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.

Rejoindre la conversation

Les membres connectés peuvent publier, suivre les mises à jour, et plus encore. Nouveau ici ? Inscrivez-vous gratuitement.

Liens rapides

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