已知中心点经纬度，求笛卡尔网格转经纬度的Python实现方案

要将以参考点为原点的笛卡尔米网格转换为WGS84经纬度，最可靠的方式是使用局部横轴墨卡托投影（适合小范围区域，误差可忽略），或者基于UTM投影进行坐标偏移转换。以下是两种可行的实现方式：

方法一：自定义局部投影（推荐） #

直接创建以参考经纬度为原点的局部横轴墨卡托投影，将网格坐标直接逆转换为经纬度：

import numpy as np
import pyproj

# 生成笛卡尔网格
x = np.linspace(-3600, 3600, 7201)
y = np.linspace(-3600, 3600, 7201)
yy, xx = np.meshgrid(y, x)

# 参考经纬度（WGS84）
referenceLat = 48.85865
referenceLon = 2.33811

# 创建自定义局部横轴墨卡托投影，原点为参考点
local_proj = pyproj.Proj(
    f'+proj=tmerc +lat_0={referenceLat} +lon_0={referenceLon} '
    '+ellps=WGS84 +units=m +no_defs'
)

# 将局部笛卡尔坐标转换为WGS84经纬度
lon_grid, lat_grid = local_proj(xx, yy, inverse=True)

# 验证中心点坐标（应与参考经纬度一致）
print(f"中心点经纬度：{lat_grid[3600, 3600]:.6f}, {lon_grid[3600, 3600]:.6f}")

方法二：UTM投影偏移转换 #

先将参考点转换为对应UTM带的坐标，再通过偏移计算实际UTM坐标，最后转回经纬度：

import numpy as np
import pyproj

# 生成笛卡尔网格
x = np.linspace(-3600, 3600, 7201)
y = np.linspace(-3600, 3600, 7201)
yy, xx = np.meshgrid(y, x)

# 参考经纬度（WGS84）
referenceLat = 48.85865
referenceLon = 2.33811

# 定义投影：WGS84经纬度（EPSG:4326）和UTM31N（巴黎所在UTM带，EPSG:32631）
wgs84 = pyproj.Proj('EPSG:4326')
utm31n = pyproj.Proj('EPSG:32631')

# 将参考点转换为UTM坐标
x0, y0 = pyproj.transform(wgs84, utm31n, referenceLon, referenceLat)

# 计算实际UTM坐标（网格偏移+参考点UTM坐标）
utm_x = xx + x0
utm_y = yy + y0

# 转换回WGS84经纬度
lon_grid, lat_grid = pyproj.transform(utm31n, wgs84, utm_x, utm_y)

# 验证中心点坐标
print(f"中心点经纬度：{lat_grid[3600, 3600]:.6f}, {lon_grid[3600, 3600]:.6f}")

关键说明 #

• 自定义局部投影适用于小范围区域（如本次7.2km见方的网格），投影误差极小，无需额外计算偏移。
• UTM方法需提前确定参考点对应的UTM带（本次参考点位于巴黎，属于UTM 31N，对应EPSG:32631），若你的参考点在其他区域，需替换对应的EPSG编码。
• 两种方法都可保证中心点(0,0)精确对应给定的参考经纬度，网格其他点的转换结果符合地理坐标逻辑。

内容的提问来源于stack exchange，提问作者Rakesh