SARScape:如何根据 PS-InSAR 结果估计速度
平均速度是SARscape PS处理的输出参数。可在PS-InSAR处理后的.csv文件中找到。
平均速度是根据速度范围参数在反演过程中找到的解。
该技术应用“蛮力”算法,通过尝试不同的速度值来找到与多时相干涉图相位的最佳拟合。
基本上,第一次反演期间求解的方程组如下:(根据 PSI 的线性模型)
Phase1 = v*Time + K_fact*Residual_Height
Phase2 = v*Time + K_fact*Residual_Height
Phase3 = v*Time + K_fact*Residual_Height
Phase4 = v*Time + K_fact*Residual_Height
Phase5 = v*Time + K_fact*Residual_Height
.
.
.
Phase(N) = v*Time + K_fact*Residual_Height
其中 N 是干涉图的数量,v 是平均速度(毫米/年),Residual_Height(米)是相对于输入 DEM 的高度修正,K_fact是一个几何参数,取决于基线、入射角和倾斜距离。
对于每个像素,找到一对(v,Residual_Height)。
与两个未知数一起生成的多时相相干性(第一次和第二次反演步骤:_cc_first,_cc)是衡量所找到的对的质量指标。
该值越高,干涉图相位与线性模型匹配的精度越高。
详细信息请参阅SARscape PS处理所依据的论文:
Ferretti, Alessandro, Claudio Prati, and Fabio Rocca. "Permanent scatterers in SAR interferometry." IEEE Transactions on geoscience and remote sensing 39.1 (2001): 8-20.