地形图生成DEM

文章目录

• 前言
• 一、地形刨面线平面图与坡度平面图
• 二、具体实现

• 1.利用地形刨面数据和matlab绘制刨面图
• 1.计算地形刨面图的坡度值并绘制坡度平面图

• 总结

前言

1.如何利用地形刨面数据和matlab进行绘制平面图？

2.如何计算地形刨面图的坡度值并绘制坡度平面图？

3.地形刨面图与坡度图有何联系？

注：（基于个人爱好，以下文章纯属个人撰写，但由于技术拙略，如有错误之处，望各位读者指正谢谢！）

一、地形刨面线平面图与坡度平面图

• 地形剖面图指沿地表某一直线方向上的垂直剖面图，以显示剖面线上断面地势起伏状况。
• 坡度平面图是采用计算刨面图相邻点之间的坡度值并绘制的二维平面图

二、具体实现

1.利用地形刨面数据和matlab绘制刨面图

某地区的刨面数据获取和txt数据源（头文件去除）

代码如下（示例）：

%利用matlab绘制抛面线以及坡度图
%刨面线绘制
A=importdata('抛面线实验数据.txt');%此处的数据可以自行去下载
xyzdata=A.data;
length_data=length(xyzdata)
max_data=max(xyzdata(:,3))
X=xyzdata(:,1);
Y=xyzdata(:,2);
Z=xyzdata(:,3);
scatter3(X,Y,Z,5,"black","filled","d")
pl=line(X,Y,Z);
pl.Color='red';
%二维平面图
x=linspace(0,max_data,length_data);
plot(x,Z)
scatter(x,Z,5,"red","filled",'d')

结果显示：

1：该地区立体三维刨面图

2：该地区二维平面刨面图

3：散点图(便于查看某些点的位置和属性)

1.计算地形刨面图的坡度值并绘制坡度平面图

代码如下（示例）：

%计算坡度
length_gd1=diff(Z);
pd=atan(length_gd1/20);%此处的20属于变动值，及点间距需结合数据手动输入
x1=linspace(1,max_data,length_data-1);
plot(x1,pd)
scatter(x1,pd,5,"red","filled",'d')

结果显示：

1.坡度平面图

2.坡度散点图

总结

• 我认为，上述操作主要有以下两个浅方面的小作用：

（1）地形形变监测：
针对一个地形发生微小形变的区域，当我们的刨面图不足以观察到细微变化时，可以通过坡度图进行观察，对数据进行拉伸和放大处理，就可以实现该区域的形变监测。

（2）选址规划：
当我们需要在一个区域选取地点修建梯田或者其他建筑时，往往要考虑到坡度的大小。因此，坡度图可以用来规划该地区的一个适应项目的选址方案。

• 另外，我还需要考虑以下几个问题。希望读者能够不吝赐教！
• 1.上述程序虽然解决了刨面平面图和坡度图的绘制，但是目前刨面图和坡度图联合分析却是少之又少，绘制了坡度图到底能够实现什么社会实际问题？
• 1.若是将刨面平面图换成地形图，又该如何绘制出其坡度图，如何用matlab进行实现？