舟曲泥石流受灾面积计算模型..docVIP

舟曲泥石流受灾面积计算模型 第七组 王政 白燚焱 邓敏 摘要 本文研究的是舟曲泥石流受灾面积计算问题。主要分两大步走：第一步，将舟曲进行区域划分；第二步，计算各个区域的面积。论文的关键在于将灾前和灾后图像进行配准。在MATLAB软件的支持下，将灾前和灾后图像进行通过实用点映射法和归一化灰度组合法配准图像后我们首先将灾前和配准后的灾后灰度图像转换为矩阵；然后将矩阵相减得到新的矩阵，矩阵中小于等于0的部分为不受灾区，其余部分为受灾区；再将此结果与原舟曲图像进行比对，将农田、城区、村庄、山体和其它划分出来；最后利用图像任意两点与Googleearth中对应的两点多次比对，得出每个点在实际位置中代表的大小为：S=36 m/点。再利用配准后的矩阵中点的个数计算出各个受灾区域的面积和总体面积。并且结合实际情况形成灾情报告。 在模型一中：由于灾前和灾后的灰度图像拍摄过程中存在天气、照度等差异，所以我们需要将灾前和灾后图像进行配准。在此模型中，我们利用实用点映射的图像配准法和归一化灰度组合法将灾前和灾后图像进行配准得到配准后的灾后图形如图三。两种方法的结合是本文的特色之处。 为了提高分析的精确度，在MATLAB软件的支持下，我们将舟曲灾前、配准后的灾后灰度图像转化为了矩阵。 在模型二中，我们首先在模型一矩阵的基础上，将矩阵三与矩阵四相减得到新的矩阵；然后分析矩阵，大于0的部分则为泥石流受灾区，其余地区不受泥石流灾害的影响；最后将泥石流受灾区与舟曲灾前图像进行对比，划分出农田、城区、山体、村庄和其它。 在模型三中，我们在舟曲图像中多次随机取两个点与Googleearth中对应点进行比对，并计算出：S=36 m/点。利用这个结果与模型二中配准后的灾后灰度图像转化的矩阵结合，求出受灾不同区域的面积和总面积。受泥石流灾害的总面积达到＝2502684。 综合前面的模型，结合实际情况，对灾区的灾情进行总结，形成报告。 关键词 区域的划分、区域面积的计算、MATLAB软件、图像转化为矩阵 、实用点映射的图像配准法、归一化灰度组合法 一、问题的重述 8月7日晚至8日凌晨，甘肃甘南州舟曲县突发特大泥石流灾害，造成重大人员伤亡。测绘部门紧急行动，服务抢先救灾。先遣队的主要任务是，抵达现场选建无人机起降场地，利用无人机对堰塞湖、河流、道路、居民区等受灾地区进行航空摄影，获取救援急需的影像资料，并现场处理、制作高分辨率影像图，为省、洲、县抗洪救灾指挥部提供决策用图。 中心科研人员利用上述两期数据对舟曲县泥石流灾害进行信息提取和对比分析，此次灾害造成的损失如下： （1）泥石流灾害i：舟曲县城主要遭受由北面罗家峪和三眼峪两条大型泥石流的冲击，造成农田及城区建筑物的严重损毁，航空遥感数据解译泥石流淹没面积的达0.44平方公里，泥石流淹没范围如附件一图中红色显示。 （2）形成0.13平方公里堰塞湖：此次泥石流灾害造成县城白龙江下游堵塞，形成堰塞湖，造成县城部分区域淹没，有航空数据解译可知洪水淹没面积共计0.13平方公里，洪水淹没范围如附件二中蓝色所示。 资料：http://www,svsn,giv,cn/article/chyw/201008/20100800071112.shtml http://www,svsn,giv,cn/article/chyw/201008/20100800071159.shtml http://www,meizhou,cn/news/chyw/1008/10/10081000194.html 通过对获取的舟曲的某一具体区域的地形图和遥感图形，分析如下问题： （1）根据要干地形图对泥石流受灾区域进行划分； （2）对受灾区域内不同类别区域（如住宅区、农田、山体等）的受灾面积的甲酸及总受灾面积计算。 二、模型的假设 1、在泥石流发生的时候不伴随其他自然灾害，如：地震； 2、在拍摄灾后图像时泥石流巨大冲击已经结束，整个地区的泥石流灾害已经趋于稳定； 3、图像的采集及时没有耽误； 4、在图像处理过程中没有造成信息的流失； 5、忽略堰塞湖所造成的灾害的面积； 6、忽略后续的泥石流，即模型以采集照片时