CN105046707B - 基于n阶多项式函数拟合海杂波的SAR图像船只检测方法 - Google Patents

基于n阶多项式函数拟合海杂波的SAR图像船只检测方法 Download PDF

Info

Publication number
CN105046707B
CN105046707B CN201510409676.0A CN201510409676A CN105046707B CN 105046707 B CN105046707 B CN 105046707B CN 201510409676 A CN201510409676 A CN 201510409676A CN 105046707 B CN105046707 B CN 105046707B
Authority
CN
China
Prior art keywords
sea clutter
fitting
sar image
value
sea
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
Application number
CN201510409676.0A
Other languages
English (en)
Other versions
CN105046707A (zh
Inventor
郎海涛
石洪基
谢添慧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Chemical Technology
Original Assignee
Beijing University of Chemical Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Chemical Technology filed Critical Beijing University of Chemical Technology
Priority to CN201510409676.0A priority Critical patent/CN105046707B/zh
Publication of CN105046707A publication Critical patent/CN105046707A/zh
Application granted granted Critical
Publication of CN105046707B publication Critical patent/CN105046707B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/0002—Inspection of images, e.g. flaw 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

Abstract

本发明公开基于n阶多项式函数拟合海杂波的SAR图像船只检测方法，主要解决传统的海杂波拟合方式应用于高分辨高海况的SAR图像时，整体检测性能差的问题。本发明采用n阶多项式函数拟合海杂波分布，综合考虑拟合性能及拟合效率。基于多项式函数拟合海杂波分布的方法相较于传统分布具有更好的拟合优度(GoF‑Goodness of Fitting)，将拟合的海杂波分布与恒虚警率(CFAR‑Constant False Alarm Rate)方法相结合，继而得到更好的目标检测效果。本发明可应用于渔业监管、船只检测等方面。

Description

基于n阶多项式函数拟合海杂波的SAR图像船只检测方法
技术领域
本发明属于遥感图像处理技术领域，具体涉及图像处理、海杂波的拟合、CFAR检测算法。可应用于渔业监管、船只监测等方面。
背景技术
在资源愈加匮乏的现代，海洋维权愈发重要。合成孔径雷达(SAR)作为一种高分辨率成像传感器，具有全天时和全天候观测的能力，已成为海洋监测的有效手段。
在SAR图像中，由于船只主要由金属制成，后向散射系数很大，对应的亮度较高；而海洋的后向散射系数很小，因此海洋在SAR图像中灰度很小。根据灰度值的不同，而将船只目标与海杂波区分开来[1](参见唐沐恩,林挺强等.遥感图像中舰船检测方法综述[J].计算机应用研究,2011(28):29-36)。但在不同的天气、风速下，海况相差很大，SAR图像海面背景与船只目标的对比度随之变化，因此呈现在SAR图像中的海杂波效果也有很大的差别，针对这些复杂情况，在检测过程中需要有一个自适应的检测方法，现应用得较普遍的则是恒虚警率(CFAR)检测方法。
这个算法的核心思想是在保证虚警率为常数的同时，根据虚警率和SAR图像海杂波的统计特性(即海洋杂波的概率密度函数(PDF-Probability Distribution Function))计算得到检测船只目标的阈值。再将所得阈值与SAR图像像素值进行比较，大于阈值的则为所求的船只目标，小于阈值的则为海洋背景。
由于使用不同的海杂波模型可得到具有不同形式和检测效果的CFAR检测器，因此该算法的关键在于海杂波PDF的建模，海杂波模型与实际海况的拟合度越高，得到的检测效果越好。目前已有多种海杂波拟合分布，如加拿大海洋监测系统(OMW)中使用的监测算法就是一个使用K分布作为海杂波分布模型的CFAR检测算法[2](参见VACHON P W；ADLAKHA P；EDEL H，el al.Canadian progress toward marine and coastal applications ofsynthetic aperture radar 2000(01))；基于高斯分布的CFAR检测算法[3](参见NovakLM.，Halversen S.D.Effects of Polarization and Resolution on SAR ATR[J].IEEEtrans.on Aerospace and Electronic Systems.1997，33(1).102-115)；基于伽马分布的CFAR算法[4](参见Principe J.C.,Radisavljevic A.Target prescreening based on aquadratic gamma discriminator[J].IEEE Transactions on Aerospace andElectronic Systems，34(3)：706-715)等。而随着SAR图像分辨率的进一步提高，以及对复杂海况下SAR图像研究的展开，这些传统海杂波分布模型的检测结果虚检率和漏检率增高，需要针对高分辨率SAR图像构建更好的拟合模型。
发明内容
本发明的目的在于针对上述已有技术的不足，提出一种基于多项式函数拟合海杂波分布的SAR图像船只检测方法。多项式函数拟合海杂波的主要优点在于对数据的拟合比较灵活，线性地计算复杂度，拟合过程简单。
实现本发明的技术方案为：输入SAR图像，基于n阶多项式函数拟合其海杂波分布，根据拟合的海杂波的概率密度分布计算累积分布密度，并结合预先设定的恒虚警率，计算检测阈值，将船只目标从海洋背景中分割出来。
其具体实现步骤包括如下：
①获得SAR图像数据值分布的概率密度函数f1(x)。SAR图像可以进行滤波处理也可以不进行滤波处理；SAR图像数据值可以取幅度值或强度值。
②对①所述的概率密度函数进行以a为底的对数运算，得到f2(x)。即f2(x)＝logaf1(x)。
对a不做特别要求，本发明推荐a＝10。
③n阶多项式函数共有n+1个未知量，为了使采取的样本足够进行函数拟合，拟在f2(x)函数中均匀采样m个数据点，m取远大于n的值，如m>10n,进一步说，m的取值范围：m>80。
采样数据点的坐标为(x1,f2(x1))、(x2,f2(x2))、…、(xm,f2(xm))。用n阶多项式函数拟合f2(x)：其中，pi为多项式函数的参数，x为图像的灰度值，8位SAR图像取值范围为[0,255]，16位SAR图像取值范围为[0,65535]。构造超定方程如下：
联立求解，得到各项系数p1,p2,…,pn,pn+1的值。得到拟合的函数f3(x)。
⑤根据拟合的海杂波的概率密度分布f4(x)计算累积分布密度F(x),
⑥当虚警率Pfa取某一确定常数时，通过Pfa＝1-F(Tn)，可以求得阈值Tn。Pfa取值取决于对最后检测结果精度的要求，一般取0.0001对于检测区域中的数据值，若大于阈值Tn，则判定为船只目标，否则为海杂波。
若需检测的SAR图像噪点过多，或对结果检测精度有着较严格的要求，可以将图像先进行滤波处理；一般情况下，为节省时间，也可不进行滤波处理。
附图说明
图1是本发明的流程示意图；
图2是本发明具体实施例中所使用的Radarsat-2高分辨SAR图像；
图3是本发明具体实施例中检测结果；
图4是本发明具体实施例中检验结果与专家解译结果对比图。
具体实施方式
参照图1，本发明的具体实施步骤如下：
步骤1：输入图像(图2)，进行归一化处理，将图像像素值归一化到[0,1]。
本发明具体实施例一中SAR图像的来源及详细信息参见下表：
获得SAR图像数据值分布的概率密度函数f1(x)。
步骤2：对f1(x)进行以10为底的对数运算，得到f2(x)，f2(x)＝log10f1(x)。
步骤3：在f2(x)中均匀采样100个数据点。采样数据点的坐标为
(x1,f2(x1))、(x2,f2(x2))、…、(x100,f2(x100))。本实施例中采用6阶多项式函数拟合f2(x)：
构造超定方程如下：
联立求解，得到各项系数p1,p2,…,p5,p6,p7的值：
p1＝1383
p2＝607.7
p3＝-8922
p4＝1861
p5＝1.42e+04
p6＝-1.083e+04
p7＝2284
得到拟合的函数f3(x)。
步骤4：本发明采用R-square评价拟合函数的拟合优度。
其中，SSR(Sum of squares of the regression)即预测数据与原始数据均值之差的平方和，公式如下:SST(Total sum of squares)即原始数据和均值之差的平方和，公式如下:R-square(确定系数)越接近1，表明方程的变量对y的解释能力越强，这个模型对数据拟合得也较好。根据上述公式求得：R-square＝0.9829。
步骤5：通过以a为底数、f3(x)为指数的指数运算，得到拟合后的海杂波概率密度函数f4(x)。即根据拟合的海杂波的概率密度分布f4(x)计算累积分布密度F(x)。
步骤6：设虚警率Pfa＝0.0001，通过Pfa＝1-F(Tn)，可以求得阈值Tn＝0.7217。对于检测区域中的数据值，若大于阈值Tn，则判定为船只目标，否则为海杂波。
步骤7：将所得的结果与专家解译结果对比(图4：正确检测目标用圆圈标注，漏检目标用三角形标注，虚警目标用正方形标注)，得知检验结果。其中，算法性能根据正确检测目标数，漏检目标数，虚警目标数进行评价，参考品质因数FoM指标参数，其定义为：其中：Ntt为检测结果中正确的检测目标数，Nfa为虚警目标数，Ngt为实际的目标数。FoM越趋近于100％，证明检测效果越好。最终检测结果为：
参照优秀的FoM参数，可以证明，该发明对船只目标的检测性能良好。

Claims (2)

1.一种基于多项式函数拟合海杂波的SAR图像船只检测方法，其特征在于步骤如下：
2.1获得SAR图像数据值分布的概率密度函数f1(x)；x为图像的灰度值，8位SAR图像取值范围为[0,255]，16位SAR图像取值范围为[0,65535]；
2.2对2.1所述的概率密度函数进行以a为底的对数运算，得到f2(x)；
2.3采用n阶多项式函数拟合f2(x)，得到拟合函数f3(x)，其中pi为拟合系数；n取大于等于4小于等于7的自然数：
2.4进行以a为底，f3(x)为指数的指数运算，a的值与步骤2.2中的相同；得到拟合后的海杂波概率密度函数f4(x)；
2.5根据拟合的海杂波的概率密度分布f4(x)计算累积分布密度当虚警率Pfa取某一确定常数时，通过Pfa＝1-F(Tn)，求得阈值Tn；Pfa取值取决于对最后检测结果精度的要求，对于检测区域中的数据值，若大于阈值Tn，则判定为船只目标，否则为海杂波。
2.根据权利要求1所述的方法，其特征在于：底数a取10。
CN201510409676.0A 2015-07-13 2015-07-13 基于n阶多项式函数拟合海杂波的SAR图像船只检测方法 Active CN105046707B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510409676.0A CN105046707B (zh) 2015-07-13 2015-07-13 基于n阶多项式函数拟合海杂波的SAR图像船只检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510409676.0A CN105046707B (zh) 2015-07-13 2015-07-13 基于n阶多项式函数拟合海杂波的SAR图像船只检测方法

Publications (2)

Publication Number Publication Date
CN105046707A CN105046707A (zh) 2015-11-11
CN105046707B true CN105046707B (zh) 2018-06-26

Family

ID=54453225

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510409676.0A Active CN105046707B (zh) 2015-07-13 2015-07-13 基于n阶多项式函数拟合海杂波的SAR图像船只检测方法

Country Status (1)

Country Link
CN (1) CN105046707B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10594439B2 (en) 2017-09-08 2020-03-17 Huawei Technologies Co., Ltd. Channel encoding method and apparatus in wireless communications to output a polar encoded bit sequence
CN115378546A (zh) * 2017-09-08 2022-11-22 华为技术有限公司 编译码方法及装置
CN112764030B (zh) * 2020-12-25 2023-09-08 中国人民解放军空军工程大学 基于电磁与统计混合模型的海面sar图像杂波生成方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102968799A (zh) * 2012-12-12 2013-03-13 北京航空航天大学 一种基于积分图像的快速acca-cfar sar图像目标检测方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102968799A (zh) * 2012-12-12 2013-03-13 北京航空航天大学 一种基于积分图像的快速acca-cfar sar图像目标检测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
水陆交界区域SAR图像舰船目标检测技术;阚瀛芝;《中国优秀硕士学位论文全文数据库 信息科技辑》;20141215;第I136-145页 *

Also Published As

Publication number Publication date
CN105046707A (zh) 2015-11-11

Similar Documents

Publication Publication Date Title
Xia et al. Infrared small target detection based on multiscale local contrast measure using local energy factor
Manessa et al. Satellite-derived bathymetry using random forest algorithm and worldview-2 imagery
CN108171193B (zh) 基于超像素局部信息度量的极化sar舰船目标检测方法
CN103971127B (zh) 一种前视雷达成像海面目标关键点检测识别方法
CN102096824A (zh) 基于选择性视觉注意机制的多光谱图像舰船检测方法
CN110889843B (zh) 基于最大稳定极值区域的sar图像舰船目标检测方法
CN108133468B (zh) 自适应参数增强和尾迹辅助检测的恒虚警率舰船检测方法
Fan et al. Floating raft aquaculture information automatic extraction based on high resolution SAR images
CN105046707B (zh) 基于n阶多项式函数拟合海杂波的SAR图像船只检测方法
CN113536963B (zh) 基于轻量化yolo网络的sar图像飞机目标检测方法
Zhang et al. Sparse coding of 2D-slice Zernike moments for SAR ATR
Yang et al. Evaluation and mitigation of rain effect on wave direction and period estimation from X-band marine radar images
CN105046706B (zh) 基于有理多项式函数拟合海杂波的sar图像船只检测方法
Wang et al. A PolSAR ship detector based on a multi-polarimetric-feature combination using visual attention
Chen et al. A novel scheme for extracting sea surface wind information from rain-contaminated x-band marine radar images
Potapov Strategic directions in synthesis of new topological radar detectors of low-contrast targets against the background of high-intensity noise from the ground, sea and precipitations
CN108710816B (zh) 一种基于概率融合的极化sar图像建筑物检测方法
Pelliza et al. Optimal Canny’s parameters regressions for coastal line detection in satellite-based SAR images
CN109031294B (zh) 基于相似性检验的极化sar舰船目标检测方法
Lang et al. A new scattering similarity based metric for ship detection in polarimetric synthetic aperture radar image
CN104992159A (zh) 基于高斯函数拟合海杂波的sar图像船只检测方法
Fein-Ashley et al. Benchmarking Deep Learning Classifiers for SAR Automatic Target Recognition
CN108830290B (zh) 基于稀疏表示和高斯分布的sar图像分类方法
Stasolla et al. Combining wavelets and mathematical morphology to detect changes in time series
Mehdizadeh et al. Exploitation of sub-look analysis and polarimetric signatures for ship detection in PolSAR data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Lang Haitao

Inventor after: Shi Hongji

Inventor after: Xie Tianhui

Inventor before: Lang Haitao

Inventor before: Xie Tianhui

GR01 Patent grant
GR01 Patent grant