6《地质灾害隐患综合遥感识别技术规程》（报批稿）.docx

ICS07.060CCSA77DZ中华人民共和国地质矿产行业标准XXTXXXXX-XXXX地质灾害隐患综合遥感识别技术规程Codeofpracticeforpotentia1.geohazardidentificationbycomprehensiveremotesensing（点击此处添加与国际标准一致性程度的标识）（报批稿）XXXX-XX-XX发布XXXX-XX-XX实施中华人民共和国自然资源部发布目次刖三I引言II1范围12规范性引用文件13术语和定义14缩略语15总则25.1 目的任务25.2 工作层次25.3 基本要求25.4 工作流程36基础资料与数据获取36.1 基础资料36.2 SAR数据36.3 光学遥感数据47地表形变特征提取47.1 地表形变InSAR数据处理方法47.2 D-InSAR数据处理57.3 StaCking-InSAR数据处理67.4 时序InSAR数据处理67.5 其他SAR/InSAR方法数据处理77.6 形变聚集区提取78地质灾害要素光学遥感解译78.1 光学遥感数据处理78.2 解译内容78.3 解译方法88.4 解译要求99地质灾害隐患综合判识与初步评价99.1 判识与评价内容99.2 综合判识99.3 风险初步评价99.4 判识与评价要求910野外查证1010.1 查证目的1010.2 查证内容1010.3 查证要求1010.4 查证方法1011成果编制1111.1 数据集编辑1111.2 图件编制1111.3 报告编写1211.4 成果提交12附录A（资料性）现有可用星载SAR传感器基本参数及应用特征表13附录B（资料性）常用InSAR数据处理方法及特点14附录C（资料性）典型地质灾害隐患变形破坏迹象遥感影像特征表15附录D（资料性）典型地质灾害遥感影像特征16附录E（资料性）地质灾害隐患综合遥感判识依据表18附录F（资料性）地质灾害隐患风险初步评价依据19附录G（规范性）地质灾害隐患综合遥感识别记录表21附录H（规范性）地质灾害隐患综合遥感识别野外查证记录表22附录I（资料性）地质灾害隐患综合遥感识别成果数据属性结构建议表23附录J（资料性）地质灾害隐患综合遥感识别成果报告提纲26参考文献27本文件按照GB/T1.12020标准化工作导则第1部分：标准化文件的结构和起草规则的规定起草。请注意本文件的某些内容可能涉及专利。本文件的发布机构不承担识别专利的责任。本文件由中华人民共和国自然资源部提出。本文件由全国自然资源与国土空间规划标准化技术委员会(SAC/TC93)归口。本文件起草单位：中国自然资源航空物探遥感中心、成都理工大学、四川省地质调查院、云南省地质环境监测院、四川省国土空间生态修复与地质灾害防治研究院、中国地质科学院岩溶地质研究所、长安大学。本文件主要起草人：王珊珊、葛大庆、刘斌、童立强、郭兆成、张玲、李为乐、梁京涛、黄成、马志刚、程洋、袁蔚林、王艳、蒋校、余琛、阎书豪、韩建云。我国地质环境与地质灾害孕灾机制复杂，传统调查方式难以有效发现隐蔽性、潜伏性地质灾害。破解复杂艰险地区重大地质灾害早期识别，动态掌控风险隐患发育状况是地质灾害调查工作的重要任务。综合应用星载InSAR地表形变监测、光学遥感解译等技术方法，以遥感观测要素与地质灾害关联特征为识别对象，在孕灾环境、致灾因子、承灾体判识分析的基础上，开展滑坡、泥石流、崩塌、地面塌陷等地质灾害隐患识别，有助于进一步摸清重大地质灾害风险隐患底数，提升地质灾害综合防治能力。重大地质灾害隐患综合遥感识别已成为地质灾害防治体系中的重要环节。为规范和指导地质灾害隐患综合遥感识别工作，制定本文件。地质灾害隐患综合遥感识别技术规程1范围本文件确立了地质灾害隐患综合遥感识别的工作程序，并规定了基础资料与数据获取、地表形变特征提取、灾害要素光学遥感解译、隐患综合判识与初步评价、野外查证、成果编制等阶段的技术要求。本文件适用于滑坡、崩塌、泥石流、地面塌陷等地质灾害隐患的综合遥感识别，其他类型地质灾害隐患遥感识别工作可参照执行。2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件；不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。GB/T39612低空数字航摄与数据处理规范CH/T30191：250001：50000光学遥感测绘卫星影像产品生产技术规范CH/T9008.3基础地理信息数字成果1：500、1：1000>1：2000数字正射影像CH/T9009.3基础地理信息数字成果1：5000、1：10000>1：25000>1：50000>1：100000数字正射影像图DZ"0261滑坡崩塌泥石流灾害调查规范（1：50000）DZ/T0284地质灾害排查规范DZ/T0438地质灾害风险调查评价规范（1:50000）3术语和定义下列术语和定义适用于本文件。3.1地质灾害隐患关联遥感特征potentia1.geohazardre1.atedremotesensingcharacteristics利用综合遥感观测手段获取的地形地貌、地表形变、变形破坏迹象、威胁对象等可直接或间接表征地质灾害形成、发育状况与潜在危害的特征。3.2地质灾害隐患综合遥感识别potentia1.geohazardidentificationbycomprehensiveremotesensing综合应用InSAR、光学遥感等对地观测技术获取地质灾害隐患关联遥感特征，结合地质灾害形成条件分析，判识地质灾害隐患的位置、类型、范围和潜在风险的过程。3.3形变聚集区deformationaccumu1.ationzone利用InSAR等技术提取的地表发生移动、变形和破坏，且与周边有明显差异的区域。3.4变形破坏迹象deformationandfai1.ure地表裂缝、岩体崩落、局部坍塌、建（构）筑物破坏等反映地质灾害孕育、发生、发展的地表异常现象和特征。4缩略语下列缩略语适用于本文件。CR：角反射器（COrnerRef1.ector）DEM：数字高程模型（Digita1.EIeVationMOde1.）D0M：数字正射影像（Digita1.C）HhophotoMap）DS-InSAR:分布式散射体干涉测量(DiStribUtedScattererInSAR)D-InSAR:合成孔径雷达差分干涉测量(Differentia1.InSAR)FFT:快速傅里叶变换(FaStFOUrierTranSfOn1.I)InSAR：合成孔径雷达干涉测量(SyntheticApertureRadarInterferometry)1.iDAR：激光雷达(1.ightDetectionandRanging)POT：像元偏移量跟踪(PiXe1.offSet-TraCking)PS-InSAR：永久散射体干涉测量(PemIanent/PersistentScattererInSAR)SAR：合成孔径雷达(SytheticApertureRadar)SBAS-InSAR:短基线集干涉测量(SmaI1.BaSeIineSUbSetS1.nSAR)S1.C：单视复数影像(Sing1.e1.ookComp1.ex)5总体要求5.1 目的任务综合利用多源遥感技术及其观测数据，开展地质灾害隐患关联遥感特征提取、解译与综合判识，获取地质灾害隐患的位置、类型、范围、形变特征和威胁对象等信息，初步评价地质灾害隐患风险等级，支撑地质灾害调查评价、监测预警、综合治理与应急防治等工作。5.2 工作层次5.2.1地质灾害隐患综合遥感识别分为一般识别与重点识别两个层次。5.2.2 一般识别按照县级行政区划为基本工作单元开展，也可按照自然地理单元开展；5.2.3 重点识别应优先选择地质灾害发育密集、地质环境条件复杂的城镇及重大工程建设规划区、人口聚集区开展。5.2.4地质环境条件复杂程度划分依据按照DZ/T0261执行。5.2.5一般识别结合工作区面积和地质灾害发育状况等确定，应采用不低于1：50(XX)比例尺；重点识别应采用不低于1：10000比例尺。5.3 基本要求5.3.1 地质灾害隐患综合遥感识别以滑坡为主，兼顾崩塌、泥石流、地面塌陷等地质灾害。5.3.2 应以SAR卫星数据、光学遥感影像、DEM为主要数据源，结合采用地形地貌、水工环地质、地质灾害、地震、气象等基础资料，通过开展地表形变特征提取、地质灾害要素光学遥感解译、地质灾害隐患综合判识与初步评价、野外查证等实现。5.3.3 3,3地表形变特征提取应利用星载InSAR相位干涉技术获取地表形变速率或形变量，确定形变地质体的形变部位、范围和强度等要素。5.3.4 地质灾害要素光学遥感解译应基于光学遥感影像与DEM数据解译获取地质灾害隐患关联的变形破坏迹象及孕灾地质环境。5.3.5 地质灾害隐患综合判识与初步评价应综合地质体的地表形变、变形破坏迹象特征及孕灾地质环境等，确认疑似地质灾害隐患的空间分布和类型，并初步评价其风险等级。5.3.6 野外查证应通过野外核查验证分析隐患识别信息的正确性，辅助识别成果修改完善。5.3.7 地质环境条件复杂区应加强无人机遥感、1.iDAR和倾斜摄影测量等技术应用，获取高精度光学遥感影像和地形数据，以提高识别效果。5.3.8 平面坐标系采用2000国家大地坐标，地图投影采用高斯-克吕格投影，1：25000及以下比例尺采用6。分带，1：10O(X)及以上比例尺采用3。分带；高程基准采用1985国家高程基准。5.3.9 应形成地质灾害隐患综合遥感识别空间数据集，编制地质灾害隐患综合遥感识别成果图件。5.4 工作流程地质灾害隐患综合遥感识别的工作流程见图1,包括基础资料与数据获取、地表形变特征提取、地质灾害要素光学遥感解译、地质灾害隐患综合判识与初步评价、野外查证和成果编制。基础资料与数据嬴)成果编制图1地质灾害隐患综合遥感识别技术流程图6基础资料与数据获取6.1基础资料6.1.1 搜集地质灾害调查、监测、防治等成果资料，尤其是工作区在册地质灾害隐患数据，用于分析工作区地质灾害空间分布状况和典型地质灾害发育特征及变化趋势。6.1.2 搜集基础地理资料，包含但不限于以下要素：县级行政区划、水系、居民点及设施、交通、管线、地貌、植被、地名及注记等。6.1.3 搜集地形资料，主要包括DEM数据、地形等高线等。若现有地形数据无法满足识别精度，应采用无人机遥感、1.iDAR、倾斜摄影测量等技术，获取高精度地形数据。6.1.4 搜集地质构造、地层岩性、气候、降雨、地震等调查成果资料，用于分析总结地质灾害发育规律和区域环境特征，基本掌握区域地质灾害的孕灾条件、诱发因素和成灾模式。6.1.5 基础资料要求现势性强，工作精度高，工作程度满足情况下资料比例尺应不低于隐患识别的比例尺。对于工作程度低的地区应收集比例尺最大或成果完成时间最新的基础资料。6.2 SAR数据6.2.1 2.1应综合工作区地形、气候以及植被、冰雪等地表覆盖条件，评估给定成像模式下的有效观测范围，优先选择干涉性能稳定，成像几何条件敏感，轨道测量精度高，基线控制满足InSAR技术要求的S1.C格式数据作为数据源。6.2.2 2.2应根据地质灾害隐患识别的工作层次、测量精度、时间长度等，结合工作区SAR数据存档情况，获取存档数据，根据需要编程定制工作周期内的与存档数据相同模式的SAR数据。当前可供选择的主要星载SAR数据源参见附录A。一般识别宜选取存档数据累积量多的数据，重点识别宜选取波段长、重复观测频率高的数据。6.2.3 2.3在数据存档或编程满足情况下应选择时间长度不短于1年，数据存档不低于20期的SAR数据。6.2.4 应在评估工作区地表形变强度、地表自然地物复杂度、地形坡度、隐患范围、时域变形特征的基础上确定SAR成像模式和参数。地表变形快的地区应选择长波长、短周期、高分辨率的SAR数据，地表变形慢的地区应选择波长较短的SAR数据。6.2.5 SAR数据波长选择时应充分考虑植被覆盖状况，植被覆盖较高的地区应优先选择长波长、穿透能力强的数据，依据植被盖程度依次选择1.波段、C波段或X波段。6.2.6 2.6SAR数据分辨率选择依据识别工作层次确定。一般识别中SAR数据分辨率应不低于20m（多视后）；重点识别中针对区域识别时SAR数据分辨率应优于IOm（多视后），针对单体隐患监测时SAR数据分辨率应优于3mo6.2.7 SAR卫星重访周期应遵循短周期原则，成像条件相近时应根据环境条件与变形特征，优先选择短周期SAR数据。6.2.8 SAR数据入射角选择应遵循SAR视线向与最大位移方向夹角最小为原则，最大程度避免因山体阴影、叠掩、透视收缩等成像畸变引起的无效观测。地形起伏大、高山峡谷地区应结合多源、多入射角、升降轨联合观测。6.2.9 SAR干涉组合应优先选择垂直基线优于1/3临界基线的像对，在数据存档满足情况下，累积数据应不低于1年，宜同时获取升轨和降轨观测数据。6.2.10 2.10同轨数据若按照单景定制，上下两景影像的重叠度应大于影像长度的10%；跨轨数据相邻两景影像重叠度应超过影像幅宽的10%o6.2.11 工作区不具备多种SAR数据可选时，应在保证获取最优性形变监测效果的前提下，适当降低卫星数据参数要求。6.3 光学遥感数据6.3.1 3.1应根据解译精度确定光学遥感数据的空间分辨率。不同解译比例尺对光学遥感数据空间分辨率的要求参照表1执行。为便于对地质灾害变形破坏迹象及相关微地貌特征进行精细解译，应优先选用空间分辨率不低于Im的光学遥感数据。表1解译比例尺对应的光学遥感数据空间分辨率列表解译比例尺空间分辨率m1：2000<0.21：5000<0.51：10000<11：50000<26.3.2 光学遥感数据时效性应不超过1年，部分地区在存档数据无法满足情况下应不超过2年，时相以覆盖雷达监测数据时间间隔为最优，应晚于具有较大影响的地震、强降雨事件发生时间。6.3.3 一般识别光学遥感数据以卫星数据为主；重点识别中若卫星数据无法满足要求时，应根据解译比例尺要求通过航空摄影测量方式获取高精度光学遥感数据。航空遥感数据采集技术要求按照GB/T39612执行。6.3.4 光学遥感数据的云、雪等覆盖率不宜大于5%,且不能覆盖重要地物。图像应符合层次丰富清晰、色调均匀、反差适中等要求，条带、噪声应尽可能少。6.3.5 针对地质灾害隐患集中发育区或者具有重大风险的地质灾害隐患，应获取多时相光学遥感数据开展对比解译，以提高解译精度。7地表形变特征提取7,1地表形变InSAR数据处理方法7.1.1 根据应用目标、工作层次与可用的SAR数据，在平衡有效测量信息、计算能力与处理效率、测量精度与覆盖完整性的基础上确定InSAR形变数据处理方法。7.1.2 一般识别中应根据SAR数据情况，选择技术复杂度低、形变信息量丰富且计算效率高的方法，在SAR数据量不足时应优先选择D-InSAR方法提取相对变形量；在SAR数据量丰富时应选择Stacking-InSAR技术提取形变速率，获取工作区地表形变分布状况。7.1.3 重点识别中应根据存档与编程数据获取情况，结合工作区环境因素和总体相干性保持能力，选择时序InSAR分析方法。区域尺度形变监测应优先选择SBAS-InSAR提取形变速率与累积形变量。单体隐患选择PS-InSAR或DS-InSAR开展长时序分析，获取相干目标时间序列。7.1.4 1,4地表大变形超出干涉相位临界条件，无法采用D-InSAR、StaCking、时序InSAR等相位干涉方法解决时，可选择高分辨率SAR影像，采用PoT方法予以求解；需开展InSAR与其他观测手段同步测量时，可布设人工角反射器，选择相位干涉或PoT方法予以处理。7.1.5 2D-InSAR数据处理7.2.1 数据预处理读入SAR卫星S1.C（1级）数据并实现数据格式转换，生成单视复数据与多视强度影像。根据卫星精密轨道数据提供情况，选择是否在S1.C影像中添加精轨数据。7.2.2 数据配准利用强度相关法进行主辅影像精确配准，条带成像模式SAR数据要求方位向和距离向误差应均小于0.1个像元，TOPS模式数据方位向频谱重叠应优于0.001个像元。SAR与DEM配准时以多视SAR强度影像为参考进行配准，根据SAR数据与DEM分辨率，对DEM进行降采样或过采样处理，并模拟成SAR影像，将模拟SAR影像与多视SAR强度影像进行空间相关计算，二者配准精度应优于0.3个像元，利用SAR坐标与DEM坐标对应查找表将DEM转换至SAR坐标系下。7.2.3 干涉相位和相干系数计算对已配准的主辅影像逐像元共辗相乘计算生成干涉相位图，依据相干系数公式并选择窗体大小逐像元计算生成相干系数图。7.2.4 差分干涉计算根据空间基线以及距离-多普勒模型和地球椭球模型计算平地相位；利用配准后DEM模拟计算地形相位。从干涉相位中去除平地和地形相位，逐像元计算生成差分干涉图，并计算生成相干系数图。7.2.5 差分干涉相位滤波利用空间自适应滤波对差分干涉图与相干系数进行空间滤波，降低噪声影响，并生成相干系数图。滤波方法根据处理范围大小与计算效率而定，一般选择GO1.dStein滤波方法。滤波窗口以5x5为基础，应根据干涉相位总体分布而扩大窗口大。局部形变区滤波选择非局域（NOn-IOCa1.）滤波方法，以保持局部形变相位细节。7.2.6 2.6相位解缠采用区域生长法或最小费用流方法对差分干涉相位进行解缠处理，得到解缠相位图。检查解缠相位图，对出现解缠错误的相位图进行重新解缠，调整参数至解缠相位连续，以保证后续处理中间结果质量。7.2.7 趋势相位去除对存在趋势相位的差分干涉图，根据干涉条纹密度分别选择FFT方法或多项式函数估计趋势相位，并从差分干涉相位或解缠相位中去除。趋势相位高频分布（条纹密集）且规则时应采用FFr方法估计空间相位分布密度，趋势相位低频分布时（条纹低于5个时）应采用多项式函数模拟。根据分布趋势选择函数项和次数，模拟趋势相位并去除，将解缠相位作为后续形变解算处理的基础。若工作区地形复杂且差分干涉图中大气相位与地形分布相关明显，应进行大气相位校正。一般采用GACOS（通用型卫星SAR大气改正，GenericAtmosphericCorrectionOn1.ineService）数据或幕-指函数进行大气相位估算，然后从差分干涉相位或其解缠结果中予以去除。幕指函数计算窗口大小应根据地形和大气相位的分布特征确定。7.2.9 视线向（1.oS）形变计算根据SAR成像几何关系将解缠相位转换至SAR视线向（1.OS）o7.2.10 地理编码利用DEM坐标系与雷达影像坐标系的转换查找表，将形变结果（形变量、形变速率或形变序列）转换到地理坐标系下。7.3Stacking-InSAR数据处理7.3.1参考影像选取分析多时相SAR数据基线分布特征，选取时间基线居中的SAR影像作为参考影像，其余SAR影像统一配准至参考影像坐标系下。7.3.2干涉像对组合利用SAR数据轨道数据计算参考影像与其余影像之间的时间基线与空间基线，生成时空基线分布图。以给定时间和空间基线阈值确定干涉相对组合方式，时空基线阈值设置应结合所用数据的相干性保持能力、工作区环境因素以及形变量值确定，一般应选择时间基线小于60天以内，空间基线优于三分之一临界基线的干涉像对生成干涉数据集。7.3.3差分干涉图集生成按照干涉像对组合进行差分干涉处理，生成差分干涉相位图集，具体按照724-728执行。7.3.4相位叠加（StaCking）处理选取解缠相位无明显跳变且连续分布的干涉相位图组合，基于最小二乘拟合估算线性形变速率，得到形变速率相位。7.3.5相位叠加结果地理编码对形变速率相位进行视线向形变计算，得到形变速率；利用DEM坐标系与SAR坐标系转换查找表，将StaCking处理所得形变速率转换到地理坐标系下。1.1.2 InSAR数据处理1.1.3 4.1时序InSAR数据预处理按照724-7.2.8生成差分干涉图集，得到时序形变差分干涉相位。7.4.2 相干目标识别相干目标包括点目标与分布式目标两类。点目标选择应采用幅度离散指数方法或相干系数阈值方法初步筛选相干目标，通过相位质量评价筛选出高相干点目标；分布式目标应依据相似性判别准则筛选出的同质像元及优化萃取后的相位，按照相位一致性准则筛选出可靠的分布式目标。7.4.3 相干目标时序分析根据先验知识给定高程改正和形变速率估计的范围，将相干像元构成DeIaUnay三角网，采用二维周期函数估计网络中任一弧段的高程误差估计值与线性速率估计值，完成相位积分并得到形变速率；利用奇异值分解方法求解非线性形变序列，将线性和非线性分量相加得到每个相干点目标的形变序列。根据大气延迟相位空间域低频和时间域高频的特征，采用时空滤波方法分离大气相位误差，一般空间滤波窗口应不大于Ikm,时域滤波窗口应不大于120天。基于相干点目标时序形变与DEM和空间基线之间的相关性，利用时域函数模型等进行DEM误差估计和校正。7.4.4 形变序列结果地理编码利用DEM坐标系与影像坐标系的转换查找表，将形变速率与形变序列转换到地理坐标系，得到时序形变数据。7.5其他SAR/InSAR方法数据处理7. 5.1POT方法对米级大变形的局部变形探测可采用分辨率优于3m的高分辨率SAR数据提取变形。选取拍摄时间早的SAR数据作为主影像，利用影像强度互相关算法精确估计主辅影像间的整体和局部偏移量，并从中扣除多项式模型拟合的轨道偏移量，得到与局部地表形变有关的偏移量残余，获得最佳变形值。8. 5.2CR-InSAR方法植被覆盖密集的地区可通过布设角反射器测量长期变形。以CR干涉相位为观测量，采用D-InSAR、时序方法等获取给定时间间隔或长时序形变序列。以CR强度峰值中心为依据检测相对位移量。9. 5.3二维形变场反演在具备SAR升、降轨观测条件下，以实际地表形变的主要方向为约束，利用升降轨1.OS向投影转换，估计水平向（东西向）、垂直向变形量，并在此基础上确定实际变形量。2. 6形变聚集区提取7. 6.1无效区域剔除综合利用SAR卫星轨道信息、DEM、SAR相干图与强度图、光学影像等数据剔除水域、阴影、叠掩等无效观测信息，生成有效观测区模板，一般保存为矢量文件。8. 6.2地表形变结果生成利用有效观测区模板对地表形变结果进行掩膜处理，得到剔除无效观测的形变结果。Stacking-InSAR形变速率结果以图像表示；时序InSAR形变结果以矢量文件表示，属性表中应包含地表形变速率，形变序列及其他过程质量控制信息。矢量文件和图像文件按照地表形变值分级设色表示。9. 6.3形变聚集区探测通过目视交互解译，或利用空间数据聚类分析方法（如阈值分割法、热点图法），或通过深度学习方法探测并识别形变聚集区。10. 6.4形变信息提取基于7.6.2生成的地表形变结果和7.6.3生成的形变聚集区文件，通过空间量算和叠加分析等，计算形变聚集区的位置、面积、形变速率等信息，并对SAR数据源、监测周期、处理方法、升/降轨道等信息进行标记。8地质灾害要素光学遥感解译11. 1光学遥感数据处理11.1.1 于光学遥感数据，通过计算机数字图像处理，提高光学遥感图像质量，制作工作区DOM,最大限度地反映地质灾害隐患及其孕灾地质环境信息。11.1.2 星光学遥感数据处理方法按照CH/T3019的相关规定执行，航空光学遥感数据处理方法按照GB/T39612执行。8.1.3制作获取的DOM成果按照比例尺应分别满足CH/T9008.3或CH/T9009.3规定的成果要求。8.2解译内容解译工作区地形地貌、地质构造、地层岩性、土地利用、人类工程活动等孕灾地质环境要素。a）地形地貌包括：1）划分工作区地貌单元，确定各地貌单元的分界线、形态和成因类型；2）陡坡、陡崖、冲沟等微地貌的位置、形态等个体特征和组合特征。b）地质构造包括：1）断层的位置、长度和延伸方向；2）褶皱的类型、长度及延伸方向；3）破碎带的性质、分布。c）地层岩性包括：D结合已有地质资料，确定地层、岩性类别及岩层产状；2）划分岩土体工程地质岩组类型。d）土地利用包括：地表水体、植被、耕地、荒坡地、城镇、交通等用地类型和分布现状。e）人类工程活动包括：工程切坡、水库库岸、露天采矿场、尾矿库、固体废物堆场等人类工程活动的类型、空间分布及其稳定性。8.2.2滑坡、崩塌、泥石流、地面塌陷等相关的变形破坏迹象解译工作区滑坡、崩塌、泥石流、地面塌陷等孕育、发生、发展有关的地表异常，包括地表裂缝、岩体崩落、局部坍塌、建（构）筑物破坏等，获取地表异常所处位置、范围、形态、长度、宽度、面积、所处部位、空间分布及其变化状况。8.2.3地质灾害解译形变聚集区和变形破坏迹象所在地质体及周边区域的地质灾害，应包含滑坡、崩塌、泥石流、地面塌陷四种类型。a）滑坡解译包括：D滑坡体所处位置、范围、形态、总体滑动方向、前后缘高程；2）滑坡体地貌部位、坡度、沟谷发育状况、植被发育状况，滑坡与重要建筑物的关系及影响程度等。b）崩塌解译包括：1）崩塌危岩体所处位置、范围、形态、总体崩塌方向、分布高程、相对高差、坡度；2）崩塌堆积体范围、面积、坡度、崩塌方向、植被类型等。c）泥石流解译包括：1）泥石流流域边界、形态、面积、主沟长度、主沟纵降比、平均坡度；2）物源区边界、水源条件、集水面积、地形坡度、岩层性质，区内植被覆盖及断裂、滑坡、崩塌、松散堆积物等不良地质现象；3）流通区边界、沟床纵横坡度、冲淤变化及泥石流痕迹，阻塞地段堆积类型，以及跌水、急弯、卡口情况等；4）堆积区边界、性质、面积，堆积扇坡降和土地覆盖等。d）地面塌陷解译包括：D地面塌陷的位置、范围、形状、面积、所在斜坡坡度；2）塌陷对地面设施的破坏程度和造成的成灾范围。8.3解译方法8.3.1 解译应在二维与三维结合的平台环境下，以交互解译为主，计算机信息提取为辅的方式进行。其中，边界勾绘和基本属性获取通过交互解译实现，位置、长度、面积、高程方向、坡度、纵降比等通过叠加DEM、DOM等数据利用空间量算和分析实现。8.3.2 孕灾地质环境解译，若前期基础地质、土地利用调查等工作程度较高，应充分利用已有资料采用补充解译和修正方式；如果基础地质、土地利用工作程度较低或无相关资料，则应开展较为详细的解译。土地利用、人类工程活动解译应采用多时相光学遥感影像对比解译相关地物的动态变化。8.3.3 变形破坏迹象解译应侧重地表异常及其变化解译。不同类型地质灾害隐患变形破坏迹象的遥感影像特征参见附录C。没有形变聚集区分布的地质体，应采用多时相光学遥感影像对比解译变形破坏迹象的变化情况，分析变形破坏迹象的稳定状态和发展趋势。8.3.4 地质灾害解译在形变聚集区和变形破坏迹象分布的地质体及周边区域开展，侧重已发生滑坡、崩塌、泥石流和地面塌陷的位置、范围、形态和发育特征的解译。不同类型地质灾害遥感影像特征参见附录D。8.4解译要求8 .4.1孕灾地质环境解译要求上图图斑面积大于4mm2的面状地物，长度大于2cm的线状地物均应解译。9 .4.2变形破坏迹象和地质灾害解译，要求按照影像原始分辨率能识别的地物均应解译。9地质灾害隐患综合判识与初步评价9-1判识与评价内容基于地表形变特征提取和地质灾害要素光学遥感解译成果，结合孕灾地质环境条件，综合分析形变聚集区或变形破坏迹象所在地质体是否具备形成崩塌、滑坡、泥石流、地面塌陷的条件和潜在威胁对象，若具备则判断为地质灾害隐患，提取其地理位置、类型、范围、运动方向、活动程度、威胁对象等信息，完成风险等级的遥感初步评价。9.2综合判识9.2.1 综合地表形变特征提取的形变聚集区，地质灾害要素光学遥感解译获取的变形破坏迹象，分析形变聚集区形变特征、变形破坏迹象发育程度、两者之间的时空关系以及孕灾地质环境和已有地质灾害发育情况，结合工作区地质灾害的形成条件，判断所在地质体是否具备形成滑坡、崩塌、泥石流、地面塌陷的条件，判断依据参考附录E执行。9.2.2 若地质体具备形成地质灾害的条件，分析其物质运移方式和影响范围，提取潜在威胁对象。潜在威胁对象包括：a）受威胁的居民点、城镇、水利水电工程、工矿企业、铁路、公路、隧道、桥梁、河流等基础设施；b）受威胁的包括耕地、园地、林地等自然资源。9.2.3 具有潜在威胁对象的地质体判断为地质灾害隐患，绘制地质灾害可能触发边界作为地质灾害隐患图斑。9.2.4 未被识别为地质灾害隐患的地质体，删除相关形变聚集区提取和地质灾害要素遥感解译成果。9. 3风险初步评价9.1.1 对综合判识确定的地质灾害隐患，利用地表形变特征、光学遥感特征和潜在威胁对象特征等进行定性评价，初步确定隐患风险等级，作为后续开展地面调查、监测、评价和风险管控的基础资料。9.1.2 风险初步评价在活动性和危害性做出定性评价的基础上，通过矩阵叠加方法实现风险等级初步划分。活动性、危害性和风险等级皆划分为“极高”“高”“中”“低”四个等级，遥感评价方法参见附录F。9.4判识与评价要求9.4.1地质灾害隐患识别成果采用面状图斑勾绘。9.4.2 应基于光学遥感数据、地表形变InSAR监测成果、地形数据和工作区在册地质灾害隐患数据等，提取地质灾害隐患的地理位置、高程、类型、是否新增、面积、运动方向、有无形变、光学遥感数据源、潜在威胁对象、风险等级初判等信息。9.4.3 应加强高位远程和链式地质灾害隐患分析，重点关注隐患源区特征、微地貌特征、地表形变和变形破坏迹象的变化发展情况，分析地质灾害的形成条件、运动过程、链生效应和可能影响范围，完成风险初步评价。9.4.4 识别过程中应填写地质灾害隐患综合遥感识别记录表，记录表格式按照附录G执行。10野外查证10.1 查证目的10.1.1 核查所识别地质灾害隐患的孕灾环境、发育特征、变形迹象、诱发因素、承灾体、发展趋势等情况。验证识别地质灾害隐患地理位置和类型的正确性。10.1.2 验证地质灾害隐患风险等级初步评价的正确性。10.1.3 修改完善地质灾害隐患识别成果。10.2 查证内容10.2.1 综合判识正确性查证核查识别地质灾害隐患所在地质体上滑坡、崩塌、泥石流、地面塌陷相关的孕灾环境、发育特征、变形迹象、威胁对象及其发展变化等情况，评价地质体上发生滑坡、崩塌、泥石流、地面塌陷等地质灾害的可能性，验证识别成果中地质灾害隐患地理位置、类型和范围是否正确。10.2.2 风险等级初步评价正确性查证分析识别地质灾害隐患所在地质体上滑坡、崩塌、泥石流、地面塌陷等的稳定性（发育程度）和发展趋势，结合威胁对象的性质和数量评价其风险等级，验证识别成果中风险等级初步评价的正确性。10.3 查证要求10.3.1 3.1野外查证点总数应不低于识别隐患数目的20%,新增地质灾害隐患识别点的查证率应不低于40%O10.3.2 针对每处查证点应填写地质灾害隐患遥感识别野外查证记录表，表格格式按照附录H执行。10.3.3 野外查证结束后应按照查证结果修改初步识别成果，使识别成果更加客观、准确。10.4 查证方法10.4.1 野外查证点选择应重点考虑以下内容：a）存疑地质灾害隐患点；b）新增地质灾害隐患点；c）与隐患库信息差异显著的已知隐患点；d）规模大或具有重大威胁的地质灾害隐患点；e）综合遥感特征突出的地质灾害隐患点。10. 4.2野外查证应携带必要的仪器工具，包括：GNSS、测距仪、卷尺、罗盘、数码照相机、摄像机、无人机摄影测量系统等。野外调查采用目视观测与现场量测相结合的方式，通过沿路观察进行记录描述，记录方式包括相机拍摄、无人机航空摄影、现场文字记录等。10.3.3 地质灾害隐患野外调查方法参照DZ/T0284执行。隐患发育特征调查中，应重视形变聚集区和变形破坏迹象解译图斑及周边地区的调查，以便与遥感识别成果对比分析。隐患风险评价按照DZ/T0438执行。10.3.4 充分重视对查证点当地群众的访问。访问内容包括查证点处是否为已发地质灾害及灾害规模、危害、可能诱发因素等信息；各类地表变形特征的发育历史，特别是区域内难以到达裂缝和垮塌的类型、规模、趋势等特征；现场调查未见明显变形特征点的灾害发育情况。10.3.5 现场照片采集应获取地质灾害隐患点整体及周边区域全景照片不少于1张，局部特征位置照片不少于3张，局部照片应详细记录拍摄位置、拍摄时间、镜头指向、变形特征等信息；照片分辨率不低于300dpi。对现场难以到达的查证点，应采用无人机航空摄影方式获取上述照片。10.3.6 隐患识别正确的判断依据应同时满足以下条件：a）野外核查确认识别的地质灾害隐患范围内具备发育地质灾害的可能，且潜在地质灾害范围大部分落在识别地质灾害隐患范围内；b）具有明确的承灾体。11成果编制11. 1数据集编辑11.1.1 地质灾害隐患综合遥感识别成果数据应以地理信息系统的空间数据集形式进行整理。11.1.2 空间数据集中应包含地表形变特征提取、地质灾害要素光学遥感解译和地质灾害隐患综合判识与初步评价的成果数据。成果数据涉及的图层（要素类）参照表2建立。表2地质灾害隐患综合遥感识别成果数据图层（要素类）列表序号编码名称图元（要素）类型约束条件来源1XBQ形变聚集区面必选地表形变特征提取2BXJX变形破坏迹象点、线、面必选光学遥感解译3YZHJ孕灾地质环境点、线、面可选光学遥感解译4DZYH地质灾害隐患面必选隐患综合判识与初步评价5YXFW潜在影响范围面可选隐患综合判识与初步评价11.1.3不同要素类（图层）的属性字段结构设置参见附录I。11.2图件编制11.2.1 应编制地质灾害隐患综合遥感识别成果图，图件内容包括：a）主图：展示识别获取的地质灾害隐患分布情况及其综合遥感特征，应包括三个层次：1）遥感底图：应采用地形阴影、高分辨率光学遥感影像、InSAR形变速率图或组合叠加。2）基础地理要素：包括村镇及以上级别地名、主要水系、重要交通干线、重要工程设施、行政区划界线等。3）地质灾害隐患要素：反映各类地质灾害隐患的位置、类型、范围、空间分布、运动方向、潜在威胁对象、风险等级初判及其形变集中分布区和变形破坏迹象等特征。b）镶图：应根