CN114609633B - 一种圆周聚束模式干涉sar测高方法 - Google Patents
一种圆周聚束模式干涉sar测高方法 Download PDFInfo
- Publication number
- CN114609633B CN114609633B CN202210263135.1A CN202210263135A CN114609633B CN 114609633 B CN114609633 B CN 114609633B CN 202210263135 A CN202210263135 A CN 202210263135A CN 114609633 B CN114609633 B CN 114609633B
- Authority
- CN
- China
- Prior art keywords
- antenna
- target
- main antenna
- radar
- distance
- 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
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/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric 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/904—SAR modes
- G01S13/9047—Doppler beam sharpening mode
-
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明属于雷达微波干涉测量技术领域,具体涉及一种圆周聚束模式干涉SAR测高方法。本发明利用合成孔径雷达常用的线性调频信号的回波信号进行分析,建立圆周聚束模式下干涉测高模型,根据空间立体几何关系构建四面体模型,推导目标高程反演公式,公式除干涉相位外,其他参数均是由雷达参数决定的常数,不涉及复杂的计算操作,实现方式也较为简单。
Description
技术领域
本发明属于干涉合成孔径雷达微波测量领域,具体涉及一种圆周聚束模式干涉SAR测高方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)是一种主动式微波遥感测量技术,作为一种高分辨率的二维成像系统,它发射宽脉冲的线性调频信号,解决了距离向分辨率与作用距离之间的矛盾,利用真实孔径的小天线相对于目标运动,等效于等间隔的天线阵元在空间上合成一个长的真实孔径天线,从而获得方位向高分辨率,随着射频器件水平及数字技术的不断提升,SAR系统功能逐步完善。例如,干涉合成孔径雷达(Interferometric Synthetic Aperture Radar,InSAR),在机载单轨测量模式下,它通过两幅天线同时观测获取地面同一目标的复影像对,根据目标与两天线位置的几何关系和复图像上产生的相位差信息可以提取地面目标地形的三维信息,将地形测绘由二维平面拓展到三维空间;视频合成孔径雷达(Video Synthetic Aperture Radar,ViSAR),通过对场景的持续观测获取目标场景的连续的图像帧,使其能检测目标场景的动态变化。为了将视频SAR与干涉SAR技术相结合,必须考虑视频SAR成像算法是在圆迹聚束模式下的,这与传统的条带模式的干涉SAR有所区别,目前还没有成熟实用的针对圆周聚束模式干涉测高模型的建立。
发明内容
针对上述问题,本发明提出了一种圆周聚束模式干涉SAR测高方法,该方法的基本思路是假设雷达平台在圆周聚束模式下沿着航迹轨道匀速飞行,先根据发射的线性调频信号,计算其回波经过匹配滤波后的回波干涉相位,考虑载机按照圆迹轨迹飞行,主副天线与空间目标的位置不在一个平面,根据目标与雷达主副天线的空间几何关系构建四面体模型,得到目标高度与主副天线与目标之间距离差的关系,进而得到目标高度与相位差的关系式,从而得到高程反演公式。
本发明的技术方案为:一种圆周聚束模式干涉SAR测高方法,用于机载双天线视频SAR干涉测量系统,包括以下步骤:
S1、计算回波相位,通常情况下,雷达平台发射具有大时带宽积的线性调频信号,在接收到回波信号后,首先对回波数据做距离向的匹配滤波,得到处理后回波信号的相位;
S2、计算干涉相位和主副天线与目标之间距离差的关系,信号共轭相乘体现在相位计算上为相位相减得到相位差,主副天线与目标点的斜距不同,从而产生不同的相位,因此得到和主副天线与目标距离差有关的相位差;
S3、根据空间几何关系,构建四面体模型,原点向上平移以目标高度为单位的点O′,主天线所在的位置点A1,副天线所在的位置点A2,目标点P,这四个点在空间组成四面体O′PA1A2,根据勾股定理与余弦定理计算目标高度和主副天线与目标之间距离差的关系;
S4、根据目标高度和主副天线与目标之间距离差的关系,主副天线与目标之间距离差和相位差的关系,消去主副天线与目标之间距离差的影响,得到目标高度与相位差的关系式,从而得到目标高程反演公式;
本发明的有益效果为,本方法提出了聚束模式的视频SAR干涉测量模型及高程反演公式,并针对相关步骤给出了具体方法。
附图说明
图1圆周聚束模式干涉测高模型;
图2空间三棱锥示意图;
具体实施方式
下面结合附图对本发明进行详细的描述。
本发明的适用于机载双天线干涉SAR系统,主要包括:
S1、计算回波相位,通常情况下,雷达平台发射具有大时带宽积的线性调频信号,在接收到回波信号后,首先对回波数据做距离向的匹配滤波,处理后回波信号的相位可写为:
其中t为方位向时间,为慢时间,fτ表示快时间频率,Rref表示场景中心和雷达天线中心之间的参考距离,Ri用于表示点目标和雷达天线中心的斜距;
如图1所示,雷达机载平台在距离地面高度为H的位置按照圆迹飞行,主天线记为A1,副天线记为A2,地面场景中心为Q,主天线与场景中心之间的参考距离为Rref1,副天线与场景中心之间的参考距离为Rref2,对于地面高度为h的点目标P,距离主天线的距离为R1,记为R,距离副天线的距离为R2,记为R+△R,α为基线B与水平方向的夹角,主天线在参考平面的投影点为原点,记为O,x轴为零多普勒位置载机飞行方向,z轴方向朝上,y轴由右手坐标系确定,β为主天线和目标连线与主天线和原点连线之间的夹角,点O′为原点O向上平移h个单位的点,设fc为雷达载频,λ为雷达波长,主天线A1和副天线A2之间的连线称为基线,基线的长度记为B。
主天线A1得到匹配滤波处理后的回波相位为:
主天线A2得到匹配滤波处理后的回波相位为:
S2、计算干涉相位和主副天线与目标之间距离差的关系,信号共轭相乘体现在相位计算上为相位相减得到相位差,主副天线与目标点的斜距不同,从而产生不同的相位,因此得到和主副天线与目标距离差有关的相位差;
将主天线A1得到匹配滤波处理后的回波相位减去主天线A2得到匹配滤波处理后的回波相位得到相位差为:
S3、根据空间几何关系,构建四面体模型,原点向上平移以目标高度为单位的点O′,主天线所在的位置点A1,副天线所在的位置点A2,目标点P,这四个点在空间组成四面体O′PA1A2,根据勾股定理与余弦定理计算目标高度和主副天线与目标之间距离差的关系;
如图2所示,原点向上平移以目标高度为单位距离的点O′,主天线所在位置点A1,副天线所在位置点A2,目标点P,这四个点在空间组成四面体O′PA1A2,具体求目标高度和主副天线与目标之间距离差的关系分为以下步骤:
第一步:求O′A2的长度:
在△O′A1A2中,根据余弦定理可得O′A2的长度:
第二步:求O′P的长度:
在△O′A1P中,因为O′P⊥O′A1,根据余弦定义得
O′P=Rsinβ
第三步:先求A2P的长度,再得到夹角β的计算式,
在△O′PA2中,O′P⊥O′A2,根据勾股定理得:
(R+△R)2=|A2P|2=|O′P|2+|O′A2|2=B2+(H-h)2-2B(H-h)cos(α+90°)+R2sin2β
变形可得:
sin2β=f(△R)=[(R+△R)2-B2+(H-h)2+2B(H-h)sinα]/R2
第四步:根据几何关系得高程反演公式为:
在直角△A1PO′中,根据余弦定理得目标高度与夹角β的关系:
h=H-Rcosβ
联立上面两式消去夹角β可得:
S4、根据目标高度和主副天线与目标之间距离差的关系,主副天线与目标之间距离差与相位差的关系,消去距离差的影响,得到目标高度与相位差的关系式;
距离差可以由相位差表示为:
将距离差表达式代入目标高度与距离差的关系式得到高程反演公式为:
高程反演公式看似复杂,其实除了相位差△φ之外,均为定值常数,在雷达干涉过程中,只要参数设定,便可由相位值得到目标的高程信息。
Claims (1)
1.一种圆周聚束模式干涉SAR测高方法,用于机载双天线视频SAR干涉测量系统,其特征在于,根据圆周聚束模式下的成像模式,建立圆周聚束模式干涉测高模型,定义雷达机载平台在距离地面高度为H的位置按照圆迹飞行,主天线记为A1,副天线记为A2,以主天线在参考平面的投影点为原点,记为O,x轴为零多普勒位置载机飞行方向,z轴方向朝上,y轴由右手坐标系确定,地面场景中心为Q,对于地面高度为h的点目标P,距离主天线的距离为R,距离副天线的距离为R+△R;所述测高方法包括:
S1、雷达机载平台接收到回波信号后,对回波数据做距离向的匹配滤波,得到相位:
其中t为方位向时间,为慢时间,fτ表示快时间频率,Rref表示地面场景中心和雷达天线中心之间的参考距离,Ri表示点目标和雷达天线中心的斜距;
主天线A1得到匹配滤波处理后的回波相位为:
主天线A2得到匹配滤波处理后的回波相位为:
其中Rref1为主天线与地面场景中心之间的参考距离,Rref2为副天线与地面场景中心之间的参考距离;
S2、将主天线A1得到匹配滤波处理后的回波相位减去主天线A2得到匹配滤波处理后的回波相位得到相位差为:
S3、根据空间几何关系,构建四面体模型,原点向上平移以目标高度为单位的点O′,主天线所在的位置点A1,副天线所在的位置点A2,目标点P,这四个点在空间组成四面体O′PA1A2,根据勾股定理与余弦定理计算目标高度和主副天线与目标之间距离差的关系;
其中,B为基线的长度,α为基线B与水平方向的夹角;
S4、根据目标高度和主副天线与目标之间距离差的关系,主副天线与目标之间距离差和相位差的关系,消去主副天线与目标之间距离差的影响,得到目标高度与相位差的关系式,从而得到目标高程反演公式,得到目标高程信息:
其中,fc为雷达载频,λ为雷达波长。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210263135.1A CN114609633B (zh) | 2022-03-17 | 2022-03-17 | 一种圆周聚束模式干涉sar测高方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210263135.1A CN114609633B (zh) | 2022-03-17 | 2022-03-17 | 一种圆周聚束模式干涉sar测高方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114609633A CN114609633A (zh) | 2022-06-10 |
CN114609633B true CN114609633B (zh) | 2023-09-01 |
Family
ID=81865084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210263135.1A Active CN114609633B (zh) | 2022-03-17 | 2022-03-17 | 一种圆周聚束模式干涉sar测高方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114609633B (zh) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
CN101551455A (zh) * | 2009-05-13 | 2009-10-07 | 西安电子科技大学 | 干涉合成孔径雷达三维地形成像系统及其高程测绘方法 |
CN103543453A (zh) * | 2013-10-28 | 2014-01-29 | 北京理工大学 | 一种地球同步轨道合成孔径雷达干涉的高程反演方法 |
CN104007439A (zh) * | 2014-05-30 | 2014-08-27 | 中国科学院电子学研究所 | 一种干涉圆迹sar高程估计处理方法 |
CN104698457A (zh) * | 2014-09-02 | 2015-06-10 | 电子科技大学 | 一种迭代曲面预测InSAR成像及高度估计方法 |
CN105929399A (zh) * | 2016-04-25 | 2016-09-07 | 电子科技大学 | 一种干涉sar数据成像及高程估计方法 |
CN106199601A (zh) * | 2016-07-01 | 2016-12-07 | 西安电子科技大学 | 基于粗数字高程图的InSAR绝对相位模糊估计方法 |
CN107607943A (zh) * | 2017-09-07 | 2018-01-19 | 西安电子科技大学 | 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法 |
CN108594222A (zh) * | 2018-03-21 | 2018-09-28 | 中国科学院电子学研究所 | 一种双频干涉合成孔径雷达的高程重建方法和装置 |
CA2968794A1 (en) * | 2017-05-31 | 2018-11-30 | Institut National D'optique | Synthetic aperture imaging assisted by three-dimensional scanning imaging for height reconstruction |
KR20190043478A (ko) * | 2017-10-18 | 2019-04-26 | 서울대학교산학협력단 | 고해상도 수치 표고 모델 생성 시스템 및 방법 |
CN112882030A (zh) * | 2021-01-12 | 2021-06-01 | 中国科学院空天信息创新研究院 | InSAR成像干涉一体化处理方法 |
CN113093184A (zh) * | 2021-03-31 | 2021-07-09 | 电子科技大学 | 一种基于视频合成孔径雷达的干涉测量方法 |
CN113640797A (zh) * | 2021-08-09 | 2021-11-12 | 北京航空航天大学 | 一种用于参考条带模式InSAR的前斜视测高方法 |
CN113835089A (zh) * | 2021-08-27 | 2021-12-24 | 中国空间技术研究院 | 一种基于差频InSAR的缓变高程反演方法 |
-
2022
- 2022-03-17 CN CN202210263135.1A patent/CN114609633B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
CN101551455A (zh) * | 2009-05-13 | 2009-10-07 | 西安电子科技大学 | 干涉合成孔径雷达三维地形成像系统及其高程测绘方法 |
CN103543453A (zh) * | 2013-10-28 | 2014-01-29 | 北京理工大学 | 一种地球同步轨道合成孔径雷达干涉的高程反演方法 |
CN104007439A (zh) * | 2014-05-30 | 2014-08-27 | 中国科学院电子学研究所 | 一种干涉圆迹sar高程估计处理方法 |
CN104698457A (zh) * | 2014-09-02 | 2015-06-10 | 电子科技大学 | 一种迭代曲面预测InSAR成像及高度估计方法 |
CN105929399A (zh) * | 2016-04-25 | 2016-09-07 | 电子科技大学 | 一种干涉sar数据成像及高程估计方法 |
CN106199601A (zh) * | 2016-07-01 | 2016-12-07 | 西安电子科技大学 | 基于粗数字高程图的InSAR绝对相位模糊估计方法 |
CA2968794A1 (en) * | 2017-05-31 | 2018-11-30 | Institut National D'optique | Synthetic aperture imaging assisted by three-dimensional scanning imaging for height reconstruction |
CN107607943A (zh) * | 2017-09-07 | 2018-01-19 | 西安电子科技大学 | 基于干涉相位辅助的延迟多普勒雷达高度表的测高方法 |
KR20190043478A (ko) * | 2017-10-18 | 2019-04-26 | 서울대학교산학협력단 | 고해상도 수치 표고 모델 생성 시스템 및 방법 |
CN108594222A (zh) * | 2018-03-21 | 2018-09-28 | 中国科学院电子学研究所 | 一种双频干涉合成孔径雷达的高程重建方法和装置 |
CN112882030A (zh) * | 2021-01-12 | 2021-06-01 | 中国科学院空天信息创新研究院 | InSAR成像干涉一体化处理方法 |
CN113093184A (zh) * | 2021-03-31 | 2021-07-09 | 电子科技大学 | 一种基于视频合成孔径雷达的干涉测量方法 |
CN113640797A (zh) * | 2021-08-09 | 2021-11-12 | 北京航空航天大学 | 一种用于参考条带模式InSAR的前斜视测高方法 |
CN113835089A (zh) * | 2021-08-27 | 2021-12-24 | 中国空间技术研究院 | 一种基于差频InSAR的缓变高程反演方法 |
Non-Patent Citations (1)
Title |
---|
"PollnSAR/TomoSAR 植被和建筑物高度反演研究";解金吾;《中国博士学位论文全文数据库 信息科技辑》(第7(2020)期);第I136-69页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114609633A (zh) | 2022-06-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5659318A (en) | Interferometric SAR processor for elevation | |
CN112305510B (zh) | 一种基于dem匹配的合成孔径雷达影像几何定标方法 | |
CN112034460B (zh) | 基于天线相位方向图补偿的圆弧孔径雷达成像方法和雷达 | |
CN109633642B (zh) | 一种太赫兹高速目标雷达成像方法 | |
CN105372657A (zh) | 基于回波数据的视频合成孔径雷达运动补偿成像方法 | |
CN110068817B (zh) | 一种基于激光测距和InSAR的地形测图方法、仪器和系统 | |
CN113093184B (zh) | 一种基于视频合成孔径雷达的干涉测量方法 | |
JP6381825B2 (ja) | 移動目標検出装置 | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN104122552B (zh) | 一种滑动式双站圆周合成孔径雷达成像方法 | |
KR100441590B1 (ko) | 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용하여지형고도를 측정하기 위한 디지털 고도모형 생성방법 | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 | |
CN107728144B (zh) | 一种基于前视双基模式的干涉sar成像方法 | |
CN111707996B (zh) | 基于改进grft-stap的geo星机sar动目标检测方法 | |
CN114609633B (zh) | 一种圆周聚束模式干涉sar测高方法 | |
CN115712095A (zh) | 基于单个角反射的sar卫星三维定位误差改正方法及系统 | |
CN116148891A (zh) | 基于顺轨干涉的星地双基海面二维流场测量方法 | |
Yu et al. | Preliminary Analgsis of Geometric Positioning Accuracy Based on Gaofen-3 Data | |
Gong et al. | High resolution 3D InISAR imaging of space targets based on PFA algorithm with single baseline | |
CN112859018A (zh) | 一种基于图像几何校正的视频sar成像方法 | |
CN115291213B (zh) | 一种基于双波束sar的动目标方位速度估计方法 | |
Lim et al. | Phase statistical model and correction in imagery of ground based synthetic aperture radar (GBSAR) for land deformation monitoring | |
CN114839632B (zh) | 一种联合非摄影测量观测条件约束的sar无控几何定标方法及其系统 | |
De Macedo et al. | Recent LC-and X-band MetaSensing airborne SAR campaigns for emerging applications | |
Wang et al. | Dual-Mode interferometric experiment using Ka DBF-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 |