UTM2LL将通用横切墨卡托(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”获取语法和帮助。
引用作为
弗朗索瓦Beauducel(2021)。LL2UTM和UTM2LL(//www.tianjin-qmedu.com/matlabcentral/fileexchange/45699-ll2utm-and-utm2ll), MATLAB中央文件交换。检索.