沿轮廓线的旗拐点

Question

我有一个三维流中线形状文件，我正试图使用Python和Python3，特别是在标高变化0.5英尺或更高的情况下，沿着ArcGIS Pro标出拐点。以下是我所说的拐点的一些例子：

​

​

我有一张桌子，上面有每一张的顶点和高度。我添加了一个名为“旗帜”的字段来存储哪些顶点是拐点。

​

​

我如何标出海拔上升或下降0.5英尺的点？

Answer 1

我不确定您的数据目前是如何在您的程序中表示的。如果可以将其转换为numpy数组，diff()将计算元素之间的离散差异。

如果列表中有值，则转换为numpy数组非常容易。

my_numpy_array = np.array(my_list)

假设数组y中有y轴值，数组x中有x轴值。

我们可以通过以下途径找到f'(x)

f1 = diff(y)/diff(x)

和f''(x)

f2 = diff(f1) / diff(x)

焊盘使用nan表示这些数组的长度，以便它们与x的长度相同。

f1_p = np.pad(f1, (0, x.shape[0]), mode='constant', constant_values=(np.nan))
f2_p = np.pad(f2, (0, x.shape[0]), mode='constant', constant_values=(np.nan))

然后您可以使用numpy的logical_and来检查f2_p == 0和f1_p > 0.5

flag = np.logical_and(f2_p == 0, f1_p > 0.5)