法律状态公告日
法律状态信息
法律状态
2019-02-12
授权
授权
2018-08-24
实质审查的生效 IPC(主分类):G06F17/50 申请日:20180110
实质审查的生效
2018-07-31
公开
公开
技术领域
本发明涉及环境监测技术,尤其涉及根据监测点-小流域-县域不同尺度土壤侵蚀监测数据之间的联系提出的一种土壤侵蚀动态监测分析方法。
背景技术
土壤侵蚀一直是全球性环境问题,土壤侵蚀造成的水土流失严重阻碍中国经济发展、破坏中国土地资源、严重影响中国生态环境。2013年公布《第一次全国水利普查水土保持情况公报》(普查基准年为2011年)结果,中国土壤侵蚀面积为294.91万Km2(未含香港、澳门特别行政区及台湾省),占中国国土面积的31.12%,但是中国的土壤侵蚀状况严峻。不同尺度的分析方法各异,目前还没有一种分析不同尺度土壤侵蚀和建立不同尺度之间联系的分析方法,因而还不能通过一种尺度的分析结果分析其余尺度的土壤侵蚀现状。
发明内容
发明目的:为了解决现有技术无法通过一种尺度的分析结果进而分析其余尺度土壤侵蚀现状的问题,本发明提供一种县域-小流域-径流小区的土壤侵蚀动态分析方法。
技术方案:本发明提供一种县域-小流域-径流小区的土壤侵蚀动态分析方法,包括以下步骤:
(1)根据县域地形图及遥感影像制作县域遥感数据土壤侵蚀分布图,分析县域土壤侵蚀模数;
(2)在县域遥感数据土壤侵蚀分布图制作基础上,截取小流域,分析小流域土壤侵蚀模数遥感数据;在小流域的出水口实地观测并记录径流携带泥沙数据,分析原始小流域土壤侵蚀模数实测数据;对比分析小流域土壤侵蚀模数遥感数据与原始小流域土壤侵蚀模数实测数据,对原始小流域土壤侵蚀模数实测数据进行修正,得到修正小流域土壤侵蚀模数实测数据;由修正小流域土壤侵蚀模数实测数据对小流域土壤侵蚀模数遥感数据进行修正;由县域包含的各修正小流域土壤侵蚀模数实测数据对县域土壤侵蚀模数遥感数据进行修正;建立县域-小流域土壤侵蚀模数修正系数表;
(3)在县域遥感数据土壤侵蚀分布图制作基础上,利用GPS实地观测记录径流小区的经纬度,分析径流小区土壤侵蚀模数遥感数据;在径流小区野外实地观测并记录监测点土壤侵蚀量数据,分析原始径流小区土壤侵蚀模数实测数据;对比分析径流小区土壤侵蚀模数遥感数据与原始径流小区土壤侵蚀模数实测数据,对原始径流小区土壤侵蚀模数实测数据进行修正,得到修正径流小区土壤侵蚀模数实测数据;由修正径流小区土壤侵蚀模数实测数据对径流小区土壤侵蚀模数遥感数据进行修正;由小流域包含的全部修正径流小区土壤侵蚀模数实测数据对小流域土壤侵蚀模数遥感数据进行修正;由县域包含的全部修正径流小区土壤侵蚀模数实测数据对县域土壤侵蚀模数遥感数据进行修正;建立县域-径流小区、小流域-径流小区土壤侵蚀模数修正系数表;
(4)县域土壤侵蚀模数遥感数据经过县域-小流域土壤侵蚀模数修正系数表修正后,依据土壤侵蚀分类分级标准SL190-2007对修正县域土壤侵蚀模数遥感数据进行分析,并制作新的县域遥感数据土壤侵蚀分布图,得出县域及小流域土壤侵蚀模数、土壤侵蚀强度等级及潜在土壤侵蚀危险区域;县域土壤侵蚀模数遥感数据经过县域-径流小区土壤侵蚀模数修正系数表修正后,依据土壤侵蚀分类分级标准SL190-2007对修正县域土壤侵蚀模数遥感数据进行分析,并制作新的县域遥感数据土壤侵蚀分布图,得出县域及小流域土壤侵蚀模数、土壤侵蚀强度等级及潜在土壤侵蚀危险区域,得出径流小区土壤侵蚀模数、土壤侵蚀强度等级;由小流域或径流小区原始土壤侵蚀模数实测数据经过修正系数表修正,依据土壤侵蚀分类分级标准SL190-2007对修正小流域实测数据或径流小区实测数据进行分析,得出县域、小流域土壤侵蚀模数、土壤侵蚀强度等级及潜在土壤侵蚀危险区域等土壤侵蚀现状,得出径流小区土壤侵蚀模数、土壤侵蚀强度等级等土壤侵蚀现状。
进一步的,步骤(2)中,采用以下方法对县域-小流域土壤侵蚀模数进行修正以及建立修正系数表:
(21)统计小流域土壤侵蚀模数遥感数据与原始小流域土壤侵蚀模数实测数据之间差值并制作原始小流域土壤侵蚀模数实测数据与差值的散点图,通过散点图确定回归方程模型,利用MATLAB软件进行曲线拟合并对回归方程模型进行误差分析,确定最佳回归方程模型,利用回归方程修正原始小流域土壤侵蚀模数实测数据,得到修正小流域土壤侵蚀模数实测数据;
(22)依据土壤侵蚀分类分级标准SL190-2007将小流域土壤侵蚀模数遥感数据与修正小流域土壤侵蚀模数实测数据按照分级指标-土壤侵蚀模数,划分微度、轻度、中度、强烈、极强烈和剧烈等六组数据,分别统计六组修正小流域土壤侵蚀模数实测数据与小流域土壤侵蚀模数遥感数据的商值,分别对六组商值做平均值处理;统计六组商值的平均值,并制作小流域土壤侵蚀分级对应商值平均值的统计表,所得统计表即为小流域土壤侵蚀模数修正系数表,所得小流域土壤侵蚀分级对应商值平均值即为小流域土壤侵蚀模数修正系数;
(23)将各小流域土壤侵蚀模数修正系数表统计合制一张土壤侵蚀模数修正系数表,即为县域-小流域土壤侵蚀模数修正系数表,表中小流域土壤侵蚀模数修正系数即为县域-小流域土壤侵蚀模数修正系数;利用Arcgis软件将各小流域土壤侵蚀模数修正遥感数据赋值替换掉各小流域土壤侵蚀模数遥感数据,制作新的县域遥感数据土壤侵蚀分布图;
进一步的,步骤(3)中,采用以下方法对县域-径流小区、小流域-径流小区土壤侵蚀模数进行修正以及建立修正系数表:(31)统计径流小区土壤侵蚀模数遥感数据与原始径流小区土壤侵蚀模数实测数据之间差值并制作原始径流小区土壤侵蚀模数实测数据与差值的散点图,通过散点图确定回归方程模型,利用MATLAB软件进行曲线拟合并对回归方程模型进行误差分析,确定最佳回归方程模型,利用回归方程修正原始径流小区土壤侵蚀模数实测数据,得到修正径流小区土壤侵蚀模数实测数据;
(32)依据土壤侵蚀分类分级标准SL190-2007将径流小区土壤侵蚀模数遥感数据与修正径流小区土壤侵蚀模数实测数据按照分级指标-土壤侵蚀模数,划分微度、轻度、中度、强烈、极强烈和剧烈等六组数据,分别统计六组修正径流小区土壤侵蚀模数实测数据与径流小区土壤侵蚀模数遥感数据的商值,分别对六组商值做平均值处理,统计六组商值的平均值,并制作径流小区土壤侵蚀分级对应商值平均值的统计表,所得统计表即为径流小区土壤侵蚀模数修正系数表,所得径流小区土壤侵蚀分级对应商值平均值即为径流小区土壤侵蚀模数修正系数;
(33)将径流小区所在小流域的土壤侵蚀模数修正系数表统计合制一张土壤侵蚀模数修正系数表,即为小流域-径流小区土壤侵蚀模数修正系数表,表中土壤侵蚀模数修正系数即为小流域-径流小区土壤侵蚀模数修正系数;
(34)将径流小区所在县域的土壤侵蚀模数修正系数表统计合制一张土壤侵蚀模数修正系数表,即为县域-径流小区土壤侵蚀模数修正系数表,表中土壤侵蚀模数修正系数即为县域-径流小区土壤侵蚀模数修正系数;利用Arcgis软件制作修正径流小区土壤侵蚀模数遥感数据分布图图层,利用Arcgis软件中图层叠加模块将DEM图层及土壤侵蚀模数修正后的分布图图层进行叠加,利用Arcgis软件中空间分析模块进行Kriging空间插值制作新的县域遥感数据土壤侵蚀分布图;
进一步的,所述步骤(1)县域土壤侵蚀分布图的制作方法为:
(11)提供县域地形图及遥感影像,根据县域地形图制作DEM图层,根据DEM图层制作降雨侵蚀力专项图、土壤可侵蚀力专项图及坡度坡长因子专项图;
(12)对遥感影像进行预处理,对预处理后的遥感影像进行土地利用现状图层制作,并提取制作归一化差分植被指数NDVI图层,制作植被覆盖因子专项图及水土保持因子专项图;
(13)将降雨侵蚀力专项图、土壤可侵蚀力专项图、坡度坡长因子专项图、植被覆盖因子专项图及水土保持因子专项图统一坐标系,通过USLE方程并将各因子专项图图层运算叠加,制作土壤侵蚀分布图,分析县域土壤侵蚀模数。
进一步的,所述小流域为最小的自然积水单元,利用Arcgis软件中Polygon Subset模块截取小流域,分析小流域土壤侵蚀模数。
进一步的,所述步骤(11)中制作DEM图层的方法为:利用R2V软件对县域地形图进行矢量数字化,利用Arcgis软件对矢量化数字图层定义地理坐标,利用Arcgis软件中3D Analyst模块制作DEM图层。
进一步的,步骤(11)中,降雨侵蚀力专项图的制作方法为:利用Arcgis软件制作气象站点分布图图层,并赋值各气象站点气象观测数据;利用Arcgis软件中图层叠加模块将DEM图层及气象站点分布图图层进行叠加;利用Arcgis软件中空间分析模块进行Kriging空间插值进行空间插值和降雨侵蚀力专项图制作;
土壤可侵蚀力专项图的制作方法为:通过土壤剖面样点的资料收集和数据采集,并对各土壤剖面样点的地形、坡度、海拔、植被、土壤信息环境条件进行坐标定位,利用Arcgis软件制作不同土壤类型坡面样点分布图,或通过中国科学院南京土壤研究所直接获取不同土壤类型坡面样点分布图;利用Arcgis软件中图层叠加模块将DEM图层及不同土壤类型剖面样点分布图图层进行叠加;利用Arcgis软件中空间分析模块进行Kriging空间插值制作土壤可侵蚀力专项图;
坡度坡长因子专项图的制作方法为:对DEM图层进行填洼处理;利用Arcgis软件中Hydrology模块及子模块Flow Accumulator对填洼处理后的DEM图层进行水流方向及汇流累计量处理,对水流方向及汇流累计量处理后的图层提取坡长并制作坡长图;利用Arcgis软件中Raster Calculator模块叠加处理坡度图和坡长图,进行坡度坡长因子专项图制作。
进一步的,所述步骤(12)中遥感影像进行预处理包括:利用Erads Imagine软件中的Image Geometric模块对遥感影像进行几何校正,利用Arcgis软件中Define Projection模块对遥感影像进行投影变换、利用Arcgis软件中Extraction模块对遥感影像进行影像拼接。
进一步的,步骤(12)中,植被覆盖因子专项图的制作方法为:利用Erads Imagine软件中的Supervised Classification模块对遥感影像进行土地利用现状图层制作;利用Arcgis软件中Spatial analyst模块提取制作归一化差分植被指数NDVI图层;利用植被指数与植被覆盖因子的数学关系对NDVI图层进行赋值的方式制作植被覆盖因子专项图;
通过资料收集,参照美国农业部水土保持措施因子表及中国水土保持措施因子值研究文献成果,汇总整理适用监测区土地利用现状的各土地利用类型水土保持措施因子值,并制作水土保持措施因子值表格;利用Erads Imagine软件中的Supervised Classification模块对遥感影像进行土地利用现状图层制作;对照水土保持措施因子值表格中的数值,对不同土地利用类型进行赋值,进行水土保持因子专项图制作。
有益效果:本发明提供的一种县域-小流域-径流小区的土壤侵蚀动态分析方法,相比较现有技术,建立不同土壤侵蚀监测尺度之间联系,既能满足大尺度数字化分析及实测结果相互校正的需求,又可以节省人力、物力。
1.本方法各监测尺度根据成熟监测方法,但又不同于传统监测方法,以土壤侵蚀模数为指标,建立县域-小流域-径流小区的土壤侵蚀模数修正系数表,可通过一个尺度的分析结果结合修正系数表修正后直接分析其余尺度土壤侵蚀现状。省去了很多人力、物力资源,操作性、实用性、数据可靠性强,避免了重复性劳动,具有较好的参考和实用价值。
2.本方法主要包括县域土壤侵蚀分析部分、流域土壤侵蚀分析部分、径流小区土壤侵蚀部分。基础在于各个尺度分析结果已经过多年数据采集,数据时间序列长,数据量大,分析时不需要重复实测、采集数据和进行实测区的重复性建设。目的在于以土壤侵蚀模数为指标,建立县域-小流域-径流小区的土壤侵蚀模数修正系数表,通过一个尺度的分析结果结合修正系数表修正后直接分析其余尺度土壤侵蚀现状。因此本监测方法落脚点在于建立县域-小流域-径流小区的土壤侵蚀模数修正系数表,仅须对修正系数表进行动态更新,就可根据需要分析出不同尺度的土壤侵蚀模数,这对于掌握省、市、县等不同行政区划或全国范围的土壤侵蚀动态监测现状来说具有重大和深远的意义。
附图说明
图1为土壤侵蚀分布图制作流程图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明。
如图1所示,县域-小流域-径流小区的土壤侵蚀动态分析方法,包括以下步骤:
第一步:根据县域地形图及遥感影像制作县域遥感数据土壤侵蚀分布图,分析县域土壤侵蚀模数。县域土壤侵蚀分布图的制作方法为:
(11)提供县域地形图及遥感影像,利用R2V软件对县域地形图进行矢量数字化,利用Arcgis软件对矢量化数字图层定义地理坐标,利用Arcgis软件中3D Analyst模块根据县域地形图制作DEM图层;
根据DEM图层,先利用Arcgis软件制作气象站点分布图图层,并赋值各气象站点气象观测数据;利用Arcgis软件中图层叠加模块将DEM图层及气象站点分布图图层进行叠加;利用Arcgis软件中空间分析模块进行Kriging空间插值进行空间插值和降雨侵蚀力(R)专项图制作;
通过土壤剖面样点的资料收集和数据采集,并对各土壤剖面样点的地形、坡度、海拔、植被、土壤信息环境条件进行坐标定位,利用Arcgis软件制作不同土壤类型坡面样点分布图,或通过中国科学院南京土壤研究所直接获取不同土壤类型坡面样点分布图;利用Arcgis软件中图层叠加模块将DEM图层及不同土壤类型剖面样点分布图图层进行叠加;利用Arcgis软件中空间分析模块进行Kriging空间插值制作土壤可侵蚀力(K)专项图;
对DEM图层进行填洼处理;利用Arcgis软件中Hydrology模块及子模块Flow Accumulator对填洼处理后的DEM图层进行水流方向及汇流累计量处理,对水流方向及汇流累计量处理后的图层提取坡长并制作坡长图;利用Arcgis软件中Raster Calculator模块叠加处理坡度图和坡长图,进行坡度坡长因子(SL)专项图制作。
(12)对遥感影像进行预处理,预处理包括:利用Erads Imagine软件中的Image Geometric模块对遥感影像进行几何校正,利用Arcgis软件中Define Projection模块对遥感影像进行投影变换、利用Arcgis软件中Extraction模块对遥感影像进行影像拼接;对预处理后的遥感影像进行土地利用现状图层制作,并提取制作归一化差分植被指数NDVI图层,制作植被覆盖因子(C)专项图及水土保持因子(P)专项图;
其中植被覆盖因子专项图的制作方法为:利用Erads Imagine软件中的Supervised Classification模块对遥感影像进行土地利用现状图层制作;利用Arcgis软件中Spatial analyst模块提取制作归一化差分植被指数NDVI图层;利用植被指数与植被覆盖因子的数学关系对NDVI图层进行赋值的方式制作植被覆盖因子专项图;
水土保持因子专项图的制作方法为:通过资料收集,参照美国农业部水土保持措施因子表及中国水土保持措施因子值研究文献成果,汇总整理适用监测区土地利用现状的各土地利用类型水土保持措施因子值,并制作水土保持措施因子值表格;利用Erads Imagine软件中的Supervised Classification模块对遥感影像进行土地利用现状图层制作;对照水土保持措施因子值表格中的数值,对不同土地利用类型进行赋值,进行水土保持因子专项图制作。
(13)将降雨侵蚀力专项图、土壤可侵蚀力专项图、坡度坡长因子专项图、植被覆盖因子专项图及水土保持因子专项图统一坐标系,通过USLE方程并将各因子专项图图层运算叠加,制作土壤侵蚀分布图,分析县域土壤侵蚀模数。
第二步:在县域遥感数据土壤侵蚀分布图制作基础上,截取小流域,分析小流域土壤侵蚀模数遥感数据,所述小流域为最小的自然积水单元,利用Arcgis软件中Polygon Subset模块截取小流域,分析小流域土壤侵蚀模数;在小流域的出水口实地观测并记录径流携带泥沙数据,分析原始小流域土壤侵蚀模数实测数据;对比分析小流域土壤侵蚀模数遥感数据与原始小流域土壤侵蚀模数实测数据,对原始小流域土壤侵蚀模数实测数据进行修正,得到修正小流域土壤侵蚀模数实测数据;由修正小流域土壤侵蚀模数实测数据对小流域土壤侵蚀模数遥感数据进行修正;由县域包含的各修正小流域土壤侵蚀模数实测数据对县域土壤侵蚀模数遥感数据进行修正;建立县域-小流域土壤侵蚀模数修正系数表;
遥感数据因卫星传感器、遥感影像处理方法等出现相应的误差,实测数据因数据采集仪器、数据采集人员采集和记录数据产生相应的误差,用实测数据修正遥感数据时不能简单的以实测数据替换遥感数据;具体采用以下方法对县域-小流域土壤侵蚀模数进行修正以及建立修正系数表:
(21)统计小流域土壤侵蚀模数遥感数据与原始小流域土壤侵蚀模数实测数据之间差值并制作原始小流域土壤侵蚀模数实测数据与差值的散点图,通过散点图确定回归方程模型,利用MATLAB软件进行曲线拟合并对回归方程模型进行误差分析,确定最佳回归方程模型,利用回归方程修正原始小流域土壤侵蚀模数实测数据,得到修正小流域土壤侵蚀模数实测数据;
(22)依据土壤侵蚀分类分级标准SL190-2007将小流域土壤侵蚀模数遥感数据与修正小流域土壤侵蚀模数实测数据按照分级指标-土壤侵蚀模数,划分微度、轻度、中度、强烈、极强烈和剧烈划六组数据,分别统计六组修正小流域土壤侵蚀模数实测数据与小流域土壤侵蚀模数遥感数据的商值,分别对六组商值做平均值处理;统计六组商值的平均值,并制作小流域土壤侵蚀分级对应商值平均值的统计表,所得统计表即为小流域土壤侵蚀模数修正系数表,所得小流域土壤侵蚀分级对应商值平均值即为小流域土壤侵蚀模数修正系数;
(23)将各小流域土壤侵蚀模数修正系数表统计合制一张土壤侵蚀模数修正系数表,即为县域-小流域土壤侵蚀模数修正系数表,表中小流域土壤侵蚀模数修正系数即为县域-小流域土壤侵蚀模数修正系数;利用Arcgis软件将各小流域土壤侵蚀模数修正遥感数据赋值替换掉各小流域土壤侵蚀模数遥感数据,制作新的县域遥感数据土壤侵蚀分布图。
第三步:在县域遥感数据土壤侵蚀分布图制作基础上,利用GPS实地观测记录径流小区的经纬度,分析径流小区土壤侵蚀模数遥感数据;在径流小区野外实地观测并记录监测点土壤侵蚀量数据,分析原始径流小区土壤侵蚀模数实测数据;对比分析径流小区土壤侵蚀模数遥感数据与原始径流小区土壤侵蚀模数实测数据,对原始径流小区土壤侵蚀模数实测数据进行修正,得到修正径流小区土壤侵蚀模数实测数据;由修正径流小区土壤侵蚀模数实测数据对径流小区土壤侵蚀模数遥感数据进行修正;由小流域包含的全部修正径流小区土壤侵蚀模数实测数据对小流域土壤侵蚀模数遥感数据进行修正;由县域包含的全部修正径流小区土壤侵蚀模数实测数据对县域土壤侵蚀模数遥感数据进行修正;建立县域-径流小区、小流域-径流小区土壤侵蚀模数修正系数表;
其中,采用以下方法对县域-径流小区、小流域-径流小区土壤侵蚀模数进行修正以及建立修正系数表:
(31)统计径流小区土壤侵蚀模数遥感数据与原始径流小区土壤侵蚀模数实测数据之间差值并制作原始径流小区土壤侵蚀模数实测数据与差值的散点图,通过散点图确定回归方程模型,利用MATLAB软件进行曲线拟合并对回归方程模型进行误差分析,确定最佳回归方程模型,利用回归方程修正原始径流小区土壤侵蚀模数实测数据,得到修正径流小区土壤侵蚀模数实测数据;
(32)依据土壤侵蚀分类分级标准SL190-2007将径流小区土壤侵蚀模数遥感数据与修正径流小区土壤侵蚀模数实测数据按照分级指标-土壤侵蚀模数,划分微度、轻度、中度、强烈、极强烈和剧烈划六组数据,分别统计六组修正径流小区土壤侵蚀模数实测数据与径流小区土壤侵蚀模数遥感数据的商值,分别对六组商值做平均值处理,统计六组商值的平均值,并制作径流小区土壤侵蚀分级对应商值平均值的统计表,所得统计表即为径流小区土壤侵蚀模数修正系数表,所得径流小区土壤侵蚀分级对应商值平均值即为径流小区土壤侵蚀模数修正系数;
(33)将径流小区所在小流域的土壤侵蚀模数修正系数表统计合制一张土壤侵蚀模数修正系数表,即为小流域-径流小区土壤侵蚀模数修正系数表,表中土壤侵蚀模数修正系数即为小流域-径流小区土壤侵蚀模数修正系数;
(34)将径流小区所在县域的土壤侵蚀模数修正系数表统计合制一张土壤侵蚀模数修正系数表,即为县域-径流小区土壤侵蚀模数修正系数表,表中土壤侵蚀模数修正系数即为县域-径流小区土壤侵蚀模数修正系数;利用Arcgis软件制作修正径流小区土壤侵蚀模数遥感数据分布图图层,利用Arcgis软件中图层叠加模块将DEM图层及土壤侵蚀模数修正后的分布图图层进行叠加,利用Arcgis软件中空间分析模块进行Kriging空间插值制作新的县域遥感数据土壤侵蚀分布图。
第四步:县域土壤侵蚀模数遥感数据经过县域-小流域土壤侵蚀模数修正系数表修正后,依据土壤侵蚀分类分级标准SL190-2007对修正县域土壤侵蚀模数遥感数据进行分析,并制作新的县域遥感数据土壤侵蚀分布图,得出县域及小流域土壤侵蚀模数、土壤侵蚀强度等级及潜在土壤侵蚀危险区域;县域土壤侵蚀模数遥感数据经过县域-径流小区土壤侵蚀模数修正系数表修正后,依据土壤侵蚀分类分级标准SL190-2007对修正县域土壤侵蚀模数遥感数据进行分析,并制作新的县域遥感数据土壤侵蚀分布图,得出县域及小流域土壤侵蚀模数、土壤侵蚀强度等级及潜在土壤侵蚀危险区域,得出径流小区土壤侵蚀模数、土壤侵蚀强度等级;由小流域或径流小区原始土壤侵蚀模数实测数据经过修正系数表修正,依据土壤侵蚀分类分级标准SL190-2007对修正小流域实测数据或径流小区实测数据进行分析,得出县域、小流域土壤侵蚀模数、土壤侵蚀强度等级及潜在土壤侵蚀危险区域等土壤侵蚀现状,得出径流小区土壤侵蚀模数、土壤侵蚀强度等级等土壤侵蚀现状。
机译: 防止流域小流域土壤侵蚀的方法
机译: 防止土壤侵蚀和径流的方法以及经过侵蚀处理的土壤。
机译: 径流分析仪和径流分析方法