我尝试在给定点周围创建正方形缓冲区,我可以创建圆形缓冲区,但不能创建正方形缓冲区
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[:]
另外,如果我有一些重叠的缓冲区,如何合并它们
在jaguardb中,它应该是 从t中选择并集(buffer(p1),buffer(p2)),其中相交(buffer(p1),buffer(p2))
相关问题 更多 >
编程相关推荐