本发明涉及到建筑物遥感图像处理技术领域,具体地说,是一种基于二进制特征分类器的建筑物提取方法。
背景技术:
随着数字城市、智慧市政等概念的提出,城市规划、用地变化监测、建筑物三维建模等业务需求飞速增长,而建筑物提取是支撑这些业务实现的核心技术。因此,如何从高分辨率遥感影像中进行准确快速的建筑物提取成为当前的热门技术领域。
目前所开展的适用于遥感影像的建筑物提取的方法研究,大部分的是以图像分割方法为基础。这类方法通常需要具有专业经验的人员进行人工操作,效率低且成本高。
由此可见,亟需一种快速、高精度的、自动建筑物提取方法。
技术实现要素:
针对现有技术的不足,本发明的目的是提供一种基于二进制特征分类器的建筑物提取方法,利用二进制特征分类器对多尺度正射影像中的超像素分割点进行分类,并根据分类信息实现对建筑物快速、高精度的、自动提取。
为达到上述目的,本发明采用的技术方案如下:
一种基于二进制特征分类器的建筑物提取方法,其关键在于包括以下步骤:
步骤1:对原始正射图像进行降采样;
步骤2:将降采样后的正射图像从rgb颜色空间转换至lab颜色空间;
步骤3:在lab颜色空间下对图像进行slic超分割处理;
步骤4:根据超分割结果标记出交叉点;
步骤5:利用二进制特征分类器对交叉点进行二分类,获得建筑物交叉点和非建筑物交叉点;
步骤6:将非建筑物交叉点所连接的图斑进行剔除,得到粗略提取的建筑物轮廓;
步骤7:设定最小建筑物区域像素阈值,并剔除错误检测的建筑物;
步骤8:获得建筑物轮廓的外接矩形,得到建筑物区域,并将其位置还原至原始正射图像上;
步骤9:在原始正射图像上,提取建筑物区域,并将其颜色空间由rgb转换为lab;
步骤10:在lab颜色空间下,对建筑物区域图像进行slic超分割处理;
步骤11:根据超分割结果标记出新交叉点;
步骤12:利用二进制特征分类器对新交叉点进行三分类,获得建筑物内部点、建筑物边缘点和非建筑物点;
步骤13:将非建筑物交叉点所连接的图斑进行剔除,得到精细提取的建筑物轮廓;
步骤14:设定最小建筑物区域像素阈值,剔除错误检测的建筑物;
步骤15:提取建筑物轮廓曲线,获得建筑物轮廓图像。
进一步的,步骤1中所述降采样为将原始正射图像的像素分辨率降采样至一米每像素。
进一步的,步骤3与步骤10中所述slic超分割处理的步骤一致。
更进一步的,步骤3中所述slic超分割处理的具体步骤为:
步骤3-1:初始化种子点,按照设定的超像素个数,在图像内均匀分配种子点,计算相邻种子点的距离l:
其中,n为图像的像素点总数,k为超分割块设定个数;
步骤3-2:在种子点的n×n邻域内重新选择种子点,通过计算该邻域内所有像素点的梯度值,将种子点移到该邻域内梯度最小的地方;
步骤3-3:在每个种子点的2l×2l的邻域内为每个像素点分配类标签,对于每个搜索到的像素点,分别计算它与种子点的距离;
步骤3-4:重复步骤3-2和步骤3-3,迭代得到最优结果。
再进一步的,步骤3-3中像素点与种子点的距离计算方法如下:
其中,d’表示种子点的距离,nc为最大颜色空间,nl是类内最大空间距离,dc表示颜色距离,dl表示空间距离,dc和dl的计算公式如下:
进一步的,步骤5中与二进制特征分类器相匹配的二进制混合描述符包含方向特征和描述符两部分,具体描述方式为:
方向特征为:
θ=atan2(m01,m10),
其中,θ为角度特征,m01和m10为每个图片分割块的矩,计算公式如下:
mpq=∑x,yxpyqi(x,y),
描述符为:
fn(p)=σ0≤j<4σ0≤i<n2n×j+iτj(p;ai,bi),
其中,fn(p)表示分割块p的描述符,τ1为纹理特征测试对,τ2为相似度特征测试对,τ3和τ4为颜色信息测试对,计算公式分别如下:
其中,pl(a)和pl(b)分别表示l通道上以a为中心和以b为中心的分割块p,pa(a)和pa(b)分别表示a通道上以a为中心和以b为中心的分割块p,pb(a)和pb(b)分别表示b通道上以a为中心和以b为中心的分割块p,s为相似度阈值,ca为a通道的颜色阈值,cb为b通道的颜色阈值。
进一步的,步骤12中与步骤5中,与所述二进制特征分类器相匹配的二进制混合描述符的描述方式一致。
本方案包括建筑物区域提取和建筑物轮廓提取两个部分。在第一个部分中,首先对原图像进行降采样处理,转换至lab颜色空间,然后在lab颜色空间之下实施slic超分割获得交叉点,利用二进制特征分类器对交叉点进行二分类获得建筑物交叉点和非建筑物交叉点,从而获得建筑物区域,实现建筑物的粗略提取。第二个部分是在第一个部分的基础之上进行的,首先将第一部分提取出的建筑物区域图转换至lab颜色空间,然后在lab颜色空间之下实施slic超分割获得新交叉点,利用二进制特征分类器对新交叉点进行三分类,获得建筑物内部点、建筑物边缘点和非建筑物点,从而得到建筑物轮廓。
本发明的显著效果是:分为建筑物区域提取和建筑物轮廓提取两个部分,分别实现了建筑物轮廓的粗提取和精提取两个过程,并利用二进制特征分类器对特征进行分类,实现了建筑物轮廓的快速准确提取;与传统方法相比,本发明精度更高、速度更快且算法过程无交互。
附图说明
图1是本发明的方法流程图;
图2是本发明方法的主要阶段图,其中(a)是原始正射图像,(b)是经步骤3对降采样后正射图像进行slic超分割的结果,(c)是经步骤8得到的建筑物区域的结果,(d)是经步骤10对建筑物区域进行超分割后的结果,(e)是采用本发明方法提取的建筑物轮廓的最终结果。
具体实施方式
下面结合附图对本发明的具体实施方式以及工作原理作进一步详细说明。
如图1所示,一种基于二进制特征分类器的建筑物提取方法,具体按照以下步骤进行:
步骤1:对原始正射图像分别进行降采样,从而将原始正射图像降采样至一米每像素;
步骤2:将降采样后的正射图像从rgb颜色空间转换至lab颜色空间,以更好地利用原图像中的纹理细节特征和颜色特征;
步骤3:在lab颜色空间下对图像进行slic(简单的线性迭代聚类)超分割处理,具体处理步骤为:
步骤3-1:初始化种子点,按照设定的超像素个数,在图像内均匀分配种子点,计算相邻种子点的距离l:
其中,n为图像的像素点总数,k为超分割块设定个数;
步骤3-2:在种子点的n×n邻域内重新选择种子点,通过计算该邻域内所有像素点的梯度值,将种子点移到该邻域内梯度最小的地方;
步骤3-3:在每个种子点的2l×2l的邻域内为每个像素点分配类标签,对于每个搜索到的像素点,分别计算它与种子点的距离;
其中,像素点与种子点的距离计算方法如下:
其中,d’表示种子点的距离,nc为最大颜色空间,用固定常数值代替,本例中取为10;nl是类内最大空间距离,定为s,适合每个聚类。dc表示颜色距离,dl表示空间距离,dc和dl的计算公式如下:
步骤3-4:重复步骤3-2和步骤3-3,迭代10次,获得迭代出的最优结果。
步骤4:根据slic超分割标记出交叉点;
步骤5:利用二进制特征分类器对交叉点进行二分类,获得建筑物交叉点和非建筑物交叉点;
本实施例中,与二进制特征分类器相匹配的二进制混合描述符包含方向特征(orientation)和描述符(descriptor)两部分:
方向特征为:
θ=atan2(m01,m10),
其中,θ为角度特征,m01和m10为每个图片分割块的矩,计算公式如下:
mpq=∑x,yxpyqi(x,y),
描述符为:
fn(p)=∑0≤j<4∑0≤i<n2n×j+iτj(p;ai,bi),
其中,fn(p)表示分割块p的描述符,τ1为纹理特征测试对,τ2为相似度特征测试对,τ3和τ4为颜色信息测试对,计算公式分别如下:
其中,pl(a)和pl(b)分别表示l通道上以a为中心和以b为中心的分割块p,pa(a)和pa(b)分别表示a通道上以a为中心和以b为中心的分割块p,pb(a)和pb(b)分别表示b通道上以a为中心和以b为中心的分割块p,s为相似度阈值,ca为a通道的颜色阈值,cb为b通道的颜色阈值,在本发明中,ca和cb均设为255。
最终得到分割块p的描述符形式定义为
其中,n取256。
步骤6:将非建筑物交叉点所连接的图斑进行剔除,得到粗略提取的建筑物轮廓;
步骤7:设定最小建筑物区域像素阈值,当建筑物区域像素和小于阈值,则被认为是非建筑物区域,以此剔除错误检测的建筑物;
步骤8:获得建筑物轮廓的外接矩形,得到建筑物区域,并将其位置还原至原始正射图像上;
步骤9:在原始正射图像上,提取建筑物区域,并将其颜色空间由rgb转换为lab;
步骤10:在lab颜色空间下,对建筑物区域图像进行slic超分割处理,分割方法见步骤3;
步骤11:根据超分割结果标记出新交叉点;
步骤12:利用二进制特征分类器对新交叉点进行三分类,获得建筑物内部点、建筑物边缘点和非建筑物点,特征描述方法见步骤5;
步骤13:将非建筑物交叉点所连接的图斑进行剔除,得到精细提取的建筑物轮廓;
步骤14:设定最小建筑物区域像素阈值,剔除错误检测的建筑物;
步骤15:提取建筑物轮廓曲线,获得建筑物轮廓图像。
本例中,在对于建筑物轮廓进行提取时,分为两个阶段,首先进行对原始正射图像中的建筑物区域进行提取,然后在建筑物区域图像的基础之上进行建筑物轮廓的提取,从而经过由粗到精的建筑物提取过程,可准确得到最终的建筑物轮廓,大大地减少误检可能。
最后,为了便于对本发明所述的处理方法有更直观的认识,图2给出了本发明方法在各阶段的步骤图示,其中(a)是原始正射图像,(b)是经步骤3对降采样后正射图像进行slic超分割的结果,(c)是经步骤8得到的建筑物区域的结果,(d)是经步骤10对建筑物区域进行超分割后的结果,(e)是采用本发明方法提取的建筑物轮廓的最终结果。