CN105866777B - 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 - Google Patents
多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 Download PDFInfo
- Publication number
- CN105866777B CN105866777B CN201610187369.7A CN201610187369A CN105866777B CN 105866777 B CN105866777 B CN 105866777B CN 201610187369 A CN201610187369 A CN 201610187369A CN 105866777 B CN105866777 B CN 105866777B
- Authority
- CN
- China
- Prior art keywords
- angle
- dimensional
- pixel
- satellite
- period
- 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/9058—Bistatic or multistatic SAR
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)
- Image Processing (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种多角度多时段导航卫星双基地PS‑InSAR三维形变反演方法,该方法通过在待监测的场景附近放置接收机,利用全向天线接收卫星直达波信号,利用喇叭天线接收场景回波信号,利用直达波信号对回波信号进行同步、成像处理获得多角度多时段的SAR图像序列,通过PS方法实现PS点的识别与形变量估计,再利用空间、时间插值实现不同角度下的形变量空间、时间的连续化,最后利用加权最小二乘实现三维形变量的反演;本发明成本低,配置灵活,可广泛应用在形变监测领域。
Description
技术领域
本发明属于雷达信号处理的技术领域,具体涉及一种多角度多时段导航卫星双基地PS-InSAR三维形变反演方法。
背景技术
基于导航卫星的双基地重轨差分干涉SAR利用导航卫星作为发射平台,地面配置接收机,对需要探测形变的区域进行连续成像、差分干涉和形变反演处理可获得监测区域的一维形变。
然而,单颗导航卫星的重轨时间较长,同时仅能实现一维形变的反演,精度受限,实用性较差。实际上,导航卫星系统由多颗卫星组成,对于同一地区始终有多颗卫星的信号连续覆盖,而且随着时间的不同,卫星相对于一个特定场景的几何关系也会发生变化,即可提供不同视角的信息。因此,该系统在提高形变监测精度方面具有很大潜力。
因此,开发一种基于多角度多时段导航卫星照射的双基地PS-InSAR的三维连续形变监测方法,对于高精度形变测量领域具有重要意义。
发明内容
有鉴于此,本发明提供了一种多角度多时段导航卫星双基地PS-InSAR三维形变反演方法,能够利用多卫星、多时段的回波数据实现照射场景的三维形变测量,其精度高、成本低且时间连续。
实现本发明的技术方案如下:
步骤一:利用不同角度不同时段下的直达波信号对相应角度相应时段的回波信号进行同步和成像处理,得到同一场景的不同角度的双基地SAR图像序列;
步骤二:对所述双基地SAR图像序列,利用PS方法实现各角度下的PS点识别与PS点的一维形变量的反演;
步骤三:根据所述PS点的一维形变量的反演结果,利用Kriging插值方法对所述双基地SAR图像序列的每个像素的一维形变量进行反演;
步骤四:设置一个统一的时间轴,利用时间序列分析方法对所述双基地SAR图像序列的每个像素的一维形变量的反演结果进行时间插值,得到同一场景的不同角度下的每个像素在统一的时间轴下的一维形变量历史;
步骤五:利用加权最小二乘估计方法,对所述每个像素在统一的时间轴下的一维型变量历史进行处理,实现每个像素的三维形变量历史的反演。
有益效果:
本发明对比已有技术,充分利用了导航卫星数目多、导航信号覆盖范围广、照射时间长的特点,与监测区域附近的接收机配置构成双基地SAR系统,可对监测区域实现时间、空间连续的三维形变测量,能够获取更加全面的形变信息。同时作为被动测量方式,本发明成本低,配置灵活,可广泛应用在形变监测领域。
附图说明
图1为本发明系统构型示意图。
图2为本发明算法总流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明方法所利用的系统构型如图1所示,不同卫星从不同角度发射导航信号到场景,场景反射回波到回波天线,同时,一个全向直达波天线接收卫星的直达波信号。对于每颗卫星,采集其不同重轨时段的信号。以接收机直达波天线作为原点O,正东为X轴,正北为Y轴,竖直向上为Z轴建立空间直角坐标系O-XYZ。如图2所示为本发明方法的流程,该方法具体步骤如下:
步骤一:利用不同角度不同时段下的直达波信号对相应角度相应时段的回波信号进行同步和成像处理,得到同一场景的不同角度下的双基地SAR图像序列。设第i颗卫星第j组时段(一组时段对应同一个角度)的Mij+1幅重轨图像序列为Qi,j:
其中,S表示卫星总数,Ki表示第i颗卫星的时段总数。qi,j,0表示第i颗卫星第j组时段的第一幅图像,qi,j,1表示第i颗卫星第j组时段的第二幅图像,依此类推。
第i颗卫星第j组时段的卫星重轨时间序列Ti,j为:
其中,ti,j,0表示第i颗卫星第j组时段的第一个重轨时刻,ti,j,1表示第i颗卫星第j组时段的第二个重轨时刻,其他依此类推。
由于不同时段,卫星相对于观测场景的观测角度不同,因此,所有卫星所有时段的观测角度总数L为:
步骤二:对第i颗卫星第j组时段的Mij+1幅重轨图像进行PS点识别与PS点的一维形变量反演,得到各个PS点沿着等效视线方向的一维形变反演结果ΔLi,j(Pi,j;ti,j,0),设为(设qi,j,0为主图像):
其中,Pi,j为第i颗卫星的第j组图像中的PS点集合。Δli,j(Pi,j;ti,j,0,ti,j,1)表示第i颗卫星的第j组图像中PS点沿着等效视线方向的ti,j,1时刻相对于ti,j,0时刻的形变反演结果。Δli,j(Pi,j;ti,j,0,ti,j,2)表示第i颗卫星的第j组图像中PS点沿着等效视线方向的ti,j,2时刻相对于ti,j,0时刻的形变反演结果。其他依此类推。
步骤三:由于不同卫星不同角度下探测到的PS点不同,无法直接计算某个像素的三维形变量,因此,根据不同角度下SAR图像中各个PS点的一维形变反演结果,利用Kriging插值方法,对每个角度下整个场景的双基地SAR图像序列的每个像素进行空间插值,得到空间插值后的每个角度下场景中每个像素的形变量矩阵:
其中,表示场景中的像素点集合。表示空间插值后的Δli,j(Pi,j;ti,j,0,ti,j,1),表示空间插值后的Δli,j(Pi,j;ti,j,0,ti,j,2),其他依此类推。
步骤四:由于卫星的重轨时刻未必完全相同,为实现场景中每个像素三维形变量历史的提取,需要将所获取的不同角度下的形变量历史的时间轴进行统一。设统一时间轴为:
t=[t0,t1,…,tN-1] (6)
其中,t0为时间轴的第1个时刻,tk=t0+k·T(k=1,2,…,N-1)为该时间轴的第k+1个时刻,其中,T表示统一时间轴的时间间隔,N表示统一时间轴的时刻总数。
那么,对进行基于时间序列分析方法的时间插值得到空间、时间插值后的形变量矩阵:
其中,表示时间插值后的形变量在t0时刻的形变量,其他依此类推。
步骤五:
设在第tm(m=0,1,…N-1)时刻,场景中的点的三维形变量为:
其中,和分别为点在第tm时刻沿着前述坐标系O-XYZ上的X、Y和Z轴方向的形变量。
则有:
其中,为噪声,标准差为为点在第i颗卫星第j组时段的形变量进行时间插值和空间插值后在时刻tm的形变量,为等效视线方向的向量,可写为:
其中:
Ψ(Sij)=[xij,yij,zij]T (11)
表示卫星i在第j组时段的孔径中心时刻的三维位置,其中xij、yij和zij分别表示卫星i在第j组时段的孔径中心时刻在前述坐标系O-XYZ上的X、Y和Z轴方向的坐标;
Ψ(R)=[xR,yR,zR]T (12)
为回波天线R的三维位置。其中xR、yR和zR分别表示回波天线R在前述坐标系O-XYZ上的X、Y和Z轴方向的坐标;
为点的三维位置。其中和分别表示点在前述坐标系O-XYZ上的X、Y和Z轴方向的坐标;
那么,点处三维形变量的加权最小二乘估计为:
其中:
Φ为权值矩阵,为:
其中,diag表示对角矩阵,即:Φ的第一行第一列为第二行第二列为其他依此类推。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.多角度多时段导航卫星双基地PS-InSAR三维形变反演方法,其特征在于,包括以下步骤:
步骤一:利用不同角度不同时段下的直达波信号对相应角度相应时段的回波信号进行同步和成像处理,得到同一场景的不同角度的双基地SAR图像序列;
步骤二:对所述双基地SAR图像序列,利用PS方法实现各角度下的PS点识别与PS点的一维形变量的反演;
步骤三:根据所述PS点的一维变形量的反演结果,利用Kriging插值方法对所述双基地SAR图像序列的每个像素的一维形变量进行反演;
步骤四:设置一个统一的时间轴,利用时间序列分析方法对所述双基地SAR图像序列的每个像素的一维形变量的反演结果进行时间插值,得到同一场景的不同角度下的每个像素在统一的时间轴下的一维形变量历史;
步骤五:利用加权最小二乘估计方法,对所述每个像素在统一的时间轴下的一维型变量历史进行处理,实现每个像素的三维形变量历史的反演。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610187369.7A CN105866777B (zh) | 2016-03-29 | 2016-03-29 | 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610187369.7A CN105866777B (zh) | 2016-03-29 | 2016-03-29 | 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105866777A CN105866777A (zh) | 2016-08-17 |
CN105866777B true CN105866777B (zh) | 2018-10-16 |
Family
ID=56625937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610187369.7A Active CN105866777B (zh) | 2016-03-29 | 2016-03-29 | 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105866777B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108507454B (zh) * | 2018-03-09 | 2019-12-03 | 北京理工大学 | 一种基于导航卫星Bi-InSAR形变反演图像提取方法 |
CN109613531B (zh) * | 2019-01-07 | 2020-10-02 | 北方工业大学 | 一种微变感知预警雷达的多阈值优化形变反演方法及系统 |
CN110888132B (zh) * | 2019-11-22 | 2022-05-20 | 深圳市城市公共安全技术研究院有限公司 | 基于InSAR监测的桥梁形变分析方法及系统 |
CN112284332B (zh) * | 2020-08-31 | 2021-10-08 | 北京四象爱数科技有限公司 | 基于高分辨率insar的高层建筑沉降监测结果三维定位方法 |
CN112986990B (zh) * | 2021-02-04 | 2023-02-17 | 中国地质大学(北京) | 一种大气相位改正方法及系统 |
CN113777606B (zh) * | 2021-08-12 | 2023-12-26 | 北京理工大学 | 分布式geo sar三维形变反演多角度选取方法及装置 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2566132T3 (es) * | 2008-07-04 | 2016-04-11 | Telespazio S.P.A. | Identificación y análisis de dispersores persistentes en series de imágenes de SAR |
CN101770027B (zh) * | 2010-02-05 | 2012-05-16 | 河海大学 | 基于InSAR与GPS数据融合的地表三维形变监测方法 |
CN102608584B (zh) * | 2012-03-19 | 2014-04-16 | 中国测绘科学研究院 | 基于多项式反演模型的时间序列InSAR形变监测方法及装置 |
JP6349937B2 (ja) * | 2014-05-09 | 2018-07-04 | 日本電気株式会社 | 変動検出装置、変動検出方法および変動検出用プログラム |
CN104091064B (zh) * | 2014-07-02 | 2017-02-22 | 北京航空航天大学 | 基于优化解空间搜索法的PS‑DInSAR地表形变测量参数估计方法 |
-
2016
- 2016-03-29 CN CN201610187369.7A patent/CN105866777B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105866777A (zh) | 2016-08-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105866777B (zh) | 多角度多时段导航卫星双基地PS-InSAR三维形变反演方法 | |
CN107917880B (zh) | 一种基于地基云图的云底高度反演方法 | |
CN105137430B (zh) | 一种前视阵列sar的回波稀疏获取及其三维成像方法 | |
CN103472450B (zh) | 基于压缩感知的非均匀空间构形分布式sar动目标三维成像方法 | |
CN105445730A (zh) | 一种基于角度分集的海洋流场反演星载sar系统及其方法 | |
CN105467370A (zh) | 一种复合双基地雷达进动目标isar图像横向定标方法 | |
CN109613583A (zh) | 基于单星与地面站测向及联合测时差的无源目标定位方法 | |
CN103698764A (zh) | 一种稀疏采样条件下的干涉合成孔径雷达成像方法 | |
CN108107429A (zh) | 基于最大似然估计的前视超分辨成像方法 | |
CN112904337A (zh) | 一种基于Offset Tracking技术的边坡形变时序监测方法 | |
CN107797110B (zh) | 综合孔径投影辐射的系统成像方法及综合孔径投影辐射计 | |
CN111025256A (zh) | 一种机载雷达的微弱生命体征信号的检测方法及系统 | |
Bishop et al. | Video-SAR using higher order Taylor terms for differential range | |
CN103091682B (zh) | 基于时频分析InISAR多动目标成像和运动轨迹重建法 | |
CN107144815B (zh) | 一种基于一维测向的三维定位方法 | |
CN107561534B (zh) | 一种基于全极化高轨sar的电离层时变tec测量方法 | |
CN110927765B (zh) | 激光雷达与卫星导航融合的目标在线定位方法 | |
CN104833971A (zh) | 基于滑动散射中心的双/多基地雷达图像ps点关联方法 | |
CN110618409B (zh) | 顾及叠掩及阴影的多通道InSAR干涉图仿真方法及系统 | |
KR102151362B1 (ko) | 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법 | |
CN110276240B (zh) | 一种sar图像建筑物墙面窗户信息提取方法 | |
Aloi et al. | Comparative analysis of single-channel direction finding algorithms for automotive applications at 2400 MHz in a complex reflecting environment | |
CN103869312A (zh) | 一种连续场景的稀疏阵列sar侧视三维成像方法 | |
CN107526079A (zh) | 一种基于l型三天线干涉处理的空间自旋目标宽带雷达三维成像方法 | |
CN106371096B (zh) | 机载双天线InSAR三维构像模型构建方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |