使用 gma 绘制全球 UTM 投影分带

步骤0：UTM 投影简介

UTM（Universal Transverse Mercator，通用横轴墨卡托）投影是一种广泛应用于全球地图制图、GIS（地理信息系统）和工程测量的横轴墨卡托投影。它采用分带投影的方式，将地球划分为多个经度带，以减少投影变形，确保每带内的地图精度。

本文所用到 UTM 投影为 WGS 84 / UTM 投影 6°分带，从180°W ~ 180°E，每分带 6°，对应的EPSG编码为 32601 ~ 32660（北半球，0° ~ 84°N）、 32701 ~ 32760（南半球，80°S~0°）。

步骤1：准备

1.1 环境

Python: 3.12
gma: 2.1.5.2（下载：https://pypi.org/project/gma/）
安装：
pip install gma==2.1.5.2

1.2 使用 gma 查看一个 UTM 分带投影信息

（以 EPSG:32651为例；Jupyter Notebook）

from gma import crs
proj = crs.SpatRef(32651)

在 Jupyter Notebook 中，运行 proj 单元格可获得如下信息：

步骤2：绘制

from gma.map import plot, inres
from gma import crs
import matplotlib.pyplot as plt

## 0 创建地图框
MapF = plot.MapFrame()
## 1 添加底图
MapL0 = MapF.AddDataSetDiscrete(inres.WorldDataSet.NaturalEarth)

## 2 添加 UTM 投影分带范围
## 2.1 WGS 84 / UTM - N
for epsgID in range(32601, 32661):
    proj = crs.SpatRef(epsgID)
    ## 获得投影范围
    proj_ft = proj.GetUseFeatureInWGS84()
    ## 添加投影范围要素
    MapLn = MapF.AddFeature(proj_ft, FaceColor = 'none', LineColor = 'red', LineWidth = 0.3)
    ## 添加标签
    MapLn.AddLabel(proj.Name.split()[-1], FontSize=4, Offset=(-5, 0), rotation=-90)
## 2.2 WGS 84 / UTM - S (重复操作)
for epsgID in range(32701, 32761):
    proj = crs.SpatRef(epsgID)
    proj_ft = proj.GetUseFeatureInWGS84()
    MapLn = MapF.AddFeature(proj_ft, FaceColor = 'none', LineColor = 'blue', LineWidth = 0.3)
    MapLn.AddLabel(proj.Name.split()[-1], FontSize=4, Offset=(-5, 0), rotation=-90)

## 3 其他的一些要素
GridLines = MapF.AddGridLines(LONRange=(-180, 180, 6), LATRange=(-90, 90, 90), LineWidth = 0.2)

MapF.Axes.text(-181, 0, 0, size = 5, ha = 'right', va = 'center')
MapF.Axes.text(-181, 84, '84°N', size = 5, ha = 'right', va = 'center')
MapF.Axes.text(-181, -80, '80°S', size = 5, ha = 'right', va = 'center')
for lon in range(-180, 181, 6):
    tx = f'{abs(lon)}°W' if lon < 0 else f'{lon}°E' if lon > 0 else lon
    MapF.Axes.text(lon, -92, tx, size = 4, ha = 'center', va = 'top', rotation=-90)
    
MapF.Axes.text(0, 95, 'WGS 84 / UTM zone', size = 6, ha = 'center', va = 'center')

plt.savefig(r"UTM 分带.jpg", bbox_inches = 'tight')

结果如下：