基于ENVI5.1的“资源三号”立体像对DEM提取介绍.docVIP

基于ENVI5.1的“资源三号”立体像对DEM提取 ? 数字高程模型(Digital Elevation Model)，简称DEM，它是用一组有序数值阵列形式表示地面高程的一种实体地面模型。DEM除了包括地面高程信息外，还可以派生地貌特性，包括坡度、坡向等，还可以计算地形特征参数，包括山峰、山脊、平原、位面、河道和沟谷等。建立DEM的方法有多种，从数据源及采集方式主要有：根据航空或航天影像，通过摄影测量途径获；野外测量或者从现有地形图上采集高程点或者等高线，后通过内插生成DEM等方法。 ??资源三号搭载了四台光学相机，包括一台地面分辨率2.1m的正视全色TDI CCD相机、地面分辨率3.5m的前视和后视全色TDI CCD相机、一台地面分辨率5.8m的正视多光谱相机。其中前正后视全色相机，推扫成像形成三线阵立体像对，可用于DEM提取。经过测试发现，前视和正视，或后视和正视，可以组成立体像对进行DEM提取，效果较好，优于前视和后视提取的DEM。 2.详细操作步骤 ? DEM Extraction工具分为9个步骤，总体上我们可分为6个步骤，如下图所示。 图1.1 DEM提取流程 图2.1 Step 1 of 9输入立体像对 第一步：输入立体像对 （1）单击File-Open，以正视和后视图像为例。打开ZY3_01a_hsnnavp_*和ZY3_01a_hsnbavp_*文件夹内的tif文件。 （2）在Toolbox中，双击/ Terrain/DEM Extraction/DEM Extraction Wizard: New。 （3）单击Select Stereo Image按钮，选择正视图像（NAVP）为Left Image，后视图像（BAVP）为Right Image。 注：1、ENVI一般能自动识别RPC文件，否则需要手动选择RPC文件。 2、左影像选择垂直观测的影像，或者观测角度小的影像。可以简单通过对比立体像对两幅影像的地面分辨率，分辨率高的当作左影像。 （4）单击Next按钮。 第二步：定义地面控制点 提供三种定义地面控制点方式：不定义（No GCPs (relative DEM values only)）、交互式定义（Define GCPs Interactively）和读取控制点文件（Read GCPs From File）。由于缺少地面控制点数据，这里选择不定义控制点，单击Next按钮。注：不定义地面控制点得到的DEM是相对高程，即以卫星地势面为基准。 图2.2 Step2 of 9定义地面控制点 图2.3 Step 4 of 9自动选择Tie点参数设置 第三步：定义连接点 （1）提供三种定义连接点方式：自动寻找、交互式手工定义和外部读取控制点文件。选择选择自动寻找，参数设置如下所示， 连接点数目（Number of Tie Points）：1000。 搜索窗口大小（Search Windows Size）：481。 移动窗口大小（Moving Window Size）：101。 最小相关系数（Minimum Correlation）：0.90。 平均高程（Region Elevation）：自动从影像读取。 是否检查连接点（Examine and Edit Tie Points）：Yes。 注：几个参数的设置说明如下说明： 搜索窗口大小（Search Windows Size）：大于或等于21的任意整数，并且必须比移动窗口大。该参数的值越大，找到匹配点的可能性也越大，但同时也要耗费更多的计算时间。大致确定搜索窗口大小的方法：在立体像对（带有粗略地理坐标）两个图像上找到一个同名点，量测两个图像上同名点的距离D（像素单位），搜索窗口最小设置为（D+1）*2。 移动窗口大小（Moving Window Size）：在搜索窗口中进行检查，寻找地形特征匹配的小区域。移动窗口大小必须是奇数。最小的移动窗口大小是5，即为5×5像素。使用较大的移动窗口将会获得更加可靠的匹配结果，但也需要更多的处理时间。移动窗口的大小跟图像空间分布率有关系，根据如下所列设置： 大于等于10 米分辨率图像，设置值的范围是9~15。 5~10 米分辨率图像，设置值的范围是11~21。 1~5米分辨率图像，设置值的范围是15~41。 小于1米分辨率图像, 设置值的范围是21~81 或者更高。 （2）单击Next按钮，进入Step 5/9步骤。 （3）本例子中，自动找到414个Tie点，Maximum Y Parallax（Y视差）=2.3394，小于10即可进行下一步处理。从点的分布来看，有部分地区没有Tie点，下面手动交互添加一些Tie点。 注：