CN110991257B - 基于特征融合与svm的极化sar溢油检测方法 - Google Patents
基于特征融合与svm的极化sar溢油检测方法 Download PDFInfo
- Publication number
- CN110991257B CN110991257B CN201911094691.5A CN201911094691A CN110991257B CN 110991257 B CN110991257 B CN 110991257B CN 201911094691 A CN201911094691 A CN 201911094691A CN 110991257 B CN110991257 B CN 110991257B
- Authority
- CN
- China
- Prior art keywords
- cnn network
- features
- sar
- pol
- oil spill
- 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
- 239000003305 oil spill Substances 0.000 title claims abstract description 47
- 238000001514 detection method Methods 0.000 title claims abstract description 37
- 230000004927 fusion Effects 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 claims abstract description 51
- 230000008569 process Effects 0.000 claims abstract description 27
- 230000009467 reduction Effects 0.000 claims abstract description 20
- 238000000605 extraction Methods 0.000 claims abstract description 15
- 230000010287 polarization Effects 0.000 claims abstract description 11
- 238000012549 training Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 21
- 238000000513 principal component analysis Methods 0.000 claims description 20
- 238000011176 pooling Methods 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 12
- 230000001427 coherent effect Effects 0.000 claims description 9
- 238000013178 mathematical model Methods 0.000 claims description 6
- 210000002569 neuron Anatomy 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000009977 dual effect Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000007781 pre-processing Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 239000000284 extract Substances 0.000 abstract description 9
- 239000012075 bio-oil Substances 0.000 abstract description 3
- 238000013527 convolutional neural network Methods 0.000 description 53
- 238000012706 support-vector machine Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 7
- 230000006872 improvement Effects 0.000 description 6
- 238000011160 research Methods 0.000 description 5
- 238000013135 deep learning Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 239000003921 oil Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 3
- 239000013535 sea water Substances 0.000 description 3
- 230000004913 activation Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000007774 longterm Effects 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
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/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
- G06F18/2135—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods based on approximation criteria, e.g. principal component analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
- G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
- G06F18/2411—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on the proximity to a decision surface, e.g. support vector machines
-
- 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
- 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
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computing Systems (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Remote Sensing (AREA)
- Astronomy & Astrophysics (AREA)
- Multimedia (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了基于特征融合与SVM的极化SAR溢油检测方法,包括以下步骤:构建Pol‑SAR数据块,然后构建溢油样本数据集;搭建CNN网络,并从Pol‑SAR数据中抽取训练样本对CNN网络模型进行训练;利用训练好的CNN网络模型,对Pol‑SAR数据进行深度特征提取,提取高层次的特征;将不同的高层次的特征进行融合,利用带有RBF核函数的SVM分类器进行分类;本发明方法避免了传统方法需要手动提取并选择特征的复杂过程,能够有效的减少相干斑噪声的影响,减少溢油检测的虚警率,提取出CNN网络的两层高层特征,通过PCA降维对两层特征进行特征融合,可以有效地区分溢油与生物油膜,提高溢油检测精度。
Description
技术领域
本发明涉及溢油检测技术领域,尤其涉及基于特征融合与SVM的极化SAR溢油检测方法。
背景技术
海上石油泄漏往往造成大规模的海洋污染,对海洋生态系统和环境、渔业、野生动物和其他社会利益造成严重危害,开展及时有效的海洋溢油监测具有重要意义。其中,合成孔径雷达(Synthetic Aperture Radar,SAR)由于不受阳光、云层的影响,且具有全天时、全天候的工作特点及大范围观测以及高分辨率的优势,已成为海洋溢油检测的主要手段;
近年来,Po1-SAR海面溢油检测已成为研究热点,许多研究证实了Pol-SAR在海面溢油检测方面的有效性。然而,在Pol-SAR图像目标检测和分类中地物特征的自动提取和选择特征是一项长期的挑战,随着深度学习在光学图像识别中的成功应用,其在SAR目标识别方面的研究工作也已经开始,基于深度学习的识别技术具有多级特征表达和非线性数据拟合的能力,它可以自动挖掘更具有鉴别性和代表性的特征,从而保证了分类性能,其中CNN能够挖掘数据中的空间相关性;
目前大多数研究都是基于目标电磁波散射特性从Pol-SAR数据中人工的提取特征并进行检测,溢油检测的准确性在很大程度上取决于提取的特征的优劣,人工提取特征复杂且耗时,需要长期的研究和丰富的Pol-SAR图像处理经验,且提取的特征可能存在信息涵盖不全面或存在冗余等现象,且不同的特征在不同的海况下性能会有较大的差异,从而导致溢油检测性能不高,因此,本发明提出基于特征融合与SVM的极化SAR溢油检测方法,以解决现有技术中的不足之处。
发明内容
针对上述问题,本发明提出基于特征融合与SVM的极化SAR溢油检测方法,避免了传统方法需要手动提取并选择特征的复杂过程,通过CNN网络自动的从原始Pol-SAR数据中提取鲁棒的有辨别性的特征,考虑了像素周边邻域信息,能够有效的减少相干斑噪声的影响,减少溢油检测的虚警率,提取出CNN网络的两层高层特征,通过PCA降维,对两层特征进行特征融合,涵盖的信息更加丰富,可以有效地区分溢油与生物油膜,提高溢油检测精度。
本发明提出基于特征融合与SVM的极化SAR溢油检测方法,包括以下步骤:
步骤一:对Pol-SAR数据进行预处理,首先构建Pol-SAR数据块,然后构建溢油样本数据集,保证Pol-SAR数据适应CNN网络的输入;
步骤二:搭建CNN网络,并从Pol-SAR数据中抽取训练样本对CNN网络模型进行训练;
步骤三:利用步骤二中训练好的CNN网络模型,先保存整个CNN网络模型,将CNN网络模型中的每一个卷积池化层当作是一个特征提取层,然后对Pol-SAR数据进行深度特征提取,获得深度特征,并引入主成分分析进行特征降维;
所述步骤三中的具体过程为:
(1)首先确定一个滑动窗口,大小为15×15,从步骤一中构建的原始的Pol-SAR数据块的左上角第一个像素(x,y),x∈m,y∈n开始,边界部分采用镜像模式填充,取出15×15×9的小数据块,输入到步骤二中训练好的CNN网络模型中;
(2)提取CNN网络第3个卷积层后并拉伸为1维的输出值,该输出是一个1×80的特征向量,记为fea_1,再初始化一个三维矩阵F1(m,n,80),将fea_1保存到相应的位置(x,y),x∈m,y∈n中;
(3)提取全连接层的输出值,记为fea_2,其大小为1×240,然后构建一个用于保存该特征的三维矩阵F2(m,n,240),将fea_2填入到相应的位置(x,y)中;
(4)对于x+1,y+1,重复上述步骤(1)-(3),直到遍历完整个图像;
(5)对于F1和F2,引入PCA分别对两个数据进行降维,取降维后的前3个主成分作为降维特征,记为PCA_fea_1,PCA_fea_2,尺寸大小均为m×n×3;
步骤四:将步骤(5)中得到的降维特征进行融合,将降维特征进行堆叠,得到融合后的特征,记为fuse_fea,尺寸大小为m×n×6,最后利用带有RBF核函数的SVM分类器进行分类。
进一步改进在于:所述步骤一中构建Pol-SAR数据块的具体过程为:Pol-SAR数据中的每个像素呈现为2×2散射矩阵,包含用于描述相干性或纯散射体的丰富的潜在信息,在多视和互易假设下,将单站Pol-SAR由对称的3×3复数相干矩阵T3进行表示,复数相干矩阵T3形式如公式(1)所示:
其中,T11,T22,T33是实数,其他矩阵元素均是复数;
将矩阵复数相干矩阵T3作为CNN的输入,让对于每个像素,极化数据可以定义为矢量t,其表达如公式(2)所示:
t=[T11,T22,T33,Re(T12),Im(T12),Re(T13),Im(T13),Re(T23),Im(T23)] (2)
其中,Re表示实部,Im表示虚部;
Pol-SAR数据为一个包含9个通道的图像数据块D(m,n,bands),其中:m表示图像的长,n表示宽度,bands表示图像的通道数,该图像数据块用于生成样本,作为CNN的输入。
进一步改进在于:所述步骤一中构建溢油样本数据集的具体过程为:对步骤一中产生的9个通道的图像数据块,从每一类地物中随机选取4000个样本点,取其周围15*15的小块作为一个样本,记为patch,其中心像素点对应的标签作为该样本的标签,记为label(i),i∈N,N为地物类别数。
进一步改进在于:所述步骤二中搭建CNN网络的具体过程为:搭建一个适合于从PolSAR数据中自动的提取深度特征的CNN网络模型,该CNN网络包含3个卷积层、两个最大池化层、一个扁平化层、一个全连接层以及连接到输出的Softmax分类器。
进一步改进在于:所述CNN网络包含的3个卷积层中前2个卷积层分别包含30和60个卷积核,大小为3×3,用于从输入数据中自动提取特征;前2个卷积层后面分别接着一个最大池化层,池大小为2×2,步长为2个像素;第3个卷积层采用了1×1卷积,卷积核的个数为20,用于删除冗余特征并提高特征的代表性,降低计算复杂度,然后对所有的卷积层均进行L2正则化,再将得到的特征图拉伸为一个1维向量,用于输入全连接网络;所述全连接层的神经元个数为240,最后一层神经元个数为待分类的个数,并在所述全连接层后增加了Dropout操作;所述Softmax分类器用于计算每个类的概率。
进一步改进在于:所述步骤二中对CNN网络模型进行训练的具体过程为:将取出的溢油训练样本输入搭建好的搭建CNN网络中,经过卷积层、最大池化层、扁平化层、一个全连接层前向传播,然后得到CNN网络的预测概率值,最后基于预测值和真实的样本标签值计算损失函数,计算公式如公式(4)所示:
然后根据BackPropagation算法对CNN网络参数进行迭代更新,在反向传播中,采用随机梯度下降算法,并在梯度下降的过程中使用Adadelta优化算法进行自适应调整学习率,最后得到一个预估结果。
进一步改进在于:所述步骤四中具体过程为:基于CNN网络模型,利用卷积池化层从Pol-SAR图像中提取深度特征,然后使用全连接层和Softmax层对特征进行分类,并在分类时引入RBF-SVM算法,对提取的深度融合特征进行分类,设T={(x1,y1),(x2,y2)...,(xn,yn)}为待分类的样本,RBF-SVM算法的数学模型如公式(5)所示:
其中,w表示法向量,决定了超平面的方向;b表示偏移量,决定了超平面与原点之间的距离;yi是正确的分类标签;C代表惩罚系数;ξi是松弛系数;
将公式(4)的数学模型变形为拉格朗日对偶问题,其表达式如公式(6)所示:
其中,αi代表拉格朗日权重;k(xi,xj)表示RBF-SVM的核函数。
本发明的有益效果为:避免了传统方法需要手动提取并选择特征的复杂过程,将卷积神经网络引入到溢油检测中,通过CNN网络自动的从原始Pol-SAR数据中提取鲁棒的有辨别性的特征,学习到的深度特征,考虑了像素周边邻域信息,能够有效的减少相干斑噪声的影响,减少溢油检测的虚警率,此外,对于一般的深度学习方法未能融合多层深度特征,本发明方法提取出CNN网络的两层高层特征,通过PCA降维,对两层特征进行特征融合,涵盖的信息更加丰富,最后将融合得到的特征输入到SVM分类器得到最终检测结果,本发明方法可以有效地区分溢油与生物油膜,提高溢油检测精度。
附图说明
图1为本发明方法流程示意图;
图2为本发明中PolSAR数据的9-通道形式结构示意图;
图3为本发明中CNN网络结构图;
图4为本发明中数据1的海洋溢油检测分类结果示意图;其中,(a)T3-SVM;(b)PF-SVM;(c)T3-CNN;(d)本发明方法;(e)解译图;(Oilspill:溢油、Seawater:海水);
图5本发明中数据2的海洋溢油检测分类结果示意图。
具体实施方式
为了加深对本发明的理解,下面将结合实施例对本发明做进一步详述,本实施例仅用于解释本发明,并不构成对本发明保护范围的限定。
根据图1、2、3、4、5所示,本实施例提出基于特征融合与SVM的极化SAR溢油检测方法,包括以下步骤:
步骤一:对Pol-SAR数据进行预处理,首先构建Pol-SAR数据块,然后构建溢油样本数据集,保证Pol-SAR数据适应CNN网络的输入;
构建Pol-SAR数据块的具体过程为:Pol-SAR数据中的每个像素呈现为2×2散射矩阵,包含用于描述相干性或纯散射体的丰富的潜在信息,在多视和互易假设下,将单站Pol-SAR由对称的3×3复数相干矩阵T3进行表示,复数相干矩阵T3形式如公式(1)所示:
其中,T11,T22,T33是实数,其他矩阵元素均是复数;
将矩阵复数相干矩阵T3作为CNN的输入,让对于每个像素,极化数据可以定义为矢量t,其表达如公式(2)所示:
t=[T11,T22,T33,Re(T12),Im(T12),Re(T13),Im(T13),Re(T23),Im(T23)] (2)
其中,Re表示实部,Im表示虚部;
Pol-SAR数据为一个包含9个通道的图像数据块D(m,n,bands),其中:m表示图像的长,n表示宽度,bands表示图像的通道数,该图像数据块用于生成样本,作为CNN的输入;
构建溢油样本数据集的具体过程为:对步骤一中产生的9个通道的图像数据块,从每一类地物中随机选取4000个样本点,取其周围15*15的小块作为一个样本,记为patch,其中心像素点对应的标签作为该样本的标签,记为label(i),i∈N,N为地物类别数;
步骤二:搭建CNN网络,并从Pol-SAR数据中抽取训练样本对CNN网络模型进行训练;
CNN的基本结构由输入层、卷积层、池化层和全连接层堆叠组成,所有层都是串联连接的,每一层的输入数据是前一层的输出数据,输入层用于接收图像数据,由于层级连接结构,CNN可以从低级特征中提取高级特征,
其中卷积层也称为特征提取层,通过对输入数据进行卷积操作得到特征图,假设第l层为卷积层,则第l层第j个特征图的计算公式如下所示:
其中,表示第l-1层第i个特征图,Mj表示输入的特征图,表示卷积操作后得到的输出特征图,表示第l-1层的第i个特征图与第l层第j个特征图之间的卷积核;表示偏差,f(·)代表非线性激活函数,如sigmoid函数和整流线性单位(ReLU);
搭建CNN网络的具体过程为:搭建一个适合于从PolSAR数据中自动的提取深度特征的CNN网络模型,该CNN网络包含3个卷积层、两个最大池化层、一个扁平化层、一个全连接层以及连接到输出的Softmax分类器;
所述CNN网络包含的3个卷积层中前2个卷积层分别包含30和60个卷积核,大小为3×3,用于从输入数据中自动提取特征;前2个卷积层后面分别接着一个最大池化层,池大小为2×2,步长为2个像素;第3个卷积层采用了1×1卷积,卷积核的个数为20,用于删除冗余特征并提高特征的代表性,降低计算复杂度,然后对所有的卷积层均进行L2正则化,再将得到的特征图拉伸为一个1维向量,用于输入全连接网络;所述全连接层的神经元个数为240,最后一层神经元个数为待分类的个数,本实施例中设置为,并在所述全连接层后增加了Dropout操作,Dropout参数设置为0.5;所述Softmax分类器用于计算每个类的概率,ReLU激活函数应用于所有的卷积层以及全连接层;
对CNN网络模型进行训练的具体过程为:将取出的溢油训练样本输入搭建好的搭建CNN网络中,经过卷积层、最大池化层、扁平化层、一个全连接层前向传播,然后得到CNN网络的预测概率值,最后基于预测值和真实的样本标签值计算损失函数,计算公式如公式(4)所示:
然后根据Back Propagation算法对CNN网络参数进行迭代更新,在反向传播中,采用随机梯度下降算法,并在梯度下降的过程中使用Adadelta优化算法进行自适应调整学习率,其中batchsize设置为64,最后得到一个预估结果
步骤三:利用步骤二中训练好的CNN网络模型,先保存整个CNN网络模型,将CNN网络模型中的每一个卷积池化层当作是一个特征提取层,然后对Pol-SAR数据进行深度特征提取,获得深度特征,并引入主成分分析进行特征降维;
(1)首先确定一个滑动窗口,大小为15×15,从步骤一中构建的原始的Pol-SAR数据块的左上角第一个像素(x,y),x∈m,y∈n开始,边界部分采用镜像模式填充,取出15×15×9的小数据块,输入到步骤二中训练好的CNN网络模型中;
(2)提取CNN网络第3个卷积层后并拉伸为1维的输出值,该输出是一个1×80的特征向量,记为fea_1,再初始化一个三维矩阵F1(m,n,80),将fea_1保存到相应的位置(x,y),x∈m,y∈n中;
(3)提取全连接层的输出值,记为fea_2,其大小为1×240,然后构建一个用于保存该特征的三维矩阵F2(m,n,240),将fea_2填入到相应的位置(x,y)中;
(4)对于x+1,y+1,重复上述步骤(1)-(3),直到遍历完整个图像;
(5)对于F1和F2,引入PCA分别对两个数据进行降维,取降维后的前3个主成分作为降维特征,记为PCA_fea_1,PCA_fea_2,尺寸大小均为m×n×3;
步骤四:将步骤(5)中得到的降维特征进行融合,将降维特征进行堆叠,得到融合后的特征,记为fuse_fea,尺寸大小为m×n×6,最后利用带有RBF核函数的SVM分类器进行分类;
基于CNN网络模型,利用卷积池化层从Pol-SAR图像中提取深度特征,然后使用全连接层和Softmax层对特征进行分类,并在分类时引入RBF-SVM算法,对提取的深度融合特征进行分类,设T={(x1,y1),(x2,y2)...,(xn,yn)}为待分类的样本,RBF-SVM算法的数学模型如公式(5)所示:
其中,w表示法向量,决定了超平面的方向;b表示偏移量,决定了超平面与原点之间的距离;yi是正确的分类标签;C代表惩罚系数;ξi是松弛系数;
将公式(4)的数学模型变形为拉格朗日对偶问题,其表达式如公式(6)所示:
其中,αi代表拉格朗日权重;k(xi,xj)表示RBF-SVM的核函数。
为了验证本发明方法的有效性,设置了对比实验;
对比实验1:基于T3原始数据的SVM分类方法(T3-SVM);
对比实验2:基于极化特征和SVM的分类方法(PF-SVM),选取了几种溢油检测中较为常用的极化特征组合作为SVM的输入,分别是Span,H,α,μ,DoP;
对比实验3:基于CNN的分类方法(T3-CNN);
以两景RADARSAT-2全极化SAR溢油数据作为实验数据1和数据2,图4和图5分别是本发明方法和其他算法的检测结果,表1和表2是精度评定结果,其中,OA(Overallaccuracy)代表总体分类精度,Kappa为Kappa系数;
表1数据1溢油检测分类精度
表2数据2溢油检测分类精度
由图4、图5以及表1和表2的结果可以比较直观的看出,本发明方法在目视效果上纯净度较高,在海面部分存在较少的杂斑,海水错分为溢油的部分较其他方法大大减少,相比于其他方法,本发明方法有效的融合了两层高层特征,并使用SVM分类器替换了CNN中较弱的softmax分类器,增强了解决非线性问题的能力;表1可以看到本发明方法OA和kappa系数均都高于其他方法,分别为99.19%和0.9587;
对于图4,主要用于验证本发明方法对于溢油与生物油膜的区分能力,本发明方法检测效果最好,OA和kappa系数最高,分别为97.56%和0.7795,与其他方法相比,不仅明显降低了海面的杂斑,而且对于生物油膜也能够很好的区分,此外还减少了生物油膜误分为溢油的现象,说明本发明方法能够有效提高溢油的分类精度。
本发明的有益效果为:避免了传统方法需要手动提取并选择特征的复杂过程,将卷积神经网络引入到溢油检测中,通过CNN网络自动的从原始Pol-SAR数据中提取鲁棒的有辨别性的特征,学习到的深度特征,考虑了像素周边邻域信息,能够有效的减少相干斑噪声的影响,减少溢油检测的虚警率,此外,对于一般的深度学习方法未能融合多层深度特征,本发明方法提取出CNN网络的两层高层特征,通过PCA降维,对两层特征进行特征融合,涵盖的信息更加丰富,最后将融合得到的特征输入到SVM分类器得到最终检测结果,本发明方法可以有效地区分溢油与生物油膜,提高溢油检测精度。
以上显示和描述了本发明的基本原理、主要特征和优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (7)
1.基于特征融合与SVM的极化SAR溢油检测方法,其特征在于,包括以下步骤:
步骤一:对Pol-SAR数据进行预处理,首先构建Pol-SAR数据块,然后构建溢油样本数据集,保证Pol-SAR数据适应CNN网络的输入;
步骤二:搭建CNN网络,并从Pol-SAR数据中抽取训练样本对CNN网络模型进行训练;
步骤三:利用步骤二中训练好的CNN网络模型,先保存整个CNN网络模型,将CNN网络模型中的每一个卷积池化层当作是一个特征提取层,然后对Pol-SAR数据进行深度特征提取,获得深度特征,并引入主成分分析进行特征降维;
所述步骤三中的具体过程为:
(1)首先确定一个滑动窗口,大小为15×15,从步骤一中构建的原始的Pol-SAR数据块的左上角第一个像素(x,y),x∈m,y∈n开始,边界部分采用镜像模式填充,取出15×15×9的小数据块,输入到步骤二中训练好的CNN网络模型中;
(2)提取CNN网络第3个卷积层后并拉伸为1维的输出值,该输出是一个1×80的特征向量,记为fea_1,再初始化一个三维矩阵F1(m,n,80),将fea_1保存到相应的位置(x,y),x∈m,y∈n中;
(3)提取全连接层的输出值,记为fea_2,其大小为1×240,然后构建一个用于保存该特征的三维矩阵F2(m,n,240),将fea_2填入到相应的位置(x,y)中;
(4)对于x+1,y+1,重复上述步骤(1)-(3),直到遍历完整个图像;
(5)对于F1和F2,引入PCA分别对两个数据进行降维,取降维后的前3个主成分作为降维特征,记为PCA_fea_1,PCA_fea_2,尺寸大小均为m×n×3;
步骤四:将步骤(5)中得到的降维特征进行融合,将降维特征进行堆叠,得到融合后的特征,记为fuse_fea,尺寸大小为m×n×6,最后利用带有RBF核函数的SVM分类器进行分类。
2.根据权利要求1所述的基于特征融合与SVM的极化SAR溢油检测方法,其特征在于:所述步骤一中构建Pol-SAR数据块的具体过程为:Pol-SAR数据中的每个像素呈现为2×2散射矩阵,包含用于描述相干性或纯散射体的丰富的潜在信息,在多视和互易假设下,将单站Pol-SAR由对称的3×3复数相干矩阵T3进行表示,复数相干矩阵T3形式如公式(1)所示:
其中,T11,T22,T33是实数,其他矩阵元素均是复数;
将矩阵复数相干矩阵T3作为CNN的输入,让对于每个像素,极化数据可以定义为矢量t,其表达如公式(2)所示:
t=[T11,T22,T33,Re(T12),Im(T12),Re(T13),Im(T13),Re(T23),Im(T23)] (2)
其中,Re表示实部,Im表示虚部;
Pol-SAR数据为一个包含9个通道的图像数据块D(m,n,bands),其中:m表示图像的长,n表示宽度,bands表示图像的通道数,该图像数据块用于生成样本,作为CNN的输入。
3.根据权利要求2所述的基于特征融合与SVM的极化SAR溢油检测方法,其特征在于:所述步骤一中构建溢油样本数据集的具体过程为:对步骤一中产生的9个通道的图像数据块,从每一类地物中随机选取4000个样本点,取其周围15*15的小块作为一个样本,记为patch,其中心像素点对应的标签作为该样本的标签,记为label(i),i∈N,N为地物类别数。
4.根据权利要求1所述的基于特征融合与SVM的极化SAR溢油检测方法,其特征在于:所述步骤二中搭建CNN网络的具体过程为:搭建一个适合于从PolSAR数据中自动的提取深度特征的CNN网络模型,该CNN网络包含3个卷积层、两个最大池化层、一个扁平化层、一个全连接层以及连接到输出的Softmax分类器。
5.根据权利要求4所述的基于特征融合与SVM的极化SAR溢油检测方法,其特征在于:所述CNN网络包含的3个卷积层中前2个卷积层分别包含30和60个卷积核,大小为3×3,用于从输入数据中自动提取特征;前2个卷积层后面分别接着一个最大池化层,池大小为2×2,步长为2个像素;第3个卷积层采用了1×1卷积,卷积核的个数为20,用于删除冗余特征并提高特征的代表性,降低计算复杂度,然后对所有的卷积层均进行L2正则化,再将得到的特征图拉伸为一个1维向量,用于输入全连接网络;所述全连接层的神经元个数为240,最后一层神经元个数为待分类的个数,并在所述全连接层后增加了Dropout操作;所述Softmax分类器用于计算每个类的概率。
6.根据权利要求4所述的基于特征融合与SVM的极化SAR溢油检测方法,其特征在于:所述步骤二中对CNN网络模型进行训练的具体过程为:将取出的溢油训练样本输入搭建好的搭建CNN网络中,经过卷积层、最大池化层、扁平化层、一个全连接层前向传播,然后得到CNN网络的预测概率值,最后基于预测值和真实的样本标签值计算损失函数,计算公式如公式(4)所示:
然后根据Back Propagation算法对CNN网络参数进行迭代更新,在反向传播中,采用随机梯度下降算法,并在梯度下降的过程中使用Adadelta优化算法进行自适应调整学习率,最后得到一个预估结果。
7.根据权利要求1所述的基于特征融合与SVM的极化SAR溢油检测方法,其特征在于:所述步骤四中具体过程为:基于CNN网络模型,利用卷积池化层从Pol-SAR图像中提取深度特征,然后使用全连接层和Softmax层对特征进行分类,并在分类时引入RBF-SVM算法,对提取的深度融合特征进行分类,设T={(x1,y1),(x2,y2)...,(xn,yn)}为待分类的样本,RBF-SVM算法的数学模型如公式(5)所示:
其中,w表示法向量,决定了超平面的方向;b表示偏移量,决定了超平面与原点之间的距离;yi是正确的分类标签;C代表惩罚系数;ξi是松弛系数;
将公式(4)的数学模型变形为拉格朗日对偶问题,其表达式如公式(6)所示:
其中,αi代表拉格朗日权重;k(xi,xj)表示RBF-SVM的核函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911094691.5A CN110991257B (zh) | 2019-11-11 | 2019-11-11 | 基于特征融合与svm的极化sar溢油检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911094691.5A CN110991257B (zh) | 2019-11-11 | 2019-11-11 | 基于特征融合与svm的极化sar溢油检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110991257A CN110991257A (zh) | 2020-04-10 |
CN110991257B true CN110991257B (zh) | 2023-04-18 |
Family
ID=70083732
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911094691.5A Active CN110991257B (zh) | 2019-11-11 | 2019-11-11 | 基于特征融合与svm的极化sar溢油检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110991257B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2022078754A (ja) * | 2020-11-13 | 2022-05-25 | 株式会社東芝 | 画像識別装置、レーダ装置、画像識別方法、およびプログラム |
JP7006875B1 (ja) * | 2021-06-18 | 2022-01-24 | 株式会社スペースシフト | 学習モデル、信号処理装置、飛翔体、及びプログラム |
CN115485742B (zh) | 2021-06-18 | 2023-06-27 | 空间转移株式会社 | 学习模型、信号处理装置、飞行物以及程序 |
CN113570589B (zh) * | 2021-08-03 | 2022-04-19 | 河海大学 | 一种基于特征融合的深度学习sar影像溢油区识别方法 |
CN114724030B (zh) * | 2022-04-06 | 2023-06-02 | 西安电子科技大学 | 基于对比学习的极化sar地物分类方法 |
CN117095308B (zh) * | 2023-10-19 | 2024-01-05 | 中国石油大学(华东) | 一种基于遗传算法的全极化sar影像海面溢油检测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109490886A (zh) * | 2018-11-28 | 2019-03-19 | 北京工业大学 | 一种极化合成孔径雷达遥感海面溢油区域精确提取方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563422B (zh) * | 2017-08-23 | 2019-08-27 | 西安电子科技大学 | 一种基于半监督卷积神经网络的极化sar分类方法 |
CN108932480B (zh) * | 2018-06-08 | 2022-03-15 | 电子科技大学 | 基于1d-cnn的分布式光纤传感信号特征学习与分类方法 |
-
2019
- 2019-11-11 CN CN201911094691.5A patent/CN110991257B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109490886A (zh) * | 2018-11-28 | 2019-03-19 | 北京工业大学 | 一种极化合成孔径雷达遥感海面溢油区域精确提取方法 |
Non-Patent Citations (1)
Title |
---|
万剑华 ; 韩仲志 ; 宋欣欣 ; 刘杰 ; .多模式融合下的海洋溢油高光谱成像油种识别方法.《发光学报》.2016,(第04期),473-480. * |
Also Published As
Publication number | Publication date |
---|---|
CN110991257A (zh) | 2020-04-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110991257B (zh) | 基于特征融合与svm的极化sar溢油检测方法 | |
CN110443143B (zh) | 多分支卷积神经网络融合的遥感图像场景分类方法 | |
CN110472627B (zh) | 一种端到端的sar图像识别方法、装置及存储介质 | |
CN106845478B (zh) | 一种字符置信度的二次车牌识别方法及装置 | |
CN106127204B (zh) | 一种全卷积神经网络的多方向水表读数区域检测算法 | |
CN112288008B (zh) | 一种基于深度学习的马赛克多光谱图像伪装目标检测方法 | |
Pally et al. | Application of image processing and convolutional neural networks for flood image classification and semantic segmentation | |
Liu et al. | Remote sensing image change detection based on information transmission and attention mechanism | |
CN111709313B (zh) | 基于局部和通道组合特征的行人重识别方法 | |
CN109359661B (zh) | 一种基于卷积神经网络的Sentinel-1雷达图像分类方法 | |
CN112597815A (zh) | 一种基于Group-G0模型的合成孔径雷达图像舰船检测方法 | |
CN113408340B (zh) | 基于增强型特征金字塔的双极化sar小型船只检测方法 | |
CN108932455B (zh) | 遥感图像场景识别方法及装置 | |
CN109543585A (zh) | 基于卷积神经网络的水下光学目标检测与识别方法 | |
CN109558803B (zh) | 基于卷积神经网络与np准则的sar目标鉴别方法 | |
CN113393457A (zh) | 一种结合残差密集块与位置注意力的无锚框目标检测方法 | |
CN111126155B (zh) | 一种基于语义约束生成对抗网络的行人再识别方法 | |
CN113128564A (zh) | 一种基于深度学习的复杂背景下典型目标检测方法及系统 | |
Yaohua et al. | A SAR oil spill image recognition method based on densenet convolutional neural network | |
CN114187506B (zh) | 视点意识的动态路由胶囊网络的遥感图像场景分类方法 | |
Zhang et al. | Nearshore vessel detection based on Scene-mask R-CNN in remote sensing image | |
CN116580322A (zh) | 一种地面背景下无人机红外小目标检测方法 | |
CN113505712B (zh) | 基于类平衡损失函数的卷积神经网络的海面溢油检测方法 | |
CN116129280B (zh) | 一种遥感影像雪检测的方法 | |
CN112508106A (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 |