DEM数字高程数据之高程网格格式转换
上一篇文章介紹了DEM數(shù)字高程數(shù)據(jù)的下載及matlab處理,這篇介紹制作固定大小的高程網(wǎng)格格式數(shù)據(jù)。
高程網(wǎng)格格式數(shù)據(jù)即為點云數(shù)據(jù),后綴名為.xyz。在有些場合需要固定大小的高程網(wǎng)格格式數(shù)據(jù),100100或者200200的,在地理空間數(shù)據(jù)云上下載的DEM數(shù)據(jù)往往數(shù)據(jù)量很大,不符合詩句需要,因此就要做裁剪,我用Global Mapper 19 制作了大小為101*101的高程網(wǎng)格格式的數(shù)據(jù):
首先打開下載的地圖,在地理空間數(shù)據(jù)云上下載90m分辨率的數(shù)據(jù),我這里下載的是SRTMDEM,數(shù)據(jù)的壓縮格式為.img,在Global Mapper中打開,如下圖所示:
下面介紹格式轉換步驟:
1:選擇【文件】->【導出】->【導出高程網(wǎng)格格式】,即可以打開如下對話框,這里分辨率一般默認就好
選定導出邊界選項卡,自行設定經(jīng)緯度范圍,由于分辨率為0.000833333333333333,如果我需要制作101101的高程網(wǎng)格數(shù)據(jù),可以設定如下顯示的經(jīng)緯度
即可截取選定區(qū)域,生成的網(wǎng)格數(shù)據(jù)為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存放經(jīng)度信息
%%yy存放緯度信息
%%zz即為需要的的網(wǎng)格矩陣,矩陣元素表示高度
顯示結果如下圖:
如有侵權,請聯(lián)系作者
總結
以上是生活随笔為你收集整理的DEM数字高程数据之高程网格格式转换的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 去除word PDF水印 去除pdf限
- 下一篇: 使用GDAL读取SRTM格式高程数据