CN112560897A - 基于船长统计模型的雷达舰船分类能力评价方法 - Google Patents

基于船长统计模型的雷达舰船分类能力评价方法 Download PDF

Info

Publication number
CN112560897A
CN112560897A CN202011334711.4A CN202011334711A CN112560897A CN 112560897 A CN112560897 A CN 112560897A CN 202011334711 A CN202011334711 A CN 202011334711A CN 112560897 A CN112560897 A CN 112560897A
Authority
CN
China
Prior art keywords
captain
statistical model
ship
sea
radar
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.)
Granted
Application number
CN202011334711.4A
Other languages
English (en)
Other versions
CN112560897B (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.)
Xidian University
Original Assignee
Xidian University
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 Xidian University filed Critical Xidian University
Priority to CN202011334711.4A priority Critical patent/CN112560897B/zh
Publication of CN112560897A publication Critical patent/CN112560897A/zh
Application granted granted Critical
Publication of CN112560897B publication Critical patent/CN112560897B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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)
  • Data Mining & Analysis (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种基于船长统计模型的雷达舰船分类能力评价方法，包括：获取某一海域的若干船长样本数据；根据船长样本数据得到该海域的船长统计模型；建立不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型；根据船长统计模型和条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线；根据依赖曲线对雷达舰船分类能力进行评价。本发明基于不同海域的船长统计模型，通过给出正确分类概率达到要求时容许的径向尺寸估计误差范围，能够更精确地分析和评价不同海用雷达和舰船径向尺寸估计方法的舰船粗分类能力。

Description

基于船长统计模型的雷达舰船分类能力评价方法
技术领域
本发明属于雷达技术领域，具体涉及一种基于船长统计模型的雷达舰船分类能力评价方法。
背景技术
对海监视雷达是现代海上防御作战系统的主要组成部分，其在海岸和海岛防御中起着重要作用。此外，对海监视雷达还可用于海关缉私监控、海上交通监视等，为切断海上走私通道提供信息引导。高分辨对海监视雷达需要对海面舰船具有“粗”分类的能力，所谓“粗”分类即把海面舰船按照其长度分为“大、中、小”三类：船长小于a米的舰船为小型船，船长介于a 米至b米之间的舰船为中型船，船长大于b米的舰船为大型船，其中a＜b，区分舰船所属类型的分界值a，b通常由工程规定。
对于工作在扫描状态下的高分辨对海监视雷达，常用的分类手段是基于舰船的高分辨距离像(High-resolution range profile，HRRP)对其径向尺寸进行估计，然后从径向尺寸和航向角得到船长的估计值，最后依据船长进行大中小分类。
然而，由于海杂波、目标姿态变化以及高分辨雷达脉冲压缩距离副瓣导致的舰船目标回波扩展等问题，径向尺寸估计均存在一定的误差，该误差的大小可能会导致将目标舰船错误分类，影响监视效果。而对于我国近海舰船粗分类的应用，目前尚没有有效的方法可以实现对不同对海监视雷达和舰船径向尺寸估计方法的舰船粗分类能力的评估，以帮助工程人员判断雷达是否满足工程对正确分类概率的需求。
发明内容
为了解决现有技术中存在的上述问题，本发明提供了一种基于船长统计模型的雷达舰船分类能力评价方法。本发明要解决的技术问题通过以下技术方案实现：
一种基于船长统计模型的雷达舰船分类能力评价方法，包括：
获取某一海域的若干船长样本数据；
根据所述船长样本数据得到该海域的船长统计模型；
建立不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型；
根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线；
根据所述依赖曲线对雷达舰船分类能力进行评价。
在本发明的一个实施例中，根据所述船长样本数据得到该海域的船长统计模型，包括：
根据所述船长样本数据计算船长样本的一、二阶原点矩；
利用所述船长样本的一、二阶原点矩计算船长统计模型的形状参数和尺度参数的估计值；
根据Weibull分布模型、所述形状参数和尺度参数的估计值得到该海域的船长统计模型。
在本发明的一个实施例中，所述船长统计模型的表达式为：
Figure RE-GDA0002940513020000021
其中，f(L；ν,σ)表示形状参数为ν、尺度参数为σ的船长统计模型在船长为L米时的概率密度，exp表示以自然对数为底的指数操作。
在本发明的一个实施例中，所述不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型表示为：
当ε＝0时，θ(L,ε)＝1；
当ε＜0时，θ(L,ε)的表达式为：
Figure RE-GDA0002940513020000031
当ε＞0时，θ(L,ε)的表达式为：
Figure RE-GDA0002940513020000032
其中，θ(L,ε)表示船长为L米的舰船被径向尺寸估计误差为ε米的雷达正确分类的概率，π表示圆周率，arccos表示反余弦操作，cos表示余弦操作，a表示小型船与中型船的分界船长，b表示中型船与大型船的分界船长。
在本发明的一个实施例中，根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线，包括：
根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与正确分类概率之间的函数关系式；
利用软件绘制所述函数关系式，得到所述径向尺寸估计误差和所述正确分类概率之间的依赖曲线。
在本发明的一个实施例中，所述径向尺寸估计误差与正确分类概率之间的函数关系式表示为：
Figure RE-GDA0002940513020000041
其中，PC(ε)表示径向尺寸估计误差为ε的雷达将某一海域舰船正确分类的概率，∫表示积分操作，∞表示无穷大，f(L；ν,σ)表示船长统计模型。
在本发明的一个实施例中，根据所述依赖曲线对雷达舰船分类能力进行评价，包括：
根据工程需求从所述依赖曲线得到径向尺寸估计误差的容许取值范围；
根据雷达的实际径向尺寸估计误差和所述容许取值范围对该雷达舰船分类能力进行评价。
本发明的另一个实施例提供了一种基于船长统计模型的雷达舰船分类能力评价方法，应用于中国近海海域，包括：
获取中国近海海域的若干船长样本数据；其中，所述中国近海海域包括渤海、黄海、东海以及南海；
根据所述船长样本数据分别计算中国近海及渤海、黄海、东海、南海船长统计模型的形状参数和尺度参数的估计值；
根据所述形状参数和尺度参数的估计值得到对应海域的船长统计模型；
建立不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型；
根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线；
根据所述依赖曲线对雷达舰船分类能力进行评价。
在本发明的一个实施例中，所述船长统计模型包括中国近海船长统计模型，所述中国近海船长统计模型的形状参数估计值为νOffshore≈0.95，尺度参数估计值为σOffshore≈62.56；
所述渤海船长统计模型的形状参数估计值为νB≈1.32，尺度参数估计值为σB≈105.69；
所述黄海船长统计模型的形状参数估计值为νH≈0.97，尺度参数估计值为σH≈72.32；
所述东海船长统计模型的形状参数估计值为νD≈0.93，尺度参数估计值为σD≈55.66；
所述南海船长统计模型的形状参数估计值为νN≈0.92，尺度参数估计值为σN≈59.00。
在本发明的一个实施例中，所述对应海域的船长统计模型包括中国近海船长统计模型、渤海船长统计模型、黄海船长统计模型、东海船长统计模型以及南海船长统计模型；其中，
所述中国近海船长统计模型的表达式为：
fOffshore(L)＝0.019L-0.052exp(-0.020L0.948)；
所述渤海船长统计模型的表达式为：
fB(L)＝0.003L0.322exp(-0.002L1.322)；
所述黄海船长统计模型的表达式为：
fH(L)＝0.015L-0.026exp(-0.015L0.974)；
所述东海船长统计模型的表达式为：
fD(L)＝0.022L0.074exp(-0.024L0.926)；
所述南海船长统计模型的表达式为：
fN(L)＝0.022L0.084exp(-0.024L0.916)。
本发明的有益效果：
本发明基于不同海域的船长统计模型，通过给出正确分类概率达到要求时容许的径向尺寸估计误差范围，能够更精确地分析和评价不同海用雷达和舰船径向尺寸估计方法的舰船粗分类能力，从而为高分辨对海监视雷达参数设置和不同径向尺寸估计方法的性能评估提供了参考和依据。
以下将结合附图及实施例对本发明做进一步详细说明。
附图说明
图1是本发明实施例提供的一种基于船长统计模型的雷达舰船分类能力评价方法流程示意图；
图2是本发明实施例提供的利用航向角余弦值取值范围分析将舰船正确分类的条件概率的示意图；
图3是本发明实施例提供的我国近海舰船AIS数据采集示意图；
图4是本发明实施例提供的采用Weibull分布模型拟合得到的中国近海及四大海域船长统计模型；
图5是本发明实施例提供的中国近海及四大海域舰船径向尺寸估计误差与正确分类概率之间的依赖曲线。
具体实施方式
下面结合具体实施例对本发明做进一步详细的描述，但本发明的实施方式不限于此。
实施例一
请参见图1，图1是本发明实施例提供的一种基于船长统计模型的雷达舰船分类能力评价方法流程示意图，包括：
S1：获取某一海域的若干船长样本数据。
具体地，采集某一指定海域内N艘舰船的AIS(Automatic IdentificationSystem，舰船自动识别系统)信息，提取其中的船长数据，记为 L1,L2,…,Li,…,LN(单位：米)，其中Li为第i艘舰船的船长，1≤Li≤400， i＝1,2,…,N。
其中，样本数据越多，即舰船的数量N越大，船长统计模型的拟合效果越好。优选的，本实施例采用N＞30000。
S2：根据船长样本数据得到该海域的船长统计模型。
在本实施例中，采用Weibull分布模型来拟合船长统计模型。
Weibull分布是指数分布的推广，通过尺度参数和形状参数联合控制概率密度曲线的形状和拖尾，能够很好地拟合船长的分布，其概率密度函数表达式如下：
Figure RE-GDA0002940513020000071
其中，f(x；ν,σ)表示形状参数为ν，尺度参数为σ的Weibull分布模型在 x点的概率密度。
Weibull分布模型一、二阶原点矩m1,m2的表达式如下：
Figure RE-GDA0002940513020000081
其中，Γ(·)表示伽马函数。
本实施例采用Weibull分布模型来拟合船长统计模型的过程具体如下：
S21：计算船长样本的一、二阶原点矩
Figure RE-GDA0002940513020000082
为：
Figure RE-GDA0002940513020000083
S22：联立Weibull分布模型一、二阶原点矩表达式，消去σ，用
Figure RE-GDA0002940513020000084
代替m1,m2，得：
Figure RE-GDA0002940513020000085
其中，
Figure RE-GDA0002940513020000086
表示该海域船长统计模型的形状参数估计值。由于上式含有伽马函数，难以获得
Figure RE-GDA0002940513020000087
的解析解，因此采用数值方法解得上式的零点，得到
Figure RE-GDA0002940513020000088
再将其代回Weibull分布模型一阶原点矩表达式，即可得到该海域船长统计模型的尺度参数估计值表达式：
Figure RE-GDA0002940513020000089
S23：根据Weibull分布模型、形状参数和尺度参数的估计值得到船长统计模型。
具体地，将
Figure RE-GDA00029405130200000810
的值代入Weibull分布模型表达式，得到该海域的船长统计模型。
S3：建立不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型。
具体地，设船长满足L≤a的舰船为小型船，满足a＜L≤b的舰船为中型船，满足L＞b的舰船为大型船，其中a＜b。
对于目标舰船，从其HRRP中提取到的径向尺寸R为船长L在雷达视线方向上的投影，舰船的航向角
Figure RE-GDA0002940513020000091
为雷达视线方向与船长方向之间的夹角，在
Figure RE-GDA0002940513020000092
服从均匀分布。三者的关系为
Figure RE-GDA0002940513020000093
设雷达的径向尺寸估计误差为ε，则船长的估计值
Figure RE-GDA0002940513020000094
为：
Figure RE-GDA0002940513020000095
其中，
Figure RE-GDA0002940513020000096
表示径向尺寸估计值。当
Figure RE-GDA0002940513020000097
时，
Figure RE-GDA0002940513020000098
舰船相对于雷达接近切向航行，径向尺寸难以有效估计船长。因此，工程通常将关注的航向角范围限制在区间
Figure RE-GDA0002940513020000099
内。
(3.1)当L≤a，径向尺寸欠估计(ε＜0)时，该小型船总是被正确判定，θ(L,ε)＝1。
(3.2)当L≤a，径向尺寸过估计(ε＞0)时，会出现将该小型船错判为中型船的情况，θ(L,ε)的定义式为：
Figure RE-GDA00029405130200000910
其中，Prob{B|A}表示在事件A发生的条件下，事件B发生的概率。
请参见图2，图2是本发明实施例提供的利用航向角余弦值取值范围分析将舰船正确分类的条件概率的示意图。
根据图2，若
Figure RE-GDA00029405130200000911
落在(1)处，即
Figure RE-GDA00029405130200000912
恒不成立。即在ε＞0，a-ε＜L≤a的情况下，θ(L,ε)＝0。
若
Figure RE-GDA0002940513020000101
落在(2)处，即
Figure RE-GDA0002940513020000102
恒成立。即在ε＞0，
Figure RE-GDA0002940513020000103
的情况下，θ(L,ε)＝1。
若
Figure RE-GDA0002940513020000104
落在(3)处，即
Figure RE-GDA0002940513020000105
当
Figure RE-GDA0002940513020000106
满足
Figure RE-GDA0002940513020000107
时
Figure RE-GDA0002940513020000108
成立，其中|·|表示取绝对值操作。由于
Figure RE-GDA0002940513020000109
在
Figure RE-GDA00029405130200001010
均匀分布，即在ε＞0，
Figure RE-GDA00029405130200001011
的情况下，
Figure RE-GDA00029405130200001012
(3.3)当L＞b，径向尺寸过估计(ε＞0)时，该大型船总是被正确判定，θ(L,ε)＝1。
(3.4)当L＞b，径向尺寸欠估计(ε＜0)时，会出现将该大型船错判为中型船的情况，θ(L,ε)的定义式为：
Figure RE-GDA00029405130200001013
根据图2，若
Figure RE-GDA00029405130200001014
落在(1)处，即
Figure RE-GDA00029405130200001015
恒不成立。即在ε＜0，b＜L＜b-ε的情况下，θ(L,ε)＝0。
若
Figure RE-GDA00029405130200001016
落在(2)处，即
Figure RE-GDA00029405130200001017
恒成立。即在ε＜0，
Figure RE-GDA00029405130200001018
的情况下，θ(L,ε)＝1。
若
Figure RE-GDA00029405130200001019
落在(3)处，即
Figure RE-GDA00029405130200001020
当
Figure RE-GDA00029405130200001021
满足
Figure RE-GDA00029405130200001022
时
Figure RE-GDA00029405130200001023
成立。又由于
Figure RE-GDA00029405130200001024
在
Figure RE-GDA00029405130200001025
均匀分布，即在ε＜0，
Figure RE-GDA0002940513020000111
的情况下，
Figure RE-GDA0002940513020000112
(3.5)当a＜L≤b，径向尺寸过估计(ε＞0)时，会出现将该中型船错判为大型船的情况，θ(L,ε)的定义式为：
Figure RE-GDA0002940513020000113
根据图2，若
Figure RE-GDA0002940513020000114
落在(1)处，即
Figure RE-GDA0002940513020000115
恒不成立。即在ε＞0，b-ε＜L＜b的情况下，θ(L,ε)＝0。
若
Figure RE-GDA0002940513020000116
落在(2)处，即
Figure RE-GDA0002940513020000117
恒成立。即在ε＞0，
Figure RE-GDA0002940513020000118
的情况下，θ(L,ε)＝1。
若
Figure RE-GDA0002940513020000119
落在(3)处，即
Figure RE-GDA00029405130200001110
当
Figure RE-GDA00029405130200001111
满足
Figure RE-GDA00029405130200001112
时
Figure RE-GDA00029405130200001113
成立。即在ε＞0，
Figure RE-GDA00029405130200001114
的情况下，由于
Figure RE-GDA00029405130200001115
在
Figure RE-GDA00029405130200001116
均匀分布，
Figure RE-GDA00029405130200001117
(3.6)当a＜L≤b，径向尺寸欠估计(ε＜0)时，会出现将该中型船错判为小型船的情况，θ(L,ε)的定义式为：
Figure RE-GDA00029405130200001118
根据图2，若
Figure RE-GDA00029405130200001119
落在(1)处，即
Figure RE-GDA00029405130200001120
恒不成立。即在ε＜0，a＜L＜a-ε的情况下，θ(L,ε)＝0。
若
Figure RE-GDA0002940513020000121
落在(2)处，即
Figure RE-GDA0002940513020000122
恒成立。即在ε＜0， a-ε＜L≤b的情况下，θ(L,ε)＝1。
若
Figure RE-GDA0002940513020000123
落在(3)处，即
Figure RE-GDA0002940513020000124
当
Figure RE-GDA0002940513020000125
满足
Figure RE-GDA0002940513020000126
时
Figure RE-GDA0002940513020000127
成立。由于
Figure RE-GDA0002940513020000128
在
Figure RE-GDA0002940513020000129
均匀分布，即在ε＜0，
Figure RE-GDA00029405130200001210
的情况下，
Figure RE-GDA00029405130200001211
综合上述六种情况，不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型可以表示为：
当ε＝0时，θ(L,ε)＝1；
当ε＜0时，θ(L,ε)的表达式为：
Figure RE-GDA00029405130200001212
当ε＞0时，θ(L,ε)的表达式为：
Figure RE-GDA0002940513020000131
S4：根据船长统计模型和条件概率计算模型得到径向尺寸估计误差与该海域所有舰船被正确分类的概率之间的依赖曲线，具体包括：
S41：根据船长统计模型和条件概率计算模型得到径向尺寸估计误差与正确分类概率之间的函数关系式，其表达式为：
Figure RE-GDA0002940513020000132
S42：利用软件绘制所述函数关系式，得到径向尺寸估计误差和正确分类概率之间的依赖曲线。
具体地，本实施例采用软件MATLAB的plot函数绘制出径向尺寸估计误差与舰船正确分类概率的依赖曲线。
S5：根据依赖曲线对雷达舰船分类能力进行评价，包括：
S51：根据工程需求从依赖曲线得到径向尺寸估计误差的容许取值范围。
S52：根据雷达的实际径向尺寸估计误差和容许取值范围对该雷达舰船分类能力进行评价。
若雷达的径向尺寸估计误差落在该容许范围之内，则正确分类概率能够满足工程需求；反之则不能满足工程需求。
本发明基于不同海域的船长统计模型，通过给出正确分类概率达到要求时容许的径向尺寸估计误差范围，能够更精确地分析和评价不同海用雷达和舰船径向尺寸估计方法的舰船粗分类能力，从而为高分辨对海监视雷达参数设置和不同径向尺寸估计方法的性能评估提供了参考和依据。
实施例二
对于我国近海舰船粗分类的应用，感兴趣的问题是径向尺寸估计误差需控制在怎样的范围内，能够满足工程对正确分类概率的需求。该问题的解决与船长的先验概率分布密切相关，然而在公开数据中，没有中国近海船长分布的数据和统计模型。
基于此，本实施例将上述实施例一提供的基于船长统计模型的雷达舰船分类能力评价方法应用在中国近海海域，得到中国近海船长统计模型，以用于更精确地分析和评价不同海用雷达和舰船径向尺寸估计方法的舰船粗分类能力。具体过程如下：
步骤1：获取中国近海海域的若干船长样本数据。
具体地，请参见图3，图3是本发明实施例提供的我国近海舰船AIS数据采集示意图。我国近海自北向南分为渤海、黄海、东海、南海四大海域。渤海-黄海分界线为辽东半岛南端老铁山角至山东半岛北岸蓬莱角的连线，黄海-东海分界线为长江口北岸启东角至韩国济州岛西南角的连线，东海-南海分界线为广东南澳岛至中国台湾南端鹅銮鼻之间的连线，如图3中(a)所示。数据采集基于船讯网提供的舰船AIS信息，包括船长、船宽、类型等静态信息和吃水、经纬度等动态信息，如图3中(c)所示。采集范围为渤海整体，南海十段线以内，黄海、东海距离领海基线200海里以内的海域(即专属经济区)。由于海上舰船数量繁多，分布密集，为便于采集，将海域划分为若干网格，以网格为单位进行采集和统计，如图3中(b)所示(以东海舟山群岛的一小片海域为例)。
此次采集时间自2020年三月中旬至三月底，为期约两周，由于采集期间舰船的位置会动态更新，为防止重复统计，以舰船首次出现所在的海域为准。
本实施例采用的大中小分类标准为船长小于69m为小型船，船长介于 69m至180m之间为中型船，船长大于180m为大型船，即a＝69m，b＝180m。本实施例共采集了30584艘舰船的AIS信息，包括东海14308艘(47％)、南海 9635艘(31％)、黄海4306艘(14％)、渤海2335艘(8％)；其中小型船22291艘 (73％)、中型船5972艘(19％)、大型船2321艘(8％)。
步骤2：根据船长样本数据分别计算中国近海及渤海、黄海、东海、南海海域船长统计模型的形状参数和尺度参数的估计值。
具体地，根据步骤1采集的样本数据得到中国近海船长统计模型的形状参数为νOffshore≈0.95，尺度参数为σOffshore≈62.56；渤海船长统计模型的形状参数为νB≈1.32，尺度参数为σB≈105.69；黄海船长统计模型的形状参数为νH≈0.97，尺度参数为σH≈72.32；东海船长统计模型的形状参数为νD≈0.93，尺度参数为σD≈55.66；南海船长统计模型的形状参数为νN≈0.92，尺度参数为σN≈59.00。
步骤3：根据所述形状参数和尺度参数的估计值得到对应海域的船长统计模型。
具体地，请参见图4，图4是本发明实施例提供的采用Weibull分布模型拟合得到的中国近海及四大海域船长统计模型，其中，
中国近海船长统计模型表达式为：
fOffshore(L)＝0.019L-0.052exp(-0.020L0.948)；
渤海船长统计模型表达式为：
fB(L)＝0.003L0.322exp(-0.002L1.322)；
黄海船长统计模型表达式为：
fH(L)＝0.015L-0.026exp(-0.015L0.974)；
东海船长统计模型表达式为：
fD(L)＝0.022L0.074exp(-0.024L0.926)；
南海船长统计模型表达式为：
fN(L)＝0.022L0.084exp(-0.024L0.916)。
步骤4：建立不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型。
步骤5：根据船长统计模型和条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线。
具体地，根据中国近海及四大海域的船长统计模型，绘制出了各海域正确分类概率与径向尺寸估计误差的依赖曲线，如图5所示。由依赖曲线可知有0≤PC≤1，当ε＝0时PC＝1，PC随着|ε|的增大而减小；曲线关于ε＝0左右不对称，且过估计对PC的影响略大于欠估计的影响，这是由于海上小型船所占比重或先验概率较大，而小型船的正确分类概率仅受到过估计的影响的缘故。
步骤6：根据所述依赖曲线对雷达舰船分类能力进行评价。
在本实施例中，设工程所需求的正确分类概率为90％，当中国近海舰船的径向尺寸估计误差落在区间(-12.6m,9.4m)内时，正确分类概率能够满足工程需求。各海域的情况因船长先验概率分布参数的差异而有所不同：渤海正确分类概率达到90％的要求最高，径向尺寸估计误差需落在区间(-8.1m, 7.3m)内；其次是黄海，径向尺寸估计误差需落在区间(-11.5m,8.9m)内；南海径向尺寸估计误差需落在区间(-13.2m,9.6m)内；东海的要求最低，径向尺寸估计误差需落在区间(-13.8m,9.7m)内。
最后根据雷达的实际径向尺寸估计误差进行评估，例如，雷达的实际径向尺寸估计误差为2m时，按照本发明的舰船分类能力评价方法，其正确分类概率达到93％，能够满足正确分类概率为90％的工程需求。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (10)

1.一种基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，包括：
获取某一海域的若干船长样本数据；
根据所述船长样本数据得到该海域的船长统计模型；
建立不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型；
根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线；
根据所述依赖曲线对雷达舰船分类能力进行评价。
2.根据权利要求1所述的基于船长统计模型的雷达舰船分类能力评价方法，根据所述船长样本数据得到该海域的船长统计模型，包括：
根据所述船长样本数据计算船长样本的一、二阶原点矩；
利用所述船长样本的一、二阶原点矩计算船长统计模型的形状参数和尺度参数的估计值；
根据Weibull分布模型、所述形状参数和尺度参数的估计值得到船长统计模型。
3.根据权利要求2所述的基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，所述船长统计模型的表达式为：
Figure FDA0002796843670000011
其中，f(L；ν,σ)表示形状参数为ν、尺度参数为σ的船长统计模型在船长为L米时的概率密度，exp表示以自然对数为底的指数操作。
4.根据权利要求1所述的基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，所述不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型表示为：
当ε＝0时，θ(L,ε)＝1；
当ε＜0时，θ(L,ε)的表达式为：
Figure FDA0002796843670000021
当ε＞0时，θ(L,ε)的表达式为：
Figure FDA0002796843670000022
其中，θ(L,ε)表示船长为L米的舰船被径向尺寸估计误差为ε米的雷达正确分类的概率，π表示圆周率，arccos表示反余弦操作，cos表示余弦操作，a表示小型船与中型船的分界船长，b表示中型船与大型船的分界船长。
5.根据权利要求1所述的基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线，包括：
根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与正确分类概率之间的函数关系式；
利用软件绘制所述函数关系式，得到所述径向尺寸估计误差和所述正确分类概率之间的依赖曲线。
6.根据权利要求5所述的基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，所述径向尺寸估计误差与正确分类概率之间的函数关系式表示为：
Figure FDA0002796843670000031
其中，PC(ε)表示径向尺寸估计误差为ε的雷达将某一海域内任一舰船正确分类的概率，∫表示积分操作，∞表示无穷大。
7.根据权利要求1所述的基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，根据所述依赖曲线对雷达舰船分类能力进行评价，包括：
根据工程需求从所述依赖曲线得到径向尺寸估计误差的容许取值范围；
根据雷达的实际径向尺寸估计误差和所述容许取值范围对该雷达舰船分类能力进行评价。
8.一种基于船长统计模型的雷达舰船分类能力评价方法，应用于中国近海海域，其特征在于，包括：
获取中国近海海域的若干船长样本数据；其中，所述中国近海海域包括渤海、黄海、东海以及南海；
根据所述船长样本数据分别计算中国近海及渤海、黄海、东海、南海船长统计模型的形状参数和尺度参数的估计值；
根据所述形状参数和尺度参数的估计值得到对应海域的船长统计模型；
建立不同径向尺寸估计误差下不同长度的舰船被正确分类的条件概率计算模型；
根据所述船长统计模型和所述条件概率计算模型得到径向尺寸估计误差与该海域任一目标舰船被正确分类的概率之间的依赖曲线；
根据所述依赖曲线对雷达舰船分类能力进行评价。
9.根据权利要求8所述的基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，所述中国近海船长统计模型的形状参数估计值为νOffshore≈0.95，尺度参数估计值为σOffshore≈62.56；
所述渤海船长统计模型的形状参数估计值为νB≈1.32，尺度参数估计值为σB≈105.69；
所述黄海船长统计模型的形状参数估计值为νH≈0.97，尺度参数估计值为σH≈72.32；
所述东海船长统计模型的形状参数估计值为νD≈0.93，尺度参数估计值为σD≈55.66；
所述南海船长统计模型的形状参数估计值为νN≈0.92，尺度参数估计值为σN≈59.00。
10.根据权利要求8所述的基于船长统计模型的雷达舰船分类能力评价方法，其特征在于，所述对应海域的船长统计模型包括中国近海船长统计模型、渤海船长统计模型、黄海船长统计模型、东海船长统计模型以及南海船长统计模型；其中，
所述中国近海船长统计模型的表达式为：
fOffshore(L)＝0.019L-0.052exp(-0.020L0.948)；
所述渤海船长统计模型的表达式为：
fB(L)＝0.003L0.322exp(-0.002L1.322)；
所述黄海船长统计模型的表达式为：
fH(L)＝0.015L-0.026exp(-0.015L0.974)；
所述东海船长统计模型的表达式为：
fD(L)＝0.022L0.074exp(-0.024L0.926)；
所述南海船长统计模型的表达式为：
fN(L)＝0.022L0.084exp(-0.024L0.916)。
CN202011334711.4A 2020-11-24 2020-11-24 基于船长统计模型的雷达舰船分类能力评价方法 Active CN112560897B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011334711.4A CN112560897B (zh) 2020-11-24 2020-11-24 基于船长统计模型的雷达舰船分类能力评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011334711.4A CN112560897B (zh) 2020-11-24 2020-11-24 基于船长统计模型的雷达舰船分类能力评价方法

Publications (2)

Publication Number Publication Date
CN112560897A true CN112560897A (zh) 2021-03-26
CN112560897B CN112560897B (zh) 2024-02-09

Family

ID=75043483

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011334711.4A Active CN112560897B (zh) 2020-11-24 2020-11-24 基于船长统计模型的雷达舰船分类能力评价方法

Country Status (1)

Country Link
CN (1) CN112560897B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016101279A1 (zh) * 2014-12-26 2016-06-30 中国海洋大学 一种合成孔径雷达图像舰船目标快速检测方法
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
CN107122764A (zh) * 2017-05-22 2017-09-01 中国人民解放军国防科学技术大学 基于KpN模型的海上舰船目标检测方法
CN107145874A (zh) * 2017-05-13 2017-09-08 复旦大学 复杂背景sar图像中的舰船目标检测与鉴别方法
CN108490415A (zh) * 2018-05-14 2018-09-04 浙江大学 一种智能sar雷达海上舰船目标识别系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016101279A1 (zh) * 2014-12-26 2016-06-30 中国海洋大学 一种合成孔径雷达图像舰船目标快速检测方法
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
CN107145874A (zh) * 2017-05-13 2017-09-08 复旦大学 复杂背景sar图像中的舰船目标检测与鉴别方法
CN107122764A (zh) * 2017-05-22 2017-09-01 中国人民解放军国防科学技术大学 基于KpN模型的海上舰船目标检测方法
CN108490415A (zh) * 2018-05-14 2018-09-04 浙江大学 一种智能sar雷达海上舰船目标识别系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
侯颖妮;杨予昊;李士国;江涛;: "基于ISAR像的舰船目标识别技术研究", 现代雷达, no. 03 *
魏存伟;刘先康;孙菲;徐冰超;: "一种基于雷达HRRP的舰船目标识别方法", 电子科技, no. 11 *
魏存伟;段发阶;刘先康;: "基于宽带雷达HRRP舰船目标长度估计算法", 系统工程与电子技术, no. 09 *

Also Published As

Publication number Publication date
CN112560897B (zh) 2024-02-09

Similar Documents

Publication Publication Date Title
CN107145874B (zh) 复杂背景sar图像中的舰船目标检测与鉴别方法
CN103971127B (zh) 一种前视雷达成像海面目标关键点检测识别方法
US20170016987A1 (en) Processing synthetic aperture radar images for ship detection
CN109711353B (zh) 一种基于机器视觉的船舶吃水线区域识别方法
US9031285B2 (en) Detection of floating objects in maritime video using a mobile camera
CN110084129B (zh) 一种基于机器视觉的河流漂浮物实时检测方法
CN111694913B (zh) 一种基于卷积自编码器的船舶ais轨迹聚类方法和装置
CN113112540B (zh) 一种在vts系统中使用ais标定cctv摄像头进行船舶图像目标定位的方法
CN110910448B (zh) 一种海上溢油漂移轨迹溯源方法
CN110765912B (zh) 一种基于统计约束和Mask R-CNN的SAR图片船舶目标检测方法
CN112699315A (zh) 基于ais数据的船舶异常行为探测方法
CN114972918B (zh) 基于集成学习与ais数据的遥感图像舰船目标识别方法
CN109284663A (zh) 一种基于正态和均匀混合分布模型的海面障碍物检测方法
CN115272876A (zh) 一种基于深度学习的遥感图像船舶目标检测方法
CN111553183A (zh) 一种船只检测模型训练方法、船只检测方法及装置
Kamirul et al. Automatic ship recognition chain on satellite multispectral imagery
CN112347827A (zh) 一种船舶水尺自动检测方法及系统
Heyn et al. A system for automated vision-based sea-ice concentration detection and floe-size distribution indication from an icebreaker
CN112560897A (zh) 基于船长统计模型的雷达舰船分类能力评价方法
CN112966656A (zh) 数据处理方法及装置
Zhang et al. FRS-Net: An efficient ship detection network for thin-cloud and FOG-covered high-resolution optical satellite imagery
CN111563428A (zh) 一种机场停机位入侵检测方法及其系统
CN113011376B (zh) 海上船舶遥感分类方法、装置、计算机设备及存储介质
CN113641954B (zh) 一种复杂海洋环境中的三维声场快速预报方法及系统
CN108627825A (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