北京师范大学地理学与遥感科学学院北京师范大学地理学与遥感科学学院第一次全国水利普查技术细则第一次全国水利普查技术细则 提提 纲纲 一、水土流失普查遥感影像处理技术细则一、水土流失普查遥感影像处理技术细则二、水土流失普查专题信息提取技术细则二、水土流失普查专题信息提取技术细则三、水土流失普查强度判别与成果处理技三、水土流失普查强度判别与成果处理技术细则术细则 一、水土流失普查遥感影像处理技术细则一、水土流失普查遥感影像处理技术细则Ø 1.1.基础数据基础数据Ø 2.2.遥感数据预处理遥感数据预处理Ø 3.3.质量控制质量控制Ø 4.4.成果要求成果要求 1. 基础数据基础数据1.1 HJ-1-A 、、B数据数据 HJ-1-A和HJ-1-B卫星数据对地刈宽为700公里、地面像元分辨率为30米、4个谱段（蓝、绿、红、近红外），重访周期为4天覆盖全国需约80景HJ-1-A、HJ-1-B星CCD数据HJ-1-A卫星卫星洞洞庭庭湖湖-HJ星星数数据据 HJ-1-A和HJ-1-B卫星CCD数据用于计算植被计算植被覆盖，修正土地利用覆盖，修正土地利用。
数据获取基于中国资源卫星应用中心的数据查询系统（快视）进行， 具体要求：具体要求：（1）2009年至2010年季度数据；（2）不同侵蚀类型调查区图像的时相一致或相近；（3）图像清晰，地物层次分明，色调均一；（4）图像没有坏行、缺带，没有条带、斑点噪声和耀斑；（5）云层覆盖少的图像（即以晴空图像为优）；（6）数据覆盖全国 1.2 MODIS数据数据 MODIS数据波段范围广，有36个波段，数据空间分辨率包括250米、500米和1000米三个尺度 普查使用MODIS NDVI数据和MODIS Landcover数据，空间分辨率1km，用于计算全国植计算全国植被盖度年内变化曲线被盖度年内变化曲线，具体要求具体要求：：（1）2000年至2010年的数据，时间周期16天；（2）为实现HJ-1-A、HJ-1-B星CCD数据产品与MODIS数据 产品匹配，选取时间要尽可能一致；（3）对于不同侵蚀类型调查区图像的时相一致或相近；（4）图像没有坏行、缺带；（5）尽量挑选云层覆盖少的图像（即以晴空图像为优）；（6）数据覆盖全国 1.3 Aqua/AMSR-E数据数据 Aqua/AMSR-E亮温数据是辐射计数据，来自美国冰雪数据中心，数据在网上免费下载（https://wist.echo.nasa.gov/api/），像元大小为0.25°，用于冻融侵蚀温度相关因子反演冻融侵蚀温度相关因子反演。
具体要求：具体要求：（1）2002年至2010年数据，时间周期2天；（2）图像没有坏行、缺带；（3）数据覆盖冻融区 2.1 几何精校正几何精校正2. 遥感数据预处理遥感数据预处理2.1.1 HJ-1-A、、HJ-1-B数据数据几几何何精精校校正正流流程程图图 l(1) 数据挑选数据挑选 数据源选取需能有效的覆盖普查范围，保证数据质量和方便色调处理具体标准有： 1）相邻区域时相一致或相近的图像； 2）图幅方正，图像清晰，地物层次分明，色调均一； 3）图像没有坏行缺带，没有条带、斑点噪声和耀斑； 4）尽量挑选云层覆盖少的图像（即以晴空图像为优）； 5）为了保证图像色调均一，尽量选用同一季节的图像l(2) 几何精纠正预处理几何精纠正预处理 系统提供的标准数据产品是经辐射校正和系统几何校正的2级产品，只进行了单景数据波段间的配准、纵横向随机条纹的基本滤除和CCD影像色调的平衡归一化校正 需对获取的2级标准产品影像进行以下处理： 1）选取其中的1、2、3和4波段，采用遥感图像处理软件如ERDAS IMAGINE、ENVI等，处理采用ERDAS IMAGINE软件实现波段合成，转换成IMG格式文件。
2）对有噪声的图像进行去噪声处理以及对模糊图像进行拉伸增强处理，使图像清晰化，调整影像的亮度、对比度和色阶比，达到图像最佳效果 l(3) 图像几何精校正图像几何精校正 图像的几何精纠正以图像的几何精纠正以1:50000和和1:100000地形图或（同等或更地形图或（同等或更高）分辨率卫星影像作为参考，人工选取控制点进行几何精纠正高）分辨率卫星影像作为参考，人工选取控制点进行几何精纠正 纠正结果文件采用正轴等面积割圆锥投影，又称亚尔勃斯纠正结果文件采用正轴等面积割圆锥投影，又称亚尔勃斯（（Albers）投影，以）投影，以25°N和和47°N两条纬线为标准纬线相割；椭球两条纬线为标准纬线相割；椭球体采用国际通用的体采用国际通用的WGS-84；数据格式为；数据格式为Geotiff 几何精校正用几何精校正用ERDAS IMAGINE软件完成纠正算法采用多软件完成纠正算法采用多项式校正，每幅影像的控制点均匀分布，选取项式校正，每幅影像的控制点均匀分布，选取20个以上控制点；个以上控制点；平原地区可采用平原地区可采用2次多项式进行几何校正，山区需采用次多项式进行几何校正，山区需采用3次多项式次多项式进行几何校正。
进行几何校正 几何精校正后的影像空间采样分辨率：几何精校正后的影像空间采样分辨率：HJ-1-A和和HJ-1-B卫星卫星CCD数据产品为数据产品为30米 几何精校正产品格式：几何精校正产品格式：Geotiff；； 几何精校正的命名规则：几何精校正的命名规则：Path-Row-卫星标识卫星标识-获取日期获取日期-ref l(4) 数据地形图分幅数据地形图分幅 数据分幅工作包括镶嵌及编辑处理、地形图分幅处理、质量数据分幅工作包括镶嵌及编辑处理、地形图分幅处理、质量控制三步控制三步 1）镶嵌及编辑处理）镶嵌及编辑处理 经几何精纠正后，以较高的精度把各景影像拼接起来，形成经几何精纠正后，以较高的精度把各景影像拼接起来，形成覆盖全国的整幅影像覆盖全国的整幅影像 2）地形图分幅处理）地形图分幅处理 采用基于网格新裁切方案，具体按采用基于网格新裁切方案，具体按1:250000比例尺要求，生比例尺要求，生成覆盖某个区域（如省为单位）的网格矢量线，利用网格线对该成覆盖某个区域（如省为单位）的网格矢量线，利用网格线对该区域图像进行自动裁切，生成一系列该比例尺的影像地图。
区域图像进行自动裁切，生成一系列该比例尺的影像地图 Aqua/AMSR-E数据校正，利用头文件信息直接将地理坐标信息赋给图像，在ENVI软件完成 2.1.2 MODIS数据数据 MODIS数据校正，利用头文件信息直接将地理坐标信息赋给图像，采用ENVI软件完成2.1.3 Aqua/AMSR-E亮温数据亮温数据 2.2 大气效应纠正大气效应纠正利用利用HJ-1-A和和HJ-1-B星的蓝光星的蓝光(Alt+1)、绿光（、绿光（Alt+2）、）、红光（红光（Alt+3）和近红外（）和近红外（Alt+4）共）共4个通道的遥感影像，个通道的遥感影像，以及与其相对应的头文件（以及与其相对应的头文件（*.xml），根据灰度和），根据灰度和RGB不同不同波段组合方式浏览波段组合方式浏览HJ-1-A和和HJ-1-B星影像，并判断是否存星影像，并判断是否存在云（或云阴影）、水体和浓密植被在云（或云阴影）、水体和浓密植被2.2.1 图像信息获取图像信息获取2.2.2 图像预处理图像预处理利用公式（1）将DN值转换为表观辐亮度 （式（式 1））是原始图像像元灰度值；是原始图像像元灰度值；是绝对辐射定标系数，可从头文件或中国资源卫星应用中心网站中获取。
是绝对辐射定标系数，可从头文件或中国资源卫星应用中心网站中获取 公式（2）将表观辐亮度转换为表观反射率（式（式 2））是表观反射率；是日—地距离纠正因子；是大气外太阳光谱辐照度；是太阳天顶角上述参数中有些可从原始数据头文件获得，上述参数中有些可从原始数据头文件获得，或或网站中获取网站中获取 图像的暗目标自动提取依据比值植被指数 (RVI)，土壤调整植被指数(SAVI)和归一化水体指数 (NDWI)的综合分析法实现，各指数的计算公式如下： l（1）最小反射率法 最小反射率法利用浓密植被和水体分别在遥感影像的蓝光、近红外通道具有非常小的反射率来自动提取暗目标并获取大气参数 技术流程技术流程: 依据比值植被指数(RVI)，土壤调整植被指数(SAVI)和归一化水体指数(NDWI)的综合分析法，通过决策树方法逐步实现对清洁水体和蓝光波段浓密植被作为图像暗目标的自动提取 2.2.3 大气效应纠正大气效应纠正 利用自主开发的大气订正软件完成利用自主开发的大气订正软件完成HJ-1-A和和HJ-1-B卫星数据的大气纠正主要采用了最小反射率法、卫星数据的大气纠正主要采用了最小反射率法、MODIS辅辅助法与气象台站数据法三种方法。
助法与气象台站数据法三种方法 图像暗目标自动提取技术流程图 最最小小反反射射率法大气纠正技术流程图 操作步骤操作步骤:1)将HJ-1-A、HJ-1-B的CCD数据产品转换为大气顶的表观辐亮度和表观反射率：2)软件建立了以气溶胶光学厚度（AOD）和太阳天顶角（ ）为索引的查找表（Look-Up-Table），气溶胶类型主要是大陆乡村型3)在实现图像暗目标自动提取后，依据查找表获取气溶胶光学厚度；4)依据气溶胶和太阳天顶角，通过查找表获取其它大气参数，进行大气纠正 l（2）MODIS辅助法 技术流程技术流程: 利用已经进行大气纠正的现成遥感图像数据产品对卫星影像进行大气纠正得到地表反射率 1）采用MODIS星上定标产品得到MODIS图像获取时太阳在大气顶的辐照度，将定标好的MODIS L1B辐亮度产品MOD02转换为大气顶的表观反射率；2) 将得到的MODIS大气顶反射率和地表反射率产品MOD09配合解得大气中气溶胶光学厚度和水汽含量等参数；3) 假设MODIS与HJ过境时大气状况没有发生变化，再利用简单的大气辐射传输模型得到HJ-1星角度上的大气反射率、大气透过率，进行大气纠正,最终得到HJ-1星的地表反射率。
大大气气效效应应纠纠正正技技术术流流程程图图 操作步骤操作步骤:1) 定标好的MODIS L1B辐亮度产品MOD02转换为大气顶的表观反射率2）MODIS大气顶反射率和地表反射率产品配合解得大气参数3）HJ-1-A、HJ-1-B卫星数据大气效应纠正 l（3）气象台站数据法 当获取与HJ-1-A、HJ-1-B卫星数据相应的气象资料数据时，可以用于提出大气纠正所需的大气参数气象数据应包括：蓝光(440nm)、绿光(550nm)、红光(660nm)和近红外(880nm)四个波段任意两波段的气溶胶数据和大气柱水汽含量(单位：g/cm2)，输入的气象数据中，AOD和CWV不能小于0，或可直接输入气象数据进行大气纠正 2.3 角度效应纠正角度效应纠正 将HJ-1-A、HJ-1-B星数据方向性地表反射率产品（30m）进行角度效应纠正，得到垂直向下观测的归一化植被指数NDVI技术流程技术流程1) HJ-1-A、HJ-1-B星数据方向性地表反射率根据NDVI的定义得到带有方向性特征的NDVI数据2) HJ-1-A、HJ-1-B星数据方向性NDVI数据在HJ-1-A、HJ-1-B星30m像元为均匀植被的假设前提下使用简单的余弦纠正初步得到垂直观测的NDVI数据。
影像角度效应纠正技术流程图 操作步骤操作步骤: HJ-1-A、HJ-1-B星数据方向性地表反射率根据NDVI的定义得到带有方向性特征的NDVI数据 3.1 质量控制内容质量控制内容 质量控制是卫星数据产品质量保证的重要途径，具体内容包括HJ-1-A、HJ-1-B卫星数据产品质量检查、几何精校正产品质量精度检查、大气纠正产品质量和精度检查 利用典型地物光谱库数据或实测地物反射率数据HJ-1-A、HJ-1-B星大气纠正产品进行检查，分析并估算大气纠正产品的精度3. 质量控制质量控制 3.2 质量控制方法质量控制方法l（（1）过程抽检）过程抽检 成立项目质量检查组，在项目执行过程中，定期开展质量抽查，发现问题及时解决，严格控制把关各个作业环节的质量，保证优质高效地完成项目l（（2）阶段成果的）阶段成果的“二检二检” 在项目实施过程中，针对阶段成果，严格执行作业员自检、作业员互检的“二检”制度 作业员自检：作业人员自己进行全面的检查工作，检查比例可根据本身的作业水平决定。
作业员互检：作业员之间的交换互检，检查比例为100%，并作详细记录l（（3）最终成果预检）最终成果预检 组织专家会同项目技术组组成专门检查小组，在上述抽检和二检的基础上，根据质量检查组检查报告，对项目的最终成果进行全面的预检，填写相应的记录表 l（1）辐射精度 中等反射率地物，HJ-1-A、HJ-1-B星大气纠正产品的误差在10%左右；l（2）空间精度 HJ-1-A、HJ-1-B星数据几何精纠正精度：平原地区小于1.5个像元,山区小于3个像元 分幅接边精度：平原地区<2.5个像元，山区<3.5个像素3.3 精度要求精度要求 4.1 数据投影数据投影 成果数据投影采用正轴等面积割圆锥投影，即亚尔勃斯（Albers）投影，空间数据参数如下：4. 成果要求成果要求投影名称Albers中央经线110°投影原地纬度0°东移假定值0北移假定值0第一标准纬线25°N第二标准纬线47°N椭球体WGS 84 HJ-1-A、HJ-1-B卫星、MODIS、Aqua/AMSR-E校正数据产品格式为Geotiff格式。
依据国家基础地理信息1：25万数据标准分幅存储影像地图标注县、乡镇地名及省级公路，按县域分幅成图4.2 数据格式与分幅数据格式与分幅 二、水土流失普查专题信息提取技术细则二、水土流失普查专题信息提取技术细则Ø（一）水力侵蚀专题信息提取技术细则（一）水力侵蚀专题信息提取技术细则Ø（二）风力侵蚀专题信息提取技术细则（二）风力侵蚀专题信息提取技术细则Ø（三）冻融侵蚀专题信息提取技术细则（三）冻融侵蚀专题信息提取技术细则 （一）水力侵蚀专题信息提取技术细则（一）水力侵蚀专题信息提取技术细则l（1）降雨侵蚀力因子降雨侵蚀力因子：是指降雨导致土壤侵蚀发生的潜在能力，用一次降雨总动能E与该次降雨最大30min雨强I30的乘积EI30表示反映了雨滴对土壤颗粒的击溅分离以及降雨形成径流对土壤冲刷的综合作用l（2）降雨侵蚀力季节分布降雨侵蚀力季节分布：是指一年中某时段降雨侵蚀力占全年降雨侵蚀力的百分比，用作权重因子计算水土保持生物措施因子值l（3）降雨侵蚀力等值线图降雨侵蚀力等值线图：空间上多年平均年降雨侵蚀力相等点的连线称为等侵蚀力线，由等侵蚀力线构成的空间等值线分布图称为降雨侵蚀力等值线图，反映了多年平均降雨侵蚀力的空间变化特征。
1. 降雨侵蚀力因子降雨侵蚀力因子1.1 术语术语 l相关数据包括： 各县建立的所属气象站1981-2010年逐日雨量电子数据（.dat格式或登记表电子文档），及其说明文件；装订成册的各县“气象数据登记表”l按降雨侵蚀力因子计算公式要求，进行以下数据处理：（1）非侵蚀性降雨剔除如果日雨量小于（不含）12mm，则将该日雨量设为02）处理后的数据按原格式（.dat）重新存入相同目录下，命名为“侵蚀性降雨”1.2 数据处理数据处理 1.3 因子计算因子计算 质量控制包括两个方面： 一是一是对降雨资料进行质量控制每省最少抽查1个县2-3年的记录，核对气象数据登记表数据及其电子数据 二是二是对降雨侵蚀力因子计算结果进行抽查在水蚀区最少抽查一个县的计算结果，与利用已有分钟降水资料计算的降雨侵蚀力因子进行对比，精度>75%1.4 质量控制质量控制1.5 成果构成成果构成（1）全国各县多年平均年降雨侵蚀力因子值2）全国降雨侵蚀力等值线图，及其30m网格的栅格数据3）全国各县多年平均24个半月降雨侵蚀力占年降雨侵蚀力百分比4）全国多年24个半月降雨侵蚀力占年降雨侵蚀力百分比等值线图，及 其30m网格的栅格数据。
l数字高程模型数字高程模型：表现某高程基准下地面高程空间分布的有序数字阵列数字高程模型可以基于地形图、地面或遥感测量等方式获取的高程数据，经内插建立l坡度因子坡度因子：指CSLE的坡度因子，定义为某一坡度土壤流失量与坡度为5.13，其它条件都一致的坡面产生的土壤流失量之比率对于面上的普查，采用每个栅格的坡度因子值参与土壤流失量的计算，对于抽样调查，将取每个评价单元坡度的中位数参与土壤流失量的计算l坡长因子坡长因子：指CSLE的坡长因子，定义为某一坡面土壤流失量与坡长为22.13m、其它条件都一致的坡面产生的土壤流失量之比率用法同坡度因子值2. 坡度坡长因子坡度坡长因子2.1 术语术语 数据处理过程包括地形图数据处理、流水线和分水线提取以及数据处理过程包括地形图数据处理、流水线和分水线提取以及DEM生成l（1）地形图数据处理地形图数据处理：需要1:50000地形图、1:250000地形图和100m分辨率DEM 对国家测绘局基础地理信息中心提供的数字化地形图或自行数字化的地形图，经过必要修改编辑和质量控制，以确保每条等高线和每个高程点上标注有合理的高程值、河流由高向低流动、河流和等高线位置关系正确。
l（2）流水线和分水线提取流水线和分水线提取：利用GIS水文地貌分析功能，通过编程自动提取流水线和分水线以100m分辨率DEM为基础，划分出各大流域的2级支流划分出的流域单元，面积控制在1*104—2*104km2平原区流域划分和河流提取结果，须与比较大比例尺地形图对照，并做出必要的修编2.2 数据处理数据处理 l（3）DEM生成生成：坡度坡长因子的提取，以1:10000和1:50000数字地形图建立的DEM为数据基础其中基于1:10000地形图的DEM适用于野外调查单元对坡度坡长因子的计算，基于1:50000地形图的DEM适用于全国坡度坡长因子的提取与计算利用经过质量控制的数字地形图（包括等高线、高程点、河流、湖泊和水库等要素）和专业软件ANUDEM，经过插值生成DEM1:10000和1:50000地形图插值，在丘陵地区和山区，分辨率设置为10m和25m，平原地区和东北漫岗丘陵区分别设置为5m和10m 2.3 因子计算因子计算 以DEM为基础提取坡度坡长因子，根据坡度坡长因子的定义和坡度坡长因子值的算法（式1-11、式1-12和式1-13），利用LS计算专用程序（基于C++语言开发）完成流域坡度坡长因子值提取与计算 。
LS计算专用程序界面 地形因子的质量控制采用分层次、分区控制法首先对DEM建立、坡度和坡长的计算分别控制；其次在主要水土流失类型区各选3-5个样点，通过实测坡度和坡长，完成对质量的评价，要求精度达到75%2.4 质量控制质量控制2.5 成果构成成果构成（1）全国DEM、坡度、坡长、坡度坡长因子栅格文件，精度为30m网格2）以1:250000地形图图幅为基础，对流域坡度坡长因子值进行重新组织（拼接和切割），完成图幅坡度坡长因子值该因子值可以1:250000地形图编号来命名 l土土壤壤可可蚀蚀性性因因子子：表征土壤被冲被蚀的难易程度，反映土壤对侵蚀外营力剥蚀和搬运的敏感性，是影响土壤侵蚀的内在因素国际上常用K表示K值的大小是由土壤性质本身所决定的，在众多的土壤性质中土壤颗粒组成和有机质对K值影响最大，这也是利用土壤性质在宏观尺度上估算土壤可蚀性的理论基础理论上的获取方法，是指标准小区上单位降雨侵蚀力引起的土壤流失量，单位为thm2 h/（hm2 MJmm），K值是进行土壤侵蚀和水土流失定量评价的重要依据3. 土壤可蚀性因子土壤可蚀性因子3.1 术语术语 完成全国水蚀区土壤可蚀性K值的计算及其分布规律，一方面需要土壤数据库，包括土壤亚类和土属空间数据和属性数据；此外，为满足本次土壤侵蚀普查的精度要求，K值的计算需要精确到全国909个土属，而且计算K值的土壤属性数据必须对应到中比例尺的水蚀区土壤类型图上。
具体的数据获取和处理方法说明如下：l（（1）土壤基本属性和空间数据获取）土壤基本属性和空间数据获取 收集全国34个省市的土种志和1:500000土壤类型图，对照中国土壤发生分类系统，查阅各土种的理化性质和机械组成主要包括土种的名称、土种所在的地点、分布面积、该土种所属的亚类和土属、表层的有机质含量(%)、粗砂2-0.2mm(%)、细砂0.2-0.02mm(%)、粉砂0.02-0.002mm(%)和粘粒<0.002mm(%)将这些数据进行预处理，建立土壤理化性质的属性数据库同时，对于纸质版的土壤类型图进行几何纠正、配准和数字化，建立全国1：50万土壤类型图图形数据库，用于土壤可蚀性的计算3.2 数据处理数据处理 l（（2）试点县的野外采样过程）试点县的野外采样过程 结合土壤侵蚀野外抽样调查工作，完成单元内土壤样品的采集工作在试点省份的试点县内，将野外调查布点图与土壤分布图进行叠加分析，取该点土壤样品，并集中到省级普查中心，进行土壤理化性质 土壤样品采集方法如下：土壤样品采集方法如下：ü1）记录样点信息记录样点信息土壤样品采样点用GPS进行空间定位，记录经纬度坐标，在信息表上记录采样点基本信息，包括地点（县、市、乡、村）、土壤名称、位置坐标、土地利用类型、地貌部位（坡度、坡向）、水土保持措施、植被盖度等。
ü2）取样取样选择面积约20cm×20cm的小区域，清除土层表面枯枝落叶，挖一个小垂直断面，用铁锹垂直取1个土柱，土柱深约20cm，重量大约500-800g土样即可，装入密封袋在山区，土壤样点主要依据地形选择：3个样点一般分别分布在坡面上部、中部和底部，同时还要兼顾阴阳坡和土地利用类型在平原区，可在调查单元内随机选择3个土壤样点，要避免3个样点之间的距离太近，最后混合装袋成一个样品ü3）填写标签与封袋填写标签与封袋用铅笔填写标签（包括采集日期、地点、土壤名称、采集人等）一式两份一个放入袋中，一个贴于土样袋外侧将密封袋中的空气慢慢挤出，封好密封袋，采样完毕ü4）采样点处理采样点处理用铁锹将取样时从样点挖掘出的土壤回填到样点，可用脚踩实，以枯枝落叶覆盖样点ü5）将野外采集的土壤样品带回实验室，进行理化分析将野外采集的土壤样品带回实验室，进行理化分析，分析试验指标包括土壤有机质含量和土壤机械组成，机械组成采用美国制标准，分5个粒级，即：2-0.1mm、0.1-0.05mm、0.05-0.02mm、0.02-0.002mm和<0.002mm l（（3）共享分析测试中心的土壤理化分析数据）共享分析测试中心的土壤理化分析数据 通过共享中国科学院南京土壤研究所分析测试中心自1980年以来所有的土壤理化分析资料，计算近年来不同土壤类型的K值，用于更新部分土壤类型的K值。
主要共享全国范围内的土壤样品的采集地点、土壤类型、地理坐标、机械组成、有机质含量等数据，用于计算土壤可蚀性K值l（（4）整理已经发表的文献资料）整理已经发表的文献资料 利用CSCD文库，查阅自1980年以来所有文献上的相关土壤类型的土壤理化分析结果，以及相关的土壤可蚀性资料，用于补充和修正利用第二次土壤普查资料获取的K值 3.3 因子计算因子计算 以全国数字化土壤类型图为蓝本，以全国第二次土壤普查的样点数据（土种）为基础，先计算每个土种的土壤可蚀性K值再利用中国土壤发生分类系统中土属与土种的对应关系，将土种的K值利用面积加权平均法归并到土属上来再将土属的K值分级后链接到全国土壤类型图上来，得到全国土壤可蚀性K值分布图，具体的流程图如下： 关于土壤可蚀性关于土壤可蚀性K值的计算，在获取了上述土壤理化分析资料后，可值的计算，在获取了上述土壤理化分析资料后，可以根据理化性质的具体情况，分别选用以根据理化性质的具体情况，分别选用Williams模型和模型和Wischmeier模型l（（1））Williams模型模型 Williams等人在EPIC（Erosion—Productivity Impact Calculator）模型中发展了土壤可蚀性因子K 值的估算方法，只需要土壤有机碳和颗粒组成资料， l（（2）） Wischmeier模型模型 如果利用第二次土壤普查资料或者获取的土壤分析资料不是美国制，如果利用第二次土壤普查资料或者获取的土壤分析资料不是美国制，在计算土壤可蚀性在计算土壤可蚀性K值前就必须进行土壤颗粒分析结果的转化，即由国际值前就必须进行土壤颗粒分析结果的转化，即由国际制转化为美国制。
制转化为美国制 因为第二次土壤普查中土壤颗粒分析标准是国际制在因为第二次土壤普查中土壤颗粒分析标准是国际制在国际制向美国制的转换中（见表国际制向美国制的转换中（见表1-4），一般采用具体用图解法即在半对），一般采用具体用图解法即在半对数纸上先画出国际制的土壤颗粒大小分配曲线，然后查出美国制某一粒径数纸上先画出国际制的土壤颗粒大小分配曲线，然后查出美国制某一粒径的百分数含量的百分数含量 此方法是在计算机还没有普及的情况下进行使用，根据这个思路，我此方法是在计算机还没有普及的情况下进行使用，根据这个思路，我们进行了拟合曲线法，得到了很好的拟合效果即应用们进行了拟合曲线法，得到了很好的拟合效果即应用Y=aX+b和和Y=aX2+bX+c进行拟合进行拟合, 得出了每个土壤类型的转换方程，其中两个方程中得出了每个土壤类型的转换方程，其中两个方程中X=ln(P), P为粒径大小为粒径大小(mm)，，Y是小于是小于P粒径的累计颗粒含量百分数粒径的累计颗粒含量百分数(%) 数据获取和处理过程中肯定存在一定的误差，都可能影响到最终成果的准确性为了尽可能提高K值的计算精度，需要在计算流程的多个关键环节进行质量控制：l（1）为保证土壤属性数据的代表性和全面性，在将土种的理化性质向土属归并时，需保证该土属内的土种数占该土属内所有土种总数的85%以上。
l（2）在进行土壤样品采集时，采样点的布设需保证新采集的土壤样品数能够覆盖水蚀区试点省份80%以上的亚类l（3）在进行土壤样品分析时，标样和样品之间的误差不能大于5%l（4）在进行1:500000土壤类型图数字化和属性赋值过程中，要保证土壤图斑空间位置和属性信息的准确性大于95%l（5）在K值计算中存在国际制向美制单位转换的问题，需要借助于转换方程完成，在计算过程中要保证转换方程的确定性系数（R2）不小于90%3.4 质量控制质量控制 l（1）建立全国1:500000土壤分布图与第二次土壤普查资料的土壤属性数据之间的关系，构建全国土壤基础属性和空间信息数据库l（2）通过采样和更新土壤理化性质，获取全国909个土属的最新可蚀性K值，构建全国土壤可蚀性因子数据库，再将K值链接到土壤类型图上，进行分级后可得到全国土壤可蚀性因子图，精度为30m网格3.5 成果构成成果构成 l（1）生物措施因子生物措施因子：一定条件下耕作农地上的土壤流失量与同等条件下连续休闲对照裸地上的土壤流失量之比，为一无量纲数，其值大小变化于01之间l（2）植被盖度植被盖度：是指植被（包括叶、茎、枝）在单位面积内植被的垂直投影面积所占百分比。
l（3）归一化植被指数（归一化植被指数（NDVI））： NDVI=(NIR-R)/ (NIR+R)，常用遥感观测的近红外波段反射率（NIR）和红光波段反射率（R）计算，可用于检测植被生长状态、植被盖度和消除部分辐射误差等4. 生物措施因子生物措施因子4.1 术语术语 在获取时间序列高时空分辨率在获取时间序列高时空分辨率NDVI数据的过程中，数据的过程中，用到了以下数据：用到了以下数据：l（1）30m左右的高空间分辨率的数据：HJ-1多光谱；l（2）全国2008年1:100000土地利用图；l（3）低空间分辨率时间序列的NDVI数据：MODIS 反射率产品（MCD43B4，Nadir BRDF-Adjusted Reflectance, NBAR），空间分辨率为1km，时间分辨率为16天，时间序列5年，用于生成可靠的NDVI 时间序列； l（4）低空间分辨率的MODIS分类产品：MOD12Q1，空间分辨率为1km，时间2004年4.2 植被盖度数据处理植被盖度数据处理4.2.1 数据说明数据说明 时间序列高分辨率NDVI计算的具体流程如下图 4.2.2 时间序列时间序列NDVI计算计算 l（（1））HJ-1数据和数据和MODIS 数据的配准数据的配准 MODIS产品的投影方式与HJ-1数据的UTM投影不同，将它们投影转为一致。
l （（2））HJ-1 NDVI分布图分布图 在对HJ-1多光谱数据大气纠正和角度订正的基础上，利用NDVI的定义（如公式所示），生成不同时相HJ-1 NDVI分布图l（（3）各个不同地类）各个不同地类MODIS NDVI时间序列数据的提取时间序列数据的提取 MODIS分类产品（MOD12Q1）中包含五大分类体系，根据应用目的选取了的植被功能分类（PFT）产品，其包含地物类别分别为：常绿针叶树、常绿阔叶树、落叶针叶树、落叶阔叶树、灌木、草地、谷类作物、阔叶作物、城镇、水体等 选取多个混合像元后，通过线性模型分解，利用最小二乘的方法便可以得到亚像元上各地类的NDVIl（（4）高空间、高时间分辨率）高空间、高时间分辨率NDVI产品的生成产品的生成 融合MODIS NDVI 和HJ-1 NDVI数据得到HJ-1空间尺度上全年半月尺度的各地物类NDVI数据产品 4.2.3 植被盖度计算植被盖度计算从NDVI提取植被盖度产品的技术流程如下图 l（（1））NDVI到植被盖度转换系数确定到植被盖度转换系数确定 以MODIS NDVI产品为辅助判断标准，同时判断每一种气候地理区划里每一地类一年中MODIS NDVI的最大和最小值。
在每一种气候地理区划里找到一年中每一30m分辨率不同植被类型的直方图中最大的2-5%左右选取为NDVImax， NDVImin的选取固定不同植被类型周边小范围空间区域内裸土地类的NDVI平均值l（（2））NDVI到植被覆盖度转换到植被覆盖度转换 结合地面调查单元的植被覆盖度实测数据，通过以上公式计算得到时间分辨率15天空间分辨率30m的植被盖度l（（3）植被盖度转换系数的确定和检验）植被盖度转换系数的确定和检验 通过空间定位得到与地面调查单元测量点对应的遥感NDVI图像像元，选用部分地面实测盖度和遥感像元NDVI数据拟合NDVImax， NDVImin与k使用另外部分数据对转换系数进行精度验证 进行植被盖度季节分布监测，获取不同地区主要植被类型的植被盖度季节分布标准曲线，这是计算生物措施B因子的一个重要基础，可以对宏观遥感调查计算生物措施B因子中，对遥感反演的植被盖度进行验证l（（1）植被盖度动态监测准备）植被盖度动态监测准备 依托各地水土保持监测站的试验小流域开展工作，依据全国水土流失类型区，选择30个小流域为样点进行植被盖度动态监测（表1-5）。
4.2.4植被盖度季节分布标准曲线植被盖度季节分布标准曲线 各小流域工作人员调查收集本流域的土地利用状况，确定本地区的土地利用类型，并填写右表（表1-6） l（（2）野外选择调查地块）野外选择调查地块 要求在正式进行野外植被盖度照相测量前，确定野外地块的类型、数量和位置每个小流域（样点）按农地、人工草地、天然草地、人工林地、果园等5种类型共选择12个地块进行观测具体方法如下：ü 1、地块的选择应该代表该流域某种土地利用类型的主要生长状况，同时兼顾组织实施的方便性，如靠近交通线；ü 2、每个地块的植被覆盖均匀，面积尽量在60×60m2以上，如果达不到可适当降低标准，但地块面积应尽可能相对大一些；ü 3、各土地利用类型地块的选择和数量分配： 农地结合该地区主要耕作制度进行，选择2种代表性作物,每种作物选择2个地块； 人工草地、天然草地、人工林地、果园等每种类型各选择2个地块，果园应该选择当地主要的果树类型； 如果设定的某种植被类型在本流域没有，则可以在其它植被类型上增添地块的数量，使得地块的总数量保持控制为12个。
a)根据手持GPS，确定每个地块的位置； b)生成流域（样点）地块布设表和空间分布图流域地块布设表记录了每个地块编号、植被类型、相对位置、经纬度坐标，给出农地植被的播种时间、出苗时间和收割时间，和地块景观照片编号，并在表后附上每个地块的景观照片；空间分布图以主要交通线、流域界线（或行政界线）和标志地名作为底图,表示该样点各地块的空间分布填表时，对于农地作物的播种、出苗和收获时间，如现阶段尚不能明确，可以在以后的植被盖度动态监测的过程中补充填写 c)确定全部地块的测量顺序 l（（3）植被盖度季节分布监测实施）植被盖度季节分布监测实施 植被盖度季节分布监测从2010年开始共实施两年，每半个月对所选择的全部地块进行一次植被盖度测量，要求两天之内完成全部地块测量，相当于一个测量周期每年测量24次，每次测量要求至少2人，为保证全年测量的顺利进行，中途最好不更换测量人员 每次测量安排在每个半月的后半段（即上半月的10号～15号期间，下半月的25～31号期间）进行，采用照相法测量植被盖度，同时测量植被高度，并填写野外照相法测量植被盖度调查表。
最后及时上传照相测量记录表和数码植被图像，根据专用软件及时计算植被盖度 l（（4）植被盖度计算）植被盖度计算 采用植被盖度自动计算系统（PCOVER）自动计算植被盖度植被盖度自动计算系统（PCOVER）（计算机软件著作权登记号：2008SR12421）是专为便携型植被盖度摄影仪所开发的配套软件，运行环境为Windows 2000或Windows XP，系统安装要求约30M硬盘空间l（（5）获取植被盖度季节分布标准曲线）获取植被盖度季节分布标准曲线 完成一年植被盖度动态监测后，对计算的植被盖度信息进行整理，完成各小流域主要植被类型的植被盖度季节分布标准曲线的绘制 4.3 生物措施因子计算生物措施因子计算 以时间分辨率以时间分辨率15天、空间分辨率天、空间分辨率30m的植被盖度影像系列为基础，的植被盖度影像系列为基础，结合时间分辨率结合时间分辨率15天、空间分辨率天、空间分辨率30m的的降雨侵蚀力比例图和分辨率的的降雨侵蚀力比例图和分辨率30m的土地利用图，进行生物措施的土地利用图，进行生物措施B因子的计算具体步骤：因子的计算。
具体步骤：l（1）结合分辨率30m的土地利用图，按农、林、草等3大类型，根据时间分辨率15天、空间分辨率30m的植被盖度影像系列分别计算每半月时段的土壤流失比率SLR（水体和建筑用地忽略），得到各半月时段的土壤流失比率SLR图l（2）根据各半月时段的土壤流失比例SLR和各半月时段降雨侵蚀力占全年的比率，计算得到生物措施B因子图： 采取两种方式对基于遥感的生物措施采取两种方式对基于遥感的生物措施B因子进行质量控制因子进行质量控制l（（1）遥感反演的植被盖度季节分布进行验证和质量控制）遥感反演的植被盖度季节分布进行验证和质量控制ü全国选择30个样点，每个样点按农地、人工草地、天然草地、人工林地、果园等选择2个地块，确定每个地块的经纬度位置ü监测获取每个地块的植被盖度季节分布数据ü从时间分辨率15天、空间分辨率30m的植被盖度影像系列上，获取30个验证样点每个地块的植被盖度季节分布数据ü针对30个样点的每个地块，对实测的地表植被盖度季节分布和遥感反演的植被盖度季节分布进行对比分析，精度>75%l（（2）生物措施）生物措施B因子进行验证和质量控制因子进行验证和质量控制ü在每个水蚀类型区，统计整理该类型区的径流小区的侵蚀资料和生物措施B因子计算成果，精度>75%。
ü根据历史径流小区计算的生物措施B因子结果与基于遥感的生物措施B因子计算进行对比分析4.4 质量控制质量控制 l（1）全国30个小流域主要植被类型的植被盖度季节分布标准曲线l（2）全国半月尺度植被覆盖度图，精度为30m网格l（3）全国半月B因子图和年B因子图，精度为30m网格4.5 成果构成成果构成 l工程措施：是指通过改变小地形（如坡改梯等平整土地的措施），拦蓄地表径流，增加土壤降雨入渗，改善农业生产条件，充分利用光、温、水土资源，建立良性生态环境，减少或防止土壤侵蚀，合理开发、利用水土资源而采取的措施l耕作措施：是以保水保土保肥为主要目的，以提高农业生产为宗旨，以犁、锄、耙等为耕（整）地农具所采取的措施l地块：是指空间上连续的、具有同一土地利用类型、同一覆盖度、同一水土保持措施的区域5. 工程措施和耕作措施因子工程措施和耕作措施因子5.1 术语术语 l（（1）野外调查数据）野外调查数据 工程措施和耕作措施因子依据《野外调查单元确定及数据采集技术规定》中相关要求完成l（（2）高分辨率影像相结合的）高分辨率影像相结合的ET专题图处理专题图处理 工程措施因子E和耕作措施因子T因子主要依托土壤侵蚀野外调查单元获得的工程措施和耕作措施数据，并结合高分辨率的CBERS全色影像，采用人机交互式解译的方法获取典型区内的水土保持工程措施和耕作措施数据，最后以土地利用为依据，按照比例平均分配的方法获取。
专题图以行政县为单位，根据该县1:100000的土地利用图，在GIS软件的支持下，修改图层的属性表，添加工程措施因子E和生物措施因子T两个属性字段，将调查单元上获得的各土地利用类型的措施因子对应赋值到土地利用的属性表上，根据工程措施因子E和生物措施因子T两个属性分别生成工程措施因子E和生物措施因子T的专题图，从而实现措施因子的空间离散5.2 数据处理数据处理 ET因子获取技术流程图 5.3 因子计算因子计算l（（1）措施数量统计）措施数量统计 在野外调查单元获取的工程措施和耕作措施的数据基础上，以各县的土地利用类型为基本计算单元，统计该县内所有野外调查单元内采取了各种工程措施和耕作措施的土地总面积，并计算其占野外调查单元内该土地利用类型的比例l（（2））ET因子计算因子计算 按照工程措施、耕作措施面积比例计算工程措施因子和耕作措施因子某种土地利用类型的ET因子即为该土地利用类型中采取的每种措施的ET因子（基于野外调查单元）与该措施在该土地利用中的比例乘积的累积 ， 注：水土保持措施未涉及地区的工程措施（E）因子值尚在收集当中。
部分地区工程措施（E）措施因子值将结合调查人员收集的径流小区资料进行计算、更新 注：水土保持措施未涉及地区的耕作措施（T）因子值尚在收集当中部分地区耕作措施（T）因子值将结合调查人员收集的径流小区资料进行更新 工程措施和耕作措施因子的质量控制采用分层次、分区控制法在主要水土流失类型区各选3-5个样点，通过实地调查检验，完成对质量的评价，要求精度达到75%5.4 质量控制质量控制5.5 成果构成成果构成（1）野外调查单元内获取的工程措施和耕作措施数据（2）全国ET因子专题图，精度为30m网格 （二）风力侵蚀专题信息提取技术细则（二）风力侵蚀专题信息提取技术细则l风力因子风力因子：是指不同等级风速和不同风向对土壤风蚀发生的潜在能力，单位m/s风速和风向采用气象台站观测数据l临界侵蚀风速临界侵蚀风速：是指风力作用于地表，产生土壤风蚀的最小风速大于和等于临界侵蚀风速的各等级风速，均为侵蚀风速考虑到流沙地表的临界侵蚀风速约5m/s（在各类地表状况中的临界侵蚀风速最小），因而在统计风力因子时，自≥5m/s的风速开始统计l各等级风速累积时间分布栅格图各等级风速累积时间分布栅格图：全年大于和等于临界侵蚀风速的各等级风速的累积时间（min）在空间上的分布，反映了各等级风速累积时间的空间变化特征，采用30m分辨率的栅格数据表达各等级风速累积时间的空间变化。
1. 风力因子风力因子1.1 术语术语 在实际操作中，不同气象站和不同时段内，可能既有逐日整点风速风向数据，又有逐日4次观测的风速风向数据为了统一数据格式，便于计算，提出以下2条技术途径l（（1）若收集到的风速和风向资料为逐日）若收集到的风速和风向资料为逐日24小时整点风速风向，则按小时整点风速风向，则按1m/s间间隔直接统计全年隔直接统计全年≥5m/s的各等级风速的累积时间统计步骤如下：的各等级风速的累积时间统计步骤如下： 1）按照给定的“逐日24小时整点风速风向统计表”完成逐日24小时整点风速风向电子版本（Excel格式） 2）统计全年5.0～5.9m/s风速发生的频数，每次代表1小时，按下式换算为该等级风速的累积时间 累积时间=频数×1hr/次×60min/hr. 为便于计算，介于5.0～5.9m/s的风速等级取平均值U＝5.5m/s，6.0～6.9m/s的风速等级取平均值U＝6.5m/s，依此类推 3）按步骤（2）方法，依次统计5.0～5.9m/s、6.0～6.9m/s、7.0～7.9m/s、……各等级风速的累积时间，最高风速统计到该气象台站记录到的最大风速（Excel格式）。
1.2 数据处理数据处理 l（2）有的气象台站如果仅有每天4次（观测时间为2时、8时、14时、20时）风速风向记录，则采用线性插值的方法，获取24小时整点风速估计值，然后按1m/s间隔统计全年≥5m/s的各等级风速的累积时间统计步骤如下： 1） 按照给定的“逐日24小时整点风速风向统计表”完成逐日4次风速风向电子版本（Excel格式）逐日4次的风速风向填写在对应的时刻（2时、8时、14时、20时），待完成步骤2和步骤3后，将各整点风速插值填入对应时刻 2） 按照线性插值方法，计算相邻两个整点时刻的风速，3）完成逐日24小时整点风速风向统计表（Excel格式）4）按照“方法一”中“3）”的相同方法，统计全年各等级风速发生的频数 1.3 因子计算因子计算 考虑到风力因子在空间上的变化是渐进的，即风速和风向在空间上是逐渐变化的，相邻气象台站之间的风速和风向变化呈逐渐过渡关系因此，风力因子计算可以在GIS软件和自主开发程序的支持下实现自动运算l（1） 逐年统计1981～2010年1～5月和10～12月的逐日逐时不同风速等级和相应风向l（2）将每个月分为上半月和下半月（每月1～15日、16～30日或31日）两个时段，统计每半月时段内各风速等级的累积时间（单位：min）和16方位风向出现频率。
例如，某气象台站某一等级的风速累积时间和16方位风向出现频率为： 累积时间（min）＝频数×1hr/次×60min/hr. 式中的各风力等级按照5.0～5.9m/s取其中值5.5m/s，6.0～6.9m/s取其中值6.5m/s，7.0～7.9m/s取其中值7.5m/s，…… 依此类推：依此类推： 风向频率（％）＝风向频率（％）＝100×某风向出现频数某风向出现频数÷（半月天数（半月天数×24）） 式中的风向频率（％），按照气象台站统一的式中的风向频率（％），按照气象台站统一的16方位风向观测数方位风向观测数据计算每一风速等级对应一组据计算每一风速等级对应一组16方位风向频率数据（共方位风向频率数据（共16个数个数据），各风向频率（％）之和等于据），各风向频率（％）之和等于100％l（（3）依据）依据“（（2））”计算结果，在计算结果，在GIS软件和自主开发程序的支软件和自主开发程序的支持下，将每个气象台站作为空间上的一个点，生成持下，将每个气象台站作为空间上的一个点，生成1981～～2010年年1～～5月和月和10～～12月，每半个月的不同风力等级和相应风向频率的月，每半个月的不同风力等级和相应风向频率的30m分辨率栅格图。
分辨率栅格图 由于每半个月的不同风力等级和相应风向频率的30m分辨率栅格图是计算机自动生成的，已经有一套成熟的算法因此，风力因子质量控制的关键环节是气象数据采集、各等级风速累积时间和16方位风向出现频率计算过程 为了保证风力因子的数据质量，对质量控制的重点环节——“水土保持普查气象数据登记表（表号：P501）”、各等级风速累积时间和16方位风向出现频率统计结果，加大抽检比例其中，“水土保持普查气象数据登记表”按照30％的比例进行抽样检查，如果1张表格中出现1个或1组数据错误，记为整张表格数据有误，须重新采集数据抽检表格中的合格表格数占抽检表格总数的95％以上，方可被认为数据合格1981～2010年1～5月和10～12月期间，每半个月的不同风力等级和相应风向频率的30m分辨率栅格图，须进行100％检查1.4 质量控制质量控制 l（1）1981～2010年1～5月和10～12月期间，每半个月的不同风力等级和相应风向频率统计表l（2）1981～2010年1～5月和10～12月期间，每半个月的不同风力等级30m分辨率栅格图1.5 成果构成成果构成 l表表土土湿湿度度因因子子：是指土壤表层2.5cm深度范围内含水率（%）对土壤风蚀的潜在影响，是反映土壤对风力侵蚀的敏感程度的重要因子。
2. 表土湿度因子表土湿度因子2.1 术语术语2.2 数据处理数据处理l（1） 利用土壤介电模型（Dobson模型），参数化的地表辐射模型（Qp模型），以及辐射传输方程（ω-τ模型）模拟各种地表条件下的微波辐射，构建地表辐射模拟数据集l（2） 利用地表冻融状态判别算法进行冻融状态的判定，屏蔽冻土区域，而只对融土区进行表土含水率的反演 l（（3）） 推导微波指数（推导微波指数（MVIs）与植被透过率之间的数学关系，并利）与植被透过率之间的数学关系，并利用模拟数据建立计算植被透过率与植被含水量的半经验模型，用模拟数据建立计算植被透过率与植被含水量的半经验模型， 消除植被对表土含水率反演的影响消除植被对表土含水率反演的影响l（（4）） 利用多通道微波数据，估算地表温度，并用以求解利用多通道微波数据，估算地表温度，并用以求解10.65GHz下的裸露土壤的水平和垂直极化发射率下的裸露土壤的水平和垂直极化发射率l（（5）） 利用利用Qp参数化模型，构建能够消除地表粗糙度影响的土壤湿参数化模型，构建能够消除地表粗糙度影响的土壤湿度指数，通过模拟数据建立湿度指数与表土含水率之间的数学关系，度指数，通过模拟数据建立湿度指数与表土含水率之间的数学关系，并利用全国实测数据点拟合方程系数，建立表土含水率反演的半经并利用全国实测数据点拟合方程系数，建立表土含水率反演的半经验模型：验模型：l（（6）） 生成表土含水率数据，并进行产品比较与验证（见表土含水生成表土含水率数据，并进行产品比较与验证（见表土含水率计算模型的反演流程图）。
率计算模型的反演流程图）2011年的年的1-5月和月和10-12月，每半月计月，每半月计算一次（每月的算一次（每月的15日、日、30日或日或31日） 表土含水率计算模型的反演流程图表土含水率计算模型的反演流程图 2.3 因子计算因子计算 表土湿度因子按以下步骤计算，并获取表土湿度因子分布图表土湿度因子按以下步骤计算，并获取表土湿度因子分布图l（1）利用土壤水分计算模型反演得到的表土含水率（P）数据l（2）根据表土含水率（P），计算表土湿度因子（W）： W＝0.0932Ln（P）-0.0864l（3）利用表土湿度计算模型和湿度因子计算方法，在GIS软件支持下，绘制表土湿度因子（W）空间等值线图，等值线差值为0.52.4 质量控制质量控制 根据近年来京津风沙源治理工程区内的测量数据，对表土含水率计算模型进行校验2.5 成果构成成果构成每半月一次表土湿度因子（W）空间等值线图，等值线差值为0.5 l地表粗糙度因子：是描述植被、微地形和农田耕作技术措施对土壤风蚀影响的定量因子，单位cm3. 地表粗糙度地表粗糙度3.1 术语术语3.2 数据处理数据处理 由于影响地表粗糙度大小的因素较多，且受时间和人力资源限制，难以进行达范围的实地观测。
所以，根据植被类型、植被盖度与平均高度、地形起伏和微地貌等调查结果，通过查表获得地表粗糙度（表2-1） 3.3 因子计算因子计算 在风力侵蚀和风水交错侵蚀区，除农田呈斑块状分布，并且地表粗在风力侵蚀和风水交错侵蚀区，除农田呈斑块状分布，并且地表粗糙度在边界发生突变以外，其他诸如草地和沙地（沙漠）等的边界都呈糙度在边界发生突变以外，其他诸如草地和沙地（沙漠）等的边界都呈逐渐过渡特征因此，在地表粗糙度计算时，依据土地利用类型图，将逐渐过渡特征因此，在地表粗糙度计算时，依据土地利用类型图，将相同耕作措施的农田粗糙度作为同一粗糙度对待，单独进行计算其他相同耕作措施的农田粗糙度作为同一粗糙度对待，单独进行计算其他类型（草地、沙地、沙漠、戈壁等）之间的地表粗糙度，按照逐渐过渡类型（草地、沙地、沙漠、戈壁等）之间的地表粗糙度，按照逐渐过渡的关系，在的关系，在GIS软件和自主开发程序的支持下实现自动运算软件和自主开发程序的支持下实现自动运算l（（1）根据）根据“水土保持普查风蚀野外调查表（表号：水土保持普查风蚀野外调查表（表号：P503））”信息，信息，通过查表（表通过查表（表2-1 不同下垫面粗糙度不同下垫面粗糙度Z0值）获得地表粗糙度。
值）获得地表粗糙度l（（2）单独生成农田）单独生成农田30m分辨率地表粗糙度空间分布栅格图分辨率地表粗糙度空间分布栅格图 l（3）将每个野外调查单元作为空间上的一个点，利用查表获得的地表粗糙度数值，在GIS软件和自主开发程序的支持下，生成覆盖风力侵蚀和风水交错侵蚀区的30m分辨率空间分布栅格图l（4）在30m分辨率栅格图上剔除农田地块，即得到除农田以外的30m分辨率地表粗糙度空间分布栅格图l（5）将农田30m分辨率地表粗糙度空间分布栅格图，与除农田以外的30m分辨率地表粗糙度空间分布栅格图合并，最终生成覆盖风力侵蚀和风水交错侵蚀区的30m分辨率地表粗糙度空间分布栅格图 对最终生成的30m分辨率地表粗糙度栅格图进行抽样检查，抽检结果与野外调查单元的实地调查数据进行对比，各项数据的相对误差＜±10％，数据质量视为合格；各项数据的相对误差＞±10％，数据质量视为不合格抽样检查样点中，数据合格样点数占抽检样点总数95％以上，数据质量视为合格；数据合格样点数占抽检样点总数95％以下，数据质量视为不合格，须要重新计算和生成30m空间分布率的地表粗糙度空间分布栅格数据 3.4 质量控制质量控制3.5 成果构成成果构成风力侵蚀和风水交错侵蚀区30m分辨率地表粗糙度空间分布栅格图。
与水力侵蚀专题信息提取技术规定中的“植被盖度因子”相同在实地调查中，需要增加记录植被类型（草本群落、灌草群落、灌木群落、乔木群落）和植被平均高度（m）调查内容和格式见“水土保持普查风蚀野外调查表”（表号：P503）4. 植被盖度因子植被盖度因子 （三）冻融侵蚀专题信息提取技术细则（三）冻融侵蚀专题信息提取技术细则l冻冻融融侵侵蚀蚀区区：：侵蚀营力以冻融作用为主的区域冻融侵蚀区具有能够产生强烈冻融作用的寒冷气候条件，冻融作用是最普遍、最主要的外营力侵蚀过程，同时有相应的冻融侵蚀地貌形态表现1. 冻融侵蚀区范围界定冻融侵蚀区范围界定1.1 术语术语1.2 数据处理数据处理 冻融侵蚀区范围界定所需数据为涉及西藏、四川、云南、青海、新疆、甘肃、内蒙古、黑龙江等省、自治区1:50000DEM数据，需对以上区域DEM数据进行拼接处理，同时需要最新全国县（市、区）一级行政区划图电子版 1.3 界定方法界定方法 我国冻融侵蚀区面积广大，不同区域由于地理位置的差异，冻融侵蚀区的下限海拔高度差异很大须根据不同位置的经纬度利用冻融侵蚀分布下界海拔公式计算其冻融侵蚀区下界海拔高度。
按该式计算的海拔高度以上的区域，且剔除冰川、湖泊、沙漠等区域即为冻融侵蚀区具体计算步骤如下：（1）用DEM提取公式所需要的纬度（X1）图层和经度（X2）图层；（2）利用公式做出冻融侵蚀区下界海拔高度（H）图层；（3）把 DEM和H图层做减法运算，获得冻融侵蚀区下界海拔以上区域图层； （4）把冻融侵蚀区下界海拔以上区域图和土地利用类型图叠加，在冻融侵蚀区下界海拔以上区域图层中剔出冰川、湖泊、沙漠等分布区域；（5）适当考虑行政边界的完整性，与全国行政区划图进行叠加，对部分区域进行修正冻融侵蚀区调查技术流程 冻融侵蚀区由DEM数据和土地利用类型图数据计算获得，尽量用比例尺较大的数据，在计算前要确保所用DEM数据和土地利用类型图数据准确可靠此方法确定冻融侵蚀区范围和用遥感调查的冻融日循环天数大于180天的范围基本一致1.4 质量控制质量控制1.5 成果构成成果构成冻融侵蚀区分布图，比例尺1：100000 l冻融日循环冻融日循环：土壤温度日最大值大于0℃而日最小值小于0℃时，即认为土壤发生夜间冻结白天消融的日冻融循环过程l年冻融日循环天数年冻融日循环天数：是指一年中冻融日循环发生的天数。
2. 年冻融日循环天数年冻融日循环天数2.1 术语术语2.2 数据处理数据处理 地表的冻融日循环通过星载被动微波辐射计AMSR-E进行监测AMSR-E数据可以通过网络从美国冰雪数据中心免费下载通过选取高时间分辨率的AMSR-E辐射计数据，对冻融日循环进行监测由于AMSR-E亮温数据产品在中低纬度地区不能达到每日的全覆盖，则需要根据拼接算法进行相邻两日亮温数据的拼接 根据拼接后得到的辐射计亮温数据，计算衡量地表温度和发射率变化的判据因子，利用判别算法对地表的冻融状态进行区分，标记冻融循环发生区域，统计一年内日循环发生次数主要步骤是：l（1）提取通道36.5V亮温以衡量地表温度的变化，计算Tb18.7H同Tb36.5V的比值用以衡量地表发射率的变化；生成地表温度因子图层和地表发射率因子图层；l（2）将提取得到以上两个图层数据代入下式式进行计算，并比较计算结果的大小，若F>T，则为冻土，反之为融土；2.3 因子计算因子计算 l（3）对计算结果进行标识，获取冻融日循环发生区域，累计年冻融日循环发生天数，统计得到全国年冻融日循环天数分布图l（4）冻融日循环因子的计算流程如下： 对日冻融循环天数因子计算的数据质量进行评估，需要地面实际测量的地表状况参数对地表冻融状况的判别结果进行对比验证。
因此，在调查中，需要对土壤和地表的实际状况进行实际测量和调查，以进一步对算法进行校正及检验，提高判别精度需要观测的地表参数包括：土壤含水量、土壤温度、土壤容重、土壤质地和地表的粗糙度 AMSR-E数据的象元面积较大（约25km *25km），因此要实现对卫星算法的验证，对地面测量的观测要求，布置数据自动采集系统对地表温度和含水量进行自动采集，以保证对卫星算法精度的验证在调查区选择3个同步区域，观测区域内垂直剖面（深度>5cm）的土壤水分（单位：%）和土壤温度（单位：℃），在每个同步区域内布设100个点，每个点上分别布设一台数据采集仪、一个土壤水分传感器和两个土壤温度传感器，分别来测定土壤0~5cm含水量、0~2cm和2~5cm的土壤温度，检验精度>75%2.4 质量控制质量控制 l（1）野外测量数据集；l（2）中国区域每日地表冻融循环动态监测结果分布图，栅格精度25Km；l（3）中国区域年冻融日循环天数分布图，比例尺1：1000002.5 成果构成成果构成 l日日均均冻冻融融相相变变水水量量：是指土体冻融循环过程中发生相变的水量相变水量增加，冻结时由于水体结冰体积增大而对土体的破坏作用增加。
在冻融交替过程中，土壤含水量的增大使得土壤抗剪切力减小日均冻融相变水量反映了土壤含水量对冻融侵蚀强度的影响3. 日均冻融相变水量日均冻融相变水量3.1 术语术语3.2 数据处理数据处理 地表实际发生的相变水量可以通过AMSR-E数据的表现出的微波辐射特征变化来反映 AMSR-E数据可以通过网络从美国冰雪数据中心免费下载通过选取高时间分辨率的AMSR-E辐射计数据，对冻融日循环进行监测由于AMSR-E亮温数据产品在中低纬度地区不能达到全覆盖，要根据拼接算法进行相邻两日亮温数据的拼接 根据拼接后的AMSR-E数据，通过相变水量算法估计地表冻融循环过程中的相变水量其计算步骤为：l（1）提取通道36.5V和10.65V的亮温，计算Tb10.65V同Tb36.5V的比值用以衡量地表发射率的变化，求取升降轨数据的准发射率的差异；l（2）将提取得到的图层数据代入下式进行相变水量的计算3.3 因子计算因子计算l（3）累计年冻融相变水总量，进行平均获取日均冻融相变水量 日冻融相变水量因子计算流程 需要开展地面观测土壤水分和温度等参数来验证冻融相变水量因子，已保证其获取精度。
对于需要开展的地面测量参数布点、观测步骤与冻融日循环天数的观测规范相同，检验精度>75%3.4 质量控制质量控制3.5 成果构成成果构成（1）野外测量数据集；（2）中国区域每日地表冻融相变水量动态监测结果分布图，栅格精度25Km；（3）中国区域年均冻融相变水量分布图，比例尺1：100000 l年均降雨量：一个区域多年的年降水量的平均值4. 年均降雨量年均降雨量4.1 术语术语4.2 数据处理数据处理 收集西藏、四川、云南、青海、新疆、甘肃、内蒙古、黑龙江等省、自治区各气象站和水文站1981-2010年的降水量资料计算多年平均年降水量，绘制降水量等值线图将等值线图空间插值成为30m空间分布率的栅格数据 4.3 因子计算因子计算•（1）利用各气象和水文台站降水量资料，计算多年平均降水量•（2）绘制降水量等值线图基于GIS软件，利用各地多年平均降水量，绘制冻融侵蚀区降水量等值线图•（3）对降水量等值线图进行空间插值，建立30m分辨率冻融侵蚀区年均降水量栅格图•（4）根据冻融侵蚀年降水量因子赋值标准将年均降水量分布图分成四级，得到冻融侵蚀年均降水量因子图计算技术流程如下图： 确保收集到的降水资料全面准确，对西藏、青海等气象台站稀少的区域可以2000至2010年间TRMM 3B42数据估算的年均降水量数据进行补充。
4.4 质量控制质量控制4.5 成果构成成果构成冻融侵蚀年均降雨量因子图，比例尺1：100000 参照水力侵蚀坡度计算技术细则5. 坡度坡度 l坡坡向向：坡向为坡面法线在水平面上的投影的方向，不同坡向的地表由于接受太阳辐射强度不同，造成坡向小气候和土壤理化性质存在差异，进而影响冻融作用的强弱程度6. 坡向坡向6.1 术语术语6.2 数据处理数据处理 需要西藏、四川、云南、青海、新疆、甘肃、内蒙古、黑龙江等省、自治区1:50000地形图数据数据处理同水蚀LS因子 6.3 因子计算因子计算•（1）利用DEM，在GIS软件中生成坡向图；•（2）根据冻融侵蚀强度分级评价指标赋值标准（见下图）将坡向图分成四级，得到坡向因子图坡向因子由DEM生成，调查技术流程如下图 冻融侵蚀强度分级评价指标赋值标准 在主要冻融侵蚀类型区交通便利地方选3-5个点，通过坡向量测进行质量检验，抽样合格率要求达到75%6.4 质量控制质量控制6.5 成果构成成果构成成果冻融区坡向图，栅格精度30m l植植被被盖盖度度：是指植被（包括叶、茎、枝）在单位面积内植被的垂直投影面积所占百分比。
7. 植被盖度植被盖度7.1 7.1 术语术语术语术语7.2 7.2 因子计算因子计算因子计算因子计算 计算冻融侵蚀区植被盖度图层，计算方法同水蚀； 根据冻融侵蚀植被盖度因子赋值标准将冻融侵蚀区植被盖度图分成四级，得到冻融侵蚀植被盖度因子图7.3 7.3 质量控制质量控制质量控制质量控制7.4 7.4 成果构成成果构成成果构成成果构成同水蚀植被盖度调查一致同水蚀植被盖度调查一致冻融侵蚀植被盖度因子图冻融侵蚀植被盖度因子图 三、水土流失普查强度判别与成果处理技术细则三、水土流失普查强度判别与成果处理技术细则Ø 1.1.水力侵蚀强度判别水力侵蚀强度判别Ø 2.2.风力侵蚀强度判别风力侵蚀强度判别Ø 3.3.冻融侵蚀强度判别冻融侵蚀强度判别Ø 4.4.成果处理技术规定成果处理技术规定 1.1 模型说明模型说明1. 水力侵蚀强度判别水力侵蚀强度判别 计算分两个层次进行：计算分两个层次进行： 一是按一定抽样密度，在全国分层抽样布设野外调查单元，计算以野外调查单元为基础的土壤水蚀模数 二是在全国划分30m30m网格，计算以该网格为基础的土壤水蚀模数。
在此基础上，基于野外调查单元计算结果对网格计算的相关参数与结果进行检验和校正 1.2 计算方法计算方法 野外调查单元布设及数据采集详见《野外调查单元的确定及数据采集技术细则》 水力侵蚀力模型中各个因子计算详见《土流失专题信息提取技术细则》中的“水力侵蚀专题信息提取技术细则”部分 水力侵蚀计算过程水力侵蚀计算过程 将计算得到的全国水蚀模数按照土壤侵蚀分类分级标准（SL190-2007，表1）进行分级，评价土壤水蚀强度 1.3 判别标准判别标准注：不同侵蚀类型区采用不同的容许土壤流失量，东北黑土区、北方土注：不同侵蚀类型区采用不同的容许土壤流失量，东北黑土区、北方土石山区容许土壤流失量为石山区容许土壤流失量为200 t/(km2·a)，南方红壤丘陵区、西南土石山，南方红壤丘陵区、西南土石山区容许土壤流失量为区容许土壤流失量为500 t/(km2·a)，西北黄土高原区容许土壤流失量为，西北黄土高原区容许土壤流失量为1000 t/(km2·a) 质量控制分两方面进行：一是通过抽样，比较土壤水蚀模数计算结果与小区观测结果的差异；二是进行专家评估。
l与小区观测结果比较：在每个水蚀类型区，最少抽取一个水土保持试验站不同土地利用或不同水土保持措施下的小区观测资料，与其所对应地区相同土地利用或水土保持措施下的计算结果进行比较l专家评估：在每个水蚀类型区，最少抽取一个省，给出不同空间分辨率，如30m30m、5km5km、县级等土壤水蚀模数分布及土壤水蚀强度分布图通过专家经验打分，对计算结果进行评估，评估精度>75% 1.4 质量控制质量控制 2.1 模型说明模型说明2. 风力侵蚀强度判别风力侵蚀强度判别风力侵蚀量计算采用国内经验模型 l（（1）计算时段）计算时段 风力侵蚀采取每半月计算一次，全年逐月累加的方法，获得风力侵和风水交错侵蚀区的风力侵蚀强度数据采取风力侵蚀和水力侵蚀累加的方法，获得风水交错侵蚀区的侵蚀强度数据考虑到每年6～9月基本没有风蚀发生，全年风力侵蚀计算时段为1～5月和10～12月，共计算16次l（（2）计算方法）计算方法 在获得植被盖度空间分布栅格数据、临界侵蚀风速、风力因子空间分布栅格数据、表土湿度因子空间等值线图、地表粗糙度空间分布栅格数据之后，将每个栅格数据（等值线）图作为一个图层，在ENVI4.5+IDL和自主开发程序wind_erosion.pro的支持下，实现计算机自动运算。
运算流程如下图 2.2 计算方法计算方法 风力侵蚀运算流程风力侵蚀运算流程 风力侵蚀强度依据风力侵蚀强度依据“土壤侵蚀分类分级标准土壤侵蚀分类分级标准（（SL190-2007））”判别（表判别（表3），分为微度、轻度、中度、），分为微度、轻度、中度、强烈、极强烈、剧烈共六级强烈、极强烈、剧烈共六级 2.3 判别标准判别标准表3 风力侵蚀强度分级标准（SL190-2007）级别平均侵蚀模数[t/(km2·a)]微度<200轻度200~2500中度2500~5000强烈5000~8000极强烈8000~15000剧烈>15000 风力侵蚀强度成果的质量控制，采取计算结果与多年平均风力侵蚀强度和野外实地测量风力侵蚀强度相对比的方法，检查计算结果是否准确具体技术途径如下：具体技术途径如下：ü计算结果与多年平均风力侵蚀强度的对比计算结果与多年平均风力侵蚀强度的对比首先，在农牧交错带和典型草原区分别选择一个县（市、旗），利用Cs137元素示踪的方法，测定近50年来的多年平均风力侵蚀强度然后，在1981～2010年期间，挑选其中降水量和风力接近于多年平均状况的年份，用本次调查使用的模型进行计算，获得模型计算出的风力侵蚀强度。
最后，将模型计算的风力侵蚀强度与Cs137元素示踪方法测定的多年平均风力侵蚀强度进行比对，检查两者的差异如果两者接近，则证明模型计算结果是可靠的；相反，则应重新调整模型参数，直至两者接近为止ü计算结果与野外实地测量风力侵蚀强度的对比计算结果与野外实地测量风力侵蚀强度的对比根据多年积累的不同土地利用类型风力侵蚀强度实地观测数据，对比实地观测点风力侵蚀强度和对应点模型计算风力侵蚀强度值，检查两者的差异如果两者接近，则证明模型计算结果是可靠的；相反，则应重新调整模型参数，直至两者接近为止2.4 质量控制质量控制 冻融侵蚀评价采用指标评价法选择年冻融日循环天数、日均冻融相变水量、年均降水量、坡度、坡向和植被盖度作为评价指标，对各指标赋予权重，计算综合评价指数，按照分级评价标准判定对冻融侵蚀强度 3. 冻融侵蚀强度判别冻融侵蚀强度判别3.1 模型说明模型说明 l根据冻融侵蚀强度分级评价指标赋值标准获得冻融侵蚀评价各因子图层；l将各因子图层代入模型计算综合评价指数图层；l对冻融侵蚀综合评价指数图层按照冻融侵蚀判别标准分级得到冻融侵蚀强度分布图3.2 计算方法计算方法3.3 判别标准判别标准冻融侵蚀强度判别标准见表 l（1）遥感反演的的冻融相变水量数据，抽样率1%，抽样合格率要求达到80%以上；l（2）“冻融侵蚀野外调查表”数据，抽样率2%，抽样合格率要求达到100%以上；l（3）整个调查抽取2%野外调查单元进行验证，合格率要求达到100% 。
3.4 质量控制质量控制 l（（1）区域内部数据处理）区域内部数据处理 以县级行政单元为区内数据平衡基本单元，以其面积为标准，进行区域内部数据平衡l（（2）区域之间数据处理）区域之间数据处理 以省级行政单元为区域数据平衡基本单元，以其面积为标准，进行区域之间数据平衡全国面积为总面积控制标准4. 成果处理技术规定成果处理技术规定4.1 数据平衡数据平衡4.2 成果汇总成果汇总各省、县土壤侵蚀强度分级面积统计 l（（1）技术依据）技术依据 《SL73.8-2001水利水电工程制图标准水土保持图》 《GB/T 20257.3-2006 国家基本比例尺地图图式第三部分：1：25000 1：50000 1：100000地形图图式》； 《GB/T 20257.4-2006国家基本比例尺地图图式第四部分：1：250000 1：500000 1：1000000地形图图式》； 《CH/T 4015-2001地图符号库建立的基本规定》； 《GB/T 18317-2001专题地图信息分类与代码》； 《GB/T 13923-2006基础地理信息要素分类与代码》。
l（（2）精度要求）精度要求 专题图精度与专题数据保持一致，虽然挂图和图集涉及到多种比例尺数据要求，但由于此次土壤侵蚀普查数据库采用统一比例尺成果，因此，考虑到工作量等原因，对于不同比例尺专题图编制，不对数据作缩编处理 基础地理要素则应根据不同比例尺图面表达的要求，进行综合取舍 4.3 普查图件编制普查图件编制 l（（3）图件编制说明）图件编制说明1）图面尺寸及比例尺设计）图面尺寸及比例尺设计 根据不同的成图范围，设计不同的制图比例尺，编制国家级、流域级、省（市）级和县级专题挂图图面尺寸采用宽不超过2米，高不超过1.8米为宜 全国专题挂图全国专题挂图比例尺采用1：400比例尺； 流域级专题挂图流域级专题挂图根据流域面积大小，采用适当的比例尺； 省级专题挂图省级专题挂图根据面积大小，可采用1:250000-1:500000比例尺； 市级专题挂图市级专题挂图根据面积大小，可采用1:50000-1:100000比例尺； 县级专题挂图县级专题挂图根据面积大小，可采用1:50000比例尺左右比例尺 2）图集）图集 以省为单位编制专题图集，每个省的图集内容包括各个县的12项专题图，图集印刷尺寸：16K，各县根据面积大小，采用适当的比例尺。
3）地图表示内容）地图表示内容 基础地理要素：经纬度坐标及网线、水系及其注记等，根据不同比例尺图面表达要求，应进行综合取舍；4）图件编制要求）图件编制要求 图外要素齐全，包括图名、数学基础说明、成图比例尺、制图单位全称、图例等 图内、外要素的颜色、图案、线型等表示符合要求，图件整饰规范、美观 内容的选取和表示要层次分明，符号、注记等正确，清晰易读 汇报完毕！汇报完毕！谢谢！谢谢！ 。