CN1808172A - 机载干涉合成孔径雷达原始回波生成方法 - Google Patents
机载干涉合成孔径雷达原始回波生成方法 Download PDFInfo
- Publication number
- CN1808172A CN1808172A CN 200510011214 CN200510011214A CN1808172A CN 1808172 A CN1808172 A CN 1808172A CN 200510011214 CN200510011214 CN 200510011214 CN 200510011214 A CN200510011214 A CN 200510011214A CN 1808172 A CN1808172 A CN 1808172A
- Authority
- CN
- China
- Prior art keywords
- fft
- original echo
- exp
- zoom factor
- multiply
- 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
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及信息获取与处理技术领域,特别是一种机载干涉合成孔径雷达(InSAR)原始回波生成方法。其方法包括:步骤S1:对等效后向散射系数γ1,2(x,)作方位向傅里叶变换(FFT),然后乘以包含散焦项和缩放因子在内的相位函数P1=;步骤S2:对步骤S1得到的结果作距离向FFT,然后乘以包含缩放因子在内的相位函数 ;步骤S3:对步骤S2得到的结果作距离向逆FFT,然后乘以包含缩放因子在内的相位函数;;步骤S4:对步骤S3得到的结果作距离向FFT,然后乘以包含空不变冲激响应函数、 距离平移因子和缩放因子在内的函数P4=G(ξ,η,γ0)exp{-jηγs1,2+j[Ω1,2(ξ)-1]Bη2};步骤S5:对步骤S4得到的结果作二维逆FFT,可以得到所需的原始回波数据h1,2(x′,′)。
Description
技术领域
本发明涉及信息获取与处理技术领域,特别是一种机载干涉合成孔径雷达(InSAR)原始回波生成方法。
背景技术
在研制InSAR(干涉合成孔径雷达)系统的过程中,为了在获得雷达数据之前分析并设计相应的InSAR处理方法,需要利用计算机仿真的方法生成处理所需的各种InSAR数据。国外许多研究室为此设计了一系列的InSAR仿真系统,其中根据文献的分类方法,这些InSAR仿真方法大致可以归结为三大类,即相干系统仿真法、非相关图像仿真法以及基于SAR图像的仿真法。
为了研究能与机载双天线InSAR的特点相结合的新型InSAR成像处理方法,有必要在获得真实的雷达数据之前利用计算机仿真的方法生成InSAR原始回波数据,因此需要采用上述的第一类仿真方法。在此类仿真方法中,就InSAR原始回波数据生成方面而言,Alberti等人是在时域内利用卷积的方法实现的,而Franceschetti等人则是在频域内利用系统冲激响应函数的方法完成的。前者由于采用时域卷积的方法,在每个方位位置上都要计算雷达照射范围内所有目标回波的相干累加,需要消耗惊人的计算时间,因此只适用于简单的稀疏点阵目标仿真,对于复杂的分布目标仿真则显得难以适从。后者由于采用频域二维FFT(傅里叶变换)的方法,无需计算每个单点目标的回波信号,因此比较适用于复杂分布目标的仿真要求。不过,由于其在生成原始回波数据时并未将雷达偏离参考轨迹的情况以及两幅单视图像失配的情况考虑在内,因此无法满足新型InSAR成像处理方法研究的要求。
一般的正侧视机载双天线InSAR系统的成像几何关系如图1所示。天线A1和A2之间的基线垂直于飞机航线,其长为B,与水平方向的夹角为α。对于场景内的目标P(x,y,z),设两天线与目标间的最短距离分别为r和r+Δr,则根据余弦定理,有
(r+Δr)2=r2+B2-2Brsin(θ-α) (1)
其中
为天线视角,H为天线高度,z为目标高度。对于机载InSAR系统,天线基线B一般都远小于作用距离r,即B<<r,因此由1式可求得
Δr≈-Bsin(θ-α) (3)
假设场景地形的平均坡度为β,根据(2)式和(3)式,可以求得不同距离处两天线间的距离差Δr与场景中心处两天线间的距离差Δr0之间的关系为
Δr=Δr0+ks(r-r0) (4)
其中
对于一般的双天线InSAR系统,两副天线A1和A2既可以采用一副天线发送、两副天线接收的工作模式(即标准模式),也可以采用两副天线轮流发送并接收的工作模式(即乒乓模式)。对于机载双天线InSAR系统,设天线的发射信号为调频率为k、脉宽为τ的线性调频信号,则经过混频解调后的两副天线的接收信号可表示为
h1,2(x′,r′)=∫∫dxdrγ(x,r)g1,2(x′-x,r′-r;r) (6)
式中下标1和2分别表示由天线A1和天线A2所接收的信号,其中
为新场景的背散射系数,其中γ0(x,r)为场景的实际背散射系数,λ为发射信号载波波长,而g1,2(x′-x,r′-r;r)则为系统的冲激响应函数。
对(6)式中的天线接收信号求二维FT,有
H1,2(ξ,η)=G(ξ,η;r0)∫∫dxdrγ(x,r)exp(-jξx-jηr)1,2(ξ,η;r) (8)
其中G(ξ,η;r0)和1,2(ξ,η;r)分别为系统冲激响应函数的空不变部分和空变部分。利用驻定相位原理,并略去无关紧要的常数系数,可以求得
其中
这里的D为真实天线孔径长度,Δf为调频信号带宽,f0为载波中心频率。而
(12)
将(12)式分别代入(8)式中,并令
进行变量代换
并化简后,可以得到天线A1和A2所接收信号的频域表达式为
其中
rs1,2=r0+mΔr0 (15)
(14)式中的积分项可以表示为以Ω1,2(ξ)η为输出变量的
的距离向FT。利用缩放原理,可以通过两次卷积和两次二次相位平移将以η为输出变量的
的距离向FT变换为其以Ω1,2(ξ)η为输出变量的距离向FT。
发明内容
本发明的目的在于提供一种机载干涉合成孔径雷达原始回波生成方法。
本发明提出一种新的机载双天线InSAR原始回波数据生成的方法。该方法根据机载双天线InSAR回波数据对的特点,通过在距离向上引入平移和伸缩因子,充分地考虑了雷达的运动参数和位置关系的影响。此外,由于采用频域二维FFT的方法生成机载InSAR原始回波数据对的计算机仿真方法结构,因此可以大大缩短程序运行的时间。
本发明所提出的InSAR原始回波数据生成方法的基本思路可以概括为:先对场景后向散射系数γ1,2(x,r)进行缩放变换,求得其缩放后的散射系数频谱Γ1,2(ξ,Ω1,2(ξ)η);然后乘以空不变的系统冲激响应函数G(ξ,η;r0),得到原始回波数据的频谱H1,2(ξ,η);最后通过二维逆FT求得原始回波数据h1,2(x′,r′)。根据这一思路,可以得到如图2所示的InSAR原始回波数据生成的实现框图,其中图中的常数B理论上应该不受任何限制。
从图中可以看出,InSAR原始回波数据生成过程可以由两次方位FT和四次距离FT以及四次复数相乘来完成。其中第一复数相乘项中的前部分对应于原始数据中与聚焦深度有关的项,而后部分则对应于与距离缩放因子有关的项。最后一个复数相乘项除了与距离缩放因子有关的项外,还包括了空不变的系统冲击函数以及距离平移因子。整个方法充分考虑了InSAR系统的聚焦深度、耦合因子等各项成像因素,具有结构清晰简洁、实现简单有效等特点。
附图说明
图1是一般的正侧视机载InSAR系统的成像几何关系图。
图2是机载干涉合成孔径雷达原始回波生成方法流程图。
具体实施方式
图2机载干涉合成孔径雷达原始回波生成方法,步骤如下:
步骤S2:对步骤S1得到的结果作距离向FFT,然后乘以包含缩放因子在内的相位函数
步骤S3:对步骤S2得到的结果作距离向逆FFT,然后乘以包含缩放因子在内的相位函数
步骤S4:对步骤S3得到的结果作距离向FFT,然后乘以包含空不变冲激响应函数、距离平移因子和缩放因子在内的函数P4=G(ξ,η;r0)exp{-jηrs1,2+j[Ω1,2(ξ)-1]Bη2};
步骤S5:对步骤S4得到的结果作二维逆FFT,可以得到所需的原始回波数据h1,2(x′,r′)。
Claims (3)
1.一种机载干涉合成孔径雷达原始回波生成方法,其特征在于,通过在距离向上引入平移和伸缩因子模拟雷达运动参数和位置关系的影响,并利用频域二维FFT方法生成机载InSAR原始回波数据对的计算机仿真方法结构。
2.一种机载干涉合成孔径雷达原始回波生成方法,其步骤是:先对场景后向散射系数γ1,2(x,r)进行缩放变换,求得其缩放后的散射系数频谱Г1,2(ξ,Ω1,2(ξ)η);然后乘以空不变的系统冲激响应函数G(ξ,η;r0),得到原始回波数据的频谱H1,2(ξ,η);最后通过二维逆FT求得原始回波数据h1,2(x′,r′)。
3.根据权利要求1或2的机载干涉合成孔径雷达原始回波生成方法,
其特征在于,包括如下步骤:
步骤S2:对步骤S1得到的结果作距离向FFT,然后乘以包含缩放因子在内的相位函数
步骤S3:对步骤S2得到的结果作距离向逆FFT,然后乘以包含缩放因子在内的相位函数
步骤S4:对步骤S3得到的结果作距离向FFT,然后乘以包含空不变冲激响应函数、距离平移因子和缩放因子在内的函数P4=G(ξ,η;r0)exp{-jηrs1,2+j[Ω1,2(ξ)-1]Bη2};
步骤S5:对步骤S4得到的结果作二维逆FFT,可以得到所需的原始回波数据h1,2(x′,r′)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200510011214XA CN100547432C (zh) | 2005-01-20 | 2005-01-20 | 机载干涉合成孔径雷达原始回波生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200510011214XA CN100547432C (zh) | 2005-01-20 | 2005-01-20 | 机载干涉合成孔径雷达原始回波生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1808172A true CN1808172A (zh) | 2006-07-26 |
CN100547432C CN100547432C (zh) | 2009-10-07 |
Family
ID=36840173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB200510011214XA Expired - Fee Related CN100547432C (zh) | 2005-01-20 | 2005-01-20 | 机载干涉合成孔径雷达原始回波生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100547432C (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101846740A (zh) * | 2010-05-17 | 2010-09-29 | 北京航空航天大学 | 一种星载sar指定纬度区域回波仿真方法 |
CN101876704A (zh) * | 2010-06-03 | 2010-11-03 | 中国人民解放军国防科学技术大学 | 干涉合成孔径雷达三维陆地场景回波仿真方法 |
CN101887121A (zh) * | 2010-06-22 | 2010-11-17 | 北京航空航天大学 | 基于半牛顿迭代法的星载干涉合成孔径雷达基线估计方法 |
CN101482616B (zh) * | 2008-08-13 | 2011-06-01 | 中国科学院电子学研究所 | 一种地形测量方法 |
CN101587500B (zh) * | 2008-05-23 | 2011-11-09 | 中国科学院电子学研究所 | 双站合成孔径雷达海面成像的计算机仿真方法 |
CN102253384A (zh) * | 2010-05-19 | 2011-11-23 | 中国科学院电子学研究所 | 一种多通道或多时相雷达成像方法 |
CN102466799A (zh) * | 2010-10-29 | 2012-05-23 | 中国科学院电子学研究所 | 基于pos运动数据的干涉sar回波数据仿真方法 |
CN102478653A (zh) * | 2010-11-22 | 2012-05-30 | 电子科技大学 | 一种基于距离分割的sar回波时频混合模拟方法 |
CN102544755A (zh) * | 2011-12-31 | 2012-07-04 | 哈尔滨工业大学 | 一种基于强散射点的均匀线阵校准方法 |
CN102830390A (zh) * | 2012-08-24 | 2012-12-19 | 电子科技大学 | 一种双基地sar回波模拟方法 |
-
2005
- 2005-01-20 CN CNB200510011214XA patent/CN100547432C/zh not_active Expired - Fee Related
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101587500B (zh) * | 2008-05-23 | 2011-11-09 | 中国科学院电子学研究所 | 双站合成孔径雷达海面成像的计算机仿真方法 |
CN101482616B (zh) * | 2008-08-13 | 2011-06-01 | 中国科学院电子学研究所 | 一种地形测量方法 |
CN101846740A (zh) * | 2010-05-17 | 2010-09-29 | 北京航空航天大学 | 一种星载sar指定纬度区域回波仿真方法 |
CN102253384B (zh) * | 2010-05-19 | 2013-06-26 | 中国科学院电子学研究所 | 一种多通道或多时相雷达成像方法 |
CN102253384A (zh) * | 2010-05-19 | 2011-11-23 | 中国科学院电子学研究所 | 一种多通道或多时相雷达成像方法 |
CN101876704B (zh) * | 2010-06-03 | 2012-01-11 | 中国人民解放军国防科学技术大学 | 干涉合成孔径雷达三维陆地场景回波仿真方法 |
CN101876704A (zh) * | 2010-06-03 | 2010-11-03 | 中国人民解放军国防科学技术大学 | 干涉合成孔径雷达三维陆地场景回波仿真方法 |
CN101887121B (zh) * | 2010-06-22 | 2012-06-27 | 北京航空航天大学 | 基于半牛顿迭代法的星载干涉合成孔径雷达基线估计方法 |
CN101887121A (zh) * | 2010-06-22 | 2010-11-17 | 北京航空航天大学 | 基于半牛顿迭代法的星载干涉合成孔径雷达基线估计方法 |
CN102466799B (zh) * | 2010-10-29 | 2013-06-26 | 中国科学院电子学研究所 | 基于pos运动数据的干涉sar回波数据仿真方法 |
CN102466799A (zh) * | 2010-10-29 | 2012-05-23 | 中国科学院电子学研究所 | 基于pos运动数据的干涉sar回波数据仿真方法 |
CN102478653A (zh) * | 2010-11-22 | 2012-05-30 | 电子科技大学 | 一种基于距离分割的sar回波时频混合模拟方法 |
CN102478653B (zh) * | 2010-11-22 | 2014-01-01 | 电子科技大学 | 一种基于距离分割的sar回波时频混合模拟方法 |
CN102544755A (zh) * | 2011-12-31 | 2012-07-04 | 哈尔滨工业大学 | 一种基于强散射点的均匀线阵校准方法 |
CN102544755B (zh) * | 2011-12-31 | 2013-12-11 | 哈尔滨工业大学 | 一种基于强散射点的均匀线阵校准方法 |
CN102830390A (zh) * | 2012-08-24 | 2012-12-19 | 电子科技大学 | 一种双基地sar回波模拟方法 |
CN102830390B (zh) * | 2012-08-24 | 2014-03-05 | 电子科技大学 | 一种双基地sar回波模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100547432C (zh) | 2009-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1808172A (zh) | 机载干涉合成孔径雷达原始回波生成方法 | |
CN111521989B (zh) | 用于雷达系统中的超分辨率的深度学习 | |
CN108107431B (zh) | 一种圆柱扫描sar三维成像快速实现方法 | |
CN101900812B (zh) | 一种圆迹合成孔径雷达的大场景极坐标格式三维成像方法 | |
US20210018593A1 (en) | Methods and systems for processing radar reflections | |
KR20210093967A (ko) | 희소 안테나 어레이로부터 형성되는 가상의 어퍼처를 사용한 잠재적으로 움직이는 참조 프레임으로부터의 오브젝트 감지 | |
EP0298112B1 (en) | Processing parameter generator for synthetic aperture radar | |
CN109738894B (zh) | 一种用于大视场合成孔径雷达的高斜视多角度成像方法 | |
WO2008073011A1 (en) | A sar radar system and a method relating thereto | |
CN110596646A (zh) | 一种基于mimo体制的提高雷达角度分辨率的布局及方法 | |
CN102812379B (zh) | 超分辨率成像雷达 | |
CN103293528B (zh) | 扫描雷达超分辨成像方法 | |
CN111145337A (zh) | 基于分辨率逼近的快速稀疏重构的线阵sar三维成像方法 | |
CN102176018A (zh) | 机械扫描雷达的多普勒波束锐化快速成像方法 | |
CN110208760A (zh) | 一种基于时域上采样的雷达回波仿真方法 | |
Qiu et al. | Bistatic SAR data processing algorithms | |
CN111722227B (zh) | 基于近似观测矩阵的聚束sar压缩感知成像方法 | |
WO2019118050A1 (en) | System and method for synthetic aperture radar image formation | |
CN101566689B (zh) | 一种消除三维合成孔径雷达图像栅瓣的方法 | |
CN110231616B (zh) | 基于北斗卫星辐射源的海面运动目标检测与定位方法 | |
Giroux et al. | T-fftradnet: Object detection with swin vision transformers from raw adc radar signals | |
Li et al. | Azimuth super-resolution for fmcw radar in autonomous driving | |
CN113608218A (zh) | 一种基于后向投影原理的频域干涉相位稀疏重构方法 | |
Chen et al. | Analysis and Compensation for Systematical Errors in Airborne Microwave Photonic SAR Imaging by 2-D Autofocus | |
CN117129994A (zh) | 基于相位补偿核gnss-sar改进后向投影成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20091007 Termination date: 20160120 |
|
EXPY | Termination of patent right or utility model |