上一篇文章介绍了DEM数字高程数据的下载及matlab处理,这篇介绍制作固定大小的高程网格格式数据。
高程网格格式数据即为点云数据,后缀名为.xyz。在有些场合需要固定大小的高程网格格式数据,100100或者200200的,在地理空间数据云上下载的DEM数据往往数据量很大,不符合诗句需要,因此就要做裁剪,我用Global Mapper 19 制作了大小为101*101的高程网格格式的数据:
首先打开下载的地图,在地理空间数据云上下载90m分辨率的数据,我这里下载的是SRTMDEM,数据的压缩格式为.img,在Global Mapper中打开,如下图所示:
下面介绍格式转换步骤:
1:选择【文件】->【导出】->【导出高程网格格式】,即可以打开如下对话框,这里分辨率一般默认就好
选定导出边界选项卡,自行设定经纬度范围,由于分辨率为0.000833333333333333,如果我需要制作101101的高程网格数据,可以设定如下显示的经纬度
即可截取选定区域,生成的网格数据为101101
处理代码如下:
data=Loadmap(‘030.xyz’);
x=data(:,1);
y=data(:,2);
z=data(:,3);
xx=zeros(101,101);
yy=zeros(101,101);
zz=zeros(101,101);
for i=1:10201
zz(i)=z(i);
xx(i)=x(i);
yy(i)=y(i);
end
Mapx = xx(1,:);
Mapy = (yy(:,1))’;
Mapz=zz;
%%xx存放经度信息
%%yy存放纬度信息
%%zz即为需要的的网格矩阵,矩阵元素表示高度
显示结果如下图:
如有侵权,请联系作者