CN108805057B - 一种基于联合显著性分析的sar图像油库区检测方法 - Google Patents
一种基于联合显著性分析的sar图像油库区检测方法 Download PDFInfo
- Publication number
- CN108805057B CN108805057B CN201810530745.7A CN201810530745A CN108805057B CN 108805057 B CN108805057 B CN 108805057B CN 201810530745 A CN201810530745 A CN 201810530745A CN 108805057 B CN108805057 B CN 108805057B
- Authority
- CN
- China
- Prior art keywords
- image
- group
- sar
- images
- map
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
- G06F18/232—Non-hierarchical techniques
- G06F18/2321—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions
- G06F18/23213—Non-hierarchical techniques using statistics or function optimisation, e.g. modelling of probability density functions with fixed number of clusters, e.g. K-means clustering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/253—Fusion techniques of extracted features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/20—Image preprocessing
- G06V10/30—Noise filtering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/46—Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
- G06V10/462—Salient features, e.g. scale invariant feature transforms [SIFT]
- G06V10/464—Salient features, e.g. scale invariant feature transforms [SIFT] using a plurality of salient features, e.g. bag-of-words [BoW] representations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/50—Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/46—Descriptors for shape, contour or point-related descriptors, e.g. scale invariant feature transform [SIFT] or bags of words [BoW]; Salient regional features
- G06V10/467—Encoded features or binary features, e.g. local binary patterns [LBP]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V2201/00—Indexing scheme relating to image or video recognition or understanding
- G06V2201/07—Target detection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Probability & Statistics with Applications (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开一种基于联合显著性分析的SAR图像油库区检测方法,属于SAR图像处理及图像识别技术领域。实施过程包括:1)利用增强的方向平滑滤波器对一组SAR图像进行降噪处理;2)提取一组SAR图像中的亮度特征、纹理特征和曲线特征;3)利用模糊C均值聚类获得簇;4)利用全局对比度计算各个簇的显著值,得到共性显著特征图;5)计算每幅图像的共生直方图,通过分析共生直方图,计算得到单图显著特征图;6)融合共性显著特征图和单图显著特征图,得到油库区显著图;7)通过最大类间方差法进行阈值分割提取油库区。与传统的方法相比,本发明实现了对SAR图像油库区的准确提取,可应用于港区建设、环境监测及石油储量分析等方面。
Description
技术领域
本发明属于遥感图像处理及图像识别技术领域,具体涉及一种基于联合显著性分析的SAR图像油库区检测方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)不受天气、光照等条件的限制,具有无可比拟的全天候、全天时等诸多优点,因此被广泛应用于船只检测、地形勘察、资源勘探、目标识别以及环境监测等各个方面。SAR***中雷达技术的不断发展使得SAR图像的分辨率日益提高,数据量不断增加,因此传统的SAR图像解译方法已经不能满足海量数据的处理要求。如何从复杂的地物场景中有效检测出SAR图像目标,如:车辆、舰船、坦克、油库等,目前已成为遥感图像处理及图像识别技术领域的研究热点。
最早开展SAR图像目标识别技术研究的是美国麻省理工学院的林肯实验室。林肯实验室的SAIP***,使用的是经典的模板匹配方法。基于模板匹配的方法简单易实现,但由于SAR图像特征常常受各种因素的影响发生显著的变化,使得模板匹配方法在面对较为复杂的目标或有较强背景杂波干扰时常常发生漏检与误检等问题,已无法满足人们的需要。
一些专家学者借鉴了光学遥感影像的目标检测方法,提出了SAR图像的目标检测模型。但是SAR图像成像机制与光学遥感影像不同,如果直接使用现有的模型来处理SAR图像,将存在以下几个主要问题:(1)由于SAR独特的相干成像原理,SAR图像中存在大量的斑点噪声,这和光学遥感影像中的加性噪声不同。斑点噪声严重地影响了SAR图像的质量,并且也给后续的目标检测带来不利影响。(2)SAR数据成像时对目标的方位角变化很敏感,方位角不同会导致成像结果产生很大的差异,SAR图像中的目标信息也会产生变化。因此光学遥感影像中常用的边界信息无法在SAR图像中发挥作用。(3)与光学遥感影像不同,SAR图像缺少光谱信息,因此一些常用的光学遥感影像特征无法应用于SAR图像。此外,SAR图像中的目标一般在图像中只占很小的比例,使得目标检测精度严重下降。例如,使用Zhang等人在文章“Regions of Interest Detection in Panchromatic Remote Sensing ImagesBased on Multiscale Feature Fusion”中提出的基于多尺度特征融合的MFF模型进行SAR图像目标检测,将无法获取精确结果。
人眼在面对外界海量的信息时处理能力是有限的,视觉注意机制是保证信息处理过程中高效率的关键,它使得人眼只对其中显著的信息做出响应而舍弃其他不重要的信息。视觉显著性就是受此原理发展起来的,它的目标是在处理图像时使计算机可以具备类似人眼视觉的主动性和选择性。由于显著性目标检测具有不依赖先验知识、识别准确率高以及能快速检测出图像中的显著目标等优势,它已经被广泛应用于自然景物图像目标检测和光学遥感影像分析领域。利用显著性检测方法获取图像中的显著目标后,其结果可用于图像重建、图像融合和变化检测等领域,并且可以大大降低计算量。例如Imamoglu等人在文章“Asaliency detection model using low-level features based on wavelettransform”中提出的基于小波变换的WT模型,能够较为准确高效的计算自然景物图像的显著图并提取显著目标。因此,将视觉显著性引入SAR图像处理,能够为高效、精确地从SAR图像中提取出目标提供新的思路。
但是,将视觉显著性引入SAR图像目标检测还有一些问题需要注意。一方面SAR图像中缺少由可见光提供的真实色彩信息,因此需要针对SAR图像寻找新的特征,另一方面一组具有相似特征的SAR图像间存在共性关系,结合多幅图像中的共性显著特征进行目标检测,不仅保留了视觉显著性分析的优点,而且能够有效降低由真实色彩缺失所带来的SAR图像目标检测精度降低问题,还能够进一步抑制背景干扰,避免误判。联合显著性分析就是指在一组具有相似特征的图像中挖掘共有的显著目标,因此本发明使用联合显著性分析对SAR图像进行目标检测以提高检测精度。
发明内容
本发明的目的在于提供一种基于联合显著性分析的SAR图像油库区检测方法,该方法用于对一组SAR图像中的油库区进行精确检测。传统的SAR图像目标检测主要借鉴光学遥感影像的目标检测方法,首先需要对SAR图像中的斑点噪声进行建模,但是对SAR图像斑点噪声的建模是一个很复杂的问题,需要考虑多方面因素构建模型。此外SAR图像的成像机理不同于光学遥感影像,针对SAR图像寻找合适的特征也是一个重要问题。因此本发明方法主要关注三个方面:
1)降低斑点噪声对检测结果的不利影响;
2)针对SAR图像缺少光谱信息问题,寻找合适的新的特征;
3)挖掘多幅图像间的共有信息弥补SAR图像颜色信息的缺失,实现更为准确的SAR图像的油库区检测。
本方法针对一组具有相似地物特征的SAR图像进行处理,首先,利用方向平滑滤波器对一组SAR图像进行降噪处理,其次,提取一组SAR图像的亮度特征、纹理特征和曲线特征,进行共性显著特征分析,得到一组SAR图像的共性显著特征图,再次,基于共生直方图,计算得到一组SAR图像中每一幅图像的单图显著特征图,然后,融合共性显著特征图和单图显著特征图,得到一组SAR图像的最终显著图,最后通过最大类间方差法对最终显著图进行阈值分割提取油库区,具体包括以下步骤:
步骤一:进行SAR图像降噪处理,即构建一个增强的方向平滑滤波器,通过该滤波器的核函数在滑动窗口中执行低通滤波,使一组SAR图像中每幅图像的斑点噪声均得到有效降低;
步骤二:在一组SAR图像中,提取每幅图像各像素的亮度特征、纹理特征和曲线特征,即对于一组SAR图像中的每一幅图像,提取图像中各像素的亮度信息作为该像素的亮度特征,对图像进行下采样,在下采样图像上使用基于邻域阈值的局部二值模式提取每个像素的纹理特征,对图像中的所有像素点进行张量投票,计算得到每个像素的曲线特征,将提取的各像素的亮度特征、纹理特征以及曲线特征作为本发明中聚类使用的共性显著特征;
步骤三:利用模糊C均值聚类算法完成聚类,即在一组SAR图像中,利用每幅图像各像素的亮度特征、纹理特征和曲线特征,构建像素的共性显著特征向量,通过模糊C均值聚类算法,将一组SAR图像中所有像素的共性显著特征向量进行聚类,得到k个簇;
步骤四:计算一组SAR图像中每幅图像的共性显著特征图,即通过测量各像素的共性显著特征向量在各个簇中的出现频率,来计算簇与簇之间的距离,然后将第i个簇中含有的像素数与总像素数相除,相除的结果定义为第i个簇的权重,其中,得到所有k个簇的权重后,令一个簇到其他簇的距离的加权和为该簇的对比度显著值,把簇的对比度显著值赋给每一个属于该簇的像素点,由此获得一组共性显著特征图;
步骤五:计算一组SAR图像中每幅图像的单图显著特征图,即在一组SAR图像的每一幅图像中,统计所有强度值为a的像素的8邻域内,强度值为b的像素的个数,其中a和b的取值范围为0至255,然后将统计得到的像素对的个数排列成一个方阵,该方阵即为图像的共生直方图,然后用共生直方图中不同数值的出现概率的负对数,来描述直方图中每个点的初始显著值,最后增强初始显著值差异,得到一组SAR图像的单图显著特征图;
步骤六:计算一组SAR图像中每幅图像的最终显著图,即对一组SAR图像的每一幅单图显著特征图进行归一化,将归一化后的单图显著特征图中每个像素的灰度值记为对应SAR图像相同位置像素的显著性权重,用该权重来增强对应SAR图像的共性显著特征图,从而获得一组SAR图像的最终显著图;
步骤七:利用最大类间方差法准确提取一组SAR图像中每幅图像的油库区,即通过最大类间方差法计算得到最终显著图的分割阈值,利用该阈值将一组SAR图像中每幅图像的最终显著图分割为一幅二值图像模板,该模板中“1”代表油库区,“0”代表非油库区,最后将二值图像模板与对应SAR图像相乘,从而得到该SAR图像的油库区提取结果。
附图说明
图1为本发明的流程图。
图2为本发明的SAR图像共性显著特征聚类结果。第一行是SAR图像,第二行是聚类结果。
图3为本发明的共性显著性分析与单图显著性分析的结果对比。第一行是共性显著特征图Sp_m;第二行是单图显著特征图Sp_s。
图4为示例SAR图像及其对应的地面实况(Ground-Truth)。第一行是示例SAR图像;第二行是地面实况图。
图5为示例图片采用本发明方法和其他方法生成的显著图的比较。(a)-(d)分别表示采用MFF,WT,CSFA和本发明方法的结果,其中CSFA是发明内容中步骤二至步骤四所提出的共性显著特征分析方法(common salient feature analysis,CSFA)。
图6为示例图片采用本发明方法和其他方法所检测到的油库区比较。(a)-(d)分别表示采用MFF,WT,CSFA和本发明方法的结果。
具体实施方式
下面结合附图对本发明做进一步详细说明。本发明的总体框架如图1所示,现介绍每一步的实现细节。
步骤一:进行SAR图像降噪处理;
本发明中使用的一组尺寸为M×N的SAR图像,SAR图像中的斑点噪声通常被建模为纯乘性噪声过程,用下式表示
v(r,c)=l(r,c)s(r,c)
其中v表示雷达测量到的值,l表示图像的实际辐射值,s表示斑点噪声。对于单视SAR图像,s是瑞利分布或均值为1的负指数分布。对于具有独立外观的多视SAR图像,s是均值为1的伽马分布。M表示图像的行数,N表示图像的列数,r、c分别表示像素的横、纵坐标,r=1,2,...,M,c=1,2,...,N。
增强的方向平滑(enhanced directional smoothing,EDS)滤波器可以降低图像中的斑点噪声,它用核函数在滑动窗口中执行低通滤波。滤波窗口的典型尺寸范围为3×3到33×33,窗口的大小通常设为奇数。较大的滤波器窗口意味着可以利用较大的图像区域进行计算,但可能需要更多的计算时间,具体取决于滤波器算法的复杂程度。如果滤波器窗口的尺寸过大,重要细节将因过度平滑而丢失。另一方面,如果滤波器窗口的尺寸太小,则不利于降低斑点噪声。根据相关文献的实验结果分析,实验中采用3×3或7×7的滤波窗通常会产生最佳结果。本发明中使用的是一组尺寸为M×N的SAR图像,利用尺寸为3×3的滤波器降噪,将降噪后的SAR图像组记为该组SAR图像的总数为Q,Ix表示一组SAR图像的第x幅,x=1、2…Q。
步骤二:提取一组SAR图像中每幅图像所有像素的亮度特征、纹理特征和曲线特征;
在本方法中,将滤波后的SAR图像的亮度特征、纹理特征以及油罐特有的曲线特征作为本发明中聚类使用的共性显著特征。
亮度特征:直接使用Ix中像素pn的亮度信息作为该像素的亮度特征Gn。pn表示图像Ix中的第n个像素,n=1,2,…,M×N。
纹理特征:纹理特征是从Ix中提取的,主要包括两步。在本发明中,首先对Ix以1/8的比例进行下采样,然后在下采样后的图像上使用基于邻域阈值的局部二值模式来获取纹理特征图,计算每个像素的纹理特征Texturen的方法如下:
Texturen=∑ss(|Gnw-Gn|-δ)
其中Gn和Gnw分别表示当前像素及其邻域像素的灰度值。本发明中,pnw表示当前像素8-领域内的像素,δ是一个阈值,将其设为40。最后将纹理特征图恢复至初始大小就可获得每个像素的纹理特征。
曲线特征:张量投票是一种可以从噪声图像中提取显著特征结构的算法,因此本发明利用张量投票来获取SAR图像中每个像素的曲线特征Curven。二维平面中张量投票的投票域有两种:棒形张量和球形张量。
首先,对图像Ix中的所有像素点进行球形张量编码,即对其编码为:
根据矩阵谱分解原理,张量Tensor可分解如下的特征值λ和对应的特征向量e的线性组合。
Tensor=(λ1-λ2)e1e1 T+λ2(e1e1 T+e2e2 T)
其中,λ1和λ2为张量的非负特征值;e1和e2为两个特征值对应的特征向量;e1e1 T和e1e1 T+e2e2 T分別表示棒形张量和球形张量;对应的(λ1-λ2)和λ2分别表示线特征的显著性和点特征的显著性。因此,利用下式即可获得每个像素的曲线特征Curven。
Curven=λ1-λ2
步骤三:利用模糊C均值聚类算法完成聚类;
经过步骤二后,我们得到了Gn,Texturen和Curven三种SAR图像的共性显著特征。考虑到每个特征对聚类结果的贡献值,本发明对每个特征赋予不同的权重,然后将其组合为一个显著特征向量Featuren。
Featuren=(0.8·Gn,0.8·Texturen,Curven)
模糊C均值聚类(Fuzzy C-means algorithm,FCM)是一种柔性划分方法,主要是度量每个数据点属于某个簇的隶属度。FCM这种对样本类别划分的不确定性是一种更为客观的描述,因此本发明采用FCM完成多幅SAR图像的共性显著特征聚类。FCM使用u函数来确定每个数据属于某个簇的程度。
其中k表示簇的个数。ui(pn)表示像素pn属于第i个簇Ci的隶属度。FCM的目标函数J如下:
公式中的m是加权指数,默认m=2。ci是簇Ci的中心。FCM的算法步骤如下:首先,初始化各个簇中心ci使其满足的约束条件;然后,更新簇中心ci并计算目标函数J。计算J和前一步获得的J之间的差值,如果差值小于迭代条件则迭代终止,否则继续更新聚类结果。当k=3时,聚类结果如图2所示。从图中可以看出油库区被准确地分到同一个簇(黑色区域),海洋被分为一类(灰色区域),道路以及附近的建筑被分为一类(白色区域)。
步骤四:计算一组SAR图像中每幅图像的共性显著特征图;
全局对比度的计算基于一个通用规则:具有独特特征的稀有群集在直观上是引人注目的并且是更加显著的。为了达到这个目的,本发明首先构造了一个特征空间,空间的三个维度分别对应三个共性显著特征:Gn,Texturen和Curven。对每个维度分别进行均匀量化,量化后的级数分别表示为BG,BT和BC,其中,本发明设置BG=BT=BC=8。
对于每个簇,通过测量各特征向量在该簇中的出现频率来计算它在特征空间中的特征直方图。基于特征直方图,簇Ci的对比度显著值S(Ci)可以被估计为在特征空间中到其他簇的距离的加权和,具体计算如下。
ω(Ci)表示簇Ci中的像素数量与所有图像的总像素数量之比。D(ci,cj)表示簇Ci与簇Cj的距离,Ci表示第i个簇,BN代表直方图中特征量化的级数,BN=BG×BT×BC。另外,fi,t表示第t个特征向量出现在簇Ci中的频率。由此可获得共性显著特征图Sp_m。
步骤五:计算一组SAR图像中每幅图像的单图显著特征图;
共生直方图是表示像素之间的共生关系的一种统计方法;它可以揭示图像灰度值的空间关系的本质。共生直方图中的元素是图像中像素的强度信息。如果图像是由灰度值相似的像素块组成,则共生直方图的主对角线上元素将是相对较大的值;反之如果图像像素灰度值在局部有波动,则共生直方图的非主对角线元素具有相对较大的值。
对于一幅M×N的图像,共生直方图COH计算如下:
COH=[coh(a,b)]
COH是一个256×256的对称矩阵。coh(a,b)就表示强度值为a的像素的8-邻域内,强度值为b的像素的个数,强度值的取值范围为0至255。
COH就是由所有的coh(a,b)构建而成,然后使用下面的公式将COH归一化到0-1。
pCOH(a,b)表示的是共生直方图中不同数值的出现概率。理论上,对于研究者想要的显著地感兴趣区域,这些区域应该具有一个较小的pCOH(a,b),相反,对于非显著的区域,pCOH(a,b)应该是较大的数值。
为了使pCOH(a,b)与强度对(a,b)的显著性的关系是一个正相关的映射,受Boltzmann熵定理的启发,用负对数关系来描述每个亮度对的显著性关系,如下公式所示:
L(a,b)=-ln(pCOH(a,b))
其中,L(a,b)可以被看作是图像的初始显著值。因此显著区域应该具有一个较大的L(a,b),与此同时,对于非显著区域而言L(a,b)应该是一个较小的数值。但是,此时显著区和非显著区的初始显著值的差别不是很大,因此考虑增强这两者之间的差异。本发明利用K-means算法来增强两者之间的差异:首先,利用k-means将两者区分开,然后分别对两类做相应计算来增加两者之间的差异。增强计算如下:
其中,es(a,b)表示的是最后增强后的显著值。y(a,b)代表k-means算法的分类结果,值为1说明在该像素在较大数值的簇里,值为0则说明该像素位于较小数值的类别中。U表示共生直方图的平均显著值,经过此操作,显著区域与非显著区域之间的差异就会被进一步加大。
单图显著特征图中每个像素的显著值都取决于它的邻域,对于像素pn,它的显著值smn计算如下:
smn=∑es(Gn,Gnw)
由此得到的显著图中包含一些道路和其他冗余信息。因此,本发明用高斯滤波器来处理显著图,在一定程度上削弱了这些非显著区域的强度,由此可获得单图显著特征图Sp_s。
步骤六:计算一组SAR图像中每幅图像的最终显著图;
从图3可以看出:第一,单图显著性分析方法在显著背景干扰区域上产生了严重的误判,而共性显著性分析则准确地抑制其显著性。第二,相比于Sp_m,Sp_s中背景干扰更加明显。第三,虽然Sp_m能有效地抑制背景干扰,但是显著性较高的油库区中的部分油罐并不完整。
综上,这两种方法的检测结果都有各自的优点,所以本发明考虑融合这两种显著特征图来去除所有多余的信息的干扰。为了避免油库区中的部分油罐不完整,本发明使用单图显著特征图Sp_s来优化共性显著特征图Sp_m,以计算最终显著图MSP。
MSP=γ×Sp_m
γ=Sp_s/255
第一个公式的作用是融合单图显著特征图和共性显著特征图获得最终显著图,这里是一种点乘运算。首先对Sp_s进行归一化获得单幅图像中内部显著性的权重指标γ,然后用γ来增强Sp_m获得最终显著图MSP,如图5所示。从图中可以看出,MSP确实明显地突显出油库区,有效地抑制了非油库区的显著性。
步骤七:利用最大类间方差法准确提取SAR图像的油库区;
通过最大类间方差法得到最终显著图的分割阈值,利用该阈值将最终显著图分割为一幅二值图像模板,用“1”代表油库区,用“0”代表非油库区。最后将二值图像模板与SAR图像进行点乘得到最终的油库区提取结果。
本发明的效果可通过以下实验结果与分析进一步说明:
1.实验数据
实验中用的SAR数据是2014年8月29日在东京地区由ALOS2卫星拍摄获取的。该SAR图像是由卫星的L波段采用HH极化方式获得,分辨率为3m。我们从实验数据中随机选取7幅图像进行结果展示,图4展示了示例SAR图像及其对应的地面实况(Ground-Truth)。
2.对比实验
为了评价本发明方法的性能,我们设计了如下的对比实验,选取了现有的具有代表性的视觉注意方法选取了多特征融合方法(MFF),小波变换方法(WT),以及步骤二至步骤四所提出的共性显著性分析方法(CSFA)与本发明方法进行性能对比。从主观上分别对比了不同方法生成的显著图和感兴趣区域图,如图5和图6所示。
从结果可以看出:除了CSFA和本发明方法,其他的方法在显著背景干扰区域上都产生了严重的误判。MFF,WT的显著性检测结果中都包含有大量的背景信息,而CSFA和本发明则对背景信息进行了有效地抑制。这是由于后者采用对所有显著特征进行聚类然后计算显著值,可以有效地将油库区和背景区域进行区分。而前四种方法主要针对单图中的显著特征进行计算,常用的亮度以及频域分析无法准确地区分与油库区表现相似的背景区域。
Claims (1)
1.一种基于联合显著性分析的SAR图像油库区检测方法,本方法针对一组具有相似地物特征的SAR图像进行处理,首先,利用方向平滑滤波器对一组SAR图像进行降噪处理,其次,提取一组SAR图像的亮度特征、纹理特征和曲线特征,进行共性显著特征分析,得到一组SAR图像的共性显著特征图,再次,基于共生直方图,计算得到一组SAR图像中每一幅图像的单图显著特征图,然后,融合共性显著特征图和单图显著特征图,得到一组SAR图像的最终显著图,最后通过最大类间方差法对最终显著图进行阈值分割提取油库区,其特征在于,包括以下步骤:
步骤一:进行SAR图像降噪处理,即构建一个增强的方向平滑滤波器,通过该滤波器的核函数在滑动窗口中执行低通滤波,使一组SAR图像中每幅图像的斑点噪声均得到有效降低;
步骤二:在一组SAR图像中,提取每幅图像各像素的亮度特征、纹理特征和曲线特征,即对于一组SAR图像中的每一幅图像,提取图像中各像素的亮度信息作为该像素的亮度特征,对图像进行下采样,在下采样图像上使用基于邻域阈值的局部二值模式提取每个像素的纹理特征,对图像中的所有像素点进行张量投票,计算得到每个像素的曲线特征,将提取的各像素的亮度特征、纹理特征以及曲线特征作为聚类使用的共性显著特征;
步骤三:利用模糊C均值聚类算法完成聚类,即在一组SAR图像中,利用每幅图像各像素的亮度特征、纹理特征和曲线特征,构建像素的共性显著特征向量,通过模糊C均值聚类算法,将一组SAR图像中所有像素的共性显著特征向量进行聚类,得到k个簇;
步骤四:计算一组SAR图像中每幅图像的共性显著特征图,即通过测量各像素的共性显著特征向量在各个簇中的出现频率,来计算簇与簇之间的距离,然后将第i个簇中含有的像素数与总像素数相除,相除的结果定义为第i个簇的权重,其中i=1、2…k,得到所有k个簇的权重后,令一个簇到其他簇的距离的加权和为该簇的对比度显著值,把簇的对比度显著值赋给每一个属于该簇的像素点,由此获得一组共性显著特征图;
步骤五:计算一组SAR图像中每幅图像的单图显著特征图,即在一组SAR图像的每一幅图像中,统计所有强度值为a的像素的8邻域内,强度值为b的像素的个数,其中a和b的取值范围为0至255,然后将统计得到的像素对(a,b)的个数排列成一个方阵,该方阵即为图像的共生直方图,然后用共生直方图中不同数值的出现概率的负对数,来描述直方图中每个点的初始显著值,最后增强初始显著值差异,得到一组SAR图像的单图显著特征图;
步骤六:计算一组SAR图像中每幅图像的最终显著图,即对一组SAR图像的每一幅单图显著特征图进行归一化,将归一化后的单图显著特征图中每个像素的灰度值记为对应SAR图像相同位置像素的显著性权重,用该权重来增强对应SAR图像的共性显著特征图,从而获得一组SAR图像的最终显著图;
步骤七:利用最大类间方差法准确提取一组SAR图像中每幅图像的油库区,即通过最大类间方差法计算得到最终显著图的分割阈值,利用该阈值将一组SAR图像中每幅图像的最终显著图分割为一幅二值图像模板,该模板中“1”代表油库区,“0”代表非油库区,最后将二值图像模板与对应SAR图像点乘,从而得到该SAR图像的油库区提取结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810530745.7A CN108805057B (zh) | 2018-05-29 | 2018-05-29 | 一种基于联合显著性分析的sar图像油库区检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810530745.7A CN108805057B (zh) | 2018-05-29 | 2018-05-29 | 一种基于联合显著性分析的sar图像油库区检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108805057A CN108805057A (zh) | 2018-11-13 |
CN108805057B true CN108805057B (zh) | 2020-11-17 |
Family
ID=64090777
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810530745.7A Active CN108805057B (zh) | 2018-05-29 | 2018-05-29 | 一种基于联合显著性分析的sar图像油库区检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108805057B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110070549A (zh) * | 2019-04-25 | 2019-07-30 | 中国石油大学(华东) | 一种基于最优尺度邻域信息的海上溢油sar图像软分割方法 |
CN110956083A (zh) * | 2019-10-21 | 2020-04-03 | 山东科技大学 | 一种基于高分四号光学卫星的渤海海冰漂移遥感探测方法 |
CN113592033B (zh) * | 2021-08-20 | 2023-09-12 | 中科星睿科技(北京)有限公司 | 油罐图像识别模型训练方法、油罐图像识别方法和装置 |
CN113920283B (zh) * | 2021-12-13 | 2022-03-08 | 中国海洋大学 | 基于聚类分析和特征过滤的红外图像尾迹检测及提取方法 |
CN117314795B (zh) * | 2023-11-30 | 2024-02-27 | 成都玖锦科技有限公司 | 一种利用背景数据的sar图像增强方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104951799A (zh) * | 2015-06-12 | 2015-09-30 | 北京理工大学 | 一种sar遥感影像溢油检测识别方法 |
CN106096505A (zh) * | 2016-05-28 | 2016-11-09 | 重庆大学 | 基于多尺度特征协同表示的sar目标识别方法 |
CN106557740A (zh) * | 2016-10-19 | 2017-04-05 | 华中科技大学 | 一种遥感图像中油库目标的识别方法 |
CN107256409A (zh) * | 2017-05-22 | 2017-10-17 | 西安电子科技大学 | 基于sae和显著性检测的高分辨sar图像变化检测方法 |
CN107369131A (zh) * | 2017-07-04 | 2017-11-21 | 华中科技大学 | 图像的显著性检测方法、装置、存储介质和处理器 |
CN107832796A (zh) * | 2017-11-17 | 2018-03-23 | 西安电子科技大学 | 基于曲线波深度阶梯网络模型的sar图像分类方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9958521B2 (en) * | 2015-07-07 | 2018-05-01 | Q Bio, Inc. | Field-invariant quantitative magnetic-resonance signatures |
US20170016987A1 (en) * | 2015-07-17 | 2017-01-19 | Her Majesty The Queen In Right Of Canada, As Represented By The Minister Of National Defence | Processing synthetic aperture radar images for ship detection |
-
2018
- 2018-05-29 CN CN201810530745.7A patent/CN108805057B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104951799A (zh) * | 2015-06-12 | 2015-09-30 | 北京理工大学 | 一种sar遥感影像溢油检测识别方法 |
CN106096505A (zh) * | 2016-05-28 | 2016-11-09 | 重庆大学 | 基于多尺度特征协同表示的sar目标识别方法 |
CN106557740A (zh) * | 2016-10-19 | 2017-04-05 | 华中科技大学 | 一种遥感图像中油库目标的识别方法 |
CN107256409A (zh) * | 2017-05-22 | 2017-10-17 | 西安电子科技大学 | 基于sae和显著性检测的高分辨sar图像变化检测方法 |
CN107369131A (zh) * | 2017-07-04 | 2017-11-21 | 华中科技大学 | 图像的显著性检测方法、装置、存储介质和处理器 |
CN107832796A (zh) * | 2017-11-17 | 2018-03-23 | 西安电子科技大学 | 基于曲线波深度阶梯网络模型的sar图像分类方法 |
Non-Patent Citations (1)
Title |
---|
Saliency detection and region of interest extraction based on multi-image common saliency analysis in satellite images;Libao Zhang.et.;《Neurocomputing》;20180329;第283卷;第150-165页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108805057A (zh) | 2018-11-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108805057B (zh) | 一种基于联合显著性分析的sar图像油库区检测方法 | |
Li et al. | A review of remote sensing image classification techniques: The role of spatio-contextual information | |
Chitradevi et al. | An overview on image processing techniques | |
Touati et al. | An energy-based model encoding nonlocal pairwise pixel interactions for multisensor change detection | |
CN106170819B (zh) | 一种合成孔径雷达图像舰船目标快速检测方法 | |
CN109035188B (zh) | 一种基于目标特征驱动的智能图像融合方法 | |
Tang et al. | Grabcut in one cut | |
CN109871902B (zh) | 一种基于超分辨率对抗生成级联网络的sar小样本识别方法 | |
Yousif et al. | Improving SAR-based urban change detection by combining MAP-MRF classifier and nonlocal means similarity weights | |
CN110309781B (zh) | 基于多尺度光谱纹理自适应融合的房屋损毁遥感识别方法 | |
CN106846322B (zh) | 基于曲线波滤波器和卷积结构学习的sar图像分割方法 | |
CN110458192B (zh) | 基于视觉显著性的高光谱遥感图像分类方法及*** | |
CN104103082A (zh) | 一种基于区域描述和先验知识的图像显著性检测方法 | |
CN111027497B (zh) | 基于高分辨率光学遥感影像的弱小目标快速检测方法 | |
Rusyn et al. | Segmentation of atmospheric clouds images obtained by remote sensing | |
CN115690086A (zh) | 一种基于对象的高分辨率遥感影像变化检测方法及*** | |
Sun et al. | Probabilistic neural network based seabed sediment recognition method for side-scan sonar imagery | |
CN108805186B (zh) | 一种基于多维显著特征聚类的sar图像圆形油库检测方法 | |
CN108875798B (zh) | 一种基于空间金字塔池化的超像素级特征提取方法 | |
Scharfenberger et al. | Image saliency detection via multi-scale statistical non-redundancy modeling | |
Cheng et al. | Tensor locality preserving projections based urban building areas extraction from high-resolution SAR images | |
Toure et al. | Coastline detection using fusion of over segmentation and distance regularization level set evolution | |
CN113963270A (zh) | 一种高分遥感影像建筑物检测方法 | |
Ko et al. | Adaptive growing and merging algorithm for image segmentation | |
Sun et al. | A New Method for Edge Detection Based on the Criterion of Separability. |
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 |