DEM建立与应用

[分享]DEM建立与应用

一、  目的

 

DEM是对地形地貌的一种离散的数字表达，是对地面特性进行空间描述的一种数字方法、途径，它的应用可遍及整个地学领域。通过对本次实习的学习，我们应：

1、              加深对DEM建立过程的原理、方法的认识；

2、              熟练掌握ARCVIEW中建立DEM、TIN的技术方法。

3、              结合实际、掌握应用DEM解决地学空间分析问题的能力。

二、  实验准备

 

1、  软件准备：Arcview

2、  数据准备：文件feapt-clip1.dbf，feapt-clip1.shp，feapt-clip1.shx，文件terlk-clip1.dbf，terlk-clip1.shp，terlk-clip1.shx，文件夹cal2和info。

三、  实验内容

 

1、  DEM及TIN的建立

(1)   由采样点数据建立表面

1)         在视图目录表中添加并激活采样点层面feapt-clip1.shp。

 

 

2)         从【Surface】菜单中选择【Interpolate Grid】命令。

 

 

3)         在出现的Output Grid Specification对话框中设定输出主题的范围、栅格单元大小及栅格行、列数。

 

 

图1.  由采样点生成的表面

 

 

4)         接下来出现的Interpolate Surface对话框中，从Method列表中选择Spline（注意：在菜单中只有IDW和Spline两种内插方法可以选择）。在Z Value Field列表中选择Elev（高程）字段，单击OK。

 

 

5)         生成新的栅格主题Surface from Soilsamp.shp（如图1）。

6)         将得到的结果层面与原始的DEM层面（Map Calculation1）相减，得到结果层面如图2所示：

图2.  Surface层面与原始DEM相减的结果及其统计值

(2)   由点、线数据生成TIN转为GRID

1)            添加并激活点层面feapt-clip1和线层面terlk-clip1。

2)            点击【Surface】菜单下的【Create TIN from features】；

3)           在“Create New TIN”对话框中定义每个主题的数据使用方式；    

 

 

在“Create New TIN”对话框中，对某一个主题需要指定高程源（Height Source），以何种表面特征输入（Input as），以及选某一个值的字段来作为属性信息（可以为None）。（见图3）

 

 

图3.  生成TIN过程中的Create New TIN对话框

 

 

4)           确定生成文件的名称及其路径，生成新的层面tin-point。（见图4）

 

 

图4.  生成的TIN

 

 

5)           点击【Theme】菜单下的【Convert to grid】，确定生成文件的名称及其路径，生成新的Grid层面。（见图5）

 

 

图5.  由TIN生成的DEM

 

 

6)           将新生成的Grid层面与原始的DEM层面（Map Calculation1）相减，得到结果层面如图6所示：

 

 

图6.  新生成的DEM层面与原始DEM相减的结果及其统计

 

 

2、  DEM的应用

 

 

(1)         地形指标的提取

 

 

1)             坡度

 

 

具体的方法步骤如下:

 

 

I、               添加Dem数据并激活它。

 

 

II、            从【Surface】菜单中选择【Derive Slope】命令。

 

 

III、       生成新的坡度主题slope of Dem。

 

 

IV、          双击左边的图例，在弹出的Legend Editor对话框中可重新调整坡度分级（如图7）。

 

 

图7.  提取坡度

 

 

2)             坡向

 

 

具体的方法步骤如下:

 

 

I、               在视图目录表中添加dem并激活它。

 

 

II、            从【Surface】菜单中选择【Derive Aspect】命令。

 

 

III、       显示并激活生成的坡向主题Aspect of Dem（如图8）。

 

 

图8.  提取坡向

 

 

3)             平面曲率

 

 

I、               激活坡向数据。

 

 

II、            从【Surface】菜单中选择【Derive Slope】命令。

 

 

III、       生成平面曲率层面Slope of Aspect。（见图9）

 

 

图9.  提取平面曲率

 

 

4)             剖面曲率

 

 

具体的方法步骤如下:

 

 

I、               激活坡度数据。

 

 

II、            从【Surface】菜单中选择【Derive Slope】命令。

 

 

III、       显示并激活生成的剖面曲率层面Slope of Slope（如图10）。

 

 

图10.  提取剖面曲率

 

 

(2)         提取等高线

 

 

具体的方法步骤如下:

 

 

I、               在视图目录表中添加dem并激活它。

 

 

II、            从【Surface】菜单中选择【Create Contours】命令。

 

 

III、          在出现的Contours Parameters对话框中输入等高距Contour interval和基础等高线的值Base Contours。

 

 

IV、          生成等高线主题Contours of Dem（如图11）。

 

 

图11.  提取等高线

 

 

(3)         地形表面的阴影图

 

 

具体的方法步骤如下:

 

 

I、               在视图目录表中添加dem并激活它。

 

 

II、            从【Surface】菜单中选择【Compute Hillshade】命令。

 

 

III、          在Compute Hillshade对话框（如图12）中，输入计算Hillshade的参数值。

 

 

IV、          

图12  Compute Hillshade对话框

生成地表阴影主题Hillshade of Dem（如图13）。

 

 

(4)        

图13.  提取地表阴影图

可视性分析

 

 

1)             通视性分析

 

 

在Arcview中，进行通视性分析有两种具体操作

 

 

第一种操作：

 

 

例如：分析某区域内S与P两点间的通视情况。

 

 

I、             添加Dem或TIN主题作为通视性分析的地形表面并激活它。

 

 

II、          从工具栏选择Line of sight工具 。

 

 

III、        

图14  Line of sight对话框

在出现的Line of Sight对话框中输入观察者Observer与目标物Target距地面的距离（如图14），单击OK。

 

 

IV、        按住鼠标左键，屏幕上将会出现十字光标，将光标从观察点S移向目标点P，然后释放光标。在观察点到目标点之间将会出现一条视线，其中可视的部分为红色，不可视的部分为绿色。并且，在Arcview窗口底部的状态栏显示了从观察点到目标点是否可视。若不可视，在视线上将会用圆点表示第一个障碍物的位置，它的xy坐标将会在状态栏显示（如图15）。

 

 

图15  S、P两点间通视情况示意

S

P

 

 

第二种操作：

 

 

I、             从【File】菜单中选择【Extensions】命令。

 

 

II、          在Extensions的对话框中选择Visibility Tools选项。

 

 

III、        在Arcview工具条中出现Line of sight工具 。

 

 

IV、        选择此工具，在Available Grids对话框中选择栅格主题，例如：Dem作为通视性分析的地形表面（如图16）。

 

 

图16  Available Grids对话框

  

 

 

图16  Available Grids对话框          

 

 

V、           在出现的Set Visibility Parameters对话框中输入观察者与目标物距地面的距离。单击OK。

 

 

VI、          按住鼠标左键，将光标从观察点A移向目标点A′，然后释放光标。在观察点到目标点之间将会出现一条视线，其中可视的部分为红色，不可视的部分为绿色（如图17）。同时，Arcview

会自动绘出A—A′两点间的通视剖面图（如图18）。

 

 

A

A′

A′

A

    

 

 

 

 

 

图 17  A -A′间的通视情况示意                           图 18   A —A′两点间的通视剖面图

 

 

                                                                         （红色实线条为可视区，绿色虚线为不可视区）

 

 

 

2)             可视区分析

 

 

具体操作如下：

 

 

I、             在视图目录表中添加Dem作为可视区分析的地形表面。

 

 

II、          创建或添加一个点主题point.shp包含观测点。

 

 

III、        同时激活Dem和point.shp主题。

 

 

IV、        从【Surface】菜单中选择【Calculate Viewshed】命令。

 

 

V、           生成可视区栅格主题visibility of point（如图19）。

 

 

图19  可视区分析

(5)         水文分析

 

 

1)   洼地填充

 

 

操作步骤如下：

 

 

I、               激活视图目录表中的Dem。

 

 

II、            从【Hydro】菜单中选择【Fill sinks】命令。

 

 

III、       在视图目录表中显示生成的Filled Dem，既为填充过洼地的Dem数据。

 

 

3)             水流方向计算

 

 

在Arcview中提取水流方向的步骤如下：

 

 

I、               从视图目录表中激活Filled Dem。

 

 

II、            从【Hydro】菜单中选择【Flow Direction】命令。

 

 

III、       显示新生成的水流流向数据Flow Direction（见图20）。

 

 

4)            

图20  提取水流方向

流水累积量

 

 

计算流水累积量的步骤如下：

 

 

I、               从视图目录表中激活Flow Direction。

 

 

II、            从【Hydro】菜单中选择【Flow Accumulation】命令。

 

 

III、       显示新生成的流水累积量数据主题Flow Accumulation。双击此主题的图例，在出现的Legend Editor对话框（见图21）中，

 

 

 

 

图21  Legend Editor对话框

 

 

IV、          单击Classify按钮，在出现的Classification对话框中选择标准偏差Standard Deviation为分类类型，单击OK。

 

 

提取的Flow Accumulation主题（见图22）显示了地面水系的分布状况。

 

 

图22  提取水流累积量

 

 

4)   水网

 

 

提取水网密度的步骤如下：

 

 

I、               从视图目录表中激活Flow Direction。

 

 

II、            从【Hydro】菜单中选择【Stream Network as Line Shape】命令。

 

 

III、          在出现的Stream Network对话框中输入计算水网密度的最小河流长度，其单位为栅格数。

 

 

IV、          在Hydro.StreamNetwork对话框中选择水流方向Flow Direction主题，单击OK。

 

 

V、             生成新的主题Stream Network Shape即为水网密度分布图。

 

 

不同的河流最短长度值所计算的水网密度如下图（见图23、图24）所示：

 

 

图23  最短河流长度为20提取的水网密度

    

图24  河流最小长度为 100米 提取的水网密度

 

 

5)   水流长度计算

 

 

提取水流长度的步骤如下：

 

 

I、               从视图目录表中激活Flow Direction。

 

 

II、            从【Hydro】菜单中选择【Flow Length】命令。

 

 

显示输出的Flow Length栅格主题（见图25）。

 

 

图25  水流长度提取

 

 

6)   流域分析

 

 

采用水流方向和流水累积量数据生成Watershed数据的步骤如下：

 

 

I、               在视图目录中添加Flow Direction和Flow Accumulation。

 

 

II、            击活Flow Accumulation。

 

 

III、          从【Hydro】菜单中选择【Properties】命令。

 

 

IV、          在出现的Hydro Properties对话框（见图26）中输入计算所需的Flow Direction和Flow Accumulation主题。

 

 

图26  Hydro Properties对话框

 

 

V、             从【Hydro】菜单中选择【Watershed】命令。

 

 

VI、          在出现的Watershed对话框中输入计算流域的面积（最少栅格数）。

 

 

VII、        生成Watersheds主题（见图27）。

 

 

图27  提取Watershed

 

 

四、  实习报告要求

 

将所做工作以幻灯片形式做以汇报，内容包括原理、过程、结果。