ISCE+Mintpy进行SBASInSAR处理(持续更新)

ISCE安装教程 q由于近期安装isce,实在特别坎坷,对电脑需要足够熟悉,需要改很多参数,才能成功安装成功。 主要安装有两种方法 1.根据isce官网在github上的链接提示https://github.com/scottyhq/isce_notes/tree/masterhttps://github.com/scottyhq/isce_notes/tree/master 这里面有mac苹果电脑安装教程,以及Ubuntu安装的教程,图片如下: 大家如果是18版本的乌班图系统可以参考,施贺青师兄写的i. 阅读详情

目录

ISCE

stackSentinel.py参数列表

ISCE处理流程

stackSentinel.py函数实现

1.边界框处理

2.数据准备与筛选

3.相邻日期对的选择

4.SLC 堆栈处理

5.相关性堆栈处理(correlationStack):

6.干涉图堆栈处理(interferogramStack):

Mintpy

SBAS流程实现

1.配准

2.干涉图生成

3.干涉图滤波

4.相位解缠

5.地理编码

ICA分解


ISCE

SBAS处理流程可以参考博主

【时序InSAR开源软件】用isce2(stackSentinel.py)+mintpy进行哨兵一号SBAS小基线集InSAR时序处理全流程-CSDN博客

【INSAR形变监测】SENTINEL-1 TOPS INSAR全流程 (ISCE2+MINTPY)_sentienl-1 形变处理-CSDN博客

ISCE主要是使用stackSentinel.py生成干涉图

源码:isce2/contrib/stack/topsStack/stackSentinel.py at main · isce-framework/isce2 · GitHub

stackSentinel.py输入输出为:

输入数据:SLC、DEM、Orbits数据

输出数据:一系列可运行文件,每个运行文件实际上有几个相互独立的命令,每个运行文件的配置文件包括执行特定命令/函数的处理选项。

stackSentinel.py

 stackSentinel.py使用示例

# interferogram workflow with 2 nearest neighbor connections (default coregistration is NESD):
stackSentinel.py -s../SLC/ -d ../../MexicoCity/demLat_N18_N20_Lon_W100_W097.dem.wgs84 -b '19 20 -99.5 -98.5' -a ../../AuxDir/ -o ../../Orbits -c 2

# interferogram workflow with all possible interferograms and coregistration with only geometry:
stackSentinel.py -s../SLC/ -d ../../MexicoCity/demLat_N18_N20_Lon_W100_W097.dem.wgs84 -b '19 20 -99.5 -98.5' -a ../../AuxDir/ -o ../../Orbits -C geometry -c all

# correlation workflow with all possible correlation pairs and coregistration with geometry:
stackSentinel.py -s../SLC/ -d ../../MexicoCity/demLat_N18_N20_Lon_W100_W097.dem.wgs84 -b '19 20 -99.5 -98.5' -a ../../AuxDir/ -o ../../Orbits -C geometry -c all -W correlation

# slc workflow that produces a coregistered stack of SLCs
stackSentinel.py -s../SLC/ -d ../../MexicoCity/demLat_N18_N20_Lon_W100_W097.dem.wgs84 -b '19 20 -99.5 -98.5' -a ../../AuxDir/ -o ../../Orbits -C NESD  -W slc

stackSentinel.py参数列表

| 参数 | 短选项 | 长选项 | 类型 | 是否必需 | 默认值 | 选项 | 含义 |
|--|--|--|--|--|--|--|
| display_help | -H | --hh | 无 | 否 | 无 | 无 | 显示详细帮助信息,调用 customArgparseAction 打印 helpstr 并退出程序 |
| slc_directory | -s | --slc_directory | str | 是 | 无 | 无 | 包含所有 Sentinel SLC 的目录 |
| orbit_directory | -o | --orbit_directory | str | 是 | 无 | 无 | 包含所有轨道数据的目录 |
| aux_directory | -a | --aux_directory | str | 是 | 无 | 无 | 包含所有辅助文件的目录 |
| working_directory | -w | --working_directory | str | 否 | ./ | 无 | 工作目录 |
| dem | -d | --dem | str | 是 | 无 | 无 | DEM 文件的路径 |
| polarization | -p | --polarization | str | 否 | vv | 无 | SAR 数据的极化方式 |
| workflow | -W | --workflow | str | 否 | interferogram | slc, correlation, interferogram, offset | InSAR 处理的工作流 |
| swath_num | -n | --swath_num | str | 否 | 1 2 3 | 无 | 要处理的条带列表 |
| bbox | -b | --bbox | str | 否 | 无 | 无 | 经纬度边界(SNWE),示例:'19 20 -99.5 -98.5' |
| exclude_dates | -x | --exclude_dates | str | 否 | 无 | 无 | 要排除处理的日期列表,示例:'20141007,20141031' |
| include_dates | -i | --include_dates | str | 否 | 无 | 无 | 要包含处理的日期列表,示例:'20141007,20141031' |
| start_date | --start_date | 无 | str | 否 | 无 | 无 | 堆栈处理的开始日期,格式应为 YYYY-MM-DD,如 2015-01-23 |
| stop_date | --stop_date | 无 | str | 否 | 无 | 无 | 堆栈处理的结束日期,格式应为 YYYY-MM-DD,如 2017-02-26 |
| coregistration | -C | --coregistration | str | 否 | NESD | geometry, NESD | 配准选项 |
| reference_date | -m | --reference_date | str | 否 | 无 | 无 | 参考采集的目录 |
| snr_misreg_threshold | --snr_misreg_threshold | 无 | str | 否 | 10 | 无 | 使用互相关估计距离失配准的 SNR 阈值 |
| esd_coherence_threshold | -e | --esd_coherence_threshold | str | 否 | 0.85 | 无 | 使用增强谱分集估计方位失配准的相干阈值 |
| num_overlap_connections | -O | --num_overlap_connections | str | 否 | 3 | 无 | 用于 NESD 计算(方位偏移失配准)的每个日期和后续日期之间重叠干涉图的数量 |
| num_connections | -c | --num_connections | str | 否 | 1 | 无 | 每个日期和后续日期之间干涉图的数量 |
| azimuth_looks | -z | --azimuth_looks | str | 否 | 3 | 无 | 干涉图多视处理时的方位视数 |
| range_looks | -r | --range_looks | str | 否 | 9 | 无 | 干涉图多视处理时的距离视数 |
| filter_strength | -f | --filter_strength | str | 否 | 0.5 | 无 | 干涉图滤波的强度 |
| unw_method | -u | --unw_method | str | 否 | snaphu | icu, snaphu | 解缠方法 |
| rmFilter | --rmFilter | 无 | bool | 否 | False | 无 | 是否生成一个额外的去除滤波效果的解缠文件 |
| param_ion | --param_ion | 无 | str | 否 | 无 | 无 | 电离层估计参数文件,若提供将进行电离层估计 |
| num_connections_ion | --num_connections_ion | 无 | str | 否 | 3 | 无 | 用于电离层估计的每个日期和后续日期之间干涉图的数量 |
| useGPU | -useGPU | --useGPU | bool | 否 | False | 无 | 当可用时允许程序使用 GPU |
| num_process | --num_proc | --num_process | int | 否 | 1 | 无 | 每个运行文件中并行运行的任务数量 |
| num_process4topo | --num_proc4topo | --num_process4topo | int | 否 | 1 | 无 | 并行进程的数量(仅针对地形) |
| text_cmd | -t | --text_cmd | str | 否 | 无 | 无 | 要添加到每个运行文件每行开头的文本命令 |
| virtual_merge | -V | --virtual_merge | str | 否 | 无 | True, False | 对合并的 SLC 和几何文件使用虚拟文件,不同工作流有不同默认值 |

其中比较重要的

-s:slc(single look complex)影像所在目录。

-d表示包含影像范围的DEM文件(WGS84),可利用dem.py命令线上下载。

-a为S1-A的的辅助文件目录,因为S1-A的SAR影像可能由不同版本的设备处理仪器(IPF)获得,而S1-B则不需要

-o为S1轨道文件的目录,缺失的文件将自动下载

-b为研究区经纬度范围

-c为SAR影像连接数,确定每景影像与其相邻日期最近的影像生成干涉图的数目,取决于时间基线,对S1而言一般取3

-u为解缠算法,一般选snaphu

-n为SAR影像的subswath的编号,从左到右为1,2,3

-f为干涉影像的滤波强度

-z为方位向的视数

-r为距离向视数,S1的SAR影像的分辨率为5 m * 20 m(距离*方位),视数越小则干涉影像目标点数(或者像元)越多,处理越费时,但是形变细节保留更多。

--reference_date为影像的参考日期。

       执行stackSentinel.py命令后会在当前目录生成一系列文件和文件夹,在runfile文件夹下有16个python脚本文件,分别代表干涉处理的步骤,依次运行脚本则得到一堆干涉影像。


注意:使用stackSentinel.py后可能会生成11个run文件,这个是由于stackSentinel.py参数不同导致的


ISCE处理流程

run_01_unpack_topo_reference:

用于配置和执行解包堆栈参考SLC的操作。它会调用unpackStackReferenceSLC方法,从safe_dict中获取相关信息,对参考Sentinel - 1 SLC数据进行解包,将数据从原始的SAFE文件格式转换为适合后续处理的格式。

run_02_unpack_secondary_slc:

负责配置和执行解包辅助SLC的任务。通过调用unpackSecondarysSLC方法,使用stackReferenceDate、secondaryDates和safe_dict中的信息,对辅助日期对应的Sentinel - 1 SLC数据进行解包处理,同样将数据转换为可进一步处理的格式。

 run_03_average_baseline:

用于配置和执行计算平均基线的操作。它调用averageBaseline方法,利用stackReferenceDate和secondaryDates相关数据,计算平均基线。基线信息在后续的干涉测量等处理中非常重要,它反映了不同SLC数据采集时卫星轨道的相对位置关系。

 run_04_extract_burst_overlaps:

由于inps.coregistration为NESD且updateStack为False,被生成并用于配置和执行提取突发重叠区域的操作。通过调用extractOverlaps方法,从解包后的SLC数据中确定不同突发之间的重叠区域,这些重叠区域的数据将用于后续的精确配准步骤。

 run_05_overlap_geo2rdr:

负责配置和执行将重叠区域从地理坐标转换为距离向坐标并计算偏移量的操作。调用geo2rdr_offset方法,针对secondaryDates对应的SLC数据,在之前提取的重叠区域基础上,进行坐标转换和偏移量计算,为后续的重采样做准备。

  run_06_overlap_resample:

用于配置和执行对重叠区域进行带载波重采样的操作。通过调用resample_with_carrier方法,以secondaryDates数据为基础,结合之前计算的偏移量,对重叠区域的数据进行重采样,使不同SLC数据在重叠区域的坐标和分辨率达到一致,便于后续处理。

 run_07_pairs_misreg:

如果updateStack为True,配置和执行针对secondaryDates的成对配准误差处理;若updateStack为False,则针对acquisitionDates。通过调用pairs_misregistration方法,利用safe_dict中的数据信息,对每对SLC数据进行配准误差的计算和分析,以确定数据之间的精确差异,为后续的时间序列处理提供基础。 

run_08_timeseries_misreg:

用于配置和执行时间序列的配准误差处理。调用timeseries_misregistration方法,综合之前计算的成对配准误差信息,对整个时间序列的SLC数据进行配准误差的进一步分析和处理,以提高数据在时间维度上的一致性和准确性。

 run_09_fullBurst_geo2rdr:

负责配置和执行对全突发数据进行从地理坐标到距离向坐标的偏移处理。调用geo2rdr_offset方法,并设置fullBurst='True',对secondaryDates对应的全突发SLC数据进行坐标转换和偏移量计算,与之前重叠区域的处理类似,但针对的是整个突发数据。

 run_10_fullBurst_resample:

用于配置和执行对全突发数据进行带载波的重采样操作。调用resample_with_carrier方法,同样设置fullBurst='True',对经过地理到距离向坐标偏移处理后的全突发secondaryDates数据进行重采样,确保全突发数据在坐标和分辨率上的一致性,满足后续处理要求。 

 run_11_extract_stack_valid_region:

负责配置和执行提取堆栈有效区域的操作。通过调用extractStackValidRegion方法,对经过前面一系列处理后的SLC堆栈数据进行分析,确定数据中的有效区域,排除无效或噪声数据区域,提高后续处理的数据质量。 

run_12_merge_reference_secondary_slc:

因为mergeSLC为True,用于配置和执行合并参考SLC和辅助SLC的操作。通过调用mergeReference方法合并参考日期的SLC数据,调用mergeSecondarySLC方法合并secondaryDates对应的辅助SLC数据,并设置virtual = 'False',表示采用实际合并而非虚拟合并的方式,将多个SLC数据合并为一个更便于后续处理的数据集。 

run_13_grid_baseline:

负责配置和执行网格化基线的操作。调用gridBaseline方法,使用stackReferenceDate和secondaryDates相关数据,将之前计算的平均基线信息进行网格化处理,将基线数据转换为适合网格化分析和后续应用的格式。

run_14_merge_burst_igram:

配置了合并突发干涉图的操作,接着调用runObj.igram_mergeBurst(acquisitionDates, safe_dict, pairs)方法,将之前生成的突发干涉图进行合并,以得到更完整或适合后续处理的干涉图数据。

run_15_filter_coherence:

用于配置对干涉图相干性进行滤波的操作,然后调用runObj.filter_coherence(pairs)方法,对由pairs确定的干涉图数据进行相干性滤波处理,以提高干涉图的质量和可靠性。

InSAR个人笔记(一)ISCE安装(利用anaconda) 一般选择Anaconda3-2023.03-Linux-x86_64.sh ,下载。安装完成,输入conda,检查conda是否添加到了环境变量中,若显示。自动加入环境变量,之后输入conda,显示。Anaconda3在ubuntu22.04中安装和使用(自用)_赈早见 琥珀川144的博客-CSDN博客。Ubuntu安装步骤(图文)_安装ubuntu_Anliy876的博客-CSDN博客。Ubuntu安装图文记录_安装ubuntu。3、通过conda安装isce2。1\ubuntu系统的安装。 阅读详情

相关推荐

InSAR学习笔记之ISCE 软件安装

InSAR学习笔记之ISCE 软件安装 ISCE是一款常用的InSAR数据处理软件,2018年更新的版本基于ubuntu18.04系统安装过程简化了很多,本文分享一下安装过程。(之前在ubuntu16上安装过,过程非常坎坷,花了很长时间才成功) 配置安装环境 sudo apt update sudo apt install -y gfortran libmotif-dev libhdf5-dev libfftw3-dev libgdal-dev scons python3 cython3 python3-s

weixin_42128764的博客 8574

ISCE开源软件Linux系统配置-小白从头到尾篇

虚拟机+Ubuntu系统安装 安装Anaconda 下载清华大学Anaconda镜像 ISCE安装 ISCEgithub地址 根据说明文件,首先需要检查基本软件状况。 2.1 安装github sudo apt install git 2.2 clone源码 git clone https://github.com/isce-framework/isce2.git 2.1...

qq_45110581的博客 2952

Isce+mintpy处理流程

配准+干涉对形成+解缠,采用stackSentinel.py命令的示例1或示例2.示例1:stackSentinel.py -s */mintpy_test/slc/ -d */mintpy_test/DEM/demLat_S01_N00_Lon_W092_W090.dem.wgs84 -b '-0.60 -0.20 -91.75 -91.35' -a */mintpy_test/orbits/ -o */mintpy_test/orbits/ -c 2 -m 20180619。

weixin_56301271的博客 2821

ISCE软件常用命令

1.下载和拼接dem mkdir DEM; cd DEM dem.py -a stitch -b 18 20 -100 -97 -r -s 1 -c rm demLat*.dem demLat*.dem.xml demLat*.dem.vrt cd .. 2.stamps预处理 stackSentinel.py -s ../SLC/ -d ../DEM/demLat_N18_N20_Lon_W100_W097.dem.wgs84 -a ../../AuxDir/ -o ../../Orbits

m0_60086135的博客 1901

isce3:InSAR科学计算环境***注意-该存储库当前是内部托管项目的镜像。 暂时未接受外部拉动请求。***

ISCE-InSAR科学计算环境 | InSAR科学计算环境(ISCE)是一个开放源代码库,用于处理星载和机载干涉式合成Kong径雷达(InSAR)数据。 该项目是框架的继承者。 这是一个彻底的重新设计,着重于改进的模块化,文档和测试驱动的开发。 它的最初开发是由美国航空航天局地球科学技术办公室(ESTO)根据高级信息系统技术(AIST)2008资助的,目前由美国国家航空航天局ISRO SAR(NISAR)项目资助。 :warning: 注意:ISCE3处于早期开发中-其功能和界面可能会发生变化 安装 入门 | 贡献 欢迎捐款。 有关入门的提示和贡献的准则,请参阅贡献的指南。 提交问题,请让我们知道您是否遇到错误。 执照 这是“按原样”提供给您的研究代码,无任何准确性保证。 自行承担风险使用。 此软件的使用由许可文件中指定的条款控制。 此软件的出口受EAR99 NLR分类条款中规定的美

别再让地图加载慢吞吞!GeoServer + GeoWebCache 缓存配置实战(附性能对比)

本文详细介绍了如何通过GeoServer和GeoWebCache优化地图服务性能,实现地图加载速度的显著提升。从缓存机制解析到实战配置流程,再到高级性能调优技巧,帮助开发者构建高效地图缓存系统,解决地图加载慢的问题。

weixin_42525369的博客 400

android 运动目标检测_MindSpore应用案例:AI对篮球运动员目标的检测

MindSpore作为一个端边云协同的开源的全场景AI框架,今年3月份开源以来,受到了业界的广泛关注和应用,欢迎大家参与开源贡献、模型众智合作、行业创新与应用、学术合作等,贡献您在云侧、端侧 (HiMindSpore 应用程序)、边侧以及安全领域的应用案例。MindSpore与您在AI领域共同成长,让AI使能千行百业,释放出强大的生产力。本文讲述的是MindSpore对篮球运动员目标的检测应用,通...

weixin_33603813的博客 824

three.js开发3D地图记录(二)

添加地图点击事件,点击区域,地图区域改变颜色并浮起来。

anwebDesigner的博客 1687

【智能体】零代码实战:Coze多智能体协作开发全流程解析

本文详细解析了Coze平台的多智能体协作开发全流程,从入门指南到实战技巧,帮助用户零代码构建高效AI团队。通过模块化设计和专业分工,多智能体系统能显著提升任务处理速度和准确率,适用于客服、内容创作、教育培训等多个场景。文章还提供了性能优化和常见问题解决方案,助力开发者快速掌握Multi-Agent技术。

weixin_29239279的博客 321

intel-VROC raid卡驱动 centos7.9

u盘进入系统,按e,在出现下面在quiet前面输入inst.dd initcall_blacklist=vmd_drv_init 后按Ctrl和x,开始引导。进入系统,按e,在出现下面在quiet前面输入inst.dd inicll_blacklist=vmd_drv_init 后按Ctrl和x,开始引导。驱动完成后,在输入 c 回车, 引导过程开始,进入安装界面。输入1 ,选择这个ISO驱动文件里的rpm包。输入1 ,选择这个ISO驱动文件里的rpm包。vi编辑 /etc/default/grub。

dacming的博客 4532

保姆级教程:手把手教你从NASA Earthdata下载ASTER L2地表温度数据(附shp文件拖拽技巧)

本文提供了一份详细的NASA Earthdata平台ASTER L2地表温度数据下载教程,涵盖账号注册、数据筛选、shp文件拖拽技巧及质量检查等关键步骤。特别适合遥感领域初学者,帮助快速掌握高效获取ASTER L2T数据的方法,适用于城市热岛效应和农业监测等研究。

weixin_30736301的博客 465

ffmpeg 4.4 cenc-aes-ctr 加解密 MP4 工程性质分析

ffmpeg 4.4 cenc-aes-ctr 加解密 MP4 工程性质分析

Dream_Dynasty的博客 3026

Windows Server 2019 服务器一天之中内存利用率越来越高直到死机,疑似内存泄漏

环境: 华为桌面云虚拟机 Windows Server 2019 版本1809 OS内部版本17763.107 问题描述: Windows Server 2019服务器 一天之中内存利用率越来越高直到死机,不断扩容内存从10G到16G一重启就正常然后慢慢又循环,连续几天都这样,疑似内存泄漏 内存泄漏是什么? 内存泄漏指由于疏忽或错误造成程序未能释放已经不再使用的内存。内存泄漏并非指内存在物理上的消失,而是应用程序分配某段内存后,由于设计错误,导致在释放该段内存之前就失去了对该段内存的控制,从而造成了内存的浪

本博客,博文仅代表个人操作经验,不能完全解决你的问题,仅供参考,佛系回复。 5384

复分析 阿尔福斯(Ahlfors) 第三版答案 英文

答案分为十章,有部分习题没有答案,全英文复分析阿尔福斯(Ahlfors) 第三版答案

NTscan扫描工具

ntscan是一个扫描器,可以扫描没有密码和弱密码的共享,先到配置文件中填写猜测的用户名和密码,扫描过程中右侧出现的是扫描日志。

廖雪峰Python教程[完整版]

廖雪峰的python教程整理的,已经排好版了,而且还是完整的,适合小白

C# WinForm 基于DirectShow.dll编写的程序,可边预览边保存视频,音频和图片

基于微软的DirectShow.dll的二次开发1.自动识别连接的摄像头跟麦克风且可以对其进行选择(默认选中第一个)2.可录制wmv和avi两种格式的视频并且都带有音频,avi可以录制的分辨率和帧数进行选择,wmv默认进行了淸晰度优化如不喜欢可去除代码中有相关提示3.可在预览点开之后进行拍照,也可边录制视频边拍照4.添加了视频和音频的压缩功能(酌情使用)这次二次开发发现C#对DirectShow的相关信息太少,这次做完也不敢藏私共享出来,还望大家指正

上一篇: c# 进阶 反射
吃西瓜的芋头
博客等级 码龄5年 38粉丝 22原创
评论 2
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付:10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值