加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

SARproZ自定义影像连接中干涉对剔除

(2018-06-11 17:50:12)
标签:

qps

干涉对剔除

分类: 软件操作及更新

针对于影像数量较少的数据集,常采用SBAS-InSAR技术进行处理,SARproZ提供了另一种解决思路——QPS技术。在经典PS-InSAR的基础上,采用不同的干涉对组合策略,将所有获得的SAR影像数据根据干涉图中相干性值作为阈值组合成若数据集。这样克服了空间去相关问题,并且提高了观测数据的时间采样率,充分使用了更多的SAR数据。例如,将12景影像数据集,根据一定得相干性阈值,组合成41个干涉对,大大提高了观测数据得时间采样率。

http://s3/mw690/006PeO3vzy7lbaOudPkd2&690

在实际处理中,存在干涉对效果不好现象,需要剔除,但软件未提供自动剔除干涉对操作,需进行手动剔除。首先,将保存的干涉对自定义连接文件(工程目录下RESULTS/MATLAB/Graph_0260loadmatlab

http://s8/mw690/006PeO3vzy7lbaQI2NNf7&690

其中,CoerGrafoImg为干涉对的相干性,第一列为第一个干涉对的相干性值。UsedImg为影像列表20180101代表影像序号1MatrIncImg为干涉对组合文件,第一列为干涉对连接序号,第二列为影像序号,(11)、(12)代表第1个干涉对连接,连接的影像为第1个影像和第2个影像。

http://s3/mw690/006PeO3vzy7lbaTqNGyb2&690

http://s6/mw690/006PeO3vzy7lbaTHHSJ05&690

http://s4/mw690/006PeO3vzy7lbaTsKwH33&690

如果想删除连接影像120180101)和影像220180112)的第1个干涉对连接,在matlab中需将MatrIncImg11=0MatrIncImg12=0赋值为0

http://s13/mw690/006PeO3vzy7lbaVCOL28c&690

http://s9/mw690/006PeO3vzy7lbb2b45af8&690

赋值完后,已将(11=1和(12=1删除。

http://s2/mw690/006PeO3vzy7lbb7OtVf11&690

由于MatrIncImgsparse矩阵,需将其转换为普通矩阵,即a=full(MatrIncImg),aMatrIncImg转换后的矩阵。

http://s13/mw690/006PeO3vzy7lbb9xCLq7c&690

可以看到a是由10-1组成的矩阵,其中行代表干涉对序列号,列代表影像序列号。例如:第2行中第1列和第3列分别为1-1,表明在第1个影像(20180101)和第3个影像(20180123)组成了第2个干涉对。由于我们删除了第1个干涉对,所以第1行全为0.

http://s13/mw690/006PeO3vzy7lbbbiEWMdc&690
删除a矩阵第1行,a变为4012列矩阵。

http://s11/mw690/006PeO3vzy7lbbedkgyaa&690

matlab中输入MatrIncImg=sparsea,a矩阵值覆盖MatrIncImg矩阵。

http://s14/mw690/006PeO3vzy7lbbgBiX34d&690

可以看到MatrIncImg也变为4012列矩阵,同时删除干涉对对应的相干性值,即将CoerGrafoImg1列删除,变为40列。

http://s12/mw690/006PeO3vzy7lbbl06AH2b&690
最后将矩阵a删除,Save Workspace保存最终的文件。

http://s13/mw690/006PeO3vzy7lbbofmDabc&690

将保存的干涉对连接文件加载,可以看到20180101-20180112干涉对连接被删除,最终干涉对连接数为40

http://s9/mw690/006PeO3vzy7lbbpP5lS38&690

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有