基于多源遥感的黄河三角洲湿地动态监测_樊彦国.pdf

2023 年第 6 期樊彦国，等:基于多源遥感的黄河三角洲湿地动态监测引文格式:樊彦国，王杰，樊博文，等 基于多源遥感的黄河三角洲湿地动态监测J 测绘通报，2023(6):27-35 DOI:10 13474/j cnki 11-22462023 0164基于多源遥感的黄河三角洲湿地动态监测樊彦国1，王杰1，樊博文2，徐梓耀1(1 中国石油大学(华东)海洋与空间信息学院，山东 青岛 266580;2 哈尔滨工程大学水声工程学院，黑龙江 哈尔滨 150001)摘要:城乡建设的逐步推进与自然生态环境保护之间的矛盾愈发激烈，我国对黄河三角洲湿地区域的发展高度重视，迫切需要对黄河三角洲近年来的动态变化信息进行监测、提取与分析。本文利用 20002021 年 8 景 Landsat 光学影像与 20152021 年 34 景Sentinel-1 降轨 VV 极化 SAR 影像，对黄河三角洲湿地区域的宏观(土地利用)、微观(地表形变)信息进行提取，并结合社会经济、自然环境数据进行驱动因素分析。结果表明:20002021 年黄河三角洲湿地区域土地利用格局未发生显著变化，但土地利用程度逐步提高;20152021 年黄河三角洲湿地区域的地表形变程度空间分布极不均匀，最大沉降速率达43.643 mm/a，最大抬升速率为 0.136 mm/a;人类活动和社会经济是引起该区域地表信息变化的主导驱动因素。关键词:黄河三角洲;多源数据协同;动态监测;土地利用;地表形变中图分类号:P237文献标识码:A文章编号:0494-0911(2023)06-0027-09Dynamic monitoring of wetlands in the Yellow River delta based onmulti-source remote sensingFAN Yanguo1，WANG Jie1，FAN Bowen2，XU Ziyao1(1 College of Oceangraphy and Space Information，China University of Petroleum(East)，Qingdao 266580，China;2 School of Underwater Acoustic Engineering，Harbin Engineering University，Harbin 150001，China)Abstract:The contradiction between the gradual urban and rural construction and natural ecological environmental protection isbecoming more and more intense，and the state and government attach great importance to the development of the Yellow River deltawetland area Therefore，it is urgent to monitor，extract and analyze the dynamic change information of the Yellow River delta in recentyears In this paper，macroscopic(land use)and microscopic(surface deformation)information of the Yellow River delta wetlandregion is extracted and analyzed by using 8-view Landsat optical images from 2000 to 2021 and 34-view Sentinel-1 downscaled VVpolarized SAR images from 2015 to 2021，combined with socio-economic and natural environment data for driving factor analysis Theresults show that the land use pattern of the Yellow River delta wetland region did not change significantly between 2000 and 2021，butthe land use degree gradually increased;the spatial distribution of the surface deformation degree of the Yellow River delta wetlandregion is extremely uneven between 2015 and 2021，with the maximum sedimentation rate reaching 43.643 mm/a and the maximumuplift rate of 0.136 mm/a;human activities and socio human activities and social economy are the dominant drivers of surfaceinformation changes in the regionKey words:the Yellow River delta;multi-source data collaboration;dynamic monitoring;land use;surface deformation近年来，我国对黄河三角洲的高质量发展做出了大量的指导意见和规划方案1-4。当前黄河三角洲地区的发展正处于机遇与挑战并存的状态。目前，该地区已经具备了高效生态经济发展的基础条件，具有在更高的起点上进一步实现可持续高质量、高效率发展和建设文明生态环境的优势条件;但是人口增长、经济发展、技术进步等因素使得人类对土地的使用范围扩大、程度加剧。因此对黄河三角洲及其周边地区的地表动态变化进行提取、监测和分析显得尤为重要，不仅方便于该地区的土地管理工作，还能够为城市的发展规划提供科学依据。卫星遥感是目前大面积区域动态变化监测分析最常用的技术手段之一。已有大量学者分别利用光学遥感、雷达遥感数据对黄河三角洲及其周边72收稿日期:2022-08-08;修回日期:2023-04-19基金项目:自主创新项目-战略专项(22CX01004A-3);国家自然科学基金青年基金(42106215)作者简介:樊彦国(1965)，男，博士，教授，主要研究方向为遥感信息智能提取。E-mail:ygfan upc edu cn通信作者:王杰。E-mail:2394834585 qq com测绘通报2023 年第 6 期地区的土地利用及其驱动因素5-8、地表形变及其诱发因素进行相关研究9-13。本文首先利用改进的监督分类对黄河三角洲及其周边地区进行土地利用分类;然后利用统计年鉴中的数据对黄河三角洲土地利用变化进行驱动力分析;最后利用 SBAS-InSAR 技术获取研究区域的地表形变信息，并探讨其诱发因素。1研究区域与数据来源1.1研究区域概况现代黄河三角洲以宁海为起点，北起套尔河口，南至支沟口，占地面积约 5400 km2。本文对涵盖黄河三角洲及其周边的河口区、垦利区、东营区等地区展开研究。1.2数据来源本文利用了 20002021 年 8 景 Landsat 系列光学遥感影像和 20152021 年 34 景 Sentinel-1SAR 遥感影像;所使用的统计数据均源自东营市统计年鉴 202114;使 用 分 辨 率 为 30 m 的ASTER GDEM V3 版 DEM 数据作为 SAR 影像处理参考。2光学遥感用于黄河三角洲湿地动态监测首先对光学遥感影像进行预处理，其次将最大似然法分类结果与 NDVI、NDBI 及 NDWI 指数进行协同判断以提高分类精度，算法流程如图 1 所示，其中，T1、T2、T3分别为不同年份植被、水体、建筑提取阈值。图 1本文监督分类算法流程2.1黄河三角洲土地利用变化特征分析监督分类与各地类空间分布中心计算结果如图 2 所示;各时期地类面积时间序列如图 3 所示;土地利用转移矩阵可视化结果如图 4 所示。由图 2 可知，20002009 年研究区植被和裸地覆盖区域分布极不均匀，在研究区域内呈现南北两端分布;20122021 年植被和裸地覆盖区域的分布逐渐转向均匀分布。滩涂的分布具有显著空间规律性，主要部分位于海陆交界处，呈现弯曲条带状。20002015 年 建 筑 区 域 分 布 较 为 离 散，20182021 年建筑区域面积明显增加，且集中分布在研究区域的中部和南部地区。水体的空间分布特征和面积无明显变化。20002021 年黄河三角洲滨海地带海岸线也发生了一定的变化，有由陆地缓慢向822023 年第 6 期樊彦国，等:基于多源遥感的黄河三角洲湿地动态监测海域扩张的趋势。图 28 个时期遥感影像监督分类及各地类空间分布重心图 3不同地类面积时间序列各地类的空间分布中心除水体外均位于陆地区域。水 体、滩 涂 区 域 重 心 位 置 较 为 固 定，但20002021 年滩涂区域重心有缓慢向东移动、靠近水体重心的趋势。20002009 年植被区域和裸地区域的重心呈南北两端分布，且两者位置有较大位移交替;20122021 年植被区域重心逐渐向研究区域的几何中心靠近，而裸地区域的重心逐渐向东南方向移动，两者之间的距离也越来越大，表明随着时间的推移，研究区域的植被分布面积不断扩大，生态环境得到明显改善。而建筑区域重心位置变化与植被、裸地重心位置相关，20002009 年建筑区域重心位置靠近植被区域重心位置，而 20092021 年建筑区域中心位置则靠近裸地区域重心位置，表明早期城镇向植被区域扩张，而在 2012 年后城镇化的扩张地区逐渐转为裸地。由图 3 可知，20002021 年研究区域水体区域面积基本保持在 9000 km2左右，面积波动幅度较小;建筑区域面积持续增加，其增长速率逐渐加快，但 2006 年建筑物面积陡增，主要是因为 2006 年选取的 Landsat 影像成像时间为 3 月 16 日，研究区域处于早春季节，裸地覆盖居多，且该影像中裸地与建筑物的区分度较小，故导致该年建筑区域的面积出现陡增现象;裸地面积呈现下降趋势，而植被面积呈现增长趋势，但 2006 和 2018 年影像成像时间分别为 3 月 16 日和 11 月 12 日，此时研究区域分别处于早春、晚秋季节，根据当地自然环境气候特征，导致植被区域在光学遥感影像中呈现出与裸地相似的特征，故这两个时相监督分类结果中出现植被减少，裸地骤增的异常现象;滩涂面积以平稳缓慢的速度较少。裸地、植被、建筑这 3 种地类的分类较易受到季节气候环境的干扰。由图 4 可知，20002018 年(2006 和 2018 年除外，因有季节影像)，建筑区域转为裸地区域的面积均小于裸地区域转为建筑区域的面积，在不考虑误差的情况下，可总体认为裸地区域向建筑区域转化，且转化速率逐渐变慢;但 20182021 年建筑区域与裸地区域的相互转换近似相等，可知此时的建筑区域面积近似达到饱和状态。裸地区域的转化起伏较大，主要向植被区域和建筑区域转换。此外，滩涂区域除了自身转换外，大部分转化为植被区域，也有少量植被区域转为建筑区域。综上所述，20002021 年各地类均有不同程度的转变，其中裸地区域转为建筑、植被区域较多，少量滩涂区域转为建筑和水体区域，水体区域在转换过程中则表现出较为稳定特征，由水体转换成其他地类的面积较少。92测绘通报2023 年第 6 期图 4土地利用面积转移矩阵由于所用的光学遥感影像分辨率为 30 m，此外，2006 与 2012 年的遥感影像是在 Landsat 7 出现故障之后拍摄的，影像上存在很多缺失的条纹，虽然使用了一定的方法对其进行修补，但还会存在一部分损失的信息。加之建筑区域与裸地的光谱特征较为相似，导致选择的训练样本无法保证完全正确，再考虑所使用的监督分类方法具有一定的局限性，因此对建筑、裸地区域的提取较为困难，导致监督分类结果中存在一定的误差，从而土地利用转移矩阵中出现部分类似于水体转建筑、建筑转植被、植被转水体等一系列不合理的土地利用转型现象。2.2黄河三角洲土地利用驱动力分析土地利用变化会受到自然环境、社会经济、人类活动等多种因素影响。为对该地区的高效生态经济发展提供决策依据，从定量角度分析该区域的土地利用结构特征和驱动因子。在参考文献 1517 的基础之上，结合实际分析，从东营市 2021 统计年鉴中选取人口数量、年均气温、生产总值等20 个驱动因子指标，并将其划分为人类活动、自然环境和社会经济 3 大类。通过相关性分析(结果如图 5 所示)选取出 15 个强相关因子，进行建模分析。水体、滩涂区域面积变化不明显，建筑和植被区域面积具有相似的特征，且两者面积与裸地面积互补。由于植被区域面积变化是黄河三角洲土地利用变化的关键，因此以植被区域面积作为因变量组 Y，以代替研究区域的土地利用变化特征。考虑到各驱动因子内部之间也可能存在相关性，为了能够在土地利用驱动模型中引入较少的自变量，需要在上述工作的基础之上，对 15 个驱动因子进行主成分分析，分析结果见表 1。然后选取前两个主成分，计算出这 15 个因子在这两个新的主成分中的载荷系数(见表 2)，并以此评价各驱动因子之间的相对重要性，从而选取主导驱动因子。图 5各驱动因子与地类之间的相关系数计算由表 1 可看出，第一主成分、第二主成分、第三主成分的方差贡献占比分别为 91.299%、6.189%、1.512%，其累计贡献为 99.000%。由表 2 可看出，在第一主成分中，这 15 个驱动因子中除了公路里程、废水排放、进出口总值之外，其余驱动因子的载荷系数均较大，尤其是社会经济类驱动因子;在第二主成分中则与人类活动相关的驱动因子，如人口数量、城镇人均居住面积、农业机械总动力及废水排放等的载荷系数较高。因此，东营市的社会经济发展和人类活动是影响黄河三角洲土地利用的主要驱动因素。032023 年第 6 期樊彦国，等:基于多源遥感的黄河三角洲湿地动态监测表 1黄河三角洲土地利用变化驱动因子主成分分析结果主成分特征值方差贡献占比/(%)累计贡献/(%)一13.69591.29991.299二0.9286.18997.488三0.2271.51299.000四0.0870.57899.577五0.0360.23899.815六0.0180.11899.933七0.0100.067100.000八8.56E-315.71E-30100.000表 2各驱动因素在第一、第二主成分中载荷系数驱动因子第一主成分第二主成分人口数量0.9640.931城镇人均居住面积0.9620.965农业机械总动力0.9750.985工业企业个数0.9340.480公路里程0.2560.988废水排放0.2000.645全市生产总值0.9800.050第一产业0.9670.129第二产业0.9690.089第三产业0.9640.220人均生产总值0.9600.017进口总值0.2470.343出口总值0.2640.128一般公共预算收入0.9660.104一般公共预算支出0.9670.102以植被区域面积作为因变量，对不显著的驱动因子进行提取，选取对土地利用影响最大的 6 个因子，分别为人口数量、全市生产总值、农业机械总动力、第二产业、第一产业、一般公共预算支出，进行多元线性回归，得到的最终模型为Y=19.269X1+3.209X2+2.658X3+0.256X4+1.586X5+0.585X6+2 308.795式中，Y 为黄河三角洲区域植被区域面积;X1 X6分别为上述 6 个驱动因子。模型的拟合度 R2为0.996，通过 1%显著检验。从建模结果可看出，这6 个主导驱动因子均与植被区域面积呈正相关，人口数量、全市生产总值及农业机械总动力对植被区域面积的影响较大。通过对上述结果进行分析，可以得到黄河三角洲地区植被区域面积变化驱动机制，即黄河三角洲土地利用机制，如下:(1)社会经济对黄河三角洲土地利用变化的驱动。社会经济发展是土地开发利用的背景条件，同时两者之间相互影响、相互制约。2000 年以来，黄河三角洲地区工业发展迅速，第二产业和旅游业的发展，极大程度上刺激了土地利用的转型。全市GDP、人均收入、农业生产总值、能源生产总值、农业生产总动力等指标均有明显上升趋势，提高了土地的生产力。从相关性分析和主成分分析的结果来看，全市区生产总值、人均生产总值、第一、第二、第三产业生产总值等社会经济因子是第一主成分的主导因子，说明在 20002021 年，社会经济发展对植被区域面积变化有着重要的驱动影响。(2)人类活动对黄河三角洲土地利用变化的驱动。人口数量的增长和人类对土地的改造、治理等干预也会造成土地利用的转型。据东营市 2021 统计年鉴，20002021 年，东营市常住人口增长量约为 40 万人，城镇化率由 48.27%上升至 71.45%，非农业人口大幅增长，这主要是由城市扩张的经济发展所导致的，城市扩张和工业发展势必导致农业用地向非农业用地的转化，因此城市化建设也是土地利用变化的一个重要驱动因素。此外，东 营 市20002021 年生活污染物排放等指标下降，这表明在过去 21 年的时间里，富裕的生活使得人们的生活、工作方式发生了明显改变。经济发展带动了农民观念意识的发展，推动了农业的发展，同时也使得黄河三角洲土地利用结构方式发生显著变化。(3)自然环境对该区域的土地利用转型作用不明显。3雷达遥感用于黄河三角洲湿地动态监测使用 共 计 34 景 Sentinel-1 SLC 数 据，进 行SBAS-InSAR 地表形变反演(如图 6 所示)。选取2018 年 8 月 20 日的 SAR 影像作为超级主影像，时间基线阈值设置为一年(366 d)，空间基线阈值设置为45%，最终生成了 161 个干涉像对，生成的时空基线连接情况如图 7 所示。图 6SBAS-InSAR 研究区域13测绘通报2023 年第 6 期图 7干涉像对时空基线连接图3.1黄河三角洲地表形变特征分析通过 SBAS-InSAR 时序分析后反演得到黄河三角洲及周边地区的地表形变结果，为分析方便，在形变结果中画出形变显著区，如图 8 所示，其中正值表示地面抬升，负值表示地面沉降。图 820152021 年 SBAS-InSAR 解算结果由图 8 可知:20152021 年黄河三角洲湿地区域的地表形变主要表现为沉降，整个研究区域的累计形变量范围为256.4428.101 mm，地面沉降平均沉降速率为 10.602 mm/a，最大沉降速率为43.343 mm/a，最大抬升速率为 0.136 mm/a;有极少地区发生轻微抬升现象。利津县和博兴县内的地表形变区域空间分布均匀，大部分属于轻微沉降和中等沉降，利津县整体的沉降速率约为 12.265 mm/a，博兴县整体的沉降速率约为 3.560 mm/a。河口区、垦利区、寿光市的沉降程度在空间分布上存在显著的差异性。河口区西部一大半形变区域为展现出轻微沉降和中等沉降的特点，沉降速率范围为 9.593 14.456 mm/a;在靠近黄海的东部地区，绝大多数的形变 为 严 重 沉 降 和 特 严 重 沉 降，沉 降 均 大 于20 mm/a。垦利区的地表沉降在东北西南方向上体现出差异，东北沿海地区大部分为水田、植被覆盖，为未解算区域;而在西南部的城区内为中等沉降，大部分地区的沉降速率 10 mm/a 以上。寿光市的地表形变差异体现在南北方向上，寿光市南部为植被山地区域，属未结算区;北部多沿海滩涂、裸地，属轻微沉降，但其中有部分区域的沉降速率大于10 mm/a。东营区、广饶县及研究区域东南角区域的大部分区域均存在形变。其中东营区北部沉降现象较为明显，有连续轻微、中等沉降，沉降速率在020 mm/a 不等，最大累计沉降量高达 67.624 mm。广饶县的北部地区均为未解算区域，地表形变主要发生 在 南 部 城 区，市 中 心 的 最 大 沉 降 速 率 达29.049 mm/a，最大累计沉降量高达 121.566 mm。研究区域东南角的地表形变则体现出空间连续均匀分布的特点，主要的形变类型为轻微抬升和轻微沉降，最大累计沉降量为 56.212 mm，最大累计抬升量为 8.101 mm。黄河三角洲湿地区域地表发生的形变在空间上整体表现为不均匀，在离海水较远的内陆地区地表形变不连续，而沿海地区的地表形变在空间分布上比较连续。中部城区主普遍发生大范围的轻微、中等程度沉降;东南部的沿海地区则产生轻微沉降和轻微抬升;西南部主要体现为轻微、中等沉降，个232023 年第 6 期樊彦国，等:基于多源遥感的黄河三角洲湿地动态监测别区域达到了严重沉降的程度;北部地区存在较多未解算区，发生地表形变的区域其形变等级为轻微沉降、中等沉降居多，但仍然存在部分地区有特别严重沉降现象。3.2黄河三角洲地表形变诱发因素分析地表形变的诱发因素可以分为自然、人为两类。一般情况下，两者共同作用更容易引发大范围的地表形变。因此，为更好地分析地表形变的特征与诱发因素，对研究区域的形变区进行划分，从而得到 11 个显著形变区，如图 8(b)中 D1D11 区所示。3.2.1自然因素为重点验证断裂带对地表形变的影响，分别在远离城区的陈南断裂、益都断裂带方向及其垂直方向上选取若干高相干点(如图 9图 10 所示)，断层数据来源于文献 1819 等，通过比较这些高相干点位之间的形变速率证明断裂带对黄河三角洲地表形变的影响。图 11 为所选取高相干点的形变速率值，对其进行样条插值，绘制该方向上的形变速率曲线。其中，图 11(a)、(b)为与断裂带垂直方向上的高相干点，伴随着相干点与断层之间空间距离的不同，其形变速率也发生着不同程度的变化，在断裂带垂直方向上，与断裂带距离越近的点，其形变速率越大;图 11(b)、(c)则为断裂带方向上的高相干点形变速率，这些点的沉降速率均大于垂直方向上高相干点的形变速率，且在一个很小的范围内上下波动。图 9黄河三角洲主要断层分布图 10高相干点空间位置分布图 11高相干点形变速率比较33测绘通报2023 年第 6 期通过对断裂带附近高相干点形变速率的分析可知，高相干点的形变速率和其他断层的空间位置存在很大的联系，表明断裂带对黄河三角洲的地表形变有较大的影响。此外，在黄河三角洲沿海地区，如图 8 中 D5、D10、D11 区的主要地表覆盖为滩涂、盐碱地等，此类土质本身质地松软，在长期海水侵蚀的作用下也会发生地面沉降。3.2.2人为因素随着社会经济的飞速发展，海量的建筑物拔地而起。建筑物具有的质量在重力作用下会对地表产生巨大的压力，从而使一些结构不稳定的土体产生压缩凝固现象，引起地面沉降。D7、D8 区分别是东营区、广饶县、淄博区的主要城市建筑区域。由图 12 可以看出，D7 区剖线上最大沉降速率为13.245 mm/a，D8 区剖线上最大沉降速率高达24.120 mm/a。位于剖线上的特征点地表形变性质均表现为沉降，而且建筑物分布密度越大的区域，对应的沉降速率越大。图 12城市建筑区光学遥感影像及其地面沉降速率长期不合理的地下水、地下资源的开采会引起其地下含水层土壤颗粒之间的压力增大，形成土层的压缩凝固现象，引发地面沉降。D3、D4 区域属临海盐田区域，这两个区域的沉降速率明显高于其他地区。地下卤水等资源的开采所引起的地面沉降速率较大，范围为 14.52 43.64 mm/a，如图 13所示。图 13地下水、地下卤水开采区光学遥感影像及其地面沉降速率油气开采会导致地下油层压力减小，从而含油区和黏土层发生地面沉降。但石油开采时会伴随地下水回注过程，回注区域会发生轻微抬升现象。如图 14 中的 D6、D9 区为油田区域，在油田开采和地下水回注的双重作用之下，发生地表形变，主要表现为缓慢沉降，但仍有部分区域会出现轻微抬升。图 14油田区光学遥感影像及其地面沉降速率可以看出，在人类活动中，对地表形变影响最大的是地下水、地下资源的开采，油气开采次之，而建筑载荷对地表形变的影响包括建筑物质量、空间分布密度等因素。4结论随着近年来社会经济的不断发展，黄河三角洲在微观、宏观尺度上发生着日新月异的变化。通过多源数据协同的方式对黄河三角洲近年来发展与变化进行监测、分析，挖掘其中的演变规律具有重要指导意义。主要结论如下:(1)多源异构遥感数据可以用于对大面积区域的动态监测与分析。(2)本文监督分类方法的分类精度均在 82%以上，其中最高分类精度为 87.46%，说明多种指数与分类结果的协同判断可以显著提高土地利用分类的精度。(3)在宏观方面，黄河三角洲 20002021 年最主要的变化体现在其土地利用层面。在此期间黄河三角洲的 5 大地类的面积变化和空间分布均有不同程度的变化，其中，裸地是植被和建筑的主要转化来源。近 21 年间，黄河三角洲的土地利用动态度逐渐增长并趋于稳定，表明人类对土地的干扰程度在不断增加，也表明土地的利用程度不断提高。人类活动(城市扩张、土地开发、工业化等)和社会发展(经济发展、产业结构、人类意识等)是影响黄河432023 年第 6 期樊彦国，等:基于多源遥感的黄河三角洲湿地动态监测三角洲土地利用变化的决定性因素;自然环境因素对黄河三角洲土地利用的驱动作用不明显。(4)在微观层面，黄河三角洲 20002021 年发生地表形变的区域在空间分布上表现出离散、不均匀的特征。黄河三角洲地区的地表形变程度当下处于中低水平，通过对地表形变显著区的分析，断裂带、人类对地下水、地下卤水及油气自然资源的开采、建筑物载荷是引起黄河三角洲地表形变的主要因素。参考文献:1山东省人民政府办公厅关于做好民政部山东省人民政府共同推进黄河三角洲高效生态经济区民政事业发展合作协议 实施工作的通知 J 山东省人民政府公报，2011(11):41-43 2苏锐 保护黄河流域生态 助推文旅事业高质量发展 N 中国文化报，2021-10-25(001)3翟俊，金点点，陈妍，等 黄河三角洲高效生态经济区生态变化、问题与对策J 资源科学，2020，42(3):517-526 4张磊，宫兆宁，王启为，等 Sentinel-2 影像多特征优选的黄河三角洲湿地信息提取 J 遥感学报，2019，23(2):313-326 5徐振田，ALI Shahzad，张莎，等 基于 Landsat 数据的黄河三角洲湿地提取及近 30 年动态研究J 海洋湖沼通报，2020(3):70-79 6侯学会，李新华，徐洪彪，等 基于时序特征的黄河三角洲土地覆盖动态及驱动力分析 J 水资源与水工程学报，2021，32(1):233-240 7高瑞，王志勇，周晓东，等 利用多时相遥感监测与分析黄河三角洲湿地变化动态J 测绘通报，2021(4):22-27 8张金芝 基于 InSAR 时序分析技术的现代黄河三角洲地面沉降监测及典型影响因子分析D 北京:中国科学院大学，2015 9刘一霖 黄河三角洲地面沉降时序 InSAR 技术监测与地下流体开采相关性分析D 青岛:中国科学院研究生院(海洋研究所)，2016 10 陈继伟，曾琪明，焦健，等 Sentinel-1 A 卫星 TOPS 模式数据的 SBAS 时序分析方法:以黄河三角洲地区为例 J 国土资源遥感，2017，29(4):82-87 11 程霞，张永红，邓敏，等 Sentinel-1A 卫星的黄河三角洲近期地表形变分析J 测绘科学，2020，45(2):43-51 12 张金盈，崔靓，刘增珉，等 利用 Sentinel-1 SAR 数据及 SBAS 技术的大区域地表形变监测J 测绘通报，2020(7):125-129 13 王辉，曾琪明，焦健，等 结合序贯平差方法监测地表形变的 InSAR 时序分析技术J 北京大学学报(自然科学版)，2021，57(2):241-249 14 东营市统计局 东 营 市 统 计 年 鉴 2021EB-OL 2022-6-1 http:-dystjj dongying gov cn-module-download-downfile jsp?classid=0filename=99aecb89c05d4c678b9bd7ad34e9be7e pdf 15 郇红艳，谭清美，朱平 城乡一体化进程中耕地利用变化的驱动因素及区域比较J 农业工程学报，2013，29(21):201-213 16 张成扬，赵智杰 近 10 年黄河三角洲土地利用/覆盖时空变化特征与驱动因素定量分析 J 北京大学学报(自然科学版)，2015，51(1):151-158 17 孙希文 区域土地利用变化与驱动力研究D 重庆:西南大学，2005 18 屈春燕 最新 1/400 万中国活动构造空间数据库的建立 J 地震地质，2008，30(1):298-304 19邓起东 中国活动构造图M 北京:地震出版社，2007(责任编辑:马保卫)53