CN112147643B - 一种gnss-r系统成像预处理的方法及装置 - Google Patents
一种gnss-r系统成像预处理的方法及装置 Download PDFInfo
- Publication number
- CN112147643B CN112147643B CN202010859336.9A CN202010859336A CN112147643B CN 112147643 B CN112147643 B CN 112147643B CN 202010859336 A CN202010859336 A CN 202010859336A CN 112147643 B CN112147643 B CN 112147643B
- Authority
- CN
- China
- Prior art keywords
- signal
- error
- phase error
- phase
- reflected
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000007781 pre-processing Methods 0.000 title claims abstract description 13
- 230000006835 compression Effects 0.000 claims abstract description 36
- 238000007906 compression Methods 0.000 claims abstract description 36
- 238000012545 processing Methods 0.000 claims description 26
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 238000012937 correction Methods 0.000 claims description 6
- 230000001902 propagating effect Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 abstract description 9
- 238000010586 diagram Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 230000001360 synchronised effect Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/14—Receivers specially adapted for specific applications
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本申请公开了一种GNSS‑R系统成像预处理的方法及装置,该方法包括:确定GNSS‑R系统中的直达信号以及反射信号,并根据预设的本地参考信号将直达信号进行压缩处理得到第一信号;根据第一信号确定相关峰信息以及根据相关峰信息对第一信号进行修正得到修正后的信号;确定GNSS‑R系统中接收机和发射机的位置信息,根据位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差;根据时延误差和相位误差重构生成反射参考信号,根据反射参考信号将反射信号进行压缩处理得到第二信号,根据相位误差对第二信号进行相位误差补偿得到第三信号,以使得基于第三信号进行成像。本申请解决了现有技术中系统成像效果差的技术问题。
Description
技术领域
本申请涉及GNSS-R成像技术领域,尤其涉及一种GNSS-R系统成像预处理的方法及装置。
背景技术
随着全球导航卫星的不断发展,全球导航卫星系统(global navigationsatellite system,GNSS)的微波遥感技术(GNSS-Reflections或GNSS-Remotesensing,GNSS-R)的应用也受到了广泛关注。基于GNSS-R的双基成像系统是一种典型的利用GNSS反射信号作为照射源的新型探测系统。该系统顾名思义是一种以全球导航卫星GPS(US)、GLONASS(RU)、Galileo(EU)和BDS(CN))作为发射平台,通过接收导航卫星的直达波和照射目标后的反射回波进行成像和定位。GNSS-R的双基成像系统作为双基成像系统,基于空基接收平台的GNSS-R的双基成像涉及包括如下两项重要的关键技术:一项是基于GNSS-R成像系统的信号同步;另一项是基于GNSS-R成像系统的成像方法。
目前,传统的基于GNSS-R成像算法中,信号同步是基于导航信号的传统同步方法,该同步方法存在如下技术问题:一方面,其同步过程包含捕获、跟踪等多个环节,同步算法流程复杂、计算复杂度高且建立同步的时间很长,无法保证高效性和实时性;另一方面,传统信号同步方法是针对PWM信号的同步方法,无法应用与基于GNSS-R成像体制的同步过程。而成像的过程推导大多基于地面固定平台,没有考虑同步残差补偿,从而导致成像效果差,甚至无法聚焦。
发明内容
本申请解决的技术问题是:针对现有技术中系统成像效果差的问题,提供了一种GNSS-R系统成像预处理的方法及装置,本申请实施例所提供的方案中,根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号,然后根据所第一信号确定相关峰信息以及根据相关峰信息对第一信号进行修正得到修正后的信号,然后根据GNSS-R系统中接收机和发射机的位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差,即根据本地参考信号和直达信号确定时延和相位同步误差;进一步,再根据时延误差和相位误差重构生成反射参考信号,最后,根据反射参考信号将反射信号进行压缩处理得到第二信号,根据相位误差对第二信号进行相位误差补偿得到第三信号,即根据相位误差对反射信号进行同步残差补偿,不仅将信号同步过程与成像过程有机结合在一起,还提高成像时的聚焦效果以及提高成像效果。
第一方面,本申请实施例提供一种GNSS-R系统成像预处理的方法,该方法包括:
确定GNSS-R系统中的直达信号以及反射信号,并根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号;
根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号;
确定所述GNSS-R系统中接收机和发射机的位置信息,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差;
根据所述时延误差和所述相位误差重构生成反射参考信号,根据所述反射参考信号将所述反射信号进行压缩处理得到第二信号,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,以使得基于所述第三信号进行成像。
本申请实施例所提供的方案中,根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号,然后根据所第一信号确定相关峰信息以及根据相关峰信息对第一信号进行修正得到修正后的信号,然后根据GNSS-R系统中接收机和发射机的位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差,即根据本地参考信号和直达信号确定由于发射与接收端使用了独立的本振,接收目标回波中错过了公共的时间参考,而可能引入时延和相位同步误差;进一步,再根据时延误差和相位误差重构生成反射参考信号,最后,根据反射参考信号将反射信号进行压缩处理得到第二信号,根据相位误差对第二信号进行相位误差补偿得到第三信号,即根据相位误差对反射信号进行同步残差补偿,不仅将信号同步过程与成像过程有机结合在一起,还提高成像时的聚焦效果以及提高成像效果。
可选地,根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号,包括:
根据所述第一信号确定任意相邻的前一时刻和后一时刻所对应的相关峰位置和相位,将所述前一时刻和后一时刻所对应的相关峰相位进行作差处理得到处理结果;
根据所述处理结果判断所述后一时刻所对应的相关峰相位是否发生变化;
若发生变化,则对所述第一信号中的相位进行修正得到所述修正后的信号。
可选地,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差,包括:
根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离;
根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差。
可选地,根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离,包括:
通过如下公式计算任一时刻所述直达信号的基线距离:
L(τ0,(x,y,z))=|T(tn)-R(tn)|+ξ(x,y,z)
其中,L(τ0,(x,y,z))表示所述基线距离;τ0表示所述任一时刻;T(tn)表示所述任一时刻发射机的位置信息;R(tn)表示所述任一时刻接收机的位置信息;ξ(x,y,z)表示所述偏置误差补偿因子。
可选地,根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差,包括:
通过如下公式计算得到所述时延误差和所述相位误差:
其中,表示所述时延误差;τp(tn)表示所述预设的直达信号的时延信息;c表示电磁波在空间传播的速度;表示所述相位误差;表示所述预设的直达信号的多普勒信息;λ表示电磁波在空间传播的波长。
可选地,根据所述时延误差和所述相位误差重构生成反射参考信号,包括:
通过如下公式重构生成反射参考信号:
其中,SR_ref(tn)表示所述反射参考信号;CB2A( )表示伪随机码。
可选地,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,包括:
通过如下公式对所述第二信号进行补偿得到所述第三信号:
其中,R′R_out(tn,τ)表示所述第三信号;RR_out(tn,τ)表示所述第二信号。
可选地,还包括:
根据预设的BP成像算法以及所述第三信号计算方位向压缩输出结果,并根据所述方位向压缩输出结果判断是否满足预设成像同步精度阈值;
若不满足,则调整所述偏置误差补偿因子重新计算出成像时反射信号压缩所对应的时延误差和相位误差,直达所述方位向压缩输出结果满足所述预设成像同步精度阈值为止。
第二方面,本申请实施例提供了一种GNSS-R系统成像预处理的装置,该装置包括:
处理单元,用于确定GNSS-R系统中的直达信号以及反射信号,并根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号;
修正单元,用于根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号;
计算单元,用于确定所述GNSS-R系统中接收机和发射机的位置信息,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差;
补偿单元,用于根据所述时延误差和所述相位误差重构生成反射参考信号,根据所述反射参考信号将所述反射信号进行压缩处理得到第二信号,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,以使得基于所述第三信号进行成像。
可选地,所述修正单元,具体用于:
根据所述第一信号确定任意相邻的前一时刻和后一时刻所对应的相关峰位置和相位,将所述前一时刻和后一时刻所对应的相关峰相位进行作差处理得到处理结果;
根据所述处理结果判断所述后一时刻所对应的相关峰相位是否发生变化;
若发生变化,则对所述第一信号中的相位进行修正得到所述修正后的信号。
可选地,所述计算单元,具体用于:
根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离;
根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差。
可选地,所述计算单元,具体用于:
通过如下公式计算任一时刻所述直达信号的基线距离:
L(τ0,(x,y,z))=|T(tn)-R(tn)|+ξ(x,y,z)
其中,L(τ0,(x,y,z))表示所述基线距离;τ0表示所述任一时刻;T(tn)表示所述任一时刻发射机的位置信息;R(tn)表示所述任一时刻接收机的位置信息;ξ(x,y,z)表示所述偏置误差补偿因子。
可选地,所述计算单元,具体用于:
通过如下公式计算得到所述时延误差和所述相位误差:
其中,表示所述时延误差;τp(tn)表示所述预设的直达信号的时延信息;c表示电磁波在空间传播的速度;表示所述相位误差;表示所述预设的直达信号的多普勒信息;λ表示电磁波在空间传播的波长。
可选地,所述补偿单元,具体用于:
通过如下公式重构生成反射参考信号:
其中,SR_ref(tn)表示所述反射参考信号;CB2A( )表示伪随机码。
可选地,所述补偿单元,具体用于:
通过如下公式对所述第二信号进行补偿得到所述第三信号:
其中,R′R_out(tn,τ)表示所述第三信号;RR_out(tn,τ)表示所述第二信号。
可选地,所述计算单元,还用于:
根据预设的BP成像算法以及所述第三信号计算方位向压缩输出结果,并根据所述方位向压缩输出结果判断是否满足预设成像同步精度阈值;
若不满足,则调整所述偏置误差补偿因子重新计算出成像时反射信号压缩所对应的时延误差和相位误差,直达所述方位向压缩输出结果满足所述预设成像同步精度阈值为止。
附图说明
图1为本申请实施例所提供的一种GNSS-R系统成像预处理的方法的流程示意图;
图2为本申请实施例所提供的一种峰值相位的差分示意图;
图3为本申请实施例所提供的一种GNSS-R系统成像预处理的装置的结构示意图。
具体实施方式
本申请实施例提供的方案中,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本申请保护的范围。
以下结合说明书附图对本申请实施例所提供的一种GNSS-R系统成像预处理的方法做进一步详细的说明,该方法具体实现方式可以包括以下步骤(方法流程如图1所示):
步骤101,确定GNSS-R系统中的直达信号以及反射信号,并根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号。
在本申请实施例所提供的方案中,在GNSS-R系统中同步和成像过程首先需要通过本地参考信号对直达信号进行匹配滤波,例如,本地参考信号可以是本地生成的P代码。由于在GNSS-R系统中导航信息的值可以为1或-1,即导航信息可以表示为一个相位项,并且导航信息的振幅AP和初始相位在以下的处理中可以忽略,考虑到所有这些因素,下面以北斗B2A信号为例,对上述确定GNSS-R系统中的直达信号以及反射信号的过程进行简要说明。
具体的,卫星发射的发射复信号通过下式表示:
SB2A(t)=CB2A[t]·exp[-j2πf0t+jπDs(t)] (1)
其中,SB2A(t)表示发射复信号;CB2A[t]表示预设的伪随机码;f0表示载波频率;Ds(t)表示导航信息的值。
进一步,根据上述公式(1)所示的发射复信号得到直达信号,以及将直达信号进行正交解调得到解调后的直达信号,具体的,解调后的直达信号通过下式表示:
其中,表示解调后的直达信号;τH(tn)表示直达信号时延;τδ(tn)表示系统同步时延误差;fd(tn)表示直达信号对应的瞬时多普勒频率;表示系统同步相位误差。
在本申请实施例所提供的方案中,确定反射信号的过程与上述确定直达信号的过程类似,在此不赘述。
进一步,在确定出直达信号之后,需要根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号,其中,第一信号是指匹配滤波后的信号。具体的,第一信号可通过如下表达式表示:
其中,RH_B2A(tn,τ)表示第一信号;SH_B2A_ref(τ)表示预设的本地参考信号;Λ[·]表示B2A码匹配滤波输出结果。
步骤102,根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号。
在本申请实施例所提供的方案中,在确定出第一信后之后,根据第一信号可以确定相关峰信息,其中,相关峰信息包括相关峰的位置以相关峰的相位。具体的,根据上述公式(3)中第一信号可以确定出相关峰的位置和相关峰的相位如下所示:
τpeak=τH(tn)+τδ(tn)
其中,τpeak表示相关峰的位置;表示相关峰的相位。
进一步,根据上述公式(4)可知,πDs(tn)这一项为±π,即意味着相位跳变将出现在导航信息被调制的地方,因此,为了避免相位跳变对第一信号的影响,需要对第一信号进行修正得到修正后的信号。具体的,在本申请实施例所提供的方案中,根据第一信号确定相关峰信息以及根据相关峰信息对第一信号进行修正得到修正后的信号方式有多种,下面以一种较佳的方式为例进行说明。
在一种可能实现的方式中,根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号,包括:根据所述第一信号确定任意相邻的前一时刻和后一时刻所对应的相关峰位置和相位,将所述前一时刻和后一时刻所对应的相关峰相位进行作差处理得到处理结果;根据所述处理结果判断所述后一时刻所对应的相关峰相位是否发生变化;若发生变化,则对所述第一信号中的相位进行修正得到所述修正后的信号。
在本申请实施例所提供的方案中,可以通过对相邻的前一时刻和后一时刻的峰值相位进行差分处理来检测导航消息是否被调制。具体的,参见图2,本申请实施例提供了一种峰值相位的差分示意图。前一时刻和后一时刻的峰值相位之间的差值可通过如下公式表示:
其中,表示前一时刻和后一时刻的峰值相位之间的差值;ΔτH(tn)表示脉冲重复间隔(Pulse Repetition Interval,PRI)的差分时延;ΔDs(tn)表示导航信息引起的相位变化;表示相位同步误差波动。
进一步,根据值判断相关峰相位是否发生变化;若发生变化,则对所述第一信号中的相位进行修正得到所述修正后的信号。具体的,修正后的信号参见下式所示:
根据上述公式(6)可确定,此时GNSS-R系统中多普勒相位项如下所示:
步骤103,确定所述GNSS-R系统中接收机和发射机的位置信息,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差。
具体的,根据发射机和接收机的位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差的方式有多种,下面以一种较佳的方式为例进行说明。
在一种可能实现的方式中,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差,包括:根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离;根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差。
具体的,直达信号的基线距离为卫星到达接收机位置的理论计算距离。
进一步,在一种可能实现的方式中,根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离,包括:
通过如下公式计算任一时刻所述直达信号的基线距离:
L(τ0,(x,y,z))=|T(tn)-R(tn)|+ξ(x,y,z)
其中,L(τ0,(x,y,z))表示所述基线距离;τ0表示所述任一时刻;T(tn)表示所述任一时刻发射机的位置信息;R(tn)表示所述任一时刻接收机的位置信息;ξ(x,y,z)表示所述偏置误差补偿因子。
进一步,在一种可能实现的方式中,根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差,包括:
通过如下公式计算得到所述时延误差和所述相位误差:
其中,表示所述时延误差;τp(tn)表示所述预设的直达信号的时延信息;c表示电磁波在空间传播的速度;表示所述相位误差;表示所述预设的直达信号的多普勒信息;λ表示电磁波在空间传播的波长。
步骤104,根据所述时延误差和所述相位误差重构生成反射参考信号,根据所述反射参考信号将所述反射信号进行压缩处理得到第二信号,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,以使得基于所述第三信号进行成像。
具体的,在计算出成像时反射信号压缩所对应的时延误差和相位误差之后,根据时延误差和相位误差重构生成反射参考信号。在本申请实施例所提供的方案中,根据时延误差和相位误差重构生成反射参考信号的方式有多种,下面以一种较佳的方式为例进行说明。
在一种可能实现的方式中,根据所述时延误差和所述相位误差重构生成反射参考信号,包括:
通过如下公式重构生成反射参考信号:
其中,SR_ref(tn)表示所述反射参考信号;CB2A( )表示伪随机码。
进一步,根据所述反射参考信号将所述反射信号进行压缩处理得到第二信号,具体的,第二信号可以通过如下公式表示:
进一步,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号。具体的,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号的方式有多种,下面以一种较佳的方式为例进行说明。
在一种可能实现的方式中,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,包括:
通过如下公式对所述第二信号进行补偿得到所述第三信号:
其中,R′R_out(tn,τ)表示所述第三信号;RR_out(tn,τ)表示所述第二信号。
根据上述第三信号的表达式可以确定反射信号的时延信息和相位信息,具体时延信息和相位信息如下所示:
其中,τp_out表示反射信号的时延信息;表示反射信号的相位信息。
根据上述公式(9)和公式(10)可以看出GNSS-R系统中同步时延误差τδ(tn)和同步相位误差已完全剔除。
进一步,为了保证处理后的反射信号满足同步成像的精度要求,在一种可能实现的方式中,本申请实施例所提供的方案,在步骤104之后,还包括:
根据预设的BP成像算法以及所述第三信号计算方位向压缩输出结果,并根据所述方位向压缩输出结果判断是否满足预设成像同步精度阈值;
若不满足,则跳转到步骤103,调整所述偏置误差补偿因子重新计算出成像时反射信号压缩所对应的时延误差和相位误差,直达所述方位向压缩输出结果满足所述预设成像同步精度阈值为止。
为了验证上述预处理后的反射信号满足成像同步精度要求,下面对得到的反射信号进行验证,下面对反射信号验证过程进行简要介绍。
首先,基于BP成像算法的点目标A累积曲线对第三信号进行积分。累积曲线由时延的理论计算确定第三信号的方位向压缩输出,具体的,第三信号的方位向压缩输出可通过下式表示:
I(A)=∫R′R_out(tn,τ)·exp[jφout(tn)]dtn (11)
其中, 表示合成孔径时间t0时点目标A的理论双基距离。
进一步,可将BP算法的方位向压缩表示为:
根据上述公式(12)可以确定出距离向与方位向模糊度函数已经分离成上述CP和sinc两个分量,故通过本申请实施例所提供的方案所获得的图像理论上具有良好的聚焦性,并且相位信息也得到了保持。
本申请实施例所提供的方案中,根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号,然后根据所第一信号确定相关峰信息以及根据相关峰信息对第一信号进行修正得到修正后的信号,然后根据GNSS-R系统中接收机和发射机的位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差,即根据本地参考信号和直达信号确定由于发射与接收端使用了独立的本振,接收目标回波中错过了公共的时间参考,而可能引入时延和相位同步误差;进一步,再根据时延误差和相位误差重构生成反射参考信号,最后,根据反射参考信号将反射信号进行压缩处理得到第二信号,根据相位误差对第二信号进行相位误差补偿得到第三信号,即根据相位误差对反射信号进行同步残差补偿,不仅将信号同步过程与成像过程有机结合在一起,还提高成像时的聚焦效果以及提高成像效果。
基于与上述图1所示的方法相同的发明构思,本申请实施例提供了一种GNSS-R系统成像预处理的装置,参见图3,该装置包括:
处理单元301,用于确定GNSS-R系统中的直达信号以及反射信号,并根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号;
修正单元302,用于根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号;
计算单元303,用于确定所述GNSS-R系统中接收机和发射机的位置信息,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差;
补偿单元304,用于根据所述时延误差和所述相位误差重构生成反射参考信号,根据所述反射参考信号将所述反射信号进行压缩处理得到第二信号,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,以使得基于所述第三信号进行成像。
可选地,所述修正单元302,具体用于:
根据所述第一信号确定任意相邻的前一时刻和后一时刻所对应的相关峰位置和相位,将所述前一时刻和后一时刻所对应的相关峰相位进行作差处理得到处理结果;
根据所述处理结果判断所述后一时刻所对应的相关峰相位是否发生变化;
若发生变化,则对所述第一信号中的相位进行修正得到所述修正后的信号。
可选地,所述计算单元303,具体用于:
根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离;
根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差。
可选地,所述计算单元303,具体用于:
通过如下公式计算任一时刻所述直达信号的基线距离:
L(τ0,(x,y,z))=|T(tn)-R(tn)|+ξ(x,y,z)
其中,L(τ0,(x,y,z))表示所述基线距离;τ0表示所述任一时刻;T(tn)表示所述任一时刻发射机的位置信息;R(tn)表示所述任一时刻接收机的位置信息;ξ(x,y,z)表示所述偏置误差补偿因子。
可选地,所述计算单元303,具体用于:
通过如下公式计算得到所述时延误差和所述相位误差:
其中,表示所述时延误差;τp(tn)表示所述预设的直达信号的时延信息;c表示电磁波在空间传播的速度;表示所述相位误差;表示所述预设的直达信号的多普勒信息;λ表示电磁波在空间传播的波长。
可选地,所述补偿单元304,具体用于:
通过如下公式重构生成反射参考信号:
其中,SR_ref(tn)表示所述反射参考信号;CB2A( )表示伪随机码。
可选地,所述补偿单元304,具体用于:
通过如下公式对所述第二信号进行补偿得到所述第三信号:
其中,R′R_out(tn,τ)表示所述第三信号;RR_out(tn,τ)表示所述第二信号。
可选地,所述计算单元303,还用于:
根据预设的BP成像算法以及所述第三信号计算方位向压缩输出结果,并根据所述方位向压缩输出结果判断是否满足预设成像同步精度阈值;
若不满足,则调整所述偏置误差补偿因子重新计算出成像时反射信号压缩所对应的时延误差和相位误差,直达所述方位向压缩输出结果满足所述预设成像同步精度阈值为止。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然,本领域的技术人员可以对本申请进行各种改动和变型而不脱离本申请的精神和范围。这样,倘若本申请的这些修改和变型属于本申请权利要求及其等同技术的范围之内,则本申请也意图包含这些改动和变型在内。
Claims (10)
1.一种GNSS-R系统成像预处理的方法,其特征在于,包括:
确定GNSS-R系统中的直达信号以及反射信号,并根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号;
根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号;
确定所述GNSS-R系统中接收机和发射机的位置信息,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差;
根据所述时延误差和所述相位误差重构生成反射参考信号,根据所述反射参考信号将所述反射信号进行压缩处理得到第二信号,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,以使得基于所述第三信号进行成像。
2.如权利要求1所述的方法,其特征在于,根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号,包括:
根据所述第一信号确定任意相邻的前一时刻和后一时刻所对应的相关峰位置和相位,将所述前一时刻和后一时刻所对应的相关峰相位进行作差处理得到处理结果;
根据所述处理结果判断所述后一时刻所对应的相关峰相位是否发生变化;
若发生变化,则对所述第一信号中的相位进行修正得到所述修正后的信号。
3.如权利要求2所述的方法,其特征在于,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差,包括:
根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离;
根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差。
4.如权利要求3所述的方法,其特征在于,根据所位置信息以及所述偏置误差补偿因子计算得到任一时刻所述直达信号的基线距离,包括:
通过如下公式计算任一时刻所述直达信号的基线距离:
L(τ0,(x,y,z))=|T(tn)-R(tn)|+ξ(x,y,z)
其中,L(τ0,(x,y,z))表示所述基线距离;τ0表示所述任一时刻;T(tn)表示所述任一时刻发射机的位置信息;R(tn)表示所述任一时刻接收机的位置信息;ξ(x,y,z)表示所述偏置误差补偿因子。
5.如权利要求4所述的方法,其特征在于,根基所述基线距离以及预设的直达信号的时延信息以及多普勒信息计算得到所述时延误差和所述相位误差,包括:
通过如下公式计算得到所述时延误差和所述相位误差:
其中,表示所述时延误差;τp(tn)表示预设的直达信号的时延信息;c表示电磁波在空间传播的速度;表示所述相位误差;表示预设的直达信号的多普勒信息;λ表示电磁波在空间传播的波长。
6.如权利要求1~5任一项所述的方法,其特征在于,根据所述时延误差和所述相位误差重构生成反射参考信号,包括:
通过如下公式重构生成反射参考信号:
其中,SR_ref(tn)表示所述反射参考信号;CB2A()表示伪随机码。
7.如权利要求6所述的方法,其特征在于,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,包括:
通过如下公式对所述第二信号进行补偿得到所述第三信号:
其中,R′R_out(tn,τ)表示所述第三信号;RR_out(tn,τ)表示所述第二信号。
8.如权利要求1~5任一项所述的方法,其特征在于,还包括:
根据预设的BP成像算法以及所述第三信号计算方位向压缩输出结果,并根据所述方位向压缩输出结果判断是否满足预设成像同步精度阈值;
若不满足,则调整所述偏置误差补偿因子重新计算出成像时反射信号压缩所对应的时延误差和相位误差,直达所述方位向压缩输出结果满足所述预设成像同步精度阈值为止。
9.一种GNSS-R系统成像预处理的装置,其特征在于,包括:
处理单元,用于确定GNSS-R系统中的直达信号以及反射信号,并根据预设的本地参考信号将所述直达信号进行压缩处理得到第一信号;
修正单元,用于根据所述第一信号确定相关峰信息以及根据所述相关峰信息对所述第一信号进行修正得到修正后的信号;
计算单元,用于确定所述GNSS-R系统中接收机和发射机的位置信息,根据所述位置信息以及预设的偏置误差补偿因子计算出成像时反射信号压缩所对应的时延误差和相位误差;
补偿单元,用于根据所述时延误差和所述相位误差重构生成反射参考信号,根据所述反射参考信号将所述反射信号进行压缩处理得到第二信号,根据所述相位误差对所述第二信号进行相位误差补偿得到第三信号,以使得基于所述第三信号进行成像。
10.如权利要求9所述的装置,其特征在于,所述修正单元,具体用于:
根据所述第一信号确定任意相邻的前一时刻和后一时刻所对应的相关峰位置和相位,将所述前一时刻和后一时刻所对应的相关峰相位进行作差处理得到处理结果;
根据所述处理结果判断所述后一时刻所对应的相关峰相位是否发生变化;
若发生变化,则对所述第一信号中的相位进行修正得到所述修正后的信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010859336.9A CN112147643B (zh) | 2020-08-24 | 2020-08-24 | 一种gnss-r系统成像预处理的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010859336.9A CN112147643B (zh) | 2020-08-24 | 2020-08-24 | 一种gnss-r系统成像预处理的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112147643A CN112147643A (zh) | 2020-12-29 |
CN112147643B true CN112147643B (zh) | 2024-02-09 |
Family
ID=73888461
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010859336.9A Active CN112147643B (zh) | 2020-08-24 | 2020-08-24 | 一种gnss-r系统成像预处理的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112147643B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116528345A (zh) * | 2022-01-21 | 2023-08-01 | 维沃移动通信有限公司 | 定时校准的方法、装置及通信设备 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1875291A (zh) * | 2003-10-28 | 2006-12-06 | 天宝导航有限公司 | 三个或更多载波的gnss信号的模糊度估计 |
CN101846746A (zh) * | 2010-03-24 | 2010-09-29 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的载波相位测高装置及其方法 |
JP2014153087A (ja) * | 2013-02-05 | 2014-08-25 | Mitsubishi Electric Corp | Gnss測位装置およびgnss測位方法 |
CN104486279A (zh) * | 2014-11-27 | 2015-04-01 | 西安空间无线电技术研究所 | 一种uqpsk微波直接调制器调制特性测试方法 |
CN105891860A (zh) * | 2016-03-25 | 2016-08-24 | 中国科学院光电研究院 | 一种基于误差分离模式的gnss区域伪距差分增强定位方法 |
WO2016145723A1 (zh) * | 2015-03-19 | 2016-09-22 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的河流流速测量方法与系统 |
CN106291610A (zh) * | 2015-06-12 | 2017-01-04 | 北京信息科技大学 | 一种用于gnss信号压缩捕获处理装置的压缩并行相关模块及其实现方法 |
CN109270507A (zh) * | 2018-11-05 | 2019-01-25 | 金华航大北斗应用技术有限公司 | Gnss-r双基sar成像机载运动误差补偿模型 |
CN109298435A (zh) * | 2017-11-02 | 2019-02-01 | 天津讯联科技有限公司 | 卫星反射信号的跟踪方法、装置及接收机 |
CN109358326A (zh) * | 2018-09-11 | 2019-02-19 | 金华航大北斗应用技术有限公司 | 基于gnss-rsar地基的数据采集成像装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8022868B2 (en) * | 2008-12-31 | 2011-09-20 | Javad Gnss, Inc. | Inter-channel bias calibration for navigation satellite systems |
US8248300B2 (en) * | 2010-06-02 | 2012-08-21 | Honeywell International Inc. | System and method of compensating for micro-jump events |
KR102065666B1 (ko) * | 2012-12-12 | 2020-02-11 | 삼성전자 주식회사 | 위성 항법 시스템의 신호 추적 방법, 신호 추적 장치 및 이를 포함하는 위성신호 수신기 |
-
2020
- 2020-08-24 CN CN202010859336.9A patent/CN112147643B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1875291A (zh) * | 2003-10-28 | 2006-12-06 | 天宝导航有限公司 | 三个或更多载波的gnss信号的模糊度估计 |
CN101846746A (zh) * | 2010-03-24 | 2010-09-29 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的载波相位测高装置及其方法 |
JP2014153087A (ja) * | 2013-02-05 | 2014-08-25 | Mitsubishi Electric Corp | Gnss測位装置およびgnss測位方法 |
CN104486279A (zh) * | 2014-11-27 | 2015-04-01 | 西安空间无线电技术研究所 | 一种uqpsk微波直接调制器调制特性测试方法 |
WO2016145723A1 (zh) * | 2015-03-19 | 2016-09-22 | 中国科学院空间科学与应用研究中心 | 一种基于gnss-r技术的河流流速测量方法与系统 |
CN106291610A (zh) * | 2015-06-12 | 2017-01-04 | 北京信息科技大学 | 一种用于gnss信号压缩捕获处理装置的压缩并行相关模块及其实现方法 |
CN105891860A (zh) * | 2016-03-25 | 2016-08-24 | 中国科学院光电研究院 | 一种基于误差分离模式的gnss区域伪距差分增强定位方法 |
CN109298435A (zh) * | 2017-11-02 | 2019-02-01 | 天津讯联科技有限公司 | 卫星反射信号的跟踪方法、装置及接收机 |
CN109358326A (zh) * | 2018-09-11 | 2019-02-19 | 金华航大北斗应用技术有限公司 | 基于gnss-rsar地基的数据采集成像装置 |
CN109270507A (zh) * | 2018-11-05 | 2019-01-25 | 金华航大北斗应用技术有限公司 | Gnss-r双基sar成像机载运动误差补偿模型 |
Non-Patent Citations (2)
Title |
---|
A Fast Algorithm for GNSS-R Reflected Signals Based on Dynamic Phase Compensation and DBZP;Yanguang Wang;《 China Satellite Navigation Conference (CSNC) 2020 Proceedings》;第242-251页 * |
GNSS-R 双基SAR 成像运动误差分析及补偿;郝令政;《导航定位学报》;第7卷(第2期);第36-43页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112147643A (zh) | 2020-12-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210215829A1 (en) | Spoofing detection in real time kinematic positioning | |
EP0965048B1 (en) | On-the-fly accuracy enhancement for civil gps receivers | |
JP5968998B2 (ja) | 測位信号検出方法、測位信号検出プログラム、測位信号受信装置、測位装置および情報機器端末 | |
US10908294B2 (en) | Detection and elimination of GNSS spoofing signals with PVT solution estimation | |
JP7377866B2 (ja) | 利用不可能なgnss測定値を再作成するための方法およびシステム | |
CA2530616A1 (en) | A hardware architecture for processing galileo alternate binary offset carrier (altboc) signals | |
WO2002063328A2 (en) | Low cost system and method for making dual band gps measurements | |
AU2013294159B2 (en) | Satellite positioning signal receiving method and device | |
CN114935770B (zh) | 一种多历元加快精密单点定位收敛速度的方法及装置 | |
CN112147643B (zh) | 一种gnss-r系统成像预处理的方法及装置 | |
NL2003860A (en) | Method and apparatus for weak data bit sync in a positioning system. | |
CN110857994A (zh) | 具有gnss伪随机噪声延迟序列生成器的gnss接收器装置 | |
CN1982913B (zh) | 定位装置、定位装置的控制方法 | |
CN105549054A (zh) | 一种基于北斗地球同步卫星的快速定位方法 | |
US20070010212A1 (en) | Apparatus for transmitting positioning signal, positioning system including the apparatus, and system for transmitting positioning signal | |
US6750814B1 (en) | Efficient algorithm for processing GPS signals | |
JP2007278708A (ja) | 衛星航法装置 | |
JP2008058146A (ja) | 測位装置、その制御方法、制御プログラム及びその記録媒体 | |
JPH06265624A (ja) | 位置計測装置 | |
CN113296137B (zh) | 干涉式形变监测方法、装置和接收机 | |
JP2001281318A (ja) | Gps受信装置における疑似雑音符号の検出方法およびその装置 | |
US20230003902A1 (en) | Method for performing a parallel search, receiver, computer program product and non-volatile storage medium | |
CN113296136B (zh) | 多目标形变监测方法、装置和接收机 | |
CN109828293B (zh) | 一种信号处理方法及其相关设备 | |
CN113296135B (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 |