isce2学习（三）下载哨兵数据、POD精轨数据、DEM数据、AUX_CAL辅助文件下载

isce2学习（三）下载哨兵数据、POD精轨数据、DEM数据、AUX_CAL辅助文件下载

0、介绍

注：ISCE支持的哨兵哨兵数据类型有：

• IW SLC——干涉宽幅模式（TOPS Mode）的斜距单视复数产品

• EW SLC——超宽幅模式（TOPS Mode）的斜距单视复数产品

• SM SLC ——条带模式的斜距单视复数产品

但是在中国境内，哨兵只提供IW模式的数据。

查看哨兵数据的网址：https://asf.alaska.edu/data-sets/sar-data-sets/sentinel-1

具体格式如下：

任何关于ISCE的问题可以到ISCE论坛进行交流，ISCE论坛网址：http://earthdef.caltech.edu/projects/isce_forum/boards

以下是用ISCE的topsStack处理哨兵数据的流程。

1 、需要的数据下载

首先建立一个工程文件夹，下面分别建立SLC文件夹，DEM文件夹，轨道数据（orbits）文件夹（ POD Precise Orbit Ephemerides即为POD精密定轨星历数据，这是最精确的轨道数据，但该数据距离GNSS下行21天之后才可以使用。每天会产生一个文件，每个文件覆盖26个小时（一整天24小时加上一天开始前1小时和一天结束后的1小时），定位精度优于5厘米。文件后缀为aux_poeorb；POD Restituted Orbit（POD回归轨道数据），这是比较精确轨道数据，文件在接收到GNSS数据的3小时内产生。该文件覆盖一个卫星的轨道，从升序节点ANX加593 OSV重叠于卫星轨道的时间跨度之前。定位精度优于10厘米，文件后缀为aux_resorb）。和AUX文件（Calibration Auxiliary Data，校准辅助数据，文件后缀为aux_cal）。其中SLC数据，AUX_CAL数据和AUX_Poeorb数据都可以在ASF官网下载。

2、SLC哨兵数据下载

到ASF官网下载ASF官网

3、orbit精密轨道文件下载

（1）eesa下载链接为：https://scihub.copernicus.eu/gnss/#/home

登陆账号：gnssguest 密码：gnssguest

下载精密轨道文件

（2）ESA上的下载地址：
• https://scihub.copernicus.eu/gnss/#/home

（3）ASF官网上的下载地址：
• https://s1qc.asf.alaska.edu/aux_poeorb/

其中第二个网站的账号密码均为：gnssguest

下面以成像时间为20160705（S1A_IW_SLC__1SDV_20160705T222043_20160705T222116_012019_0128FB_1232.SAFE）数据为例，其对应的轨道数据如下：

如上图所示：可看到每个轨道数据文件都有3个日期，1:精密轨道数据发布日期，2：数据成像时间的前一天，3：数据成像时间的后一天，网站右边查询过滤器的日期为成像时间的前一天，因此成像时间为20160705的数据文件对应轨道文件为S1A_OPER_AUX_POEORB_OPOD_20160725T121611_V20160704T225943_20160706T005943.EOF。点击文件用浏览器下载即可。

一般选择第三个网址进行下载https://s1qc.asf.alaska.edu/，

进去之后选择aux_poeorb, 页面显示如下：

文件名称露出来的时间是精密轨道数据发布日期，鼠标放在链接上，查看完整的文件名称，选择需要的数据，点击后出现一个网页，等待网页完全加载后，另存为，另存为之后浏览器回自动加.txt后缀，删除.txt后缀，保持后缀为.EOF即可。

以上为手动下载精轨文件自己添加到文件夹orbits里，也可以直接利用python脚本fetchOrbit.py进行下。但是欧空局在2020年2月改变了哨兵一号轨道的位置。因此需要使用ISCE2 GitHub上一个新版本的“fetchOrbit.py”，才能自动下载文件，它拥有这个新站点。

4、DEM数据下载

(1) *先设置**Earthdata 帐户**（**此处设置一次便可）*

您需要拥有 urs.earthdata.nasa.gov 用户名和密码，并且需要将LPDAAC应用程序包含在您的帐户中。在 Earthdata 帐户中将LPDAAC 应用程序添加到您的帐户：在此处登录到 earthdata： https://urs.earthdata.nasa.gov/home - 在配置文件上单击我的应用程序 - 单击"添加更多应用程序" - 搜索"LP DAAC" - 选择"LP DAAC 数据池"和"LP DAAC OpenDAP"并批准。https://lpdaac.usgs.gov/

最后得到如下

*（2）建立.netrc文件*

在DEM数据下载之前，先要在home目录下建立.netrc文件，

先进入home目录：

cd /home

在home目录下建立.netrc文件

 sudo touch .netrc

文件中存入的内容如下

进入如下界面，按a开始编辑，将

machine urs.earthdata.nasa.gov
login xxxx
password xxxxxx

添加进去，按esc 退出，按 ：wq退出

dem下载 运行：

dem.py -a stitch -b 34 38 119 122 -r -s 1 -c

（出现没有指令-c的错误，则是网址被墙，没下载完成，-c指令为剪裁并保存，输入的bbox值必须是整数）
下载DEM数据,并且会按照输入的bbox值(矩形框的位置信息青岛的bbox: 34 38 119 122

崂山区位置信息bbox:36 37 120 121)截取DEM：

运行成功之后显示的如上图，运行之后生成的结果如下图。

DEM拼接，然后执行：

rm demLat*.dem demLat*.dem.xml demLat*.dem.vrt

最后剩余文件：

拼接得到的DEM 文件也可以用gdal_translate 转换格式用到gmt中绘制地形图

gdal_translate -of GMT -b 1 demLat_N36_N37_Lon_E120_E121.dem.wgs84 dem.grd

拼接得到的DEM 文件也可以用gdal_translate 转换格式用到gmt中绘制地形图

gdal_translate -of GMT -b 1 demLat_N36_N37_Lon_E120_E121.dem.wgs84 dem.grd

5、AUX_CAL 数据下载

下载地址为：https://qc.sentinel1.copernicus.eu/aux_cal/
AUX_CAL辅助文件主要是用来天线方向图校正，IPF(Instrument Processing Facility)仪器处理设备如果是同一个版本的话则不需要辅助文件。主要使用的是2015年3月之前，也就是说当2015年3月之后和之前的图像生成干涉图时会导致相位偏移的变化，因此需要此文件进行校正。但是还是需要文佳夹Aux存在，后期命令必须。