这一段正要用到这个函数，欢迎大家共同学习。【数学建模&MATLAB|利用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）。
所有区域的命名由其左上角的坐标所决 定。
LatitudeLongitudeElevation
TileMinimumMaximumMinimumMaximumMinimumMaximumMeanStd.Dev.
-----------------------------------------------------------------------
W180N904090-180-14016098448482
W140N904090-140-10014635730596
W100N904090-100-6012416333280
W060N904090-60-20139401624933
也就 是由纬度的最大值和经度的最小值组成，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]的图像为
