引用DEM的矩阵

Question

我已按照Tile E的步骤，从NOAA / NGDC-GLOBE DEM（tileinfo＆amp; metadata）读到this blog。

我的目标略有不同，即找到特定坐标的W，SW和NW的最高点。所以在我得到该区域的矩阵dem之后，我创建了一个引用矩阵：

Rdem=makerefmat('RasterSize',[6000 10800],'Latlim',[0 50],'Lonlim',[-180 -90]);

当我使用ltln2val(dem,Rdem,40,-105)时，我得到的值为-500，表示海洋。这个位置在科罗拉多州。有人可以解释我哪里出错吗？

由于 多米尼克

Answer 1

makerefmat假设dem矩阵从底部开始而不是顶部，这对于许多映射函数是常见的。但是因为Tile是使用泛型函数（而不是映射工具箱）读入的，所以它就像普通矩阵一样。 makerefmat可以选择切换为地理配准索引行的方向。 Rdem的正确语法是：

Rdem=makerefmat('RasterSize',[6000 10800],'Latlim',[0 50],'Lonlim',[-180 -90],'ColumnsStartFrom','north');