利用ISCE2处理时序COSMO-SkyMed数据

一、COSMO-SkyMed数据

COSMO-SkyMed卫星是意大利航天局和意大利国防部共同研发的COSMO-SkyMed高分辨率雷达卫星星座的第二颗卫星。COSMO-SkyMed卫星是X波段合成孔径雷达（SAR）卫星，波长约3.11cm。

二、ISCE2 时序处理

ISCE2（InSAR Scientific Computing Environment version 2）是常用的开源InSAR数据处理软件，其为Sentinel-1、ALOS、TerraSAR-X等多个卫星提供了D InSAR、时序InSAR的处理流程。
ISCE2中，单独的一个名为Stack Processors的模块是专门用于处理时序数据的。该模块为Sentinel-1的tops数据和ALOS数据提供了单独的处理模块，而其他聚束模式的数据都用stripmapStack模块进行处理。
ISCE2 GitHub 地址：https://github.com/isce-framework/isce2
ISCE2 时序处理说明：https://github.com/isce-framework/isce2/tree/main/contrib/stack

三、过程记录

1. 使用dem.py准备DEM

dem.py -a stitch -b -37 -31 -72 -69 -r -s 1 -c

3. 将所有数据解压到单独文件中，这里将这些文件保存到rawdata文件夹中
4. 使用prepRawCSK.py脚本生成转换命令run_unPackCSK

prepRawCSK.py -i ./rawdata/ -o SLC

4. 这里直接用生成的文件有问题，我将其进行了批量修改

unpackFrame_CSK.py -i /home/ISCE/rawdata/20190904 -o /home/ISCE/SLC/20190904

然后，执行执行

chmod +x run_unPackCSK
bash run_unPackCSK

程序将把原始文件转换为ISCE格式的SLC文件

5. 使用stackStripMap.py脚本生成运行文件

 stackStripMap.py -s ./SLC/ -m 20190904 -d ./demLat_N39_N40_Lon_E116_E117.dem.wgs84 -w ./work -t 150 -b 1000 -u snaphu -W interferogram --nofocus

具体参数含义可以通过stackStripMap.py -h查看

6. 运行脚本

chmod +x ./work/run_files/
bash /home/ISCE/work/run_files/run_01_reference
bash /home/ISCE/work/run_files/run_02_focus_split
bash /home/ISCE/work/run_files/run_03_geo2rdr_coarseResamp
bash /home/ISCE/work/run_files/run_04_refineSecondaryTiming

前四步问题不大，一步一步执行，第四步耗时较长，建议后台运行。

7. 修改第四步的结果
   第四步做完之后，在 work/refineSecondaryTiming/pairs/ 目录下将生成所有干涉对的中间结果
   但是这个结果中可能存在部分错误的结果
   错误的结果中misreg.dir文件只有一行，而正常的文件有四行
   将所有的misreg.dir文件的大小输出出来就可以手动找到出问题的文件

cd work/refineSecondaryTiming/pairs/
ll *_*/misreg.dir > pair_list
## 找到错误的结果
cat pair_list | while read l; do rm -rf ${l} ; done
rm -rf pair_list

8. 继续执行其他的命令

bash /home/ISCE/work/run_files/run_05_invertMisreg
bash /home/ISCE/work/run_files/run_06_fineResamp
bash /home/ISCE/work/run_files/run_07_grid_baseline
bash /home/ISCE/work/run_files/run_08_igram

第八步执行完之后，可以在work/Igrams/ 文件夹中看到所有的干涉对包括其解缠结果
第八步比较耗时，建议后台执行

P.S. 在GitHub的ISCE2的Issues中有各种各样的问题和回答，可以在这里找到很多处理方法

ISCE2 Issues Website