1

我一直在研究这个问题一段时间，在 ESRI 论坛页面或我编写的一些 FORTRAN 三角剖分脚本上没有发现任何乐趣。

我有两个 .csv 文件，其中包含数百个 xy 点数据。这些点代表潮间带范围的高端和低端。高点和低点彼此平行，我想创建多边形条，将每英尺的四个点连接成单独的多边形。取决于高点和低点之间的距离，多边形的高度将为 x。下面的链接显示了两张图片来说明我的意思：

http://forums.arcgis.com/threads/39757-Feature-to-Line..?p=135880&posted=1#post135880

主要问题是编写多边形脚本以在角落正确形成。我知道当绕弯移动时，你不能有一个底部直径为 1 英尺，顶部直径为 1 英尺的多边形。但这只是我在尝试解决这个问题时遇到的许多问题之一......

任何帮助将不胜感激，谢谢

4

1 回答 1

0

这应该适用于插值（比如潮汐线是距离岸上某个控制点的 x，y 距离）：

import math

example_triangle = [[0,0], [5,5], [0,5], [0,0]]
high_tide_line = [[0, 0], [5., 1.5], [10., 3.2], [20., 1.], [30., 4.5], [80.,2.], [80,0]]
low_tide_line = [[0, 10.], [5., 11.5], [10., 13.2], [20., 11.], [30., 14.5], [80., 12.], [80, 10]]

def points_from_geom(geom):
    idx = 0
    line_lengths = []
    unit_vectors = []
    interpolated_points = []
    while idx < (len(geom) - 1):
        dy, dx = ((geom[idx+1][1] - geom[idx][1]), (geom[idx+1][0] - geom[idx][0]))
        line_lengths.append(math.sqrt(dy**2 + dx**2))
        try:
            angle = math.atan(dy/dx)
            unit_vectors.append([math.cos(angle)*cmp(dx, 0),
                math.sin(angle)*cmp(dy, 0)])
        except ZeroDivisionError:
            if geom[idx+1][1] < geom[idx][1]:
                direction = [0, -1]
            else:
                direction = [0, 1]
            unit_vectors.append(direction)
        idx += 1

    for i, length in enumerate(line_lengths):
        inter = 0
        while inter <= length:
            interpolated_points.append([geom[i][0] + unit_vectors[i][0]*inter,\
                geom[i][1] + unit_vectors[i][1]*inter])
            inter += .3048 # a ft in proper units ;)

    return interpolated_points

ln1 = points_from_geom(example_triangle)
ln2 = points_from_geom(high_tide_line)
ln3 = points_from_geom(low_tide_line)

print ln1, ln2, ln3

有很多不同的方法来做多边形。我可能会尝试确定参考海岸线，然后以固定间隔或在每个线段的中间取垂直线，然后在相邻海岸线上的交点处制作多边形或边界框。顺便说一句，在实际应用程序中，您需要确保cmp不在 (0,0) 上运行。

于 2011-09-23T11:09:09.290 回答