利用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]的图像为