融合机载激光扫描数据的航空影像水体特征提取培训课件

融合机载激光扫描数据的航空影像水体特征提取研究目的影像水体提取在环境遥感、环境监测等领域应用非常广泛，尤其在环境保护、洪水分析、水资源调查中，影像的水体提取技术为这些应用提供了强有力的保障。目前的影像水体特征提取主要分为：目视人工解译----适用于各类影像基于光谱信息的解译----适用于多光谱、高光谱影像单波段法多波段法基于图像纹理信息的解译----适用于各类影像航空影像的水体特征提取是当前遥感影像水体特征提取的薄弱环节，提取的效率和精度都有待提高。2融合机载激光扫描数据的航空影像水体特征提取第一部分:机载激光扫描（LightDetectionAndRanging，LiDAR）3融合机载激光扫描数据的航空影像水体特征提取机载激光扫描原理机载激光扫描是在飞机上搭载激光扫描设备，由GPS提供的飞机位置数据、INS惯性导航系统提供的飞机姿态数据以及激光发射器获取激光测距数据获得一个大范围的带状区域内的地物点云数据。激光扫描的设备组成部分和激光测距原理4融合机载激光扫描数据的航空影像水体特征提取机载激光扫描的外、内业处理过程外业采集的数据:

飞机上安装的GPS实时观测数据地面基站的GPS实时观测数据飞机的INS系统数据激光测距数据航空影像数据处理过程:

数据分离差分GPS处理差分GPS和IMU数据集成集成数据的内插和平滑点云数据生成航空影像内、外定向最终数据激光扫描点云数据飞机的实时定位、定姿数据航空影像5融合机载激光扫描数据的航空影像水体特征提取点云（PointsCloud）点云是以一个整体的状态存在，能够很清晰地看到其表面的各种细节形态的大量三维点数据。点云数据在逆向工程、三维重建中具有广泛的用途。点云数据一般以三维坐标X,Y,Z的方式给出，在一些仪器采集的数据中，还可能包含回光强度值信息和激光点的颜色信息。在同济大学校园里采集的点云数据6融合机载激光扫描数据的航空影像水体特征提取研究过程和步骤方法：人工目视判读方法获取同名点进行匹配输出：影像向激光扫描数据转换的坐标参数点云与影像匹配机载激光扫描数据分析狄洛尼三角网和水体三角网水体三角网栅格化、缓冲航空影像与缓冲栅格的运算水体边缘的精化三角网的直方图统计分析带坐标框架的航空影像均值偏移向量法(MeanShift)影像分类8融合机载激光扫描数据的航空影像水体特征提取试验数据所在位置实验数据来自美国俄亥俄州北部，数据覆盖范围约为2.3平方公里。大致位于北纬41度44分，西经81度17分http://metadataexplorer.gis.state.oh.us/metadataexplorer/explorer.jsp

9融合机载激光扫描数据的航空影像水体特征提取实验数据机载激光扫描点云数据：共439，996点，点的平均密度约为0.02/平方英尺。由Leica的ALS系统采集。航空影像数据：真彩色航空影像，5000*5000分辨率，像素分辨率为约为1英尺。由普通CCD相机采集。10融合机载激光扫描数据的航空影像水体特征提取点云数据与航空影像的匹配将激光点云数据进行按高程和回光强度值内插。采用人工判读的方法，按房屋和道路对象提取同名点。各提取了14个同名点，并进行6参数的迭代平差坐标转换，转换结果如下：当同名点选房屋时，转换参数为：2299963.64240717765001.8608268831.007301949161960.00765490911180677-4.4084076080253e-005-0.999834732412637X方向的残差中误差为2.35英尺，Y方向的残差中误差为1.97英尺当同名点选道路时，转换参数为：2300000.52919327765005.0717043550.9989745181146650.00114835895237775-0.000278379912288074-1.00039224778675X方向的残差中误差为2.77英尺，Y方向的残差中误差为3.09英尺结论：采用同名点判读的方法可以实现影像向点云坐标系统的匹配；房屋同名点的匹配精度略优于道路同名点；X方向和Y方向的匹配中误差约为2-3英尺。对于影像的水体提取，该匹配精度已经可以完成相关的工作。11融合机载激光扫描数据的航空影像水体特征提取融合激光扫描数据的航空影像水体特征提取原理当前主流的LiDAR系统采用的波段均为近红外的波段，如：Leica公司的ALS系统和ALTM系统采用1.064微米波段IGI公司的LiteMapper系统采用1.550微米波段此外，还有1.047微米也常用于机载LiDAR系统测深LiDAR系统还采用透水性较好的蓝绿激光波段，如532纳米常见地物的在近红外波段的光谱特性水体土壤植被/Geography/DeLiberty/Geog474/spectral_signatures.gif12融合机载激光扫描数据的航空影像水体特征提取各类地物的狄洛尼三角网从激光扫描点云数据和航空影像对比可以发现，航空影像的水体所在位置，在激光数据中出现了大量的缺失；而在其他区域，水体的数据分布基本均匀。为了定量的获取激光点云数据中的数据缺失部分，采用基于四方边缘结构的(Quad-Edge)逐点插入算法，对激光点云数据构建狄洛尼三角网。数据缺失区域的三角网形状将有别于数据正常区域。草地、道路、水体、建筑物、树木、山体的三角网比较13融合机载激光扫描数据的航空影像水体特征提取水体三角网通过各类地物的三角网比较可以发现，由于水体部分缺少激光点，其三角网形状与周围存在较大差异。其他地物的三角网在平面位置和大小形状方面差异不大。因此，采用平面面积比较的方法对三角网进行统计分析。最小面积0.003最大面积155410平均面积28.018从面积统计表可以发现，三角网的最大最小面积差异很大。为了更好的分析面积的分布情况，取面积小于100的三角形进行直方图统计。并以2倍的平均面积作为阈值进行水体三角网的提取。14融合机载激光扫描数据的航空影像水体特征提取水体三角网的栅格化、缓冲水体三角网是矢量结构的数据，无法直接用于GIS的空间分析，因此，在ArcGIS软件中，将矢量格式的三角网数据通过栅格化的方式转换为栅格数据。缓冲分析是将数据按照一定的距离产生缓冲区。我们对栅格数据采用三角网平均边长（8.8英尺）的缓冲分析。缓冲分析的主要目的是：消除匹配误差带来的影像，在影像和点云匹配过程中，控制点的残差中误差为2-3英尺，8.8英尺也相当于3倍的控制点中误差。能消除匹配误差带来的影像。缓冲分析过程15融合机载激光扫描数据的航空影像水体特征提取影像与栅格数据的栅格运算1001在栅格运算前，先将栅格数据进行重分类，水体三角网所在的栅格点赋值为1，其他点赋值为0。栅格运算类似于矩阵的点运算，即对应像素进行相应的运算。栅格乘法运算25512233128*16融合机载激光扫描数据的航空影像水体特征提取粗水体的精化粗水体：由于激光点反应的位置不一定就是水体的边缘，此外，在前面的处理中，我们还对水体三角网进行了缓冲，因此，得到的水体边缘还不能精确的反应影像中的水体边缘特征。改进的MeanShift算法用于影像分类传统的MeanShift分类算法是针对二维或三维的信息的，在粗水体信息中，存在像素的位置信息X、Y以及像素的属性光谱信息R、G、B，因此，需要对MeanShift算法的距离函数进行改进。MeanShift算法的阈值也是值得探讨的问题。在实验中，我们截取了其中部分数据作为试验，分别选取了0.04,0.05,0.06,0.07,0.08,0.09,0.10,0.15,0.20和0.30作为阈值进行处理，对比得到的结果发现阈值是0.06时，提取效果最佳。因此，对所有的粗水体数据，我们采用0.06作为阈值进行水体的精化。17融合机载激光扫描数据的航空影像水体特征提取精水体特征18融合机载激光扫描数据的航空影像水体特征提取结果比较本文方法人工判读解译ERDAS监督分类19融合机载激光扫描数据的航空影像水体特征提取结论人工判读、监督分类以及本方法都可以实现对水体特征的提取。采用人工判读的提取可以充分考虑操作员对地域的熟悉程度，结合操作员的经验，从全局来看提取结果较为准确。对于细小的水体，人工判读的方法也可以进行有效提取。然而人工判读的结果在边缘部分提取精度不高。采用监督分类的提取只考虑样本的光谱信息，因此同谱异质以及同质异普现象极易对提取结果产生影响。本方法与人工判读的结果相比，提取的水体边缘精度高，而在细小水体上无法实现准确的分类。然而，采用本文的方法可以自动实现水体的提取，减少人工操作带来的大量工作和个体差异。与监督分类的结果相比，本方法既考虑了光谱信息，也充分融合了LIDAR的空间信息以及激光的光谱特性，提取精度明显优于监督分类。在边缘部分，由于狄洛尼三角网边缘是凸多边形的性质，导致提取结果在影像边缘部分与实际不符。而人工分类、监督分类在边缘部分不存在这种情况。20融合机载激光扫描数据的航空影像水体特征提取后续的研究狄洛尼三角网的凸性质，使得提