利用MATLAB绘制GTOPO30的DEM地形图

这一段正要用到这个函数，欢迎大家共同学习。

GTOPO30 是一个全球的数字高程模型，它覆盖西经180度至东经180度，南纬90至北纬90度的所有区域。它的分辨率为30秒（即0.00833333度），生成 一个21,600（行）×43,200（列）的DEM。其高程值范围在-407至8752米。在数字高程模型中，海洋地区被指定为-9999 。 低地沿海一带，海拔至少有1米。由于图象为栅格数据，在影像上小于一个像素的点被忽略（即面积小于1平方千米的小岛）。

为了便于数据的分 发，GTOPO30 被划分为33个小的区域（tiles）。从南纬60度至北纬90度，西经180度至东经180度的地区被划分为大小为50（纬度）×40（经度）的27块 区域，南极洲（南纬90至南纬60度，西经180至东经180度的区域）被划分为6个区域（30×60）。

所有区域的命名由其左上角的坐标所决 定。

Latitude Longitude Elevation

Tile Minimum Maximum Minimum Maximum Minimum Maximum Mean Std.Dev.

------- ---------------- ---------------- --------------------------------

W180N90 40 90 -180 -140 1 6098 448 482

W140N90 40 90 -140 -100 1 4635 730 596

W100N90 40 90 -100 -60 1 2416 333 280

W060N90 40 90 -60 -20 1 3940 1624 933

也就 是由纬度的最大值和经度的最小值组成，EW,NS由符号决定。

命名的图示如下：



该数据可以在ftp://edclxs71.cr.usgs.gov/pub/data/gtopo30/global中下载，为USGS匿名服务器。

一、解压缩

将数据解压出来。一般来说，每个包都含有8个文件，后缀分别是DEM, DMW, GIF, HDR, PRJ, SCH, SRC, STX，绘图时是需要DEM文件，但建议其它格式文件也都保留，另外的一些绘图工具要用到。

二、绘图

MATLB继承的gtopo30函数可以直接将图像绘制，现绘制W100N40的区域的全图

[datagrid refvec] = gtopo30('W100N40', 2, [24 31], [-90 -78]); %gather topological data for 1' resolution
datagridt = flipud(datagrid);
image(datagridt);

flipud函数是用作翻转

使用image函数可以直接可视化图像，则区域[24 31], [-90 -78]的图像为