CN113158806B - 一种基于OTD_Loglogistic的SAR数据海洋目标检测方法 - Google Patents

一种基于OTD_Loglogistic的SAR数据海洋目标检测方法 Download PDF

Info

Publication number
CN113158806B
CN113158806B CN202110311348.2A CN202110311348A CN113158806B CN 113158806 B CN113158806 B CN 113158806B CN 202110311348 A CN202110311348 A CN 202110311348A CN 113158806 B CN113158806 B CN 113158806B
Authority
CN
China
Prior art keywords
target
model
ocean
detection
logistic
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
CN202110311348.2A
Other languages
English (en)
Other versions
CN113158806A (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.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and 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 Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN202110311348.2A priority Critical patent/CN113158806B/zh
Publication of CN113158806A publication Critical patent/CN113158806A/zh
Application granted granted Critical
Publication of CN113158806B publication Critical patent/CN113158806B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
    • 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/04—Architecture, e.g. interconnection topology
    • G06N3/047—Probabilistic or stochastic 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/04—Architecture, e.g. interconnection topology
    • G06N3/048—Activation functions
    • 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
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Computational Linguistics (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Astronomy & Astrophysics (AREA)
  • Remote Sensing (AREA)
  • Multimedia (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种基于OTD_Loglogistic的SAR数据海洋目标检测方法，属于海洋目标检测领域，采用基于深度学习模型的初检和基于Loglogistic的CFAR精检相结合的检测方法，首先构建了OceanTDA9轻量级深度学习模型，基于该深度学习模型进行海洋目标初检，然后在此基础上采用基于Loglogistic模型的CFAR方法进行海洋目标精细检测，提取海洋目标特征。本发明将基于深度学习的海洋目标初检和基于Loglogistic的CFAR方法结合，克服了滑动窗口CFAR需要逐一统计图像中的像素计算耗费长、无法对图像边缘的像素进行处理造成边缘处漏检等缺陷，提高了海洋目标的检测效率和精度。

Description

一种基于OTD_Loglogistic的SAR数据海洋目标检测方法
技术领域
本发明属于海洋目标检测领域，具体涉及一种基于OTD_Loglogistic的SAR数据海洋目标检测方法。
背景技术
近年来，神经网络已应用在海洋目标检测上，但在应用时深度神经网络存在维数灾难的缺陷，会降低检测速度。平静的海面海洋杂波可以近似运用高斯模型建模，然而在风浪、潮汐等复杂海况下，SAR后向散射概率分布出现长尾，不再符合高斯分布。传统CFAR算法基于SAR数据中海杂波服从高斯分布的假定，但在复杂海况中，海杂波分布极其复杂，此假定不具有科学性和普遍性。而且，CFAR需要逐一统计图像中的像素，计算耗费长，滑动窗口无法对图像边缘的像素进行处理，会造成边缘处漏检的缺陷。
发明内容
为了解决上述问题，本发明提出了一种基于OTD_Loglogistic的SAR数据海洋目标检测方法，采用基于深度学习模型的初检和基于Loglogistic的CFAR精检相结合的方法，克服了现有技术的不足，提高了海洋目标检测效率和精度。
本发明的技术方案如下：
一种基于OTD_Loglogistic的SAR数据海洋目标检测方法，采用基于深度学习模型的初检和基于Loglogistic的CFAR精检相结合的检测方法，首先，构建OceanTDA9轻量级深度学习模型，采用所构建的轻量级深度学习模型进行海洋目标初检，在此基础上采用基于Loglogistic模型的CFAR方法进行海洋目标精细检测，提取海洋目标特征。
优选地，OceanTDA9模型包含4个卷积层、1个卷积组和3个全连接层；每个卷积层的组织形式是Convolution2D-ReLU-Dropout-Maxpooling；卷积组的组织形式是(Convolution2D-ReLU-Dropout)*2-Maxpooling；3个全连接层中，前两个全连接层的组织形式是Dense-ReLU-Dropout，最后一个全连接层的组织形式是Dense。
优选地，OceanTDA9模型的训练过程为：
S101.对海洋目标样本特征数据x进行归一化处理，采用的归一化公式如下：
Figure BDA0002989844640000011
其中，μ为特征的期望，σ为方差；
S102.配置训练优化算法、损失函数、监控参数；
S103.加载数据集，验证数据的合法性；
S104.将数据集分为多个周期，每个周期分为多个批次；
S105.将每个周期按批次打混索引，进行批次循环；
S106.计算交叉熵，采用反向传播算法和所改进的梯度下降算法WinR-Adagrad，以0.01的学习速率不断地修改变量以最小化交叉熵；
S107.对批次指定的样本数进行训练，学习权重、偏置，计算损失、精度直到一个周期的批次都循环结束，再进入下一个周期进行批循环；
S108.直到循环完所有周期，保存模型及训练损失、精度。
优选地，OceanTDA9模型损失函数采用目标分类和模型预测分类之间的交叉熵，其公式为：
Figure BDA0002989844640000021
式中，y-i是第i个样本的输入值，hθ(xi)是第i个样本x的输出值，参数θ0的初值设置为0.1，θ1的初值设置为标准差为0.1的正态分布浮点数。
优选地，海洋目标精细检测采用Loglogistic概率密度函数来拟合复杂海况的海杂波分布，突出海杂波的长尾特征，在此基础上采用CFAR方法进行海洋目标特征提取。
优选地，海洋目标精细检测的原理和过程为：
S201.计算初检保存在数组x中每个像素的频数绘制直方图；
S202.提取直方图中最小像素到最大像素间有效像素及对应的频数作为拟合数据并保存到数组x2、y2中；
S203.调用Loglogistic复杂海况的海洋杂波概率密度函数式
Figure BDA0002989844640000022
拟合x2、y2，其中，x>0，α>0，β>0；
求出拟合曲线的α、β参数；
S204.计算拟合优度，包括卡方和绝对值误差；
S205.计算拟合曲线有效像素的积分
Figure BDA0002989844640000023
当积分值达到阈值时，输出临界像素Xtb；
S206.调用目标检测模型PixelsVal＞Xtb，将疑似目标保存到数组row2中；
S207.对数组row2图像窗口中疑似目标像素聚簇，并保存各个目标簇的参数，包括目标簇中心坐标、长度、宽度、像素总数及倾斜角。
本发明所带来的有益技术效果：
通过构建的轻量级卷积神经网络——OceanTDA9，对海洋硬目标进行初检，克服了深度神经网络的维数灾难，提高了海洋图像中目标检测的速度；对复杂海况海杂波进行建模，采用Loglogistic分布拟合复杂海况海杂波的分布，拟合效果更好，突出了海杂波的长尾分布特征，并基于Loglogistic模型采用CFAR进行海洋目标特征提取，提高了精度；
采用基于深度学习的海洋目标初检和基于Loglogistic的CFAR方法精简相结合的方式，克服了滑动窗口CFAR需要逐一统计图像中的像素计算耗费长、无法对图像边缘的像素进行处理造成边缘处漏检等缺陷，提高了海洋目标的检测效率和精度。
附图说明
图1是本发明轻量级深度学习模型OceanTDA9的结构图；
图2是本发明OceanTDA9模型训练流程图；
图3是本发明海洋目标精检流程图；
图4是本发明基于OTD_Loglogistic的海洋目标检测流程图；
图5是本发明多核并行OTD_Loglogistic目标并行检测CPU-t图。
具体实施方式
下面结合附图以及具体实施方式对本发明作进一步详细说明：
一、轻量级深度学习模型构建
1、OceanTDAx系列模型构建
在构建轻量级深度学习模型时，首先设计了轻量级卷积神经网络模型OceanTDAx系列，OceanTDAx系列包括4种模型，分别为OceanTDA2、OceanTDA4、OceanTDA9、OceanTDA16。初步实验OceanTDA9模型检测效果最好。
OceanTDA9模型的结构如图1所示，包含4个卷积层、1个卷积组和3个全连接层，前4个卷积层分别为Conv2D_1、Conv2D_2、Conv2D_3、Conv2D_4，每个卷积的形式一样，都是Convolution2D-ReLU-Dropout-Maxpooling；中间1个卷积组是Conv2D_g，组织形式是(Convolution2D-ReLU-Dropout)*2-Maxpooling；最后3个是密集全连接Dense层，分别为Dense_1、Dense_2、dense_3，其中前两个Dense，每组都是Dense-ReLU-Dropout，最后一个全连接Dense层仅有Dense。
卷积过程是使用一个卷积核W，在每层像素矩阵X上不断按滑动步长stride扫描下去，每次扫到的数值X会和卷积核W中对应位置的数值进行相乘，然后与偏值b相加求和，生成一个新的矩阵XW+b。卷积核W相当于卷积操作中的一个过滤器，用于提取图像的特征，特征提取后会得到一个特征图。卷积核W里面的每个值就是训练模型过程中的神经元参数，即权重，训练开始时赋予W随机的初始值，训练网络过程中，网络会通过后向传播不断更新这些参数值，通过loss损失函数来评估参数值，直到取得最佳的参数值。OceanTDA9模型开始输入的是长28个像素、宽28个像素的单通道图像，卷积核的大小是3×3，滑动步长stride的大小为2。
OceanTDA9模型中的卷积层和全连接层的激活函数都使用ReLU(RectifiedLinear Units，整流线性单元)激活函数，ReLU激活函数应用非饱和激活函数f(x)＝max(0,x)，增加了决策函数和整个网络的非线性特性，而不影响卷积层的感受野。其他函数也可以用来增加非线性，例如饱和双曲正切函数f(x)＝tanh(x)，f(x)＝|tanh(x)|和sigmoid函数f(x)＝(1+e-x)-1。与其它激活函数相比，采用ReLU激活函数的训练神经网络速度快几倍，而且对泛化精度没有显著影响。
OceanTDA9模型在训练过程中，按照一定的比例将网络中的神经元进行丢弃，可以防止模型训练过拟合的情况。OceanTDA9模型中Dropout都设置成0.2。
卷积操作后提取到的特征信息，相邻区域会有相似特征信息，如果全部保留这些特征信息会存在信息冗余，增加计算难度。通过池化可将这些相似特征信息相互替代，不断减小数据的空间大小，参数的数量和计算量会有相应的下降，这在一定程度上控制了过拟合。池化操作相当于降维操作，有最大池化(Max pooling)和平均池化，根据海洋目标检测特点，OceanTDA9模型采用核尺寸为2×2，滑动步长stride为2的Max pooling，每次池化操作后，矩阵的长宽都降低一半。
Flatten将池化后的数据拉开，变成一维向量来表示，方便输入到全连接网络。OceanTDA9模型卷积、池化后，经Flatten平整化后为一维512个特征的向量，进入全连接层。
OceanTDA9模型最后做3层Dense全连接层，特征由512个神经元维压缩到64个，ReLU又接Dropout 0.2过渡，并再次用一个包含64个神经元的全连接Dense作为缓冲后进入全连接层Dense进一步压缩到2个神经元输入损失函数Loss层中的Softmax进行分类。
损失函数层(Loss Layer)用于决定训练过程如何来“惩罚”网络的预测结果和真实结果之间的差异，它通常是网络的最后一层。不同类型的任务采用不同的损失函数，SoftMax交叉熵损失函数常常被用于在K个类别中选出一个，而Sigmoid交叉熵损失函数常常用于多个独立的二分类问题，欧几里德损失函数常常用于结果取值范围为任意实数的问题。OceanTDA9模型确定损失函数采用目标分类和模型预测分类之间的交叉熵，其公式为：
Figure BDA0002989844640000051
式中，y-i是第i个样本的输入值，hθ(xi)是第i个样本x的输出值，即yi。参数θ0的初值设置为0.1，θ1的初值设置为标准差为0.1的正态分布浮点数。
对OceanTDAx系列的4种模型进行对比实验，分别得到损失与批次的相关性曲线、精度与批次的相关性曲线。OceanTDA9模型训练次数为275000次，每次100个样本图像，模型训练到150000次左右接近拐点，随后损失下降缓慢，此时损失是0.02左右，耗时13982秒，最后50次模型损失均值是0.0131，标准差是0.00007；OceanTDA16模型训练次数也是275000次，其它2个模型的训练次数设置为550000次。模型训练开始时波动较大，随着时间推移，模型损失波动逐渐减少。从损失与批次相关性曲线可以看出损失下降由慢到快依次是2层、4层、16层和9层的OceanTDAx模型，综合效果表现最好的是9层模型OceanTDA9。再比较模型精度，从OceanTDAx系列模型训练精度相关性曲线中可以看出9层模型OceanTDA9训练精度为0.9973，训练精度最好。综合损失和精度，模型OceanTDA9检测效果最好，所以本发明最终选择OceanTDA9模型对海洋目标进行初检。
2、OceanTDA9模型训练
首先，要对样本数据x进行特征归一化处理。样本特征的取值范围不同，这可能导致迭代较慢，为了避免此弊端，需要对特征数据进行归一化处理。采用的归一化公式如下：
Figure BDA0002989844640000052
其中，μ为特征的期望，σ为方差。这样特征的新期望为0，新方差为1，迭代次数可以大大加快。
然后，本发明已确定训练的损失函数为目标分类和模型预测分类之间的交叉熵，采用反向传播算法和所改进的梯度下降算法WinR-Adagrad，以0.01的学习速率不断地修改变量以最小化交叉熵。模型具体训练流程如图2所示，训练的关键点在于存在回调函数。回调函数callbacks.on_train_begin()在训练开始时调用，指定logs保存的目录，callbacks.on_epoch_begin()在每个周期开始时调用，记录本次epoch的历史信息，在每个周期的循环中，主要是对每个epoch进行循环，其中核心针对每个batch的处理，按照batch批次打混索引，得到一个批次的索引后，进入batch批次循环。回调函数callbacks.on_batch_begin()的logs包含size，即当前batch的样本数，batch批次循环对指定的样本数进行训练，学习权重、偏置，计算损失、精度等。
模型训练好后，就可调用evaluate()方法在测试数据集上测试评估模型的性能。
loss,accuracy＝model.evaluate(X_test,y_test)
evaluate()方法验证数据的有效性后，首先采用tf.argmax()预测正确的标签，tf.argmax()能返回在一个张量里沿着某条轴的最高条目的索引值，例如tf.argmax(y_ocean,1)是模型认为每个输入最有可能对应的那些标签，而tf.argmax(y_,1)代表正确的标签。然后用tf.equal来检测模型预测是否与真实标签匹配，一般函数会返回一组布尔值，为了确定正确预测项的比例，可以把布尔值转换成浮点数，然后取平均值。损失采用分类交叉熵的平均值做为每个批次的总损失。
correct_prediction＝tf.equal(tf.argmax(y_ocean,1),tf.argmax(y_,1))
accuracy＝tf.reduce_mean(tf.cast(correct_prediction,tf.float32))
categorical_crossentropy＝tf.reduce_mean(tf.nn.softmax_cross_entropy_with_logits(labels＝y_,logits＝y_ocean))
evaluate()返回模型配置时需要返回的损失loss和精度accuracy，可将这些简要信息保存到文件中，通过可视化软件进行可视化，以便利用这些简要信息进一步评估模型。同时训练脚本会为所有学习变量计算其移动均值，评估脚本则直接将所有学习到的模型参数替换成对应的移动均值，这一替代方式可以在评估过程中改善模型的性能。
采用改进优化算法WinR-Adagrad的9层OceanTDA9模型，学习速率设置为0.01，间隔距离tn设置为10，修正系数delta设置为0.6，训练8250次后的模型，对测试数据集的10000个样本图像进行测试，损失为0.0766，精度为99.83％，用时3.345秒。
模型训练好后，通过Gfile()方法将模型保存成pb(Protocol Buffers)格式的文件，便于海洋硬目标检测时调用。
二、基于OceanTDA9模型的海洋目标初检
为了进一步验证本发明OceanTDA9模型的优越性，分别对OceanTDA9、OceanTDvgg、OceanTDcnn模型进行海洋目标初检实验，得到损失_批次、精度_批次的关系曲线。OceanTDA9模型训练结束时损失是0.0625，精度为0.9996，耗时687秒，最后50次模型损失均值是0.2091，标准差是0.00029，精度平均值为0.9214，标准差为0.00015；OceanTDvgg模型训练结束时损失是0.0563，精度为0.9985，耗时1244秒，最后50次模型损失均值是0.2101，标准差是0.00032，精度平均值为0.9212，标准差为0.00015；OceanTDcnn模型训练结束时损失是0.0896，精度为0.9898，耗时540秒，最后50次模型损失均值是0.1930，标准差是0.00027，精度平均值为0.9284，标准差为0.00019。
所构建的OceanTDA9模型训练后共检测到钻井平台疑似目标90个。目标漏检数为0个，检出率为100％，虚警率为10.9％，用时2.36秒，SAR图像检测能力约58.5km2/s。满足无漏检并保持10％左右的虚警率的检测需求。
三、基于Loglogistic分布的海洋背景建模
平静的海面海洋杂波可以近似用高斯模型建模，然而在风浪、潮汐等复杂海况下，SAR后向散射概率分布出现长尾，不再符合高斯分布。对研究海域海洋杂波分布情况进行实验和分析，得到复杂海况的海洋杂波概率密度分布。鉴于海洋杂波的概率密度函数具有较重的长尾，本发明首次将Loglogistic(对数逻辑)分布函数用于拟合复杂海况海洋杂波。Loglogistic的概率密度函数(Probability Density Function，PDF)如下式所示，x>0，α>0，β>0。
Figure BDA0002989844640000071
对数逻辑分布的累积分布函数(Cumulative Distribution Function，CDF)如下式所示：
Figure BDA0002989844640000072
实验表明，Loglogistic概率分布与复杂海况的海洋杂波概率密度拟合优度较高。所以，本发明采用对数逻辑分布对复杂海况的海洋杂波进行建模，基于此进行海洋目标的提取。
采用深度学习模型对28*28像素区域进行10次检测，得到疑似目标和非疑似目标，非疑似目标可视为背景。然后，统计直方图及得到Loglogistic拟合曲线，经Loglogistic统计模型拟合的海洋杂波得到最佳拟合曲线，在像素值233前积分达99.5％，极值像素255的概率值为0.0395，像素总数为28*28*8，图像的分辨率是13.89米*13.89米，极值像素255可达13.89米*13.89米*28*28*8*0.0395＝47797.77m2。所以，仅从概率分布上很难确定47797.77m2的极值像素255是否为目标。
通过多次实验验证，取α＝98.151，β＝3.6501的对数逻辑模型作为海洋杂波拟合模型，此时拟合优度卡方值ks＝0.5908，p＝1，可以认为观测值和预期值是拟合的。对应的标准差σ＝0.0032，绝对误差ME＝0.00198。
实验研究表明，采用tb＝GeneralizedExtremeValuepdffitokss(labb,labtb,rowb)或tb＝mb+σab*t都是基于统计模型，统计模型缺乏空间特征，必须结合聚类等操作，通过空间连通性来确定疑似像素点是否为目标点。最终，对海洋杂波长度介于50.0-600.0米的簇聚类，得到检测因子的检测结果。
四、基于Loglogistic模型的海洋目标精检测
为了改进双参数CFAR逐窗口检测、计算量大的缺陷，采用基于深度学习OceanTDA9模型进行初检，再采用Loglogistic模型进行目标精检测，提取海洋目标特征。计算初检保存在数组x中每个像素的频数绘制直方图，提取直方图中最小像素到最大像素间有效像素及对应的频数做为拟合数据并保存到数组x2、y2中，调用Loglogistic复杂海况的海洋杂波概率密度函数式(3)拟合x2、y2，求出拟合曲线的α、β参数值并计算卡方、绝对值误差等拟合优度。利用下式计算拟合曲线有效像素的积分，当积分值达到阈值时，输出临界像素Xtb。
Figure BDA0002989844640000081
调用目标检测模型PixelsVal＞Xtb，将3*3图像窗口中疑似目标保存到数组row2中。对row2中的28*3行、28*3列图像窗口中疑似目标像素聚簇并保存各个目标簇的参数，包括目标簇的中心坐标等。具体检测流程如图3所示：(1)对行和列均为28*3的图像窗口中疑似目标像素聚簇并保存各个簇的参数；(2)当满足最小宽度和最大长度纵横比时，目标直线拟合提取目标长度、宽度及倾斜角；(3)当满足实际纵横比时进行目标椭圆拟合；(4)保存目标中心坐标、长度、宽度、像素总数及倾斜角；(5)当所有疑似目标检测完毕后结束。
采用Loglogistic模型共提取海洋目标12个，得到编号为7032的提取目标的相关参数可视化结果。提取的钻井平台长13.89*36.0米，宽13.89*15.0米，面积约13.89*13.89*302平方米。该图像所在九宫格的像素均值为116.75，有3个疑似目标子图像，6个海杂波子图像，海杂波像素平均值为116.66，标准差为64.01，基于Loglogistic的海洋目标检测的分割阈值为232，九宫格中像素值大于分割阈值的点数百分比为8.89％。提取目标中心轴拟合直线表达式为y＝0.01x+32.50，九宫格中矩形拟合中心(50.00,33.00)，长36.0像素，宽15.0像素，斜率-0.019，长宽比2.39。待检图像中椭圆拟合参数中心(2206.0,341.0)，长轴中心线斜率-0.019，长半轴17.99像素，短半轴7.51像素，椭圆包含像素总数302.0个。
五、初检和精检结合的海洋目标检测实验
采用OceanTDA9初检和基于Loglogistic的CFAR精检相结合的检测方法——OTD_Loglogistic进行海洋目标检测实验，实验采用并行分布式构架进行，流程如图4所示。图5是双机16核CPU并行计算，0-15核CPU随机分配1-3个任务，最先完成任务的CPU用时3秒，最后完成任务的CPU用时约33秒。主节点的CPU显示第一幅疑似目标的时间是在程序运行3秒或2秒时显示的，其它工作节点第一幅疑似目标分别是在程序运行2、1、2、1或2时显示的。
双机16核OTD_Loglogistic海洋目标参数并行提取各节点CPU任务执行状态如表1所示。
表1双机多核OTD_Loglogistic CPU任务执行状态
Figure BDA0002989844640000082
按照上述OTD_Loglogistic海洋目标检测方法，对研究海域采用并行算法进行实验，得到28*28图像中的非疑似目标、疑似目标和目标。聚类后采用28*28图像中长度大于250米、小于1500米，目标纵横比<＝10去虚警。
本发明通过采用OceanTDA9初检和基于Loglogistic的CFAR精检相结合的检测方法，克服了滑动窗口CFAR需要逐一统计图像中的像素计算耗费长、无法对图像边缘的像素进行处理造成边缘处漏检等缺陷，可以有效提高海洋目标的检测效率和精度。
当然，上述说明并非是对本发明的限制，本发明也并不仅限于上述举例，本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换，也应属于本发明的保护范围。

Claims (4)

1.一种基于OTD_Loglogistic的SAR数据海洋目标检测方法，其特征在于，采用基于深度学习模型的初检和基于Loglogistic的CFAR精检相结合的检测方法，首先，构建OceanTDA9轻量级深度学习模型，采用所构建的轻量级深度学习模型进行海洋目标初检，在此基础上采用基于Loglogistic模型的CFAR方法进行海洋目标精细检测，提取海洋目标特征；
OceanTDA9模型包含4个卷积层、1个卷积组和3个全连接层；每个卷积层的组织形式是Convolution2D-ReLU-Dropout-Maxpooling；卷积组的组织形式是(Convolution2D-ReLU-Dropout)*2-Maxpooling；3个全连接层中，前两个全连接层的组织形式是Dense-ReLU-Dropout，最后一个全连接层的组织形式是Dense；
海洋目标精细检测的原理和过程为：
S201.计算初检保存在数组x中每个像素的频数绘制直方图；
S202.提取直方图中最小像素到最大像素间有效像素及对应的频数作为拟合数据并保存到数组x2、y2中；
S203.调用Loglogistic复杂海况的海洋杂波概率密度函数式
Figure FDA0003569914470000011
拟合x2、y2，其中，x>0，α>0，β>0；
求出拟合曲线的α、β参数；
S204.计算拟合优度，包括卡方和绝对值误差；
S205.计算拟合曲线有效像素的积分
Figure FDA0003569914470000012
当积分值达到阈值时，输出临界像素Xtb；
S206.调用目标检测模型PixelsVal＞Xtb，将疑似目标保存到数组row2中；
S207.对数组row2图像窗口中疑似目标像素聚簇，并保存各个目标簇的参数，包括目标簇中心坐标、长度、宽度、像素总数及倾斜角。
2.根据权利要求1所述的一种基于OTD_Loglogistic的SAR数据海洋目标检测方法，其特征在于，OceanTDA9模型的训练过程为：
S101.对海洋目标样本特征数据x进行归一化处理，采用的归一化公式如下：
Figure FDA0003569914470000013
其中，μ为特征的期望，σ为方差；
S102.配置训练优化算法、损失函数、监控参数；
S103.加载数据集，验证数据的合法性；
S104.将数据集分为多个周期，每个周期分为多个批次；
S105.将每个周期按批次打混索引，进行批次循环；
S106.计算交叉熵，采用反向传播算法和所改进的梯度下降算法WinR-Adagrad，以0.01的学习速率不断地修改变量以最小化交叉熵；
S107.对批次指定的样本数进行训练，学习权重、偏置，计算损失、精度直到一个周期的批次都循环结束，再进入下一个周期进行批循环；
S108.直到循环完所有周期，保存模型及训练损失、精度。
3.根据权利要求2所述的一种基于OTD_Loglogistic的SAR数据海洋目标检测方法，其特征在于，OceanTDA9模型损失函数采用目标分类和模型预测分类之间的交叉熵，其公式为：
Figure FDA0003569914470000021
式中，y-i是第i个样本的输入值，hθ(xi)是第i个样本x的输出值，参数θ0的初值设置为0.1，θ1的初值设置为标准差为0.1的正态分布浮点数。
4.根据权利要求1所述的一种基于OTD_Loglogistic的SAR数据海洋目标检测方法，其特征在于，海洋目标精细检测采用Loglogistic概率密度函数来拟合复杂海况的海杂波分布，突出海杂波的长尾特征，在此基础上采用CFAR方法进行海洋目标特征提取。
CN202110311348.2A 2021-03-24 2021-03-24 一种基于OTD_Loglogistic的SAR数据海洋目标检测方法 Active CN113158806B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110311348.2A CN113158806B (zh) 2021-03-24 2021-03-24 一种基于OTD_Loglogistic的SAR数据海洋目标检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110311348.2A CN113158806B (zh) 2021-03-24 2021-03-24 一种基于OTD_Loglogistic的SAR数据海洋目标检测方法

Publications (2)

Publication Number Publication Date
CN113158806A CN113158806A (zh) 2021-07-23
CN113158806B true CN113158806B (zh) 2022-05-17

Family

ID=76888400

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110311348.2A Active CN113158806B (zh) 2021-03-24 2021-03-24 一种基于OTD_Loglogistic的SAR数据海洋目标检测方法

Country Status (1)

Country Link
CN (1) CN113158806B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114359720B (zh) * 2021-12-23 2024-04-26 湖南国科轩宇信息科技有限公司 一种基于卫星光学图像的海上目标检测方法、系统及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145874A (zh) * 2017-05-13 2017-09-08 复旦大学 复杂背景sar图像中的舰船目标检测与鉴别方法
CN108830242A (zh) * 2018-06-22 2018-11-16 北京航空航天大学 基于卷积神经网络的sar图像海洋目标分类检测方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106156744B (zh) * 2016-07-11 2019-01-29 西安电子科技大学 基于cfar检测与深度学习的sar目标检测方法
CN107563411B (zh) * 2017-08-07 2020-11-24 西安电子科技大学 基于深度学习的在线sar目标检测方法
CN107748361A (zh) * 2017-09-22 2018-03-02 合肥工业大学 基于截断杂波统计的sar图像双参数cfar检测方法
CN109766811B (zh) * 2018-12-31 2023-05-02 复旦大学 一种星载sar图像中海面船只的端到端检测与识别方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107145874A (zh) * 2017-05-13 2017-09-08 复旦大学 复杂背景sar图像中的舰船目标检测与鉴别方法
CN108830242A (zh) * 2018-06-22 2018-11-16 北京航空航天大学 基于卷积神经网络的sar图像海洋目标分类检测方法

Also Published As

Publication number Publication date
CN113158806A (zh) 2021-07-23

Similar Documents

Publication Publication Date Title
CN110598029B (zh) 基于注意力转移机制的细粒度图像分类方法
Ströfer et al. Data-driven, physics-based feature extraction from fluid flow fields
Springenberg et al. Improving deep neural networks with probabilistic maxout units
CN111028327B (zh) 一种三维点云的处理方法、装置及设备
CN109101897A (zh) 水下机器人的目标检测方法、系统及相关设备
CN111914728B (zh) 高光谱遥感影像半监督分类方法、装置及存储介质
CN110765865B (zh) 基于改进的yolo算法的水下目标检测方法
CN109242019B (zh) 一种水面光学小目标快速检测与跟踪方法
CN112115916A (zh) 域适应Faster R-CNN半监督SAR检测方法
CN113344045B (zh) 一种结合hog特征提高sar船只分类精度的方法
He et al. SAR target recognition and unsupervised detection based on convolutional neural network
US20230267703A1 (en) Hierarchical constraint (hc)-based method and system for classifying fine-grained graptolite images
CN113723572B (zh) 船只目标识别方法、计算机系统及程序产品、存储介质
CN116580322A (zh) 一种地面背景下无人机红外小目标检测方法
CN113158806B (zh) 一种基于OTD_Loglogistic的SAR数据海洋目标检测方法
CN116363439A (zh) 一种基于多头自注意力的点云分类方法、装置及设备
CN113436125B (zh) 基于风格迁移的侧扫声呐仿真图像生成方法、装置及设备
CN113096080B (zh) 图像分析方法及系统
Yang et al. An improved algorithm for the detection of fastening targets based on machine vision
CN110738166A (zh) 基于PCNN和PCANet的渔政监控系统红外目标识别方法及存储介质
CN116205918B (zh) 基于图卷积的多模态融合半导体检测方法、装置及介质
CN116758411A (zh) 一种基于遥感图像逐像素处理的舰船小目标检测方法
Wang et al. FPA-DNN: a forward propagation acceleration based deep neural network for ship detection
Yu et al. A lightweight ship detection method in optical remote sensing image under cloud interference
CN115393631A (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