CN114415181B - 一种合成孔径雷达的原始回波生成方法和装置 - Google Patents

一种合成孔径雷达的原始回波生成方法和装置 Download PDF

Info

Publication number
CN114415181B
CN114415181B CN202210327750.4A CN202210327750A CN114415181B CN 114415181 B CN114415181 B CN 114415181B CN 202210327750 A CN202210327750 A CN 202210327750A CN 114415181 B CN114415181 B CN 114415181B
Authority
CN
China
Prior art keywords
target
azimuth direction
echo data
processed
sea surface
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
CN202210327750.4A
Other languages
English (en)
Other versions
CN114415181A (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 CN202210327750.4A priority Critical patent/CN114415181B/zh
Publication of CN114415181A publication Critical patent/CN114415181A/zh
Application granted granted Critical
Publication of CN114415181B publication Critical patent/CN114415181B/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/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
    • 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
    • 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/415—Identification of targets based on measurements of movement associated with the target
    • 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)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种合成孔径雷达的原始回波生成方法和装置，涉及海洋工程的技术领域，包括：获取待处理海域的海洋场景仿真参数；获取待处理海域的目标方位向的仿真时刻，基于海洋场景仿真参数和目标方位向的仿真时刻，计算出时变传播系数，以及基于时变传播系数，计算出目标参数；基于目标参数和散射模型，构建目标方位向的目标矩阵；基于目标矩阵，对目标方位向的目标回波强度值进行合并，得到目标方位向的原始回波数据；重复执行上述步骤，直至得到处理海域的所有方位向的原始回波数据，并基于所有方位向的原始回波数据，确定出待处理海域的原始回波数据，解决了现有的原始回波生成方法难以生成高精度原始回波数据的技术问题。

Description

一种合成孔径雷达的原始回波生成方法和装置
技术领域
本发明涉及海洋工程的技术领域，尤其是涉及一种合成孔径雷达的原始回波生成方法和装置。
背景技术
合成孔径雷达采用相干成像的方式获取目标信息，是对陆对海观测领域的重要传感器。在海洋学及海洋观测领域，SAR更是应用广泛的传感器，由于合成孔径雷达图像中包含十分丰富的海洋信息，以及其全天候、全天时对海面进行连续观测的特点，使得其在海洋理论与应用研究领域发挥着越来越重要的作用。合成孔径雷达原始回波数据的仿真，有助于评估载荷性能及改进并优化成像处理算法，是合成孔径雷达进行海洋工程探测领域中的一项重要的处理过程。
与陆地固定目标的合成孔径雷达仿真不同，海表面随着时间不停变化。而合成孔径雷达是通过相干积分方式成像的，在整个积分时间内，海面运动会造成海面后向散射是分辨率单元内散射小散射元的运动效应在积分时间内的积分关系。传统的海面仿真成像，多为基于静止的二维或三维海面，在相干积分时间内，没有考虑海面运动效应，这就会使得最终仿真的原始回波数据与真实雷达工作情形存在偏差，对载荷性能评估及算法验证产生影响。
针对上述问题，还未提出有效的解决方案。
发明内容
有鉴于此，本发明的目的在于提供一种合成孔径雷达的原始回波生成方法和装置，以缓解了现有的原始回波生成方法难以生成高精度原始回波数据的技术问题。
第一方面，本发明实施例提供了一种合成孔径雷达的原始回波生成方法，包括：获取步骤，获取待处理海域的海洋场景仿真参数；计算步骤，获取所述待处理海域的目标方位向的仿真时刻，基于所述海洋场景仿真参数和所述目标方位向的仿真时刻，计算出时变传播系数，以及基于所述时变传播系数，计算出目标参数，其中，所述目标参数包括：海表面高度和海表面斜率，所述目标方位向为所述待处理海域的任意一个方位向；构建步骤，基于所述目标参数和散射模型，构建所述目标方位向的目标矩阵，其中，所述目标矩阵包括：所述目标方位向的后向散射系数矩阵和所述后向散射系数矩阵对应的斜距矩阵；合并步骤，基于所述目标矩阵，对所述目标方位向的目标回波强度值进行合并，得到所述目标方位向的原始回波数据，其中，所述目标回波强度为所述目标方位向的回波强度值中斜距相同的回波强度值；重复执行所述计算步骤、所述构建步骤和所述合并步骤，直至得到所述待处理海域的所有方位向的原始回波数据，并基于所述所有方位向的原始回波数据，确定出所述待处理海域的原始回波数据。
进一步地，所述海洋场景仿真参数包括：海面尺度、海面分辨率、风速和风向。
进一步地，所述时变传播系数的计算公式为
Figure P_220329090442515_515605001
，其中，
Figure P_220329090442547_547335002
为所述时变传播系数，
Figure P_220329090442562_562949003
和
Figure P_220329090442578_578624004
为服从高斯分布的随机数，
Figure P_220329090442611_611321005
为海面海浪谱，
Figure P_220329090442626_626959006
为所述目标方位向的仿真时刻，
Figure P_220329090442642_642588007
为复数，波数向量
Figure P_220329090442673_673841008
，
Figure P_220329090442689_689468009
，
Figure P_220329090442720_720699010
，
Figure P_220329090442736_736352011
，
Figure P_220329090442751_751962012
为重力加速度常数，
Figure P_220329090442784_784109013
为所述风速，
Figure P_220329090442800_800247014
为所述风向，所述海面尺度为
Figure P_220329090442815_815920015
、所述海面分辨率为
Figure P_220329090442847_847125016
。
进一步地，所述海表面高度的计算公式为
Figure P_220329090442862_862763001
；
所述海表面斜率为
Figure P_220329090442894_894010001
，所述海表面斜率的计算公式
Figure P_220329090442925_925251002
，
Figure P_220329090442940_940898003
；其中，
Figure P_220329090442972_972135004
为取实部，
Figure P_220329090442989_989169005
为傅里叶逆变换。
进一步地，所述的散射模型为复合Brag后向散射系数计算模型。
进一步地，基于所述目标矩阵，对所述目标方位向的目标回波强度值进行合并，得到所述目标方位向的原始回波数据，包括：利用sinc函数构造插值核函数；利用所述插值核函数、所述后向散射系数矩阵和所述斜距矩阵，对所述目标回波强度值进行合并处理，得到所述目标方位向的原始回波数据。
进一步地，基于所述所有方位向的原始回波数据，确定出所述待处理海域的原始回波数据，包括：利用加噪处理公式，对所述所有方位向的原始回波数据进行加噪处理，得到所述待处理海域的原始回波数据，其中，所述加噪处理公式
Figure P_220329090443004_004866001
，
Figure P_220329090443021_021018002
为所述待处理海域的原始回波数据，
Figure P_220329090443051_051728003
为所有方位向的原始回波数据，
Figure P_220329090443067_067829004
为等效噪声系数，
Figure P_220329090443082_082993005
和
Figure P_220329090443099_099119006
为服从高斯分布的随机数。
进一步地，所述待处理海域的所有方位向的数量的计算公式为：
Figure P_220329090443129_129848001
；其中，
Figure P_220329090443145_145487002
为所述待处理海域的所有方位向的数量，
Figure P_220329090443161_161586003
为四舍五入取整运算，
Figure P_220329090443177_177220004
为最近斜距，
Figure P_220329090443212_212357005
为天线长度，
Figure P_220329090443227_227989006
为所述待处理海域的方位长度，
Figure P_220329090443243_243618007
为地速，
Figure P_220329090443259_259243008
为合成孔径雷达的脉冲重复频率；所述目标方位向的仿真时刻的计算公式为：
Figure P_220329090443290_290498009
，其中，
Figure P_220329090443306_306132010
为所述目标方位向时间，
Figure P_220329090443321_321759011
为所述目标方位向的序号，取值范围为1至
Figure P_220329090443352_352990012
。
第二方面，本发明实施例还提供了一种合成孔径雷达的原始回波生成装置，包括：所述获取单元，用于获取待处理海域的海洋场景仿真参数；所述计算单元，获取所述待处理海域的目标方位向的仿真时刻，基于所述海洋场景仿真参数和所述目标方位向的仿真时刻，计算出时变传播系数，以及基于所述时变传播系数，计算出目标参数，其中，所述目标参数包括：海表面高度和海表面斜率，所述目标方位向为所述待处理海域的任意一个方位向；所述构建单元，基于所述目标参数和散射模型，构建所述目标方位向的目标矩阵，其中，所述目标矩阵包括：所述目标方位向的后向散射系数矩阵和所述后向散射系数矩阵对应的斜距矩阵；所述合并单元，基于所述目标矩阵，对所述目标方位向的目标回波强度值进行合并，得到所述目标方位向的原始回波数据，其中，所述目标回波强度为所述目标方位向的回波强度值中斜距相同的回波强度值；所述处理单元，用于控制所述计算单元、所述构建单元和所述合并单元重复工作，直至得到所述待处理海域的所有方位向的原始回波数据，并基于所述所有方位向的原始回波数据，确定出所述待处理海域的原始回波数据。
第三方面，本发明实施例还提供了一种电子设备，包括存储器以及处理器，所述存储器用于存储支持处理器执行上述第一方面中所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
第四方面，本发明实施例还提供了一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序。
在本发明实施例中，通过计算时变传播系数，计算海表面高度及海表面斜率；基于海表面高度及海表面斜率，结合散射模型计算当前方位向下场景后向散射系数；基于场景后向散射系数，实现相同斜距合并处理，得到单测线的原始回波数据；遍历所有方位向位置，同步更新二维动态海面，最终生成整个场景的原始回波数据，达到了基于时变海面的合成孔径雷达原始回波数据生成的目的，解决现有技术中因缺少对海面时变因素考量导致仿真结果不精确的技术问题，从而实现了能够基于时变动态海洋的原始回波数据仿真的技术效果。
本发明的其他特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附附图，作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案，下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种合成孔径雷达的原始回波生成方法的流程图 ；
图2为本发明实施例提供的一种合成孔径雷达的原始回波生成装置的示意图 ；
图3为本发明实施例提供的一种电子设备的示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一：
根据本发明实施例，提供了一种合成孔径雷达的原始回波生成方法的实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的一种合成孔径雷达的原始回波生成方法的流程图，如图1所示，该方法包括如下步骤：
步骤S102，获取步骤，获取待处理海域的海洋场景仿真参数；
需要说明的是，海洋场景仿真参数包括：海面尺度、海面分辨率、风速和风向。
步骤S104，计算步骤，获取所述待处理海域的目标方位向的仿真时刻，基于所述海洋场景仿真参数和所述目标方位向的仿真时刻，计算出时变传播系数，以及基于所述时变传播系数，计算出目标参数，其中，所述目标参数包括：海表面高度和海表面斜率，所述目标方位向为所述待处理海域的任意一个方位向；
步骤S106，构建步骤，基于所述目标参数和散射模型，构建所述目标方位向的目标矩阵，其中，所述目标矩阵包括：所述目标方位向的后向散射系数矩阵和所述后向散射系数矩阵对应的斜距矩阵；
需要说明的是，优选的，上述的散射模型采用复合Brag后向散射系数计算模型。
步骤S108，合并步骤，基于所述目标矩阵，对所述目标方位向的目标回波强度值进行合并，得到所述目标方位向的原始回波数据，其中，所述目标回波强度为所述目标方位向的回波强度值中斜距相同的回波强度值；
步骤S110，重复执行所述计算步骤、所述构建步骤和所述合并步骤，直至得到所述待处理海域的所有方位向的原始回波数据，并基于所述所有方位向的原始回波数据，确定出所述待处理海域的原始回波数据。
在本发明实施例中，通过计算时变传播系数，计算海表面高度及海表面斜率；基于海表面高度及海表面斜率，结合散射模型计算当前方位向下场景后向散射系数；基于场景后向散射系数，实现相同斜距合并处理，得到单测线的原始回波数据；遍历所有方位向位置，同步更新二维动态海面，最终生成整个场景的原始回波数据，达到了基于时变海面的合成孔径雷达原始回波数据生成的目的，解决现有技术中因缺少对海面时变因素考量导致仿真结果不精确的技术问题，从而实现了能够基于时变动态海洋的原始回波数据仿真的技术效果。
在本发明实施例中变传播系数的计算公式为
Figure P_220329090443368_368665001
，其中，
Figure P_220329090443385_385197002
为所述时变传播系数，
Figure P_220329090443416_416989003
和
Figure P_220329090443432_432581004
为服从高斯分布的随机数，
Figure P_220329090443448_448242005
为海面海浪谱，
Figure P_220329090443463_463849006
为所述目标方位向的仿真时刻，
Figure P_220329090443495_495073007
为复数，波数向量
Figure P_220329090443510_510747008
，
Figure P_220329090443541_541965009
，
Figure P_220329090443557_557570010
，
Figure P_220329090443573_573235011
，
Figure P_220329090443606_606413012
为重力加速度常数，
Figure P_220329090443622_622041013
为所述风速，
Figure P_220329090443637_637659014
为所述风向，所述海面尺度为
Figure P_220329090443653_653288015
、所述海面分辨率为
Figure P_220329090443684_684568016
。所述海表面高度的计算公式为
Figure M_220329090443699_699701001
；
所述海表面高度的计算公式为
Figure P_220329090443747_747064001
；
所述海表面斜率为
Figure P_220329090443778_778321001
，所述海表面斜率的计算公式
Figure P_220329090443795_795350002
，
Figure P_220329090443827_827158003
；
其中，
Figure P_220329090443842_842771001
为取实部，
Figure P_220329090443874_874024002
为傅里叶逆变换。
在本发明实施例中，步骤S106包括如下步骤：
利用所述目标参数，结合外部输入的雷达参数，其中，雷达参数包括入射角，方位角，工作频率，极化方式。
利用复合Brag后向散射系数计算模型，计算目标方位向的后向散射系数，并根据目标方位向的后向散射系数构建目标方位向的目标矩阵。
在本发明实施例中，步骤S108包括如下步骤：
步骤S11，利用sinc函数构造插值核函数；
步骤S12，利用所述插值核函数、所述后向散射系数矩阵和所述斜距矩阵，对所述目标回波强度值进行合并处理，得到所述目标方位向的原始回波数据。
插值核函数采用
Figure P_220329090443889_889649001
函数，原始长度为
Figure P_220329090443920_920916002
，过采样因子为
Figure P_220329090443936_936515003
，生成的
Figure P_220329090443967_967774004
序列为
Figure P_220329090443984_984325005
；
（2）合并处理计算方法为：
Figure P_220329090444000_000457001
其中，
Figure P_220329090444047_047343001
为斜距总点数，
Figure P_220329090444062_062982002
为未过采样的
Figure P_220329090444078_078588003
函数长度，
Figure P_220329090444094_094219004
为所述后向散射系数矩阵，
Figure P_220329090444125_125444005
为与所述后向散射系数矩阵对应的斜距矩阵，
Figure P_220329090444141_141102006
为向下取整运算，
Figure P_220329090444156_156692007
为四舍五入取整运算，
Figure P_220329090444172_172339008
为取余运算，
Figure P_220329090444206_206029009
，
Figure P_220329090444221_221649010
，
Figure P_220329090444237_237281011
为中间变量，
Figure P_220329090444268_268523012
为目标方位向的原始回波数据。
在本发明实施例中，在得到所有方位向的原始回波数据之后，需要利用加噪处理公式，对所述所有方位向的原始回波数据进行加噪处理，得到所述待处理海域的原始回波数据，其中，所述加噪处理公式
Figure P_220329090444284_284144001
，
Figure P_220329090444299_299763002
为所述待处理海域的原始回波数据，
Figure P_220329090444331_331046003
为所有方位向的原始回波数据，
Figure P_220329090444346_346634004
为等效噪声系数，
Figure P_220329090444362_362273005
和
Figure P_220329090444377_377902006
为服从高斯分布的随机数。
需要说明的是，在第一次执行获取步骤之前，需要计算出待处理海域的所有方位向的数量和目标方位向的仿真时刻，待处理海域的所有方位向的数量的计算公式为：
Figure P_220329090444398_398365001
；
其中，
Figure P_220329090444430_430167001
为所述待处理海域的所有方位向的数量，
Figure P_220329090444445_445769002
为四舍五入取整运算，
Figure P_220329090444461_461379003
为最近斜距，
Figure P_220329090444477_477008004
为天线长度，
Figure P_220329090444508_508285005
为所述待处理海域的方位长度，
Figure P_220329090444523_523934006
为地速，
Figure P_220329090444539_539530007
为合成孔径雷达的脉冲重复频率；
所述目标方位向的仿真时刻的计算公式为：
Figure P_220329090444555_555234001
，其中，
Figure P_220329090444586_586411002
为所述目标方位向时间，
Figure P_220329090444603_603452003
为所述目标方位向的序号，取值范围为1至
Figure P_220329090444619_619621004
。
另外，还需要说明的是，在每次获取到目标参数之后，可以基于目标参数，构建二维动态海面。
实施例二：
本发明实施例还提供了一种合成孔径雷达的原始回波生成装置，该装置用于执行本发明实施例上述内容所提供的合成孔径雷达的原始回波生成方法，以下是本发明实施例提供的装置的具体介绍。
如图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 (8)

1.一种合成孔径雷达的原始回波生成方法，其特征在于，包括：
获取步骤，获取待处理海域的海洋场景仿真参数；
计算步骤，获取所述待处理海域的目标方位向的仿真时刻，基于所述海洋场景仿真参数和所述目标方位向的仿真时刻，计算出时变传播系数，以及基于所述时变传播系数，计算出目标参数，其中，所述目标参数包括：海表面高度和海表面斜率，所述目标方位向为所述待处理海域的任意一个方位向；
构建步骤，基于所述目标参数和散射模型，构建所述目标方位向的目标矩阵，其中，所述目标矩阵包括：所述目标方位向的后向散射系数矩阵和所述后向散射系数矩阵对应的斜距矩阵，所述的散射模型为复合Brag后向散射系数计算模型；
合并步骤，基于所述目标矩阵，对所述目标方位向的目标回波强度值进行合并，得到所述目标方位向的原始回波数据，其中，所述目标回波强度为所述目标方位向的回波强度值中斜距相同的回波强度值；
重复执行所述计算步骤、所述构建步骤和所述合并步骤，直至得到所述待处理海域的所有方位向的原始回波数据，并基于所述所有方位向的原始回波数据，确定出所述待处理海域的原始回波数据；
其中，所述海洋场景仿真参数包括：海面尺度、海面分辨率、风速和风向；
所述时变传播系数的计算公式为
Figure P_220512170029088_088698001
，其中，
Figure P_220512170029119_119946002
为所述时变传播系数，
Figure P_220512170029151_151237003
和
Figure P_220512170029166_166862004
为服从高斯分布的随机数，
Figure P_220512170029198_198072005
为海面海浪谱，
Figure P_220512170029248_248869006
为所述目标方位向的仿真时刻，
Figure P_220512170029280_280099007
为复数，波数向量
Figure P_220512170029311_311365008
，
Figure P_220512170029326_326992009
，
Figure P_220512170029358_358233010
，
Figure P_220512170029389_389480011
，
Figure P_220512170029424_424636012
为重力加速度常数，
Figure P_220512170029455_455905013
为所述风速，
Figure P_220512170029487_487146014
为所述风向，所述海面尺度为
Figure P_220512170029502_502817015
、所述海面分辨率为
Figure P_220512170029534_534028016
。
2.根据权利要求1所述的方法，其特征在于，
所述海表面高度的计算公式为
Figure P_220512170029580_580900001
；
所述海表面斜率为
Figure P_220512170029613_613578001
，所述海表面斜率的计算公式为
Figure P_220512170029645_645353002
，
Figure P_220512170029676_676607003
；
其中，
Figure P_220512170029707_707842001
为取实部，
Figure P_220512170029739_739088002
为傅里叶逆变换。
3.根据权利要求1所述的方法，其特征在于，基于所述目标矩阵，对所述目标方位向的目标回波强度值进行合并，得到所述目标方位向的原始回波数据，包括：
利用sinc函数构造插值核函数；
利用所述插值核函数、所述后向散射系数矩阵和所述斜距矩阵，对所述目标回波强度值进行合并处理，得到所述目标方位向的原始回波数据。
4.根据权利要求1所述的方法，其特征在于，基于所述所有方位向的原始回波数据，确定出所述待处理海域的原始回波数据，包括：
利用加噪处理公式，对所述所有方位向的原始回波数据进行加噪处理，得到所述待处理海域的原始回波数据，其中，所述加噪处理公式
Figure P_220512170029770_770353001
，
Figure P_220512170029803_803535002
为所述待处理海域的原始回波数据，
Figure P_220512170029944_944669003
为所有方位向的原始回波数据，
Figure P_220512170030164_164382004
为等效噪声系数，
Figure P_220512170030308_308439005
和
Figure P_220512170030336_336270006
为服从高斯分布的随机数。
5.根据权利要求1所述的方法，其特征在于，
所述待处理海域的所有方位向的数量的计算公式为：
Figure P_220512170030367_367512001
；
其中，
Figure P_220512170030398_398767001
为所述待处理海域的所有方位向的数量，
Figure P_220512170030416_416304002
为四舍五入取整运算，
Figure P_220512170030448_448084003
为最近斜距，
Figure P_220512170030479_479333004
为天线长度，
Figure P_220512170030494_494949005
为所述待处理海域的方位长度，
Figure P_220512170030648_648783006
为地速，
Figure P_220512170030805_805012007
为合成孔径雷达的脉冲重复频率；
所述目标方位向的仿真时刻的计算公式为：
Figure P_220512170030858_858770001
，其中，
Figure P_220512170031061_061328002
为所述目标方位向时间，
Figure P_220512170031329_329936003
为所述目标方位向的序号，取值范围为1至
Figure P_220512170031509_509118004
。
6.一种合成孔径雷达的原始回波生成装置，其特征在于，包括：获取单元，计算单元，构建单元，合并单元和处理单元，其中，
所述获取单元，用于获取待处理海域的海洋场景仿真参数；
所述计算单元，获取所述待处理海域的目标方位向的仿真时刻，基于所述海洋场景仿真参数和所述目标方位向的仿真时刻，计算出时变传播系数，以及基于所述时变传播系数，计算出目标参数，其中，所述目标参数包括：海表面高度和海表面斜率，所述目标方位向为所述待处理海域的任意一个方位向；
所述构建单元，基于所述目标参数和散射模型，构建所述目标方位向的目标矩阵，其中，所述目标矩阵包括：所述目标方位向的后向散射系数矩阵和所述后向散射系数矩阵对应的斜距矩阵，所述的散射模型为复合Brag后向散射系数计算模型；
所述合并单元，基于所述目标矩阵，对所述目标方位向的目标回波强度值进行合并，得到所述目标方位向的原始回波数据，其中，所述目标回波强度为所述目标方位向的回波强度值中斜距相同的回波强度值；
所述处理单元，用于控制所述计算单元、所述构建单元和所述合并单元重复工作，直至得到所述待处理海域的所有方位向的原始回波数据，并基于所述所有方位向的原始回波数据，确定出所述待处理海域的原始回波数据；
其中，所述海洋场景仿真参数包括：海面尺度、海面分辨率、风速和风向；
所述时变传播系数的计算公式为
Figure P_220512170031524_524749001
，其中，
Figure P_220512170031555_555996002
为所述时变传播系数，
Figure P_220512170031587_587259003
和
Figure P_220512170031604_604297004
为服从高斯分布的随机数，
Figure P_220512170031636_636056005
为海面海浪谱，
Figure P_220512170031651_651680006
为所述目标方位向的仿真时刻，
Figure P_220512170031682_682934007
为复数，波数向量
Figure P_220512170031714_714185008
，
Figure P_220512170031792_792320009
，
Figure P_220512170031827_827947010
，
Figure P_220512170031859_859223011
，
Figure P_220512170031890_890475012
为重力加速度常数，
Figure P_220512170031906_906093013
为所述风速，
Figure P_220512170031937_937333014
为所述风向，所述海面尺度为
Figure P_220512170031968_968585015
、所述海面分辨率为
Figure P_220512170031999_999827016
。
7.一种电子设备，其特征在于，包括存储器以及处理器，所述存储器用于存储支持处理器执行权利要求1至5任一项所述方法的程序，所述处理器被配置为用于执行所述存储器中存储的程序。
8.一种计算机可读存储介质，计算机可读存储介质上存储有计算机程序，其特征在于，计算机程序被处理器运行时执行上述权利要求1至5任一项所述方法的步骤。
CN202210327750.4A 2022-03-31 2022-03-31 一种合成孔径雷达的原始回波生成方法和装置 Active CN114415181B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210327750.4A CN114415181B (zh) 2022-03-31 2022-03-31 一种合成孔径雷达的原始回波生成方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210327750.4A CN114415181B (zh) 2022-03-31 2022-03-31 一种合成孔径雷达的原始回波生成方法和装置

Publications (2)

Publication Number Publication Date
CN114415181A CN114415181A (zh) 2022-04-29
CN114415181B true CN114415181B (zh) 2022-07-12

Family

ID=81263684

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210327750.4A Active CN114415181B (zh) 2022-03-31 2022-03-31 一种合成孔径雷达的原始回波生成方法和装置

Country Status (1)

Country Link
CN (1) CN114415181B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115629364B (zh) * 2022-12-22 2023-03-28 中国海洋大学 一种面向动态海面的星载小角度sar海况偏差仿真方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102144173B (zh) * 2008-05-29 2014-08-20 电视广播有限公司 合成孔径雷达成像的海域中的目标检测
CN105445711B (zh) * 2015-11-27 2017-08-01 南京信息工程大学 一种基于逆Omega‑K算法的海面要素SAR原始数据仿真方法
CN107229051B (zh) * 2017-05-26 2020-05-01 西安电子科技大学 基于gpu的视频sar回波仿真并行实现方法
CN107271977B (zh) * 2017-07-25 2020-04-24 哈尔滨工业大学 基于移动激励源fdtd算法的高精度sar回波仿真方法
CN107656253B (zh) * 2017-08-25 2019-12-13 北京航空航天大学 电磁涡旋合成孔径雷达回波信号仿真方法及装置
CN110208764A (zh) * 2019-05-05 2019-09-06 南京航空航天大学 基于电磁散射计算的动态海洋大型场景回波仿真方法
CN110456348B (zh) * 2019-08-19 2020-08-25 中国石油大学(华东) 多视向sar海浪谱数据融合的海浪截断波长补偿方法
CN110988878B (zh) * 2019-11-28 2023-04-07 中国海洋大学 一种基于rd算法的sar海浪成像仿真方法
CN111522007A (zh) * 2020-07-06 2020-08-11 航天宏图信息技术股份有限公司 真实场景与目标仿真融合的sar成像仿真方法和系统
CN114047511B (zh) * 2021-11-02 2022-04-12 中国海洋大学 一种基于csa算法的时变海面机载sar成像仿真方法

Also Published As

Publication number Publication date
CN114415181A (zh) 2022-04-29

Similar Documents

Publication Publication Date Title
Sutherland et al. Evaluation of coastal area modelling systems at an estuary mouth
CN106872958B (zh) 基于线性融合的雷达目标自适应检测方法
JP2012533051A (ja) 同一領域上で取得されたｓａｒ画像から得られるインターフェログラムのフィルタリング処理方法
KR102319145B1 (ko) 고해상도 해양 데이터를 생성하기 위한 방법 및 이를 이용한 장치
CN114415181B (zh) 一种合成孔径雷达的原始回波生成方法和装置
US20190227164A1 (en) Automotive Testing Method, System and Computer Program Product
CN111611731B (zh) 卫星数据融合方法、装置及电子设备
CN103325092B (zh) 一种生成二维相位解缠质量图的方法及装置
Qi et al. Phase-resolved wave field simulation calibration of sea surface reconstruction using noncoherent marine radar
CN115079109A (zh) 一种合成孔径雷达成像质量指标评价系统
CN111487621A (zh) 一种基于雷达图像的海表流场反演方法及电子设备
Hammer et al. Coherent simulation of SAR images
CN116680994B (zh) 一种基于激光雷达的气溶胶追踪和风场反演方法及系统
CN114488107B (zh) 一种海杂波时空分布及影响分级产品制作的方法及装置
CN115877345B (zh) 一种风廓线雷达缺测数据的补齐方法和装置
CN114814779B (zh) 浮标涌浪波高观测数据误差评估方法、系统、设备及介质
CN107271981B (zh) 一种目标电磁散射数据模型的建立方法
CN113064130B (zh) 粒子谱分布的确定方法、装置、存储介质及程序产品
Lu et al. Research on rainfall identification based on the echo differential value from X-band navigation radar image
CN112733801B (zh) 海洋内波发生频率的计算方法、装置和电子设备
Chen et al. Sand ripple characterization using an extended synthetic aperture sonar model and parallel sampling method
Di Martino et al. Physical models for SAR speckle simulation
CN116359877B (zh) 一种基于风廓线雷达的大雾天气监测方法和装置
Daniel et al. Millimetre-wave and sub-THz maritime radar scene simulation
CN112394355B (zh) 圆周扫描地基sar的三维成像旁瓣抑制方法及装置

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