GISer用于坐标系选择的懒癌神器:utm

本篇介绍一个Python包:utm用于坐标系选择。

1 GIS坐标系

地理信息系统（Geographic Information System，简称GIS）中的坐标系是用于在地球表面定位和表示空间数据的数学模型和参考框架。坐标系的选择和应用对于准确地存储、分析和展示地理空间数据至关重要（来自ChatGPT）。

其实从业务角度关键是区分地理坐标系和投影坐标系。从最直观的角度来说，前者单位是经纬度，后者单位通常是米。因此在计算一些几何指标如距离和面积是必须要转成投影坐标计算。

2 UTM坐标系

UTM（Universal Transverse Mercator，通用横轴墨卡托）坐标系是一种广泛使用的投影坐标系统，用于将地球表面的三维地理坐标（经度和纬度）转换为二维平面坐标（东（Easting）和北（Northing））。它由国际地理测量和地图制图组织（International Association of Geodesy）制定，旨在提供高精度的地图和空间数据表示（来自ChatGPT）。

UTM是我们比较常用的一个坐标系，但是这个坐标系比较适用于小区域研究，因为它在全球是分带的，不同区域需要根据带号区分。那么问题来了，我们怎么知道哪个区域用哪个带号呢？

3 UTM坐标系组成(来自ChatGPT)

UTM坐标系的组成

• 投影带（Zone）：全球共60个投影带，每带宽6度，经度从西经180度开始编号，东经180度结束。每带编号从1到60，东经0度位于第31带。
• 东移量（Easting）：为了避免负数，UTM坐标系在每个投影带的中央经线设置了一个500,000米的假东移量。实际东移量为相对于中央经线的距离。
• 北移量（Northing）：在北半球，赤道设为0米，向北增加。在南半球，赤道设为10,000,000米，向南减少，以避免负数。

这是简单的原理。

4 本篇重点:utm包

之前在找对应坐标系的时候发现了一个Python包：:utm，专门用于自动选择合适的UTM分带坐标系，直接用pip安装即可。

pip install utm

utm github链接

具体使用方法为导入包，给定经纬度即可，这里假设以东经119，北纬40为例。

import utm
utm.from_latlon(40,119)

所以应该选择50N作为带号。

参考链接：

Auto-select suitable UTM Zone based on Grid Intersection