1.一种面向黄土填土沉降区的边界识别与提取方法，其特征在于，包括以下步骤：S1、获取原始SAR影像数据并进行预处理，使用时间序列InSAR技术获取影像覆盖范围内黄土填土区的地表不均匀沉降数据；

S2、根据获得的地面沉降空间分布结果将沉降区划分为m个区域，其中，m为沉降漏斗个数；再按照划分的区域导出沉降区的沉降点的基本属性表，并以此建立地面沉降点数据库，将沉降区域进行分类管理；

S3、在填土沉降区的基础地表形变数据库基础之上，统计出每个沉降区域内沉降速率相同或相近的沉降观测点，将不同区域内所有等值点通过连接形成等值线；

S4、设定沉降速率和累计沉降量阈值，将大于一定阈值的沉降点视为位于填土区内的沉降点，相应的沉降点构成的等值线即可视为该区域内的填土区边界。

2.根据权利要求1所述的一种面向黄土填土沉降区的边界识别与提取方法，其特征在于，步骤S1包括：

在影像预处理阶段，首先利用外部精密轨道数据完成轨道校正；其次选定主影像后利用外部数字高程模型数据完成地形校正；再将副影像配准到主影像几何，通过干涉生成干涉图。

3.根据权利要求2所述的一种面向黄土填土沉降区的边界识别与提取方法，其特征在于，步骤S1还包括：

获得影像覆盖范围内的黄土填土区地面沉降分布情况：首先设定时空基线阈值将预处理获得的干涉图构成干涉对基线网络，然后设定时空相干性阈值筛选相干点目标，以去除相干性较差的点；最后通过相位解缠得到干涉相位时间序列结果，通过网络反演得到研究区内的相应沉降结果信息。

4.根据权利要求1所述的一种面向黄土填土沉降区的边界识别与提取方法，其特征在于，基本属性表包括形变速率表和累计形变量表；形变速率表中包含沉降观测点的坐标、对应的形变速率信息；累计形变量表中包含沉降观测点的坐标、SAR影像时间序列中各个时间段的累计形变量信息。

5.根据权利要求4所述的一种面向黄土填土沉降区的边界识别与提取方法，其特征在于，每个时间段的沉降观测点都被赋予唯一的ID进行存储。