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

3D編集で維持したいPolyline Featureクラスがあります

MathieuBoonen

<\/HEAD>

皆さん、こんにちは、<\/P>

ガイダンスか奇跡を求めています..<\/P>

<\/P>

最初に作成したときはZ対応のPolyline featureクラスがありますが、それを調整して新しいフィーチャを追加するには、更新されたZ対応のPolylineフィーチャで別のFeatureクラスを作成するしかないようです。そして、静的なデータセットを維持するためにデータベース内のデータを削除して置き換える必要があります。<\/P>

ArcSceneで3D編集を試みましたが、非常に限られた成功しか得られていません。例えば、15頂点、250メートルの線をデジタイズするのに約10分かかります（表示パフォーマンスが遅いです）。そして、多くの場合、線に沿った頂点のZ属性が入力されません。成功したのは一度だけです。（RasterからTINや他のTIFなど、さまざまなElevation Sourcesを試しました）。<\/P>

同じセッションを再起動して翌日にArcSceneドキュメントを再読み込みし、新しいフィーチャを追加すると、また10分かかり、Z対応が失われています。<\/P>

<\/P>

Polyline Layerは標高が3D Surfaceから派生していることを認識していますが、Z値は入力されません。<\/P>

実際にはArcMapでZ対応フィーチャを編集したいです。なぜなら、それは別のリレーショナルデータベースにリンクされているので、ポリラインを半分に分割すると、そのIDもリレーショナルデータベースで適切な新しいレコードと共に更新する必要があるからです。<\/P>

3Dライン自体は、線の勾配を表面上で計算しようとしているため、開始点と終了点のZ値だけが知られていれば十分です。<\/P>

<\/P>

私にはこれは基本的な3D編集のように思えますが、簡単かつ迅速にできないことに驚いています。私が望んでいるのは、表面から補間されたZ値がPolyline自体に保存されることだけです。<\/P>

<\/P>

なぜPolylineフィーチャを変更するのか？それは勾配を決定し、終点間の12％勾配など物理的制約内に保つという意思決定プロセスに関わるからです。XY平面でポリラインを上下に移動するとZ値も変わる必要があるため、この編集プロセス中にできるだけリアルタイムでこれを反映しつつ、リレーショナルデータベースへのリンクも維持する必要があります。<\/P>

これを簡単に達成する方法がわからず困っています。一回のセッションで50〜100回の変更があり得て、一つの勾配変更がネットワーク内の次のポリライン勾配に影響します。<\/P>

<\/P>

ありがとうございます<\/P>

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

タグ付けされた投稿をもっと見る

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.

会話に参加する

サインインしたメンバーは投稿、更新のフォローなどができます。初めてですか？無料アカウントを登録してください。

クイックリンク

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