袁煜伟, 李志伟, 母岷政. 中国第一代全球大气再分析数据在InSAR大气改正中的应用研究:以南加州地区为例[J]. 武汉大学学报 ( 信息科学版). DOI: 10.13203/j.whugis20230154
引用本文: 袁煜伟, 李志伟, 母岷政. 中国第一代全球大气再分析数据在InSAR大气改正中的应用研究:以南加州地区为例[J]. 武汉大学学报 ( 信息科学版). DOI: 10.13203/j.whugis20230154
YUAN Yuwei, LI Zhiwei, MU Minzheng. Application of China's First Generation Global Atmospheric Reanalysis Data in InSAR Atmospheric Correction: A Case Study of Southern California[J]. Geomatics and Information Science of Wuhan University. DOI: 10.13203/j.whugis20230154
Citation: YUAN Yuwei, LI Zhiwei, MU Minzheng. Application of China's First Generation Global Atmospheric Reanalysis Data in InSAR Atmospheric Correction: A Case Study of Southern California[J]. Geomatics and Information Science of Wuhan University. DOI: 10.13203/j.whugis20230154

中国第一代全球大气再分析数据在InSAR大气改正中的应用研究:以南加州地区为例

Application of China's First Generation Global Atmospheric Reanalysis Data in InSAR Atmospheric Correction: A Case Study of Southern California

  • 摘要: 大气延迟是合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)地表形变监测的主要误差源之一。提出采用中国第一代全球大气再分析业务系统(CRA40)数据改正InSAR对流层延迟,并以南加州SBAS-InSAR地表变形监测为例,进行了详细分析与评估。首先提出了利用该产品计算对流层延迟的方法,通过顾及大气参数垂直分层及水平变化的物理特性,分别对该产品原始气象参数进行垂直向和水平向的插值,沿卫星视线方向积分计算大气延迟,通过与原始干涉图差分,得到大气改正结果;其次将结果从干涉图标准差、空间相关性和相位-高程相关系数三个方面进行误差分析;最后利用Sentinel-1数据进行南加州形变测量实验,将本文方法与经典SBAS-InSAR结果对比,使用GPS数据进行验证。结果表明:(1)校正后的干涉图平均标准差减小了34.7%,65%的干涉图平均标准差减小了20%以上; (2)干涉图空间结构函数的期望平方差显著性下降表明该产品能有效抑制长波大气; (3)相关性的拟合系数变化表明该产品能够有效降低高程影响带来的大气垂直分层分量; (4)说明利用分辨率较低的全球大气模型数据进行大气改正时考虑对流层空间变化的有效性。所得结果说明了CRA产品在InSAR大气校正上的可行性。

     

    Abstract: Objectives: Atmospheric delay is one of the main error sources ininterferometric synthetic aperture radar ( InSAR) surface deformation monitoring. Using the global meteorological model to correct the tropospheric delay has been successfully applied in recent years. In this paper, an InSAR atmospheric correction method is presented, which introduces China's first Generation Global Atmospheric reanalysis data (CRA). Taking into account the physical characteristics of tropospheric meteorological parameters, the method was tested in Los Angeles, Southern California, and evaluated quantitatively. The results show the feasibility of CRA products in InSAR atmospheric correction. Methods: The meteorological parameters such as temperature, water vapor and air pressure are obtained by CRA, taking into account the spatial variation characteristics of atmospheric parameters, according to the vertical stratification characteristics of meteorological data, the piecewise function is used to interpolate in the vertical direction. By considering the spatial variability of the atmosphere, the interpolation is carried out in the horizontal plane. According to the results, the atmospheric delay is calculated along the satellite line of sight. The atmospheric correction result is obtained by difference with the original interferogram. Verify and analyze the results, first, the error analysis is carried out from three aspects: Interferogram standard deviation, spatial correlation and phaseelevation correlation coefficient; second, it is compared with the results of GACOS and PyAPS methods; Finally, 78 interferograms obtained from 40 scenes Sentinel-1 data of Southern California based on short baseline principle are used to carry out sequential deformation measurement experiments, the proposed method is compared with classical SBAS-InSAR results, and verified by GPS data. Results: The results show that: (1) The average standard deviation of the corrected interferogram is reduced by 34.7%, The average standard deviation of 65% interferograms has been reduced by more than 20%., especially for the interferograms which are seriously affected by the atmosphere. It is better than GACOS and PYAPS with an average improvement of 30%. (2) The significant decrease of the expected square variance of the spatial structure function of the interferogram shows that the product can effectively suppress the long-wave atmosphere. (3) The change of fitting coefficient of correlation shows that the product can effectively reduce the atmospheric vertical stratification component caused by the influence of elevation. Depending on the degree of atmospheric influence, CRA can improve the vertical stratification component of the atmosphere, ranging from 20% to 60% or more. (4) It is effective to take the spatial variation of convection into account when establishing a high-resolution InSAR tropospheric delay map using low-resolution global atmospheric model data. Conclusions: This paper confirms that CRA can correct and improve the overall accuracy of deformation monitoring. Through quantitative evaluation, researchers can fully understand the correction effect and performance of this product, which is helpful to promote the application and development of this product. However, due to the differences of external data, follow-up research can combine GNSS data with CRA solution to further improve the accuracy of InSAR atmospheric correction.

     

/

返回文章
返回