CN114488107A - 一种海杂波时空分布及影响分级产品制作的方法及装置 - Google Patents

一种海杂波时空分布及影响分级产品制作的方法及装置 Download PDF

Info

Publication number
CN114488107A
CN114488107A CN202210382406.5A CN202210382406A CN114488107A CN 114488107 A CN114488107 A CN 114488107A CN 202210382406 A CN202210382406 A CN 202210382406A CN 114488107 A CN114488107 A CN 114488107A
Authority
CN
China
Prior art keywords
sea clutter
sea
time distribution
calculation model
space
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
CN202210382406.5A
Other languages
English (en)
Other versions
CN114488107B (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.)
Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Aerospace Hongtu Information Technology Co Ltd
Original Assignee
Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Aerospace Hongtu Information Technology Co Ltd
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 Southern Marine Science and Engineering Guangdong Laboratory Guangzhou, Aerospace Hongtu Information Technology Co Ltd filed Critical Southern Marine Science and Engineering Guangdong Laboratory Guangzhou
Priority to CN202210382406.5A priority Critical patent/CN114488107B/zh
Publication of CN114488107A publication Critical patent/CN114488107A/zh
Application granted granted Critical
Publication of CN114488107B publication Critical patent/CN114488107B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06—Systems determining position data of a target
    • 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/93—Radar or analogous systems specially adapted for specific applications for anti-collision purposes
    • G01S13/937—Radar or analogous systems specially adapted for specific applications for anti-collision purposes of marine craft
    • 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/95—Radar or analogous systems specially adapted for specific applications for meteorological use
    • G01S13/953—Radar or analogous systems specially adapted for specific applications for meteorological use mounted on aircraft
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40—Means for monitoring or calibrating
    • G01S7/4052—Means for monitoring or calibrating by simulation of echoes
    • 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
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G01S7/414—Discriminating targets with respect to background clutter
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • 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)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Electromagnetism (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Ocean & Marine Engineering (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种海杂波时空分布及影响分级产品制作的方法及装置，涉及海洋工程的技术领域，包括：获取目标海域的目标数据和目标海域的实测海杂波数据；基于实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；基于目标数据、校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出海杂波时空分布产品；基于的海杂波时空分布产品，确定出海杂波时空分布产品的分级阈值；基于分级阈值，对海杂波时空分布产品进行分级，得到的海杂波影响分级产品，解决了难以对海杂波进行大尺度分布趋势分析和影响效能评估的技术问题。

Description

一种海杂波时空分布及影响分级产品制作的方法及装置
技术领域
本发明涉及海洋工程的技术领域，尤其是涉及一种海杂波时空分布及影响分级产品制作的方法及装置。
背景技术
对于海面目标探测雷达而言，海表面对发射信号的后向散射常常严重地限制雷达对舰船、飞机、导弹、导航浮标和其他和海表面同在一个雷达分辨单元内的目标的探测能力。这个干扰信号通常称为海杂波。海杂波除了与海洋环境参数（如海面风场、海浪、潮汐等）相关，还与雷达参数（频率、极化、观测方法、擦地角）密切相关。为了更好了解这种海洋干扰信息，就需要建立海杂波信号与雷达参数和海洋环境参数的关系。
目前针对海杂波特性的研究，多采用试验方法，通过岸基雷达、机载航飞试验的方式获得海面回波测量数据。该方法只能获取某区域或很小范围内的海杂波结果，缺少对大尺度海域海杂波分布趋势的分析。此外，目前海杂波后向散射系数计算模型，可以分为两类，一类为半经验的海杂波模型GTI，TSC，HYB，NRL；一类为解析近似模型，如KA，SPM，复合BRAG模型。各个模型的适用范围有限，单一模型无法更好用于计算多种的应用场景。
针对上述问题，还未提出有效的解决方案。
发明内容
有鉴于此，本发明的目的在于提供一种海杂波时空分布及影响分级产品制作的方法及装置，以缓解了现有技术难以对海杂波进行大尺度分布趋势分析和影响效能评估的技术问题。
第一方面，本发明实施例提供了一种海杂波时空分布及影响分级产品制作的方法，包括：获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品。
进一步地，所述海况数据包括：海面风速和海面风向；所述雷达观测参数包括观测频率、观测方位角、擦地角和极化方式。
进一步地，基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型，包括：基于所述实测海杂波数据，确定出海杂波强度的计算参数的观测值和所述观测值对应的海况数据，其中，所述海杂波强度的计算参数包括：观测频率、擦地角、观测相对风向和极化方式；基于所述观测值对应的海况数据和所述海杂波半经验计算模型，计算出所述海杂波强度的计算参数的理论值；基于所述观测值和所述理论值，对所述海杂波半经验计算模型进行偏差校正，得到所述校正后的海杂波半经验计算模型。
进一步地，基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品，包括：确定出所述目标数据中的第一擦地角和第二擦地角，其中，所述第一擦地角大于或等于预设阈值，所述第二擦地角小于所述预设阈值；基于所述第一擦地角对应的目标数据和所述复合布拉格后向散射计算模型，确定出所述第一擦地角的海杂波时空分布产品；基于所述第二擦地角对应的目标数据和所述校正后的海杂波半经验计算模型，确定出所述第二擦地角的海杂波时空分布产品。
进一步地，基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值，包括：计算出所述海杂波时空分布产品的累计概率密度分布；基于所述累计概率密度分布，确定出累计概率分割点；基于所述累计概率分割点，确定出所述海杂波时空分布产品的分级阈值。
进一步地，基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品，包括：确定出所述海杂波时空分布产品的海杂波影响强度值；基于所述海杂波影响强度值和所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品。
第二方面，本发明实施例还提供了一种海杂波时空分布及影响分级产品制作的系统，包括：获取单元，校正单元，海杂波时空分布产品生成单元，确定单元和海杂波影响分级产品生成单元，其中，所述获取单元，用于获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；所述校正单元，用于基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；所述海杂波时空分布产品生成单元，用于基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；所述确定单元，用于基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；所述海杂波影响分级产品生成单元，用于基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品。
进一步地，所述海况数据包括：海面风速和海面风向；所述雷达观测参数包括观测频率、观测方位角、擦地角和极化方式。第三方面，本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述第一方面中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
第四方面，本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序。
在本发明实施例中，通过获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品，达到了生成海杂波时空分布产品和海杂波影响分级产品的目的，进而解决了现有技术难以对海杂波进行大尺度分布趋势分析和影响效能评估的技术问题，从而实现了对海杂波进行大尺度分布趋势分析和影响效能评估的技术效果。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种海杂波时空分布及影响分级产品制作的方法的流程图；
图2为本发明实施例提供的一种海杂波时空分布及影响分级产品制作的系统的示意图；
图3为本发明实施例提供的一种电子设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一：
根据本发明实施例，提供了一种海杂波时空分布及影响分级产品制作的方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的一种海杂波时空分布及影响分级产品制作的方法的流程图，如图1所示，该方法包括如下步骤：
步骤S102，获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；
需要说明的是，上述的海况数据包括：海面风速和海面风向；上述的雷达观测参数包括观测频率、观测方位角、擦地角和极化方式。
步骤S104，基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；
步骤S106，基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；
步骤S108，基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；
步骤S110，基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品。
在本发明实施例中，通过获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品，达到了生成海杂波时空分布产品和海杂波影响分级产品的目的，进而解决了现有技术难以对海杂波进行大尺度分布趋势分析和影响效能评估的技术问题，从而实现了对海杂波进行大尺度分布趋势分析和影响效能评估的技术效果。
在本发明实施例中，步骤S104包括如下步骤：
步骤S11，基于所述实测海杂波数据，确定出海杂波强度的计算参数的观测值和所述观测值对应的海况数据，其中，所述海杂波强度的计算参数包括：观测频率、擦地角、观测相对风向和极化方式；
步骤S12， 基于所述观测值对应的海况数据和所述海杂波半经验计算模型，计算出所述海杂波强度的计算参数的理论值；
步骤S13，基于所述观测值和所述理论值，对所述海杂波半经验计算模型进行偏差校正，得到所述校正后的海杂波半经验计算模型。
需要说明的是，观测相对风向等于雷达观测方位角与海面风向的夹角，一般分为三类，0°代表顺风观测，90°代表测风观测，180°代表逆风观测。擦地角，一般可分为1°,10°,30°,60°；雷达观测参数分为L，S，C，Ku波段；极化方式一般分为HH和VV两类。
在本发明实施例中，优选的，实测海杂波数据选择Nathansan海杂波数据表，该数据表包含了七个波段、七种擦地角、七种海况、两种极化方式所对应的后向散射系数测量值，是目前为止海杂波散射率比较完整的公开数据库。该数据表不区分具体的观测风向，记录了在顺风、逆风和测风观测下海杂波测量的平均值。半经验海杂波模型选择TechnologyService Corporation模型（TSC）。TSC模型是在总结实测平均后向散射系数数据的基础上提出的，主要包含三个影响因子，即小擦地角、风速、风向因子；同时考虑了异常传播情况，如蒸发波导的影响，所该模型得到的散射系数变化较平缓。在未知传播条件时，TSC模型能更好地描述海面状况，预测后向散射系数。
所述观测值和所述理论值之间的偏差校正量的计算公式如下：
Figure P_220413090802873_873047001
；
式中，
Figure P_220413090802935_935553001
为观测频率，
Figure P_220413090802951_951178002
为擦地角，
Figure P_220413090802966_966808003
为极化方式，
Figure P_220413090802998_998061004
为顺风、逆风、测风对应的相对观测，取值为0，90°和180°，
Figure P_220413090803013_013707005
为计算模型，
Figure P_220413090803029_029310006
为Nathansan数据表，
Figure P_220413090803060_060557007
计算出来的在该组频率、入射角及极化下计算出的偏差校正量。
TSM半经验海杂波模型的计算公式如下:
Figure P_220413090803076_076182001
。
式中，
Figure P_220413090803107_107469001
为HH极化方式得到的后向散射系数，
Figure P_220413090803138_138682002
为VV极化方式得到的后向散射系数，
Figure P_220413090803154_154315003
为擦地角，单位为弧度，
Figure P_220413090803169_169937004
为小擦地角计算因子，
Figure P_220413090803201_201210005
为风向因子，
Figure P_220413090803216_216812006
为风速因子，
Figure P_220413090803232_232420007
为波长，单位为英尺，
Figure P_220413090803253_253394008
海面高度标准差，单位为英尺，
Figure P_220413090803285_285150009
为观测频率。
通过观测值与理论值相减得到偏差校正量，并将偏差校正量作为误差项补偿在海杂波半经验计算模型中，得到校正后的海杂波半经验计算模型。
在本发明实施例中，步骤S106包括如下步骤：
步骤S21，确定出所述目标数据中的第一擦地角和第二擦地角，其中，所述第一擦地角大于或等于预设阈值，所述第二擦地角小于所述预设阈值；
步骤S22，基于所述第一擦地角对应的目标数据和所述复合布拉格后向散射计算模型，确定出所述第一擦地角的海杂波时空分布产品；
步骤S23，基于所述第二擦地角对应的目标数据和所述校正后的海杂波半经验计算模型，确定出所述第二擦地角的海杂波时空分布产品。
在本发明实施例中，首先对目标数据中的擦地角的大小进行判断，进而选取对应的计算模型。如果擦地角小于10°（预设阈值），则该擦地角为第一擦地角，选用校正后的海杂波半经验计算模型，并利用校正后的海杂波半经验计算模型和第一擦地角对应的目标数据，确定出第一擦地角的海杂波时空分布产品。
如果擦地角大于或等于10°（预设阈值），则该擦地角为第二擦地角，选用复合布拉格后向散射计算模型，并利用复合布拉格后向散射计算模型和第二擦地角对应的目标数据，确定出第二擦地角的海杂波时空分布产品。
需要说明的是，校正后的海杂波半经验计算模型
Figure P_220413090803300_300802001
计算公式如下描述：
Figure P_220413090803316_316413001
；
式中，
Figure P_220413090803347_347675001
为校正后的海杂波半经验计算模型；
Figure P_220413090803363_363299002
为海杂波半经验计算模型，
Figure P_220413090803394_394565003
为补偿偏差。
复合布拉格散射模型通过将表面粗糙度视为大小两种粗糙度的叠加，将传统布拉格散射计算的散射系数在大粗糙度的斜率分布上作集平均。
复合布拉格散射模型的计算公式如下描述：
Figure P_220413090803410_410164001
；
式中，
Figure P_220413090803442_442362001
雷达观测波数，
Figure P_220413090803474_474184002
为雷达波入射角，
Figure P_220413090803489_489744003
为海面波浪的波数谱，
Figure P_220413090803505_505380004
相对观测方向，
Figure P_220413090803521_521015005
为归一化的后向散射系数，下标
Figure P_220413090803552_552263006
和
Figure P_220413090803567_567900007
表示极化（H或V），
Figure P_220413090803583_583532008
为本地入射角，
Figure P_220413090803614_614772009
和
Figure P_220413090803630_630385010
分别为x和y方向海面坡面的倾角。
在本发明实施例中，步骤S108包括如下步骤：
步骤S31，计算出所述海杂波时空分布产品的累计概率密度分布；
步骤S32，基于所述累计概率密度分布，确定出累计概率分割点；
步骤S33，基于所述累计概率分割点，确定出所述海杂波时空分布产品的分级阈值。
优选的，上述的预设时间一般为12个月。
在本发明实施例中，首先，计算的海杂波时空分布产品的累计概率密度分布。
然后，确定出累计概率分割点，在本发明实施例中，累计概率分割点为两个，分别为0.5和0.9。
最后，根据累计概率分割点，确定出海杂波时空分布产品的分级阈值。
在累计概率分割点为0.5和0.9时分割阈值A和B分别为-33.1079和-30.2807。
在本发明实施例中，步骤S110包括如下步骤：
步骤S41，确定出所述海杂波时空分布产品的海杂波影响强度值；
步骤S42，基于所述海杂波影响强度值和所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品。
在本发明实施例中利用的海杂波时空分布产品，在的海杂波时空分布产品中，定义海杂波影响强度值小于阈值A的海杂波时空分布产品对应的影响区域，确定为弱影响区域；将海杂波影响强度值大于阈值B的海杂波时空分布产品对应的影响区域，确定为强影响区域；其余海杂波时空分布产品对应的影响区域，确定为中等影响区域。
本发明实施例，突破传统海杂波分析多基于就位探测数据的局限，同时提出了多入射角的海杂波组合计算模型和基于幅度分布特性的影响分级模型，有利于进一步用于海杂波大尺度分布趋势分析及海杂波影响效能评估。
实施例二：
本发明实施例还提供了一种海杂波时空分布及影响分级产品制作的系统，该海杂波时空分布及影响分级产品制作的系统用于执行本发明实施例上述内容所提供的海杂波时空分布及影响分级产品制作的方法，以下是本发明实施例提供的海杂波时空分布及影响分级产品制作的系统的具体介绍。
如图2所示，图2为上述海杂波时空分布及影响分级产品制作的系统的示意图，该海杂波时空分布及影响分级产品制作的系统包括：获取单元10，校正单元20，海杂波时空分布产品生成单元30，确定单元40和海杂波影响分级产品生成单元50。
所述获取单元10，用于获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；
所述校正单元20，用于基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；
所述海杂波时空分布产品生成单元30，用于基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；
所述确定单元40，用于基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；
所述海杂波影响分级产品生成单元50，用于基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品。
在本发明实施例中，通过获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品，达到了生成海杂波时空分布产品和海杂波影响分级产品的目的，进而解决了现有技术难以对海杂波进行大尺度分布趋势分析和影响效能评估的技术问题，从而实现了对海杂波进行大尺度分布趋势分析和影响效能评估的技术效果。
优选的，所述海况数据包括：海面风速和海面风向；所述雷达观测参数包括观测频率、观测方位角、擦地角和极化方式。
实施例三：
本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述实施例一中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
参见图3，本发明实施例还提供一种电子设备100，包括：处理器60，存储器61，总线62和通信接口63，所述处理器60、通信接口63和存储器61通过总线62连接；处理器60用于执行存储器61中存储的可执行模块，例如计算机程序。
其中，存储器61可能包含高速随机存取存储器（RAM，Random Access Memory），也可能还包括非不稳定的存储器（non-volatile memory），例如至少一个磁盘存储器。通过至少一个通信接口63（可以是有线或者无线）实现该系统网元与至少一个其他网元之间的通信连接，可以使用互联网，广域网，本地网，城域网等。
总线62可以是ISA总线、PCI总线或EISA总线等。所述总线可以分为地址总线、数据总线、控制总线等。为便于表示，图3中仅用一个双向箭头表示，但并不表示仅有一根总线或一种类型的总线。
其中，存储器61用于存储程序，所述处理器60在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的装置所执行的方法可以应用于处理器60中，或者由处理器60实现。
处理器60可能是一种集成电路芯片，具有信号的处理能力。在实现过程中，上述方法的各步骤可以通过处理器60中的硬件的集成逻辑电路或者软件形式的指令完成。上述的处理器60可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(Digital SignalProcessing，简称DSP)、专用集成电路(Application Specific Integrated Circuit，简称ASIC)、现成可编程门阵列(Field-Programmable Gate Array，简称FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本发明实施例所公开的方法的步骤可以直接体现为硬件译码处理器执行完成，或者用译码处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器，闪存、只读存储器，可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器61，处理器60读取存储器61中的信息，结合其硬件完成上述方法的步骤。
实施例四：
本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，计算机程序被处理器运行时执行上述实施例一中所述方法的步骤。
另外，在本发明实施例的描述中，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。
在本申请所提供的几个实施例中，应该理解到，所揭露的系统、装置和方法，可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，又例如，多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些通信接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

Claims (10)

1.一种海杂波时空分布及影响分级产品制作的方法，其特征在于，包括：
获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；
基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；
基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；
基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；
基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到海杂波影响分级产品。
2.根据权利要求1所述的方法，其特征在于，
所述海况数据包括：海面风速和海面风向；
所述雷达观测参数包括观测频率、观测方位角、擦地角和极化方式。
3.根据权利要求2所述的方法，其特征在于，基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型，包括：
基于所述实测海杂波数据，确定出海杂波强度的计算参数的观测值和所述观测值对应的海况数据，其中，所述海杂波强度的计算参数包括：观测频率、擦地角、观测相对风向和极化方式；
基于所述观测值对应的海况数据和所述海杂波半经验计算模型，计算出所述海杂波强度的计算参数的理论值；
基于所述观测值和所述理论值，对所述海杂波半经验计算模型进行偏差校正，得到所述校正后的海杂波半经验计算模型。
4.根据权利要求1所述的方法，其特征在于，基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品，包括：
确定出所述目标数据中的第一擦地角和第二擦地角，其中，所述第一擦地角大于或等于预设阈值，所述第二擦地角小于所述预设阈值；
基于所述第一擦地角对应的目标数据和所述复合布拉格后向散射计算模型，确定出所述第一擦地角的海杂波时空分布产品；
基于所述第二擦地角对应的目标数据和所述校正后的海杂波半经验计算模型，确定出所述第二擦地角的海杂波时空分布产品。
5.根据权利要求1所述的方法，其特征在于，基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值，包括：
计算出所述海杂波时空分布产品的累计概率密度分布；
基于所述累计概率密度分布，确定出累计概率分割点；
基于所述累计概率分割点，确定出所述海杂波时空分布产品的分级阈值。
6.根据权利要求1所述的方法，其特征在于，基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品，包括：
确定出所述海杂波时空分布产品的海杂波影响强度值；
基于所述海杂波影响强度值和所述分级阈值，对所述海杂波时空分布产品进行分级，得到所述海杂波影响分级产品。
7.一种海杂波时空分布及影响分级产品制作的系统，其特征在于， 包括：获取单元，校正单元，海杂波时空分布产品生成单元，确定单元和海杂波影响分级产品生成单元，其中，
所述获取单元，用于获取目标海域的目标数据和所述目标海域的实测海杂波数据，其中，所述目标数据包括：海况数据，雷达观测参数和所述目标海域的范围数据；
所述校正单元，用于基于所述实测海杂波数据，对海杂波半经验计算模型进行偏差校正，得到校正后的海杂波半经验计算模型；
所述海杂波时空分布产品生成单元，用于基于所述目标数据、所述校正后的海杂波半经验计算模型和复合布拉格后向散射计算模型，确定出所述海杂波时空分布产品；
所述确定单元，用于基于所述海杂波时空分布产品，确定出所述海杂波时空分布产品的分级阈值；
所述海杂波影响分级产品生成单元，用于基于所述分级阈值，对所述海杂波时空分布产品进行分级，得到海杂波影响分级产品。
8.根据权利要求7所述的系统，其特征在于，
所述海况数据包括：海面风速和海面风向；
所述雷达观测参数包括观测频率、观测方位角、擦地角和极化方式。
9.一种电子设备，其特征在于，包括存储器以及处理器，所述存储器用于存储支持处理器执行权利要求1至6任一项所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
10.一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，其特征在于，计算机程序被处理器运行时执行上述权利要求1至6任一项所述方法的步骤。
CN202210382406.5A 2022-04-13 2022-04-13 一种海杂波时空分布及影响分级产品制作的方法及装置 Active CN114488107B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210382406.5A CN114488107B (zh) 2022-04-13 2022-04-13 一种海杂波时空分布及影响分级产品制作的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210382406.5A CN114488107B (zh) 2022-04-13 2022-04-13 一种海杂波时空分布及影响分级产品制作的方法及装置

Publications (2)

Publication Number Publication Date
CN114488107A true CN114488107A (zh) 2022-05-13
CN114488107B CN114488107B (zh) 2022-07-19

Family

ID=81488554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210382406.5A Active CN114488107B (zh) 2022-04-13 2022-04-13 一种海杂波时空分布及影响分级产品制作的方法及装置

Country Status (1)

Country Link
CN (1) CN114488107B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116976091A (zh) * 2023-06-29 2023-10-31 中国空间技术研究院 适用于微波遥感卫星的海浪环境因素影响分析方法

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060238412A1 (en) * 2005-04-21 2006-10-26 Blunt Shannon D Method and apparatus for detecting slow-moving targets in high-resolution sea clutter
JP2011117899A (ja) * 2009-12-07 2011-06-16 Mitsubishi Electric Corp レーダ装置
CN107490790A (zh) * 2017-10-10 2017-12-19 北京航空航天大学 一种连续多脉冲相参海杂波的仿真方法
CN107678003A (zh) * 2017-09-15 2018-02-09 国家海洋局第海洋研究所 一种地波雷达海杂波背景下的目标检测方法和装置
CN109143182A (zh) * 2018-09-10 2019-01-04 南京莱斯电子设备有限公司 一种多尺度的雷达分区图快速处理方法
CN110441752A (zh) * 2019-07-30 2019-11-12 中国人民解放军海军航空大学 海上目标检测方法及装置
CN110907907A (zh) * 2019-10-19 2020-03-24 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种海杂波多普勒谱特性分析对比方法
CN111308433A (zh) * 2020-03-26 2020-06-19 北京航空航天大学 一种包含纹理信息的海杂波仿真方法
CN113065221A (zh) * 2021-02-10 2021-07-02 北京理工大学 大擦地角下海杂波模型、模型构建系统以及模拟仿真方法
CN113189561A (zh) * 2021-06-16 2021-07-30 中国人民解放军火箭军工程大学 一种海杂波参数估计方法、系统、设备及存储介质
CN113281709A (zh) * 2021-04-21 2021-08-20 中国海洋大学 基于区域耦合预报系统的雷达性能评估方法
CN113505463A (zh) * 2021-02-10 2021-10-15 北京理工大学 海杂波与角反复合模型、构建系统及模拟仿真方法

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060238412A1 (en) * 2005-04-21 2006-10-26 Blunt Shannon D Method and apparatus for detecting slow-moving targets in high-resolution sea clutter
JP2011117899A (ja) * 2009-12-07 2011-06-16 Mitsubishi Electric Corp レーダ装置
CN107678003A (zh) * 2017-09-15 2018-02-09 国家海洋局第海洋研究所 一种地波雷达海杂波背景下的目标检测方法和装置
CN107490790A (zh) * 2017-10-10 2017-12-19 北京航空航天大学 一种连续多脉冲相参海杂波的仿真方法
CN109143182A (zh) * 2018-09-10 2019-01-04 南京莱斯电子设备有限公司 一种多尺度的雷达分区图快速处理方法
CN110441752A (zh) * 2019-07-30 2019-11-12 中国人民解放军海军航空大学 海上目标检测方法及装置
CN110907907A (zh) * 2019-10-19 2020-03-24 中国电波传播研究所（中国电子科技集团公司第二十二研究所） 一种海杂波多普勒谱特性分析对比方法
CN111308433A (zh) * 2020-03-26 2020-06-19 北京航空航天大学 一种包含纹理信息的海杂波仿真方法
CN113065221A (zh) * 2021-02-10 2021-07-02 北京理工大学 大擦地角下海杂波模型、模型构建系统以及模拟仿真方法
CN113505463A (zh) * 2021-02-10 2021-10-15 北京理工大学 海杂波与角反复合模型、构建系统及模拟仿真方法
CN113281709A (zh) * 2021-04-21 2021-08-20 中国海洋大学 基于区域耦合预报系统的雷达性能评估方法
CN113189561A (zh) * 2021-06-16 2021-07-30 中国人民解放军火箭军工程大学 一种海杂波参数估计方法、系统、设备及存储介质

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
ANTHONY FICHE等: "Statistical analysis of low grazing angle high resolution X-band SAR sea clutter", 《2014 INTERNATIONAL RADAR CONFERENCE》 *
LUKE ROSENBERG等: "Modeling the Statistics of Microwave Radar Sea Clutter", 《IEEE AEROSPACE AND ELECTRONIC SYSTEMS MAGAZINE》 *
YIWEN WEI等: "Simulation and analysis of the clutter from time varying sea with breaking waves", 《2015 ASIA-PACIFIC MICROWAVE CONFERENCE》 *
丁昊等: "海杂波特性认知研究进展与展望", 《雷达学报》 *
吴涛: "不同海域多频段海杂波特性差异与主要影响因素分析", 《中国优秀博硕士学位论文全文数据库(博士)信息科技辑(月刊)》 *
张勇: "小擦地角海面散射与海杂波研究", 《中国优秀博硕士学位论文全文数据库(硕士)信息科技辑(月刊)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116976091A (zh) * 2023-06-29 2023-10-31 中国空间技术研究院 适用于微波遥感卫星的海浪环境因素影响分析方法
CN116976091B (zh) * 2023-06-29 2024-03-01 中国空间技术研究院 适用于微波遥感卫星的海浪环境因素影响分析方法

Also Published As

Publication number Publication date
CN114488107B (zh) 2022-07-19

Similar Documents

Publication Publication Date Title
US9097805B2 (en) Systems and methods for calibrating dual polarization radar systems
Cho et al. Identification and removal of ground echoes and anomalous propagation using the characteristics of radar echoes
US9261583B2 (en) Systems and methods for calibrating dual polarization radar systems
CN110095760A (zh) 一种气象雷达的测试装置和方法
CN103064073A (zh) 一种基于频率捷变改变雷达目标特性的方法
Richardson et al. Bragg scatter detection by the WSR-88D. Part I: Algorithm development
CN114488107B (zh) 一种海杂波时空分布及影响分级产品制作的方法及装置
CN113281709B (zh) 基于区域耦合预报系统的雷达性能评估方法
CN102405422B (zh) 在声雷达系统中检测和补偿降水的方法
Myagkov et al. Evaluation of the reflectivity calibration of W-band radars based on observations in rain
CN114397636A (zh) 一种地基雷达反射率因子均一性评估方法、系统及设备
Anderson et al. Analysis of ASCAT ocean backscatter measurement noise
Berne et al. Influence of the vertical profile of reflectivity on radar-estimated rain rates at short time steps
Wyatt et al. Factors affecting the accuracy of SHOWEX HF radar wave measurements
Kuster et al. Using ZDR columns in forecaster conceptual models and warning decision-making
Huuskonen et al. Radar performance monitoring using the angular width of the solar image
Asai et al. Validation of X-band multiparameter phased-array weather radar by comparing data from Doppler weather radar with a parabolic dish antenna
Unal High-resolution raindrop size distribution retrieval based on the Doppler spectrum in the case of slant profiling radar
CN109541565B (zh) 一种雷达回波强度均一性检测方法及系统
CN115546526B (zh) 三维点云聚类方法、装置及存储介质
CN114814779B (zh) 浮标涌浪波高观测数据误差评估方法、系统、设备及介质
Sun et al. Remote sensing of ocean surface wind direction with shipborne high frequency surface wave radar
CN112906579B (zh) 基于k均值聚类与svm的海杂波弱目标分类方法及系统
CN113219466A (zh) 海面的风速确定方法、装置、电子设备及存储介质
Østenstad et al. Evaporation duct height climatology for Norwegian waters using hindcast data

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