CN114387329B - 基于高分遥感影像的建筑物轮廓递进式规则化方法 - Google Patents
基于高分遥感影像的建筑物轮廓递进式规则化方法 Download PDFInfo
- Publication number
- CN114387329B CN114387329B CN202210046818.1A CN202210046818A CN114387329B CN 114387329 B CN114387329 B CN 114387329B CN 202210046818 A CN202210046818 A CN 202210046818A CN 114387329 B CN114387329 B CN 114387329B
- Authority
- CN
- China
- Prior art keywords
- building
- contour
- corner
- outline
- points
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000000750 progressive effect Effects 0.000 title claims abstract description 10
- 238000001514 detection method Methods 0.000 claims abstract description 27
- 238000005457 optimization Methods 0.000 claims abstract description 20
- 238000012216 screening Methods 0.000 claims abstract description 16
- 230000001788 irregular Effects 0.000 claims abstract description 11
- 230000007246 mechanism Effects 0.000 claims abstract description 7
- 238000000605 extraction Methods 0.000 claims description 20
- 239000011159 matrix material Substances 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000005316 response function Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 238000013527 convolutional neural network Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000007430 reference method Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 244000291564 Allium cepa Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000012946 outsourcing Methods 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20112—Image segmentation details
- G06T2207/20164—Salient point detection; Corner detection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Mathematical Optimization (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Computational Linguistics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Molecular Biology (AREA)
- Evolutionary Biology (AREA)
- Geometry (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Health & Medical Sciences (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于高分遥感影像的建筑物轮廓递进式规则化方法，其利用改进的Harris角点检测算法对轮廓进行角点提取，再通过角点筛选机制剔除无用角点，顺序拟合保留的角点集，实现轮廓初步规整优化；然后利用基于Frechet距离的最小面积外接矩形算法对建筑物轮廓边界进行优化；得到整体规则、局部不规则的建筑物轮廓；最后再通过Shi‑Tomasi算法对局部不规整且呈锯齿状的建筑物轮廓边缘进行深度规则化。本发明的规则化总体精度达到了85.36％，相较于初始轮廓提高了13.17％。表明本方法适用于建筑物的轮廓规则化，有效提高了建筑物轮廓边缘的表达精度，能够准确地适应建筑物轮廓的细节变化。
Description
技术领域
本发明涉及建筑物轮廓规则化算法领域，具体涉及一种基于高分遥感影像的建筑物轮廓递进式规则化方法。
背景技术
基于高分辨率遥感影像的建筑物轮廓规则化方法是目前的热点研究方向之一，主要解决建筑物轮廓提取过程中存在的形状不规则、边缘锯齿状、不能准确的体现出建筑物轮廓形状和大小，无法保证提取的建筑物轮廓与原建筑物轮廓一致性等问题。所以，进一步将提取的建筑物轮廓整体规则化，使其与真实建筑物轮廓保持一致性，对地区的发展和规划有着重要意义。
当前轮廓规则化方法主要有以下三类：第一类是针对建筑物轮廓存在的纹理信息和空间形状特征的差异性，对图像形态特征的修复来实现建筑物轮廓进行规则化，如王伟玺提出一种基于栅格填充的规则化方法，利用轮廓的整体信息进行规则化，并用图像处理中的腐蚀、膨胀算法进行优化，使用图像像素二值化提取优化后的建筑物轮廓。黄金库针对建筑物的不同形态，分别采用手扶、自动跟踪数字化方法来规则化建筑物轮廓。此类方法以轮廓个体为单位进行整体运算，容易造成细节损失，同时会引起错误规则化的问题。第二类是针对建筑物轮廓边缘转折点的提取，拟合连接得到规则化的建筑物轮廓。郭珍珍[i]采用改进的管子算法确定轮廓线的关键点，拟合直线的交点作为新的关键点，连接各个关键点后构成的轮廓线能粗略呈现建筑物的形状，再用自适应的强制正交规则化算法得到规则化的轮廓线。该类方法得到的建筑物轮廓规则化效果较好，但是在一些局部区域仍然存在锯齿状情况，不适用于密集建筑物轮廓规则化。第三类方法是基于深度学习进行建筑物轮廓规则化。丁亚洲提出了一种基于多星形约束的图割和轮廓规则化的交互式半自动提取高分影像上直角建筑物的方法。Shiqing Wei提出了一个自动建筑足迹提取框架，该框架由基于卷积神经网络(CNN)的分割和经验多边形正则化组成，该正则化将分割地图转换为结构化的单个建筑多边形，试图用算法取代测绘领域中涉及的建筑足迹的部分手动描绘。黄小赛提出了一种基于卷积神经网络的集成方法，包括建筑物定位、形状判断、形状匹配等过程。但该类方法实现复杂、步骤繁多，且得到的建筑物轮廓存在较严重的锯齿状情况，不能取得理想的规则化效果。
发明内容
为解决上述问题，本发明提供了一种基于高分遥感影像的建筑物轮廓递进式规则化方法。
为实现上述目的，本发明采取的技术方案为：
基于高分遥感影像的建筑物轮廓递进式规则化方法，包括如下步骤：
S1、在提取的原始建筑物轮廓基础上，利用改进的Harris角点检测算法对轮廓进行角点提取，再通过角点筛选机制剔除无用角点，顺序拟合保留的角点集，实现轮廓初步规整优化；
S2、利用基于Frechet距离的最小面积外接矩形对拟合连接后的建筑物轮廓进行边缘线段优化，将建筑物轮廓线段与最小面积外接矩形线段进行离散等分，计算得出各个等分点对应的最短距离dmin，设置距离阈值δ，判断建筑物轮廓线段等分点坐标是否替换为最小面积外接矩形边界等分点坐标，再顺序拟合保留下来的离散等分点，得到初步规则化的建筑物轮廓；
S3、利用Shi-Tomasi算法对不规则的局部区域依次进行角点检测、筛选、拟合，进行深度规则化。
进一步地，所述步骤S1包括如下步骤：
S11、在提取的原始建筑物轮廓基础上，利用改进的Harris角点检测算法对轮廓进行角点提取；具体的：
首先计算出图像中的灰度变化E(u，v)：
E(u，v)＝∑w(x，y)[I(x+u，y+v)-I(x，y)]2 (1)
式中，(u，v)表示的是窗口偏移量，w(x，y)是移动的窗口函数，I(x+u，y+v)是平移后的图像灰度，I(x，y)是图像灰度；
I(x+u，y+v)＝I(x，y)+Ixu+Iyv+O(u2，v2) (2)
转化得到：
对于局部微小的窗口移动量[u，v]，可以近似得到：
其中，M是为梯度的协方差矩阵，由图像导数可得：
协方差矩阵M的特征值分析：
其中，λ1，λ2是M的两个特征值，由此得到定义角点响应函数CRF：
R＝detM-k[trace(M)]2 (8)
式中，detM＝λ1λ2，trace(M)＝λ1+λ2，k为经验常数，取值范围为[0.04，0.06]；
基于协方差矩阵M矩阵保存候选角点位置，初值设置为0，角点值设置为1，当角点(i，j)八邻域的“相似度”参数在中心点与领域其他八个点的像素值之差在(-t，+t)之间，确认它们为相似点，且相似点不在候选角点中；
S12、将检测得到的角点集合排序，利用角点筛选机制决定当前角点是否保留；
S13、剔除无关角点后，顺序拟合各个角点得到初始规则化后的建筑物轮廓。
本发明的规则化总体精度达到了85.36％，相较于初始轮廓提高了13.17％。表明本方法适用于建筑物的轮廓规则化，有效提高了建筑物轮廓边缘的表达精度，能够准确地适应建筑物轮廓的细节变化。
附图说明
图1(a)Harris角点检测算法；(b)改进的Harris角点检测算法
图2为角点间夹角示意图。
图3为初始轮廓边缘规则化；
图中：(a)初始建筑物轮廓提取结果；(b)初始轮廓的角点检测；(c)角点筛选拟合后的建筑物轮廓；(d)拟合后的轮廓外接矩形
图4为凸包坐标旋转示意图。
图5为使用不同外接矩形的多边形拟合轮廓结果；
图中：(a)多边形拟合效果；(b)最小边长外接矩形结果；(c)最小面积外接矩形结果；(d)最小外包矩形结果。
图6为Frechet距离规整轮廓流程示意图；
图中：(a)初始建筑物轮廓提取结果；(b)计算边界上各个等分离散点的欧氏距离；(c)离散等分点的最短距离集合；(d)初步规则化后效果。
图7为初步规则化过程；
图中：(a)初始矩形建筑物轮廓提取结果；(b)角点筛选拟合后结果；(c)矩形轮廓边界规则化；(d)矩形建筑物轮廓初步规则化结果。
图8：(a)初始非矩形建筑物轮廓提取结果；(b)非矩形轮廓初始规则化结果。
图9为深度规则化过程；
图中：(a)初始矩形建筑物轮廓提取结果；(b)初始规则化结果；(c)深度局部区域角点检测；(d)深度局部区域角点拟合效果。
图10为本发明提出的规则化方法实现过程；
图中：(a)建筑物轮廓真值；(b)提取的建筑物轮廓初始值；(c)角点检测；(d)角点拟合结果；(e)轮廓边界优化；(f)规则化结果。
图11为本发明实施例基于高分遥感影像的建筑物轮廓递进式规则化方法的流程图。
图12为规则化算法结果比较；
图中：(a)建筑物轮廓真值；(b)初始轮廓提取结果；(c)文献[2]方法的建筑物轮廓优化结果；(d)文献[1]方法的建筑物轮廓优化结果；(e)手动规则化建筑物轮廓结果；(f)本文算法规则化结果。
图13为规则化算法结果比较；
图中：(a)建筑物轮廓真值；(b)初始轮廓提取结果；(c)文献[18]方法的建筑物轮廓优化结果；(d)文献[17]方法的建筑物轮廓优化结果；(e)手动规则化建筑物轮廓结果；(f)本文方法的建筑物轮廓规则化结果。
具体实施方式
为了使本发明的目的及优点更加清楚明白，以下结合实施例对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明是实施例提供了一种基于高分遥感影像的建筑物轮廓递进式规则化方法，包括如下步骤：
S1、在提取的原始建筑物轮廓基础上，利用改进的Harris角点检测算法对轮廓进行角点提取，再通过角点筛选机制剔除无用角点，顺序拟合保留的角点集，实现轮廓初步规整优化；
S2、利用基于Frechet距离的最小面积外接矩形对拟合连接后的建筑物轮廓进行边缘线段优化，将建筑物轮廓线段与最小面积外接矩形线段进行离散等分，计算得出各个等分点对应的最短距离dmin，设置距离阈值δ，判断建筑物轮廓线段等分点坐标是否替换为最小面积外接矩形边界等分点坐标，再顺序拟合保留下来的离散等分点，得到初步规则化的建筑物轮廓；
S3、利用Shi-Tomasi算法对不规则的局部区域依次进行角点检测、筛选、拟合，进行深度规则化。
轮廓初步规则化
本实施例中，采用改进的的Harris算法对建筑物轮廓依次进行角点提取、剔除、拟合，使得建筑物轮廓边界清晰。改进的的Harris算法原理是寻找图像边缘曲线中曲率极大值的点，对于一个灰度图像I，将窗口w在I中移动，同时设置参数t为八邻域的“相似度”参数，当中心点与邻域其他八个点的像素值只差在(-t，+t)之间确认它们为相似角点，因为感兴趣的角点只出现在边界上，所以不用全部检测图像每个点，这提高了算法的时间复杂度，具体的，包括如下步骤：
计算出图像中的灰度变化E(u，v)：
E(u，v)＝∑w(x，y)[I(x+u，y+v)-I(x，y)]2 (1)
式中，(u，v)表示的是窗口偏移量，w(x，y)是移动的窗口函数，I(x+u，y+v)是平移后的图像灰度，I(x，y)是图像灰度；
I(x+u，y+v)＝I(x，y)+Ixu+Iyv+O(u2，v2) (2)
转化得到：
对于局部微小的窗口移动量[u，v]，可以近似得到：
其中，M是为梯度的协方差矩阵，由图像导数可得：
协方差矩阵M的特征值分析：
其中，λ1，λ2是M的两个特征值，由此得到定义角点响应函数CRF：
R＝detM-k[trace(M)]2 (8)
式中，detM＝λ1λ2，trace(M)＝λ1+λ2，k为经验常数，取值范围为[0.04，0.06]；
协方差矩阵M矩阵用来保存候选角点位置，初值设置为0，角点值设置为1，当角点(i，j)八邻域的“相似度”参数在中心点与领域其他八个点的像素值之差在(-t，+t)之间，确认它们为相似点，且相似点不在候选角点中。
而改进的Harris算法并没有全部检测图像的每个点，而是除去了边界上boundary个像素(最佳取值为4)，提高了轮廓的角点检测效率，角点检测结果如图1所示，通过对比可知，改进的Harris角点检测算法在检测角点个数和检测准确率上都要明显优于Harris算法，其中在检测轮廓效果大致相同的情况下，改进的Harris角点检测算法检测的无用角点和错误角点个数更少，算法耗时更少，且检测准确度提高了11.09％，如表1所示。
表1 Harris角点检测算法精度对比
本实施例中，将检测得到的角点集合排序，利用角点筛选机制决定当前角点是否保留。如图2所示，以Pi为起始点，计算Pi-1、Pi、Pi+1三点连接形成的夹角α，设置夹角角度阈值β(-75°，75°)决定角点是否保留。设Pi-1Pi线段斜率为k1，PiPi+1线段斜率为k2，则可以得出Pi-1PiPi+1三点之间的夹角α。若三点两线段的角度α绝对值小于β，则默认该Pi-1、Pi、Pi+1点在同一直线上，剔除Pi点；否则，该点为边缘转折点，保留Pi点。依次迭代各个建筑物轮廓提取的边缘点，剔除无关角点，保留剩下的角点集合。夹角计算公式为：
α＝arctan[(k1-k2)/(1+k1k2)] (9)
剔除无关角点后，顺序拟合各个角点得到初始规则化后的建筑物轮廓，如图3(c)，相较于初始建筑物轮廓，进行角点筛选拟合后的轮廓边缘更加清晰、规则，有利于接下来建筑物轮廓的边界优化以及深度规则化。
轮廓边缘优化
如图3(d)所示，建筑物轮廓边缘初步规则化后，轮廓边界和转角会存在一些缺失和不规整，而利用基于Frechet距离的最小面积外接矩形算法优化后，解决了建筑物轮廓边缘存在部分角落缺失和边界不规整等细节问题。Frechet距离能够准确地衡量建筑物拟合轮廓与最小面积外接矩形存在的边界距离差异，以最小面积外接矩形边界为参考边界，Q1为起始点，顺序计算外接矩形边界与建筑物轮廓边缘的离散等分点的欧式距离，设定距离阈值δ来决定外接矩形边界离散等分点的取舍，最后得到初步规则化的建筑物轮廓。
(1)最小面积外接矩形
1)选取所得凸包中一条边作为起始边，以该边左端点为中心旋转凸包，使得该边平行于X轴，如图4所示，计算其最小外接矩形，记录最小外接矩形的坐标和旋转角度。
2)依次选取其他边，按照步骤1)记录最小外接矩形的坐标和旋转角度。
3)比较所有最小外接矩形的面积，找出其中最小的外接矩形，按照其所对应的旋转角度，以该边的左端点为圆心顺时针旋转即为所求的最小面积外接矩形。
(2)离散Frechet距离算法
建筑物轮廓经常存在转折弯曲等局部变化特征，经多边形拟合出来的轮廓可以减少轮廓点的同时保留细节特征。为了更准确地保留和完善建筑物的轮廓特征，引入了离散Frechet距离算法，以准确地衡量建筑物的拟合轮廓与最小面积外接矩形存在的距离差异，作为评判拟合边界是否合适的标准。离散Frechet距离算法原理是基于最小面积外接矩形的轮廓为P且长度为N，初步规则化建筑物轮廓为Q且长度为M，如图所示，将P和Q的各个线段边界集合等分为n、m个等分离散点，两个轨迹各个等分点所对应的欧式距离即为边界等分点的相似度，通过对相似度与距离阈值δ的比较，决定建筑物边界上等分点是否保留。两者位置的描述可以用一个t变量的连续递增函数来刻画，用α(t)来表示最小面积外接矩形起始点位置描述函数，用β(t)表示初步规则化建筑物轮廓起始点位置描述函数。将变量t约束到区间[0，1]内，那么有α(0)＝0，α(1)＝N，β(0)＝0，β(1)＝M。P(α(t))和Q(β(t))分别表示t时刻两点在各自轮廓轨迹上的位置，两点之间的距离会随着α(t)和β(t)函数本身的不同和变量t的变化而不同。Frechet距离数学表达式如下：
δF(P，Q)＝minα[0，1]→[0，N]，β[0，1]→[0，M]{max∈[0，1]d(P(α(t))，Q(β(t)))}(10)
其中，d(α(t)，β(t))为整个过程中Pi点到Qi点在t时刻的欧氏距离，即轮廓边界点相似度。首先对建筑物边缘轮廓点和最小面积外接矩形轮廓点构成的线段集合依次顺时针排序，如图6(a)所示，建筑物轮廓线段集合为{P1P2，P2P3，…，Pn-1Pn}，最小面积外接矩形线段集合为{Q1Q2，Q2Q3，Q3Q4，Q4Q1}，依次判断建筑物轮廓各个线段集合与对应的矩形边缘的线段关系。例如，在Q1Q2段内，对应P1P2，P2P3，P3P4，P4P5，P5P6，P6P7六个线段，其中线段P2P3、P5P6垂直于线段Q1Q2，则线段P2P3、P5P6上的等分离散点不与矩形边缘等分离散点做距离计算，保留P2P3和P5P6线段位置；若线段不与Q1Q2呈垂直关系，则计算其等分离散点间的欧氏距离。如图6(b)所示，将P3P4、Q1Q2分别被划分为n、m个等分离散点，线段Q1Q2长度通常远大于对应的建筑物轮廓线段长度，所以n取值区间[25，30]，m取值区间[40，45]。如图6(b)所示，以建筑物轮廓线段P3P4上的等分离散点Wi为例，分别拟合连接Q1Q2线段上的Y1，Y2，…，Yj，得到最短距离H(Wi，Yj)作为该离散点与最小面积外接矩形边缘点的最短距离Hij，进而依次计算P3P4线段上所有离散等分点集合{W1，W2，…，Wi}与与最小面积外接矩形边缘所有离散等分点集合{Y1，Y2，…，Yj}的最短距离集合{d1，d2，…，dn}。距离阈值δ决定建筑物轮廓上的各个离散等分点是否保留，如图6(c)。距离阈值δ公式为：
δ＝(Sbuild/Srect)×(Lmin/2) (11)
式中：Sbuild为当前建筑物轮廓的面积；Srect为基于最小面积外接矩形的面积；Lmin为基于最小面积外接矩形的最短边长度。当建筑物存在内凹或者复杂拐角时，如果直接使用Lmin的值，δ的值就会偏大，无法保留拐角细节，甚至会破坏建筑物轮廓形状，因此使用Lmin/2更符合优化要求，不至于破坏建筑物轮廓形状，而Sbuild与Srect的比值越大则说明建筑物真实轮廓越接近矩形。若两点距离di＜δ，则视为该轮廓线段与建筑物基于最小面积外接矩形的轮廓相似度较高，此时将建筑物轮廓等分离散点Wi的坐标替换为其最小面积外接矩形上最短距离对应等分离散点Yi的坐标。若欧氏距离di＞δ，则将欧氏距离大于δ的建筑物轮廓等分点Wi的坐标保留。依次完成距离阈值δ控制下的各项等分离散点的规整，进一步规则化建筑物轮廓。
如图7所示，初始轮廓优化方法直接适用于规则化边缘轮廓为矩形的所有轮廓建筑物，且整体性与建筑物真值基本保持一致。但是大多数建筑物轮廓为不规则的、非矩形的，如图8所示，若存在非矩形建筑物轮廓，则初始规则化方法无法完整的优化局部锯齿状边缘，需要进一步深度规则化其局部轮廓。
局部轮廓深度规则化
(1)局部锯齿状区域的角点检测：
如图9(b)，在对其他复杂多边形形状的建筑物轮廓提取中，初步规则化后的建筑物轮廓边界仍然会出现锯齿状或者区域性不规则，需要对残存的多齿状边缘深度规则化，要求局部小范围角点检测且时间复杂度要求不高，用Shi-Tomasi算法对区域性的多齿状边缘进行角点检测、剔除以及拟合，得到完整规则化的建筑物轮廓。
(2)锯齿状区域角点的筛选、拟合：
将角点进行排序，利用角度阈值β对角点筛选，删除无用的角点，最后顺序拟合保留的角点，得到完整规则化的建筑物轮廓。为验证该方法的有效性，选取乌鲁木齐天山区建筑物轮廓进行检验，步骤实现如图11所示。
验证与分析
实用性验证
本发明采用改进的Harris算法对初始建筑物轮廓进行角点检测、筛选、拟合；然后利用基于Frechet距离的最小面积外接矩形算法对建筑物轮廓边界进行优化；得到整体规则、局部不规则的建筑物轮廓；最后再通过Shi-Tomasi算法对局部不规整且呈锯齿状的建筑物轮廓边缘进行深度规则化，结果与建筑物原始轮廓基本一致。为了更清晰地表示本发明方法与其他方法的对比结果，选取新疆乌鲁木齐市天山区为数据原始影像以及初始轮廓提取结果，以文献[1](Fischler M A，Bolles R C.Random sample consensus：a paradigmfor model fitting with applications to image analysis and automatedcartography[J].Communications of the ACM，1981，24(6)：381-395.)、文献[2](王杰茜，冯德俊，陈建飞.对比Harris算子和Susan算子的建筑物边界规则化方法[J].测绘通报，2020(04)：11-15.[19]Sampath A，Shan J.Building boundary tracing andregularization from airborne LiDAR point clouds[J].PhotogrammetricEngineering&Remote Sensing，2007，73(7)：805-812.)的方法和手动规则化为参照方法，文献[1]是通过引入卷积层特征金字塔的多尺度聚合，建立了一种尺度鲁棒的全卷积网络(FCN)。采用两种后处理策略对FCN分割图进行细化，对细化后的分割图进行矢量化和多边形化。且提出一种由粗调整和细调整组成的多边形正则化算法[19]，将初始多边形转换为结构化足迹。该多边形正则化算法在不同建筑风格、图像分辨率甚至低质量分割的挑战性情况下具有鲁棒性。文献[2]对粗提取并预处理后的建筑物边界采用Harris算子进行角点检测、排序，然后通过算法对其进行规则化边界拟合处理得到接近建筑物实际边界的规则化边界。该方法的边界拟合效果较为依赖建筑物边界角点检测结果，且规则化后的结果会出现较多的边缘毛刺突起。手动绘制建筑轮廓方法是使用arcgis对导入的建筑物轮廓提取的初始影像进行描点连线绘制，以起点为绘制点，直线连接轮廓形状的各个转折点，最后得到趋近于建筑物轮廓真值的结果，虽然该方法规则化的建筑物轮廓结果趋近于建筑物轮廓真实值，但由于手动绘制的限制性因素，该方法仅局限于对局部建筑物轮廓的规则化，不能应用于大范围的轮廓规则化，没有实际应用性。本发明方法利用基于Frechet距离的最小面积外接矩形和对遗漏和缺失的建筑物局部轮廓信息就行了较好的还原和补充，该方法规则化结果趋近于手动绘制的结果，优于文献[1]方法以及文献[2]方法结果，体现出了本发明方法的有效性和高效性，如图12和图13所示。
2.2对比分析
为了准确的阐述本发明算法的精确性和高效性，利用二分类评价体系中的完整度(CM)、正确率(CR)、综合值(F1)以及总体精度(0A)，从像素的角度对提取结果进行精度评定。选取图13中建筑物轮廓规则化结果进行论证比较，表2给出三种建筑物轮廓规则化的精度结果，以看到本发明方法相较于其他两类方法，轮廓规则化后精度均有较大提升。本发明方法相较于初始轮廓，综合值和总体精度分别提高了8.18％和13.17％，相较于文献[1]和文献[2]两种优化方法，综合值分别提高了1.43％和3.14％，总体精度分别提高了5.65％和7.77％。
表2.不同建筑物轮廓规则化方法精度对比
本发明方法利用基于最小面积外接矩形的规整操作有效地还原了缺失的部分，使得轮廓优化结果更接近原始建筑物形状。此外，针对复杂轮廓的局部细节优化问题，本发明方法提取角点后依据角度阈值选取角点，保留了建筑物轮廓的细节部分，使复杂建筑物整体形状更加规整，能有效应用于影像建筑物形状排列复杂、周围地物干扰多的场景。总之，本发明方法深度改善了建筑物结果的规整性，综合值和总体精度均优于初始提取结果，通过与两种轮廓优化参照方法的对比，说明本发明通过递进式规则化，取得了明显的效果，进一步提高了建筑物轮廓的表达精度。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以作出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.基于高分遥感影像的建筑物轮廓递进式规则化方法，其特征在于：包括如下步骤：
S1、在提取的原始建筑物轮廓基础上，利用改进的Harris角点检测算法对轮廓进行角点提取，再通过角点筛选机制剔除无用角点，顺序拟合保留的角点集，实现轮廓初步规整优化；
S2、利用基于Frechet距离的最小面积外接矩形对拟合连接后的建筑物轮廓进行边缘线段优化，将建筑物轮廓线段与最小面积外接矩形线段进行离散等分，计算得出各个等分点对应的最短距离dmin，设置距离阈值δ，判断建筑物轮廓线段等分点坐标是否替换为最小面积外接矩形边界等分点坐标，再顺序拟合保留下来的离散等分点，得到初步规则化的建筑物轮廓；
S3、利用Shi-Tomasi算法对不规则的局部区域依次进行角点检测、筛选、拟合，进行深度规则化；
所述步骤S1包括如下步骤：
S11、在提取的原始建筑物轮廓基础上，利用改进的Harris角点检测算法对轮廓进行角点提取；具体的：
首先计算出图像中的灰度变化E(u，v)：
E(u，v)＝∑w(x，y)[I(x+u，y+v)-I(x，y)]2 (1)
式中，(u，v)表示的是窗口偏移量，w(x，y)是移动的窗口函数，I(x+u，y+v)是平移后的图像灰度，I(x，y)是图像灰度；
I(x+u，y+v)＝I(x，y)+Ixu+Iyv+O(u2，v2) (2)
转化得到：
对于局部微小的窗口移动量[u，v]，可以近似得到：
其中，M是为梯度的协方差矩阵，由图像导数可得：
协方差矩阵M的特征值分析：
其中，λ1，λ2是M的两个特征值，由此得到定义角点响应函数CRF：
R＝detM-k[trace(M)]2 (8)
式中，detM＝λ1λ2，trace(M)＝λ1+λ2，k为经验常数，取值范围为[0.04，0.06]；
基于协方差矩阵M矩阵保存候选角点位置，初值设置为0，角点值设置为1，当角点(i，j)八邻域的“相似度”参数在中心点与领域其他八个点的像素值之差在(-t，+t)之间，确认它们为相似点，且相似点不在候选角点中；
S12、将检测得到的角点集合排序，利用角点筛选机制决定当前角点是否保留；
S13、剔除无关角点后，顺序拟合各个角点得到初始规则化后的建筑物轮廓。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210046818.1A CN114387329B (zh) | 2022-01-16 | 2022-01-16 | 基于高分遥感影像的建筑物轮廓递进式规则化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210046818.1A CN114387329B (zh) | 2022-01-16 | 2022-01-16 | 基于高分遥感影像的建筑物轮廓递进式规则化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114387329A CN114387329A (zh) | 2022-04-22 |
| CN114387329B true CN114387329B (zh) | 2024-03-22 |
Family
ID=81202227
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210046818.1A Active CN114387329B (zh) | 2022-01-16 | 2022-01-16 | 基于高分遥感影像的建筑物轮廓递进式规则化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114387329B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114898119B (zh) * | 2022-07-08 | 2022-11-01 | 浙江大华技术股份有限公司 | 一种建筑物轮廓绘制方法、装置、设备及介质 |
| CN117115161B (zh) * | 2023-10-24 | 2024-01-02 | 四川新康意众申新材料有限公司 | 一种塑料缺陷检查方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2787856A1 (en) * | 2012-05-12 | 2013-11-12 | University Of Florida Research Foundation, Inc. | Systems and methods for estimating the geographic location at which image data was captured |
| CN109903304A (zh) * | 2019-02-25 | 2019-06-18 | 武汉大学 | 一种基于卷积神经元网络和多边形规则化的建筑物轮廓自动提取算法 |
| CN112489185A (zh) * | 2019-08-20 | 2021-03-12 | 黎欧思照明(上海)有限公司 | 一种基于空间数据采集的集成灯光建模方法 |
| CN112487537A (zh) * | 2020-12-08 | 2021-03-12 | 亿景智联(北京)科技有限公司 | 一种基于全卷积神经网络的建筑面多级优化提取方法 |
-
2022
- 2022-01-16 CN CN202210046818.1A patent/CN114387329B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2787856A1 (en) * | 2012-05-12 | 2013-11-12 | University Of Florida Research Foundation, Inc. | Systems and methods for estimating the geographic location at which image data was captured |
| CN109903304A (zh) * | 2019-02-25 | 2019-06-18 | 武汉大学 | 一种基于卷积神经元网络和多边形规则化的建筑物轮廓自动提取算法 |
| CN112489185A (zh) * | 2019-08-20 | 2021-03-12 | 黎欧思照明(上海)有限公司 | 一种基于空间数据采集的集成灯光建模方法 |
| CN112487537A (zh) * | 2020-12-08 | 2021-03-12 | 亿景智联(北京)科技有限公司 | 一种基于全卷积神经网络的建筑面多级优化提取方法 |
Non-Patent Citations (2)
| Title |
|---|
| "双激光短程机动车测速仪研究";马龙;《工程科技Ⅱ辑》;第1-6章 * |
| "高分遥感影像建筑物轮廓的逐级优化方法";常京新 等;《中国激光》;第47卷(第10期);第１０１０００２－１至１０１０００２－１４页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114387329A (zh) | 2022-04-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110443836B (zh) | 一种基于平面特征的点云数据自动配准方法及装置 | |
| CN114418957B (zh) | 基于机器人视觉的全局和局部二值模式图像裂缝分割方法 | |
| CN114387329B (zh) | 基于高分遥感影像的建筑物轮廓递进式规则化方法 | |
| US9123112B2 (en) | Method for the pre-processing of a three-dimensional image of the surface of a tyre for use in the inspection of said surface | |
| Xu et al. | Reconstruction of scaffolds from a photogrammetric point cloud of construction sites using a novel 3D local feature descriptor | |
| CN110222642B (zh) | 一种基于全局图聚类的平面建筑构件点云轮廓提取方法 | |
| CN109685080B (zh) | 基于霍夫变换与区域生长的多尺度平面提取方法 | |
| CN106981077B (zh) | 基于dce和lss的红外图像和可见光图像配准方法 | |
| CN110910387B (zh) | 一种基于显著性分析的点云建筑物立面窗户提取方法 | |
| CN101976504B (zh) | 一种基于颜色空间信息的多车辆视频跟踪方法 | |
| CN111508073B (zh) | 一种三维建筑模型屋顶轮廓线的提取方法 | |
| CN112164145B (zh) | 基于点云数据快速提取室内三维线段结构的方法 | |
| CN111860501B (zh) | 基于形状匹配的高铁高度调整杆脱出故障图像识别方法 | |
| CN103727930A (zh) | 一种基于边缘匹配的激光测距仪与相机相对位姿标定方法 | |
| CN110428425B (zh) | 一种基于海岸线矢量数据的sar图像海陆分离方法 | |
| CN106096497B (zh) | 一种针对多元遥感数据的房屋矢量化方法 | |
| CN116740072B (zh) | 基于机器视觉的道路表面缺陷检测方法及系统 | |
| CN111783722B (zh) | 一种激光点云的车道线提取方法和电子设备 | |
| CN107610174B (zh) | 一种鲁棒的基于深度信息的平面检测方法及系统 | |
| Wang | Automatic extraction of building outline from high resolution aerial imagery | |
| CN112489207B (zh) | 一种空间约束的密集匹配点云平面基元提取方法 | |
| CN102679871B (zh) | 亚像素精度工业物体快速检测方法 | |
| CN114937211B (zh) | 基于不变对象的遥感影像变化检测方法 | |
| Lin et al. | Unsupervised window extraction from photogrammetric point clouds with thermal attributes | |
| Kang et al. | Line feature extraction from RGB laser point cloud |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |