Loading [MathJax]/jax/output/SVG/jax.js

基于导航卫星的干涉SAR数据采集策略优选方法分析

张凌志 刘飞峰 胡程

张凌志, 刘飞峰, 胡程. 基于导航卫星的干涉SAR数据采集策略优选方法分析[J]. 雷达学报, 2019, 8(5): 624–630. doi: 10.12000/JR19065
引用本文: 张凌志, 刘飞峰, 胡程. 基于导航卫星的干涉SAR数据采集策略优选方法分析[J]. 雷达学报, 2019, 8(5): 624–630. doi: 10.12000/JR19065
ZHANG Lingzhi, LIU Feifeng, and HU Cheng. Optimization method and analysis of data acquisition strategy based on interference SAR with GNSS transmitters[J]. Journal of Radars, 2019, 8(5): 624–630. doi: 10.12000/JR19065
Citation: ZHANG Lingzhi, LIU Feifeng, and HU Cheng. Optimization method and analysis of data acquisition strategy based on interference SAR with GNSS transmitters[J]. Journal of Radars, 2019, 8(5): 624–630. doi: 10.12000/JR19065

基于导航卫星的干涉SAR数据采集策略优选方法分析

DOI: 10.12000/JR19065 CSTR: 32380.14.JR19065
基金项目: 国家自然科学基金项目(61601032, 61625103)
详细信息
    作者简介:

    张凌志(1994–),男,四川成都人,硕士研究生。2017年在东北大学计算机科学与工程学院获得学士学位,现为北京理工大学信息与电子学院硕士研究生。主要研究方向为SS-BiSAR信号处理,差分干涉SAR等。E-mail: zlz_bit@outlook.com

    刘飞峰(1983–),男,陕西渭南人,博士,副教授。2012年在北京理工大学信息与电子学院获得博士学位,现担任北京理工大学信息与电子学院副教授。主要研究方向为双/多基地SAR信号处理、GEOSAR成像算法和SS-BiSAR目标检测等。E-mail: feifengliu_bit@bit.edu.cn

    胡 程(1981–),男,湖南岳阳人,博士,教授。2009年在北京理工大学信息与电子学院获得博士学位,现担任北京理工大学信息与电子学院教授,博士生导师,主要研究方向为地球同步轨道SAR、双基地/前向散射雷达信号处理和昆虫雷达等。E-mail: hucheng.bit@gmail.com

    通讯作者:

    刘飞峰 feifengliu_bit@bit.edu.cn

  • 责任主编:张晓玲 Corresponding Editor: ZHANG Xiaoling
  • 中图分类号: TP79

Optimization Method and Analysis of Data Acquisition Strategy Based on Interference SAR with GNSS Transmitters

Funds: The National Natural Science Foundation of China (61601032, 61625103)
More Information
  • 摘要: 基于导航卫星的干涉SAR(GNSS-InSAR)使用在轨导航卫星作为照射源,近地面部署接收机,利用导航卫星的星座特性以及重轨特性,可实现区域性的连续观测。对于场景1维/3维形变反演而言,需要连续时间的数据采集,由于导航卫星并非严格意义上的重轨,且重轨时间具有不确定性,原始数据冗余度高,数据对齐时截取量大,数据有效性低。该文针对GNSS-InSAR场景数据采集时间精确性问题,提出了一种重轨数据采集优化模型,该方法通过实际轨迹与TLE预测轨迹相结合的方式,通过空间相干系数的滑窗轨迹对齐,以获取相邻天导航卫星重轨时间间隔,实现精确的GNSS-InSAR数据采集,在降低原始数据冗余度下,保证数据的有效合成孔径时间。实测数据表明所提方法的有效性。

     

  • 基于导航卫星的双基地SAR(Bistatic Synthetic Aperture Radar based on Global Navigation Satellite System, GNSS-BSAR)是空-地双基地SAR中一种典型的应用[1],使用在轨的导航卫星作为发射源,地面部署接收机(地基、车载、机载)构成双基地SAR系统[2]。由于导航星座的日趋完善,其全球覆盖性以及重轨特性所带来的优势是其他照射源暂时所不能替代的,其中以地基接收机为主的导航卫星干涉合成孔径雷达(Interference Synthetic Aperture Radar based on the Global Navigation Satellite System, GNSS-InSAR)在场景形变监测领域有着广阔的应用前景[3],成为了近年来研究热点。

    在GNSS-BSAR系统成像方面,已有研究者分别使用不同的导航卫星星座进行了成像验证,包括了北斗[4,5]、GPS[6]、格洛纳斯[7]、伽利略[8]。除此以外文献[9]还提出了多角度融合方法以增强图像信噪比。在形变监测方面,来自伯明翰大学的学者们[10]使用直达波天线,配合长约50 m的线缆构建了理想点目标,并使用格洛纳斯作为发射源,首次实现了精度约为1 cm的1维形变反演结果。该实验初步验证了GNSS-InSAR应用于形变监测的可行性。为了进一步验证场景形变监测的可能,2017年文献[11]通过对接收机进行高精度移位来模拟场景建筑形变,成功反演出了形变,精度约为1 cm。在3维形变方面,2018年北京理工大学的技术团队[12]通过人为构建转发器,进行了精度可控的强点目标形变模拟,使用我国的北斗IGSO卫星,成功实现了精度优于5 mm的3维形变反演,这些验证性实验充分表明了GNSS-InSAR应用于场景形变检测的可能。

    若要实现GNSS-InSAR场景的3维形变反演,需要同时至少3颗卫星从不同角度照射场景。由于GNSS-InSAR系统的拓扑高度非对称性以及导航信号的窄带特性[13],加上导航卫星的重轨并非是严格意义上的重轨,除了不可避免的空间基线外,重轨时间也并非严格一致,因此在实际数据采集中,需要对系统构型以及数据采集时间进行严格的优化设计。文献[14]提出了一种联合优化方法,解决了面向大场景下的多星多角度构型优化问题,配合多个接收机实现综合分辨性能优异的大场景成像。文献[15]提出了空间去相干的理论描述框架,表明了空间去相干在GNSS-InSAR中的必要性,但未对数据采集时间进行说明。从当前实际情况出发,不精确的数据采集时间可能会造成存储资源浪费,空间去相干导致的数据截取进一步降低了数据有效性。具体如图1所示:

    图  1  数据截取与有效数据示意图
    Figure  1.  Effective data interception diagram

    针对上述问题,本文提出了一种GNSS-InSAR场景连续数据采集优化方法,通过结合当前数据的卫星轨迹和两行星历数据文件(STK Two-Line Element sets, TLE)预测轨迹,基于相干系数轨迹对齐,获取卫星重轨时间间隔,得到最优的数据采集策略,从源头上降低数据的空间去相干性,提升所采集数据有效性,节约存储资源。在第2部分对GNSS-InSAR场景数据采集优化方法进行了详细介绍。第3部分针对提出的方法进行了实验设计,开展了实测数据采集,并针对采集的数据进行了初步分析。第4部分对全文进行总结。

    对于GNSS-InSAR图像而言,经过保相成像处理后,场景中任意一点(x0,y0)的像素信息分别对应分辨单元内所有散射体回波的相参叠加,可建模为

    s(x0,y0)=f(x,y;t)exp[j2πλr(x,y;P)]W(xx0,yy0;P)dxdy+n(x0,y0;t) (1)

    其中,f(x,y;t)为时间t下的地表散射系数,P为对应的合成孔径中心位置矢量,W(x,y;P)表示系统的点扩散函数(PSF), n(x,y;t)为图像的加性噪声。对于SAR图像的同名点像素,其相干系数可表示为[16]

    ρ=sm(x0,y0)ss(x0,y0)dxdysm(x0,y0)sm(x0,y0)dxdyss(x0,y0)ss(x0,y0)dxdy (2)

    其中,下标m表示主图像,s表示辅图像。根据柯西不等式可以判断:0ρ1,当ρ=0时表示同名点完全不相干,当ρ=1时,同名点完全相干。

    将点目标像素模型式(1)带入到式(2)并化简得到

    ρ=sav(x,y)exp[j2πλ(r(x,y;Pm)r(x,y;Ps))]|W(xx0,yy0;Pm)|2dxdy(sm(x,y)exp[j2πλr(x,y;Pm)]|W(xx0,yy0;Pm)|2dxdy+nm)×(ss(x,y)exp[+j2πλr(x,y;Ps)]|W(xx0,yy0;Pm)|2dxdy+ns) (3)

    式(3)的推导使用了如下近似:

    (1) 由于导航卫星的高轨道特性,使得W(xx0,yy0;Pm)W(xx0,yy0;Ps)成立;

    (2) 相邻两天的噪声相干系数为0,即

    n(x0,y0;tm)×n(x0,y0;ts)=0 (4)

    (3) 相邻两天的目标散射系数为sav(x,y),即

    sav(x,y)f(x,y;tm)f(x,y;ts) (5)

    对式(3)中的相干系数ρ进一步分解得到

    ρ=ρth×ρti×ρsp (6)

    其中,热噪声相干系数ρth与时间相干系数ρti分别由系统与实际目标决定。

    对于PS点[17]而言,地表散射系数相对稳定,不随时间变化,同时为了便于后续分析,假定散射系数为1得到空间相干系数ρsp的简化式为

    ρsp=exp[j2πλ(r(x,y;Pm)r(x,y;Ps))]|W(xx0,yy0;Pm)|2dxdy(|W(xx0,yy0;Pm)|2dxdy)(|W(xx0,yy0;Pm)|2dxdy) (7)

    从式(7)推导结果可以知道,空间基线主要是影响r(x,y;P)从而导致空间去相干。

    导航卫星的重轨时间并非稳定不变,因此需要对数据采集时间进行有效预测,从源头上降低空间去相干,提高数据有效性。

    假定主图像数据采集时间为tm,该采集时间可以通过文献[14]中的广义优化模型进行求解,辅图像数据采集时间为ts=tm+Δt, Δt为时间间隔,那么最优化数据采集模型可通过式(7)推导而来

    Δt=argmax{|˜W(x,y;Pm)|2exp[j2πλ(r(x,y;Pm)r(x,y;P(tm+Δt)))]dxdy} (8)

    其中,˜W(x,y;Pm)tm下等效归一化PSF, P(tm+Δt)Δt时间偏置下得合成孔径中心位置矢量。

    第1天数据采集需要进行实验设计以确定最优数据采集时刻,往后的重轨天数据采集可以根据数据采集优化模型,同时结合星历文件进行预测。整体的预测流程如图2所示,n为任意一天采集的数据,k为重轨天数间隔。

    图  2  GNSS-InSAR数据采集时间优化流程
    Figure  2.  Time optimization process of GNSS-InSAR data acquisition

    实际卫星位置对应的实际时间设为tn,经过模型优化得到的时间偏差为Δt,那么第n+k天对应的实际数据采集时间可表示为

    tn+k=tn+Δt (9)

    对于固定场景的形变监测,首次数据采集的时候需要严格设计系统构型,使分辨率达到最优化。本次实验接收机部署在北京理工大学信息科学试验楼楼顶西北角,实施监测场景位于西偏北30°。使用理论分辨率计算公式[18]对该场景进行分辨率设计。仿真参数具体参见表1

    表  1  数据采集试验仿真参数
    Table  1.  Data acquisition test simulation parameters
    参数
    照射源北斗 IGSO1~5
    PRF1000 Hz
    带宽10.23 MHz
    合成孔径时间600 s
    TLE文件更新日期2019年4月29日
    预定数据采集日期2019年4月30日
    下载: 导出CSV 
    | 显示表格

    以分辨单元面积作为判定依据,得到预定采集日期当天全时段下各个卫星在预定场景下所能得到的分辨单元面积如图3所示。

    图  3  全时段下场景分辨单元面积
    Figure  3.  Scene resolution unit area in full time

    为了实现3维形变反演,需要同一时间下有3颗卫星对场景进行照射。图3中10点前后与17点前后满足当前场景上空有3颗IGSO卫星可见的条件。更进一步,为了使分辨单元面积达到最优,可以得到具体的数据采集时间。具体如图4红框标注,分别是9点30分前后与17点30分前后。

    图  4  首次数据采集时间设计结果
    Figure  4.  Design results of first data acquisition time

    为了配合实验,在场景布置转发器,整体的系统构型如图5所示。

    图  5  GNSS-InSAR场景3维形变反演实验拓扑构型设计结果
    Figure  5.  GNSS-InSAR topological configuration design results of 3D deformation retrieval experiment

    以2019年4月30日采集的实测数据作为第n天数据,对于北斗的IGSO而言,重轨时间约为1天,即m=1,同时下载当天最新的TLE文件。以IGSO1为例,结合图2进行详细说明:

    (1) 使用实测数据的直达波进行卫星位置解算,同时根据TLE文件推算当天和相邻天的卫星轨迹。经过相干系数轨迹匹配之后,得到的轨迹如图6所示。

    图  6  对齐后的TLE卫星轨迹与实测数据卫星轨迹
    Figure  6.  Aligned TLE satellite trajectory and measured data satellite trajectory

    (2) 以匹配得到的TLE卫星轨迹作为参考,对重轨天的TLE卫星轨迹进行数据采集优化模型求解,系统的PSF与优化模型仿真结果分别如图7图8所示。

    图  7  场景[–147, 20, 0]处理论PSF
    Figure  7.  Theoretical PSF in scene at position of [–147, 20, 0]
    图  8  数据采集优化模型仿真结果
    Figure  8.  Simulation results of data acquisition optimization model

    图8的结果分析可知,第1个峰值点为其本身,由于空间基线为0,相干系数为1。第2个峰值点相干系数为0.999644,满足除了第1个峰值点外相邻天相干系数最大值条件,因此第2个峰值点就是最佳重轨时的空间相干系数。此时经过模型优化得到的时间间隔为:Δt=86163s=23h56min3s,结合第1天的实测数据轨迹对应的时间t1=9h26min0s,第2天准确的数据采集时间为:t2=9h22min3s

    为了说明优化结果的正确性,在实验场景中放置转发器模拟理想点目标(图5),同时按照优化后的时间进行5月1日数据采集。实际采集时间为9h21min53s,总采集时间约650 s。相邻两天的空间相干系数轨迹匹配结果如图9所示。

    图  9  实测数据重轨空间相干系数
    Figure  9.  The spatial coherence coefficient of measured data

    图9中峰值点位置来看,重轨数据采集优化模型得到的结果和实际结果相吻合。为了进一步说明,图10给出了IGSO1卫星实测数据成像结果。

    图  10  场景成像结果
    Figure  10.  Imaging results of scene

    对相邻两天的图像相干系数进行求解,得到图11所示结果。在同一坐标系下,仿真目标位于[–147, 20, 0],空间相干系数为0.999644;转发器位于[–147, 20, 0],相干系数为0.9996;两者的相干系数基本保持一致。

    图  11  相干系数结果
    Figure  11.  Coherence coefficient result

    图9图10的结果表明经过数据采集优化模型后得到的时间间隔与实际卫星轨迹的重轨时间相互吻合,在保证600 s预期合成孔径时间下,可以最大限度减少数据采集时间,节约存储资源。同时避免后期由于数据对齐带来的数据有效性降低问题。

    在GNSS-InSAR场景1维/3维形变反演应用中,针对由于导航卫星重轨时间的非严格一致性与有效数据截取带来的数据冗余,数据有效性低等问题,本文提出了一种面向GNSS-InSAR场景数据采集的优化模型,采用实测数据与TLE文件相结合,根据当天数据采集时间,预测相邻天重轨时间,从而实现精确的数据采集。实测数据验证结果表明了数据采集时间优化模型的正确性。该方法的提出有利于GNSS-InSAR场景1维/3维形变反演实验的开展,在降低原始数据冗余度基础上,保证了有效数据时间长度大于预期合成孔径时间。

  • 图  1  数据截取与有效数据示意图

    Figure  1.  Effective data interception diagram

    图  2  GNSS-InSAR数据采集时间优化流程

    Figure  2.  Time optimization process of GNSS-InSAR data acquisition

    图  3  全时段下场景分辨单元面积

    Figure  3.  Scene resolution unit area in full time

    图  4  首次数据采集时间设计结果

    Figure  4.  Design results of first data acquisition time

    图  5  GNSS-InSAR场景3维形变反演实验拓扑构型设计结果

    Figure  5.  GNSS-InSAR topological configuration design results of 3D deformation retrieval experiment

    图  6  对齐后的TLE卫星轨迹与实测数据卫星轨迹

    Figure  6.  Aligned TLE satellite trajectory and measured data satellite trajectory

    图  7  场景[–147, 20, 0]处理论PSF

    Figure  7.  Theoretical PSF in scene at position of [–147, 20, 0]

    图  8  数据采集优化模型仿真结果

    Figure  8.  Simulation results of data acquisition optimization model

    图  9  实测数据重轨空间相干系数

    Figure  9.  The spatial coherence coefficient of measured data

    图  10  场景成像结果

    Figure  10.  Imaging results of scene

    图  11  相干系数结果

    Figure  11.  Coherence coefficient result

    表  1  数据采集试验仿真参数

    Table  1.   Data acquisition test simulation parameters

    参数
    照射源北斗 IGSO1~5
    PRF1000 Hz
    带宽10.23 MHz
    合成孔径时间600 s
    TLE文件更新日期2019年4月29日
    预定数据采集日期2019年4月30日
    下载: 导出CSV
  • [1] CHERNIAKOV M. Space-surface bistatic synthetic aperture radar—prospective and problems[C]. Proceedings of RADAR 2002, Edinburgh, UK, 2002.
    [2] ANTONIOU M and CHERNIAKOV M. Experimental demonstration of passive GNSS-based SAR imaging modes[C]. Proceedings of IET International Radar Conference 2013, Xi’an, China, 2013.
    [3] ZENG Tao, ZHANG Tian, TIAN Weiming, et al. A novel subsidence monitoring technique based on space-surface bistatic differential interferometry using GNSS as transmitters[J]. Science China Information Sciences, 2015, 58(6): 1–16.
    [4] ZENG Tao, ZHANG Tian, TIAN Weiming, et al. Bistatic SAR imaging processing and experiment results using BeiDou-2/Compass-2 as illuminator of opportunity and a fixed receiver[C]. Proceedings of the 2015 IEEE 5th Asia-Pacific Conference on Synthetic Aperture Radar, Singapore, Singapore, 2015.
    [5] FAN Xuezhen, LIU Feifeng, ZHANG Tian, et al. Passive SAR with GNSS transmitters: Latest results and research progress[C]. Proceedings of 2017 IEEE International Geoscience and Remote Sensing Symposium, Fort Worth, USA, 2017.
    [6] ANTONIOU M, CHERNIAKOV M, and MA Hui. Space-surface bistatic synthetic aperture radar with navigation satellite transmissions: A review[J]. Science China Information Sciences, 2015, 58(6): 1–20.
    [7] LIU Feifeng, ANTONIOU M, ZENG Zhangfan, et al. Point spread function analysis for BSAR with GNSS transmitters and long dwell times: Theory and experimental confirmation[J]. IEEE Geoscience and Remote Sensing Letters, 2013, 10(4): 781–785. doi: 10.1109/LGRS.2012.2223655
    [8] MA Hui, ANTONIOU M, and CHERNIAKOV M. Passive GNSS-based SAR resolution improvement using joint Galileo E5 signals[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(8): 1640–1644. doi: 10.1109/LGRS.2015.2417594
    [9] ZENG Tao, ZHANG Tian, TIAN Weiming, et al. Space-surface bistatic SAR image enhancement based on repeat-pass coherent fusion with Beidou-2/Compass-2 as illuminators[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(12): 1832–1836. doi: 10.1109/LGRS.2016.2614337
    [10] ANTONIOU M, LIU F, ZENG Z, et al. Coherent change detection using GNSS-based passive SAR: First experimental results[C]. Proceedings of 2012 IET International Conference on Radar Systems, Glasgow, UK, 2012.
    [11] TZAGKAS D, ANTONIOU M, and CHERNIAKOV M. Coherent change detection experiments with GNSS-based passive SAR[C]. Proceedings of 2016 European Radar Conference, London, UK, 2016.
    [12] LIU Feifeng, FAN Xuezhen, ZHANG Tian, et al. GNSS-based SAR interferometry for 3-D deformation retrieval: Algorithms and feasibility study[J]. IEEE Transactions on Geoscience and Remote Sensing, 2018, 56(10): 5736–5748.
    [13] DONG Xichao, HU Cheng, TIAN Weiming, et al. Feasibility study of inclined geosynchronous SAR focusing using Beidou IGSO signals[J]. Science China Information Sciences, 2016, 59(12): 129302. doi: 10.1007/s11432-016-5524-x
    [14] LIU Feifeng, FAN Xuezhen, ZHANG Lingzhi, et al. GNSS-based SAR for urban area imaging: Topology optimization and experimental confirmation[J]. International Journal of Remote Sensing, 2019, 40(12): 4668–4682. doi: 10.1080/01431161.2019.1569790
    [15] ZHANG Qilei, ANTONIOU M, CHANG Wenge, et al. Spatial decorrelation in GNSS-based SAR coherent change detection[J]. IEEE Transactions on Geoscience and Remote Sensing, 2015, 53(1): 219–228. doi: 10.1109/TGRS.2014.2321145
    [16] ZEBKER H A and VILLASENOR J. Decorrelation in interferometric radar echoes[J]. IEEE Transactions on Geoscience and Remote Sensing, 1992, 30(5): 950–959. doi: 10.1109/36.175330
    [17] FERRETTI A, PRATI C, and ROCCA F. Permanent scatterers in SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(1): 8–20. doi: 10.1109/36.898661
    [18] ZENG Tao, CHERNIAKOV M, and LONG Teng. Generalized approach to resolution analysis in BSAR[J]. IEEE Transactions on Aerospace and Electronic Systems, 2005, 41(2): 461–474. doi: 10.1109/TAES.2005.1468741
  • 加载中
图(11) / 表(1)
计量
  • 文章访问数: 
  • HTML全文浏览量: 
  • PDF下载量: 
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-07-04
  • 修回日期:  2019-08-12
  • 网络出版日期:  2019-10-01

目录

/

返回文章
返回