1、win10中安装linux运行环境,相比虚拟机和双系统,该方法更有效,在windows应用商店搜索ubuntu,这里选择ubuntu20.
该运行环境不需要挂接硬盘,通过cd /mnt/c,就可以访问c盘。
2、ubuntu中安装Anaconda3环境,鉴于安装速度,可以把ubuntu和Anaconda3的安装源设置为国内镜像,教程:https://mirrors.tuna.tsinghua.edu.cn/help/anaconda/,在Anaconda3中创建python3环境,isce2为python开发。安装教程参考:https://zhuanlan.zhihu.com/p/269183148。可用python --version查看python是否创建。
3、安装isce2,这里通过在线安装,用过本地安装,发现有些插件版本不匹配会导致安装失败,推荐在线安装,保证好网速,一次不成功重复多次。安装教程也参考:https://zhuanlan.zhihu.com/p/269183148。安装后,用topsApp.py -h测试是否成功。
4、利用isce2进行干涉,获取地表形变,首先准备三个文件,注意isce2.5中主从影像不用master、slave表示,而用reference、secondary表示。其中dem需要在线下载,使用函数dem.py,需设置账号到配置文件,
machine urs.earthdata.nasa.gov
login ######
password ######
教程参考:
https://blog.csdn.net/ZB18810657133/article/details/116271087?utm_medium=distribute.pc_relevant.none-task-blog-2%7Edefault%7ECTRLIST%7Edefault-1.no_search_link&depth_1-utm_source=distribute.pc_relevant.none-task-blog-2%7Edefault%7ECTRLIST%7Edefault-1.no_search_link,
https://blog.csdn.net/weixin_42128764/article/details/108602917。文件如下:
topsApp.xml
SENTINEL1 [3] 7 3 [30, 31, 110, 111] reference.xml secondary.xml True snaphu_mcf dem/demLat_N30_N31_Lon_E110_E111.dem.wgs84 dem/demLat_N30_N31_Lon_E110_E111.dem.wgs84
reference.xml
/insar/raw/S1A_IW_SLC__1SDV_20210315T103600_20210315T103628_037008_045ADC_F43C.SAFE reference orbit
secondary.xml
/raw/S1A_IW_SLC__1SDV_20210420T103601_20210420T103629_037533_046D04_EA5F.SAFE secondary orbit
然后在Ubuntu中输入命令:topsApp.py topsApp.xml,开始运行程序,结果会得到多个文件夹.
其中merged为最终结果所在。形变结果为文件filt_topophase.unw.geo.vrt,将波段切回2波段即可看到地表形变结果,单位为mm。
其中强度图为topophase.cor.geo.vrt,类似光学影像中全色波段,分辨率20米。
5、isce2+stamps是主流的开源insar处理模式,前者生成干涉对,后者进行时序分析,可配合使用。



