19、地理空间数据处理：pyproj与栅格数据全解析

地理空间数据处理：pyproj与栅格数据全解析

1. pyproj空间参考的使用

在地理空间数据处理中，pyproj是一个强大的工具，它与PROJ.4配合使用，和OSR不同，它处理的是坐标值列表，这些列表可以是Python列表、元组、数组、NumPy数组或标量。如果你有一个包含坐标的文本文件，pyproj模块中的函数是将其转换到其他坐标系的理想选择。你可以在 这里 找到pyproj模块的在线文档和下载资源。

1.1 空间参考系统间的坐标转换

使用pyproj在空间参考系统之间转换坐标有几种不同的方法。可以使用Proj类在地理坐标和投影坐标之间进行转换，也可以使用模块级的transform函数在两个空间参考系统之间进行转换。

示例：将埃菲尔铁塔的经纬度坐标转换为UTM 31N带坐标

import pyproj
utm_proj = pyproj.Proj('+proj=utm +zone=31 +ellps=WGS84')
x, y = utm_proj(2.294694, 48.858093)
print(x, y)

这里将单个x和y坐标传递给utm_proj，它会返回一个x和一个y。你也可以传递x值列表和y值列表（其中x[i]和y[i]是一个坐标对），然后会返回两个列表。

若要从投影坐标转换回地理坐标，将可选的inverse参数设置为True并传入UT