[InSAR] 第一篇博客,总结这次项目实验吧

呼呼,做了四天终于做完了80多景影像,总共150多个结果,第一次做项目,总结一下这次所学!
实验总内容:DGB地区14-18年所有数据两两/三三差分干涉
数据:Sentinal-1A/B

#1数据下载
第一次下这么多景哨兵数据,一开始先使用freedownload进行下载,下载完14-15的1B数据全部解压失败(原因都是条带的tiff文件损坏)。之后经过波哥指点,将失败数据全部转用渣雷下载。最后终于数据全部搞定。

#2聚焦轨道矫正
接着批处理脚本S1_SLC_batch_pro对.safe文件进行聚焦并轨道矫正。
参数只需轨道文件orbit即可,此时影像压缩文件不需要解压。
注:批处理脚本会导致寻找轨道文件失败,矫正失败文件都存在noorbit中,需要批处理完成后使用独立脚本对失败文件进行轨道矫正。

#3数据裁剪
根据强度影像获取研究区大概位置后,首先使用独立脚本对swath和burst进行裁剪。要确保每个swath中至少存在两个swath,才可保证后续的配准工作。这里的难点就是对不同的frame和swath间进行裁剪,需要进行一些实验。有时同一个frame中会存在10个burst的情况,需要加大注意,不可直接使用循环进行批裁剪。

#4配准
使用S1_coreg_batch_pro进行批处理
S1_coreg_batch_pro ‘master’ ‘dem’ - - - - - - - -
配准可能导致失败,先看配准质量文件,检查方位向距离向的配准标准差
如果没有质量文件生成,则会在diff.int里发生跳变,进一步检查。确定最终的文件

#5行列裁剪
接着进行行列号的裁剪,进一步缩小区域。
SLC_copy_all
还应保留一些相对稳定的区域(如果存在的话),这样方便之后解缠参考点的选取。

#6DInSAR处理
最后进行DInSAR处理,实验选区相干掩模阈值为0.3。已是最低,不可再低。

#循环shell
for ifg in cat para.txt
do
M_P=echo $ifg | cut -c 1-8
S_P=echo $ifg | cut -c 9-16
(#SLC_copy)
(#DInSAR)
done

#获取参数shell
cat csv | cut -c 73-80 | sed ‘s/[[:space:]]//g’ > para1.txt
cat frame | sed ‘s/[[:space:]]//g’ > para2.txt
paste -d_ para1.txt para2.txt > para.txt
rm para1.txt
rm para2.txt

放一点渣渣结果:
[InSAR] 第一篇博客,总结这次项目实验吧_第1张图片
[InSAR] 第一篇博客,总结这次项目实验吧_第2张图片
两个图的colorbar还没来得及调。

呼呼,瞎G2写了一点点总结,这次也算是更加熟练了吧。之前跑过StaMPS,IPTA,SBAS和少景的DInSAR处理,但都只是停留在实验的层面上,并不能扩充到生产的水平,这次也算是熟悉了一下生产的操作,尽管最后结果一般(aaaa主要是山区相干性太差了)。还是学到了挺多的!!不管是从批处理(裁剪 配准 ),还是到之后问题影像的检查等等等等。
还有就是干涉命令 SLC_intf (rsf 以及 azsf)这两个参数还是默认为好,改成0后干涉图会出现一道一道的暗条纹,具体原因不明。

最后在总结一个学到的命令:
mk_kml $(编码好的影像) dem_seg_par .kml #用来生成kml

最后最后最后 IG牛逼!!!!

你可能感兴趣的:(InSAR)