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

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.mzph.cn/news/112432.shtml

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

相关文章

多测师肖sir_高级金牌讲师___ui自动化之selenium001

一、认识selenium (1)selenium是什么? a、selenium是python中的一个第三方库 b、Selenium是一个应用于web应用程序的测试工具,支持多平台,多浏览器,多语言去实现ui自动化测试,我们现在讲的Sel…

Atlassian Confluence OGNL表达式注入RCE CVE-2021-26084

影响版本 All 4.x.x versions All 5.x.x versions All 6.0.x versions All 6.1.x versions All 6.2.x versions All 6.3.x versions All 6.4.x versions All 6.5.x versions All 6.6.x versions All 6.7.x versions All 6.8.x versions All 6.9.x versions All 6.1…

RK3568笔记四:基于TensorFlow花卉图像分类部署

若该文为原创文章,转载请注明原文出处。 基于正点原子的ATK-DLRK3568部署测试。 花卉图像分类任务,使用使用 tf.keras.Sequential 模型,简单构建模型,然后转换成 RKNN 模型部署到ATK-DLRK3568板子上。 在 PC 使用 Windows 系统…

使用telegram机器人发送通知

文章目录 背景1 创建机器人2 与机器人的会话3 调用API让机器人发送消息 背景 在训练深度学习模型时,除了粗略估计外,很难预测训练何时结束。此外,我们可能还想随时随地查看训练情况,如果每次都需要登录回服务器的话并不方便。因此…

wordpress网站部署了ssl证书之后就排版混乱了

刚给自己的小网站部署了SSL证书,之后就发现https访问主页竟然乱套了。在手机上访问却是正常的。 直接上解决方案: 编辑网站根目录下的wp-config.php文件 在自定义文本处添加以下代码: if ($_SERVER[HTTP_X_FORWARDED_PROTO] https) $_SE…

PHP-FFMpeg 操作音视频

✨ 目录 🎈 安装PHP-FFMpeg🎈 视频中提取一张图片🎈 视频中提取多张图片🎈 调整视频大小🎈 视频添加水印🎈 生成音频波形🎈 音频转换🎈 给音频添加元数据🎈 拼接多个音视…

利用ArcGIS获取每一个冰川的中心位置经纬度坐标:要素转点和要素折点转点的区别

问题概述:下图是天山地区的冰川的分布,我们可以看到每一条冰川是一个面要素,要求得到每一个冰川(面要素)的中心经纬度坐标。 1.采用要素转点功能 选择工具箱的【数据管理工具】-【要素】-【要素转点】。完成之后再采用…

计算机基础知识36

数据库数据的演变史 ATM:1. 把数据都存在了文件中,文件名不规范 kevin|123 kevin123 kevin*123 2. 存储数据的文件越来越多,放在db文件夹,占用空间,查询存储不方便,速度慢 # 数据库软件能解…

lnmp架构部署Discuz论坛并配置重定向转发

lnmp架构部署Discuz论坛并配置重定向转发 文章目录 lnmp架构部署Discuz论坛并配置重定向转发环境说明部署Discuz论坛系统下载Discuz论坛系统代码包,官网地址如下:部署Discuz论坛系统步骤:解压安装Discuz源码包配置虚拟主机进入Discuz安装界面…

Janus: 逆向思维,以数据为中心的MoE训练范式

文章链接:Janus: A Unified Distributed Training Framework for Sparse Mixture-of-Experts Models 发表会议: ACM SIGCOMM 2023 (计算机网络顶会) 目录 1.背景介绍all-to-allData-centric Paradigm 2.内容摘要关键技术Janus细粒度任务调度拓扑感知优先级策略预取…

Android推送问题排查

针对MobPush智能推送服务在使用过程中可能出现的问题,本文为各位开发者们带来了针对MobPush安卓端推送问题的解决办法。 TCP在线推送排查 排查TCP在线收不到推送时,我们先通过客户端的RegistrationId接口获取设备的唯一标识 示例: MobPush…

代码随想录算法训练营第五十六天 | 1143.最长公共子序列、1035.不相交的线 、53. 最大子序和 动态规划

1143.最长公共子序列 视频讲解:动态规划子序列问题经典题目 | LeetCode:1143.最长公共子序列_哔哩哔哩_bilibili 代码随想录 (1)代码 1035.不相交的线 视频讲解:动态规划之子序列问题,换汤不换药 | Leet…

高数二阶导数例子

例子: 当 x 0 x 0 x0时 f(x) 1,当 x ≠ 0 x \ne 0 x0时, f ( x ) sin ⁡ x x f(x) \frac{\sin x}{x} f(x)xsinx​, 求f’(x)。 解: 二阶导数是导数的导数,按照定义为: f ′ ′ ( x ) lim ⁡ x →…

基于Java的垃圾分类管理系统设计与实现(源码+lw+部署文档+讲解等)

文章目录 前言具体实现截图论文参考详细视频演示为什么选择我自己的网站自己的小程序(小蔡coding) 代码参考数据库参考源码获取 前言 💗博主介绍:✌全网粉丝10W,CSDN特邀作者、博客专家、CSDN新星计划导师、全栈领域优质创作者&am…

安卓逆向之抽象函数public abstract的hook定位处理

一、抽象类调用过程 1、定义接口 public abstract class AbsNetHelper{public abstract int X(int x,int y); }2、其他类调用 public class f extends AbsNetHelper {public abstract int X(int x,int y){ return xy;} }二、如何HOOK 想hook调用AbsNetHelper类的X抽…

linux uboot kenerl filesystem关系

每次烧写都是uboot,image,roots,为什么不整一套,为什么分开 什么是uboot 什么是内核,三者为什么分开 如果你上网去搜一下,你会发现这样的答案 另注:本文含有大量的错误和意测成份&#xff0c…

Linux系统编程04

进程的概念 进程(动态)是一个正在运行的程序(静态) 多道程序设计缺点: (1)缺乏隔离,各个程序之间可以直接访问,使用对方的数据 (2)内存使用率低&a…

香港服务器在国内访问太慢怎么能提高?

​  一直以来,全球化业务需求的增长是跟随着蓬勃向上的互联网而发展的。有了网络,海外贸易就在鼠标的轻点中完成。而IDC市场中的香港服务器也因为免备案政策的特性,开始逐渐成为企业想要跨越地域壁垒而考虑的对象。但在使用过程中&#xff…

新手如何找到Docker容器(redis)中的持久化文件?

具体步骤 要查看Docker容器的dump.rdb和appendonly.aof文件(如果启用了AOF持久化)的位置,我们需要知道容器中Redis配置文件的内容或者容器的数据卷的挂载位置。 这里是一般步骤: 查找容器的数据卷挂载位置 使用docker inspect命令…

内存爆炸、CPU100%问题定位

目录 一、内存爆炸相关1、关于-Xms(最小堆内存)和-Xmx(最大堆内存)2、JVM初始化时申请实际物理内存3、OutOfMemory问题排查(1) 堆内存溢出排查(2) 堆外内存溢出排查 二、CPU 100及死锁问题定位1、CPU 100问题排查(1) 找到程序对应进程号(2) 查找进程对应的线程编号(3)查看线程堆…