Python, GIS и Fortran: Попытка создать несколько полигонов из данных точек xy - PullRequest
Новая
       35

Python, GIS и Fortran: Попытка создать несколько полигонов из данных точек xy

1 голос
/ 22 сентября 2011

Я работал над этой проблемой некоторое время и не нашел радости на странице форума ESRI или с каким-то написанным мной сценарием триангуляции FORTRAN.

У меня есть два файла .csv с сотнями данных точек xy. Эти точки представляют верхний и нижний пределы литорального диапазона. Верхняя и нижняя точки проходят параллельно друг другу, и я хочу создать осколки полигонов, которые соединяют четыре из этих точек каждой ногой в отдельные многоугольники. Высота многоугольника будет равна x в зависимости от расстояния между верхней и нижней точками. Ниже ссылка показывает два изображения, которые иллюстрируют, что я имею в виду:

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

Основная проблема заключалась в написании сценариев для правильного формирования углов в углах. Я понимаю, что вы не можете иметь многоугольник диаметром 1 фут в нижней части и диаметром 1 фут в верхней части при движении вокруг поворота. Но это только одна из многих проблем, с которыми я столкнулся при попытке решить эту проблему ...

Любая помощь будет принята с благодарностью, Спасибо

1 Ответ

0 голосов
/ 23 сентября 2011

Это должно работать для интерполяции (скажем, линия прилива 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).

...