【问题标题】:Is there an easy way to create square buffers around point and if they intersect, merge them?有没有一种简单的方法可以在点周围创建方形缓冲区,如果它们相交,合并它们?
【发布时间】:2019-03-10 00:17:07
【问题描述】:

我正在尝试在给定点周围创建方形缓冲区,我能够创建圆形缓冲区但不能创建方形缓冲区。

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[:]

另外,如果我有一些重叠的缓冲区,我该如何合并它们?

【问题讨论】:

    标签: python buffer coordinates geospatial


    【解决方案1】:

    在 jaguardb 中,应该是简单的 select union(buffer(p1), buffer(p2)) from t where intersect(buffer(p1), buffer(p2));

    【讨论】:

      猜你喜欢
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 2022-10-15
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      • 1970-01-01
      相关资源
      最近更新 更多