UTM2LL将Universal横向墨卡托(UTM)东/北坐标转换为纬度/经度。
LL2UTM将纬度/经度坐标转换为UTM。
这两个函数都使用精确的公式(毫米精度)，可能的用户定义数据(WGS84是默认值)，并且都是向量化的(代码中没有循环)。这意味着大量的点矩阵，比如整个DEM网格，可以非常快地转换。
(需要readhgt示例。作者的函数):
X = readhgt(36:38 12:15,“合并”,“作物”,(36.5,38.5,12.2,16),“阴谋”);
(经度、纬度)= meshgrid (X.lon X.lat);
[x, y,区]= ll2utm(纬度、经度);做这份工作!
z =双(X.z);z(z==-32768 | z<0) = NaN;
数字
pcolor (x, y, z);阴影平;抓住
轮廓(x, y, z,(0,0),“w”)
推迟;轴相等;轴紧
包含(“东(m)”);ylabel(北(m))
标题(sprintf('Sicily - UTM zone %d WGS84'，zone))
加载西西里岛SRTM全分辨率DEM (2400x4500网格)，将其转换为UTM，并用pcolor和等高线绘制结果。要制作一个常规的UTM网格，可以用griddata函数对x和y进行插值。

请参阅“doc ll2utm”和“doc utm2ll”获取语法和帮助。