😏★,°:.☆(￣▽￣)/$:.°★ 😏
这篇文章主要介绍pyproj地理计算库应用。
欢迎来到我的博客，一起学习，共同进步。
喜欢的朋友可以关注一下，下次更新不迷路🥞
1. pyproj介绍
pyproj 是 Python 中一个非常有用的库，用于执行地理坐标系统的转换和地图投影。它是基于 proj 库的 Python 接口，因此继承了 proj 库的强大功能，但以 Python 的形式提供了更友好的接口。
2. 环境配置
安装：pip install pyproj
3. 使用说明
经纬度转UTM坐标：
import pyproj
# 创建UTM坐标系对象
utm_crs = pyproj.CRS.from_epsg(32650) # EPSG 32650对应的是WGS84的UTM Zone 50N
# 创建经纬度坐标系对象
lat_lon_crs = pyproj.CRS('EPSG:4326') # WGS84经纬度坐标系
# 创建转换器
transformer = pyproj.Transformer.from_crs(lat_lon_crs, utm_crs, always_xy=True)
# 定义经纬度坐标
lon, lat = 120.30550, 31.47758
# 将经纬度坐标转换成UTM坐标
x, y = transformer.transform(lon, lat)
print("UTM东坐标：", x)
print("UTM北坐标：", y)
UTM坐标转经纬度：
import pyproj
def utm_to_latlon(easting, northing, zone_number, zone_letter):
proj = pyproj.Proj(proj='utm', zone=zone_number, ellps='WGS84', north=bool(zone_letter >= 'N'))
lon, lat = proj(easting, northing, inverse=True)
return lat, lon
# 示例 UTM 坐标：Zone 18N, Easting: 500000, Northing: 4500000
easting = 464101.487
northing = 4378816.449
zone_number = 49
zone_letter = 'N'
latitude, longitude = utm_to_latlon(easting, northing, zone_number, zone_letter)
print("经度（Longitude）：", longitude)
print("纬度（Latitude）：", latitude)
计算两点之间距离和方位角：
from pyproj import Geod
# 创建 Geod 对象
geod = Geod(ellps="WGS84")
# 定义两个点的经纬度
lon1, lat1 = 12.4924, 41.8902 # 罗马斗兽场
lon2, lat2 = 2.2945, 48.8584 # 巴黎埃菲尔铁塔
# 计算距离和方位角
azimuth1, azimuth2, distance = geod.inv(lon1, lat1, lon2, lat2)
print(f"Distance: {distance} meters")
print(f"Initial Azimuth: {azimuth1} degrees")
print(f"Final Azimuth: {azimuth2} degrees")
以上。
