CN110187347B - 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法 - Google Patents

一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法 Download PDF

Info

Publication number
CN110187347B
CN110187347B CN201910558301.9A CN201910558301A CN110187347B CN 110187347 B CN110187347 B CN 110187347B CN 201910558301 A CN201910558301 A CN 201910558301A CN 110187347 B CN110187347 B CN 110187347B
Authority
CN
China
Prior art keywords
azimuth
frequency
representing
distance
channel
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
CN201910558301.9A
Other languages
English (en)
Other versions
CN110187347A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201910558301.9A priority Critical patent/CN110187347B/zh
Publication of CN110187347A publication Critical patent/CN110187347A/zh
Application granted granted Critical
Publication of CN110187347B publication Critical patent/CN110187347B/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
    • G01S13/904—SAR modes
    • G01S13/9058—Bistatic or multistatic SAR
    • 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/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
    • G01S13/9004—SAR image acquisition techniques
    • G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth

Landscapes

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

Abstract

本发明提供了一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法，属于合成孔径雷达技术领域。本发明首先进行多通道同时进行TOPS扫描录取回波，并对接收到的多通道回波进行处理，去除欠采样模糊，解决旋转模糊，拼接得到大幅宽的成像结果。本发明针对GEO照射源对地面的大范围覆盖，在机载接收站采用TOPS模式对多个测绘带进行扫描录取回波，充分利用发射站的照射范围，扩大GEO星机双基SAR的成像幅宽。本发明特征在于对双基构型中的单基TOPS回波进行处理时，结合多通道重构技术与TOPS处理方法同时解决SAR回波频谱中的欠采样模糊和旋转模糊。

Description

一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法
技术领域
本发明属于合成孔径雷达技术领域，特别涉及一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法。
背景技术
合成孔径雷达(SAR)是一种高分辨成像雷达，具有穿透性强，能全天时、全天候的特点，目前技术已经相较成熟，被广泛应用在地球遥感、资源勘探、地形测绘等领域。
地球同步轨道合成孔径雷达(Geosynchronous Synthetic Aperture Radar,GEOSAR)的轨道高度在36000km左右的同步轨道，其波束对地面的覆盖能达到280km，有着重访周期短，可观测范围广的优势，可在大范围内保持对指定区域的持续照射。
地球同步轨道星机双基合成孔径雷达(Geosynchronous Spaceborne-AirborneBistatic SAR,GEO BiSAR)以GEO SAR的卫星作为照射源，机载平台作为接收站，是一种新的对地成像的双基SAR架构。由于其收发分置，具有双基SAR系统配置灵活的特点；并且机载接收站不主动辐射能量，隐蔽性好，发射站在地球同步轨道，生存力强，具有长期、稳定、可靠的成像能力。同时，这种“远发近收”的构型可以获得比单基GEO SAR更高的成像信噪比，并且引入高转角速度的机载平台可以在短时间内积累多普勒带宽保证成像分辨率。
然而，机载接收站受限于高度与最小天线单元面积限制，接收波束覆盖范围一般只能达到几公里级别(3km左右)，这与GEO SAR的波束覆盖范围(280km左右)相差甚远。因此，为了充分利用发射站的覆盖范围，扩大接收站的观测范围，本发明采用TOPS模式对地面进行多测绘带扫描接收。但是由于一般的GEO-BiSAR的点目标回波多普勒带宽已经超过了发射站的脉冲重复频率(PRF)，会导致频谱因欠采样而发生混叠，带来欠采样模糊，现有技术中，针对GEO BiSAR的欠采样模糊的多通道技术被提出。而在TOPS模式下接收到的回波频谱会存在多普勒质心空变，不同方位点目标的频谱占据着不同频段，导致整个场景内的多普勒带宽远远超出PRF；另外，还有针对一般TOPS回波去除旋转模糊的方法。但是这类方法都针对的都是PRF大于点目标多普勒带宽的情形，一旦PRF低于点目标多普勒带宽，方法将失效，无法进行不模糊的成像。
发明内容
本发明的目的在于解决上述问题，提出了一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法，使用GEO卫星作为照射源，在机载接收平台上沿航迹布置多个通道同时进行TOPS扫描接收回波，然后对接收到的多个通道的回波进行处理，去除欠采样模糊，再解决旋转模糊，抑制了GEO星机双基SAR在采用TOPS模式接收回波时的频谱模糊，实现了GEO星机双基SAR回波的无模糊成像。
一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法，包括以下步骤：
S1、N个通道同时对子测绘带录取回波，第m个通道的基带回波信号为Sechom(τ,η)，其中，τ和η分别表示距离向时间和方位向时间；
S2、对多通道的基带回波信号进行多通道重构，得到无欠采样模糊的回波频谱；
S3、对所述回波频谱进行去旋转模糊；
S4、对回波信号进行聚焦处理；
S5、进行去方位时域混叠，得到所述子测绘带的一次burst图像；
S6、将得到的多个子测绘带的burst图像进行拼接，获得宽测绘带SAR图像，完成大幅宽成像。
进一步地，所述步骤S1包括：
发射线性调频信号，N个通道同时对子测绘带录取回波，进行解调得到基带回波信号，第m个通道的基带回波信号为
Figure GDA0002129451410000031
其中，wr(·)和wa(·)分别表示距离向和方位向的窗函数，rect(·)表示矩形窗函数，Rm(η)表示第m个通道的双基距离和历史，
Figure GDA0002129451410000032
表示波束中心穿越目标的时刻，旋转系数
Figure GDA0002129451410000033
τ和η分别表示距离向时间和方位向时间，c表示光速，j表示虚数单位，r0表示目标点到接收站航迹的垂直距离，Td表示波束驻留时间，Tb表示一次burst的时间，Kr表示发射信号的调频率，vf表示波束脚印的速度，vR表示平台速度，ωrot表示扫描旋转角速度，λ表示中心波长，
Figure GDA0002129451410000034
c和f0分别表示光速和SAR发射脉冲载频。
进一步地，其特征在于，所述步骤S2包括：
S21、对所述步骤S1得到的每个通道的基带回波信号进行质心去斜，即相位相乘
Sderampm(τ,η)＝Sechom(τ,η)exp(-jπKdcη2)
其中，Kdc为多普勒质心随方位时间变化的斜率，
Figure GDA0002129451410000035
S22、构造第m个通道相对于参考通道的差异函数
Figure GDA0002129451410000041
其中，Δxm表示第m个通道到参考通道的距离，fη表示方位频率，k1T表示发射站到目标点的斜距一阶系数；
S23、将Hm(fη)在频域平移PRF的整数倍的结果作为元素构成系统矩阵
Figure GDA0002129451410000042
其中，N表示通道数，PRF表示脉冲重复频率；
S24、获取所述系统矩阵的逆，得到重构矩阵
Figure GDA0002129451410000043
S25、将所述步骤S21得到的Sderampm(τ,η)进行方位向快速傅里叶变换
SRDm(τ,fη)＝FFTaz{Sderampm(τ,η)}
其中，FFTaz{·}表示方位向快速傅里叶变换运算；
S26、将N个通道的基带回波信号经过所述步骤S25的结果排列成矩阵形式
SRD(τ,fη)＝[SRD1(τ,fη) SRD2(τ,fη) … SRDN(τ,fη)]
S27、将SRD(τ,fη)与所述步骤S24中的P(τ,fη)相乘，完成多通道重构，得到无欠采样模糊的回波频谱
SPRE(τ,fη)＝SRD(τ,fη)·P(τ,fη)
＝(U1(τ,fη) U2(τ,fη+PRF) … UN(τ,fη+(N-1)PRF))
其中，Uk(τ,fη),k＝1,2,…N表示重构后的第k个子带频谱，重构后的信号频谱由N个频谱拼接而成，重构后的fη由fη∈[-PRF/2,PRF/2]转换为fη∈[-N·PRF/2,N·PRF/2]，等效方位采样率为重构前的N倍，抑制欠采样模糊。
进一步地，所述步骤S3包括：
S31、对SPRE(τ,fη)进行方位向快速傅里叶变换
Sderamp(τ,η)＝IFFTaz{SPRF(τ,fη)}
S32、对Sderamp(τ,η)进行方位向快速傅里叶变换
Sp1(τ,η1)＝IFFTaz{Sderamp(τ,η)}
其中，η1表示所述步骤S32的方位向傅里叶变换后的方位时间，η1∈[-0.5N·PRF/Kdc,0.5N·PRF/Kdc]；
S33、对Sp1(τ,η1)进行方位向相位相乘
Figure GDA0002129451410000051
S34、对Sp2(τ,η1)进行方位向快速傅里叶变换
Sp3(τ,fη1)＝FFTaz{Sp2(τ,η1)}
其中，fη1表示η1对应的方位频率，fη1∈[-KdcTb/2,KdcTb/2]；
S35、对Sp3(τ,fη1)进行相位相乘
Figure GDA0002129451410000052
进一步地，所述步骤S4包括：
S41、对Secho(τ,fη1)进行距离向快速傅里叶变换
S2df(fτ,fη1)＝FFTra{Secho(τ,fη1)}
其中，FFTra{·}表示距离向快速傅里叶变换运算，fτ表示距离频率；
S42、对S2df(fτ,fη1)进行相位相乘
Figure GDA0002129451410000061
其中，RT0ref表示burst中心时刻的场景中心到发射站的斜距，RR0ref表示burst中心时刻的场景中心到接收站的斜距；
Figure GDA0002129451410000062
Rb0＝RR0ref+RT0ref，其中，k1T和k2T分别表示发射站到目标点的斜距一阶和二阶系数，θst表示接收站的初始斜视角；
S43、对所述步骤S42中相位相乘的结果在距离频率向进行Stolt插值，Stolt插值过程中距离频率与原距离频率fτ的映射关系为
Figure GDA0002129451410000063
其中，fτ'表示Stolt后的距离频率，
Figure GDA0002129451410000064
RT0表示场景内任意点到发射站的距离，RR0表示场景内任意点到接收站的距离，θe表示场景内任意点的等效斜视角，
Figure GDA0002129451410000065
S44、对所述步骤S43的结果在方位频率向进行Stolt插值，Stolt插值过程中方位频率与原方位频率fη1的映射关系为
Figure GDA0002129451410000066
其中，f′η1表示Stolt插值后的方位频率，
Figure GDA0002129451410000071
Figure GDA0002129451410000072
(xc,yc)表示场景中心点的坐标，(x,y)表示场景内任意点的坐标；
S45、对所述步骤S44的Stolt插值的结果进行距离向快速傅里叶变换
Sfocus(τ',f′η1)＝IFFTra{S'2df(f′τ,f′η1)}
其中，τ'表示fτ'对应的距离向时间，S'2df(f′τ,f′η1)表示所述步骤S44的Stolt插值的结果。
进一步地，所述步骤S5包括：
S51、对Sfocus(τ',f′η1)进行相位相乘
Figure GDA0002129451410000073
其中，K'dc表示所述步骤S45完成后的信号多普勒质心变化斜率，
Figure GDA0002129451410000074
Rr0表示接收站到场景中心的最近斜距；
S52、对Spost1(τ',f′η1)进行方位向快速傅里叶变换
Spost1(τ',η′1)＝IFFTaz{Spost1(τ',f′η1)}
S53、对Spost1(τ',η′1)进行相位相乘
Spost2(τ',η′1)＝Spost1(τ',η′1)exp(jπK'dcη′1 2)
S54、对Spost2(τ',η′1)进行方位向快速傅里叶变换
Spost3(τ',η2)＝FFTaz{Spost2(τ',η′1)}
其中，η2表示此处方位向快速傅里叶变换的方位时间，η2∈[-0.5TbKdc/K'dc,0.5TbKdc/K'dc]；
S55、对Spost3(τ',η2)进行相位相乘
Figure GDA0002129451410000081
得到子测绘带的一次burst图像。
进一步地，所述步骤S6包括：
依次对各子测绘带执行步骤S1-S5，得到各子测绘带的burst图像，将得到的多个burst图像进行拼接，获得宽测绘带SAR图像，完成大幅宽成像。
进一步地，所述步骤S1之前，包括：
初始化系统参数，包括脉冲重复频率、方位向burst采样点数、通道间隔、通道个数、扫描旋转角速度、机载平台速度及地球同步轨道卫星轨道参数。
本发明的有益效果：本发明提供了一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法，本发明根据GEO星机双基SAR的成像特点，采用在接收站配置多通道并且进行TOPS扫描多个距离测绘带的手段来录取回波，一方面可以利用多通道技术抑制回波频谱中的欠采样模糊，一方面用Two-step和Modified two-step方法处理TOPS模式下的回波，再结合多个测绘带的数据聚焦出多个测绘带的SAR图像，将其拼接后即可得到大幅宽的成像结果。由于本发明使用了多通道技术，对系统的PRF要求较低，可以大大减轻GEO SAR发射功率与数据存储的负担。而TOPS接收模式对多个距离向的测绘带进行扫描，大幅提高了机载接收站的成像范围。
附图说明
图1为本发明实施例的GEO星机双基SAR回波录取几何示意图。
图2为本发明实施例的流程图。
图3为本发明实施例的点目标分布示意图。
图4为本发明实施例的大幅宽SAR成像结果图。
图5为图4中点目标A的等高图。
图6为图4中点目标B的等高图。
图7为图4中点目标C的等高图。
具体实施方式
下面结合附图对本发明的实施例做进一步的说明。
本发明提出了一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法，地球同步轨道合成孔径雷达(Geosynchronous Synthetic Aperture Radar,GEO SAR)，为一种将SAR有效载荷放置在地球同步轨道卫星上的主动式遥感传感器。而GEO星机双基SAR(GEO-BiSAR)，即地球同步轨道星机双基合成孔径雷达，是双基地SAR的一种，是一种以地球同步轨道卫星为发射站，接收站是机载平台的星机双基地SAR。
本实施例中，GEO星机双基SAR几何结构如图1所示，系统参数如下表1所示。
表1 GEO星机双基SAR系统参数表
Figure GDA0002129451410000091
Figure GDA0002129451410000101
请参阅图2，本发明通过以下步骤实现：
S1、N个通道同时对子测绘带录取回波，第m个通道的基带回波信号为Sechom(τ,η)，其中，τ和η分别表示距离向时间和方位向时间。
本实施例中，发射线性调频信号，N＝4个通道同时对子测绘带进行TOPS扫描录取回波，进行解调得到基带回波信号，第m个通道的基带回波信号为
Figure GDA0002129451410000102
其中，m∈(1,2,3,4)，wr(·)和wa(·)分别表示距离向和方位向的窗函数，rect(·)表示矩形窗函数，Rm(η)表示第m个通道的双基距离和历史，
Figure GDA0002129451410000103
表示波束中心穿越目标的时刻，旋转系数
Figure GDA0002129451410000104
τ和η分别表示距离向时间和方位向时间，c表示光速，j表示虚数单位，r0表示目标点到接收站航迹的垂直距离，Td表示波束驻留时间，Tb表示一次burst的时间，Kr表示发射信号的调频率，vf表示波束脚印的速度，vR表示平台速度，ωrot表示扫描旋转角速度，λ表示中心波长，
Figure GDA0002129451410000111
c和f0分别表示光速和SAR发射脉冲载频。
S2、对多通道的基带回波信号进行多通道重构，得到无欠采样模糊的回波频谱。
本实施例中，步骤S2通过以下子步骤实现：
S21、对步骤S1得到的每个通道的基带回波信号进行质心去斜，即相位相乘
Sderampm(τ,η)＝Sechom(τ,η)exp(-jπKdcη2)
其中，Kdc为多普勒质心随方位时间变化的斜率，
Figure GDA0002129451410000112
S22、构造第m个通道相对于参考通道的差异函数
Figure GDA0002129451410000113
其中，Δxm表示第m个通道到参考通道的距离，fη表示方位频率，k1T表示发射站到目标点的斜距一阶估计系数。
S23、将Hm(fη)在频域平移PRF的整数倍的结果作为元素构成系统矩阵
Figure GDA0002129451410000114
其中，N表示通道数，PRF表示脉冲重复频率。
S24、获取S23中系统矩阵的逆，得到重构矩阵
Figure GDA0002129451410000115
S25、将步骤S21得到的Sderampm(τ,η)进行方位向快速傅里叶变换
SRDm(τ,fη)＝FFTaz{Sderampm(τ,η)}
其中，FFTaz{·}表示方位向快速傅里叶变换运算。
S26、将N＝4个通道的基带回波信号经过步骤S25的结果排列成矩阵形式
SRD(τ,fη)＝[SRD1(τ,fη) SRD2(τ,fη) … SRD4(τ,fη)]
S27、将S26的结果SRD(τ,fη)与步骤S24中的P(τ,fη)相乘，完成多通道重构，得到无欠采样模糊的回波频谱
SPRE(τ,fη)＝SRD(τ,fη)·P(τ,fη)
＝(U1(τ,fη) U2(τ,fη+PRF) … U4(τ,fη+3PRF))
其中，Uk(τ,fη),k＝1,2,…4表示重构后的第k个子带频谱，重构后的信号频谱由N＝4个频谱拼接而成，重构后的fη由fη∈[-PRF/2,PRF/2]转换为了fη∈[-4·PRF/2,4·PRF/2]，等效方位采样率为重构前的N＝4倍，欠采样模糊得到了有效抑制。
S3、对回波频谱进行去旋转模糊。
本实施例中，步骤S3通过以下子步骤实现：
S31、对S27的结果SPRE(τ,fη)进行方位向快速傅里叶变换
Sderamp(τ,η)＝IFFTaz{SPRF(τ,fη)}
S32、对S31的结果Sderamp(τ,η)进行方位向快速傅里叶变换
Sp1(τ,η1)＝IFFTaz{Sderamp(τ,η)}
其中，η1表示该步骤S32的方位向傅里叶变换后的方位时间，η1∈[-0.5N·PRF/Kdc,0.5N·PRF/Kdc]；
S33、对S32的结果Sp1(τ,η1)进行方位向相位相乘
Figure GDA0002129451410000121
S34、对S33的结果Sp2(τ,η1)进行方位向快速傅里叶变换
Sp3(τ,fη1)＝FFTaz{Sp2(τ,η1)}
其中，fη1表示η1对应的方位频率，fη1∈[-KdcTb/2,KdcTb/2]，Tb表示一次burst的时间。
S35、对S34的结果Sp3(τ,fη1)进行相位相乘
Figure GDA0002129451410000131
S4、对回波信号进行聚焦处理。
本实施例中，步骤S4通过以下子步骤实现：
S41、对S35中的相位相乘结果Secho(τ,fη1)进行距离向快速傅里叶变换
S2df(fτ,fη1)＝FFTra{Secho(τ,fη1)}
其中，FFTra{·}表示距离向快速傅里叶变换运算，fτ表示距离频率。
S42、对S41的结果S2df(fτ,fη1)进行相位相乘
Figure GDA0002129451410000132
其中，RT0ref表示burst中心时刻的场景中心到发射站的斜距，RR0ref表示burst中心时刻的场景中心到接收站的斜距；
Figure GDA0002129451410000133
Rb0＝RR0ref+RT0ref，其中，k1T和k2T分别表示发射站到目标点的斜距一阶和二阶估计系数，θst表示接收站的初始斜视角。
S43、对S42中相位相乘的结果在距离频率向进行Stolt插值，Stolt插值过程中新距离频率与原距离频率fτ的映射关系为
Figure GDA0002129451410000134
其中，fτ'表示Stolt后的新距离频率，
Figure GDA0002129451410000141
RT0表示场景内任意点到发射站的距离，RR0表示场景内任意点到接收站的距离，θe表示场景内任意点的等效斜视角，
Figure GDA0002129451410000142
S44、对步骤S43的结果在方位频率向进行Stolt插值，Stolt插值过程中新方位频率与原方位频率fη1的映射关系为
Figure GDA0002129451410000143
其中，f′η1表示Stolt插值后的新方位频率，
Figure GDA0002129451410000144
Figure GDA0002129451410000145
(xc,yc)表示场景中心点的坐标，(x,y)表示场景内任意点的坐标。
S45、对步骤S44的Stolt插值的结果进行距离向快速傅里叶变换
Sfocus(τ',f′η1)＝IFFTra{S'2df(f′τ,f′η1)}
其中，τ'表示fτ'对应的距离向时间，S'2df(f′τ,f′η1)表示步骤S44的Stolt插值的结果。
S5、进行去方位时域混叠，得到子测绘带的一次burst图像。
本实施例中，步骤S5通过以下子步骤实现：
S51、对S45的结果Sfocus(τ',f′η1)进行相位相乘
Figure GDA0002129451410000146
其中，K'dc表示步骤S45完成后信号的多普勒质心变化斜率，
Figure GDA0002129451410000151
Rr0表示接收站到场景中心的最近斜距；
S52、对S51的结果Spost1(τ',f′η1)进行方位向快速傅里叶变换
Spost1(τ',η′1)＝IFFTaz{Spost1(τ',f′η1)}
S53、对S52的结果Spost1(τ',η′1)进行相位相乘
Figure GDA0002129451410000152
S54、对S53的结果Spost2(τ',η′1)进行方位向快速傅里叶变换
Spost3(τ',η2)＝FFTaz{Spost2(τ',η′1)}
其中，η2表示此处方位向快速傅里叶变换的方位时间，η2∈[-0.5TbKdc/K'dc,0.5TbKdc/K'dc]。
S55、对S54的结果Spost3(τ',η2)进行相位相乘
Figure GDA0002129451410000153
得到子测绘带的一次burst图像。
S6、将得到的多个子测绘带的burst图像进行拼接，获得宽测绘带SAR图像，完成大幅宽成像。
本实施例中，依次对各子测绘带执行步骤S1-S5，得到各子测绘带的burst图像，将得到的多个burst图像进行拼接，获得宽测绘带SAR图像，完成大幅宽成像。
图3为点目标分布示意图。如图4-7所示为本发明实施例中点目标的结果图。图4为3个测绘带拼接而成的大幅宽SAR图像；图5为图4中点目标A的等高图；图6为图4中点目标B的等高图；图7为图4中点目标C的等高图。
本领域的普通技术人员将会意识到，这里的实施例是为了帮助读者理解本发明的原理，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (7)

1.一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法，其特征在于，包括以下步骤：
S1、N个通道同时对子测绘带录取回波，第m个通道的基带回波信号为Sechom(τ,η)，其中，τ和η分别表示距离向时间和方位向时间；
S2、对多通道的基带回波信号进行多通道重构，得到无欠采样模糊的回波频谱；所述步骤S2包括：
S21、对所述步骤S1得到的每个通道的基带回波信号进行质心去斜，即相位相乘
Figure FDA0002755253710000011
其中，Kdc为多普勒质心随方位时间变化的斜率，
Figure FDA0002755253710000012
vR表示平台速度，ωrot表示扫描旋转角速度，λ表示中心波长；
S22、构造第m个通道相对于参考通道的差异函数
Figure FDA0002755253710000013
其中，Δxm表示第m个通道到参考通道的距离，fη表示方位频率，k1T表示发射站到目标点的斜距一阶系数；
S23、将Hm(fη)在频域平移PRF的整数倍的结果作为元素构成系统矩阵
Figure FDA0002755253710000014
其中，N表示通道数，PRF表示脉冲重复频率；
S24、获取所述系统矩阵的逆，得到重构矩阵
Figure FDA0002755253710000021
S25、将所述步骤S21得到的Sderampm(τ,η)进行方位向快速傅里叶变换
SRDm(τ,fη)＝FFTaz{Sderampm(τ,η)}
其中，FFTaz{·}表示方位向快速傅里叶变换运算；
S26、将N个通道的基带回波信号经过所述步骤S25的结果排列成矩阵形式
SRD(τ,fη)＝[SRD1(τ,fη) SRD2(τ,fη)…SRDN(τ,fη)]
S27、将SRD(τ,fη)与所述步骤S24中的P(fη)相乘，完成多通道重构，得到无欠采样模糊的回波频谱
SPRE(τ,fη)＝SRD(τ,fη)·P(fη)
＝(U1(τ,fη) U2(τ,fη+PRF)…UN(τ,fη+(N-1)PRF))
其中，Uk(τ,fη),k＝1,2,…N表示重构后的第k个子带频谱，重构后的信号频谱由N个频谱拼接而成，重构后的fη由fη∈[-PRF/2,PRF/2]转换为fη∈[-N·PRF/2,N·PRF/2]，等效方位采样率为重构前的N倍，抑制欠采样模糊；
S3、对所述回波频谱进行去旋转模糊；
S4、对回波信号进行聚焦处理；
S5、进行去方位时域混叠，得到所述子测绘带的一次burst图像；
S6、将得到的多个子测绘带的burst图像进行拼接，获得宽测绘带SAR图像，完成大幅宽成像。
2.如权利要求1所述的地球同步轨道星机双基合成孔径雷达大幅宽成像方法，其特征在于，所述步骤S1包括：
发射线性调频信号，N个通道同时对子测绘带录取回波，进行解调得到基带回波信号，第m个通道的基带回波信号为
Figure FDA0002755253710000031
其中，wr(·)和wa(·)分别表示距离向和方位向的窗函数，rect(·)表示矩形窗函数，Rm(η)表示第m个通道的双基距离和历史，
Figure FDA0002755253710000032
表示波束中心穿越目标的时刻，旋转系数
Figure FDA0002755253710000033
τ和η分别表示距离向时间和方位向时间，c表示光速，j表示虚数单位，r0表示目标点到接收站航迹的垂直距离，Td表示波束驻留时间，Tb表示一次burst的时间，Kr表示发射信号的调频率，vf表示波束脚印的速度，
Figure FDA0002755253710000034
c和f0分别表示光速和SAR发射脉冲载频。
3.如权利要求2所述的地球同步轨道星机双基合成孔径雷达大幅宽成像方法，其特征在于，所述步骤S3包括：
S31、对SPRE(τ,fη)进行方位向快速傅里叶逆变换
Sderamp(τ,η)＝IFFTaz{SPRF(τ,fη)}
S32、对Sderamp(τ,η)进行方位向快速傅里叶逆变换
Sp1(τ,η1)＝IFFTaz{Sderamp(τ,η)}
其中，η1表示所述步骤S32的方位向傅里叶逆变换后的方位时间，η1∈[-0.5N·PRF/Kdc,0.5N·PRF/Kdc]；
S33、对Sp1(τ,η1)进行方位向相位相乘
Figure FDA0002755253710000035
S34、对Sp2(τ,η1)进行方位向快速傅里叶变换
Sp3(τ,fη1)＝FFTaz{Sp2(τ,η1)}
其中，fη1表示η1对应的方位频率，fη1∈[-KdcTb/2,KdcTb/2]；
S35、对Sp3(τ,fη1)进行相位相乘
Figure FDA0002755253710000041
4.如权利要求3所述的地球同步轨道星机双基合成孔径雷达大幅宽成像方法，其特征在于，所述步骤S4包括：
S41、对Secho(τ,fη1)进行距离向快速傅里叶变换
S2df(fτ,fη1)＝FFTra{Secho(τ,fη1)}
其中，FFTra{·}表示距离向快速傅里叶变换运算，fτ表示距离频率；
S42、对S2df(fτ,fη1)进行相位相乘
Figure FDA0002755253710000042
其中，RT0ref表示burst中心时刻的场景中心到发射站的斜距，RR0ref表示burst中心时刻的场景中心到接收站的斜距；
Figure FDA0002755253710000043
Rb0＝RR0ref+RT0ref，其中，k1T和k2T分别表示发射站到目标点的斜距一阶和二阶系数，θst表示接收站的初始斜视角；
S43、对所述步骤S42中相位相乘的结果在距离频率向进行Stolt插值，Stolt插值过程中距离频率与原距离频率fτ的映射关系为
Figure FDA0002755253710000051
其中，f′τ表示Stolt后的距离频率，
Figure FDA0002755253710000052
RT0表示场景内任意点到发射站的距离，RR0表示场景内任意点到接收站的距离，θe表示场景内任意点的等效斜视角，
Figure FDA0002755253710000053
S44、对所述步骤S43的结果在方位频率向进行Stolt插值，Stolt插值过程中方位频率与原方位频率fη1的映射关系为
Figure FDA0002755253710000054
其中，f′η1表示Stolt插值后的方位频率，
Figure FDA0002755253710000055
Figure FDA0002755253710000056
(xc,yc)表示场景中心点的坐标，(x,y)表示场景内任意点的坐标；
S45、对所述步骤S44的Stolt插值的结果进行距离向快速傅里叶逆变换
Sfocus(τ',f′η1)＝IFFTra{S′2df(f′τ,f′η1)}
其中，τ'表示fτ'对应的距离向时间，S′2df(f′τ,f′η1)表示所述步骤S44的Stolt插值的结果。
5.如权利要求4所述的地球同步轨道星机双基合成孔径雷达大幅宽成像方法，其特征在于，所述步骤S5包括：
S51、对Sfocus(τ',f′η1)进行相位相乘
Figure FDA0002755253710000061
其中，K′dc表示所述步骤S45完成后的信号多普勒质心变化斜率，
Figure FDA0002755253710000062
Rr0表示接收站到场景中心的最近斜距；
S52、对Spost1(τ',f′η1)进行方位向快速傅里叶逆变换
Spost1(τ',η′1)＝IFFTaz{Spost1(τ',f′η1)}
S53、对Spost1(τ',η′1)进行相位相乘
Spost2(τ',η′1)＝Spost1(τ',η′1)exp(jπK′dcη′1 2)
S54、对Spost2(τ',η′1)进行方位向快速傅里叶变换
Spost3(τ',η2)＝FFTaz{Spost2(τ',η′1)}
其中，η2表示此处方位向快速傅里叶变换的方位时间，η2∈[-0.5TbKdc/K′dc,0.5TbKdc/K′dc]；
S55、对Spost3(τ',η2)进行相位相乘
Figure FDA0002755253710000063
得到子测绘带的一次burst图像。
6.如权利要求5所述的地球同步轨道星机双基合成孔径雷达大幅宽成像方法，其特征在于，所述步骤S6包括：
依次对各子测绘带执行步骤S1-S5，得到各子测绘带的burst图像，将得到的多个burst图像进行拼接，获得宽测绘带SAR图像，完成大幅宽成像。
7.如权利要求1-6任一项所述的地球同步轨道星机双基合成孔径雷达大幅宽成像方法，其特征在于，所述步骤S1之前，包括：
初始化系统参数，包括脉冲重复频率、方位向burst采样点数、通道间隔、通道个数、扫描旋转角速度、机载平台速度及地球同步轨道卫星轨道参数。
CN201910558301.9A 2019-06-26 2019-06-26 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法 Active CN110187347B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910558301.9A CN110187347B (zh) 2019-06-26 2019-06-26 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910558301.9A CN110187347B (zh) 2019-06-26 2019-06-26 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法

Publications (2)

Publication Number Publication Date
CN110187347A CN110187347A (zh) 2019-08-30
CN110187347B true CN110187347B (zh) 2021-01-12

Family

ID=67723413

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910558301.9A Active CN110187347B (zh) 2019-06-26 2019-06-26 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法

Country Status (1)

Country Link
CN (1) CN110187347B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111665506B (zh) * 2020-06-22 2022-07-05 西安空间无线电技术研究所 基于去斜处理的高分辨率sar变重频均匀化重采样方法
CN111880181B (zh) * 2020-07-31 2021-11-30 北京环境特性研究所 一种微弱痕迹处理方法、装置、检测方法及装置
CN112379379A (zh) * 2020-11-10 2021-02-19 河南大学 Tops sar系统的通道相位偏差估计方法和装置
CN113238229B (zh) * 2021-05-25 2022-06-24 电子科技大学 一种geo星机双基sar无模糊成像方法
CN113933805B (zh) * 2021-12-18 2022-03-22 中国科学院空天信息创新研究院 一种星机混合模式成像同步误差估计方法
CN114325705B (zh) * 2021-12-31 2023-09-22 电子科技大学 一种高低轨双基地合成孔径雷达频域快速成像方法
CN114998176B (zh) * 2022-06-24 2024-08-09 北京市遥感信息研究所 一种大幅宽卫星遥感影像处理方法及装置
CN117169887B (zh) * 2023-11-03 2024-04-19 武汉能钠智能装备技术股份有限公司 基于方向确定的sar地面移动目标定位方法
CN118068330B (zh) * 2024-04-19 2024-06-18 中国科学院空天信息创新研究院 一种基于通道对消的多通道sar系统方位模糊抑制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135100A (zh) * 2013-01-31 2013-06-05 西安电子科技大学 同轨双基sar的动目标参数估计方法
CN104808204A (zh) * 2015-04-13 2015-07-29 电子科技大学 一种固定发射站双基前视sar动目标检测方法与成像方法
CN104898120A (zh) * 2015-06-12 2015-09-09 西安电子科技大学 一种基于回波模拟的双基前视高机动平台sar成像方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7701380B2 (en) * 2007-03-07 2010-04-20 Chirp Corporation Beam phase modulation for improved synthetic aperture detection and estimation
US9335409B2 (en) * 2013-03-20 2016-05-10 Raytheon Company Bistatic inverse synthetic aperture radar imaging
DE102013214676A1 (de) * 2013-07-26 2015-01-29 Deutsches Zentrum für Luft- und Raumfahrt e.V. Verfahren zur Korrektur der Zeit- und Phasenreferenzen von nicht-synchronen SAR-Daten
CN103454632B (zh) * 2013-08-24 2015-05-20 西安电子科技大学 一站固定式调频连续波双基地sar成像方法
CN103412310B (zh) * 2013-08-26 2015-07-15 电子科技大学 双基地前视合成孔径雷达地面动目标检测方法与成像方法
CN108039980B (zh) * 2018-02-07 2019-10-25 电子科技大学 Geo星-机多基sar拓扑结构设计方法
CN109001733B (zh) * 2018-06-15 2020-04-28 电子科技大学 一种geo星机双基合成孔径雷达频谱模糊抑制方法
CN109164449A (zh) * 2018-09-20 2019-01-08 北京空间飞行器总体设计部 一种高低轨双基sar斜距确定方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135100A (zh) * 2013-01-31 2013-06-05 西安电子科技大学 同轨双基sar的动目标参数估计方法
CN104808204A (zh) * 2015-04-13 2015-07-29 电子科技大学 一种固定发射站双基前视sar动目标检测方法与成像方法
CN104898120A (zh) * 2015-06-12 2015-09-09 西安电子科技大学 一种基于回波模拟的双基前视高机动平台sar成像方法

Also Published As

Publication number Publication date
CN110187347A (zh) 2019-08-30

Similar Documents

Publication Publication Date Title
CN110187347B (zh) 一种地球同步轨道星机双基合成孔径雷达大幅宽成像方法
Li et al. A frequency-domain imaging algorithm for highly squinted SAR mounted on maneuvering platforms with nonlinear trajectory
Pastina et al. Multistatic and MIMO distributed ISAR for enhanced cross-range resolution of rotating targets
CN103744068B (zh) 双通道调频连续波sar系统的动目标检测成像方法
CN105445704B (zh) 一种sar图像中的雷达动目标抑制方法
CN106932778B (zh) 方位多通道fmcw sar滑动聚束成像方法
CN114910911B (zh) 一种基于多相位中心重构的星载多基sar成像方法
CN102288964A (zh) 一种星载高分辨率合成孔径雷达的成像处理方法
CN114545411B (zh) 一种基于工程实现的极坐标格式多模高分辨sar成像方法
Fang et al. Improved passive SAR imaging with DVB-T transmissions
CN111505639A (zh) 基于变重频采样模式的合成孔径雷达宽幅稀疏成像方法
CN104898107A (zh) 一种多发多收合成孔径激光雷达信号处理方法
CN107102328B (zh) 基于fpga的实时成像信号处理方法及fpga
Kraus et al. Addressing the terrain topography in distributed SAR imaging
CN105182335B (zh) 基于奇异值分解的地球同步轨道sar成像方法
Dogan et al. Experimental demonstration of a novel end-to-end SAR range ambiguity suppression method
Bezvesilniy et al. Synthetic aperture radar systems for small aircrafts: Data processing approaches
CN110308445B (zh) 一种基于车载数字阵调频连续波雷达的成像方法
US20230213648A1 (en) Radar processing method and system
CN115407335A (zh) 基于OFDM-chirp信号和DBF处理的MIMO-SAR解模糊方法及装置
CN104076361B (zh) 无人机载战场监视雷达超分辨广域成像方法
Li et al. High squint multichannel SAR imaging algorithm for high speed maneuvering platforms with small-aperture
Rommel Development, implementation, and analysis of a multiple-input multiple-output concept for spaceborne high-resolution wide-swath synthetic aperture radar
Kim SAR video generation of MIMO video SAR with beat frequency division FMCW
Hu et al. A signal reconstruction algorithm for distributed multichannel squint missle-borne 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