Matlab使用etopo在线地形数据绘制中国区域DEM地形图

以下是使用MATLAB绘制中国区域DEM地形图的过程和代码示例：

1. 首先，需要从etopo网站下载中国区域的地形数据。进入etopo网站（https://www.ngdc.noaa.gov/mgg/global/etopo5.HTML），找到“Download Global Relief Data”一栏，在“Area”一栏中选择“Custom”, 然后在“Bounding Box”一栏中输入中国区域的经纬度范围（例如：东经73°至135°，北纬18°至54°），在“Data Format”一栏中选择“GMT Grid Format”，最后点击“Submit”按钮进行下载。

2. 下载完成后，将下载的文件解压缩，得到两个文件：一个是地形数据的二进制文件（.grd），另一个是元数据文件（.gmt）。

3. 在MATLAB中读取地形数据文件。可以使用MATLAB自带的函数`readgeoraster`，代码如下：

   [Z, R] = readgeoraster('filename.grd');
 

   其中，`filename.grd`是地形数据文件的路径和文件名。`Z`是DEM数据矩阵，`R`是参考对象，包含地图的元数据。

4. 将DEM数据绘制成地形图。可以使用MATLAB自带的函数`geoshow`，代码如下：

   figure
   geoshow(Z, R, 'DisplayType', 'surface');
   demcmap(Z);
   colorbar
 

   其中，`Z`和`R`是从地形数据文件中读取的DEM数据矩阵和参考对象。`'DisplayType', 'surface'`表示以表面的形式显示DEM数据。`demcmap(Z)`表示使用DEM颜色映射来显示DEM数据。`colorbar`用于添加颜色条。

完整代码示例：

% 读取DEM数据
[Z, R] = readgeoraster('filename.grd');

% 绘制地形图
figure
geoshow(Z, R, 'DisplayType', 'surface');
demcmap(Z);
colorbar