分为两部分:1是下载数据,2是命令显示
1 下载数据
进入http://www.mathworks.com/support/tech-notes/2100/2101.html,在其中找到如下位置
点击上图中右下角的[1],进入如下界面
注意上面这个地图中有很多网格,其中一个格子就表示一个区的数据集,我选的美国西南角加州所在格,出现如下界面
这个地区所对应的数据文件为w140n40,就是上图中我用红笔标出来的位置,点击就可以下载了。下载后的文件是一个压缩包w140n40.tar,解压后有一个W140N40.DEM文件就是地形数据了,记下这个文件所在的目录。
2. 命令显示
现在有了数据就可以画图显示了,下面命令行中注意w140n40数据的文件名,不需要后缀.DEM了。
clear
latlim = [ 33.10 34.7]; %这二行是设置经纬度的范围,记得要在上面格子的范围内,至于跨格的,我还没试过
lonlim = [-118.0 -117.0];
gtopo30s(latlim, lonlim); %截取指定经纬度中的数据
[Z,refvec] = gtopo30('C:\w140n40\W140N40',1,latlim,lonlim); %用gtopo30函数读取30秒的数据
worldmap(Z,refvec); %产生空白底图
geoshow(Z,refvec,'DisplayType','surface') % 画图