大地坐标与横轴墨卡托投影（UTM）坐标互转原理

一、大地坐标系（Geographic Coordinate System）

以经纬度表示位置，使用椭球体模型（如WGS84）。 纬度 φ（-90°~90°，赤道为0°）； 经度 λ（-180°~180°，本初子午线为0°）。

WGS84椭球参数：

二、UTM坐标系（Universal Transverse Mercator）

投影类型：横轴墨卡托投影（高斯-克吕格投影的全球版）。
分区规则：

• 全球分60个纵向带（每6°经度一个带）；
• 带号计算：$$\text{zone}=\lfloor \frac{\lambda +180}{6}\rfloor +1$$

坐标表示：

• 东距（Easting）：中央经线为500,000米（虚假东移）；
• 北距（Northing）：赤道为0米（南半球为10,000,000米虚假北移）。

三、大地坐标 → UTM坐标（正算）

1. 计算中央经线：$${\lambda }_0=(\text{zone}\times 6-183)$$
2. 转为弧度制：
   $${\varphi }_{\text{rad}}=\varphi \cdot \frac{\pi }{180},\Delta \lambda =(\lambda -{\lambda }_0)\cdot \frac{\pi }{180}$$
3. 卯酉圈曲率半径：$$N=\frac{a}{\sqrt{1-e^2{\sin }^2{\varphi }_{\text{rad}}}}$$
4. 投影公式：
   $$\begin{array}{rl}\text{Easting}& =k_{0}N\cos {\varphi }_{\text{rad}}\cdot \Delta \lambda \cdot \left[1+\frac{\Delta {\lambda }^2}{6}{\cos }^2{\varphi }_{\text{rad}}(1-{\tan }^2{\varphi }_{\text{rad}}+{\eta }^2)\right]+500000\\ \text{Northing}& =k_0\left[{\int }_0^{\varphi }Md\varphi +\frac{\Delta {\lambda }^2}{2}N\sin {\varphi }_{\text{rad}}\cos {\varphi }_{\text{rad}}\right]\end{array}$$
   其中 ${\eta }^2=\frac{e^2}{1-e^2}{\cos }^2{\varphi }_{\text{rad}}$，$k_0=0.9996$ 为比例因子。

四、UTM坐标 → 大地坐标（反算）

1. 去除虚假位移：
   $$E^{\prime }=\text{Easting}-500000$$
   $$N^{\prime }=\{\begin{array}{ll}\text{Northing}& (\text{北半球})\\ \text{Northing}-10^7& (\text{南半球})\end{array}$$
2. 计算起点纬度（迭代）：
   $$\mu =\frac{N^{\prime }}{k_{0}a(1-e^2/4-3e^4/64-5e^6/256)}$$
3. 纬度迭代解算：
   $$\begin{array}{rl}{\phi }_{\text{foot}}=\mu & +\left(\frac{3e_1}{2}-\frac{27e_1^3}{32}\right)\sin (2\mu )\\ & +\left(\frac{21e_1^2}{16}-\frac{55e_1^4}{32}\right)\sin (4\mu )\\ & +\left(\frac{151e_1^3}{96}\right)\sin (6\mu )\end{array}$$
4. 最终纬度/经度计算：
   $$\begin{array}{rl}\rho & =\frac{a(1-e^2)}{(1-e^2{\sin }^2{\phi }_{\text{foot}}{)}^{1.5}}\\ \nu & =\frac{a}{\sqrt{1-e^2{\sin }^2{\phi }_{\text{foot}}}}\\ \phi & ={\phi }_{\text{foot}}-\frac{E^{\prime }\tan {\phi }_{\text{foot}}}{\rho \nu }\left[1-\frac{E^{\prime 2}}{6{\nu }^2}\left(3{\tan }^2{\phi }_{\text{foot}}(1+{\eta }^2)-{\eta }^2\right)\right]\\ \lambda & ={\lambda }_0+\arctan \left(\frac{\sinh (E^{\prime }/(k_0\nu ))}{\cos {\phi }_{\text{foot}}}\right)\end{array}$$

五、特殊处理

1. 跨越带处理：
   a. 边界点需计算两带坐标；
   b. 带号重定义：zone $$=\text{zone}\pm 1$$
2. 极区转换： UTM仅覆盖80°S–84°N，极区使用UPS投影。
3. 中国特殊情况： 国内采用高斯-克吕格投影（3°分带）； 带号计算：$n=\lfloor \frac{L}{3}+0.5\rfloor$，中央子午线 $L_0=3n$。