CN111157991A - 一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法 - Google Patents
一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法 Download PDFInfo
- Publication number
- CN111157991A CN111157991A CN201911424889.5A CN201911424889A CN111157991A CN 111157991 A CN111157991 A CN 111157991A CN 201911424889 A CN201911424889 A CN 201911424889A CN 111157991 A CN111157991 A CN 111157991A
- Authority
- CN
- China
- Prior art keywords
- chebyshev polynomial
- dimensional
- distance
- sar
- 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.)
- Pending
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 39
- 238000001228 spectrum Methods 0.000 claims abstract description 43
- 230000014509 gene expression Effects 0.000 claims abstract description 39
- 230000008878 coupling Effects 0.000 claims abstract description 22
- 238000010168 coupling process Methods 0.000 claims abstract description 22
- 238000005859 coupling reaction Methods 0.000 claims abstract description 22
- 238000007906 compression Methods 0.000 claims abstract description 17
- 238000013508 migration Methods 0.000 claims abstract description 15
- 230000005012 migration Effects 0.000 claims abstract description 15
- 230000006835 compression Effects 0.000 claims abstract description 14
- 238000012937 correction Methods 0.000 claims abstract description 12
- 230000001131 transforming effect Effects 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 claims description 26
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 230000003595 spectral effect Effects 0.000 claims description 4
- 238000006467 substitution reaction Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 21
- 230000000694 effects Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 5
- 238000004088 simulation Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 208000004350 Strabismus Diseases 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000007983 Tris buffer Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
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/9004—SAR image acquisition techniques
- G01S13/9011—SAR image acquisition techniques with frequency domain processing of the SAR signals in azimuth
-
- 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/9019—Auto-focussing of the SAR signals
-
- 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/9021—SAR image post-processing techniques
- G01S13/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
-
- 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/9089—SAR having an irregular aperture
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
本发明公开了一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega‑k成像方法,包括:依据曲线SAR***的几何模型和运动方程,建立曲线运动轨迹模式下的斜距表达式;利用Chebyshev多项式对所述斜距表达式进行四阶近似;将曲线SAR回波变换到二维频域得到二维频谱;再利用Chebyshev多项式对二维频谱的耦合相位项进行三阶展开;最后依据基于Chebyshev多项式的二维频谱进行相位补偿,并完成距离徙动校正、距离向聚焦和方位压缩,变换到二维时域得到聚焦图像。本发明能够解决曲线运动轨迹SAR的成像问题并提高其成像性能。
Description
技术领域
本发明涉及雷达信号处理领域,特别是指一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法。
背景技术
随着日益复杂的外界环境和多样化的应用需求,曲线运动轨迹SAR成像的研究具有特殊的应用价值。曲线运动轨迹SAR克服了机载直线SAR***因成像特性的限制而存在的前视盲区,能灵活地选择成像区域。但复杂的距离历程给曲线SAR成像处理带来了困难。
已有文献提出了曲线运动轨迹SAR的信号处理方法。例如,Tang Shi-yang,LiuJun-bin,An Hong-yang等在研究星载SAR时,用二阶或三阶Taylor级数来近似曲线SAR的斜距历程,但较大的近似误差影响了回波二维频谱的精确性。M.Bao,Xing Meng-dao,LiZhen-yu分别提出了改进的斜距模型、改进的双曲斜距方程,这些斜距的数学表达式比较简单,但准确度不够。Wang Peng-bo,Liu Wei等提出改进的等效斜视距离模型,只考虑了两维加速度,不能用来表达具有三维加速度的曲线运动轨迹SAR的斜距历程。上述文献中,成像算法处理的核心是将二维频谱关于距离向频率作泰勒展开,由此产生的相位误差与泰勒展开的阶数密切相关,并且在合成孔径边缘,相位误差快速增长。使曲线SAR的成像性能下降。
发明内容
本发明的主要目的在于克服现有技术中的上述缺陷,提出一种能够解决提高场景目标成像精度问题的基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法。
基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法,包括如下步骤:
依据曲线SAR***的几何模型和运动方程建立曲线运动轨迹模式下的斜距表达式;
对所述斜距表达式用Chebyshev多项式进行四阶近似,并按方位慢时间的幂级数排列,得到基于Chebyshev多项式的斜距表达式;
将所述基于Chebyshev多项式的斜距表达式对应的曲线SAR回波信号变换到二维频域得到二维频谱;
将所述二维频谱的相位分为距离和方位相互耦合相位项和非耦合相位项,并利用Chebyshev多项式对耦合相位项进行三阶近似,得到基于Chebyshev多项式的二维频谱;
依据所述基于Chebyshev多项式的二维频谱进行相位补偿,并完成距离徙动校正、距离向聚焦和方位压缩,然后变换到二维时域得到聚焦图像。
优选的,依据曲线SAR***的几何模型和运动方程建立曲线运动轨迹模式下的斜距表达式,所述斜距表达式如下:
其中,Vx,Vy,Vz是平台的三维速度,Ax,Ay,Az是三维加速度,(Xs,Ys,0)是目标位置,Hs是雷达的初始飞行高度,tv是方位向慢时间。
优选的,对所述斜距表达式用Chebyshev多项式进行四阶近似,并按方位慢时间的幂级数排列,具体包括:
将方位慢时间进行归一化处理,如下:
其中,Tsyn是合成孔径时间,
将x代入所述斜距表达式,对所述斜距表达式用Chebyshev多项式进行四阶近似后,得到的斜距表达式为:
其中,Ti(x)是Chebyshev多项式,αci是Chebyshev系数,n=4,
优选的,将所述基于Chebyshev多项式的斜距表达式对应的曲线SAR回波信号变换到二维频域得到二维频谱,所述二维频谱表示如下:
二维频谱的相位表示如下:
其中,χc0=2βc0,χc1=2βc1,χc2=2βc2,χc3=2βc3,χc4=2βc4,Kr是发射信号的调频斜率,c是光速,fc为载波频率。
优选的,将所述二维频谱的相位分为距离和方位相互耦合相位项和非耦合相位项,并利用Chebyshev多项式对耦合相位进行三阶近似,具体包括:
将二维频谱的相位分成关于距离和方位相互耦合相位项和非耦合相位项,所述耦合相位项表示如下:
其中,Br是发射信号的带宽。
利用Chebyshev多项式对耦合相位关于y进行三阶展开,得到:
其中,切比雪夫多项式分别为T0(y)=1,T1(y)=y,T2(y)=2y2-1,T3(y)=4y3-3y,T4(y)=8y4-8y2+1,δci是切比雪夫系数且有
此处第一项是固定相位项;第二项只与fν有关,对应于方位压缩;第三项是fr的一次项,对应于距离徙动;第四项是fr 2的二次项,对应于距离聚焦;第五项是距离和方位的高阶耦合项。
优选的,依据所述基于Chebyshev多项式的二维频谱进行相位补偿,并完成距离徙动校正、距离向聚焦和方位压缩,然后变换到二维时域得到聚焦图像,具体包括:
(2)在二维频域中,将完成了固定相位补偿的信号与Hdecou,Che(fr,fv)相乘进行高阶相位补偿,Hdecou,Che(fr,fv)具体如下:
(3)在二维频域中,将完成了高阶相位补偿的信号与Hrcm,Che(fr,fv)相乘进行距离徙动校正,Hrcm,Che(fr,fv)具体如下:
(5)在二维频域中,将所述完成了高阶相位补偿和距离徙动校正的信号与Hrc+src,Che(fr,fv)相乘进行距离压缩和二次距离压缩,实现距离向聚焦,Hrc+src,Che(fr,fv)具体如下:
(6)在二维频域中,将完成距离向聚焦处理的信号与Haz,Che(fr,fv)相乘进行方位压缩处理,Haz,Che(fr,fv)具体如下:
(7)将方位压缩后的信号变换到二维时域,得到聚焦图像。
由上述对本发明的描述可知,与现有技术相比,本发明具有如下有益效果:
(1)利用Chebyshev多项式对曲线轨迹模式下的斜距表达式进行四阶近似,比传统的泰勒近似精度高,误差小;
(2)在二维频域,再次利用Chebyshev多项式对耦合相位项进行分解,减小了解耦合的误差,提高了相位精度;
(3)Chebyshev逼近函数的数值方法,效率高,适用于工程实现;
(4)采用Chebyshev多项式对斜距和二维频谱近似,减小相位误差的同时提高了相位补偿函数的精度,从而改善了曲线SAR的成像性能。
附图说明
图1是本发明的扩展的Omega-k成像算法流程图;
图2是本发明实施提供的曲线运动轨迹SAR***的几何模型图;
图3是本发明实施提供的成像区域多点目标等高线图和亮度图;
图4是本发明实施提供的场景中心点目标P0成像效果图;
图中:(a)本发明算法;(b)基于Taylor分解的二维频域算法;(c)距离向冲激响应比较;(d)方位向冲激响应比较。
图5是本发明实施提供的场景边缘点目标P2成像效果图;
图中:(a)本发明算法;(b)基于Taylor分解的二维频域算法;(c)距离向冲激响应比较;(d)方位向冲激响应比较。
具体实施方式
以下通过具体实施方式对本发明作进一步的描述。
如图1,本发明一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k的成像方法,具体实施步骤如下:
S101,机载SAR平台由于加速度的存在而产生曲线运动,几何模型如图2所示,初始时刻位于点Q(0,0,Hs)处的SAR载体以初始速度(Vx,Vy,Vz),加速度(Ax,Ay,Az)沿Y轴方向飞行,tv时刻到达T(xn,yn,zn)点。假设观测场景内点目标P坐标为(Xs,Ys,0)。依据运动方程,瞬时斜距表达式为:
S102,用Chebyshev多项式对瞬时斜距进行四阶近似,得到曲线运动轨迹模式下的斜距表达式。具体操作如下:
Rc(x)=αc0T0(x)+αc1T1(x)+αc2T2(x)+αc3T3(x)
其中,切比雪夫多项式T0(x)=1,T1(x)=x,T2(x)=2x2-1,T3(x)=4x3-3xT4(x)=8x4-8x2+1,αci是切比雪夫系数,其计算公式如下:
S103,将斜距对应的曲线SAR回波变换到二维频域得到二维频谱。
具体的,SAR发射线性调频信号,包括散射系数在内的总增益为1,解调后回波信号表达式为:
其中,tr是距离向时间,wr(·)和wv(·)分别是距离向和方位向的时域包络,Kr是发射信号的调频斜率,λ是发射信号的波长,c是光速。
根据驻留相位原理,将点目标的回波信号进行距离向傅里叶变换,得到的表达式为:
其中,Wr(·)为距离向的谱包络,fc为载波频率,fr为距离向频率变量。
将所述距离频域的回波信号进行方位向傅里叶变换,根据级数反演法,保留四次相位得到二维频谱表达式为
其中,fν为方位向频率,Wv(·)为方位向谱包络,χc0=2βc0,χc1=2βc1,χc2=2βc2,χc3=2βc3,χc4=2βc4。
S104,利用Chebyshev多项式展开二维频谱
由于回波的二维频谱存在距离和方位的耦合,需要进行解耦合。因此首先依据二维频谱与fr和fv的关系,将其分成耦合项和非耦合项,具体如下:
在SAR***中满足fc□fr,Chebyshev多项式具有正交性、独立性、收敛性等特点,是连续函数的多项式形式的最佳一致逼近。利用切比雪夫多项式来描述SAR***的斜距和二维频谱,可以提高目标的成像性能,并减少数学和计算的复杂性。具体操作步骤如下:
其中,Br是发射信号的带宽。
第二步,将上式按照Chebyshev多项式进行三阶展开,得到耦合相位项的表达式为:
其中,切比雪夫多项式分别为T0(y)=1,T1(y)=y,T2(y)=2y2-1,T3(y)=4y3-3y,T4(y)=8y4-8y2+1,δci是切比雪夫系数且有
此处第一项是固定相位项;第二项只与fv有关,对应于方位压缩;第三项是fr的一次项,对应于距离徙动;第四项是fr 2的二次项,对应于距离聚焦;第五项是距离和方位的高阶耦合项。
S105,依据基于Chebyshev多项式的二维频谱进行相位补偿,距离徙动校正和距离向聚焦处理。具体包括:
(3)在二维频域中,将所述完成了固定相位补偿的信号与Hdecou,Che(fr,fv)相乘进行高阶相位补偿,Hdecou,Che(fr,fv)具体如下:
(4)在二维频域中,将所述完成了相位补偿的信号与Hrcm,Che(fr,fv)相乘进行距离徙动校正,Hrcm,Che(fr,fv)具体如下:
(5)在二维频域中,将所述完成了相位补偿和距离徙动校正的信号与Hrc+src,Che(fr,fv)相乘进行距离压缩和二次距离压缩,实现距离向聚焦,Hrc+src,Che(fr,fv)具体如下:
S106,在二维频域中,将完成距离向聚焦处理的信号与Haz,Che(fr,fv)相乘进行方位压缩处理,Haz,Che(fr,fv)具体如下:
S107,将方位压缩后的信号变换到二维时域,得到聚焦图像。
如图2是本发明的曲线运动轨迹SAR***的几何模型图。
以下将通过仿真实验对本发明方法加以说明。
仿真参数如表1所示。成像场景中设置了5*5的点阵目标。每个点目标间距是100m,中心点目标是P0,另外选取点目标P1,P2,P3,P4做成像质量分析。
表1扩展的Omega-k算法仿真参数
图3(a)是采用本发明方法对场景区域多点目标成像的高线图和亮度图,25个点目标清晰地排列,说明本文算法对具有三维速度和三维加速度的曲线运动轨迹SAR能有效成像。图3(b)是采用基于Taylor分解的二维频域算法对同样参数下的曲线SAR成像的高线图和亮度图,从整体效果看,两种算法的亮度图差别微小,几乎无法分辨。但本发明算法得到的高线图更加清晰,成像质量更好。
图4和图5分别是场景中心点目标P0和边缘点目标P2的细节效果图。由仿真图可知,场景中心点P0和边缘点P2的成像效果,几乎一致,这说明成像算法对场景各点的聚焦性能是一致的。由图4(a)和图5(a)可知,基于本发明算法得到的点目标的高线图,二维主旁瓣非常清晰,不存在耦合现象;而基于Taylor分解的成像算法对点目标成像得到的高线图,距离向聚焦效果较差,出现散焦和主旁瓣不对称现象,如图4(b)和图5(b)所示。此外,从图4(d)和图5(d)看到,所有目标的方位冲激响应在两种算法下几乎重合。但是采用对比算法得到的距离向冲激响应与本文算法相比,中心点目标的旁瓣明显升高,而边缘点在旁瓣升高的同时,稍微往左偏移。聚焦性能明显下降,如图4(c)和图5(c)所示。
为了进一步定量分析本发明算法的有效性,采用峰值旁瓣比(Peak side loberatio,PSLR)和积分旁瓣比(Integral side lobe ratio,ISLR)来评价本发明算法对点目标P0,P1,P2,P3,P4的成像质量,并与基于Taylor分解的二维频域成像算法比较,如表2所示。本发明方法所得到的性能参数更接近理论值。
表2本发明算法成像质量评价指标
上述仅为本发明的具体实施方式,但本发明的设计构思并不局限于此,凡利用此构思对本发明进行非实质性的改动,均应属于侵犯本发明保护范围的行为。
Claims (6)
1.一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法,其特征在于,包括如下步骤:
依据曲线SAR***的几何模型和运动方程建立曲线运动轨迹模式下的斜距表达式;
对所述斜距表达式用Chebyshev多项式进行四阶近似,并按方位慢时间的幂级数排列,得到基于Chebyshev多项式的斜距表达式;
将所述基于Chebyshev多项式的斜距表达式对应的曲线SAR回波信号变换到二维频域得到二维频谱;
将所述二维频谱的相位分为距离和方位相互耦合相位项和非耦合相位项,并利用Chebyshev多项式对耦合相位项进行三阶近似,得到基于Chebyshev多项式的二维频谱;
依据所述基于Chebyshev多项式的二维频谱进行相位补偿,并完成距离徙动校正、距离向聚焦和方位压缩,然后变换到二维时域得到聚焦图像。
3.根据权利要求2所述的基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法,其特征在于,对所述斜距表达式用Chebyshev多项式进行四阶近似,并按方位慢时间的幂级数排列,具体包括:
将方位慢时间进行归一化处理,如下:
其中,Tsyn是合成孔径时间,将x代入所述斜距表达式,对所述斜距表达式用Chebyshev多项式进行四阶近似后,得到的斜距表达式为:
其中,Ti(x)是Chebyshev多项式,αci是Chebyshev系数,n=4,
5.根据权利要求4所述的基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法,其特征在于,将所述二维频谱的相位分为距离和方位相互耦合相位项和非耦合相位项,并利用Chebyshev多项式对耦合相位进行三阶近似,具体包括:
将二维频谱的相位分成关于距离和方位相互耦合相位项和非耦合相位项,所述耦合相位项表示如下:
其中,Br是发射信号的带宽;
利用Chebyshev多项式对耦合相位关于y进行三阶展开,得到:
其中,切比雪夫多项式分别为T0(y)=1,T1(y)=y,T2(y)=2y2-1,T3(y)=4y3-3y,T4(y)=8y4-8y2+1,δci是切比雪夫系数且有
此处第一项是固定相位项;第二项只与fν有关,对应于方位压缩;第三项是fr的一次项,对应于距离徙动;第四项是fr 2的二次项,对应于距离聚焦;第五项是距离和方位的高阶耦合项。
6.根据权利要求5所述的基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法,其特征在于,依据所述基于Chebyshev多项式的二维频谱进行相位补偿,并完成距离徙动校正、距离向聚焦和方位压缩,然后变换到二维时域得到聚焦图像,具体包括:
Hfp,Che=exp(-jθfp)
(2)在二维频域中,将完成了固定相位补偿的信号与Hdecou,Che(fr,fv)相乘进行高阶相位补偿,Hdecou,Che(fr,fv)具体如下:
Hdecou,Che(fr,fν)=exp[-jθcou(fν)fr 3]
(3)在二维频域中,将完成了高阶相位补偿的信号与Hrcm,Che(fr,fv)相乘进行距离徙动校正,Hrcm,Che(fr,fv)具体如下:
Hrcm,Che(fr,fv)=exp[-jθrcm(fv)fr]
(5)在二维频域中,将所述完成了高阶相位补偿和距离徙动校正的信号与Hrc+src,Che(fr,fv)相乘进行距离压缩和二次距离压缩,实现距离向聚焦,Hrc+src,Che(fr,fv)具体如下:
(6)在二维频域中,将完成距离向聚焦处理的信号与Haz,Che(fr,fv)相乘进行方位压缩处理,Haz,Che(fr,fv)具体如下:
Haz,Che(fr,fv)=exp[-jθaz(fv)]
(7)将方位压缩后的信号变换到二维时域,得到聚焦图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911424889.5A CN111157991A (zh) | 2019-12-31 | 2019-12-31 | 一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911424889.5A CN111157991A (zh) | 2019-12-31 | 2019-12-31 | 一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111157991A true CN111157991A (zh) | 2020-05-15 |
Family
ID=70560761
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911424889.5A Pending CN111157991A (zh) | 2019-12-31 | 2019-12-31 | 一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111157991A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6178141B1 (en) * | 1996-11-20 | 2001-01-23 | Gte Internetworking Incorporated | Acoustic counter-sniper system |
US20070143734A1 (en) * | 2005-12-16 | 2007-06-21 | Uzhakov Sergey V | Method and system for improving aerial image simulation speeds |
US8094907B1 (en) * | 2007-05-03 | 2012-01-10 | University Of Virginia | System, method and computer program product for fast conjugate phase reconstruction based on polynomial approximation |
CN106249237A (zh) * | 2016-07-19 | 2016-12-21 | 西安电子科技大学 | 一种曲线轨迹下大斜视sar频域成像方法 |
US20170131397A1 (en) * | 2014-06-27 | 2017-05-11 | Robert Bosch Gmbh | Method for locating an object using an fmcw radar |
CN109557542A (zh) * | 2018-11-28 | 2019-04-02 | 江苏天域航空科技集团股份有限公司 | 俯冲模式双基地前视sar成像方法 |
CN110208799A (zh) * | 2019-06-19 | 2019-09-06 | 华侨大学 | 基于Legendre正交分解的曲线运动轨迹SAR波数域成像方法 |
-
2019
- 2019-12-31 CN CN201911424889.5A patent/CN111157991A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6178141B1 (en) * | 1996-11-20 | 2001-01-23 | Gte Internetworking Incorporated | Acoustic counter-sniper system |
US20070143734A1 (en) * | 2005-12-16 | 2007-06-21 | Uzhakov Sergey V | Method and system for improving aerial image simulation speeds |
US8094907B1 (en) * | 2007-05-03 | 2012-01-10 | University Of Virginia | System, method and computer program product for fast conjugate phase reconstruction based on polynomial approximation |
US20170131397A1 (en) * | 2014-06-27 | 2017-05-11 | Robert Bosch Gmbh | Method for locating an object using an fmcw radar |
CN106249237A (zh) * | 2016-07-19 | 2016-12-21 | 西安电子科技大学 | 一种曲线轨迹下大斜视sar频域成像方法 |
CN109557542A (zh) * | 2018-11-28 | 2019-04-02 | 江苏天域航空科技集团股份有限公司 | 俯冲模式双基地前视sar成像方法 |
CN110208799A (zh) * | 2019-06-19 | 2019-09-06 | 华侨大学 | 基于Legendre正交分解的曲线运动轨迹SAR波数域成像方法 |
Non-Patent Citations (1)
Title |
---|
陈麒 等: "基于切比雪夫多项式的双基地前视SAR成像算法", 《控制与决策》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108710111B (zh) | 一种机载双基前视sar方位相位的二维空变校正方法 | |
CN109738894B (zh) | 一种用于大视场合成孔径雷达的高斜视多角度成像方法 | |
CN102879784B (zh) | 四种模式的合成孔径雷达统一成像方法 | |
CN109870686B (zh) | 一种基于改进斜距模型的曲线轨迹sar斜视成像方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN110673144B (zh) | 一种基于时间变标的子孔径大斜视sar成像处理方法 | |
CN110208799B (zh) | 基于Legendre正交分解的曲线运动轨迹SAR波数域成像方法 | |
CN109188432B (zh) | 一种平行双基聚束sar快速bp成像方法 | |
CN106597442B (zh) | 一种方位多通道脉内聚束sar成像方法 | |
CN111649803B (zh) | 基于垂直线性阵列的三维雷达物位计及其设计方法 | |
CN111856461A (zh) | 基于改进pfa的聚束sar成像方法及其dsp实现 | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN111208514B (zh) | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 | |
CN109856636B (zh) | 曲线合成孔径雷达自适应三维成像方法 | |
CN109597076B (zh) | 用于地基合成孔径雷达的数据处理方法及装置 | |
CN102608576A (zh) | 一种用于大前斜合成孔径雷达回波图像的几何校正方法 | |
CN107942295A (zh) | 一种前视阵列sar***的稀疏天线 | |
CN109765556A (zh) | 一种基于级数反演的双基sar快速几何校正方法及装置 | |
CN108008387B (zh) | 一种机载阵列天线下视三维成像方法 | |
CN102788978A (zh) | 一种斜视星机双基地合成孔径雷达成像方法 | |
CN111157991A (zh) | 一种基于Chebyshev多项式的曲线运动轨迹SAR的Omega-k成像方法 | |
CN112649808A (zh) | 基于移变构型的双基前视sar波数域成像方法 | |
CN112285707A (zh) | 基于gps导航信号的无源多基地高分辨成像方法 | |
CN110579762A (zh) | 一种太赫兹圆迹sar快速后向投影成像方法 | |
CN116136595A (zh) | 基于两级尺度微调的协同探测双基前视sar成像处理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |