CN115880575B - 结合变化信息和建筑特征的遥感影像新增建筑提取方法 - Google Patents
结合变化信息和建筑特征的遥感影像新增建筑提取方法 Download PDFInfo
- Publication number
- CN115880575B CN115880575B CN202211314630.7A CN202211314630A CN115880575B CN 115880575 B CN115880575 B CN 115880575B CN 202211314630 A CN202211314630 A CN 202211314630A CN 115880575 B CN115880575 B CN 115880575B
- Authority
- CN
- China
- Prior art keywords
- building
- remote sensing
- newly added
- sensing image
- change information
- 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
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种结合变化信息和建筑特征的遥感影像新增建筑提取方法,属于遥感图像处理领域。其包括:分割后时相遥感影像获得超像素图;提取线段得到线段集合;计算超像素的边缘线段垂直性;计算建筑物线段指数;提取影像基于对象的形态学特征;计算建筑物强度;计算前后时相遥感影像的纹理特征;计算前后时相遥感影像的纹理特征差和光谱特征差;求纹理和光谱特征的变化强度,获得变化强度;获取新增建筑对象斑块;采用形状指数进行后处理,获取最终的新增建筑对象。本发明充分利用光谱和纹理特征提取变化信息,并通过边缘结构特征和形态学特征提取建筑特性,为大面积遥感影像实现新增建筑识别提供了可靠的技术支撑。
Description
技术领域
本发明属于遥感图像处理技术领域,具体涉及一种结合变化信息和建筑特征的遥感影像新增建筑提取方法。
背景技术
随着国民经济的发展和人口增加,我国工业化和城市化水平逐步提高,建设用地不断增加,耕地急剧减少,因此,及时掌握建设用地、耕地的更新状况,尤其是违法建筑物的情况,成为保障社会经济协调持续发展的重要任务。近年来,高分辨率遥感发展迅速,由于其可以提供地物更加详细的光谱和纹理信息,为建筑物等小尺度目标的变化检测提供了可能。传统的变化检测方法应用于高分辨率遥感影像中新增建筑物检测面临以下两点问题:一方面,基于像素级光谱信息的变化检测不能够很好的获取高分影像中地物的纹理和几何信息,提取精度较低;另一方面,高分影像中丰富的光谱和纹理信息为提取建筑物带来一定的困难,表现为其光谱与道路、裸土等地物相似性较高,结果中噪声较多。
虽然许多学者对高分辨率影像中建筑提取做了研究,但是这些研究使用的特征较为单一,无法很好体现建筑物特征及其变化信息;其次,一部分研究采用监督分类的方法,方法中面临样本难选、人工干预程度高的问题,给大面积、自动化、快速的提取变化建筑带来困难。
发明内容
本发明的目的在于提供一种结合变化信息和建筑特征的遥感影像新增建筑提取方法,其充分利用光谱和纹理特征提取变化信息,并通过边缘结构特征和形态学特征提取建筑特性,为大面积遥感影像实现新增建筑识别提供可靠的技术支撑。
本发明采用的技术方案为:
一种结合变化信息和建筑特征的遥感影像新增建筑提取方法,包括以下步骤:
步骤8,分别计算前后时相遥感影像的纹理特征差DT和光谱特征差DS;
步骤9,分别求纹理和光谱特征的变化强度IT和IS,通过直接求和获得变化强度CI;
步骤11,计算新增建筑强度图NBI的分割阈值T,对NBI二值化后获取新增建筑对象斑块NB;
步骤12,对每一新增建筑对象斑块NB分别计算形状指数GI,过滤出GI超过阈值的斑块NB,得到最终提取的新增建筑对象NBpost。
式中,F(Δβjk)表示当两条线段lj、lk与超像素Si相交时的垂直性值,Δβjk表示线段之间的夹角,θ表示倾斜角差容许误差。
进一步地,步骤4中建筑物线段指数BLISi采用如下公式计算:
其中,N()表示线段包含的像素个数,N(Si∩lm)表示线段lm与超像素Si相交的像素个数。
式中,MBI(x,y)表示形态学建筑指数MBI图的(x,y)处的象元值,p和q分别表示与超像素Si等面积的矩形的长和宽。
进一步地,步骤7中,采用GLCM的方差测度作为纹理特征,其计算方法如下:
其中,d∈{1,2},P(i,j)代表灰度级别i和j在给定方向和距离情况下的出现的可能性,μ是GLCM的均值,N代表灰度级。
进一步地,步骤8中,纹理特征差DT和光谱特征差DS的计算方法如下:
式中,S表示光谱波段,T1和T2分别表示遥感影像的两个时间,abs()表示求绝对值。
进一步地,步骤9中,纹理和光谱特征的变化强度IT和IS的计算方式如下:
其中,上标T表示转置运算,XDT表示纹理特征差DT的特征矩阵,行为像素数,列为波段数,XDS表示纹理特征差DS的特征矩阵,行为像素数,列为波段数;
变化强度CI采用如下方式计算:
CI=Nor(IT)+Nor(IS)
其中,Nor()表示对变化特征进行最大值最小值归一化处理。
进一步地,步骤10中,依据分割结果IMseg求对象级的变化强度CIseg的计算方式如下:
新增建筑强度图NBI的计算方式如下:
其中,Nor()表示对变化信息和建筑特征进行最大值最小值归一化处理。
进一步地,步骤11中,分割阈值T的计算方式如下:
T=ms+a×std
其中,ms和std分别表示统计NBI中所有对象的均值和标准差,a为调节因子;
对NBI二值化后获取新增建筑对象斑块NB,具体方式如下:
式中,T是基于所有对象隶属度确定的阈值,NB表示新增建筑物提取结果,像元灰度值为1表示新增建筑物,像元值为0表示背景值;
步骤12中,形状指数GI的计算方法为:
其中,GI表示形状指数,是基于连通域计算的特征,P表示该连通域的周长,A表示连通域的面积;
最终提取的新增建筑对象NBpost的计算方式如下:
TGI=a×std(GI)+mean(GI)
其中,GIj表示针对第j个斑块NB计算的形状指数,TGI表示形状指数阈值,像元灰度值为1表示新增建筑物,像元值为0表示背景值;a表示调节系数,std()表示求图像灰度值的均方差,mean()表示求图像灰度值的均值。
本发明具有如下有益效果:
(1)本发明提出了基于变化强度和后时相建筑特征的新增建筑提取方法,在不需要训练样本的情况下,可以有效提取多时相遥感图像中的新增建筑。
(2)应用本发明方法可以有效提取融合纹理和光谱信息的变化信息,为其他土地类型的变化检测提供关键的技术支撑。
总之,本发明采用了全新的非监督识别方法,结合了光谱、纹理、光谱空间和边缘结构特征提取建筑特征及其变化信息,有效地实现了像素级变化信息和对象级目标信息的融合,为自动高效的新增建筑提取提供重要的技术途径。
附图说明
图1为新增建筑提取方法的原理示意图。
具体实施方式
下面结合附图及具体实施方式对本发明做进一步详细的说明。
一种结合变化信息和建筑特征的遥感影像新增建筑提取方法,包括以下步骤:
Lset={l1,l2,l3...lm,m=1,...,n}
式中,F(Δβjk)表示当两条线段lj、lk与超像素Si相交时其垂直性值,Δβjk表示线段之间的夹角,θ表示倾斜角差容许误差。
其中,N()表示线段包含的像素个数,N(Si∩lm)表示lm与超像素Si相交的像素个数。
式中,MBI(x,y)表示MBI图的(x,y)处的象元值,p和q分别表示与超像素Si等面积的矩形的长和宽。
其中,d∈{1,2},P(i,j)代表灰度级别i和j在给定方向和距离情况下的出现的可能性,μ是GLCM的均值,N代表灰度级。
步骤8,分别计算前后时相遥感影像的纹理特征差DT和光谱特征差DS,方法如下:
式中,S表示光谱波段,abs()表示求加绝对值;
步骤9,分别求纹理和光谱特征的变化强度IT和IS,其计算方式如下:
其中,*T表示转置运算,XDT表示纹理特征差DT的特征矩阵,行为像素数,列为波段数,XDS表示纹理特征差DS的特征矩阵,行为像素数,列为波段数;
通过对归一化后的IT和IS直接求和,获得变化强度CI,
CI=Nor(IT)+Nor(IS)
其中,Nor()表示对变化特征归一化处理。
步骤10,依据分割结果IMseg求对象级的变化强度CIseg的计算方式如下:
式中,CI(x,y)表示CI图中(x,y)处的象元值,pc和qc分别表示与超像素Si等面积的矩形的长和宽。
其中,Nor()表示对变化信息和建筑特征进行最大值最小值归一化处理;
步骤11,计算新增建筑强度图NBI的分割阈值T,其计算方式如下:
T=ms+a×std
其中,ms和std分别统计NBI中所有对象的均值和标准差,a为调节因子;
然后对NBI二值化后获取新增建筑对象斑块NB,判断条件如下:
式中,T是基于所有对象隶属度确定的阈值,NB表示最终的新增建筑物提取结果,像元灰度值为1表示新增建筑物,像元值为0表示背景值。
步骤12,采用形状指数GI进行后处理,可以过滤道路等长宽比较大的斑块,形状指数GI计算方式如下:
其中,GI表示形状指数,是基于连通域计算的特征,P表示该连通域的周长,A表示连通域的面积;
最后提取的新增建筑对象NBpost,计算方式如下:
其中,GIj表示第j个连通域计算的形状指数,TGI表示形状指数阈值。
总之,本发明充分利用光谱和纹理特征提取变化信息,并通过边缘结构特征和形态学特征提取建筑特性,可以有效提取高分辨率遥感影像中的新增建筑物。该方法无需选择训练样本,无需人工干预识别,为便捷、快速、精确和大范围的新增建筑物提取提供了重要的技术支撑。
Claims (10)
1.一种结合变化信息和建筑特征的遥感影像新增建筑提取方法,其特征在于,包括以下步骤:
步骤8,分别计算前后时相遥感影像的纹理特征差DT和光谱特征差DS;
步骤9,分别求纹理和光谱特征的变化强度IT和IS,通过直接求和获得变化强度CI;
步骤11,计算新增建筑强度图NBI的分割阈值T,对NBI二值化后获取新增建筑对象斑块NB;
步骤12,对每一新增建筑对象斑块NB分别计算形状指数GI,过滤出GI超过阈值的斑块NB,得到最终提取的新增建筑对象NBpost。
10.根据权利要求9所述的结合变化信息和建筑特征的遥感影像新增建筑提取方法,其特征在于,步骤11中,分割阈值T的计算方式如下:
T=ms+a×std
其中,ms和std分别表示统计NBI中所有对象的均值和标准差,a为调节因子;
对NBI二值化后获取新增建筑对象斑块NB,具体方式如下:
式中,T是基于所有对象隶属度确定的阈值,NB表示新增建筑物提取结果,像元灰度值为1表示新增建筑物,像元值为0表示背景值;
步骤12中,形状指数GI的计算方法为:
其中,GI表示形状指数,是基于连通域计算的特征,P表示该连通域的周长,A表示连通域的面积;
最终提取的新增建筑对象NBpost的计算方式如下:
TGI=a×std(GI)+mean(GI)
其中,GIj表示针对第j个斑块NB计算的形状指数,TGI表示形状指数阈值,像元灰度值为1表示新增建筑物,像元值为0表示背景值;a表示调节系数,std()表示求图像灰度值的均方差,mean()表示求图像灰度值的均值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211314630.7A CN115880575B (zh) | 2022-10-26 | 2022-10-26 | 结合变化信息和建筑特征的遥感影像新增建筑提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211314630.7A CN115880575B (zh) | 2022-10-26 | 2022-10-26 | 结合变化信息和建筑特征的遥感影像新增建筑提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115880575A CN115880575A (zh) | 2023-03-31 |
CN115880575B true CN115880575B (zh) | 2023-05-16 |
Family
ID=85758935
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211314630.7A Active CN115880575B (zh) | 2022-10-26 | 2022-10-26 | 结合变化信息和建筑特征的遥感影像新增建筑提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115880575B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111339948A (zh) * | 2020-02-25 | 2020-06-26 | 武汉大学 | 一种高分辨率遥感影像新增建筑自动识别方法 |
CN113191213A (zh) * | 2021-04-12 | 2021-07-30 | 桂林电子科技大学 | 一种高分辨率遥感影像新增建筑物检测方法 |
CN115131676A (zh) * | 2022-06-28 | 2022-09-30 | 南京信息工程大学 | 一种高分遥感影像建筑物变化的检测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8548248B2 (en) * | 2011-03-01 | 2013-10-01 | Mda Information Systems, Inc. | Correlated land change system and method |
WO2020164092A1 (zh) * | 2019-02-15 | 2020-08-20 | 深圳市大疆创新科技有限公司 | 图像处理方法、设备、可移动平台、无人机及存储介质 |
-
2022
- 2022-10-26 CN CN202211314630.7A patent/CN115880575B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111339948A (zh) * | 2020-02-25 | 2020-06-26 | 武汉大学 | 一种高分辨率遥感影像新增建筑自动识别方法 |
CN113191213A (zh) * | 2021-04-12 | 2021-07-30 | 桂林电子科技大学 | 一种高分辨率遥感影像新增建筑物检测方法 |
CN115131676A (zh) * | 2022-06-28 | 2022-09-30 | 南京信息工程大学 | 一种高分遥感影像建筑物变化的检测方法 |
Non-Patent Citations (2)
Title |
---|
《利用结构面对象的高分遥感影像建筑物变化检测》;刘海飞 等;《测绘通报》;全文 * |
Sejung Jung等.《Change Detection of Building Objects in High-Resolution Single-Sensor and Multi-Sensor Imagery Considering the Sun and Sensor's Elevation and Azimuth Angles》.《remote sensing》.2021,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN115880575A (zh) | 2023-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Three-dimensional convolutional neural network model for tree species classification using airborne hyperspectral images | |
CN110544251B (zh) | 基于多迁移学习模型融合的大坝裂缝检测方法 | |
Yang et al. | Hyper-temporal remote sensing data in bare soil period and terrain attributes for digital soil mapping in the Black soil regions of China | |
CN103559500B (zh) | 一种基于光谱与纹理特征的多光谱遥感图像地物分类方法 | |
CN110390267A (zh) | 一种基于高分遥感影像的山地景观建筑提取方法和设备 | |
CN103578110A (zh) | 基于灰度共生矩阵的多波段高分辨率遥感影像分割方法 | |
Duan et al. | Using remote sensing to identify soil types based on multiscale image texture features | |
CN107564017B (zh) | 一种城市高分遥感影像阴影检测及分割方法 | |
CN103971115A (zh) | 一种基于NDVI和PanTex指数的高分辨率遥感影像新增建设用地图斑自动提取方法 | |
CN105608458A (zh) | 一种高分辨率遥感影像建筑物提取方法 | |
Cheng et al. | Multi-scale Feature Fusion and Transformer Network for urban green space segmentation from high-resolution remote sensing images | |
CN105512622B (zh) | 一种基于图分割和监督学习的可见光遥感图像海陆分割方法 | |
CN106228130A (zh) | 基于模糊自编码网络的遥感影像云检测方法 | |
CN117115652A (zh) | 一种基于可见光谱视觉遥感的林地退化目标检测识别方法 | |
CN111882573A (zh) | 一种基于高分辨率影像数据的耕地地块提取方法及系统 | |
CN114241333A (zh) | 一种基于多源时序遥感影像精准识别新生滑坡区域的方法 | |
CN116682024A (zh) | 一种基于四波段遥感影像的快速云检测方法 | |
CN107657246B (zh) | 一种基于多尺度滤波建筑指数的遥感影像建筑物检测方法 | |
CN117292272A (zh) | 结合高分遥感影像和开源dem识别大范围耕作措施的方法 | |
CN101894368B (zh) | 基于测地距离的半监督sar图像水域分割方法 | |
Guo et al. | Mapping and classification of the liaohe estuary wetland based on the combination of object-oriented and temporal features | |
CN115880575B (zh) | 结合变化信息和建筑特征的遥感影像新增建筑提取方法 | |
Du et al. | A classification method of building structures based on multi-feature fusion of UAV remote sensing images | |
CN111476182A (zh) | 一种基于多源数据和多向纹理滤波分析的梯田提取方法 | |
CN115512159A (zh) | 面向对象的高分辨率遥感影像地表覆盖分类方法及系统 |
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 |