CN111679258A - 一种雷达回波数据处理方法及装置、设备 - Google Patents

一种雷达回波数据处理方法及装置、设备 Download PDF

Info

Publication number
CN111679258A
CN111679258A CN202010421677.8A CN202010421677A CN111679258A CN 111679258 A CN111679258 A CN 111679258A CN 202010421677 A CN202010421677 A CN 202010421677A CN 111679258 A CN111679258 A CN 111679258A
Authority
CN
China
Prior art keywords
data
radar
faraday rotation
echo
matrix
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
CN202010421677.8A
Other languages
English (en)
Other versions
CN111679258B (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.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
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 Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202010421677.8A priority Critical patent/CN111679258B/zh
Publication of CN111679258A publication Critical patent/CN111679258A/zh
Application granted granted Critical
Publication of CN111679258B publication Critical patent/CN111679258B/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
    • 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
    • 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
    • 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

本发明实施例公开了一种雷达数据处理方法及装置、设备，其中，所述方法包括，根据雷达原始回波数据获取电离层中的电子总量TEC；根据TEC确定雷达原始回波数据对应的法拉第旋转角数据；法拉第旋转角数据随频率非线性变化；根据法拉第旋转角数据确定PD的影响量化值；若影响量化值大于或等于预设阈值，则根据法拉第旋转角数据对雷达原始回波数据进行校正。根据本申请实施例提供的雷达数据处理方法及装置、设备，可以对雷达原始回波数据中PD影响进行校正；提高校正效果及校正准确性。

Description

一种雷达回波数据处理方法及装置、设备
技术领域
本发明实施例涉及雷达信号处理技术领域，尤其涉及一种雷达回波数据方 法及装置、设备。
背景技术
星载P波段合成孔径雷达(SAR，Synthetic Aperture Radar)由于兼具星载 SAR全天时、全天候、大范围测绘的优势和P波段穿透性强、对森林生物量敏 感的特性，在隐蔽目标侦察、植被资源勘测等军事、民用方面具有独特的优势 和巨大的需求。
对星载SAR来说，电离层色散、闪烁、法拉第旋转(FR，Faraday Rotation)、 极化色散(PD，Polarimetric Dispersion)等传播效应会在电磁波的幅度、相位、 极化等方面引入误差，严重影响星载P波段SAR的性能(图像分辨率、峰值旁 瓣比、积分旁瓣比等)，制约了星载P波段宽带SAR***的发展。
传统技术中对雷达回波数据中的电离层法拉第旋转影响的校正还集中在 SAR成像后对FR旋转矩阵的估计与校正，没有涉及到极化色散影响的校正， 极化色散仍是影响星载P波段SAR的图像性能的技术瓶颈。
发明内容
本发明实施例期望提供一种雷达回波数据处理方法及装置、设备，能够对 雷达回波数据中的极化色散影响进行校正。
本发明的技术方案是这样实现的：
一种雷达回波数据处理方法，包括：根据雷达原始回波数据获取电离层中 的电子总量TEC；根据所述TEC确定所述雷达原始回波数据对应的法拉第旋转 角数据；所述法拉第旋转角数据随频率非线性变化；根据所述法拉第旋转角数 据确定PD的影响量化值；若所述影响量化值大于或等于预设阈值，则根据所 述法拉第旋转角数据对所述雷达原始回波数据进行校正。
一种雷达回波数据处理装置，装置包括：第一获取模块，用于根据雷达原 始回波数据获取电离层中的电子总量TEC，并根据电子总量确定雷达原始回波 数据对应的法拉第旋转角数据；第一确定模块，用于根据法拉第旋转角数据确 定影响量化值；校正模块，用于若影响量化值大于或等于预设阈值，则根据法 拉第旋转角数据对雷达原始回波数据进行校正。
一种雷达回波数据处理设备，包括：存储器，用于存储可执行数据指令； 处理器，用于执行存储器中存储的可执行指令时，实现上述的雷达回波数据处 理方法。
一种计算机可读存储介质，存储有可执行指令，用于引起处理器执行时， 实现上述的雷达回波数据处理方法。
在本申请实施例提供的雷达数据处理方法中，通过根据雷达原始回波数据 获取电离层中的电子总量TEC；根据TEC确定雷达原始回波数据对应的法拉第 旋转角数据；法拉第旋转角数据随频率非线性变化；根据法拉第旋转角数据确 定PD的影响量化值；若影响量化值大于或等于预设阈值，则根据法拉第旋转 角数据对雷达原始回波数据进行校正。根据本申请实施例提供的雷达数据处理 方法，由于在估计法拉第旋转角数据时考虑了PD对雷达原始回波数据的影响， 将该法拉第旋转角数据定义为随频率非线性变化的数据，并非一个固定值，在 进行后续校正时，可以根据该法拉第旋转角数据对雷达原始回波数据中PD影响进行校正；由于在对雷达原始回波数据进行校正之前，先根据预设阈值对PD 的影响量化值进行判断，可以保证只在PD影响较大时对雷达原始回波数据进 行校正，避免了PD影响较小时对雷达原始回波数据进行反向校正的问题。
附图说明
图1为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图2为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图3为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图4为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图5为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图6为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图7为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图8为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示意 图；
图9a至图9k为本发明实施例提供的雷达回波数据处理方法的SAR点目标 图像示意图；
图10为本发明实施例提供的雷达回波数据处理方法的一个可选的流程示 意图；
图11为本发明实施例提供的雷达回波数据处理装置的部分结构示意图；
图12为本发明实施例提供的雷达回波数据处理设备的部分结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清 楚、完整地描述。
图1是本发明实施例提供的方法的一个可选的流程示意图，将结合图1示 出的步骤进行说明。
在S102中，根据雷达原始回波数据获取电离层中的电子总量TEC。
其中，星载合成孔径雷达(Synthetic Aperture Radar,SAR)可全天时、全天 候获取大范围、高分辨图像，在军事侦察和经济建设中有着重要的应用价值。 电离层色散、闪烁和法拉第旋转(Faraday Rotation,FR)效应会在电磁波信号 幅度、相位、极化等方面引入误差，造成星载SAR***成像、定位、极化测 量等方面性能下降。本申请中获取的雷达原始回波数据为未经过处理的雷达回 波数据，由于在传输过程中受到电离层色散、闪烁和FR效应的影响，需要对 其进行校正。
其中，TEC为电离层电子浓度总含量，又称电离层电子浓度柱含量、积分 含量等，是一个非常重要的电离层参量，对电离层物理的理论研究及电离层电 波传播的应用研究均具有十分重要的意义。理论上，TEC的空间分布及时间变 化，反映了电离层的主要特性，因此通过探测与分析电离层TEC参量，可以研 究电离层不同时空尺度的分布与变化特性，如电离层扰动，电离层的周日、逐 日变化，电离层年度变化，以及电离层的长期变化等。
在一个实施例中，可以通过专用仪器设备直接测量该雷达成像时电离层中 的电子总量TEC含量，此种方法由于采用了专用仪器设备，其测量精度较高， 但时效性较差，有测量延时，不能实时反应雷达信号穿过电离层时刻电离层对 电磁波的影响。
在另一个实施例中，可以通过雷达原始回波数据估计TEC，在实施时，可 以通过Bickel and Bates方法、Freeman方法、Qi and Jin方法、Chen and Quegan 方法。这些方法可以从雷达原始回波数据中估计电离层TEC对雷达发射电磁波 产生的平均FRA，再利用FRA值反推得到电离层TEC。这种方法比仪器测量 法的测量精度低，但是能实时反应雷达回波信号穿过电离层时刻电离层对电磁 波的影响。
在另一个实施例中，可以通过将雷达原始回波数据按频率分成N个子带数 据，再利用每个子带数据分别估计子带TEC，将得到的N个子带TEC的统计 平均值作为雷达成像时电离层中的电子总量TEC估计结果，这种方法会提高 TEC估计精度。
在S104中，根据TEC确定雷达原始回波数据对应的法拉第旋转角数据； 法拉第旋转角数据随频率非线性变化。
在本实施例中，根据该雷达原始回波数据估算的电离层TEC可以确定该雷 达原始回波数据对应的法拉第旋转角度，由于本申请考虑了PD对雷达原始回 波数据的影响，因此，法拉第旋转角数据时随着频率非线性变化的。
在一个实施例中，需要获取SAR***成像时的磁信夹角、下视角及地磁场 的磁感应强度，并根据该雷达原始回波数据的TEC、SAR***成像时的磁信夹 角、下视角及地磁场的磁感应强度确定该法拉第旋转角数据。其中，磁信夹角 为磁场与雷达回波信号之间的夹角。
在实施时，可以建立该法拉第旋转角数据的表达式(1-1)，该法拉第旋转 角数据的表达式与TEC、地磁场的磁感应强度、磁信夹角、下视角有关：
Figure BDA0002497175870000041
其中，Ω(f)为法拉第旋转角数据，f为信号频率，B0为地磁场的磁感应强 度，α为下视角，θ为磁信夹角。
在S106中，根据法拉第旋转角数据确定PD的影响量化值。
在实施时，可以先获取理想状态下的发射信号，根据该理想状态下的发射 信号得到理想状态下的回波信号，在该理想状态下的回波信号的每个频点中加 入对应的法拉第旋转角，得到预测的受到PD影响的雷达回波数据，进而可以 得到对应的受到PD影响的点目标图像；再根据上述理想状态下的回波信号得 到理想的点目标图像；通过比对该受到PD影响的点目标图像及该理想的点目 标图像，即可得到PD的影响量化值。
在S108中，若影响量化值大于或等于预设阈值，则根据法拉第旋转角数据 对雷达原始回波数据进行校正。
其中，经申请人研究发现，当某些磁信夹角、下视角、与不同的TEC组合 时，PD和FRA对雷达原始回波数据的影响较小，即对SAR成像的影响很小， 此时不需要校正FRA影响。根据SAR***设计需求，当影响量化值小于该预 设阈值时，表示此时雷达原始回波数据中受到PD的影响较小，可以不对该雷 达原始回波数据进行校正。当影响量化值大于或等于预设阈值，则根据法拉第 旋转角数据对雷达原始回波数据进行校正。
在本申请实施例提供的雷达数据处理方法中，根据雷达原始回波数据获取 电离层中的电子总量TEC；根据TEC确定雷达原始回波数据对应的法拉第旋转 角数据；法拉第旋转角数据随频率非线性变化；根据法拉第旋转角数据确定PD 的影响量化值；若影响量化值大于或等于预设阈值，则根据法拉第旋转角数据 对雷达原始回波数据进行校正。根据本申请实施例提供的雷达数据处理方法， 由于在估计法拉第旋转角数据时考虑了PD对雷达原始回波数据的影响，将该 法拉第旋转角数据定义为随频率非线性变化的数据，并非一个固定值，在进行 后续校正时，可以根据该法拉第旋转角数据对雷达原始回波数据中PD影响进行校正；由于在对雷达原始回波数据进行校正之前，先根据预设阈值对PD的 影响量化值进行判断，可以保证只在PD影响较大时对雷达原始回波数据进行 校正，避免了PD影响较小时对雷达原始回波数据进行反向校正的问题。
在一些实施例中，参见图2，图2是本发明实施例提供的方法的一个可选 的流程示意图，图1示出的S102可以通过S1022至S1028实现，将结合各步 骤进行说明。
在S1022中，根据雷达原始回波数据的总带宽及子带数N，将雷达原始回 波数据划分为多个子带数据。
在本实施例中，先确定该雷达原始回波数据的总带宽，并获取需要划分的 子带数N。其中，该子带数N可以预设在***中，也可以通过交互设备接收操 作人员的实时输入。该子带数N一般为大于1的整数，随着子带数N的增加， 后续的计算量也随之增加，但会提高TEC的估计精度。
在S1024中，分别确定各子带数据对应的FRA。
在本实施例中，为了提高TEC的估计精度，会先将该原始雷达数据进行成 像得到原始图像数据，并在原始图像数据中确定目标均匀分布的区域，后续计 算中，只对该区域内的图像数据进行计算。其中，可以通过现有的图像识别算 法获取该原始图像数据中目标均匀分布的区域，也可以接收操作人员在该原始 图像数据上的框选指令确定该目标均匀分布的区域。
在一个实施例中，可以通过下式(1-2)确定各子带数据对应的FRA：
Figure BDA0002497175870000061
其中，Ωi为第i个子带数据对应的FRA，(·)*表示共轭，Z12、Z21可以通 过下式得到公式(1-3)：
Figure BDA0002497175870000062
其中，j为
Figure BDA0002497175870000063
为第i个子带数据对应的hh极化方式下的图像数据，
Figure BDA0002497175870000064
为第i个子带数据对应的hv极化方式下的图像数据，
Figure BDA0002497175870000065
为第i个子带数据对应的vh 极化方式下的图像数据，
Figure BDA0002497175870000066
为第i个子带数据对应的vv极化方式下的图像数据。
在S1026中，根据各子带数据对应的FRA，反演得到各子带数据对应的子 带TEC。
在本实施例中，各子带数据对应的子带TEC与各子带数据对应的FRA、各 子带数据对的中心频率、下视角和地磁场的磁感应强度相关，在一个实施例中， 各子带数据对应的子带TEC可以由下式(1-4)确定：
Figure BDA0002497175870000071
其中，TECi为第i个子带数据对应的子带TEC，α为下视角，B0为地磁场 的磁感应强度，fi为第i个子带数据对应的中心频率，θ为磁信夹角。其中， B0＝(4e-5)～(5e-5)，单位特斯拉。
在S1028中，将各子带TEC的平均值作为TEC。
在本实施例中，可以通过下式(1-5)确定TEC：
Figure BDA0002497175870000072
其中，UTEC为雷达成像时电离层中的电子总量TEC，N为子带数。
在本申请实施例提供的雷达数据处理方法中，根据雷达原始回波数据的总 带宽及子带数N，将雷达原始回波数据划分为多个子带数据；分别确定各子带 数据对应的FRA；根据各子带数据对应的FRA，反演得到各子带数据对应的子 带TEC；将各子带TEC的平均值作为雷达成像时电离层中的电子总量TEC。 根据本申请实施例提供的雷达数据处理方法，由于将雷达原始回波数据划分为 N个子带数据，并分别计算各子带数据对应的FRA，可以准确的对各个频率下 的FRA，得到的多个子带数据对应的FRA可以准确的反应PD影响下的FRA的非线性变化情况；由于分别计算各子带数据对应的TEC，进一步的提升了对 雷达原始回波数据中TEC预测的准确度。
在一些实施例中，基于图1，该影响量化值包括多种极化方式中每一极化 方式对应的分辨率展宽比。图1中的S108具体为：若多种极化方式中每一极化 方式对应的分辨率展宽比中存在至少一个分辨率展宽比大于或等于预设阈值， 根据法拉第旋转角数据对雷达原始回波数据进行校正。
在本实施例中，该多种极化方式可以包括四种极化方式：HH、VV、HV、 VH。雷达发射的能量脉冲的电场矢量，可以在垂直或水平面内被偏振。无论哪 个波长，雷达信号可以发射水平(H)或者垂直(V)电场矢量，接收水平(H) 或者垂直(V)或者两者的回波信号。其中，HH、VV为同极化，HV、VH为 交叉极化。
由于每种极化方式对应的雷达回波信号受到PD影响的大小不同，因此， 对应的，各极化方式对应的分辨率展宽比也会存在差异。也就是说，对同一个 雷达原始回波数据，在同极化下的雷达回波数据对应的分辨率展宽比，与在交 叉极化下的雷达回波数据对应的分辨率展宽比可能均大于该预设阈值，也可能 均小于或等于该预设阈值，还有可能在同极化下的雷达回波数据对应的分辨率 展宽比大于该预设阈值，且在交叉极化下的雷达回波数据对应的分辨率展宽比 小于或等于该预设阈值，还有可能在同极化下的雷达回波数据对应的分辨率展 宽比小于或等于该预设阈值，且在交叉极化下的雷达回波数据对应的分辨率展 宽比大于该预设阈值。因此，需要对每一极化方式下的雷达回波数据进行分别计算，以得到各极化方式下对应的分辨率展宽比。在实施中，该预设阈值可以 设置为1.03。
在本实施例中，只要存在一种极化方式下的分辨率展宽比大于或等于该预 设阈值时，就对根据法拉第旋转角数据对雷达原始回波数据进行校正。
在一些实施例中，参见图3，图3是本发明实施例提供的方法的一个可选 的流程示意图，基于图1及上述实施例，图1示出的S106可以通过S1062至 S1068实现，将结合各步骤进行说明。
在S1062中，根据法拉第旋转角数据确定法拉第旋转矩阵。
在本实施例中，已经确定法拉第旋转角数据为与TEC、地磁场的磁感应强 度、磁信夹角、下视角和频率有关表达式(1-6)：
Figure BDA0002497175870000081
根据下式(1-7)确定所述法拉第旋转矩阵：
Figure BDA0002497175870000082
其中，Ff为法拉第旋转矩阵。
在S1064中，根据特定的理想回波信号及法拉第旋转矩阵，生成受PD影 响的点目标图像数据。
在本实施例中，需要获取预先设置好的理想发射信号，该理想发射信号可 以由下式(8)确定：
Figure BDA0002497175870000083
其中，f0为中心频率，k为线性调频斜率，
Figure BDA0002497175870000084
T为脉冲重复 周期，rect(·)为宽度和幅度都为1的矩形窗函数。
任意时刻t的信号频率由下式(1-9)确定：
Figure BDA0002497175870000091
在获取到该理想发射信号后，根据特定的时延τ及该特定的理想发射信号即 可得到该特定的理想回波信号。
在本实施例中，根据该法拉第旋转矩阵，可以在该理想回波信号中每一个 频点上加入对应频点的法拉第旋转角，此时，可以预估受到PD影响的理想回 波信号，在对其进行成像处理后，可以得到受PD影响的点目标图像数据。
在S1066中，根据理想回波信号，生成理想的点目标图像数据。
在本实施例中，直接对该理想回波信号进行成像处理，生成理想的点目标 图像数据，作为对照图像。
在S1068中，将受PD影响的点目标图像数据及理想点的目标图像数据进 行比较，得到影响量化值。
在本实施例中，通过比较受PD影响的点目标图像数据与理想点的目标图 像数据，可以得到PD对点目标图像数据的影响量化值，其中，若该影响量化 值越大，表示受PD影响的点目标图像数据与理想点的目标图像数据之间的差 异越大，即PD对点目标图像数据的影响越高；若该影响量化值越小，表示受 PD影响的点目标图像数据与理想点的目标图像数据之间的差异越小，即PD对 点目标图像数据的影响越小。
在本申请实施例提供的雷达数据处理方法中，根据法拉第旋转角数据确定 法拉第旋转矩阵；根据特定的理想回波信号及法拉第旋转矩阵，生成受PD影 响的点目标图像数据；根据理想回波信号，生成理想的点目标图像数据；将受 PD影响的点目标图像数据及理想点的目标图像数据进行比较，得到影响量化值。 由于将根据理想回波信号生成的理想的点目标图像数据作为基本判断图像，去 衡量根据特定的理想回波信号及法拉第旋转矩阵生成的受PD影响的点目标图 像数据与自身的差异，得到影响量化值，可以准确的反应出该法拉第旋转矩阵 对点目标图像生成时的影响大小，即判断PD影响大小，进而为后续的步骤提 供判断依据，间接的提升了准确性。
在一些实施例中，参见图4，图4是本发明实施例提供的方法的一个可选 的流程示意图，基于图3及上述实施例，图3示出的S1064可以通过S10642 至S10646实现，将结合各步骤进行说明。
在S10642中，根据特定的理想回波信号建立理想回波散射矩阵；理想回波 信号包括各极化方式对应的理想回波信号。
在本实施例中，该理想回波信号包括多种极化方式对应的理想回波信号， 其中，可以包括HH、VV、HV、VH四种极化方式下的理想回波信号，根据各 极化方式下的理想回波信号可以建立对应的理想回波散射矩阵，如下式(1-10)。
Figure BDA0002497175870000101
其中，Ss为理想回波散射矩阵，
Figure BDA0002497175870000102
为HH极化方式对应的理想回波信号，
Figure BDA0002497175870000103
为HV极化方式对应的理想回波信号，
Figure BDA0002497175870000104
为VH极化方式对应的理想回波 信号，
Figure BDA0002497175870000105
为VV极化方式对应的理想回波信号。
在S10644中，根据理想回波散射矩阵及法拉第旋转矩阵，确定受到PD影 响的散射矩阵。
在本实施例中，由于在雷达发射过程和雷达接收过程中，均会受到PD影 响，因此，需要对计算两次法拉第旋转矩阵的影响。因此，受到PD影响的散 射矩阵可以由下式(1-11)确定：
Sf＝FfSsFf (1-11)；
其中，Sf为受到PD影响的散射矩阵，Ff为法拉第旋转矩阵，Ss为理想回 波散射矩阵。
在S10646中，对受到PD影响的散射矩阵进行距离向快速傅里叶变换，得 到受PD影响的点目标图像数据。
在本申请实施例提供的雷达数据处理方法中，根据特定的理想回波信号建 立理想回波散射矩阵；理想回波信号包括各极化方式对应的理想回波信号；根 据理想回波散射矩阵及法拉第旋转矩阵，确定受到PD影响的散射矩阵；对受 到PD影响的散射矩阵进行距离向快速傅里叶变换，得到受PD影响的点目标图 像数据。由于根据各极化方式对应的理想回波信号建立理想回波散射矩阵，并 根据法拉第旋转矩阵可以快速得到对应的受到PD影响的散射矩阵，进而获得 受PD影响的点目标图像数据，为后续与理想的点目标图像数据得比较提供数 据基础。
在一些实施例中，参见图5，图5是本发明实施例提供的方法的一个可选 的流程示意图，基于图3及上述实施例，图3示出的S1068可以通过S10682 至S10686实现，将结合各步骤进行说明。
在S10682中，获取受PD影响的点目标图像数据对应的第一3dB宽度；第 一3dB宽度包括各极化方式对应的第一3dB宽度。
在实施时，可以通过获取受PD影响的点目标图像数据及理想的点目标图 像数据中的图像特征进行比较，在本实施例中，通过获取点目标图像数据中的 3dB宽度作为图像特征。
其中，得到的受PD影响的点目标图像数据在不同极化方式下，对应的3dB 宽度也会存在差异，因此，需要先分别获取上述各极化方式对应的点目标图像 数据，并进一步分别获取上述各极化方式对应的第一3dB宽度。
在S10684中，获取理想的点目标图像数据对应的第二3dB宽度；第二3dB 宽度包括各极化方式对应的第二3dB宽度。
在本实施例中，和上述描述相同，理想的点目标图像数据也包括不同极化 方式对应的点目标图像数据，分别获取各极化方式对应的第二3dB宽度。
在S10686中，根据各极化方式对应的第一3dB宽度及第二3dB宽度，分 别确定各极化方式对应的分辨率展宽比。
在实施例中，在得到各极化方式对应的第一3dB宽度及第二3dB宽度之后， 可以分别确定各极化方式对应的分辨率展宽比，其中，对于每一极化方式，其 对应的分辨率展宽比为第一3dB宽度与第二3dB宽度的比值。在实施时，将各 极化方式对应的分辨率展宽比作为该影响量化值。
在本申请实施例提供的雷达数据处理方法中，获取受PD影响的点目标图 像数据对应的第一3dB宽度；第一3dB宽度包括各极化方式对应的第一3dB宽 度；获取理想的点目标图像数据对应的第二3dB宽度；第二3dB宽度包括各极 化方式对应的第二3dB宽度；根据各极化方式对应的第一3dB宽度及第二3dB 宽度，分别确定各极化方式对应的分辨率展宽比。由于获取3dB宽度作为各点 目标图像数据的图像特征，通过获取受PD影响的点目标图像数据及理想的点 目标图像数据中的3dB宽度之间的比值作为影响量化值，可以准确的量化两个 点目标图像数据之间的差异，因此可以准确的反应出该法拉第旋转矩阵对点目 标图像生成时的影响大小，即判断PD影响大小，进而为后续的步骤提供判断 依据，间接的提升了校正的准确性。
在一些实施例中，参见图6，图6是本发明实施例提供的方法的一个可选 的流程示意图，基于图1及上述实施例，图1中根据所述法拉第旋转角数据对 所述雷达原始回波数据进行校正的步骤，可以通过S1082至S10810实现，将 结合各步骤进行说明。
在S1082中，获取发射失真矩阵及接收失真矩阵。
在本实施例中，由于在SAR***的发射和接收过程中，会存在一定程度上 的失真，因此，需要先获取发射时的发射失真矩阵R及接收时的接收失真矩阵 T。在实施时，可以通过SAR***外定标得到该失真矩阵。其中，若R、T未 知，则可以将R、T均预设为
Figure BDA0002497175870000121
在S1084中，根据发射失真矩阵及法拉第旋转矩阵，确定频域发射校正矩 阵。
在本实施例中，可以先获取时域发射校正矩阵CR的表达式(1-12)。
Figure BDA0002497175870000122
其中，Ff为法拉第旋转矩阵，R为发射失真矩阵。
之后，对该时域发射校正矩阵CR进行时频变换，得到频域发射校正矩阵
Figure BDA0002497175870000123
在S1086中，根据接收失真矩阵及法拉第旋转矩阵，确定频域接收校正矩 阵。
在本实施例中，可以先获取时域接收校正矩阵CT的表达式(1-13)。
Figure BDA0002497175870000124
其中，Ff为法拉第旋转矩阵，T为接收失真矩阵。
之后，对该时域发射校正矩阵CR进行时频变换，得到频域接收校正矩阵
Figure BDA0002497175870000125
在S1088中，根据频域发射校正矩阵及频域接收校正矩阵，对频域的雷达 原始回波数据进行校正，得到距离向频域数据。
在本实施例中，对雷达原始回波数据进行时频变化，可以得到对应的频域 的雷达原始回波数据，其中，对该雷达原始回波数据进行距离向快速傅里叶变 换，可以得到频域的雷达原始回波数据。
根据时域相乘对应频域卷积的原理，可以对频域的雷达原始回波数据进行 校正，得到距离向频域数据，即通过下式(1-14)确定。
Figure BDA0002497175870000131
其中，
Figure BDA0002497175870000132
为频域发射校正矩阵，MFFT为频域的雷达原始回波数据，
Figure BDA0002497175870000134
为 频域接收校正矩阵，
Figure BDA0002497175870000135
为卷积。
在S10810中，对距离向频域数据进行距离向反向快速傅里叶变换，得到时 域散射矩阵；时域散射矩阵包括各极化方式对应的校正回波数据。
在本实施例中，该距离向频域数据已经为校正后的回波数据，但需要对其 进行频域到时域的变换，在实施时，可以采用距离向反向快速傅里叶变换对该 距离向频域数据进行时频变换，得到时域散射矩阵。
在该时域散射矩阵中，同样包括上述多种极化方式对应的校正回波数据， 也就是说，该时域散射矩阵包括HH极化方式对应的校正回波数据、HV极化 方式对应的校正回波数据、VH极化方式对应的校正回波数据、VV极化方式对 应的校正回波数据。
在本申请实施例提供的雷达数据处理方法中，由于获取了发射失真矩阵及 接收失真矩阵，并根据法拉第旋转矩阵建立了对应的频域发射校正矩阵及频域 接收校正矩阵，由于考虑到了SAR***在发射和接收时的失真情况，使得本实 施例中的校正过程更加准确，得到的时域散射矩阵与理想的回波信号更加接近， 提高了校正效果。
在一些实施例中，参见图7，图7是本发明实施例提供的方法的一个可选 的流程示意图，基于图6及上述实施例，在步骤在S10810之后，还可以执行 S10812至S10816。
在S10812中，将分辨率展宽比小于预设阈值的至少一种极化方式，确定为 无PD极化方式。
在本实施例中，根据上述实施例中的描述，每种极化方式对应的雷达回波 信号受到PD影响的大小不同，因此，对应的，各极化方式对应的分辨率展宽 比也会存在差异，在完成校正过程，得到各极化方式对应的校正回波数据后， 还需要判断各极化方式中分辨率展宽比是否存在小于预设阈值的极化方式，并 将小于预设阈值的至少一种极化方式确定为无PD极化方式。其中，该无PD极 化方式表示PD对该极化方式下的雷达回波数据影响较小，无需对其校正处理。
在S10814中，根据雷达原始回波数据，确定各无PD极化方式对应的原始 回波数据。
在本实施例中，若存在至少一种无PD极化方式，根据雷达原始回波数据， 提取该雷达原始回波数据中各无PD极化方式对应的原始回波数据。
在S10816中，将时域散射矩阵中各无PD极化方式对应的校正回波数据， 替换为各无PD极化方式对应原始回波数据。
在本实施例中，由于该无PD极化方式中PD对该极化方式下的雷达回波数 据影响较小，对其进行校正时会出现反向校正的现象，因此该无PD极化方式 对应的校正回波数据可能与理想的回波数据相差较大，需要将时域散射矩阵中 各无PD极化方式对应的校正回波数据，替换为各无PD极化方式对应原始回波 数据，此时得到的时域散射矩阵为更加接近理想的回波数据。
在本申请实施例提供的雷达数据处理方法中，由于将时域散射矩阵中各无 PD极化方式对应的校正回波数据，替换为各无PD极化方式对应原始回波数据， 可以避免各极化方式中对PD影响较小的无PD极化方式的回波数据的反向校正， 进一步提升了校正准确性和校正效果。
下面，将说明本发明实施例在一个实际的应用场景中的示例性应用。
本发明实施例中，利用极化合成孔径雷达获取的全极化图像数据，FAR影 响校正。利用极化SAR数据估计电离层TEC含量。评估电离层对SAR成像的 分辨率展宽比是否可忽略。当FRA及PD对SAR成像质量的影响不可忽略时， 对SAR测量数据进行FRA影响的校正处理。
图8是本发明实施例提供的方法的一个可选的流程示意图，在本实施例中， 提供了一种法拉第旋转极化色散对星载低频全极化SAR成像影响的校正方法， 该方法主要包括：
在S801中，通过PolSAR回波数据计算成像时的TEC含量；
在S802中，计算上述TEC含量下，磁信夹角、及下视角下的FRA(此FRA 为受PD影响的FRA)；
在S803中，仿真此FRA及PD下的SAR点目标图像；
在S804中，评估上述FRA及PD对SAR成像的影响程度；
在S805中，当FRA及PD对SAR成像质量的影响不可忽略时，对SAR 测量数据进行FRA影响的校正处理。
以下分别对各步骤进行具体描述：
对于上述S801：通过PolSAR回波数据计算成像时的TEC含量。
目前，电离层TEC测量方法主要有利用专用仪器设备直接测量和从SAR 回波数据中估计两类。其中专用仪器设备测量精度高，但是时效性差，有测量 延时，不能实时反应雷达信号穿过电离层时刻电离层对电磁波的影响。从SAR 回波数据中估计电离层TEC的主要方法有：Bickel and Bates方法、Freeman方 法、Qi and Jin方法、Chen and Quegan方法。这些方法可以从极化SAR数据中 估计电离层TEC对雷达发射电磁波产生的平均FRA，再利用FRA值反推得到 电离层TEC。这种方法比仪器测量法的测量精度低，但是能实时反应雷达信号 穿过电离层时刻电离层对电磁波的影响。从SAR回波数据中估计电离层TEC 的方法可以得到电离层TEC对雷达发射电磁波产生的平均FAR，前期研究结果 表明，电离层TEC为50TECU时，P波段中心频率550MHz，宽带250MHz时， FRA变化范围可以达到90°，因此，常规的直接从SAR回波数据中估计电离 层TEC的误差很大。
因此，本项目考虑将SAR回波数据按频率分成N个子带，再利用每个子 带分别估计电离层TEC，将得到的N个电离层TEC值的统计平均值作为电离 层TEC估计结果，这种方法会提高TEC估计精度。
在一些实施例中，上述S801可以包括：
步骤(11)将SAR原始回波数据按频率划分为N(N大于1)个子带。每 个子带的频率带宽为总带宽/N。N越大则后续计算量越大，TEC含量估计精度 将越高。
步骤(12)在图像中选择分布目标的均匀场景区域，计算每个子带对应频 点的FRA值Ωi。
其中，步骤(12)可采用下式(2-1)计算：
Figure BDA0002497175870000151
其中，Z12、Z21由下式(2-2)计算。
Figure BDA0002497175870000152
这里，
Figure BDA0002497175870000153
是HH、HV、VH、VV四种极化方式下的图像数据。
步骤(13)得到FRA后，可由下式(2-3)反演每个子带TEC值。
Figure BDA0002497175870000161
其中f为频率，α为雷达下视角。B0是地磁场的磁感应强度， B0＝(4e-5)～(5e-5)，单位特斯拉。
步骤(14)平均所有频点的TEC，作为成像区域TEC含量估计值UTEC。
设每个子带估计得到的TEC值为TECi，则UTEC的表达式(2-4)为：
Figure BDA0002497175870000163
在一些实施例中，对于上述S802：计算上述TEC含量下，磁信夹角、及 下视角下的FRA(此FRA为受PD影响的FRA)。上述S802可以包括：
步骤(21)确定***成像时的磁信夹角θ、下视角α。
步骤(22)根据公式计算TEC含量下，磁信夹角θ、及下视角α下的FRA， FRA是随频率非线性变化的。
Ω(f)为与频率f相关的FRA。
Figure BDA0002497175870000164
法拉第旋转角与路径TEC含量、信号频率f、地磁场的磁感应强度B0以及 磁场与雷达信号夹角θ、雷达下视角α有关。
TEC表示传播路径上的电子总量。
在一些实施例中，对于上述S803：仿真此FRA及PD下的SAR点目标图 像。上述S803可以包括：
步骤(31)仿真理想情况下的SAR点目标距离向线性调频回波信号。s(t) 是发射的线性调频信号，
Figure BDA0002497175870000165
其中，f0为中心频率，k为线性调频斜率，任意时刻t的信号频率为：
Figure BDA0002497175870000166
步骤(32)在SAR回波信号中每个频率点上加入对应频点的FRA。对全 极化SAR***，回波信号包含如矩阵S中HH、HV、VH、VV四种极化的回 波。
Figure BDA0002497175870000171
其中，
Sf＝FfSsFf (2-9)；
Figure BDA0002497175870000172
步骤(33)进行距离向脉冲压缩(FFT)，得到受PD影响的点目标图像数 据。
请参照图9a至图9k，其示出了在下视角为30度，磁信夹角为60度时， 此FRA及PD下的SAR点目标图像示意图，其中，图9a为同极化下、0TECU 的SAR点目标图像示意图，图9b为同极化下、10TECU的SAR点目标图像示 意图，图9c为同极化下、20TECU的SAR点目标图像示意图，图9d为同极化 下、30TECU的SAR点目标图像示意图，图9e为同极化下、40TECU的SAR点目标图像示意图，图9f为同极化下、50TECU的SAR点目标图像示意图， 图9g为交叉极化下、10TECU的SAR点目标图像示意图，图9h为交叉极化下、 20TECU的SAR点目标图像示意图，图9i为交叉极化下、30TECU的SAR点 目标图像示意图，图9j为交叉极化下、40TECU的SAR点目标图像示意图，图 9k为交叉极化下、50TECU的SAR点目标图像示意图。
在一些实施例中，对于上述S804：评估上述FRA及PD对SAR成像的影 响程度。上述S804可以包括：
步骤(41)计算S803中四种极化的点目标的3dB宽度。
步骤(42)仿真理想情况下的SAR点目标距离向线性调频回波信号，并进 行距离向脉冲压缩(FFT)，得到理想的点目标图像。
步骤(43)计算理想点目标的3dB宽度。
步骤(44)计算受FRA影响的点目标的分辨率展宽比。
在一些实施例中，对于上述S805：当FRA及PD对SAR成像质量的影响 不可忽略时，对SAR测量数据进行FRA影响的校正处理。上述S805可以包括：
步骤(51)研究发现，当某些磁信夹角、雷达信号下视角、与TECU组合 时，FRA对SAR成像的影响很小，此时不需要校正FRA影响。通常，根据SAR ***设计需求，当分辨率展宽比优于某值(根据经验，本发明设计为1.03)时， 可以不校正FRA的影响。
表1给出了中心频率550MHz，***信号带宽300M时，入射角30°时， 不同TECU和分辨率展宽比优于1.03所对应的磁信夹角。也就是说，当入射角 为30°时，当TEC含量和磁信夹角在表中范围时，不校正FRA的影响。
表1线极化工作方式，分辨率展宽比优于1.03所对应的磁信夹角
Figure BDA0002497175870000181
步骤(52)当点目标图像的分辨率展宽比大于1.03时，校正FRA对SAR 测量数据的影响。
全极化SAR***存在发射和接收失真时的目标的后向散射矩阵如下：
Figure BDA0002497175870000182
其中，R、T分别为发射和接收失真矩阵，可由***外定标得到。当R、T 未知时，可假设
Figure BDA0002497175870000183
分别设发射和接收校正矩阵如下：
Figure BDA0002497175870000191
Figure BDA0002497175870000192
对发射和接收校正矩阵进行傅里叶变换(FFT)处理，变换到频域。
Figure BDA0002497175870000193
Figure BDA0002497175870000194
根据时域相乘对应频域卷积的原理，得到FRA影响校正后的距离向频域散 射矩阵。
Figure BDA0002497175870000195
将校正后的距离向频域数据进行距离向IFFT，得到FRA影响校正后的时 域散射矩阵。
S＝IFFT(SFFT) (2-17)；
式(2-17)即为FRA影响校正后的全极化四通道(HH、HV、VH、VV) 时域后向散射回波数据。
图10是本发明实施例提供的方法的一个可选的流程示意图，本实施例中， 可得校正FRA对SAR成像的影响流程包括：
在S1001中，由上式(2-5)得到的TEC值计算所有频点对应的FRA值。
在S1002中，由上式(2-10)计算所有频点对应的FR矩阵。
在S1003中，由上式(2-12)和(2-13)得到发射校正矩阵和接收校正矩 阵。因为TEC对每个频点引起的法拉第旋转角不同，因而，需要在频域针对每 个频点补偿FRA。
在S1004中，将发射校正矩阵和接收校正矩阵变换到频域，如上式(2-14) 和(2-15)。
在S1005中，将SAR原始回波数据进行距离向FFT，将数据变换的距离频 域。
在S1006中，根据公式(2-16)，将发射校正矩阵、FR矩阵和接收校正矩 阵卷积，得到FRA影响校正后的距离向频域散射矩阵。
在S1007中，将校正后的距离向频域数据进行距离向IFFT，得到FRA影 响校正后的时域散射矩阵，即上式(2-17)。
本发明实施例还提供雷达回波数据处理装置，用于实施上述的雷达回波数 据处理方法。图11为本发明实施例提供的雷达回波数据处理装置的部分结构示 意图。如图11所示，所述雷达回波数据处理装置包括：第一获取模块11、第 一确定模块12、第二确定模块13和校正模块14；其中：
第一获取模块11，用于根据雷达原始回波数据获取电离层中的电子总量 TEC；第一确定模块12，用于根据所述TEC确定所述雷达原始回波数据对应的 法拉第旋转角数据；所述法拉第旋转角数据随频率非线性变化；第二确定模块 13，用于根据所述法拉第旋转角数据确定PD的影响量化值；校正模块14，用 于若所述影响量化值大于或等于预设阈值，则根据所述法拉第旋转角数据对所 述雷达原始回波数据进行校正。
在一些实施例中，所述第一获取模块11，还用于根据雷达原始回波数据的 总带宽及子带数N，将雷达原始回波数据划分为多个子带数据；分别确定各子 带数据对应的FRA；根据各子带数据对应的FRA，反演得到各子带数据对应的 子带TEC；将各子带TEC的平均值作为TEC。
在一些实施例中，所述第一确定模块12，还用于获取磁信夹角、下视角和 地磁感应强度；根据TEC、磁信夹角、下视角和地磁感应强度，确定法拉第旋 转角数据。
在一些实施例中，影响量化值包括多种极化方式中每一极化方式对应的分 辨率展宽比，所述校正模块14，还用于若多种极化方式中每一极化方式对应的 分辨率展宽比中存在至少一个分辨率展宽比大于或等于预设阈值，根据法拉第 旋转角数据对雷达原始回波数据进行校正。
在一些实施例中，所述第二确定模块13，还用于根据法拉第旋转角数据确 定法拉第旋转矩阵；根据特定的理想回波信号及法拉第旋转矩阵，生成受PD 影响的点目标图像数据；根据理想回波信号，生成理想的点目标图像数据；将 受PD影响的点目标图像数据及理想点的目标图像数据进行比较，得到影响量 化值。
在一些实施例中，所述第二确定模块13，还用于根据特定的理想回波信号 建立理想回波散射矩阵；理想回波信号包括各极化方式对应的理想回波信号； 根据理想回波散射矩阵及法拉第旋转矩阵，确定受到PD影响的散射矩阵；对 受到PD影响的散射矩阵进行距离向快速傅里叶变换，得到受PD影响的点目标 图像数据。
在一些实施例中，所述第二确定模块13，还用于获取受PD影响的点目标 图像数据对应的第一3dB宽度；第一3dB宽度包括各极化方式对应的第一3dB 宽度；获取理想的点目标图像数据对应的第二3dB宽度；第二3dB宽度包括各 极化方式对应的第二3dB宽度；根据各极化方式对应的第一3dB宽度及第二 3dB宽度，分别确定各极化方式对应的分辨率展宽比。
在一些实施例中，所述校正模块14，还用于获取发射失真矩阵及接收失真 矩阵；根据发射失真矩阵及法拉第旋转矩阵，确定频域发射校正矩阵；根据接 收失真矩阵及法拉第旋转矩阵，确定频域接收校正矩阵；根据频域发射校正矩 阵及频域接收校正矩阵，对频域的雷达原始回波数据进行校正，得到距离向频 域数据；对距离向频域数据进行距离向反向快速傅里叶变换，得到时域散射矩 阵；时域散射矩阵包括各极化方式对应的校正回波数据。
在一些实施例中，所述校正模块14，还用于将分辨率展宽比小于预设阈值 的至少一种极化方式，确定为无PD极化方式；根据雷达原始回波数据，确定 各无PD极化方式对应的原始回波数据；将时域散射矩阵中各无PD极化方式对 应的校正回波数据，替换为各无PD极化方式对应原始回波数据。
以上装置实施例的描述，与上述方法实施例的描述是类似的，具有同方法 实施例相似的有益效果。对于本申请装置实施例中未披露的技术细节，请参照 本申请方法实施例的描述而理解。
需要说明的是，本申请实施例中，如果以软件功能模块的形式实现上述的 雷达回波数据处理方法，并作为独立的产品销售或使用时，也可以存储在一个 计算机可读取存储介质中。基于这样的理解，本申请实施例的技术方案本质上 或者说对相关技术做出贡献的部分可以以软件产品的形式体现出来，该计算机 软件产品存储在一个存储介质中，包括若干指令用以使得一台雷达回波数据处 理设备(例如可以是计算机设备等)执行本申请各个实施例所述方法的全部或 部分。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ReadOnly Memory， ROM)、磁碟或者光盘等各种可以存储程序代码的介质。这样，本申请实施例不限制于任何特定的硬件和软件结合。
图12为本发明实施例提供的雷达回波数据处理设备的部分结构示意图。如 图12所示，雷达回波数据处理设备2包括：存储器21和处理器22，存储器21 与处理器22通过总线23连接；存储器21，用于存储可执行数据指令；处理器 22，用于执行所述存储器中存储的可执行指令时，实现上述的雷达回波数据处 理方法。
本发明实施例还提供一种计算机可读存储介质，存储有可执行指令，用于 引起处理器执行时，实现上述方法实施例中的雷达回波数据处理方法。
本领域内的技术人员应明白，本发明的实施例可提供为方法、***、或计 算机程序产品。因此，本发明可采用硬件实施例、软件实施例、或结合软件和 硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算 机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储 器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品 的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方 框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结 合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或 其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可 编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个 流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备 以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的 指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流 程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使 得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理， 从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程 或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范 围。

Claims (11)

1.一种雷达回波数据处理方法，其特征在于，包括：
根据雷达原始回波数据获取电离层中的电子总量TEC；
根据所述TEC确定所述雷达原始回波数据对应的法拉第旋转角数据；所述法拉第旋转角数据随频率非线性变化；
根据所述法拉第旋转角数据确定极化色散PD的影响量化值；
若所述影响量化值大于或等于预设阈值，则根据所述法拉第旋转角数据对所述雷达原始回波数据进行校正。
2.根据权利要求1所述的方法，其特征在于，所述根据雷达原始回波数据获取电离层中的电子总量TEC，包括：
根据所述雷达原始回波数据的总带宽及子带数N，将所述雷达原始回波数据划分为多个子带数据；
分别确定各所述子带数据对应的法拉第旋转角FRA；
根据各所述子带数据对应的FRA，反演得到各所述子带数据对应的子带TEC；
将各所述子带TEC的平均值作为所述TEC。
3.根据权利要求1所述的方法，其特征在于，所述根据所述TEC确定所述雷达原始回波数据对应的法拉第旋转角数据，包括：
获取磁信夹角、下视角和地磁感应强度；
根据所述TEC、所述磁信夹角、所述下视角和所述地磁感应强度，确定所述法拉第旋转角数据。
4.根据权利要求1所述的方法，其特征在于，所述影响量化值包括多种极化方式中每一极化方式对应的分辨率展宽比；
所述若所述影响量化值大于或等于预设阈值，则根据所述法拉第旋转角数据对所述雷达原始回波数据进行校正，包括：
若所述多种极化方式中每一极化方式对应的分辨率展宽比中存在至少一个分辨率展宽比大于或等于所述预设阈值，根据所述法拉第旋转角数据对所述雷达原始回波数据进行校正。
5.根据权利要求4所述的方法，其特征在于，所述根据所述法拉第旋转角数据确定影响量化值，包括：
根据所述法拉第旋转角数据确定法拉第旋转矩阵；
根据特定的理想回波信号及所述法拉第旋转矩阵，生成受PD影响的点目标图像数据；
根据所述理想回波信号，生成理想的点目标图像数据；
将所述受PD影响的点目标图像数据及所述理想点的目标图像数据进行比较，得到所述影响量化值。
6.根据权利要求5所述的方法，其特征在于，所述根据特定的理想回波信号及所述法拉第旋转矩阵，生成受PD影响的点目标图像数据，包括：
根据特定的理想回波信号建立理想回波散射矩阵；所述理想回波信号包括各所述极化方式对应的理想回波信号；
根据所述理想回波散射矩阵及所述法拉第旋转矩阵，确定受到PD影响的散射矩阵；
对所述受到PD影响的散射矩阵进行距离向快速傅里叶变换，得到受PD影响的点目标图像数据。
7.根据权利要求5所述的方法，其特征在于，所述将所述受PD影响的点目标图像数据及所述理想点的目标图像数据进行比较，得到所述影响量化值理想点的目标图像数据，包括：
获取所述受PD影响的点目标图像数据对应的第一3dB宽度；所述第一3dB宽度包括各所述极化方式对应的第一3dB宽度；
获取所述理想的点目标图像数据对应的第二3dB宽度；所述第二3dB宽度包括各所述极化方式对应的第二3dB宽度；
根据各所述极化方式对应的第一3dB宽度及第二3dB宽度，分别确定各所述极化方式对应的分辨率展宽比。
8.根据权利要求5所述的方法，其特征在于，所述根据所述法拉第旋转角数据对所述雷达原始回波数据进行校正，包括：
获取发射失真矩阵及接收失真矩阵；
根据所述发射失真矩阵及所述法拉第旋转矩阵，确定频域发射校正矩阵；
根据所述接收失真矩阵及所述法拉第旋转矩阵，确定频域接收校正矩阵；
根据所述频域发射校正矩阵及所述频域接收校正矩阵，对频域的雷达原始回波数据进行校正，得到距离向频域数据；
对所述距离向频域数据进行距离向反向快速傅里叶变换，得到时域散射矩阵；所述时域散射矩阵包括各所述极化方式对应的校正回波数据。
9.根据权利要求8所述的方法，其特征在于，所述方法还包括：
将分辨率展宽比小于所述预设阈值的至少一种极化方式，确定为无PD极化方式；
根据所述雷达原始回波数据，确定各所述无PD极化方式对应的原始回波数据；
将时域散射矩阵中各所述无PD极化方式对应的校正回波数据，替换为各所述无PD极化方式对应原始回波数据。
10.一种雷达回波数据处理装置，其特征在于，所述装置包括：
第一获取模块，用于根据雷达原始回波数据获取电离层中的电子总量TEC；
第一确定模块，用于根据所述TEC确定所述雷达原始回波数据对应的法拉第旋转角数据；所述法拉第旋转角数据随频率非线性变化；
第二确定模块，用于根据所述法拉第旋转角数据确定PD的影响量化值；
校正模块，用于若所述影响量化值大于或等于预设阈值，则根据所述法拉第旋转角数据对所述雷达原始回波数据进行校正。
11.一种雷达回波数据处理设备，其特征在于，包括：
存储器，用于存储可执行数据指令；
处理器，用于执行所述存储器中存储的可执行指令时，实现权利要求1至10任一项所述的雷达回波数据处理方法。
CN202010421677.8A 2020-05-18 2020-05-18 一种雷达回波数据处理方法及装置、设备 Active CN111679258B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010421677.8A CN111679258B (zh) 2020-05-18 2020-05-18 一种雷达回波数据处理方法及装置、设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010421677.8A CN111679258B (zh) 2020-05-18 2020-05-18 一种雷达回波数据处理方法及装置、设备

Publications (2)

Publication Number Publication Date
CN111679258A true CN111679258A (zh) 2020-09-18
CN111679258B CN111679258B (zh) 2023-09-15

Family

ID=72434028

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010421677.8A Active CN111679258B (zh) 2020-05-18 2020-05-18 一种雷达回波数据处理方法及装置、设备

Country Status (1)

Country Link
CN (1) CN111679258B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5313210A (en) * 1993-02-23 1994-05-17 Ball Corporation Polarimetric radar signal mapping process
US6163295A (en) * 1999-04-09 2000-12-19 Vsis, Inc. Ionospheric correction for single frequency GPS receivers using two satellites
US6914553B1 (en) * 2004-11-09 2005-07-05 Harris Corporation Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the Faraday rotation, and associated methods
US6919839B1 (en) * 2004-11-09 2005-07-19 Harris Corporation Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the group delay, and associated methods
CN101551450A (zh) * 2009-05-06 2009-10-07 北京航空航天大学 一种星载极化sar法拉第旋转效应校正平台的构建方法
CN103792521A (zh) * 2014-01-17 2014-05-14 西安空间无线电技术研究所 一种受Faraday旋转变化影响的雷达目标后向散射模拟方法
CA2957541A1 (en) * 2014-08-08 2016-02-11 Urthecast Corp. Apparatus and methods for quad-polarized synthetic aperture radar
CN106154268A (zh) * 2016-09-06 2016-11-23 合肥工业大学 一种基于子带分割的低波段isar电离层校正方法
CN107526079A (zh) * 2016-06-22 2017-12-29 中国解放军空军工程大学 一种基于l型三天线干涉处理的空间自旋目标宽带雷达三维成像方法
CN107544068A (zh) * 2017-07-14 2018-01-05 电子科技大学 一种基于频域bp的图像域宽带合成方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5313210A (en) * 1993-02-23 1994-05-17 Ball Corporation Polarimetric radar signal mapping process
US6163295A (en) * 1999-04-09 2000-12-19 Vsis, Inc. Ionospheric correction for single frequency GPS receivers using two satellites
US6914553B1 (en) * 2004-11-09 2005-07-05 Harris Corporation Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the Faraday rotation, and associated methods
US6919839B1 (en) * 2004-11-09 2005-07-19 Harris Corporation Synthetic aperture radar (SAR) compensating for ionospheric distortion based upon measurement of the group delay, and associated methods
CN101551450A (zh) * 2009-05-06 2009-10-07 北京航空航天大学 一种星载极化sar法拉第旋转效应校正平台的构建方法
CN103792521A (zh) * 2014-01-17 2014-05-14 西安空间无线电技术研究所 一种受Faraday旋转变化影响的雷达目标后向散射模拟方法
CA2957541A1 (en) * 2014-08-08 2016-02-11 Urthecast Corp. Apparatus and methods for quad-polarized synthetic aperture radar
CN107526079A (zh) * 2016-06-22 2017-12-29 中国解放军空军工程大学 一种基于l型三天线干涉处理的空间自旋目标宽带雷达三维成像方法
CN106154268A (zh) * 2016-09-06 2016-11-23 合肥工业大学 一种基于子带分割的低波段isar电离层校正方法
CN107544068A (zh) * 2017-07-14 2018-01-05 电子科技大学 一种基于频域bp的图像域宽带合成方法

Also Published As

Publication number Publication date
CN111679258B (zh) 2023-09-15

Similar Documents

Publication Publication Date Title
Garren et al. Full-polarization matched-illumination for target detection and identification
Ainsworth et al. Orientation angle preserving a posteriori polarimetric SAR calibration
Mohammadzadeh et al. Maximum entropy-based interference-plus-noise covariance matrix reconstruction for robust adaptive beamforming
Laviada et al. Phaseless synthetic aperture radar with efficient sampling for broadband near-field imaging: Theory and validation
US20090315762A1 (en) Adaptive specific differential phase in dual-polarization radar
Zhang et al. Transient nearfield acoustic holography based on an interpolated time-domain equivalent source method
CN108398669B (zh) 一种基于无需预延迟处理的空时宽带自适应单脉冲测角方法
JP2001188084A (ja) 宇宙機器搭載用レーダ
CN109725297B (zh) 一种基于回波极化度的有源转发式干扰识别方法
Ahmad et al. A beamforming approach to stepped-frequency synthetic aperture through-the-wall radar imaging
CN110109077B (zh) 基于时间反演的mimo雷达相干角度估计方法
Bekasiewicz et al. Time-gating method with automatic calibration for accurate measurements of electrically small antenna radiation patterns in non-anechoic environments
Foster et al. Intrinsic instrumental polarization and high-precision pulsar timing
JP2010085164A (ja) 画像レーダ装置
CN108181624B (zh) 一种差分计算成像装置及方法
CN111679258B (zh) 一种雷达回波数据处理方法及装置、设备
CN108931761A (zh) 基于星载合成孔径雷达的电离层定标方法及***
Gu et al. Two-dimensional sparse synthetic aperture radar imaging method with stepped-frequency waveform
US11269052B2 (en) Signal processing method
Yuan et al. Orbital-angular-momentum-based electromagnetic vortex imaging by least-squares method
CN116500612A (zh) 一种极化信息驱动的微波关联成像方法
You et al. A novel DOA estimation algorithm of cyclostationary signal based on UCA in impulsive noise
US11194014B1 (en) System, method and apparatus for recovering polarization radar data
Case et al. Correcting mutual coupling and poor isolation for real-time 2-D microwave imaging systems
CN113985407B (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