问题描述
我有一个简单的问题,但我找不到答案 我正在寻找与纬度或经度相关的“minimum_rotated_rectangle”多边形的主角度
df4.minimum_rotated_rectangle
有人有这个库存吗 提前致谢
解决方法
这是一个函数,它采用 minimum_rotated_rectangle
多边形并根据较长的边计算其方位角 (0-180)。
def _azimuth(point1,point2):
"""azimuth between 2 points (interval 0 - 180)"""
import numpy as np
angle = np.arctan2(point2[0] - point1[0],point2[1] - point1[1])
return np.degrees(angle) if angle > 0 else np.degrees(angle) + 180
def _dist(a,b):
"""distance between points"""
import math
return math.hypot(b[0] - a[0],b[1] - a[1])
def azimuth(mrr):
"""azimuth of minimum_rotated_rectangle"""
bbox = list(mrr.exterior.coords)
axis1 = _dist(bbox[0],bbox[3])
axis2 = _dist(bbox[0],bbox[1])
if axis1 <= axis2:
az = _azimuth(bbox[0],bbox[1])
else:
az = _azimuth(bbox[0],bbox[3])
return az
示例:
import geopandas as gpd
df = gpd.read_file(gpd.datasets.get_path('nybb'))
# single geometry
mrr = df.geometry.iloc[0].minimum_rotated_rectangle
azimuth(mrr)
# 66.65508762854085
# whole dataframe
mrrs = df.geometry.apply(lambda geom: geom.minimum_rotated_rectangle)
df['az'] = mrrs.apply(azimuth)
ax = df.plot('az',legend=True)
mrrs.boundary.plot(ax=ax)