CN108021752A - 一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置 - Google Patents
一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置 Download PDFInfo
- Publication number
- CN108021752A CN108021752A CN201711268476.3A CN201711268476A CN108021752A CN 108021752 A CN108021752 A CN 108021752A CN 201711268476 A CN201711268476 A CN 201711268476A CN 108021752 A CN108021752 A CN 108021752A
- Authority
- CN
- China
- Prior art keywords
- soil
- erosion
- modulus
- plot
- measured
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
Abstract
本发明涉及水土保持与水土流失预测领域,具体涉及一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置,是为了解决现有技术的水土流失强度的预测方法计算繁琐,会降低计算程序的执行效率,并且会加大基础资料的搜集难度的缺点而提出的,具体包括:在地理信息平台Arcgis中获取待测算区域的侵蚀图斑背景要素;背景要素包括面积、坡度、植被、土地利用以及土壤;根据所述侵蚀图斑背景要素,使用土壤流失方程对每一个地块计算待测算区域的土壤侵蚀模数;对步骤二中计算得到的各个地块的土壤侵蚀模数进行加权平均,得到全流域土壤侵蚀模数,并将所述全流域土壤侵蚀模数作为待测算区域的土壤侵蚀程度。本发明适用于水土流失预测。
Description
技术领域
本发明涉及水土保持与水土流失预测领域,具体涉及一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置。
背景技术
土壤侵蚀模数是直观反映水土流失程度强弱的一个非常重要的指标,也是水土保持措施设计标准的重要依据。由于土壤侵蚀模数影响因子的复杂性以及由此导致的动态可变性,使得探究各种土壤侵蚀类型的侵蚀模数一直成为水土保持研究者孜孜以求的目标。尤其是近年来,随着开发建设项目水土保持方案编制工作的全面展开,对开发建设项目建设过程中产生的土壤侵蚀量的推算是水保持方案编制中的一项非常重要的基础性工作和必不可少的内容。目前,对土壤侵蚀模数的推算方法有很多种,其中使用土壤流失方程进行计算的方法中,方程的条件因子的获得非常繁琐,计算量很大,在实际的设计工作中很难实现,严重影响了工作效率。因为一方面认为方程考虑到的条件因子越多、不同地貌不同气候降水区域的取值差异越大,则得到的结果越准确,而这种复杂的计算往往会明显增加计算程序的执行时间,也会加大基础资料的搜集难度(即公式中出现的因子都应该对其进行基础资料的搜集整理)。
发明内容
本发明的目的是为了解决现有技术的水土流失强度的计算方法计算繁琐,会降低计算程序的执行效率,并且会加大基础资料的搜集难度的缺点,而提出一种应用于黑土区小流域侵蚀模数的计算方法,包括:
步骤一、在地理信息平台Arcgis中获取待测算区域的侵蚀图斑背景要素;背景要素包括面积、坡度、植被、土地利用类型、土壤种类等。
步骤二、根据所述侵蚀图斑背景要素,使用土壤流失方程对流域范围内各种土地利用类型的土壤侵蚀模数进行推导。
步骤三、对步骤二中计算得到的各个地类的土壤侵蚀模数进行加权平均,得到全流域土壤侵蚀模数,并将所述全流域土壤侵蚀模数作为待测算区域的土壤侵蚀程度。
本发明还提供一种应用于黑土区小流域侵蚀模数的土壤侵蚀图生成装置,包括:
背景要素获取单元,用于通过API接口在地理信息平台Arcgis中获取待测算区域的侵蚀图斑背景要素;背景要素包括面积、坡度、植被、土地利用以及土壤。
侵蚀模数计算单元,用于根据所述侵蚀图斑背景要素,使用土壤流失方程对每一个地块计算待测算区域的土壤侵蚀模数。
侵蚀程度确定单元,用于对侵蚀模数计算单元得到的各个地块的土壤侵蚀模数进行加权平均,得到全流域土壤侵蚀模数,并将所述全流域土壤侵蚀模数作为待测算区域的土壤侵蚀程度。
侵蚀图生成单元,用于使用土壤侵蚀程度对待测算土壤的图形区域进行标注,并将包含待测算土壤图形的土壤侵蚀图进行打印。
本发明的有益效果为:1、公式的复杂度较低,计算程序的执行效率较高;2、通过各地类面积与全流域的面积百分比进行具体的计算,不会增加额外的基础资料搜集难度,使用成本较低。
附图说明
图1为本发明的一种应用于东北黑土区小流域土壤侵蚀模数的计算方法的流程图。
具体实施方式
具体实施方式一:本实施方式的一种应用于东北黑土区小流域土壤侵蚀模数的计算方法,如图1所示,包括:
步骤一、在地理信息平台Arcgis中获取待测算区域的侵蚀图斑背景要素;背景要素包括面积、坡度、植被、土地利用以及土壤。
步骤二、根据所述侵蚀图斑背景要素,使用土壤流失方程对流域范围内各种土地利用类型的土壤侵蚀模数进行推导。
步骤三、对步骤二中计算得到的各个地类的土壤侵蚀模数进行加权平均,得到全流域土壤侵蚀模数,并将所述全流域土壤侵蚀模数作为该流域的土壤侵蚀强度。
具体实施方式二:本实施方式与具体实施方式一不同的是:土壤流失方程为:
Meti=RiKiLiSiCiPi
式中,Meti为第i个地块的土壤流失量,Ri为第i个地块的降雨侵蚀力因子,Ki为第i个地块的土壤可蚀性因子,Li为第i个地块的坡长因子,Si为第i个地块的坡度因子,Ci为第i个地块的植被覆盖和经营管理因子,Pi为第i个地块的水土保持措施因子。上述因子是从侵蚀图斑背景要素中得到的。
其它步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一或二不同的是:步骤三具体为:
使用公式
进行计算,其中Me为全流域的土壤侵蚀模数;Wti为相应的侵蚀强度地块面积占全流域的面积百分比;将全流域的土壤侵蚀模数Me作为待测算区域的土壤侵蚀模数的计算结果。
由此可见,本发明提供的计算方法简便可行,可操作性很强,只需要使用面积的比重进行计算,不会增加基础资料的搜集难度。且经过多年的实践证明,该方法科学有效,运用该方法进行的水土保持规划设计、方案编制等工作成果,均得到国家级以及省级专家的一致认可。
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式提供一种应用于黑土区小流域侵蚀模数的土壤侵蚀图生成装置,包括:
背景要素获取单元,用于通过API接口在地理信息平台Arcgis中获取待测算区域的侵蚀图斑背景要素;背景要素包括面积、坡度、植被、土地利用以及土壤。即本实施方式可以通过地理信息平台Arcgis提供的接口获取必要的信息,用于后续的运算。
侵蚀模数计算单元,用于根据所述侵蚀图斑背景要素,使用土壤流失方程对每一个地块计算待测算区域的土壤侵蚀模数。
侵蚀程度确定单元,用于对侵蚀模数计算单元得到的各个地块的土壤侵蚀模数进行加权平均,得到全流域土壤侵蚀模数,并将所述全流域土壤侵蚀模数作为待测算区域的土壤侵蚀程度。
侵蚀图生成单元,用于使用土壤侵蚀程度对待测算土壤的图形区域进行标注,并将包含待测算土壤图形的土壤侵蚀图进行打印。即本实施方式可以将测算完的数据以一定的形式在图中标注出来,然后将这张图打印出来。标注的方式可以是很多种,例如直接将数值写在图中相应图形区域处,或者通过颜色深浅来表示侵蚀程度的大小等。
本实施方式是具体实施方式一对应的程序模块,原理与具体实施方式一相同,此处不再详述。
具体实施方式五:本实施方式与具体实施方式四不同的是:土壤流失方程为:
Meti=RiKiLiSiCiPi
式中,Meti为第i个地块的土壤流失量,Ri为第i个地块的降雨侵蚀力因子,Ki为第i个地块的土壤可蚀性因子,Li为第i个地块的坡长因子,Si为第i个地块的坡度因子,Ci为第i个地块的植被覆盖和经营管理因子,Pi为第i个地块的水土保持措施因子。
本实施方式是具体实施方式二对应的程序模块,原理与具体实施方式二相同,此处不再详述。
其它步骤及参数与具体实施方式四相同。
具体实施方式六:本实施方式与具体实施方式四或五不同的是:
侵蚀程度确定单元具体通过公式
进行计算,其中Me为全流域的土壤侵蚀模数;Wti为第i个地块占全流域的面积百分比;将全流域的土壤侵蚀模数Me作为作为待测算区域的土壤侵蚀程度。
本实施方式是具体实施方式三对应的程序模块,原理与具体实施方式三相同,此处不再详述。
其它步骤及参数与具体实施方式四或五相同。
具体实施方式七:本实施方式与具体实施方式四至六之一不同的是:
侵蚀图生成单元中,使用土壤侵蚀程度对待测算土壤的图形区域进行标注具体为:
将土壤侵蚀程度的数值显示在待测算土壤的图形区域内;或
使用不同深浅的颜色对待测算图形区域进行填充,土壤侵蚀程度越高,颜色越深。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (7)
1.一种应用于黑土区小流域侵蚀模数的计算方法,其特征在于,包括:
步骤一、在地理信息平台Arcgis中获取待测算区域的侵蚀图斑背景要素;
步骤二、根据所述侵蚀图斑背景要素,使用土壤流失方程对每一个地块计算待测算区域的土壤侵蚀模数;
步骤三、对步骤二中计算得到的各个地块的土壤侵蚀模数进行加权平均,得到全流域土壤侵蚀模数,并将所述全流域土壤侵蚀模数作为待测算区域的土壤侵蚀程度。
2.根据权利要求1所述的应用于黑土区小流域侵蚀模数的计算方法,其特征在于,土壤流失方程为:
Meti=RiKiLiSiCiPi
式中,Meti为第i个地块的土壤流失量,Ri为第i个地块的降雨侵蚀力因子,Ki为第i个地块的土壤可蚀性因子,Li为第i个地块的坡长因子,Si为第i个地块的坡度因子,Ci为第i个地块的植被覆盖和经营管理因子,Pi为第i个地块的水土保持措施因子。
3.根据权利要求1或2所述的应用于黑土区小流域侵蚀模数的计算方法,其特征在于,步骤三具体为:
使用公式
<mrow>
<mi>M</mi>
<mi>e</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>Met</mi>
<mi>i</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>Wt</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
进行计算,其中Me为全流域的土壤侵蚀模数;Wti为第i个地块占全流域的面积百分比;将全流域的土壤侵蚀模数Me作为作为待测算区域的土壤侵蚀程度。
4.一种应用于黑土区小流域侵蚀模数的土壤侵蚀图生成装置,其特征在于,包括:
背景要素获取单元,用于通过API接口在地理信息平台Arcgis中获取待测算区域的侵蚀图斑背景要素;背景要素包括面积、坡度、植被、土地利用以及土壤;
侵蚀模数计算单元,用于根据所述侵蚀图斑背景要素,使用土壤流失方程对每一个地块计算待测算区域的土壤侵蚀模数;
侵蚀程度确定单元,用于对侵蚀模数计算单元得到的各个地块的土壤侵蚀模数进行加权平均,得到全流域土壤侵蚀模数,并将所述全流域土壤侵蚀模数作为待测算区域的土壤侵蚀程度;
侵蚀图生成单元,用于使用土壤侵蚀程度对待测算土壤的图形区域进行标注,并将包含待测算土壤图形的土壤侵蚀图进行打印。
5.根据权利要求4所述的应用于黑土区小流域侵蚀模数的土壤侵蚀图生成装置,其特征在于,土壤流失方程为:
Meti=RiKiLiSiCiPi
式中,Meti为第i个地块的土壤流失量,Ri为第i个地块的降雨侵蚀力因子,Ki为第i个地块的土壤可蚀性因子,Li为第i个地块的坡长因子,Si为第i个地块的坡度因子,Ci为第i个地块的植被覆盖和经营管理因子,Pi为第i个地块的水土保持措施因子。
6.根据权利要求4所述的应用于黑土区小流域侵蚀模数的土壤侵蚀图生成装置,其特征在于,
侵蚀程度确定单元具体通过公式
<mrow>
<mi>M</mi>
<mi>e</mi>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>n</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>Met</mi>
<mi>i</mi>
</msub>
<mo>&times;</mo>
<msub>
<mi>Wt</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
进行计算,其中Me为全流域的土壤侵蚀模数;Wti为第i个地块占全流域的面积百分比;将全流域的土壤侵蚀模数Me作为作为待测算区域的土壤侵蚀程度。
7.根据权利要求4所述的应用于黑土区小流域侵蚀模数的土壤侵蚀图生成装置,其特征在于,侵蚀图生成单元中,使用土壤侵蚀程度对待测算土壤的图形区域进行标注具体为:
将土壤侵蚀程度的数值显示在待测算土壤的图形区域内;或
使用不同深浅的颜色对待测算图形区域进行填充,土壤侵蚀程度越高,颜色越深。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711268476.3A CN108021752A (zh) | 2017-12-05 | 2017-12-05 | 一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711268476.3A CN108021752A (zh) | 2017-12-05 | 2017-12-05 | 一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108021752A true CN108021752A (zh) | 2018-05-11 |
Family
ID=62078501
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711268476.3A Pending CN108021752A (zh) | 2017-12-05 | 2017-12-05 | 一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108021752A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859866A (zh) * | 2020-07-21 | 2020-10-30 | 中国水利水电科学研究院 | 一种反映水土流失不同强度间转变的方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101216481A (zh) * | 2008-01-04 | 2008-07-09 | 广东省生态环境与土壤研究所 | 一种反映区域土壤自然侵蚀程度的方法 |
CN103940974A (zh) * | 2014-02-19 | 2014-07-23 | 西北农林科技大学 | 基于gis的中尺度流域土壤侵蚀时空动态分析方法 |
CN104699962A (zh) * | 2015-02-14 | 2015-06-10 | 广东省生态环境与土壤研究所 | 一种土壤侵蚀模数计算方法 |
CN105160465A (zh) * | 2015-08-25 | 2015-12-16 | 四川省宇行规划设计有限公司 | 一种基于3s技术的退耕还林方法 |
KR101789120B1 (ko) * | 2016-06-23 | 2017-11-15 | (주)해동기술개발공사 | 모바일 gis를 이용한 필지별 토사재해 정보 시스템 |
-
2017
- 2017-12-05 CN CN201711268476.3A patent/CN108021752A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101216481A (zh) * | 2008-01-04 | 2008-07-09 | 广东省生态环境与土壤研究所 | 一种反映区域土壤自然侵蚀程度的方法 |
CN103940974A (zh) * | 2014-02-19 | 2014-07-23 | 西北农林科技大学 | 基于gis的中尺度流域土壤侵蚀时空动态分析方法 |
CN104699962A (zh) * | 2015-02-14 | 2015-06-10 | 广东省生态环境与土壤研究所 | 一种土壤侵蚀模数计算方法 |
CN105160465A (zh) * | 2015-08-25 | 2015-12-16 | 四川省宇行规划设计有限公司 | 一种基于3s技术的退耕还林方法 |
KR101789120B1 (ko) * | 2016-06-23 | 2017-11-15 | (주)해동기술개발공사 | 모바일 gis를 이용한 필지별 토사재해 정보 시스템 |
Non-Patent Citations (1)
Title |
---|
陈玉东 等: "基于GIS和USLE的龙墩水库小流域土壤侵蚀评估研究", 《土壤》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111859866A (zh) * | 2020-07-21 | 2020-10-30 | 中国水利水电科学研究院 | 一种反映水土流失不同强度间转变的方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi et al. | The effect of land use/cover change on surface runoff in Shenzhen region, China | |
Melesse et al. | Storm runoff prediction based on a spatially distributed travel time method utilizing remote sensing and GIS 1 | |
Bahrawi et al. | Soil erosion estimation using remote sensing techniques in Wadi Yalamlam Basin, Saudi Arabia | |
Heng et al. | Using artificial neural network to estimate sediment load in ungauged catchments of the Tonle Sap River Basin, Cambodia | |
Svoray et al. | Catchment scale analysis of the effect of topography, tillage direction and unpaved roads on ephemeral gully incision | |
Fan et al. | Impact of C factor of USLE technique on the accuracy of soil erosion modeling in elevated mountainous area (case study: the Tibetan plateau) | |
Santillan et al. | Beyond flood hazard maps: Detailed flood characterization with remote sensing, gis and 2D modelling | |
Kumar | Assessing the impact of climate change on groundwater resources | |
Benzer | Using the geographical information system and remote sensing techniques for soil erosion assessment | |
Gurjar et al. | Assessment of the natural flow regime and its variability in a tributary of Ganga River: Impact of land use and land cover change | |
Zhang | Using LiDAR-DEM based rapid flood inundation modelling framework to map floodplain inundation extent and depth | |
CN108021752A (zh) | 一种应用于黑土区小流域侵蚀模数的计算方法及土壤侵蚀图生成装置 | |
Pawattana et al. | Analytical hierarchical process (AHP)-based flood water retention planning in Thailand | |
Wickenkamp et al. | A multiscale approach to predicting soil erosion on cropland using empirical and physically based soil erosion models in a geographic information system | |
Aiello et al. | Modelling spatially–distributed soil erosion through remotely–sensed data and GIS | |
Wilkinson | Testing the capability of a sediment budget model for targeting remediation measures to reduce suspended-sediment yield | |
Shahzeidi et al. | Assessment rate of soil erosion by GIS (case study Varmishgan, Iran) | |
Singh et al. | Mapping and prediction of surface run-off using SCS-CN method | |
Hu et al. | A GIS-based modeling approach for fast assessment of soil erosion by water at regional scale, Loess Plateau of China | |
Ganapathi et al. | Assessing the impact of land use and land cover changes on the water balances in an urbanized peninsular region of India | |
Alvarez et al. | A coastal erosion risk assessment framework | |
Zhu et al. | Urban flood damage estimation using GIS and remote sensing | |
Damian et al. | Erosion assessment modeling using the SATEEC GIS model on the Prislop catchment | |
Chen | Assessment and improvement of TELEMAC-2D routines for urban flood simulation | |
Zhang | Evaluation of Sediment Yield and Risk Management of Sediment Disaster With the Natural Resources Conservation Service-Curve Number Method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 150080 No. 78 Yanxing Road, Nangang District, Harbin City, Heilongjiang Province Applicant after: Heilongjiang Prov. Water Conservancy Science Inst. Address before: 150070 No. 50 Jingjiangxi Road, Daoli District, Harbin City, Heilongjiang Province Applicant before: HEILONGJIANG PROVINCIAL INSTITUTE OF SOIL SCIENCE |
|
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180511 |