Sarscape5.6版本中导入外部控制点、写入精密轨道文件与GACOS用于大气相位

news2024/11/7 18:05:50

SARscape中导入外部GCP点用于轨道精炼

https://www.cnblogs.com/enviidl/p/16524645.html

在SAR处理时,有时会加入GCP点文件,SAR处理中用到的控制点分为两类:用于校正地理位置的几何控制点(Geometry GCP)和用于轨道精炼的控制点(Refinement GCP)。

轨道精炼控制点在以下处理步骤使用:

InSAR/DInSAR处理工作流:/SARscape/Interferometry/Phase Processing/4 - Refinement and Re-flattening
立体SAR处理工作流:/SARscape/Interferometry/Stereo-Radargrammetry/2 - Shift Refinement and Re-flattening
MAI处理工作流:/SARscape/Interferometry/MAI Processing/2 - MAI Refinement and Re-flattening
SBAS处理工作流:/SARscape/Interferometric Stacking/SBAS/3 - Refinement and Re-Flattening
轨道精炼GCP点的选择常规方法是在轨道精炼的步骤,输入文件和DEM以及参考文件,手动在输入文件上选择GCP点。

有些情况下,用户希望能从已有的底图上(光学底图或GoogleEarth上)选择GCP,或者通过已有点坐标来生成GCP文件。比如做地表形变测量时,用户基于对工程区已有的经验,在光学底图上能确定稳定的点的位置,那么在光学底图上找到的GCP点如何导入SARscape中用于轨道精炼呢?本文以SBAS工作流,从Google Earth上选择GCP为例,介绍从已有底图上选的GCP点用于轨道精炼方法。

以下的操作是在SBAS工作流已经完成了连接图生成、干涉工作流之后进行:

第一步

,打开工程区的kml文件,在Google Earth上,使用工具,在工程区内定位若干个GCP点(点要位于没有形变的位置、最好是人工地物上,相干性高),将位置另存为kmz文件。
在这里插入图片描述

第二步

:在ArcGIS中,将kmz文件转为shapefile文件,工具是kmz to layer,再将layer另存为shapefile文件。

注:此方法仅供参考,把kmz转成shapefile即可。
在这里插入图片描述
注:以上两步是在Google Earth上选GCP点的操作,如果用户要在地理坐标的光学底图上选GCP点,直接在ENVI或ArcGIS中生成shape矢量点图层即可。如果有已知的GCP点的xy坐标,直接用点坐标文件生成shape文件即可。

第三步:

将上一步得到的地理坐标系的shp,在SARscape转为SAR坐标系的shapefile文件。打开/SARscape/Basic/Intensity Processing/Geocoding/Map to SAR Shape Conversion工具。

Input Files:

  • Input File:选择上一步准备好的地理坐标系的shapefile文件;

  • Input Reference File:选择SAR坐标系的参考数据,(在接下来的轨道精炼处理中,GCP点要在哪个图像上使用,就选哪个),此处选择干涉工作流中生成的某一对像对的fint文件作为参考数据;

  • 在这里插入图片描述
    Optional Files:输入用于几何校正的GCP文件,此处不输入。
    DEM/Cartographic System:输入参考DEM或者选择坐标系,作为地理坐标系的参考,此处输入参考DEM文件。
    在这里插入图片描述
    在这里插入图片描述

第四步:

在轨道精炼一步,选择GCP的界面上,导入上一步生成的斜距shapefile文件,即可进行轨道精炼。

在轨道精炼数据输入面板上,输入工程文件,点击Create GCP,选择input file和参考DEM文件,点击Next,
在这里插入图片描述
图 点击Create GCP

在GCP生成面板上上,点击Load GCPs,选择上一步生成的SAR坐标系的shapefile文件,导入GCP点,可检查点在参考图像上的位置,如果有不在图像范围的GCP点,删除即可。导入的GCP点号有乱码的情况,忽略即可。也可以手动再添加新的GCP点。
在这里插入图片描述

图 点击Load GCPs导入SAR坐标系的GCP点shapefile文件
在这里插入图片描述

注:有的版本不支持shp中的中文属性内容,可将shp中的相应中文字段删掉即可,再使用,如下图:

图 编辑shp点的属性表
切换到Cartographic System,选择地理坐标系:GEO-GLOBAL-GEO-WGS84
在这里插入图片描述
切换到Export面板,设置输出xml点文件的路径,点击Finish。
在这里插入图片描述
在面板上确认输入生成的控制点文件之后,进行轨道精炼操作即可。
在这里插入图片描述

SARscape导入哨兵1数据自动读取轨道文件

https://www.cnblogs.com/enviidl/p/16333401.html

在SARscape中,哨兵数据导入之前,将数据的轨道文件放置于系统指定的目录下,软件在导入数据时将自动读取相应数据的轨道文件,不需要手动依次输入。

本文在SARscape5.6.2版本中操作,其他版本操作相同。

操作方法如下:

第一步,

打开系统参数/SARscape/Preferences/Preferences common,在系统参数面板中,Directories and batch file name选项卡中,Sentinel-1 auxiliary directory指定一个用来存放哨兵数据轨道文件的路径,注意文件夹均为英文命名。点击OK。

在这里插入图片描述
图 在系统参数中指定哨兵数据轨道文件的存放路径

第二步,

在该路径下新建文件夹,命名为:AUX_POEORB,将下载好的哨兵数据精密轨道文件直接放在该路径下。该文件夹中可以放所有时间的精密轨道文件,程序会自动挑选。

注:1、如果要使用到回归轨道文件,可以建一个名为AUX_RESORB的文件夹,将回归轨道文件放于该文件夹中,如果不使用回归轨道文件,可以不建此文件夹;

2、程序会首先在AUX_POEORB文件下中查找精密轨道文件,如有,就直接使用,如没有,再会在AUX_RESORB文件夹中查找回归轨道文件,如果两个文件夹中都没有查找到相应轨道文件,程序会使用哨兵1数据元数据自带的轨道信息进行导入;

3、如果数据要进行干涉处理,导入时建议使用精密轨道文件。

在这里插入图片描述
图 将下载好的精密轨道文件放在系统参数指定的位置

第三步,

打开菜单/SARscape/Import Data/SAR Spaceborne/Single Sensor/SENTINEL-1,在数据Input File List面板上,输入需导入的哨兵1数据,可输入一景,也可以同时输入多景(5.6.2版本直接输入压缩包即可)。

在这里插入图片描述
图 哨兵数据导入界面

执行哨兵1数据导入之后,如果没有找到对应数据的精轨文件或回归轨道文件,会有相应的提示,如下图。

在这里插入图片描述
图 未找到轨道数据的提示

注:如果使用的是5.2.1版本,输入数据面板中,下面输入轨道文件的地方空着即可,如下图:

在这里插入图片描述
图 5.2.1版本的哨兵数据输入界面

SARScape使用GACOS数据

https://blog.csdn.net/qq_41159191/article/details/127453080

1 GACOS数据下载

GACOS官网

http://www.gacos.net/

在这里插入图片描述
Time of insterest(in UTC)时间设置
这个在数据名称里面可以看到,如果数据的条带一致,那么每日拍摄时间应该也是一样的,如下图红框所示,T后面的数字就是HHMMSS,即10时18分43秒
在这里插入图片描述
选择二进制文件Binary grid
在这里插入图片描述
提交之后,如果没有出错,会弹出这个界面
在这里插入图片描述

2 GACOS数据处理

2.1 查看数据

下载好的GACOS数据为tar.gz格式
在这里插入图片描述
解压数据并查看Readme
GACOS tropospheric delay maps are given in a grid binary format (4-byte float little endian, naming
convention YYYYMMDD.ztd).

2.2 Import GACOS

工具路径/SARscape/Import Data/Other Format/GACOS,打开后会弹出help界面
在这里插入图片描述
弹出提醒,点击是,然后就可以看到Import Gacos了
在这里插入图片描述
输入GACOS数据,记得先解压,因为导入时候只能选择.ztd和.ztd.tif格式的数据,而下载好的数据是tar.gz格式的压缩包文件
在这里插入图片描述
在这里插入图片描述
快视图和导出设置
在这里插入图片描述
在这里插入图片描述
这一步很快,立即完成
在这里插入图片描述

2.3 GACOS用于大气相位延迟校正

工具地址:/SARscape/Interferometry/Interferometric Tools/Atmospheric Phase Delay Correction

输入干涉图fint和斜距dem(srdem)
在这里插入图片描述
输入主master从slave影像的GACOS数据,一定要看准日期
在这里插入图片描述
在Parameters中的Atmosphere External Sensors(大气外部传感器)中选择GACOS
在这里插入图片描述
导出数据
在这里插入图片描述

2.4 GACOS用于高程校正

高程校正,只需要在optional files中选择解缠文件
在这里插入图片描述
在parameters中勾选atmosphere height correlation flag(大气高度相关标志)为True
通过选择此标志,算法会估算并从每个干涉图中删除大气的高度相关成分,使用“高度相关窗口大小 [m]”(大气部分)参数来定义过滤器的大小。
在这里插入图片描述

3 GACOS在SBAS-InSAR中的使用

工具:/SARscape/Interferometric Stacking/SBAS/2 - Interferometric Process

使用方法:
将之前import好的GACOS数据导入对应的列表,在SBAS Interferometric Process的位置为Optional Water Vapour File List
在这里插入图片描述
如果要使用GACOS数据,请务必添加以下引用:

Yu, C., Li, Z., Penna, N. T., & Crippa, P. (2018). Generic atmospheric correction model for Interferometric Synthetic Aperture Radar observations. Journal of Geophysical Research: Solid Earth, 123(10), 9202-9222.

Yu, C., Li, Z., & Penna, N. T. (2018). Interferometric synthetic aperture radar atmospheric correction using a GPS-based iterative tropospheric decomposition model. Remote Sensing of Environment, 204, 109-121.

Yu, C., Penna, N. T., & Li, Z. (2017). Generation of real‐time mode high‐resolution water vapor fields from GPS observations. Journal of Geophysical Research: Atmospheres, 122(3), 2008-2025.

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1114797.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

C++类和对象(三) (this指针)

this指针 1 this指针的引出 我们先来定义一个日期类 Date class Date { public:void Init(int year, int month, int day){_year year;_month month;_day day;}void Print(){cout << _year << "-" << _month << "-" << …

【电子通识】USB接口三大类型图册

基本概念 不同时期的USB接口有不同的类型&#xff0c;USB接口分为插头和插座&#xff1a; 插头&#xff0c;plug&#xff0c;对应的也叫公口&#xff0c;即插别人的。 插座&#xff0c;receptacle&#xff0c;对应也叫做母口&#xff0c;即被插的。 USB的接口类型&#xff0…

《动手学深度学习 Pytorch版》 9.6 编码器-解码器架构

为了处理这种长度可变的输入和输出&#xff0c; 可以设计一个包含两个主要组件的编码器-解码器&#xff08;encoder-decoder&#xff09;架构&#xff1a; 编码器&#xff08;encoder&#xff09;&#xff1a;它接受一个长度可变的序列作为输入&#xff0c;并将其转换为具有固定…

模拟器-雷电-使用adb push或adb pull操作文件

一、环境 windows 10 雷电模拟器4.0.83 二、问题 有时候我们会需要往模拟器拷贝文件或者复制文件到我的电脑 三、方法 1、获取root权限 adb root adb remount 有可能遇到【daemon not running; starting now at tcp:5037】的报错 查看端口占用进程&#xff1a;netstat -…

vue2升级到vue2.7

vue2升级到vue2.7 小小的改进,大大的提升 只需要简单修改,开发体验得到大大提升. 为什么要升级Vue2.7 不能拒绝的理由: 组合式 API(解决mixins问题:命名冲突,隐式依赖)单文件组件内的 <script setup>语法模板表达式中支持 ESNext 语法(可选链:?.、空值合并:??)单文…

域渗透05-协议(NTLM)

前言&#xff1a; 首先在域内环境中主要有如下6种身份认证协议&#xff0c;最常用的还是Kerberos&#xff0c;NTLM和LDAP三种&#xff1a; Kerberos&#xff1a;Kerberos是一种网络身份验证协议&#xff0c;用于验证用户和服务之间的身份。它通过使用票据和密钥来实现认证&…

Linux 指令学习

Linux 指令学习 以此为记录&#xff0c;也方便自己日后查看回顾&#xff01; Linux命令基础格式 无论是什么命令&#xff0c;用于什么用途&#xff0c;在Linux中&#xff0c;命令有其通用的格式&#xff1a; command&#xff1a; 命令本身 options&#xff1a;[可选&#xf…

AD9371 官方例程

文章目录 前言一、HDL方面1. ZYNQ 核根据ZCU106平台修改&#xff08;**参考UG1244 ZCU106 Evaluation Board**&#xff09;&#xff0c;尤其注意**DDR**的配置&#xff08;**参考美光 MTA4ATF51264HZ**&#xff09;2.dacfifo 深度 要修改 &#xff0c;地址位宽 小于等于16&…

Windows 11 2023 Moment 4 更新内容详解

Windows 11 2023 Moment 4 更新内容详解 编辑&#xff1a;简简单单 Online zuozuo 地址&#xff1a;https://blog.csdn.net/qq_15071263 前言 主要的更新内容清单 Copilot内置应用资源管理器照片应用截图工具画图与视频剪辑Windows 备份系统界面开始菜单快速设置系统托盘任务…

【python】进程和线程

进程和线程 今天我们使用的计算机早已进入多CPU或多核时代&#xff0c;而我们使用的操作系统都是支持“多任务”的操作系统&#xff0c;这使得我们可以同时运行多个程序&#xff0c;也可以将一个程序分解为若干个相对独立的子任务&#xff0c;让多个子任务并发的执行&#xff…

优测云测试平台 | 有效的单元测试

一、前言 本文作者提出了一种评价单元测试用例的质量的思路&#xff0c;即判断用例是否达到测试的“四大目标”。掌握识别好的用例的能力&#xff0c;可以帮助我们高效地写出高质量的测试用例。 评判冰箱的好坏&#xff0c;并不需要有制造一台冰箱的能力。在开始写测试用例之…

TikTok Shop印尼站关停,跨境卖家该如何自救

10.4日&#xff0c;TikTok Shop宣布关停印尼站&#xff0c;无疑是一个巨大的炸弹投在整个跨境圈&#xff0c;让许多跨境卖家以及当地本土的卖家、品牌措手不及。 这对TikTok的东南亚市场而言&#xff0c;无疑是一次沉重的打击。 TikTok在东南亚的商业布局&#xff0c;印尼是其…

小A对我说,他现在快想钱想疯了…

昨天晚上11点&#xff0c;小A给我打电话 小A问&#xff1a;橙哥&#xff0c;有没有赚钱的事做&#xff1f; 他说&#xff0c;实在不想上班了&#xff0c;每天起早贪黑挤地铁 辛辛苦苦干一个月&#xff0c;到手工资三四千块钱&#xff0c;房租一交&#xff0c;日常开支一花&a…

Jmeter安装(快速入门)

大家好我是苏麟今天简单聊一下Jmeter . Jmeter Jmeter官网 : Apache JMeter - Download Apache JMeter 安装Jmeter Jmeter依赖于JDK&#xff0c;所以必须确保当前计算机上已经安装了JDK&#xff0c;并且配置了环境变量。 下载 : 官网下载 : 网盘自取 : 链接&#xff1a;…

【LeetCode】57. 插入区间

1 问题 给你一个 无重叠的 &#xff0c;按照区间起始端点排序的区间列表。 在列表中插入一个新的区间&#xff0c;你需要确保列表中的区间仍然有序且不重叠&#xff08;如果有必要的话&#xff0c;可以合并区间&#xff09;。 示例 1&#xff1a; 输入&#xff1a;interval…

国产开发板上打造开源ThingsBoard工业网关--基于米尔芯驰MYD-JD9X开发板

本篇测评由面包板论坛的优秀测评者“JerryZhen”提供。 本文将介绍基于米尔电子MYD-JD9X开发板打造成开源的Thingsboard网关。 Thingsboard网关是一个开源的软件网关&#xff0c;采用python作为开发语言&#xff0c;可以部署在任何支持 python 运行环境的主机上&#xff0c;灵…

前端将二维数组转换成树形格式数组

原始数据&#xff1a; 最终得到的数据&#xff1a; 具体实现代码&#xff1a;倒推&#xff08;最后一个数据&#xff09;的思路得到数组 function childData(list,i){const arr []//最后一个数组if(i<list.length){list[i].forEach(item>{const obj {name:,prop:,chil…

美摄视频特效SDK,视频特效技术解决方案

视频内容已经成为企业传播信息、提升品牌形象的重要手段。然而&#xff0c;普通的视频内容往往难以吸引观众的注意力&#xff0c;而富有创意的视频特效则能够大大提升视频的观赏性和吸引力。为了帮助企业轻松实现视频特效的添加&#xff0c;我们推出了美摄视频特效SDK&#xff…

风口下的低代码,怎么样了?

目录 一、前言 二、为什么低代码平台会受到资本的重视&#xff1f; 三、低代码的应用场景 1.基于表单/引擎驱动的模式 2.基于aPaaS平台的模式 四、低代码能解决什么问题&#xff1f; &#xff08;1&#xff09;运维工作流方向 &#xff08;2&#xff09;运维开发方向 五、低代码…

ArcGIS笔记10_如何创建渔网?

本文目录 前言Step 1 确定渔网的精度单位Step 2 有底图时创建渔网的操作 前言 ArcGIS中的渔网是一个很好用的工具&#xff0c;它可以创建出规规整整的小格子&#xff0c;每个小格子都对应一个标注点&#xff0c;可以将原本散乱的数据规整化&#xff0c;如下图&#xff1a; Ste…