专利名称：从高空间分辨率光学影像中提取倒塌房屋的方法及装置的制作方法
技术领域：
本发明涉及模式识别领域，尤其涉及一种从高空间分辨率光学影像中提取倒塌房屋的方法及装置。
背景技术：
高分辨率光学影像(高分影像)通常是指空间分辨率小于或等于I米的可见光彩色三波段或者全色单波段光学影像。因受到外力影响出现结构性粉碎性损坏而倒塌的房屋，在高分辨率光学影像中会呈现出明显的视觉特征。灾后救援是当前应对地震灾害的主要手段。灾后倒塌房屋的数量及其分布信息对灾后科学救援起着关键作用。近年来，国内外利用灾后单时相高分影像进行倒塌房屋快速提取已有初步尝试。该方法不受灾前数据缺失的制约，同时，能够很大程度减轻人机交互工作量。该方法的核心是对灾后高分影像中倒塌房屋与其它地物在高频区域的空间分布差异进行建模，但现有方法很难兼顾计算效率和计算精度。
发明内容
有鉴于此，本发明实施例提供了一种从高空间分辨率光学影像中提取倒塌房屋的方法及装置，目的在于解决现有的从高分影像中提取倒塌房屋的方法不能兼顾计算效率与计算精度的问题。为了实现上述目的，本发明实施例提供了以下技术方案:一种从高分辨率光学影像中提取倒塌房屋的方法，包括:使用结构元素，提取待处理光学影像中的方向形态纹理信息；依据所述方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像；利用预设的第一阈值分割所述方向不变纹理特征影像，得到第一分割影像，利用预设的第二阈值分割所述方向不变纹理特征影像，得到第二分割影像；将所述第一分割影像及所述第二分割影像进行基于数学形态学的重建，从而获取倒塌房屋图像。优选地，所述使用结构元素，提取待处理光学影像中的方向形态纹理信息包括:分别使用预设数量不同方向的、最小单元的线状结构元素对待处理光学影像进行高帽变换，从而得到所述预设数量的不同方向的形态纹理信息影像。优选地，所述依据所述方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像包括:对所述方向形态纹理信息影像进行平滑滤波；依据滤波后的方向形态纹理信息影像，选取满足预设条件的方向形态纹理信息影像；将所述满足预设条件的方向纹理信息影像确定为所述待处理光学影像的方向不变形态纹理特征影像。优选地，所述将所述第一分割影像及所述第二分割影像进行基于数学形态学的重建包括:将所述第一分割影像作为种子，将所述第二分割影像作为掩膜，进行基于数学形态学的重建计算，所述第一阈值大于所述第二阈值。优选地，所述方法还包括:依据所述待处理光学影像的类型，将所述待处理光学影像与所述倒塌房屋图像进行融合可视化。一种从高分辨率光学影像中提取倒塌房屋的装置，包括:提取模块，用于使用结构元素，提取待处理光学影像中的方向形态纹理信息；计算模块，用于依据所述方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像；分割模块，用于利用预设的第一阈值分割所述方向不变纹理特征影像，得到第一分割影像，利用预设的第二阈值分割所述方向不变纹理特征影像，得到第二分割影像；重建模块，用于将所述第一分割影像及所述第二分割影像进行基于数学形态学的重建，从而获取倒塌房屋图像。优选地，所述提取模块包括:
高帽变换单元，用于分别使用预设数量不同方向的、最小单元的线状结构元素对待处理光学影像进行高帽变换，从而得到所述预设数量的不同方向的形态纹理信息影像。优选地，所述计算模块包括:滤波单元，用于对所述方向形态纹理信息影像进行平滑滤波；选取单元，用于依据滤波后的方向形态纹理信息影像，选取满足预设条件的方向形态纹理彳目息影像；确定单元，用于将所述满足预设条件的方向纹理信息影像确定为所述待处理光学影像的方向不变形态纹理特征影像。优选地，所述重建模块包括:重建单元，用于将所述第一分割影像作为种子，将所述第二分割影像作为掩膜，进行基于数学形态学的重建计算，所述第一阈值大于所述第二阈值。优选地，所述装置还包括:可视化模块，用于依据所述待处理光学影像的类型，将所述待处理光学影像与所述倒塌房屋图像进行融合可视化。本发明实施例所述的从高分辨率光学影像中提取倒塌房屋的方法及装置，使用结构元素，提取待处理光学影像中的方向纹理信息，并依据方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像，再对方向不变形态纹理特征影像进行阈值分割后重建，从而获取待处理光学影像中的倒塌房屋图像，因为利用方向形态纹理信息计算得到方向不变形态纹理特征影像，所以，能够充分刻画倒塌房屋的结构方向信息，从而提高提取结果的准确性，同时，又因为提取方向形态纹理信息的方法基于数据形态学，因此，相比与利用灰度共生矩阵的方法，减小了计算量，由此可见，本发明实施例所述的方法及装置，能够兼顾计算效率与计算精度。
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本发明实施例公开的一种从高分辨率光学影像中提取倒塌房屋的方法的流程图；图2为本发明实施例公开的又一种从高分辨率光学影像中提取倒塌房屋的方法的流程图；图3为本发明实施例公开的一种从高分辨率光学影像中提取倒塌房屋的方法中的结构元素的示意图；图4为本发明实施例公开的又一种从高分辨率光学影像中提取倒塌房屋的方法的流程图；图5为本发明实施例公开的一种从高分辨率光学影像中提取倒塌房屋的装置的结构示意图。
具体实施例方式本发明实施例公开了一种从高分辨率光学影像中提取倒塌房屋的方法及装置，其核心发明点在于，使用结构元素提取待处理光学影像中的方向纹理信息，并依据提取出的方向纹理信息，计算待处理光学影像中的方向不变形态纹理特征影像，因为方向不变纹理特征影像能够更为细致地刻画倒塌房屋的结构方向信息，因此，在此基础上对图像进行阈值分割后重建，能够得到精确的倒塌房屋信息，并且，基于数学形态学的方向纹理信息提取方法与传统的基于灰度共生矩阵的方法而言，能够大大减小计算量，从而提高倒塌房屋提取的效率。下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。本发明实施例公开了一种从高分辨率光学影像中提取倒塌房屋的方法，如图1所示，包括:SlOl:使用结构元素，提取待处理光学影像中的方向形态纹理信息；所述结构元素，是基于数学形态学的图像处理方法中常用的处理工具，具体表现为矩阵的形式。通常，首先将灾后高分影像在计算机上进行降采样可视化显示，可以利用已有的基础地理数据确定居民区的大致范围，也可以通过在降分辨率后的小图像上目视确定居民区的大致范围。获取高分辨率影像中居民区自然聚集区域后，按照次序将其编号，本实施例中所述的待处理光学影像是指编号后的居民区自然聚集区域。需要说明的是，当一个居民区自然聚集区域内的影像受到地形阴影、云、雾霾等成像条件影响出现明显的对比度差异时，需要按照明显的对比度边界将居民区进行再次分块，则再次分块后的子影像作为待处理的光学影像。即本实施例中所述的待处理的光学影像为分块后的、满足处理条件的、包含居民区的子影像。S102:依据所述方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像；S103:利用预设的第一阈值分割所述方向不变纹理特征影像，得到第一分割影像，利用预设的第二阈值分割所述方向不变纹理特征影像，得到第二分割影像；S104:将所述第一分割影像及所述第二分割影像进行基于数学形态学的重建，从而获取倒塌房屋图像。本实施例所述的方法，通过提取出待处理光学影像中的方向形态纹理信息，进而计算出待处理光学影响的方向不变纹理特征影像，因为方向不变纹理特征影像能够更细致地刻画倒塌房屋的结构信息，因此，后续基于方向不变纹理特征影像提取出的倒塌房屋的准确性更高，并且，由于待处理光学影响的方向不变纹理特征影像的计算基于数学形态学，因此能够提高计算的速度，从而提到倒塌房屋提取的效率。本发明实施例公开的又一种从高分辨率光学影像中提取倒塌房屋的方法，应用于以下场景:针对<lm空间分辨率影像数据进行处理，提取结构性粉碎性损坏而倒塌的房屋。模型假设影像中倒塌房屋包含精细的纹理结构，且房屋倒塌后形成目标面积所占像元尺寸不小于5X5个像元。要求影像像元尺寸与空间分辨率的比例是1:1，不处理经过缩放后像元尺寸大于空间分辨率的影像，因为，利用常用缩放技术进行像元尺寸放大影响最小纹理信息的尺度，且不增加信息量。如图2所示，包括:S201:分别使用预设数量不同方向的、最小单元的线状结构元素对待处理光学影像进行高帽变换，从而得到所述预设数量的不同方向的形态纹理信息影像；本实施例中，优选结构元素的数量为4个，结构元素的形状为线状，结构元素的长度的3个像素，中心点为第二个像素，即中间的像素。四组线状结构元素的方向分别为0°，45° ,90°和135°，结构元素如图3所示，图3中(a)、(b)、(c)、(d)分别表示0°，45°，90°和135°的结构元素，其中黑色部分表示各个结构元素中的中心点，每一个方格表示各个结构元素中的像素，即每个结构元素中包含3个像素。本实施例中，采用最小长度的线状结构元素的目的在于，假设绝大多数倒塌房屋在影像中具有精细的纹理结构，这个假设已在实际应用中得到验证。一种极端情况，通常出现在空间分辨率小于0.2m时，由于空间分辨率很高，也即像元尺寸很大，造成倒塌房屋最小纹理单元大于3个像元，此时，需要增加结构元素尺度至与倒塌房屋最小纹理相匹配的尺度。实际应用中，为了快速获取大面积区域灾情信息，这种数据获取方式不多，存在这种情况时，建议空间分辨率小于0.2m的影像可以利用像元平均方法降分辨率重采样到0.2m再进行处理，基本不影响倒塌房屋信息的提取精度，同时，提高计算效率。数学形态学中的高帽变换不依赖于局部或全局统计特性，能够将不同方向的运算及不同像素位置的运算并行处理，极大地提高了计算效率，使用上述四组结构元素对待处理光学影像进行高帽变换的具体过程如以下公式表示:10=Tophat (I, SE (I, 3, O))；
I45=Tophat(I, SE (I, 3, 45))；I135=Tophat (I, SE (I, 3，135))；I90=Tophat (I, SE (I, 3, 90)) 其中，其中，I表示原始影像，Tophat为高帽变换，SE表示形态变换所用的结构元素，1,3,0/45/135/90分别表不结构兀素的宽、长、角度。I。、145、1135、I90分别表不生成的方向形态特征影像。通常，高帽变换通过利用原始影像与其开运算结果进行差值运算实现，其中，开运算可以由任一种具有方向性的低通滤波替代，这里不做限定。需要说明的是，本实施例中优选4组结构元素，在实际应用中，结构元素的数量还可以为8个、16个等，也可以为2个，因此，本发明实施例中并不做限定，选定4个结构元素能够较好地兼顾计算效率与对方向纹理信息的刻画精度。每一个方向结构元 素高帽变换的结构都是一幅特征图像，代表该方向上空间尺度小于或等于结构元素长度的细节信息。需要强调的是，方向纹理信息是高分影像中区分倒塌房屋与非倒塌房屋信息的一个重要标志。这主要是由于正常居民区房屋大多属于人工目标，具有相对明显的几何和方向性，但是，受损后的建筑在影像中的结构方向性不明显，因此，方向纹理信息有助于对倒塌房屋进行更准确的刻画。S202:对所述方向形态纹理信息影像进行平滑滤波；本实施例中，优选使用低通均值滤波对方向形态纹理信息影像进行滤波，如以下公式所示:F0=mean (10, size)；F45=mean (I45, size)；F135=mean (I135, size)；F90=mean (I90, size) 其中，F表示滤波后的特征图像，mean代表均值低通滤波，size表示滤波窗口大小。均值滤波所采用的窗口大小需要依据图像空间分辨率和兴趣房屋类型大小共同决定，兴趣建筑实际尺寸B固定，图像分辨率S可变，因此，空间分辨率S越高，窗口 W越大。呈线性比例变化，B=W*S。需要说明，多方向形态纹理特征对窗口大小具有一定鲁棒性，利用最小尺寸以上的窗口进行滤波能够较好增强大面积倒塌房屋，同时不明显损失小尺寸的倒塌房屋信息。S203:依据滤波后的方向形态纹理信息影像，选取满足预设条件的方向形态纹理息影像;S204:将所述满足预设条件的方向形态纹理信息影像确定为所述待处理光学影像的方向不变形态纹理特征影像；本实施例中，可以选择方向形态纹理信息影像中特征值最小的方向形态纹理信息影像作为待处理光学影像的方向不变形态纹理特征影像，即可以如以下公式所示确定方向不变形态纹理特征影像: 0IF=minmum (F0, F45, F135, F90)。
其中，minmum表示取最小值，OIF表示方向不变形态纹理特征影像。需要说明的是，除了使用上述选择最小值的方法确定方向不变形态纹理特征影像，也可以将所有滤波后的方向形态纹理信息影像均值影像确定为方向不变形态纹理特征影像。S205:利用预设的第一阈值分割所述方向不变纹理特征影像，得到第一分割影像，利用预设的第二阈值分割所述方向不变纹理特征影像，得到第二分割影像；本实施例中，第一阈值与第二阈值可以通过人工交互的方式获取，具体的获取过程包括:将原始高分辨率光学影像和所述的方向不变纹理特征影像进行几何联接显示，交互观察方向不变纹理特征影像中高亮度区域及其在原始高分辨率光学影像中对应区域，利用多边形工具选择倒塌房屋区域作为样本信息，在不同的位置选择大于5个的样本区域。统计旋转不变纹理特征影像的直方图，并从大到小进行二值分割，逐一判断所选的样本区域与分割后目标图像的关系，当被分割出来的目标与每一个样本区域的交集都达到对应样本区域面积的10%时，该阈值作为第一阈值；进一步减小阈值，当被分割出来的目标与每一个样本区域的交集都达到对应样本区域面积的90%时，该阈值作为第二阈值。需要说明的是，旋转不变纹理特征影像与原始高分辨影像具有相同的大小，可以直接在原始高分辨率影像中获取倒塌房屋区域，再将其转绘到旋转不变纹理特征影像对应位置上，这里交互利用二者进行样本区域选取，是因为旋转不变纹理特征影像中高亮度信息能够协助更快地定位到倒塌房屋区域。其中，样本区域个数越多结果更加稳定，实际应用中个数大于5块后，结果变化不明显。第一阈值获取参数可在10%±5%之间浮动，第二阈值获取参数可在90%±5%之间浮动，不同阈值获取参数对获取的倒塌房屋目标边界有一定影响，但对于减灾应用中最关心的检测和定位倒塌房屋目标的结果影响不明显。利用第一阈值与第二阈值多方向不变形态纹理特征影像进行分割的原理可以由下述公式表示:
权利要求
1.一种从高分辨率光学影像中提取倒塌房屋的方法，其特征在于，包括: 使用结构元素，提取待处理光学影像中的方向形态纹理信息； 依据所述方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像； 利用预设的第一阈值分割所述方向不变纹理特征影像，得到第一分割影像，利用预设的第二阈值分割所述方向不变纹理特征影像，得到第二分割影像； 将所述第一分割影像及所述第二分割影像进行基于数学形态学的重建，从而获取倒塌房屋图像。
2.根据权利要求1所述的方法，其特征在于，所述使用结构元素，提取待处理光学影像中的方向形态纹理信息包括: 分别使用预设数量不同方向的、最小单元的线状结构元素对待处理光学影像进行高帽变换，从而得到所述预设数量的不同方向的形态纹理信息影像。
3.根据权利要求2所述的方法，其特征在于，所述依据所述方向形态纹理信息，计算所述待处理光学影像的方向不变 形态纹理特征影像包括: 对所述方向形态纹理信息影像进行平滑滤波； 依据滤波后的方向形态纹理信息影像，选取满足预设条件的方向形态纹理信息影像； 将所述满足预设条件的方向纹理信息影像确定为所述待处理光学影像的方向不变形态纹理特征影像。
4.根据权利要求3所述的方法，其特征在于，所述将所述第一分割影像及所述第二分割影像进行基于数学形态学的重建包括: 将所述第一分割影像作为种子，将所述第二分割影像作为掩膜，进行基于数学形态学的重建计算，所述第一阈值大于所述第二阈值。
5.根据权利要求1所述的方法，其特征在于，还包括: 依据所述待处理光学影像的类型，将所述待处理光学影像与所述倒塌房屋图像进行融合可视化。
6.一种从高分辨率光学影像中提取倒塌房屋的装置，其特征在于，包括: 提取模块，用于使用结构元素，提取待处理光学影像中的方向形态纹理信息； 计算模块，用于依据所述方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像； 分割模块，用于利用预设的第一阈值分割所述方向不变纹理特征影像，得到第一分割影像，利用预设的第二阈值分割所述方向不变纹理特征影像，得到第二分割影像； 重建模块，用于将所述第一分割影像及所述第二分割影像进行基于数学形态学的重建，从而获取倒塌房屋图像。
7.根据权利要求6所述的装置，其特征在于，所述提取模块包括: 高帽变换单元，用于分别使用预设数量不同方向的、最小单元的线状结构元素对待处理光学影像进行高帽变换，从而得到所述预设数量的不同方向的形态纹理信息影像。
8.根据权利要求6所述的装置，其特征在于，所述计算模块包括: 滤波单元，用于对所述方向形态纹理信息影像进行平滑滤波； 选取单元，用于依据滤波后的方向形态纹理信息影像，选取满足预设条件的方向形态纹理 目息影像； 确定单元，用于将所述满足预设条件的方向纹理信息影像确定为所述待处理光学影像的方向不变形态纹理特征影像。
9.根据权利要求6所述的装置，其特征在于，所述重建模块包括: 重建单元，用于将所述第一分割影像作为种子，将所述第二分割影像作为掩膜，进行基于数学形态学的重建计算，所述第一阈值大于所述第二阈值。
10.根据权利要求6所述的装置，其特征在于，还包括: 可视化模块，用于依据所述待处理光学影像的类型，将所述待处理光学影像与所述倒塌房屋图像进行融合可视 化。
全文摘要
本发明实施例提供了一种从高分辨率光学影像中提取倒塌房屋的方法及装置，使用结构元素，提取待处理光学影像中的方向纹理信息，并依据方向形态纹理信息，计算所述待处理光学影像的方向不变形态纹理特征影像，再对方向不变形态纹理特征影像进行阈值分割后重建，从而获取待处理光学影像中的倒塌房屋图像，因为利用方向形态纹理信息计算得到方向不变形态纹理特征影像，所以，能够充分刻画倒塌房屋的结构方向信息，从而提高提取结果的准确性，同时，又因为提取方向形态纹理信息的方法基于数据形态学，因此，相比与利用灰度共生矩阵的方法，减小了计算量，由此可见，本发明实施例所述的方法及装置，能够兼顾计算效率与计算精度。
文档编号G06K9/00GK103218597SQ201310080188
公开日2013年7月24日 申请日期2013年3月13日 优先权日2013年3月13日
发明者李利伟, 张兵, 高连如 申请人:中国科学院对地观测与数字地球科学中心