CN115079110A - 火星车雷达信号的校正方法、装置、设备及介质 - Google Patents

火星车雷达信号的校正方法、装置、设备及介质 Download PDF

Info

Publication number
CN115079110A
CN115079110A CN202210776198.7A CN202210776198A CN115079110A CN 115079110 A CN115079110 A CN 115079110A CN 202210776198 A CN202210776198 A CN 202210776198A CN 115079110 A CN115079110 A CN 115079110A
Authority
CN
China
Prior art keywords
time
radar data
data
corrected
phase
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
CN202210776198.7A
Other languages
English (en)
Other versions
CN115079110B (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.)
National Astronomical Observatories of CAS
Original Assignee
National Astronomical Observatories 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 National Astronomical Observatories of CAS filed Critical National Astronomical Observatories of CAS
Priority to CN202210776198.7A priority Critical patent/CN115079110B/zh
Publication of CN115079110A publication Critical patent/CN115079110A/zh
Application granted granted Critical
Publication of CN115079110B publication Critical patent/CN115079110B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/40Means for monitoring or calibrating

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

本公开提供一种火星车雷达信号的校正方法、装置、设备介质,方法包括:对原始雷达数据中的每一道数据在时间域上进行插值,并对插值后的信号取绝对值,得到多道高时间分辨率雷达数据;从多道高时间分辨率雷达数据任意确定一道高时间分辨率雷达数据作为第一参考数据;将其他多道高时间分辨率雷达数据与第一参考数据进行时间域交叉相关分析,得到每一道高时间分辨率雷达数据对应的时间域相关曲线;利用每一道高时间分辨率雷达数据的时间域相关曲线的峰值对应的时间偏移对每一道高时间分辨率雷达数据的时间进行校正,并对时间校正后的信号进行相位域分析实现相位校正。该方法解决了火星车下传雷达信号出现跳变和断层等不连续性的问题。

Description

火星车雷达信号的校正方法、装置、设备及介质
技术领域
本公开涉及数据处理技术领域,尤其涉及一种火星车雷达信号的校正方法、装置、设备及介质。
背景技术
火星车次表层探测雷达(RoSPR,Rover Subsurface Penetrating Radar)是一种基于火星车平台的高分辨率火星次表层地质结构探测雷达,有低频和高频两个通道,其中,低频主要用于探测火星表层以下大约100m以内的土壤和水冰结构。
火星车次表层探测雷达的低频通道使用线性调频信号作为其发射信号,并通过基于匹配滤波的脉冲压缩方法来将回波信号转换成可以识别的脉冲,经过星上处理后,雷达回波信号被变频到零中频,以复数的时间序列形式下传至地球。在火星上实际工作的时候,受环境和系统状态的影响,其信号出现时间和相位跳变,使得下传的雷达信号矩阵出现跳变和断层等不连续性,影响到后期的处理和信号解译工作。
发明内容
鉴于上述技术问题,本公开第一方面提供一种火星车雷达信号的校正方法,包括:对原始雷达数据中的每一道数据在时间域上进行插值处理,得到多道高时间分辨率雷达数据;从多道高时间分辨率雷达数据任意确定一道高时间分辨率雷达数据作为第一参考数据;将多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与第一参考数据进行时间域交叉相关分析,得到每一道高时间分辨率雷达数据对应的时间域相关曲线;将每一道高时间分辨率雷达数据的时间域相关曲线的峰值对应的时间偏移确定为该高时间分辨率雷达数据相对于第一参考数据的时间跳变量;基于时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据。
根据本公开实施例,校正方法还包括:从多道时间校正后的雷达数据中任意确定一道时间校正后的雷达数据作为第二参考数据;将多道时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,得到每一道时间校正后的雷达数据对应的相位域相关曲线;将每一道时间校正后的雷达数据对应的相位域相关曲线的实部最大值所对应的相位值确定为该道时间校正后的雷达数据相对于第二参考数据的相位跳变;基于相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据。
根据本公开实施例,对原始雷达数据中的每一道数据在时间域上进行插值处理,具体包括:采用香农插值法对原始雷达数据中的每一道数据在时间域上进行插值处理。
根据本公开实施例,将多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与第一参考数据进行时间域的交叉相关分析,具体包括:根据
Figure BDA0003722596230000021
对高时间分辨率雷达数据与第一参考数据进行时间域的交叉相关分析,其中,Ri[Δm]表示第i道高时间分辨率雷达数据的时间域相关曲线,Sref[m]为第一参考数据,Si,hires[m]为第i道高时间分辨率雷达数据,m为高时间分辨率雷达数据对应的点数,Δm为时间域相关曲线的峰值对应的时间偏移量,||表示取绝对值。
根据本公开实施例,基于时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据,具体包括:根据
S′i,hires[m]=Si,hires[m+Δm]
对每一道高时间分辨率雷达数据的时间进行校正,其中,Si,hires[m]为第i道高时间分辨率雷达数据,S′i,hires[m]为第i道高时间分辨率雷达数据经时间校正后的雷达数据。
根据本公开实施例,将多道时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,具体包括:根据
Figure BDA0003722596230000031
对时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,其中,fi[φ]为第i道时间校正后的雷达数据的相位域相关曲线,S′i[n]为第i道时间校正后的雷达数据,S′ref[n]为第二参考数据,φi为表示第i道时间校正后的雷达数据相对于第二参考数据的相位跳变,*上标代表取共轭。
根据本公开实施例,基于相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据,具体包括:
S″i[n]=S′i[n].e
对每一道时间校正后的雷达数据的相位进行校正,其中,S′i[n]为第i道时间校正后的雷达数据,Si[n]为时间校正后的雷达数据经相位校正后的雷达数据。
本公开第一方面提供一种火星车雷达信号的校正装置,包括:插值模块,用于对原始雷达数据中的每一道数据在时间域上进行插值处理,得到多道高时间分辨率雷达数据;第一确定模块,用于从多道高时间分辨率雷达数据任意确定一道高时间分辨率雷达数据作为第一参考数据;第一相关分析模块,用于将多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与第一参考数据进行时间域的交叉相关分析,得到每一道高时间分辨率雷达数据对应的时间域相关曲线;第二确定模块,用于将每一道高时间分辨率雷达数据的时间域相关曲线的峰值对应的时间偏移确定为该高时间分辨率雷达数据相对于第一参考数据的时间跳变量;第一校正模块,用于基于时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据;第三确定模块,用于从多道时间校正后的雷达数据中任意确定一道时间校正后的雷达数据作为第二参考数据;第一相关分析模块,用于将多道时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,得到每一道时间校正后的雷达数据对应的相位域相关曲线;第四确定模块,用于将每一道时间校正后的雷达数据对应的相位域相关曲线的实部最大值所对应的相位值确定为该道时间校正后的雷达数据相对于第二参考数据的相位跳变;第二校正模块,用于基于相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据。
本公开第三方面还提供一种电子设备,包括:一个或多个处理器;存储器,用于存储一个或多个程序,其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现上述所述的方法。
本公开第四方面还提供一种计算机可读存储介质,其特征在于,其上存储有可执行指令,该指令被处理器执行时使处理器实现上述所述的方法。
根据本公开实施例提供的火星车雷达信号的校正方法、装置、设备及介质,至少能够实现以下技术效果:
通过直接从火星车雷达数据中确定参考数据,以参考数据为基准对雷达数据进行时间和相位域的交叉相关分析,来弥补雷达系统受环境影响而出现的时间相位跳变情况。在经过时间相位校正后,雷达信号在行走路线方向上的不连续性得到修复,为后续雷达信号处理和解译具有重要意义。
附图说明
通过以下参照附图对本公开实施例的描述,本公开的上述以及其他目的、特征和优点将更为清楚,在附图中:
图1示意性示出了根据本公开一实施例的火星车雷达信号的校正方法的流程图。
图2示意性示出了根据本公开实施例的经过插值后的火星车雷达信号的实部和虚部的曲线图。
图3示意性示出了根据本公开另一实施例的火星车雷达信号的校正方法的流程图。
图4示意性示出了根据本公开实施例的火星车雷达信号经时间校正前后的信号曲线对比图。
图5示意性示出了根据本公开实施例的火星车雷达信号经时间和相位校正前后的信号曲线对比图。
图6示意性示出了根据本公开实施例的火星车雷达信号校正装置的框图。
图7示意性示出了根据本发明实施例的适于实现上文描述的方法的电子设备的框图。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本公开进一步详细说明。显然,所描述的实施例是本公开一部分实施例,而不是全部的实施例。基于本公开中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本公开保护的范围。
在此使用的术语仅仅是为了描述具体实施例,而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在,但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。
在本公开中,除非另有明确的规定和限定,术语“安装”、“相连”“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接或可以互相通讯;可以是直接连接,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本公开中的具体含义。
在本公开的描述中,需要理解的是,术语“纵向”、“长度”、“周向”、“前”、“后”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本公开和简化描述,而不是指示或暗示所指的子系统或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本公开的限制。
贯穿附图,相同的元素由相同或相近的附图标记来表示。可能导致本公开的理解造成混淆时,将省略常规结构或构造。并且图中各部件的形状、尺寸、位置关系不反映真实大小、比例和实际位置关系。另外,在本公开中,不应将位于括号之间的任何参考符号构造成对本公开的限制。
类似地,为了精简本公开并帮助理解各个公开方面中的一个或多个,在上面对本公开示例性实施例的描述中,本公开的各个特征有时被一起分到单个实施例、图或者对其描述中。参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本公开的至少一个实施例或示例中。本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或者多个实施例或示例中以合适的方式结合。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。因此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本公开的描述中,“多个”的含义是至少两个,例如两个、三个等,除非另有明确具体的限定。
针对现有技术的不足,本公开实施例提供一种火星车雷达信号的校正方法,可用于对火星车次表层探测雷达的低频通道数据的时间跳变进行估计和校正,以便进行后续的雷达信号处理和解译的工作。需要说明的是,本公开实施例提供的方法不仅限于应用于火星车雷达信号的校正,对于其他需要进行时间和相位校正的雷达信号也适用,具体可以根据实际应用需求选择,本公开不做限制。下面结合具体的实施例进行详细介绍。
图1示意性示出了根据本公开一实施例的火星车雷达信号的校正方法的流程图。
如图1所示,该火星车雷达信号的校正方法例如可以包括操作S101~操作S104。
在操作S101,对原始雷达数据中的每一道数据在时间域上进行插值处理,得到多道高时间分辨率雷达数据。
由于火星车雷达低频数据的时间跳变很难恰好为低频数据的采样间隔的整数倍,因此会存在所谓“亚采样间隔”的时间跳变,影响时间跳变的定位精度。可以通过插值提高火星车雷达低频数据的时间分辨率,从而提升影响时间跳变的定位精度。
在本公开一实施例中,可以采用香农插值法对原始雷达数据中的每一道数据在时间域上进行插值处理。
具体地,为了便于后续计算,可以先将原始的雷达数据转化为矩阵,对于雷达数据矩阵中的第i道数据进行香农插值,记该道原始雷达数据为Si[n],n∈[1,N],则可以根据下式:
Figure BDA0003722596230000071
进行香农插值。其中,Si(t)是理论上香农插值可以重建得到的连续信号,T是原始雷达数据的采样时间间隔,sin c为归一化辛格函数。
图2示意性示出了根据本公开实施例的经过插值后的低频通道雷达信号的实部和虚部的曲线图。
如图2所示,在插值时,原始信号对应的时间序列为t=[0,T,...,(N-1)·T],对信号进行预设倍数的插值,输出新的时间序列变为原始时间序列的预设倍数分之一,例如,对信号进行100倍插值,新的时间序列为t=[0,T/100,...,(NT-T)/100]。通过香农插值以后,输出数据记为Si,hires[m],m∈[1,100N]。
又由于火星车数据中不仅存在时间跳变,还存在相位跳变,因此,可以通过取绝对值的方法可以消除相位的影响,从而精准确定出时间跳变量。示例性的,取绝对值时,对雷达数据矩阵中的每一个元素,取其绝对值,记为|Si,hires[m]|,组合即可得到多道高时间分辨率雷达数据对应的矩阵形式。
在操作S102,从多道高时间分辨率雷达数据任意确定一道高时间分辨率雷达数据作为第一参考数据。
在本公开实施例中,理论上可以任意从数据中选取一道数据作为第一参考信号,记为|Sref[m]|。例如,可以选取多道高时间分辨率雷达数据中的第一道高时间分辨率雷达数据作为参考信号,记为|Sref[m]|=|S1,hires[m]|。
在操作S103,将多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与第一参考数据进行时间域的交叉相关分析,得到每一道高时间分辨率雷达数据对应的时间域相关曲线。
可以根据
Figure BDA0003722596230000081
对高时间分辨率雷达数据与第一参考数据进行时间域的交叉相关分析,其中,Ri[Δm]表示第i道高时间分辨率雷达数据的时间域相关曲线,Sref[m]为参考信号,Si,hires[m]为第i道高时间分辨率雷达数据,m为高时间分辨率雷达数据对应的点数,Δm为时间域相关曲线的峰值对应的时间偏移量,||表示取绝对值。
在操作S104,将每一道高时间分辨率雷达数据的时间域相关曲线的峰值对应的时间偏移确定为该高时间分辨率雷达数据相对于第一参考数据的时间跳变量。
在本公开实施例中,可以取时间域相关曲线的峰值对应的时间偏移量Δm为时间分辨率雷达数据相对于第一参考数据的时间跳变量。
在操作S105,基于时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据。
可以根据
S′i,hires[m]=Si,hires[m+Δm]
对每一道高时间分辨率雷达数据的时间进行校正,其中,Si,hires[m]为第i道高时间分辨率雷达数据,S′i,hires[m]为第i道高时间分辨率雷达数据经时间校正后的雷达数据。当m≤0或m>100N时,信号值为0。
之后,对S′i,hires[m]进行降采样,回复原始采样间隔,得到原始采样间隔的时间校正的雷达数据,记为Si[n],降采样的方式例如可以是n倍抽取。例如,取S′i,hires[m]的第[100,200,,...,100N]点组成原始采样间隔的时间校正后的信号S′i[n]。
根据本公开实施例提供的校正方法,以参考数据为基准对雷达数据进行时间的交叉相关分析,来弥补雷达系统受环境影响而出现的时间跳变情况。
图3示意性示出了根据本公开另一实施例的火星车雷达信号的校正方法的流程图。
如图3所示,本实施例提供的校正方法与图1所示的校正方法相比,在操作S105之后,还可以包括操作S106~操作S109。
在操作S106,从多道时间校正后的雷达数据中任意确定一道时间校正后的雷达数据作为第二参考数据。
在本公开实施例中,理论上可以任意从时间校正后的雷达数据中选取一道数据作为第二参考信号S′ref[n],记为,例如,仍然选取第一道数据作为第二参考数据S′ref[n]=S1[n]。
在操作S107,将多道时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,得到每一道时间校正后的雷达数据对应的相位域相关曲线。
在本公开实施例中,位域的交叉相关分析具体可以包括:根据
Figure BDA0003722596230000091
对时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,其中,fi[φ]为第i道时间校正后的雷达数据的相位域相关曲线,S′i[n]为第i道时间校正后的雷达数据,S′ref[n]为第二参考数据,φi为表示第i道时间校正后的雷达数据相对于第二参考数据的相位跳变,*上标代表取共轭。
在操作S108,将每一道时间校正后的雷达数据对应的相位域相关曲线的实部最大值所对应的相位值确定为该道时间校正后的雷达数据相对于第二参考数据的相位跳变。
在本公开实施例中,可以取时间域相关曲线的实部的最大值所对应的相位值φi作为第i道时间校正后的雷达数据相对于第二参考数据的相位跳变。
在操作S109,基于相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据。
在本公开实施例中,相位校正的具体过程可以为:
根据
S″i[n]=S′i[n]·e
对每一道时间校正后的雷达数据的相位进行校正,其中,S′i[n]为第i道时间校正后的雷达数据,S″i[n]为时间校正后的雷达数据经相位校正后的雷达数据。
根据本公开实施例提供的校正方法,在时间校正的基础上,以参考数据为基准对雷达数据进行相位域的交叉相关分析,进一步弥补雷达系统受环境影响而出现的相位跳变情况。
为了验证本公开实施例提供的校正方法的可行性及准确性,本公开实施提供了雷达数据的时间校正及相位校正结果。
图4示意性示出了根据本公开实施例的低频通道雷达信号经时间校正前后的信号曲线对比图。
图5示意性示出了根据本公开实施例的低频通道雷达信号经时间和相位校正前后的信号曲线对比图。
如图4和图5所示,采用本公开实施例提供的校正方法,以参考数据为基准对雷达数据进行时间和相位域的交叉相关分析,来弥补雷达系统受环境影响而出现的时间相位跳变情况,在经过时间相位校正后,雷达信号在行走路线方向上的不连续性得到修复,保证了雷达信号的精度,对后续雷达信号处理和解译具有重要意义。
图6示意性示出了根据本公开实施例的低频通道雷达信号校正装置的框图。
如图6所示,低频通道雷达信号校正装置600可以包括插值模块610、第一确定模块620、第一相关分析模块630、第二确定模块640、第一校正模块650、第三确定模块660、第一相关分析模块670、第四确定模块680以及第二校正模块690。
插值模块610,用于对原始雷达数据中的每一道数据在时间域上进行插值处理,得到多道高时间分辨率雷达数据。
第一确定模块620,用于从多道高时间分辨率雷达数据任意确定一道高时间分辨率雷达数据作为第一参考数据。
第一相关分析模块630,用于将多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与第一参考数据进行时间域的交叉相关分析,得到每一道高时间分辨率雷达数据对应的时间域相关曲线。
第二确定模块640,用于将每一道高时间分辨率雷达数据的时间域相关曲线的峰值对应的时间偏移确定为该高时间分辨率雷达数据相对于第一参考数据的时间跳变量。
第一校正模块650,用于基于时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据。
第三确定模块660,用于从多道时间校正后的雷达数据中任意确定一道时间校正后的雷达数据作为第二参考数据。
第一相关分析模块670,用于将多道时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,得到每一道时间校正后的雷达数据对应的相位域相关曲线。
第四确定模块680,用于将每一道时间校正后的雷达数据对应的相位域相关曲线的实部最大值所对应的相位值确定为该道时间校正后的雷达数据相对于第二参考数据的相位跳变。
第二校正模块690,用于基于相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据。
根据本发明的实施例的模块、子模块、单元、子单元中的任意多个、或其中任意多个的至少部分功能可以在一个模块中实现。根据本发明实施例的模块、子模块、单元、子单元中的任意一个或多个可以被拆分成多个模块来实现。根据本发明实施例的模块、子模块、单元、子单元中的任意一个或多个可以至少被部分地实现为硬件电路,例如现场可编程门阵列(FPGA)、可编程逻辑阵列(PLA)、片上系统、基板上的系统、封装上的系统、专用集成电路(ASIC),或可以通过对电路进行集成或封装的任何其他的合理方式的硬件或固件来实现,或以软件、硬件以及固件三种实现方式中任意一种或以其中任意几种的适当组合来实现。或者,根据本发明实施例的模块、子模块、单元、子单元中的一个或多个可以至少被部分地实现为计算机程序模块,当该计算机程序模块被运行时,可以执行相应的功能。
例如,插值模块610、第一确定模块620、第一相关分析模块630、第二确定模块640、第一校正模块650、第三确定模块660、第一相关分析模块670、第四确定模块680以及第二校正模块690中的任意多个可以合并在一个模块/单元/子单元中实现,或者其中的任意一个模块/单元/子单元可以被拆分成多个模块/单元/子单元。或者,这些模块/单元/子单元中的一个或多个模块/单元/子单元的至少部分功能可以与其他模块/单元/子单元的至少部分功能相结合,并在一个模块/单元/子单元中实现。根据本发明的实施例插值模块610、第一确定模块620、第一相关分析模块630、第二确定模块640、第一校正模块650、第三确定模块660、第一相关分析模块670、第四确定模块680以及第二校正模块690中的至少一个可以至少被部分地实现为硬件电路,例如现场可编程门阵列(FPGA)、可编程逻辑阵列(PLA)、片上系统、基板上的系统、封装上的系统、专用集成电路(ASIC),或可以通过对电路进行集成或封装的任何其他的合理方式等硬件或固件来实现,或以软件、硬件以及固件三种实现方式中任意一种或以其中任意几种的适当组合来实现。或者,插值模块610、第一确定模块620、第一相关分析模块630、第二确定模块640、第一校正模块650、第三确定模块660、第一相关分析模块670、第四确定模块680以及第二校正模块690中的至少一个可以至少被部分地实现为计算机程序模块,当该计算机程序模块被运行时,可以执行相应的功能。
需要说明的是,本发明的实施例的低频通道雷达信号校正装置部分与本发明的实施例的低频通道雷达信号校正方法部分是相对应的,其具体实施细节及带来的技术效果也是相同的,在此不再赘述。
图7示意性示出了根据本发明实施例的适于实现上文描述的方法的电子设备的框图。图7示出的电子设备仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。
如图7所示,根据本发明实施例的电子设备700包括处理器701,其可以根据存储在只读存储器(ROM)702中的程序或者从存储部分708加载到随机访问存储器(RAM)703中的程序而执行各种适当的动作和处理。处理器701例如可以包括通用微处理器(例如CPU)、指令集处理器和/或相关芯片组和/或专用微处理器(例如,专用集成电路(ASIC)),等等。处理器701还可以包括用于缓存用途的板载存储器。处理器701可以包括用于执行根据本发明实施例的方法流程的不同动作的单一处理单元或者是多个处理单元。
在RAM 703中,存储有电子设备700操作所需的各种程序和数据。处理器701、ROM702以及RAM703通过总线704彼此相连。处理器701通过执行ROM 702和/或RAM703中的程序来执行根据本发明实施例的方法流程的各种操作。需要注意,所述程序也可以存储在除ROM702和RAM 703以外的一个或多个存储器中。处理器701也可以通过执行存储在所述一个或多个存储器中的程序来执行根据本发明实施例的方法流程的各种操作。
根据本发明的实施例,电子设备700还可以包括输入/输出(I/O)接口705,输入/输出(I/O)接口705也连接至总线704。电子设备700还可以包括连接至I/O接口705的以下部件中的一项或多项:包括键盘、鼠标等的输入部分706;包括诸如阴极射线管(CRT)、液晶显示器(LCD)等以及扬声器等的输出部分707;包括硬盘等的存储部分708;以及包括诸如LAN卡、调制解调器等的网络接口卡的通信部分709。通信部分709经由诸如因特网的网络执行通信处理。驱动器710也根据需要连接至I/O接口705。可拆卸介质711,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器710上,以便于从其上读出的计算机程序根据需要被安装入存储部分708。
根据本发明的实施例,根据本发明实施例的方法流程可以被实现为计算机软件程序。例如,本发明的实施例包括一种计算机程序产品,其包括承载在计算机可读存储介质上的计算机程序,该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分709从网络上被下载和安装,和/或从可拆卸介质711被安装。在该计算机程序被处理器701执行时,执行本发明实施例的系统中限定的上述功能。根据本发明的实施例,上文描述的系统、设备、装置、模块、单元等可以通过计算机程序模块来实现。
本发明还提供了一种计算机可读存储介质,该计算机可读存储介质可以是上述实施例中描述的设备/装置/系统中所包含的;也可以是单独存在,而未装配入该设备/装置/系统中。上述计算机可读存储介质承载有一个或者多个程序,当上述一个或者多个程序被执行时,实现根据本发明实施例的方法。
根据本发明的实施例,计算机可读存储介质可以是非易失性的计算机可读存储介质。例如可以包括但不限于:便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本发明中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
例如,根据本发明的实施例,计算机可读存储介质可以包括上文描述的ROM702和/或RAM 703和/或ROM 702和RAM 703以外的一个或多个存储器。
附图中的流程图和框图,图示了按照本发明各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,上述模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图或流程图中的每个方框、以及框图或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。本领域技术人员可以理解,本发明的各个实施例中记载的特征可以进行多种组合和/或结合,即使这样的组合或结合没有明确记载于本发明中。特别地,在不脱离本发明精神和教导的情况下,本发明的各个实施例记载的特征可以进行多种组合和/或结合。所有这些组合和/或结合均落入本发明的范围。

Claims (10)

1.一种火星车雷达信号的校正方法,其特征在于,包括:
对原始雷达数据中的每一道数据在时间域上进行插值处理,得到多道高时间分辨率雷达数据;
从所述多道高时间分辨率雷达数据任意确定一道高时间分辨率雷达数据作为第一参考数据;
将所述多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与所述第一参考数据进行时间域的交叉相关分析,得到每一道高时间分辨率雷达数据对应的时间域相关曲线;
将每一道高时间分辨率雷达数据的时间域相关曲线的峰值对应的时间偏移确定为该高时间分辨率雷达数据相对于第一参考数据的时间跳变量;
基于所述时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据。
2.根据权利要求1所述的火星车雷达信号的校正方法,其特征在于,所述校正方法还包括:
从多道所述时间校正后的雷达数据中任意确定一道时间校正后的雷达数据作为第二参考数据;
将多道所述时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与所述第二参考数据进行相位域的交叉相关分析,得到每一道时间校正后的雷达数据对应的相位域相关曲线;
将每一道时间校正后的雷达数据对应的相位域相关曲线的实部最大值所对应的相位值确定为该道时间校正后的雷达数据相对于第二参考数据的相位跳变;
基于所述相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据。
3.根据权利要求1或2所述的火星车雷达信号的校正方法,其特征在于,所述对原始雷达数据中的每一道数据在时间域上进行插值处理,具体包括:
采用香农插值法对原始雷达数据中的每一道数据在时间域上进行插值处理。
4.根据权利要求1或2所述的火星车雷达信号的校正方法,其特征在于,所述将所述多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与所述第一参考数据进行时间域的交叉相关分析,具体包括:
根据
Figure FDA0003722596220000021
对高时间分辨率雷达数据与第一参考数据进行时间域的交叉相关分析,其中,Ri[Δm]表示第i道高时间分辨率雷达数据的时间域相关曲线,Sref[m]为第一参考数据,Si,hires[m]为第i道高时间分辨率雷达数据,m为高时间分辨率雷达数据对应的点数,Δm为时间域相关曲线的峰值对应的时间偏移量,||表示取绝对值。
5.根据权利要求4所述的火星车雷达信号的校正方法,其特征在于,所述基于所述时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据,具体包括:
根据
Si,hires[m]=Si,hires[m+Δm]
对每一道高时间分辨率雷达数据的时间进行校正,其中,Si,hires[m]为第i道高时间分辨率雷达数据,S′i,hires[m]为第i道高时间分辨率雷达数据经时间校正后的雷达数据。
6.根据权利要求2所述的火星车雷达信号的校正方法,其特征在于,所述将多道所述时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与所述第二参考数据进行相位域的交叉相关分析,具体包括:
根据
Figure FDA0003722596220000022
对时间校正后的雷达数据与第二参考数据进行相位域的交叉相关分析,其中,fi[φ]为第i道时间校正后的雷达数据的相位域相关曲线,Si′[n]为第i道时间校正后的雷达数据,S′ref[n]为所述第二参考数据,φi为表示第i道时间校正后的雷达数据相对于第二参考数据的相位跳变,*上标代表取共轭。
7.根据权利要求6所述的火星车雷达信号的校正方法,其特征在于,所述基于所述相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据,具体包括:
根据
Si″[n]=S′i[n]·e
对每一道时间校正后的雷达数据的相位进行校正,其中,Si′[n]为第i道时间校正后的雷达数据,Si[n]为时间校正后的雷达数据经相位校正后的雷达数据。
8.一种火星车雷达信号的校正装置,其特征在于,包括:
插值模块,用于对原始雷达数据中的每一道数据在时间域上进行插值处理,得到多道高时间分辨率雷达数据;
第一确定模块,用于从所述多道高时间分辨率雷达数据任意确定一道高时间分辨率雷达数据作为第一参考数据;
第一相关分析模块,用于将所述多道高时间分辨率雷达数据中除第一参考数据外的其他高时间分辨率雷达数据与所述第一参考数据进行时间域的交叉相关分析,得到每一道高时间分辨率雷达数据对应的时间域相关曲线;
第二确定模块,用于将每一道高时间分辨率雷达数据的时间域相关曲线的峰值对应的时间偏移确定为该高时间分辨率雷达数据相对于第一参考数据的时间跳变量;
第一校正模块,用于基于所述时间跳变量对每一道高时间分辨率雷达数据的时间进行校正,得到时间校正后的雷达数据;
第三确定模块,用于从多道所述时间校正后的雷达数据中任意确定一道时间校正后的雷达数据作为第二参考数据;
第一相关分析模块,用于将多道所述时间校正后的雷达数据中除第二参考数据外的其他时间校正后的雷达数据与所述第二参考数据进行相位域的交叉相关分析,得到每一道时间校正后的雷达数据对应的相位域相关曲线;
第四确定模块,用于将每一道时间校正后的雷达数据对应的相位域相关曲线的实部最大值所对应的相位值确定为该道时间校正后的雷达数据相对于第二参考数据的相位跳变;
第二校正模块,用于基于所述相位跳变对每一道时间校正后的雷达数据进行相位校正,得到相位校正后的雷达数据。
9.一种电子设备,其特征在于,包括:
一个或多个处理器;
存储器,用于存储一个或多个程序,
其中,当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现权利要求1至7中任一项所述的方法。
10.一种计算机可读存储介质,其特征在于,其上存储有可执行指令,该指令被处理器执行时使处理器实现权利要求1至7中任一项所述的方法。
CN202210776198.7A 2022-06-30 2022-06-30 火星车雷达信号的校正方法、装置、设备及介质 Active CN115079110B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210776198.7A CN115079110B (zh) 2022-06-30 2022-06-30 火星车雷达信号的校正方法、装置、设备及介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210776198.7A CN115079110B (zh) 2022-06-30 2022-06-30 火星车雷达信号的校正方法、装置、设备及介质

Publications (2)

Publication Number Publication Date
CN115079110A true CN115079110A (zh) 2022-09-20
CN115079110B CN115079110B (zh) 2025-09-02

Family

ID=83256872

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210776198.7A Active CN115079110B (zh) 2022-06-30 2022-06-30 火星车雷达信号的校正方法、装置、设备及介质

Country Status (1)

Country Link
CN (1) CN115079110B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115825896A (zh) * 2022-11-28 2023-03-21 中国科学院国家天文台 校准火星车次表层探测雷达高频极化通道幅度的方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030210181A1 (en) * 2002-05-13 2003-11-13 James R. Hager Methods and apparatus for radar data processing
CN101344391A (zh) * 2008-07-18 2009-01-14 北京工业大学 基于全功能太阳罗盘的月球车位姿自主确定方法
CN106646383A (zh) * 2016-09-28 2017-05-10 中国科学院电子学研究所 一种对月壤结构探测仪数据去除背景的方法
CN110376623A (zh) * 2019-07-08 2019-10-25 中国空间技术研究院 星载gnss-r镜面反射点海洋潮汐修正定位方法和系统
CN111474524A (zh) * 2020-04-22 2020-07-31 安徽华可智能科技有限公司 一种雷达干扰装备干扰效果监测与决策支持系统
CN113567940A (zh) * 2021-07-23 2021-10-29 中山大学 基于长时间积累的宽带雷达系统通道误差估计方法、计算机装置和存储介质
CN113777564A (zh) * 2020-06-09 2021-12-10 三星电子株式会社 用于通过校正相位失真来处理雷达信号的方法和设备
US20220146664A1 (en) * 2019-03-28 2022-05-12 Sony Semiconductor Solutions Corporation Signal processing device, signal processing method, program, and information processing device
CN114594471A (zh) * 2022-02-23 2022-06-07 中国科学院国家天文台 雷达高度的测量方法和雷达高度的测量装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030210181A1 (en) * 2002-05-13 2003-11-13 James R. Hager Methods and apparatus for radar data processing
CN101344391A (zh) * 2008-07-18 2009-01-14 北京工业大学 基于全功能太阳罗盘的月球车位姿自主确定方法
CN106646383A (zh) * 2016-09-28 2017-05-10 中国科学院电子学研究所 一种对月壤结构探测仪数据去除背景的方法
US20220146664A1 (en) * 2019-03-28 2022-05-12 Sony Semiconductor Solutions Corporation Signal processing device, signal processing method, program, and information processing device
CN110376623A (zh) * 2019-07-08 2019-10-25 中国空间技术研究院 星载gnss-r镜面反射点海洋潮汐修正定位方法和系统
CN111474524A (zh) * 2020-04-22 2020-07-31 安徽华可智能科技有限公司 一种雷达干扰装备干扰效果监测与决策支持系统
CN113777564A (zh) * 2020-06-09 2021-12-10 三星电子株式会社 用于通过校正相位失真来处理雷达信号的方法和设备
CN113567940A (zh) * 2021-07-23 2021-10-29 中山大学 基于长时间积累的宽带雷达系统通道误差估计方法、计算机装置和存储介质
CN114594471A (zh) * 2022-02-23 2022-06-07 中国科学院国家天文台 雷达高度的测量方法和雷达高度的测量装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
苏宏艳;朱淮城;: "毫米波宽带雷达相位误差校正方法研究", 系统工程与电子技术, vol. 30, no. 02, 29 February 2008 (2008-02-29), pages 253 - 256 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115825896A (zh) * 2022-11-28 2023-03-21 中国科学院国家天文台 校准火星车次表层探测雷达高频极化通道幅度的方法

Also Published As

Publication number Publication date
CN115079110B (zh) 2025-09-02

Similar Documents

Publication Publication Date Title
Deller et al. DiFX-2: a more flexible, efficient, robust, and powerful software correlator
CN110018455B (zh) 星载sar成像接收通道间幅相一致性误差定标测量方法
CN104237886B (zh) 一种高精度合成孔径雷达成像方法
CN104316923A (zh) 针对合成孔径雷达bp成像的自聚焦方法
CN114240805B (zh) 一种多角度sar动态成像检测方法及装置
JP2011169869A (ja) レーダ信号処理装置
CN111751822B (zh) 一种时间和相位同步方法、系统、设备和介质
CN115079110A (zh) 火星车雷达信号的校正方法、装置、设备及介质
CN101984363A (zh) 一种步进调频体制超高分辨率sar成像方法
US20220155475A1 (en) Irregular optimized acquisition method, device, apparatus and medium for seismic data
CN110045373B (zh) 机载多通道sar成像处理方法、装置和计算机设备
CN110346797A (zh) 一种二维逆合成孔径雷达成像方法及系统
CN106054152B (zh) 基于逆扩展Omega-K算法的非理想轨迹SAR回波获取方法
JP6246338B2 (ja) 測角装置及び測角方法
JP2010060448A (ja) レーダ画像再生装置
JP2003512862A (ja) 過渡分光分析のためのデジタルフィルタのプリチャージ
CN114879174A (zh) 基于雷达的多目标测距方法以及形变监测方法
CN110426708A (zh) 一种基于方位向多通道的星载sar保相成像方法
CN104049254A (zh) 一种高分辨率扫描合成孔径雷达自聚焦处理方法和装置
CN105549010A (zh) 频域合成孔径雷达成像方法
CN117761393B (zh) 一种时域信号的获取方法及装置
CN113466859A (zh) 基于快速相位插值的自旋空间碎片目标isar二维成像方法
CN118962675A (zh) 一种基于步进频去斜和自适应相位补偿的雷达成像方法
CN114594471B (zh) 雷达高度的测量方法和雷达高度的测量装置
CN116381686A (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