地形分析精要.pptVIP

刘海新 河北工程大学资源学院 电话 ENVI/IDL 10-地形分析 ENVI/IDL 数字高程模型(Digital Elevation Model)，简称DEM，它是用一组有序数值阵列形式表示地面高程的一种实体地面模型。DEM 除了包括地面高程信息外，还可以派生地貌特性，包括坡度、坡向等，还可以计算地形特征参数，包括山峰、山脊、平原、位面、河道和沟谷等。 建立DEM 的方法有多种，从数据源及采集方式主要有：根据航空或航天影像，通过摄影测量途径获；野外测量或者从现有地形图上采集高程点或者等高线，后通过内插生成DEM 等方法。 地形分析-概述 ENVI/IDL 本节内容用到DEM EXtraction（DEM自动提取）扩展模块，该模块可以对从ALOS PRISM, ASTER, CARTOSAT-1, FORMOSAT-2, GeoEye-1, IKONOS,KOMPSAT-2, OrbView-3, QuickBird, WorldView-1, SPOT 1-5 以及航空影像立体像对中提取DEM。 DEM EXtraction模块除了自动提取DEM向导外，还有三个DEM工具：DEM编辑工具、立体3D量测工具和核线图像3D光标工具。 DEM编辑工具可以对DEM进行相应的编辑处理 立体3D量测工具可以从立体相对中量测一个点的高程信息 核线图像3D光标工具可以再3D立体视图环境中基于已存在的核线立体相对图像做3D量测。 DEM EXtraction模块菜单命令及其功能见课本P234表11.1 10.1立体像对DEM自动提取 ENVI/IDL 1、立体像对提取DEM向导 基于立体像对的DEM 提取流程图如下图所示，总体上分为六步： 10.1立体像对DEM自动提取 本次将以ASTER level 1级数据为数据源，ASTER level 1数据是以HDF-EOS格式提供，有两个文件，分别为.hdf和.met，包括了一个标准垂视（3N）和后视（3B）图像。 ENVI/IDL 第一步：输入立体像对 （1） 单击主菜单-File-OpenExternal File→EOS→ASTER，选择AST_L1A。hdf。 （2） 单击主菜单-Topographic-DEM Extraction- DEM Extraction Wizard-New，单击Select Stereo Image 按钮，选择ASTER VNIR band3N作为左影像（left image）， ASTER VNIR band3B为右影像（right image）；一般推荐垂直获得的图像为左视图，非垂直方向为右视图。 （3）系统会自动根据自带的星历参数文件获取图像区域的最大高程和最小高程，也可以根据已知信息手动录入。 （4）单击Next 按钮，进入第二步。 10.1立体像对DEM自动提取 ENVI/IDL 第二步：定义地面控制点 提供三种定义地面控制点方式： 不定义：无控制点，获得的DEM是相对高程 交互式定义：分别定为同名点，然后找到相应的x（E）和Y（N）以及高程，操作方法如几何校正。 读取控制点文件：从外部读取控制点。 这里选择不定义控制点，单击Next 按钮，进入第三步。 10.1立体像对DEM自动提取 ENVI/IDL 第三步：定义连接点 （1）提供三种定义连接点方式： 自动寻找：定义参数如下 连接点数目（Number of Tie Points）：40。 搜索窗口大小（Search Windows Size）：101。一定要比移动窗口大。范围在大于等于21的任意整数，参数值越大，找到匹配点的可能性就越高，但计算量也就增大。 移动窗口大小（Moving Window Size）：15。在搜索窗口中进行检查，寻找地形特征匹配的小区域，必须是奇数，最小值为5，使用较大的移动窗口将会获得更加可靠的匹配结果，但处理时间更长。移动窗口的大小与空间分辨率有关，一般情况是大于等于10m的图像，设置范围为9-15,5-10m分辨率的为11-21,1-5m分辨率的为15-41，小于1m分辨率的为21-81或者更高。 10.1立体像对DEM自动提取 ENVI/IDL 平均高程（Region Elevation）：自动从影像读取，根据RPC自动计算。 是否检查连接点（Examine and Edit Tie Points）：Yes 交互式手工定义 人工交互式定义连接点，至少需要9个连接点，定义方式是通过在左右两幅图中来选取同名点 外部读取控制点文件 读取外部连接点文件，文件的后缀为.pts 本例中选择选择自动寻找，点击下一步，进入查看/添加/编辑连接点步骤，单击Show Table，选择Sort Table By Error，误差大的点排在前面，逐个