ISCE2StaMPS 全过程
最编程
2024-03-19 14:23:31
...
1.准备文件夹
准备四个文件夹:AxuDir 、DEM 、Orbits 、Process、 SLC
2.准备 dem
cd DEM
dem.py -a stitch -b 38 39 116 118 -r -s 1 -c
rm demLat*.dem demLat*dem.xml demLat*.dem.vrt
cd ..
3.预处理
#stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N36_N39_Lon_E100_E104.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '37 38.5 100.3 102' -W slc
#stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N36_N39_Lon_E099_E103.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '36.43 38.4 99.43 102.59' -W slc
#stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N36_N39_Lon_E099_E103.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '36.44 38.4 99.44 102.57' -W slc
cd Process
stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N38_N39_Lon_E116_E118.dem.wgs84 -a ../AuxDir/ -o ../Orbits -b '38.62 39.23 117.1 117.6' -W slc -z 5 -r 20
生成文件夹run_files 、configs 文件 SAFE_files.txt
4.给run_files赋予权限
sudo chmod 777 -R run_files/
cd run_files/
5.对每个文件依次进行处理(如果第一步报错,可能需要把DEM复制到该文件夹下)
./run_01_unpack_topo_reference
./run_02_unpack_secondary_slc
./run_03_average_baseline
./run_04_extract_burst_overlaps
./run_05_overlap_geo2rdr
./run_06_overlap_resample
./run_07_pairs_misreg
./run_08_timeseries_misreg
./run_09_fullBurst_geo2rdr
./run_10_fullBurst_resample
./run_11_extract_stack_valid_region
./run_12_merge_reference_secondary_slc
./run_13_grid_baseline
6.最终生成的所需文件都在merged文件夹下
包括baselines文件夹 geom_reference文件夹 SLC文件夹
7.生成input_file文本文件内容
source_data slc_stack
slc_stack_path /media/estar/es/INSAR/ISCE/Process/merged/SLC
slc_stack_reference 20170320
slc_stack_geom_path /media/estar/es/INSAR/ISCE/Process/merged/geom_reference
slc_stack_baseline_path /media/estar/es/INSAR/ISCE/Process/merged/baselines
range_looks 40
azimuth_looks 10
aspect_ratio 4
lambda 0.056
slc_suffix .full
geom_suffix .full
######## 如果做sbas,在isce2/contrib/timeseries/prepStackToStaMPS/bin/make_single_reference_stack_isce中第二行 setgenerate_ifgs=n,将y改成n
在主文件夹 (与SLC同级别)下运行
make_single_reference_stack_isce # 这步操作会将input_file变为 INSAR_20*** 文件夹
SBAS 处理
一、在INSAR_20*** 文件夹下打开终端
mt_extract_info_isce
二、在此文件夹下打开matlab,在matlab中输入ps_load_info
ps_load_info
生成loading info into matlab
SERTPARM: heading = -10.3323
SERTPARM: lambda = 0.056
三、输入sb_find(0.4, 48, 150) 0.4代表相干性, 48代表最大时间基线48天 150代表最大空间基线150m
sb_find(0.4, 48, 150)
RHO_MIN:最小相干性 (0.50)
DDIFF_MAX: 总去相关的天数(1500)
BDIFF_MAX:以m为单位的相关基线(1070)
四、在INSAR_20*** 文件夹下打开终端
make_small_baselines_isce
五、做多视,在INSAR_20*** 文件夹中放入多视文件 prep_ml_isce.m 然后在matlab中敲入下面的命令
# mt_prep_isce脚本下载地址`https://pan.baidu.com/s/1KWBl18m_xp2xFA8M-Pnehw?pwd=cmih`
prep_ml_isce('401') #此处还存在疑问
修改mt_prep_isce的161行代码,加入SMALL_BASELINES(不做多视就不需要)
七、cd MultiProc文件夹,打开终端,输入
mt_prep_isce 0.6 3 3 50 200
八、在MultiProc文件夹中打开matlab,做sbas
setparm('small_baseline_flag','y');
setparm('unwrap_method','3D_QUICK');
stamps(1,8);
推荐阅读
-
客户关系管理项目管理方案:全过程项目管理、降低成本和控制风险
-
iOS 应用程序钩子反向全过程
-
大力发展装配式建筑、智能化施工技术,推进建筑全过程标准化
-
iOS 越狱全过程,解决 cydia 或 sileo 无法联网问题
-
SBM 模型分析全过程
-
从开发、部署到网站备案全过程的私人工作经验
-
什么是装配式建筑?在一篇文章中了解集成数字交付的全过程
-
广联达:"数字建筑 "将推动建筑业向现代工业化水平迈进--一是全过程、全要素、全参与方的数字化。"数字建筑 "整合了人员、流程、数据、技术和业务系统,对建筑从规划设计到施工建设、运营维护的全生命周期进行管理。 二是数字化、在线化、智能化。这也是数字化建筑的三大典型特征。其中,数字化是基础,在线化是关键,智能化是目标。 三是新设计、新建设、新运维。试想,未来通过全数字化样板设计实现个性化最优方案,通过工业化施工提高效率精益求精,通过智能化运维提升建筑品质低碳宜居,将推动建筑业向现代工业化水平迈进。 广联达的一批标杆项目和应用案例备受关注。
-
打造项目 "全过程 "数字化,为建筑业转型赋能
-
解释手机通信的全过程 #科学 #知识 #手机