时序地表沉降数据是区域地质灾害监测、工程安全评估及地表形变演化研究的重要指标。本研究利用欧空局哨兵一号（Sentinel-1）SAR 数据影像分辨率为 5m（距离向）×20m（方位向），结合数字高程模型（DEM）、轨道辅助数据等数据源，基于 SBAS-InSAR 技术对多期影像进行时序干涉处理，生成研究区 2018-2023 年 的时序地表沉降数据产品。本数据具有以下优势：时空连续性强，克服传统 InSAR 技术因时空基线过大导致的失相干问题。
", "ds_source": "Sentinel-1A 影像从欧洲航空局下载。DEM数据DEM是美国地质调查局（USGS）发布的30m较高分辨率SRTM数据，下载地址为：https://dwtkns.com/srtm30m/。 精密定轨星历数据。下载前在筛选时间的过程中需要注意文件命名格式，文件大小一般为4.4M。下载地址为：https://s1qc.asf.alaska.edu/aux_poeorb/。
", "ds_process_way": "（1）对 Sentinel-1A 等 SAR 影像进行裁剪、镶嵌，将 DEM 数据转换为适配格式，并导入精密轨道数据以减小轨道误差影响；（2）根据时间和空间基线阈值（如时间基线≤90 天、空间基线≤临界基线的 45%），生成干涉像对组合，确保影像间相干性；（3）对干涉像对进行多视处理，采用 Goldstein 滤波抑制噪声，通过三角网格最小费用流法（Delaunay MCF）进行相位解缠，获取差分干涉图和解缠结果；（4）第一次反演优化相位并估算初始形变速率，第二次反演通过大气滤波（高通窗口 365 天、低通窗口 1600 米）去除大气相位影响，得到精确的线性形变速率和时序形变结果；（5）将斜距坐标系下的结果转换为 WGS84 地理坐标系，输出分辨率与干涉处理一致
", "ds_quality": "基线长度设置为 200 d，空间基线阈值设置为理论空间基线的 50%。后续操作中将相干性较低的干涉对进行移除，这样即确保了每一景影像至少连接了 4 个干涉对，有足够的数据参与计算，又避免了干涉对完全失相干所导致的形变精度降低
", "ds_acq_start_time": "2018-01-12 00:00:00", "ds_acq_end_time": "2023-04-04 00:00:00", "ds_acq_place": "木里煤田四号坑", "ds_acq_lon_east": 99.15416666666667, "ds_acq_lat_south": 38.12777777777778, "ds_acq_lon_west": 99.12277777777777, "ds_acq_lat_north": 38.14361111111111, "ds_acq_alt_low": null, "ds_acq_alt_high": null, "ds_share_type": "apply-access", "ds_total_size": 20515476, "ds_files_count": 0, "ds_format": "*..hdr,.*sml", "ds_space_res": "", "ds_time_res": "日", "ds_coordinate": "无", "ds_projection": "", "ds_thumbnail": "05e0ee43-9b9b-42be-aaf4-9726b88e89a8.jpeg", "ds_thumb_from": 0, "ds_ref_way": "", "paper_ref_way": "", "ds_ref_instruction": "", "ds_from_station": "", "organization_id": "5b99d600-008a-4069-8fc3-7adb9c3f2f8b", "ds_serv_man": "徐培耘", "ds_serv_phone": "13259922729", "ds_serv_mail": "xupy@xust.edu.cn", "doi_value": "", "subject_codes": [ "170.45" ], "quality_level": 0, "publish_time": "2026-06-10 10:03:27", "last_updated": "2026-06-10 10:03:27", "protected": false, "protected_to": "2027-08-20 00:00:00", "lang": "zh", "cstr": "", "i18n": { "en": { "title": "Surface subsidence data set of No. 4 pit in Muli coalfield, Qinghai Province (2018-2023)", "ds_format": "*..hdr,.*sml", "ds_source": "Sentinel-1A images were downloaded from the European Aviation Agency. DEM data DEM is 30m high-resolution SRTM data released by the U.S. Geological Survey (USGS). The download address is https://dwtkns.com/srtm30m/. Precision orbit determination ephemeris data. You need to pay attention to the file naming format during the filtering time before downloading. The file size is generally 4.4M. The download address is: https://s1qc.asf.alaska.edu/aux_poeorb/.
", "ds_quality": "The baseline length was set to 200 days, and the spatial baseline threshold was set to 50% of the theoretical spatial baseline. In subsequent operations, interference pairs with low coherence are removed, which ensures that each image is connected to at least 4 interference pairs, sufficient data is available to participate in the calculation, and avoids deformation caused by complete loss of interference pairs. Accuracy reduction
", "ds_ref_way": "", "ds_abstract": "Time-series surface subsidence data is an important indicator for regional geological disaster monitoring, engineering safety assessment and surface deformation evolution research. This study uses ESA Sentinel-1 SAR data image resolution of 5m (range direction) ×20m (azimuth direction), combined with digital elevation model (DEM), orbit auxiliary data and other data sources, based on SBAS-InSAR technology, multi-period images are subjected to time-series interference processing to generate time-series surface subsidence data products in the study area from 2018 to 2023. This data has the following advantages: strong spatiotemporal continuity, overcoming the loss of coherence caused by the large spatiotemporal baseline of traditional InSAR technology.
", "ds_time_res": "", "ds_acq_place": "Muli Coalfield Pit No. 4", "ds_space_res": "", "ds_projection": "", "ds_process_way": "(1) Clipping and mosaing SAR images such as Sentinel-1A, converting DEM data into an adapted format, and importing precision orbit data to reduce the impact of orbit errors;(2) Based on temporal and spatial baseline thresholds (For example, the temporal baseline is ≤90 days, and the spatial baseline is ≤ 45% of the critical baseline), generate a combination of interference image pairs to ensure coherence between images;(3) Multi-view processing the interference image pairs, and Goldstein filtering is used to suppress noise. Through triangular grid minimum cost flow method (Delaunay MCF) performs phase unwrapping to obtain differential interferogram and unwrapping results;(4) The first inversion optimizes the phase and estimates the initial deformation rate, and the second inversion passes atmospheric filtering (High-pass window 365 days, low-pass window 1600 meters) Remove the atmospheric phase effect and obtain accurate linear deformation rate and time-series deformation results;(5) Convert the results in oblique coordinate system to WGS84 geographical coordinate system, and output The resolution is consistent with interference processing
", "ds_ref_instruction": "" } }, "submit_center_id": "ncdc", "data_level": 0, "recommendation_value": 0, "license_type": "https://creativecommons.org/licenses/by/4.0/", "doi_reg_from": "reg_local", "cstr_reg_from": "reg_local", "doi_not_reg_reason": null, "cstr_not_reg_reason": null, "is_paper_in_submitting": false, "belong_to_nieer": false, "ds_topic_tags": [ "地表沉降", "Sentinel-1", "SBAS-InSAR", "长时序" ], "ds_subject_tags": [ "地理学" ], "ds_class_tags": [], "ds_locus_tags": [ "青海省", "木里煤田四号坑" ], "ds_time_tags": [ 2018, 2019, 2020, 2021, 2022, 2023 ], "ds_contributors": [ { "true_name": "徐培耘", "email": "xupy@xust.edu.cn", "work_for": "西安科技大学", "country": "中国" }, { "true_name": "王锴", "email": "3114412346@qq.com", "work_for": "西安科技大学", "country": "中国" } ], "ds_meta_authors": [ { "true_name": "王锴", "email": "3114412346@qq.com", "work_for": "西安科技大学", "country": "中国" } ], "ds_managers": [ { "true_name": "徐培耘", "email": "xupy@xust.edu.cn", "work_for": "西安科技大学", "country": "中国" } ], "category": "其他" }