Есть ли простой способ создать квадратные буферы вокруг точки и, если они пересекаются, объединить их? - PullRequest
0 голосов
/ 10 марта 2019

Я пытаюсь создать квадратные буферы вокруг заданных точек, я могу создавать круглые буферы, но не квадратные.

from shapely.ops import transform
from shapely.geometry import Point
    local_azimuthal_projection = "+proj=aeqd +R=6371000 +units=m +lat_0={} +lon_0={}".format(lat, lon)
                    wgs84_to_aeqd = partial(
                        pyproj.transform,
                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),
                        pyproj.Proj(local_azimuthal_projection),
                    )
                    aeqd_to_wgs84 = partial(
                        pyproj.transform,
                        pyproj.Proj(local_azimuthal_projection),
                        pyproj.Proj('+proj=longlat +datum=WGS84 +no_defs'),
                    )

                    point_transformed = transform(wgs84_to_aeqd, Point(float(lon), float(lat)))
                    buffer = point_transformed.buffer(0.5*1000)
                    buffered_geom = transform(aeqd_to_wgs84, buffer).exterior.coords[:]

Кроме того, если у меня есть несколько буферов, которые перекрываются, как я могу их объединить

1 Ответ

0 голосов
/ 10 марта 2019

В jaguardb это должно быть просто Выберите объединение (буфер (p1), буфер (p2)) из t, где пересекаются (буфер (p1), буфер (p2));

...