《Python地理空间分析指南 第2版》学习笔记-5.2 方位计算

第5章 Python与地理信息系统

5.2方位计算

需要知道某直线起点和终点构成的方位信息。正北方向为0，顺时针转。

from math import atan2, cos, sin, degrees

lon1 = -90.21
lat1 = 32.31
lon2 = -88.95
lat2 = 30.43
angle = atan2(cos(lat1) * sin(lat2) - sin(lat1) *
              cos(lat2) * cos(lon2 - lon1), sin(lon2 - lon1) * cos(lat2))
bearing = (degrees(angle) + 360) % 360
print(bearing)

308.7992752836875，
1在2的西北方向

---

总结

《Python地理空间分析指南 第2版》学习笔记，仅供学习，如有侵权请联系删除。