ISCE2023学习（0.7）数据获取

博主在这里记录自学ISCE2 相关的workshop的学习情况

课程官网：https://www.earthscope.org/event/2023-insar-isce-short-course/
课程表以及课程内容详细介绍：https://docs.google.com/document/d/116EXYa5UTSz8LOp5oett-0MpVxBlFcfZ4Snb38JcCxA/edit
计划使用课程所提供的notebook及youtube教学视频

---

前言

本博客仅作为个人学习笔记用途，欢迎和各位前辈交流ISCE2学习

本节课主要包括三个文件

• Data_Access_Accounts.ipynb（创建账号）
• GRFN.ipynb （没看明白，给出的NISAR链接已失效）
• SSARA.ipynb （一个集成多方数据的获取方法）

介绍了几个常用的库以及所需要的类似通行证的帐号，包括UNAVCO/WInSAR，ASF（Alaska Satellite Facility ），ESA EO（欧空局Earth Observation ）。

---

一、创建账号

1. 首先在NASA EarthData创建账号，EarthData优先级应该是比较高的，ASF可以直接使用EarthData进行登录。EarthData注册起来比较方便，中间没太费功夫。

   需要注意的是，后边有很多软件，如ISCE2等在运行过程中会查看隐藏文件\.netrc，该文件储存了用户的EarthData帐号信息，因此在注册完成之后，我们在notebook内运行代码创建.netrc文件。

注意不要直接复制运行
将myUsername替换为你的用户名，myPassword替换为密码后运行

!echo "machine urs.earthdata.nasa.gov login myUsername password myPassword" > $HOME/.netrc

---

2. 注册WInSAR帐号，需要提供相关的机构信息并使用机构内邮箱进行注册。

3. 注册欧空局帐号，应该没有机构限制但网速较慢，流程也更繁琐些，可以接受。

4. 注册OpenTopography并获取api
   登陆后点击下方图片中按钮即可获得一个api key
   
   
   同样的，这里需要创建一个.topoapi文件来让系统默认从你的api key中下载数据，注意把myAPIkey替换为自己的api key

!echo "myAPIkey" > $HOME/.topoapi

---

为保护个人隐私安全，课本中还对两个文件进行了chmod权限修改，仅自己可以访问

!chmod 600 $HOME/.topoapi
!chmod 600 $HOME/.netrc

二、SSARA使用

SSARA是一个集合的数据库，能够用于检索ASF、ESA、WInSAR的部分雷达影像数据，目前看来教程中比较推崇这一方式，因为可以满足直接在CLI界面中进行检索以及下载（也提供了GUI界面下载的方法），同时能够找到足够满足学员需求的数据。

SSARA的核心相当于是一个python脚本文件，运行起来就能用。

1.下载SSARA

2023年UNAVCO合并到EarthScope中去，所以部分UNAVCO的旧网站可能已经不再使用或者维护中了，课本中提供的gitlab链接在维护中，下载不了。SSARA gitlab链接在下：
Seamless SAR Archvie (SSARA) Client project

git clone https://gitlab.com/earthscope/public/sar/ssara_client.git $HOME/work/SSARA

下载好之后SSARA就在系统的$HOME/work/SSARA中。修改图片中的文件：

unavuser 对应 WInSAR帐号
asfuser 对应 ASF帐号
eossouser 对应 ESA帐号

直接在notebook中运行查看结果：

!$HOME/work/SSARA/ssara_federated_query.py

不含内容的query会直接给出参数提示以及说明书

-p --platform 指定平台 （ALOS，ENVISAT，ERS-2等等）
-a 绝对轨道号
-r 相对轨道号 目前尚不清楚其定义
-i --intersectWith 几何相交 ，后跟WKT格式的字符串，如
–intersectsWith=“POLYGON((-156.42224080873712 20.480539890660793,-156.42224080873712 18.7086915511088,-154.42712388945802 18.7086915511088,-154.42712388945802 20.480539890660793,-156.42224080873712 20.480539890660793))”
-f 图框编号，目前尚不清楚其定义
-s -e 起止日期，格式为2008-01-15

---

更多参数可以直接运行查看

可以按照给出的参数提示进行需求定制，进而运行
运行py文件进行查询工作：

!$HOME/work/SSARA/ssara_federated_query.py --platform=ENVISAT --relativeOrbit=84 --frame=2943,2961 --start=2010-03-27 --end=2010-05-03 --collection="WInSAR ESA" --print 

这是一个搜索命令，使用--print打印得到结果

Running SSARA API Query:  https://web-services.unavco.org/brokered/ssara/api/sar/search?platform=ENVISAT&relativeOrbit=84&frame=2943%2C2961&start=2010-03-27&end=2010-05-03&collectionName=WInSAR+ESA&processingLevel=L0%2CL1.0%2CSLC
SSARA API query: 16.821295 seconds
###########################
ASF: The read operation timed out
###########################
Found 4 scenes
curl -X GET "https://portal.opentopography.org/API/globaldem?north=33.625395&south=31.259826&east=-114.288114&west=-116.189303&outputFormat=GTiff&demtype=SRTMGL1_E" -H "accept: */*" -o dem.wgs84.tif
WInSAR ESA,ENVISAT,42222,2010-03-28T17:50:04.291Z,2010-03-28T17:50:20.571Z,84,2943,2943,IM,S2,D,R,VV,https://imaging.unavco.org/data/sar/lts/winsar/ENV2/84/2943/ASA_IM__0CNPDE20100328_175004_000000162088_00084_42222_9504.N1
WInSAR ESA,ENVISAT,42222,2010-03-28T17:50:19.380Z,2010-03-28T17:50:35.660Z,84,2961,2961,IM,S2,D,R,VV,https://imaging.unavco.org/data/sar/lts/winsar/ENV2/84/2961/ASA_IM__0CNPDE20100328_175019_000000162088_00084_42222_9504.N1
WInSAR ESA,ENVISAT,42723,2010-05-02T17:50:01.820Z,2010-05-02T17:50:18.100Z,84,2943,2943,IM,S2,D,R,VV,https://imaging.unavco.org/data/sar/lts/winsar/ENV2/84/2943/ASA_IM__0CNPDE20100502_175001_000000172089_00084_42723_0354.N1
WInSAR ESA,ENVISAT,42723,2010-05-02T17:50:16.910Z,2010-05-02T17:50:33.190Z,84,2961,2961,IM,S2,D,R,VV,https://imaging.unavco.org/data/sar/lts/winsar/ENV2/84/2961/ASA_IM__0CNPDE20100502_175016_000000172089_00084_42723_0354.N1

注意，这里出现的curl命令是SSARA自动从openTopography中获取的你的研究区域的dem，可以将其单独复制到命令行进行下载，注意在尾部添加自己的API，如上方的输出，我们可以运行：

!curl -X GET "https://portal.opentopography.org/API/globaldem?north=33.625395&south=31.259826&east=-114.288114&west=-116.189303&outputFormat=GTiff&demtype=SRTMGL1_E&API_Key=APIAPI" -H "accept: */*" -o dem.wgs84.tif

注意把APIAPI替换为自己的api key即可运行下载

---

2023/11/15最近使用过程中出现了一些bug，ASF数据总是time out，目前尚未找到解决方案，可能是ASF官网维护中
2023/11/18问题已解决，就是因为ASF官方在维护数据，当前完全正常使用

这种方法整体是比较方便的，因为可以搜索到足够宽泛的数据来源，但是可能会出现小bug。

2.使用SSARA GUI界面

SSARA除了命令行下载模式之外还给出了网站支持GUI浏览：SSARA GUI 和其他数据网站功能基本相似。