Python、GIS 和 Fortran：尝试从 xy 点数据创建多个多边形

Question

我已经研究这个问题有一段时间了，但在 ESRI 论坛页面或我编写的一些 FORTRAN 三角测量脚本上没有找到任何乐趣。

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

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

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

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

Answer 1

这应该适用于插值（假设潮汐线是距岸上某个控制点的 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) 上运行。