【Python】应用：pyproj地理计算库应用

😏*★,°* :.☆(￣▽￣)/$:.°★ 😏

这篇文章主要介绍pyproj地理计算库应用。
无专精则不能成，无涉猎则不能通。------梁启超

欢迎来到我的博客，一起学习，共同进步。

喜欢的朋友可以关注一下，下次更新不迷路🥞

文章目录

- [:smirk:1. pyproj介绍](#:smirk:1. pyproj介绍)
- [:blush:2. 环境配置](#:blush:2. 环境配置)
- [:satisfied:3. 使用说明](#:satisfied:3. 使用说明)

😏1. pyproj介绍

pyproj 是 Python 中一个非常有用的库，用于执行地理坐标系统的转换和地图投影。它是基于 proj 库的 Python 接口，因此继承了 proj 库的强大功能，但以 Python 的形式提供了更友好的接口。

😊2. 环境配置

安装：pip install pyproj

😆3. 使用说明

经纬度转UTM坐标：

python 复制代码

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坐标转经纬度：

python 复制代码

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)

计算两点之间距离和方位角：

python 复制代码

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")

以上。