基于高分辨率SAR图像的城市高层建筑物高度提取方法及装置与流程

本发明涉及建筑物高度提取
技术领域：
，特别涉及一种基于高分辨率sar图像的城市高层建筑物高度提取方法。
背景技术：
：合成孔径雷达(syntheticapertureradar，sar)作为一种主动微波遥感手段，具有全天时、全天候、可穿透云层快速成像的能力，弥补了光学遥感的不足，在应急侦察、灾害监测等需要快速响应的城市遥感应用中发挥着重要作用，而城市建筑物是城市目标监测的主要对象。高分辨率是sar技术发展的主要趋势之一。加拿大的radarsat-2(2007)、意大利的cosmo-skymed(2007-2008)、德国的terrasar-x(2007)和tandem-x(2010)、中国的高分三号卫星(2016)等多极化、高分辨率卫星图像分辨率已突破到1m，其中terrasar-x于2013年推出的staringspotlight(st)成像模式更是将图像的方位向分辨率提高到0.24m。米级、亚米级高分辨率sar图像的广泛获取为研究建筑物高度提取的新方法、新理论提供了数据支撑。然而，高分辨率sar图像在带来建筑物目标丰富细节的同时，也带来了独立建筑物目标成像特征破碎和成像特征多变等诸多建筑物信息提取难题。同时sar固有的几何畸变和斑噪影响，以及复杂的城市背景干扰等因素，也进一步增加了建筑物目标在sar图像上的复杂性。按照使用的数据类型的不同，目前主要的建筑物高度提取方法可以归纳为：(1)干涉测量法，即利用干涉sar(insar)数据提取建筑物目标高度信息，如tomosar技术。(2)多视向立体测量法。利用同侧或异侧进行多视向立体测量，同侧立体测量可以利用不同入射角成像实现，异侧立体测量可以利用升降轨成像实现。这两种方法各有利弊，异侧立体成像通常具有较大的视差，对于提高建筑物高度估算精度是有利的，但视角差别越大，目标在sar图像上的差异越大，给目标识别造成的困难也就越大。(3)与光学图像融合法。不同数据源遥感图像往往可以实现信息互补，sar图像中建筑物的显著性线状特征，较为容易确定潜在目标，却无法得到建筑物的完整轮廓；然而，光学图像上建筑物有完整的屋顶结构有利于确定完整边界信息，但建筑物的检测难度较大，特别是对于建筑物分布密集的城区。综合利用建筑物在不同数据类型图像中的特点来提高建筑物目标识别和三维重建的精度是一种有效手段。(4)基于单景高分辨率sar图像的建筑物高度提取方法。单景高分辨率sar图像上，建筑物高度的反演主要利用图像中灰度的层次变化关系以及阴影信息。对于立于平坦地面、并且分散的独立建筑物，利用完整的阴影信息在估计出可能的房顶结构后，可较完美地反算出建筑物高度。现有方法中，利用阴影或者叠掩信息来反演建筑物高度面临以下问题。首先，需要确保提取完整的阴影或者叠掩范围，导致该方法适用于稀疏建筑物，不能应用于建筑物密集的城区，因为建筑物的阴影或者叠掩通常受到临近目标的干扰而不完整；其次，阴影区域较难与道路、广场等等暗目标区分；最后，人工建筑物墙面的不均匀性导致sar图像中叠掩区的灰度一致性较差，导致叠掩区的检测精度低，从而影响高度反演精度。技术实现要素：有鉴于克服现有技术中由于高分辨率特性带来的叠掩面破碎特性，使建筑物高度反演精度不高的问题，本发明提出了一种基于高分辨率sar图像的城市高层建筑物高度提取方法及装置。本发明提供了一种基于高分辨率sar图像的城市高层建筑物高度提取方法，包括：对原始高分辨率sar图像采用非局部滤波，对滤波后sar图像分割获取匀质区域，然后从中提取阴影区和高亮区，并利用阴影区和高亮区在距离向的位置关系确定建筑物候选区；针对独立高亮的建筑物候选区提取其最小外接矩形，计算矩形区域内距离向上每一行的亮线长度并记录其所在行位置，并基于该长度和行位置属性利用聚类算法对其进行聚类；获取聚类后的每个类别所对应的二值区域，基于叠掩区平行四边形几何形状的限制，利用仿射变换和最小二乘方法对其进行拟合，获取优化后的建筑物候选区；基于叠掩区的像元长度、像元采样间隔和雷达入射角计算建筑物高度。优选地，所述聚类算法采用ap聚类算法，包括以下步骤：统计最小外接矩形区域内每一行上的亮线长度为l，行坐标为r，则每一行构建为一个数据点x(l,r)，最终获取数据集{x1,x2...,xn}，其中，l刻画的是叠掩区长度的相似性，r刻画的是距离上的邻近性，令s是一个刻画点之间相似度的矩阵，s(i,j)表示xi与xj之间的相似度；.采用ap算法进行迭代计算，以更新两个吸引信息矩阵和归属信息矩阵，所述吸引信息矩阵和归属信息矩阵存在交替。优选地，所述吸引信息矩阵r中的每个元素r(i,k)描述的是数据对象k适合作为数据对象i的聚类中心的程度，表示从i到k的消息；所述归属信息矩阵a中的每个元素a(i,k)描述的是数据对象i选择数据对象k作为其聚类中心的适合程度，表示从k到i的消息。优选地，(1)所述吸引信息矩阵的初始值为0，迭代计算方式如下：其中s(i,k)定义为s(i,k)＝-α·(xi[l]-xk[l])2-β·(xi[r]-xk[r])2当α与β都为1时，s(i,k)为xi与xk之间的欧式距离；(2)所述归属信息矩阵的初始值为0，迭代计算方式如下：在进行上述迭代时，如果这些决策经过若干次迭代之后保持不变、或者算法执行超过设定的迭代次数、或者一个小区域内的关于样本点的决策经过数次迭代后保持不变，则算法结束。优选地，在ap算法更新信息时引入了衰减系数λ(0＜λ＜1)以避免震荡，每条信息被设置为它前次迭代更新值的λ倍加上本次信息更新值的1-λ倍，即：rt+1(i,k)＝(1-λ)rt+1(i,k)+λrt(i,k)at+1(i,k)＝(1-λ)at+1(i,k)+λat(i,k)。优选地，所述仿射变换包括：获取该二值区域的最小外接矩形，基于该矩形边界通过仿射变换最大限度匹配二值区域的边界；令矩形边界点为e0，(x0,y0)∈e0；二值区域真实边界为ep，(xp,yp)∈ep，仿射变化得到的边界ef，(xf,yf)∈ef，则仿射变换的公式为：其中，sx为水平方向的拉伸因子，tx为水平方向的偏移量，α为错切参数，且优选地，所述拟合包括：为获取ef与ep的最佳匹配，将(xp,yp)作为(xf,yf)的替代值代入上述仿射变换公式，基于最小二乘方法获取sx，tx和α三个参数，计算公式为：其中n为ef中参与参数计算的总像元个数。优选地，所述拟合还包括：首先，需要建立(xp,yp)与(xf,yf)之间的一一对应关系；叠掩区形成的平行四边形的水平边界上的点不参与计算，其余边界点可根据yp＝yf＝y0建立稳定的一一对应关系，并参与仿射变换参数的解算。优选地，基于叠掩区的像元长度、像元采样间隔和雷达入射角计算建筑物高度包括：在获取叠掩区平行四边形后，计算其距离向像元长度为n，配合图像距离向采样间隔δ，雷达入射角θ，建立建筑物高度反演公式为：本发明还公开一种基于高分辨率sar图像的城市高层建筑物高度提取装置，包括：选取区域模块，其配置为对原始高分辨率sar图像采用非局部滤波，对滤波后sar图像分割获取匀质区域，然后从中提取阴影区和高亮区，并利用阴影区和高亮区在距离向的位置关系确定建筑物候选区；聚类分析模块，其配置为针对独立高亮的建筑物候选区提取其最小外接矩形，计算矩形区域内距离向上每一行的亮线长度并记录其所在行位置，并基于该长度和行位置属性利用聚类算法对其进行聚类；建筑物叠掩区提取模块，其配置为获取聚类后的每个类别所对应的二值区域，基于叠掩区平行四边形几何形状的限制，利用仿射变换和最小二乘方法对其进行拟合，获取优化后的建筑物候选区；计算模块，其配置为基于叠掩区的像元长度、像元采样间隔和雷达入射角计算建筑物高度。本发明在分割和亮、暗像元提取的基础上初步获取建筑物叠掩区，通过距离向长度聚类和仿射变换获取完整的建筑物叠掩区范围，弥补了高分辨率特性带来的叠掩面破碎特性，提高建筑物高度反演精度。附图说明图1a为本实施例所使用的terrasar-xstaringspotlight(st)图像；图1b为本实施例所使用的光学图像；图2为本发明所提出的一种高分辨率sar图像的城市高层建筑物高度提取方法的流程图；图3为测试图像的滤波结果图；图4为测试图像的分割结果图；图5为测试图像的建筑物候选区提取结果；图6为某一独立高亮建筑物候选区的最小外接矩形提取结果；图7为图6所示区域的聚类结果；图8a和图8b为图6所示区域的优化后的平行四边形叠掩区；图9为测试图像的整体叠掩区提取结果。具体实施方式为使本领域技术人员更好的理解本发明的技术方案，下面结合附图和具体实施方式对本发明作详细说明。本发明面向单景高分辨率sar图像的建筑物高度提取，且面向城市高大建筑，主要利用阴影、叠掩和二次散射亮线来反演建筑物高度本发明针对分米级高分辨率sar图像，解决了利用单一二次散射亮线提取建筑物高度偏小以及高层建筑物叠掩区相对破碎的问题；利用局部的叠掩区平行四边形拟合和优化，建筑物高度提取误差小，为城市3d重建提供依据。下面结合附图和实施例对本发明做进一步的详细说明。本实施例采用的测试图像如图1a和图1b所示。其中，图1a为terrasar-xstaringspotlight图像，斜距图像，距离向采样间隔0.45m，雷达中心入射角32.9°；图1b为对应的光学图像。图中主要包含a、b、c、d四种类型建筑，其高度已实测获取。在利用本发明中的方法进行城市高层建筑物高度提取时，如图2所示，主要包括以下几个步骤：步骤一：对原始高分辨率sar图像采用非局部滤波sarblock-matching3-d(sar-bm3-d)方法，滤波结果如图3所示；对滤波后sar图像采用meanshift分割，结果如图4所示；基于图像先验知识，设置简单的灰度阈值为50和100，分别提取阴影区和高亮区，并利用二者在距离向的位置关系确定建筑物候选区(阴影区存在于高亮叠掩区的远距端)，结果如图5所示。步骤二：针对每一个独立高亮建筑物候选区，提取其最小外接矩形，如图6所示，计算矩形区域内距离向上每一行的亮线长度并记录其所在行位置。基于该长度和位置属性，利用ap聚类算法对其进行聚类。ap聚类算法不需要事先确定类别的个数，基于数据点间的迭代计算来传递信息，其算法过程如下：(1)数据点的构建统计最小外接矩形区域内每一行上的亮线长度为l，行坐标为r，则每一行构建为一个数据点x(l,r)，最终获取数据集{x1,x2...,xn}，其中，l刻画的是叠掩区长度的相似性，r刻画的是距离上的邻近性，令s是一个刻画点之间相似度的矩阵，s(i,j)表示xi与xj之间的相似度。(2)ap算法进行迭代计算，以更新两个矩阵(吸引信息矩阵和归属信息矩阵)，且两个矩阵存在交替。吸引信息矩阵r：r(i,k)描述的是数据对象k适合作为数据对象i的聚类中心的程度，表示从i到k的消息。归属信息矩阵a：a(i,k)描述的是数据对象i选择数据对象k作为其聚类中心的适合程度，表示从k到i的消息。两个矩阵的初始值为0，迭代计算方式为：吸引信息迭代其中s(i,k)在本发明中定义为：s(i,k)＝-α·(xi[l]-xk[l])2-β·(xi[r]-xk[r])2本实施例中，α＝0.5，β＝1。归属信息迭代对以上步骤进行迭代，本实施例中设置迭代次数少于1000次。为了避免震荡，ap算法更新信息时引入了衰减系数λ(0＜λ＜1)，本实施例中设置其值为0.9。每条信息被设置为它前次迭代更新值的λ倍加上本次信息更新值的1-λ倍，即：rt+1(i,k)＝(1-λ)rt+1(i,k)+λrt(i,k)at+1(i,k)＝(1-λ)at+1(i,k)+λat(i,k)按照上述计算过程，针对图6的聚类结果如图7所示，共得到9个聚类中心。步骤三：聚类后，获取每个类别所对应的二值区域，基于叠掩区平行四边形几何形状的限制，利用仿射变换和最小二乘方法对其进行拟合，获取优化后的建筑物候选区。(1)仿射变换首先获取该二值区域的最小外接矩形，基于该矩形边界，通过仿射变换最大限度匹配二值区域的边界。令矩形边界点为e0，(x0,y0)∈e0；二值区域真实边界为ep，(xp,yp)∈ep，仿射变化得到的边界ef，(xf,yf)∈ef，则仿射变换的公式为：其中，sx为水平方向的拉伸因子，tx为水平方向的偏移量，α为错切参数，且(2)最小二乘拟合为获取ef与ep的最佳匹配，将(xp,yp)作为(xf,yf)的替代值代入上述仿射变换公式，基于最小二乘方法获取sx，tx和α三个参数。计算公式为：其中n为ef中参与参数计算的总像元个数。首先，需要建立(xp,yp)与(xf,yf)之间的一一对应关系。由于叠掩区形成的平行四边形必定存在两条水平的边界，水平边界上的点难以获取稳定的对应关系，因此不参与计算。其余边界点可根据yp＝yf＝y0这一特征建立稳定的一一对应关系，并参与仿射变换参数的解算。基于上述计算过程，获取图6所示区域的优化后的叠掩区域如图8a和图8b中平行四边形所示，整体图像的优化结果如图9所示。步骤四：基于叠掩区像元长度、像元采样间隔和雷达入射角计算建筑物高度。获取叠掩区平行四边形后，计算其距离向像元长度为n。配合图像距离向采样间隔δ，雷达入射角θ，建立建筑物高度反演公式为：本实施例中a、b、c、d四类建筑物的高度计算结果如表1所示，参考实测建筑物高度，本发明的建筑物高度提取平均误差为2m。表1建筑物abcd实验结果(m)45.519.548.720.6实测结果(m)48.72247.822本发明还公开一种基于高分辨率sar图像的城市高层建筑物高度提取装置，包括：选取区域模块，其配置为对原始高分辨率sar图像采用非局部滤波，对滤波后sar图像分割获取匀质区域，然后从中提取阴影区和高亮区，并利用阴影区和高亮区在距离向的位置关系确定建筑物候选区；聚类分析模块，其配置为针对独立高亮的建筑物候选区提取其最小外接矩形，计算矩形区域内距离向上每一行的亮线长度并记录其所在行位置，并基于该长度和行位置属性利用聚类算法对其进行聚类；建筑物叠掩区提取模块，其配置为获取聚类后的每个类别所对应的二值区域，基于叠掩区平行四边形几何形状的限制，利用仿射变换和最小二乘方法对其进行拟合，获取优化后的建筑物候选区；计算模块，其配置为基于叠掩区的像元长度、像元采样间隔和雷达入射角计算建筑物高度。应当理解的是，本领域技术人员结合上文的详细介绍，可以得到本发明的实施例的基于高分辨率sar图像的城市高层建筑物高度提取装置的模块化架构，也可以依据本领域技术常识对本发明的实施例的基于高分辨率sar图像的城市高层建筑物高度提取装置做不超出本发明所要求保护范围内的变化扩展，因此，本申请不再在说明书附图中以附图形式赘述相应的模块化架构。本发明的优点在于：(1)本发明提出一种高分辨率sar图像的城市高层建筑物高度提取方法，该方法利用聚类分析初步确定高层建筑物叠掩区范围；可实现不同叠掩面的分离。(2)本发明提出一种高分辨率sar图像的城市高层建筑物高度提取方法，该方法精度高。利用局部的叠掩区平行四边形拟合和优化，克服图像高分辨率特性带来的叠掩面匀质性较低问题，建筑物高度提取误差小。(3)本发明提出一种高分辨率sar图像的城市高层建筑物高度提取方法，该方法在综合利用图像亮度、几何形状等信息对建筑物高度提取方面进行了开创性的工作。以上实施例仅为本发明的示例性实施例，不用于限制本发明，本发明的保护范围由权利要求书限定。本领域技术人员可以在本发明的实质和保护范围内，对本发明做出各种修改或等同替换，这种修改或等同替换也应视为落在本发明的保护范围内。当前第1页12