% 读取DEM高程数据
demFile = '你的DEM文件.tif';  % 替换为你的DEM文件路径
[demData, R] = geotiffread(demFile);

% 输入的经纬度
lat = 40.7128; % 纬度，例如纽约
lon = -74.0060; % 经度，例如纽约

% 将经纬度转换为像素坐标
[row, col] = latlon2pix(R, lat, lon);

% 确保行列号在数据范围内
if row >= 1 && row <= size(demData, 1) && col >= 1 && col <= size(demData, 2)
    % 获取对应的高度值
    height = demData(round(row), round(col));
    disp(['该位置的高度为: ', num2str(height), ' 米']);
else
    disp('经纬度超出DEM数据范围');
end
  • 1.
  • 2.
  • 3.
  • 4.
  • 5.
  • 6.
  • 7.
  • 8.
  • 9.
  • 10.
  • 11.
  • 12.
  • 13.
  • 14.
  • 15.
  • 16.
  • 17.
  • 18.
  • 19.