CN105866775B - 一种基于全极化合成孔径雷达图像的海面溢油检测方法 - Google Patents
一种基于全极化合成孔径雷达图像的海面溢油检测方法 Download PDFInfo
- Publication number
- CN105866775B CN105866775B CN201610183160.3A CN201610183160A CN105866775B CN 105866775 B CN105866775 B CN 105866775B CN 201610183160 A CN201610183160 A CN 201610183160A CN 105866775 B CN105866775 B CN 105866775B
- Authority
- CN
- China
- Prior art keywords
- polarization
- covariance matrix
- sea surface
- synthetic aperture
- image
- 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.)
- Expired - Fee Related
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 37
- 230000010287 polarization Effects 0.000 claims abstract description 68
- 239000011159 matrix material Substances 0.000 claims abstract description 39
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000012549 training Methods 0.000 claims abstract description 18
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 16
- 239000013598 vector Substances 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 230000000877 morphologic effect Effects 0.000 claims abstract description 6
- 239000003305 oil spill Substances 0.000 claims description 33
- 239000003921 oil Substances 0.000 claims description 29
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 239000002480 mineral oil Substances 0.000 claims description 7
- 235000010446 mineral oil Nutrition 0.000 claims description 7
- 238000007781 pre-processing Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- NJPPVKZQTLUDBO-UHFFFAOYSA-N novaluron Chemical compound C1=C(Cl)C(OC(F)(F)C(OC(F)(F)F)F)=CC=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F NJPPVKZQTLUDBO-UHFFFAOYSA-N 0.000 claims description 4
- 238000007670 refining Methods 0.000 claims description 4
- 238000012795 verification Methods 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- ODKSFYDXXFIFQN-UHFFFAOYSA-M argininate Chemical compound [O-]C(=O)C(N)CCCNC(N)=N ODKSFYDXXFIFQN-UHFFFAOYSA-M 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000017105 transposition Effects 0.000 claims description 3
- 230000002336 repolarization Effects 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 11
- 230000008569 process Effects 0.000 abstract description 3
- 239000000284 extract Substances 0.000 abstract description 2
- 238000005259 measurement Methods 0.000 abstract 1
- 238000010200 validation analysis Methods 0.000 abstract 1
- 206010027146 Melanoderma Diseases 0.000 description 5
- 238000000605 extraction Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000012075 bio-oil Substances 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000007635 classification algorithm Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000013524 data verification Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000009364 mariculture Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000010913 used oil Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/904—SAR modes
- G01S13/9076—Polarimetric features in SAR
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Artificial Intelligence (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于全极化合成孔径雷达(SAR)图像的海面溢油检测方法，包括以下步骤：对需要分析的极化SAR数据进行预处理，得到全极化SAR协方差矩阵；对协方差矩阵进行精制极化Lee滤波；基于滤波后的协方差矩阵或通过其计算出的Stokes向量，提取六个极化特征，构成特征组合；将有地面验证标签信息的训练样本特征输入极大似然分类器，对分类器进行训练，并优化其参数；将极化特征作为输入，利用极大似然分类器对油膜进行检测与分类；对分类结果进行基于形态学处理，利用实测的验证信息对分类精度进行评估。本发明能够提高海面溢油检测与分类方法的性能，促进全极化SAR在海面溢油监测等实际工程问题中的应用。
Description
技术领域
本发明属于海洋遥感技术领域，特别涉及一种监测海面溢油的雷达卫星遥感方法。
背景技术
海洋溢油是目前最为常见的海洋污染之一，它使海洋环境和生态***受到破坏，威胁海洋生物的生长。溢油事故会使渔业及海产养殖业蒙巨大经济损失，甚至还会通过食物链的传导危害人类的健康。及时、准确的海面溢油监测和早期预警，可以使事故在影响扩大之前得到及早处置，有效地降低事故带来的负面影响。实效性强的溢油监测手段能够为追究肇事船只和企业的责任提供依据，从而为维持海岸带及海洋生态***的健康和可持续发展提供支持。
传统的溢油检测主要依赖于飞机和船舶的海上定点巡逻监视等方法，这些方法不仅费用昂贵，需要大量人力物力，而且效率较低，受到天气、时间等自然条件的约束。与地基、海基观测等传统方式相比，航空、航天遥感极大地扩大了海面监视的范围。随着科技的进步，目前航空及卫星遥感已经成为海面及岸带遥感监测的有效手段。按照传感器的工作波段，可将遥感分类为，紫外遥感，可见光遥感，红外遥感以及微波遥感。其中微波遥感是海洋遥感的最佳手段之一，其主要原因是微波可以穿过大气窗口，对海洋进行全天时、全天候的观测，不受云雨光照等自然条件影响。目前，微波遥感手段中应用最为广泛的技术就是合成孔径雷达(SAR)，由于其具有高分辨率、宽覆盖等优势，已成为对地观测领域的前沿技术，越来越多地被用在全球环境监测中。
利用SAR可以精密地测量海面的粗糙度。由于海面油膜抑制海面毛细波，削弱对雷达后向散射起主要贡献的布拉格散射，因此在SAR图像中形成“暗区”，从而被检测出来。但是，由于SAR的成像原理，图像中会难以避免地存在相干斑噪声，影响油膜区域的提取。更为严重的问题是：一些自然现象，例如海面禁风区、生物油膜、淡水层、波浪阴影等也会在SAR图像中产生相似的“暗区”，为SAR溢油检测的实际应用带来一定的困难。
早期的SAR溢油检测主要基于单极化图像，尽管利用灰度，纹理以及形状等信息可以排除一部分似然物，但这些方法往往难以排除一些虚警。为了解决这一问题，近些年越来越多的研究考虑利用多极化SAR来进行溢油监测。这是因为在不同的极化信号照射下地物会产生不同特性的回波，而通过这些极化信息可以提高目标识别的准确度。已经有许多研究发现，溢油与某些似然物，特别是生物油膜引起的暗斑，具有非常不同的极化散射特征。因此通过利用新一代的SAR卫星提供多极化的探测能力，溢油检测的错误率可以被显著降低。有研究表明，同极化相位差(CPD)能够有效地区分溢油与生物油膜，而从极化相干矩阵中提取出的一些参数如极化熵、alpha角等也能提高溢油分类的性能。极化SAR的应用为油膜分类提供了新的特征集，提高了溢油检测的准确度。
目前基于SAR图像的海面溢油检测算法主要分为三步，即黑斑检测、特征提取、以及油膜与似然物的分类。在黑斑检测上主要有基于全局或者局部自适应的阈值法，以及基于能量泛函的水平集法等。基于对检测出的黑斑的分析，可以提取出绝对灰度、对比度等特征，长宽比、复杂度等形态学参数，以及纹理等特征。利用这些特征，可以利用特征级的监督与非监督分类将矿物油与生物油膜等似然物进行分类。目前常用的分类器有非监督的K均值法，监督的支持向量机(SVM)以及人工神经网络(ANN)等。最近的一些研究证实，在分类中引入极化特征可以有效地提高油膜分类的精度，降低误分的概率。传统的海面溢油检测算法的有效性虽然已被许多研究与实际应用所证实，但是仍然存在较大的改进空间。广泛使用的溢油检测算法大都分具有以上三步的原因是，在溢油与似然物的分类过程中需要提取特征级的信息。而极化特征的引入使得黑斑检测的步骤变得并非必要，因为极化特征提供了对于油膜与似然物的充足区分度。在基于极化特征的油膜检测算法中，基于灰度的黑斑检测或者分割有可能忽视某些非典型或者面积较小的矿物油膜，形成漏警。同时，三步法处理流程较长，需要提取特征级的信息，需要较大的计算量，这给大范围海面溢油近实时检测带来一定的挑战。
发明内容
针对目前利用合成孔径雷达(SAR)卫星图像检测海面溢油存在的似然物难以区分，算法复杂需要处理时间长等技术问题，本发明的目的是提供基于全极化合成孔径雷达图像的海面溢油检测方法，以提高海面溢油检测与分类方法的性能，便于全极化SAR在海面溢油监测等实际工程问题中的应用。
为实现上述目的，本发明采用的技术方案为：
一种基于全极化合成孔径雷达图像的海面溢油检测方法，包括以下步骤：
步骤1：从全极化合成孔径雷达图像中获取需要处理的极化SAR数据，对需要处理的极化SAR数据进行预处理，得到全极化SAR协方差矩阵(C)；
步骤2：对上述步骤得到的协方差矩阵进行精制极化Lee滤波；
步骤3：基于滤波后的协方差矩阵或通过其计算出的Stokes向量，提取六个极化特征，构成特征组合；
步骤4：将有地面验证信息的训练样本特征输入极大似然分类器，对分类器进行训练，并优化其参数；
步骤5：将需要进行溢油检测图像的极化特征作为输入，利用极大似然分类器对油膜进行检测与分类；
步骤6：对分类结果进行基于形态学的处理，去除面积过小的区域及填补闭合区域中的孔洞，利用实测的验证信息对分类精度进行评估。
进一步的，步骤1中，极化SAR数据进行预处理包括格式转换、定标。
进一步的，步骤1中，若得到的是四通道SAR复极化散射矩阵(S)，则通过下述方法转化为协方差矩阵：
其中下标i,j分别代表发射和接收的极化方式，h为水平极化，v为垂直极化；
代表散射向量(假设对于单站***，互易定理成立，即Shv＝Svh)，T代表矩阵转置；
其中上标*T代表共轭转置，〈〉表示用一定大小的窗口(例如5×5)对数据进行空域平均。
进一步的，步骤2中，对协方差矩阵进行精制极化Lee滤波算法包括如下步骤：(1)在SPAN图像上进行边缘模板匹配，选取方向窗口；(2)在该方向窗口内应用局部统计滤波器对协方差矩阵进行滤波处理。
进一步的，步骤3中，基于滤波后的协方差矩阵或通过其计算出的Stokes向量，提取如下六个极化特征构成特征组合：
a)极化度：
其中，g*为Stokes向量(*＝0,1,2,3)；
b)一致性系数：
c)同极化相位差的标准差：
其中arg()表示求复数角度；
d)极化熵：
其中Pi由下式求得：
其中，λi(i＝1,2,3)是极化相干矩阵T的特征值；
e)归一化极化基座高度
f)相关系数
进一步的，步骤4中，采用监督分类的方法，基于极化特征将油膜与海面及其似然物进行分类；训练样本及分类结果都只包含两类，即矿物油膜和非油膜；分类操作前需要获取训练样本，样本应在与待处理图像相似的海况下获得，包含不同的海面目标类型，并经人工解译或者地面调查信息进行验证；对包含训练样本的数据进行1-3步操作，提取极化特征，将这些特征输入极大似然(ML)分类器，对分类器进行训练，并优化其参数。
进一步的，步骤6中，通过形态学处理去除面积过小的区域已及填补闭合区域中的孔洞：利用圆形模板，对分类结果的二值图像先进行闭运算，再进行开运算，再去除面积过小的区域。
本发明的有益效果是：
(1)本发明利用全极化合成孔径雷达图像所提供的海面场景信息，通过基于训练样本的监督分类对海面油膜进行检测，与传统基于单极化图像的海面溢油检测算法相比，能更好地区分油膜和似然物，具有更高的检测性能。
(2)本发明利用极化滤波器对全极化SAR海面场景图像进行预处理，极大地抑制了斑点噪声，解决了其对海面油膜检测与分类算法性能影响的问题。
(3)本发明省略了传统三步法溢油检测流程中的黑斑提取步骤，将似然油膜的提取与油膜和似然物的分类同时进行，提高了算法的效率，缩短了处理时间，同时降低了漏警率。
(4)本发明基于极化SAR图像的协方差矩阵与Stokes向量，选取具有最佳分类性能的6个特征作为特征组合，很好地平衡了分类算法性能与特征维数复杂度的问题，进一步提高了算法的效率。
附图说明
图1为本发明的方法流程图；
图2为实例区域的功率灰度图像，其中，左侧图像中包含生物油膜，右侧图像中包含矿物油；
图3为极化滤波后图像；
图4为极化特征图(归一化显示):a极化基座，b极化熵，c极化度，d相关系数，e一致性系数f同极化相位差的标准差；
图5为极大似然分类结果示意图；
图6为最终油膜检测结果示意图。
具体实施方式
本发明中提出了一种利用合成孔径雷达监测海面溢油的遥感图像快速处理方法，以NASA航天飞机雷达实验中获取的SIR-C图像为例，结合实际数据验证方法的有效性。
具体处理步骤如下：
1、数据获取与预处理
从全极化合成孔径雷达图像中获取需要处理的极化SAR数据，对需要处理的极化SAR数据进行格式转换、定标等预处理(如需要)，得到全极化SAR协方差矩阵(C)。若得到的是复极化散射矩阵(S)，则通过以下方法将其转化为协方差矩阵：
其中下标i,j分别代表发射和接收的极化方式，h为水平极化，v为垂直极化。
代表散射向量(假设对于单站***，互易定理成立，即Shv＝Svh)，T代表矩阵转置；
其中上标*T代表共轭转置，〈〉表示用一定大小的窗口(例如5×5)对数据进行空域平均。
图2是实验所用SAR数据的功率灰度像，其中，左侧图像中包含生物油膜，右侧图像中包含矿物油。
2、精制极化Lee滤波
对协方差矩阵进行精制极化Lee滤波算法主要包括如下步骤：1.在极化总功率SPAN图像上进行边缘模板匹配，选取方向窗口；2.在该方向窗口内应用局部统计滤波器对协方差矩阵进行滤波处理。极化滤波结果见图3。具体方法可参见论文：
Lee J S,Grunes M R,Grandi G De.Polarimetric SA R speckle filteringand its implication for classifica tion[C]//IEEE Trans.on Geosci.RemoteSensing,1999,37(5):2363-2373.
3、极化特征提取
基于滤波后的协方差矩阵或通过其计算出的Stokes向量，提取如下六个极化特征构成特征组合。极化特征的示意图见图4。
a)极化度(Degree of polarization)
极化向量表示为：
这里因此Stokes向量表示为：
进而计算出极化度：
其中，g*为Stokes向量(*＝0,1,2,3)；
b)一致性系数(Conformity coefficient)
c)同极化相位差的标准差(Std.of Co-polarized Phase Difference)
其中arg()表示求复数角度；
d)极化熵(Entropy)
其中Pi由下式求得：
其中，λi(i＝1,2,3)是极化相干矩阵T的特征值；
e)归一化极化基座高度(Normalized Pedestal height)
f)相关系数(Correlation coefficient)
4、极大似然分类
采用监督分类的方法，基于极化特征将油膜与海面及其似然物进行分类。训练样本及分类结果都只包含两类，即矿物油膜和非油膜。分类操作前需要获取一定量的训练样本，样本应在与待处理图像相似的海况下获得，包含不同的海面目标类型，并经人工解译或者地面调查信息进行验证。对包含训练样本的数据进行1-3步操作，提取极化特征，将这些特征输入极大似然(ML)分类器，对分类器进行训练，并优化其参数。
然后从要进行溢油检测图像中提取极化特征，输入极大似然分类器进行油膜检测与分类。分类结果如图5。
根据贝叶斯(Bayes)准则，事件ω的后验概率P(i|ω)，由伴随其发生的事件的先验概率P(i)及其发生的条件概率P(ω|i)计算出：
在本方法中，ω代表图像中提取出的特征组合，即特征向量，而i代表类。根据极大似然分类法则，像元x被分为i类的判决准则是：
X∈i，若对于所有的j≠i有P(i|ω)＞P(j|ω)
5、形态学后处理
最后通过形态学处理去除面积过小的区域已及填补闭合区域中的孔洞：利用圆形模板，对分类结果的二值图像先进行闭运算，再进行开运算，再去除面积过小的区域。最终油膜检测结果如图6。圆形模板以及小面积区域的半径可根据实际情况，由实验或者经验决定。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种基于全极化合成孔径雷达图像的海面溢油检测方法，其特征在于：包括以下步骤：
步骤1：从全极化合成孔径雷达图像中获取需要处理的极化SAR数据，对需要处理的极化SAR数据进行预处理，得到全极化SAR协方差矩阵(C)；
步骤2：对上述步骤得到的协方差矩阵进行精制极化Lee滤波；
步骤3：基于上述步骤滤波后的协方差矩阵或通过其计算出的Stokes向量，提取如下六个极化特征构成特征组合：
a)极化度：
<mrow>
<mi>P</mi>
<mo>=</mo>
<mfrac>
<msqrt>
<mrow>
<msubsup>
<mi>g</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>g</mi>
<mn>2</mn>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>g</mi>
<mn>3</mn>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
<msubsup>
<mi>g</mi>
<mn>0</mn>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
其中，g*为Stokes向量，其中，*＝0,1,2,3；
b)一致性系数：
<mrow>
<mi>C</mi>
<mi>o</mi>
<mi>n</mi>
<mi>f</mi>
<mo>&cong;</mo>
<mfrac>
<mrow>
<mn>2</mn>
<mrow>
<mo>(</mo>
<mi>Re</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
<mo>*</mo>
</msubsup>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>(</mo>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<mn>2</mn>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>)</mo>
</mrow>
</mfrac>
</mrow>
c)同极化相位差的标准差：
<mrow>
<mi>C</mi>
<mi>P</mi>
<mi>D</mi>
<mo>=</mo>
<mi>arg</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>arg</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>arg</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>H</mi>
<mi>H</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>V</mi>
<mi>V</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
其中arg()表示求复数角度；
d)极化熵：
<mrow>
<mi>H</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mn>3</mn>
</munderover>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<msub>
<mi>log</mi>
<mn>3</mn>
</msub>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
</mrow>
其中Pi由下式求得：
<mrow>
<msub>
<mi>P</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mi>j</mi>
</munder>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
</mrow>
其中，λi是极化相干矩阵T的特征值，其中，i＝1,2,3；
e)归一化极化基座高度
<mrow>
<mi>N</mi>
<mi>P</mi>
<mi>H</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>min</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>&lambda;</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<msub>
<mi>&lambda;</mi>
<mn>3</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
</mrow>
f)相关系数
<mrow>
<msub>
<mi>&rho;</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
<mi>l</mi>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>|</mo>
<mfrac>
<mrow>
<mo><</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>></mo>
</mrow>
<mrow>
<mo><</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>></mo>
<mo><</mo>
<msubsup>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
<mn>2</mn>
</msubsup>
<mo>></mo>
</mrow>
</mfrac>
<mo>|</mo>
</mrow>
<mo>;</mo>
</mrow>
步骤4：将有地面验证信息的训练样本特征输入极大似然分类器，对分类器进行训练，并优化其参数；
步骤5：将需要进行溢油检测图像的极化特征作为输入，利用极大似然分类器对油膜进行检测与分类；
步骤6：对分类结果进行基于形态学的处理，去除面积过小的区域及填补闭合区域中的孔洞，利用实测的验证信息对分类精度进行评估。
2.如权利要求1所述的基于全极化合成孔径雷达图像的海面溢油检测方法，其特征在于：步骤1中，极化SAR数据进行预处理包括格式转换、定标。
3.如权利要求1所述的基于全极化合成孔径雷达图像的海面溢油检测方法，其特征在于：步骤1中，若得到的是四通道SAR复极化散射矩阵(S)，则通过下述方法转化为协方差矩阵：
<mrow>
<mi>S</mi>
<mo>=</mo>
<mfenced open = "(" close = ")">
<mtable>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>h</mi>
</mrow>
</msub>
</mtd>
<mtd>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中下标i,j分别代表发射和接收的极化方式，h为水平极化，v为垂直极化；
<mrow>
<mover>
<mi>k</mi>
<mo>&RightArrow;</mo>
</mover>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>,</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>&rsqb;</mo>
</mrow>
<mi>T</mi>
</msup>
</mrow>
代表散射向量，假设对于单站***，互易定理成立，即Shv＝Svh，T代表矩阵转置；
<mrow>
<mi>C</mi>
<mo>=</mo>
<mo><</mo>
<mover>
<mi>k</mi>
<mo>&RightArrow;</mo>
</mover>
<msup>
<mover>
<mi>k</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mo>*</mo>
<mi>T</mi>
</mrow>
</msup>
<mo>></mo>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo><</mo>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>></mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo><</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>></mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo><</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>></mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo><</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>h</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>></mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo><</mo>
<mn>2</mn>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>></mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo><</mo>
<mn>2</mn>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>></mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo><</mo>
<mn>2</mn>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>></mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo><</mo>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
<msubsup>
<mi>S</mi>
<mrow>
<mi>h</mi>
<mi>v</mi>
</mrow>
<mo>*</mo>
</msubsup>
<mo>></mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo><</mo>
<msup>
<mrow>
<mo>|</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>v</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>></mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中上标*T代表共轭转置，<>表示用窗口对数据进行空域平均。
4.如权利要求1所述的基于全极化合成孔径雷达图像的海面溢油检测方法，其特征在于：步骤2中，对协方差矩阵进行精制极化Lee滤波算法包括如下步骤：(1)在极化总功率SPAN图像上进行边缘模板匹配，选取方向窗口；(2)在该方向窗口内应用局部统计滤波器对协方差矩阵进行滤波处理。
5.如权利要求1所述的基于全极化合成孔径雷达图像的海面溢油检测方法，其特征在于：步骤4中，采用监督分类的方法，基于极化特征将油膜与海面及其似然物进行分类；训练样本及分类结果都只包含两类，即矿物油膜和非油膜；分类操作前需要获取训练样本，样本应在与待处理图像相似的海况下获得，包含不同的海面目标类型，并经人工解译或者地面调查信息进行验证；对包含训练样本的数据进行1-3步操作，提取极化特征，将这些特征输入极大似然(ML)分类器，对分类器进行训练，并优化其参数。
6.如权利要求1所述的基于全极化合成孔径雷达图像的海面溢油检测方法，其特征在于：步骤6中，通过形态学处理去除面积过小的区域已及填补闭合区域中的孔洞：利用圆形模板，对分类结果的二值图像先进行闭运算，再进行开运算，再去除面积过小的区域。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610183160.3A CN105866775B (zh) | 2016-03-28 | 2016-03-28 | 一种基于全极化合成孔径雷达图像的海面溢油检测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610183160.3A CN105866775B (zh) | 2016-03-28 | 2016-03-28 | 一种基于全极化合成孔径雷达图像的海面溢油检测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105866775A CN105866775A (zh) | 2016-08-17 |
| CN105866775B true CN105866775B (zh) | 2018-01-02 |
Family
ID=56626058
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610183160.3A Expired - Fee Related CN105866775B (zh) | 2016-03-28 | 2016-03-28 | 一种基于全极化合成孔径雷达图像的海面溢油检测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105866775B (zh) |
Families Citing this family (14)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106291555A (zh) * | 2016-08-24 | 2017-01-04 | 武汉大学 | 一种基于c波段极化sar影像的海面溢油探测方法 |
| CN106845347A (zh) * | 2016-12-19 | 2017-06-13 | 航天恒星科技有限公司 | 海上油膜识别提取方法及*** |
| CN106815559B (zh) * | 2016-12-21 | 2020-04-14 | 中国科学院深圳先进技术研究院 | 一种利用sar数据监测蚝排区域的方法及装置、用户设备 |
| CN107101948B (zh) * | 2017-06-06 | 2019-12-27 | 南京理工大学 | 一种基于偏振吸收特性的海面溢油检测方法 |
| CN107609577B (zh) * | 2017-08-23 | 2020-05-01 | 中国地质大学（武汉） | 一种利用随机森林的极化sar海面油膜提取方法 |
| CN108256471B (zh) * | 2018-01-16 | 2020-08-14 | 北方工业大学 | 一种基于双层分类网的高分辨率合成孔径雷达船只识别方法 |
| CN108960300B (zh) * | 2018-06-20 | 2021-03-02 | 北京工业大学 | 一种基于深度神经网络的城市土地利用信息分析方法 |
| CN109490886A (zh) * | 2018-11-28 | 2019-03-19 | 北京工业大学 | 一种极化合成孔径雷达遥感海面溢油区域精确提取方法 |
| CN109615595B (zh) * | 2018-12-03 | 2019-07-23 | 中国石油大学(华东) | 一种基于双边滤波的水平集sar溢油提取方法 |
| CN109829858B (zh) * | 2019-01-29 | 2022-10-18 | 大连海事大学 | 一种基于局部自适应阈值的船载雷达图像溢油监测方法 |
| CN110188628A (zh) * | 2019-05-13 | 2019-08-30 | 北京遥感设备研究所 | 一种基于深度学习模型的rd图像舰船目标识别方法 |
| CN111025291B (zh) * | 2019-11-06 | 2021-07-16 | 中国石油大学(华东) | 一种基于全极化sar新特征的海洋溢油检测方法 |
| CN113378766B (zh) * | 2021-06-25 | 2022-04-05 | 南通大学 | 一种基于合成孔径雷达的海上大规模风力发电站监测*** |
| CN113516657B (zh) * | 2021-09-14 | 2021-12-10 | 中国石油大学(华东) | 一种基于自适应权重的全极化sar影像海面溢油检测方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6646593B1 (en) * | 2002-01-08 | 2003-11-11 | Science Applications International Corporation | Process for mapping multiple-bounce ghosting artifacts from radar imaging data |
| CN102096070B (zh) * | 2010-11-16 | 2013-04-03 | 苏州科技学院 | 一种全极化合成孔径雷达海面溢油检测方法 |
| CN102999897A (zh) * | 2011-09-19 | 2013-03-27 | 香港中文大学 | 基于sar图像检测海面溢油的方法和装置 |
-
2016
- 2016-03-28 CN CN201610183160.3A patent/CN105866775B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN105866775A (zh) | 2016-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105866775B (zh) | 一种基于全极化合成孔径雷达图像的海面溢油检测方法 | |
| Yang et al. | Ship detection from optical satellite images based on sea surface analysis | |
| Singha et al. | Satellite oil spill detection using artificial neural networks | |
| Zhu et al. | Oil spill contextual and boundary-supervised detection network based on marine SAR images | |
| Mattyus | Near real-time automatic marine vessel detection on optical satellite images | |
| Fan et al. | Oil spill monitoring based on SAR remote sensing imagery | |
| CN109427055B (zh) | 基于视觉注意机制和信息熵的遥感图像海面舰船检测方法 | |
| Fan et al. | Floating raft aquaculture information automatic extraction based on high resolution SAR images | |
| Park et al. | Oil spill detection from PlanetScope satellite image: application to Oil spill Accident near Ras Al Zour area, Kuwait in august 2017 | |
| Hou et al. | End-to-end automatic ship detection and recognition in high-resolution Gaofen-3 spaceborne SAR images | |
| Akbari et al. | Iceberg detection in open and ice-infested waters using C-band polarimetric synthetic aperture radar | |
| Skrunes et al. | Characterization of low-backscatter ocean features in dual-copolarization SAR using log-cumulants | |
| He et al. | Ship detection without sea-land segmentation for large-scale high-resolution optical satellite images | |
| Ajadi et al. | Oil spill detection in synthetic aperture radar images using Lipschitz-regularity and multiscale techniques | |
| Li et al. | Detection and monitoring of oil spills using moderate/high-resolution remote sensing images | |
| Chen et al. | A novel scheme for extracting sea surface wind information from rain-contaminated x-band marine radar images | |
| Makedonas et al. | Vessel classification in cosmo-skymed SAR data using hierarchical feature selection | |
| Wang et al. | An improved semantic segmentation model based on SVM for marine oil spill detection using SAR image | |
| Li et al. | Detection of oil spills based on gray level co-occurrence matrix and support vector machine | |
| Wei et al. | Texture feature analysis in oil spill monitoring by SAR image | |
| CN111767806B (zh) | 基于Attribute的极窄脉冲雷达船只目标鉴别方法 | |
| Akbari et al. | Iceberg detection in open water and sea ice using C-band radar polarimetry | |
| Solberg et al. | Algorithms for oil spill detection in Radarsat and ENVISAT SAR images | |
| Yang et al. | A novel oil spill detection method from synthetic aperture radar imageries via a bidimensional empirical mode decomposition | |
| Gao | Design and implementation of marine automatic target recognition system based on visible remote sensing images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP02 | Change in the address of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: 210000 No. 219 Ning six road, Jiangbei new district, Nanjing, Jiangsu Patentee after: Nanjing University of Information Science and Technology Address before: The Olympic Avenue in Jianye District of Nanjing city of Jiangsu Province, No. 69 210019 Patentee before: Nanjing University of Information Science and Technology |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180102 |