Sentinel-1 SAR数据准备

这篇笔记记录博士论文期间Sentinel-1A雷达数据的准备工作， 包括Google Earth Engine中下载并与处理的基于已获取的GNSS点处的SAR数据，研究区的SAR数据下载等。

1. 下载小尺度SAR数据

可获取的GNSS点数据

东部形变数据服务中心网站

  • 在其网站数据中选择4个站点，分别为AHBB, JSLY, SDJX, SDLY进行实验。

选取的GNSS点

Name lat lon BOX1 BOX2 BOX3 BOX4 Name in Chinese
1 HECX 116.93 38.46 116.9234 38.4546 116.9401 38.4633 沧州
2 SDJX 116.35 35.43 116.3337 35.4152 116.3749 35.4401 济宁
3 XIAA 108.99 34.18 108.9522 34.1565 109.0355 34.1968 西安
4 AHBB 117.30 32.90 117.2899 32.8936 117.3106 32.9054 蚌埠

stackSentinel.py -s $SLC_PATH -w $INSAR_WORK_DIR -a $AUX_PATH -d $DEM_PATH/demLat_N38_N39_Lon_E116_E117.dem.wgs84 -o $ORBIT_PATH -n 2 -b ‘32.80 32.91 117.28 117.32’ --useGPU -C geometry -c 3 -z 1 -r 4 -f 0.6 --num_proc 4

ASF下载burstSAR数据

  • 由于GEE并不支持直接下载SLC数据，所以只能选择在ASF网站上先有选择的下载所需GNSS点处的部分SAR数据，再进行裁减等操作。

Sentinel-1 SLC data cannot currently be ingested, as Earth Engine does not support images with complex values due to inability to average them during pyramiding without losing phase information.

  • 连云港附近JSLY点的SAR数据使用aria2c下载：
1
2
3
4
5
6
7
8
cd ~/Documents/00_data/sentinel/JSLY
aria2c --http-user=mengyuchile --http-passwd=5212109Mengle --check-certificate=false --max-concurrent-downloads=2 --allow-piece-length-change=false -M asf-datapool-results-2025-07-16_03-57-07.metalink

aria2c --http-user=mengyuchile --http-passwd=5212109@MENGle --check-certificate=false --max-concurrent-downloads=2 --allow-piece-length-change=false -M asf-datapool-results-2025-07-16_03-57-07.metalink

aria2c --http-user=mengyuchile --http-passwd=5212109Mengle -c --continue=true --check-certificate=false -x 2 -s 36 -M asf-datapool-results-2025-07-15_03-33-33.metalink

aria2c --http-user=mengyuchile --http-passwd=5212109Mengle --check-certificate=false --allow-piece-length-change=false -M asf-datapool-results-2025-07-16_06-04-35.metalink
1
2
3
4
5
6
7
8
# 下载失败的数据
6f8a65|ERR |   1.7MiB/s|/home/yuchi/Documents/00_data/sentinel/JSLY/S1A_IW_SLC__1SDV_20231104T100431_20231104T100458_051066_062861_3C82.zip
b1425e|ERR |   1.4MiB/s|/home/yuchi/Documents/00_data/sentinel/JSLY/S1A_IW_SLC__1SDV_20231023T100431_20231023T100451_050891_062262_A1B5.zip
12f3ba|ERR |       0B/s|/home/yuchi/Documents/00_data/sentinel/JSLY/S1A_IW_SLC__1SDV_20231011T100430_20231011T100457_050716_061C6C_E547.zip
14361c|ERR |       0B/s|/home/yuchi/Documents/00_data/sentinel/JSLY/S1A_IW_SLC__1SDV_20230929T100431_20230929T100451_050541_06166A_5C06.zip
fd0d6c|ERR |       0B/s|/home/yuchi/Documents/00_data/sentinel/JSLY/S1A_IW_SLC__1SDV_20230917T100431_20230917T100450_050366_06106D_9A3F.zip
a16a2c|ERR |       0B/s|/home/yuchi/Documents/00_data/sentinel/JSLY/S1A_IW_SLC__1SDV_20230905T100430_20230905T100450_050191_060A81_D64A.zip

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
Failures: 17 files
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20230731T100543_20230731T100610_049666_05F8E6_9783/IW1/VV/0.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20211208T100528_20211208T100556_040916_04DBDE_EEF4/IW1/VV/3.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20211126T100529_20211126T100556_040741_04D5C6_4C10/IW1/VV/3.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1B_IW_SLC__1SDV_20210629T100444_20210629T100519_027570_034A8A_86A4/IW1/VV/2.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210717T100531_20210717T100558_038816_049483_AB7D/IW1/VV/0.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210717T100531_20210717T100558_038816_049483_AB7D/IW1/VV/1.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210705T100530_20210705T100557_038641_048F47_B75E/IW1/VV/0.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210705T100530_20210705T100557_038641_048F47_B75E/IW1/VV/1.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1B_IW_SLC__1SDV_20210629T100444_20210629T100519_027570_034A8A_86A4/IW1/VV/1.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1B_IW_SLC__1SDV_20210723T100448_20210723T100516_027920_0354D7_FF8E/IW1/VV/0.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210623T100529_20210623T100556_038466_0489FC_ED15/IW1/VV/1.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210106T100526_20210106T100553_036016_043851_6E86/IW1/VV/1.tiff
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20230731T100543_20230731T100610_049666_05F8E6_9783/IW1/VV/0.xml
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20221203T100539_20221203T100602_046166_0586F1_5712/IW1/VV/0.xml
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210822T100533_20210822T100600_039341_04A572_0643/IW1/VV/1.xml
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210211T100525_20210211T100552_036541_044A97_AA08/IW1/VV/0.xml
        - https://sentinel1-burst.asf.alaska.edu/S1A_IW_SLC__1SDV_20210130T100525_20210130T100552_036366_04448A_0373/IW1/VV/1.xml
Average Rate: N/A (not enough data for calculation)

  • 使用ASF网站提供的python脚本下载burst数据：

asf python脚本下载

因自带python脚本是单线程，下载速度太慢，使用gemini优化代码，实现多线程下载，详细修改后代码见博客。

1
2
cd ~/Documents/00_data/sentinel/HECX1
python asf_sentinel_burst_download.py

burst2safe下载ASF burst数据

  1. 安装配置burst2safe

经过测试，直接使用conda insatll下载的burst2safe有一个文件运行报错，建议直接下载开发者模式。

1
2
3
4
5
git clone https://github.com/ASFHyP3/burst2safe.git
cd burst2safe
mamba env create -f environment.yml
mamba activate burst2safe
python -m pip install -e .

以GNSS点HECX为例：

1
2
conda activate burst2safe
burst2safe S1_146285_IW1_20210106T100527_VV_6E86-BURST S1_146286_IW1_20210106T100530_VV_6E86-BURST
  1. 在明确知道自己需要的SLC burst数据的name字段的情况下，可使用burst2safe下载ASF burst数据并直接合成为SNAP, ISCE等能处理的SAFE格式数据。
1
burst2safe S1_136231_IW2_20200604T022312_VV_7C85-BURST S1_136232_IW2_20200604T022315_VV_7C85-BURST

burst2safe结果

  1. 也可以直接根据研究区下载多幅SLC burst数据：
1
2
3
4
5
6
7
8
9
10
11
12
13

# XIAA
burst2stack --rel-orbit 84 --start-date 2021-01-01 --end-date 2024-01-01 --extent 108.9522 34.1565 109.0355 34.1968 

# SDJX
burst2stack --rel-orbit 142 --start-date 2021-01-01 --end-date 2024-01-01 --extent 116.3337 35.4152 116.3749 35.4401 

# AHBB
burst2stack --rel-orbit 142 --start-date 2021-01-01 --end-date 2024-01-01 --extent 117.2899 32.8936 117.3106 32.9054 

# HECX
burst2stack --rel-orbit 69 --start-date 2021-01-01 --end-date 2024-01-01 --extent 116.9234 38.4546 116.9401 38.4633

burst2stack结果

burst2stack命令中的–rel-orbit可以在任意一幅.xml文件中搜索relativeOrbitNumber字段来获取，比如<safe:relativeOrbitNumber type="start">142</safe:relativeOrbitNumber>

  1. 也可以直接从在ASF网站上下载sentinel-burst数据，每一幅一影像包含.tiff与.xml文件，使用local2safe直接合成SAFE文件。
  • local2safe 需要设置 json_tree_path 文件:
1
2
3
4
5
6
7
8
# 激活conda环境
conda activate burst2safe

# 例如进入文件夹
cd ~/Documents/00_data/sentinel/AHBB

# 创建json文件
gedit local2safe.json
1
2
3
4
S1_146285_IW1_20210319T100526_VV_4ECC-BURST.tiff
S1_146286_IW1_20210319T100528_VV_4ECC-BURST.tiff
S1_146285_IW1_20210331T100526_VV_52D5-BURST.tiff
S1_146286_IW1_20210331T100529_VV_52D5-BURST.tiff
  • 并在文件中编辑：
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
{
    "S1A_IW_SL1__1_DV_20210118T100525_20210118T100558_036191_043E75_642B.SAFE":{
        "IW1":{
            "VV":{
                "1":{
                    "DATA":"S1_146286_IW1_20210118T100529_VV_4C21-BURST.tiff",
                    "METADATA":"S1_146286_IW1_20210118T100529_VV_4C21-BURST.xml"
                }
            }
        }
    },
    "S1A_IW_SL1__1_DV_20210118T100525_20210118T100558_036191_043E75_642B.SAFE":{
        "IW1":{
            "VV":{
                "0":{
                    "DATA":"S1_146285_IW1_20210118T100527_VV_4C21-BURST.tiff",
                    "METADATA":"S1_146285_IW1_20210118T100527_VV_4C21-BURST.xml"
                }
            }
        }
    }
}
  • 编辑好json文件后在bash中运行：
1
2
3
# 运行local2safe

local2safe local2safe.json

local2safe结果

!注意：使用local2safe只会生成一个SAFE文件，待解决。

  • 在output_dir中得到S1A_IW_SLC__1SSV_20200604T022312_20200604T022315_032861_03CE65_C78F.SAFE, 大功告成！
Data Preparation
#Sentinel-1 #SAR Data #Google Earth Engine #ASF #GNSS
Sentinel-1 SAR数据准备
https://mengyuchi.gitlab.io/2025/07/14/Sentinel-1-SAR数据准备/
Author
Yuchi Meng
Posted on
July 14, 2025
Licensed under