CN108332719A - 坡耕地土壤耕作侵蚀速率快速估算方法 - Google Patents
坡耕地土壤耕作侵蚀速率快速估算方法 Download PDFInfo
- Publication number
- CN108332719A CN108332719A CN201810012961.2A CN201810012961A CN108332719A CN 108332719 A CN108332719 A CN 108332719A CN 201810012961 A CN201810012961 A CN 201810012961A CN 108332719 A CN108332719 A CN 108332719A
- Authority
- CN
- China
- Prior art keywords
- soil
- sloping upland
- farming
- periods
- erosion rate
- 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
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
- G01C11/06—Interpretation of pictures by comparison of two or more pictures of the same area
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Image Processing (AREA)
Abstract
一种坡耕地土壤耕作侵蚀速率快速估算方法,包括:利用无人机获取坡耕地在不同耕作时期前后的影像数据;根据不同耕作时期前后的高密度点云数据生成不同耕作时期前后的高精度数字高程数据及正射影像;利用正射影像进行目视解译,以确定土壤净侵蚀区边界;根据不同耕作时期前后的高精度数字高程数据利用土壤净侵蚀区边界,对不同耕作时期前后的高精度数字高程数据进行剪裁;利用ArcGIS软件中提供的填挖方工具,计算不同耕作时期内,坡耕地耕作侵蚀过程中的土壤净流失体积;采用模型法通过耕作前后地形变化进行计算,得出模拟耕作试验中土壤流失体积,结合土壤容重和耕作次数推算耕作侵蚀速率。
Description
技术领域
本发明涉及一种快速估算坡耕地土壤耕作侵蚀速率的方法,具体涉及利用无人机影像生成高密度点云数据和正射影像数据,监测坡耕地径流小区在耕作前后,坡顶净侵蚀区地形变化的方法,用于估算坡耕地土壤耕作侵蚀速率。
背景技术
紫色土一方面具有较好的母质肥力性成为宝贵的农业耕作土壤,同时受限于其有机质含量低、容重高、分散率高、团聚体含量少等缺点,使紫色土区成为仅次于黄土地区土壤流失最严重的地区之一。紫色土坡耕地因强烈的人为耕作活动,导致坡耕地土层变薄,养分流失加剧,生产力低下,严重阻碍了当地农业的可持续发展。因此,开展紫色土坡耕地耕作侵蚀研究,对研究坡耕地地形演化过程、土壤质量改变、土壤养分转移与重分布、土壤有机质含量及土壤碳储量变化等有重要影响。
耕作侵蚀是发生在坡耕地上的一种重要侵蚀类型。耕作侵蚀是在坡耕地景观内,由于农耕工具和重力作用而使土壤发生向下坡运动与向上坡、运动(依赖于耕作方向),导致净余土壤量向下坡传输、堆积,重新分配的过程目前针对紫色土耕作侵蚀已经进行了大量的研究,结果表明耕作侵蚀部位主要集中在坡顶凸起部位,经过多年长期耕作后,坡顶地形发生改变,坡顶土壤沿耕作方向发生位移,形成土壤净侵蚀区。目前耕作侵蚀速率的计算方法主要有物理示踪法、137Cs示踪剂法、模型法和梯级法。但上述诸多方法存在采集数据方法繁琐、数据精度不高、对原始地形破坏严重等缺陷。
目前,针对坡耕地地区,如何快速精确计算坡耕地在不同耕作年限条件下的土壤侵蚀速率,尚无高效可靠的方法。因此,需要发明一种通过使用无人机近景摄影测量,获取坡耕地耕作前后地形照片,生成正射影像和高程点云数据,插值得到较为精确的两组DEM数据,结合目视解译,确定耕作侵蚀边界,计算土壤净侵蚀体积,以此来计算土壤耕作侵蚀速率的方法。
发明内容
本发明技术的目的是提供一种适合于坡耕地地区,结合无人机影像和高精度DEM数据,推算坡耕地土壤耕作侵蚀速率的方法,并产生一定的经济、社会和生态效益,属于水土保持学领域。
本发明中的基于近景摄影测量的坡耕地土壤耕作侵蚀速率快速估算方法,包括:
步骤1,利用无人机获取坡耕地在不同耕作时期前后的影像数据;
步骤2,使用RTK-GPS获取坡耕地的地面控制点;
步骤3,根据所述地面控制点分别对所述不同耕作时期前后的影像数据进行对齐处理,然后分别提取高密度点云数据;
步骤4,根据不同耕作时期前后的的高密度点云数据生成不同耕作时期前后的高精度数字高程数据及正射影像;
步骤5,通过ArcGIS软件利用正射影像进行目视解译,以确定土壤净侵蚀区边界;
步骤6,根据不同耕作时期前后的所述高精度数字高程数据利用土壤净侵蚀区边界,对不同耕作时期前后的高精度数字高程数据进行剪裁,以使二者的数据面积严格相等,具有相同的行列号;
步骤7,利用ArcGIS软件中提供的填挖方工具,计算不同耕作时期内,坡耕地耕作侵蚀过程中的土壤净流失体积;
步骤8,采用模型法通过耕作前后地形变化进行计算,得出模拟耕作试验中土壤流失体积,结合土壤容重和耕作次数推算耕作侵蚀速率。
优选地,所述影像数据包括RTK-GPS地理位置信息。
优选地,所述影像数据间的重叠度应在40%左右,便于后期影像拼接及数据生成。
优选地,所述控制点可选择坡耕地道路交叉点、人为设置的地面标识点等在影像上易分辨且较精细的特征点。
优选地,步骤8采用下式推算耕作侵蚀速率:
其中,Rt为坡面耕作侵蚀模数(t/(hm2·tillage pass));Vt为耕作后坡耕地净流失土壤体积(m3);ρb为土壤容重(1 311kg/m3);A为坡耕地面积(m2);N为耕作次数。
本发明以水土保持学和土壤学原理为依据,结合人工目视解译,快速计算坡耕地耕作后土壤净流失体积,并以此推算坡耕地耕作侵蚀速率,只需要使用无人机获取坡耕地耕作前后地形影像数据及地面控制点坐标,通过软件即可完成影像自动配准,生成高密度点云数据和正射影像。配合人工目视判读勾画和数据空间分析,即可计算任意耕作时间内坡耕地土壤耕作侵蚀速率,操作简单,数据精度较高,未对原始地形进行人为扰动,具有较强的可操作性。
附图说明
图1是本发明的流程图。
具体实施方式
下面,以一个具体的实施方式对本发明中快速估算坡耕地土壤耕作侵蚀速率的方法进行详细说明:
首先,利用无人机获取坡耕地在不同耕作时期前后的影像数据。针对坡耕地坡度较陡、传统实地测量方法繁琐,人工测量对地形破坏严重的特点,利用无人机合理规划航迹,获取坡耕地清晰无云的影像,影像需要包含地面控制点位置信息,影像间重叠度应在40%左右,便于后期影像拼接及数据生成。
其次,使用RTK-GPS获取坡耕地的地面控制点。使用RTK-GPS在坡耕地测得地面实际控制点坐标,控制点一般选择在耕地道路交叉点、人为设置地面标识点等在影像上易分辨且较精细的特征点。
再次,提取高密度点云数据,生成正射影像和高精度DEM数据。使用AgisoftPhotoscan软件加载无人机拍摄的照片,导入RTK-GPS测得的控制点,软件自动将照片对齐,生成高密度点云数据和正射影像数据。
最后,计算土壤耕作侵蚀流失体积。将耕作前后生成的坡耕地高程点云数据,通过空间插值生成高精度数字高程数据(DEM前,DEM后)。通过在ArcGIS 10.3软件中对正射影像进行目视解译,确定土壤净侵蚀区边界。在ArcGIS 10.3中使用土壤净侵蚀区边界对DEM前和DEM后进行剪裁,使两期DEM数据面积严格相等,具有相同的行列号。利用ArcGIS 10.3中提供的填挖方工具,计算在不同耕作时期内,坡耕地耕作侵蚀过程中的土壤净流失体积(Vt)。
采用模型法通过耕作前后地形变化计算得出模拟耕作试验中土壤流失体积,结合土壤容重和耕作次数推算耕作侵蚀速率:
其中,Rt为坡面耕作侵蚀模数(t/(hm2·tillage pass));Vt为耕作后坡耕地净流失土壤体积(m3);ρb为土壤容重(1 311kg/m3);A为坡耕地面积(m2);N为耕作次数。
本发明的有益效果是:
1、科学性
以水土保持学和土壤学原理为依据,针对传统计算耕作侵蚀率中数据采集方法繁琐、数据精度不高,对坡耕地原始地形扰动较大等不足,提出了应用无人机通过近景摄影测量方法,快速获取坡耕地在任意耕作期限前后地形影像,生成高密度点云数据和正射影像数据,插值生成高精度DEM数据,结合人工目视解译,快速计算坡耕地耕作后土壤净流失体积,并以此推算坡耕地耕作侵蚀速率,研究结果表明该方法科学有效。
2、可操作性强
本发明只需要使用无人机获取坡耕地耕作前后地形影像数据及地面控制点坐标,通过软件即可完成影像自动配准,生成高密度点云数据和正射影像。配合人工目视判读勾画和数据空间分析,即可计算任意耕作时间内坡耕地土壤耕作侵蚀速率,操作简单,数据精度较高,未对原始地形进行人为扰动,具有较强的可操作性。
Claims (5)
1.一种基于近景摄影测量的坡耕地土壤耕作侵蚀速率快速估算方法,其特征在于,包括:
步骤1,利用无人机获取坡耕地在不同耕作时期前后的影像数据;
步骤2,使用RTK-GPS获取坡耕地的地面控制点;
步骤3,根据所述地面控制点分别对所述不同耕作时期前后的影像数据进行对齐处理,然后分别提取高密度点云数据;
步骤4,根据不同耕作时期前后的的高密度点云数据生成不同耕作时期前后的高精度数字高程数据及正射影像;
步骤5,通过ArcGIS软件利用正射影像进行目视解译,以确定土壤净侵蚀区边界;
步骤6,根据不同耕作时期前后的所述高精度数字高程数据利用土壤净侵蚀区边界,对不同耕作时期前后的高精度数字高程数据进行剪裁,以使二者的数据面积严格相等,具有相同的行列号;
步骤7,利用ArcGIS软件中提供的填挖方工具,计算不同耕作时期内,坡耕地耕作侵蚀过程中的土壤净流失体积;
步骤8,采用模型法通过耕作前后地形变化进行计算,得出模拟耕作试验中土壤流失体积,结合土壤容重和耕作次数推算耕作侵蚀速率。
2.根据权利要求1所述的基于近景摄影测量的坡耕地土壤耕作侵蚀速率快速估算方法,其特征在于,所述影像数据包括RTK-GPS地理位置信息。
3.根据权利要求2所述的基于近景摄影测量的坡耕地土壤耕作侵蚀速率快速估算方法,其特征在于,所述影像数据间的重叠度应在40%左右,便于后期影像拼接及数据生成。
4.根据权利要求2所述的基于近景摄影测量的坡耕地土壤耕作侵蚀速率快速估算方法,其特征在于,所述控制点可选择坡耕地道路交叉点、人为设置的地面标识点等在影像上易分辨且较精细的特征点。
5.根据权利要求2所述的基于近景摄影测量的坡耕地土壤耕作侵蚀速率快速估算方法,其特征在于,步骤8采用下式推算耕作侵蚀速率:
其中,Rt为坡面耕作侵蚀模数(t/(hm2·tillage pass));Vt为耕作后坡耕地净流失土壤体积(m3);ρb为土壤容重(1 311kg/m3);A为坡耕地面积(m2);N为耕作次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810012961.2A CN108332719A (zh) | 2018-01-06 | 2018-01-06 | 坡耕地土壤耕作侵蚀速率快速估算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810012961.2A CN108332719A (zh) | 2018-01-06 | 2018-01-06 | 坡耕地土壤耕作侵蚀速率快速估算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108332719A true CN108332719A (zh) | 2018-07-27 |
Family
ID=62924753
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810012961.2A Pending CN108332719A (zh) | 2018-01-06 | 2018-01-06 | 坡耕地土壤耕作侵蚀速率快速估算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108332719A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109033459A (zh) * | 2018-08-30 | 2018-12-18 | 中国科学院地理科学与资源研究所 | 一种土壤容重数据空间格网化构建方法 |
CN110243347A (zh) * | 2019-06-25 | 2019-09-17 | 西华师范大学 | 一种基于近景摄影测量的侵蚀沟监测方法 |
CN110375668A (zh) * | 2019-07-08 | 2019-10-25 | 西北农林科技大学 | 基于点云数据的黄土坡面微地形曲面重构方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006024686A1 (es) * | 2004-07-26 | 2006-03-09 | Consejo Superior De Investigaciones Científicas | Procedimiento para la discriminación y cuantificación de olivar con cubiertas vegetales mediante teledetección con sensores de alta resolución espacial |
CN201740777U (zh) * | 2010-04-06 | 2011-02-09 | 长江水利委员会长江科学院 | 坡面土壤侵蚀快速测定系统 |
CN102749071A (zh) * | 2012-04-24 | 2012-10-24 | 北京林业大学 | 一种基于无人机航空摄影监测土壤侵蚀的方法 |
CN105004725A (zh) * | 2015-08-04 | 2015-10-28 | 珠江水利委员会珠江水利科学研究院 | 一种水土保持综合治理土壤侵蚀变化量实时定量监测方法 |
CN105205855A (zh) * | 2015-09-14 | 2015-12-30 | 西北农林科技大学 | 一种水土流失情况测定的方法 |
CN106815432A (zh) * | 2017-01-17 | 2017-06-09 | 中国科学院、水利部成都山地灾害与环境研究所 | 矿山排土场边坡土壤侵蚀速率估算方法 |
-
2018
- 2018-01-06 CN CN201810012961.2A patent/CN108332719A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006024686A1 (es) * | 2004-07-26 | 2006-03-09 | Consejo Superior De Investigaciones Científicas | Procedimiento para la discriminación y cuantificación de olivar con cubiertas vegetales mediante teledetección con sensores de alta resolución espacial |
CN201740777U (zh) * | 2010-04-06 | 2011-02-09 | 长江水利委员会长江科学院 | 坡面土壤侵蚀快速测定系统 |
CN102749071A (zh) * | 2012-04-24 | 2012-10-24 | 北京林业大学 | 一种基于无人机航空摄影监测土壤侵蚀的方法 |
CN105004725A (zh) * | 2015-08-04 | 2015-10-28 | 珠江水利委员会珠江水利科学研究院 | 一种水土保持综合治理土壤侵蚀变化量实时定量监测方法 |
CN105205855A (zh) * | 2015-09-14 | 2015-12-30 | 西北农林科技大学 | 一种水土流失情况测定的方法 |
CN106815432A (zh) * | 2017-01-17 | 2017-06-09 | 中国科学院、水利部成都山地灾害与环境研究所 | 矿山排土场边坡土壤侵蚀速率估算方法 |
Non-Patent Citations (1)
Title |
---|
杨超,等: "基于无人机影像快速估算矿山排土场边坡土壤侵蚀速率的方法", 《水土保持通报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109033459A (zh) * | 2018-08-30 | 2018-12-18 | 中国科学院地理科学与资源研究所 | 一种土壤容重数据空间格网化构建方法 |
CN109033459B (zh) * | 2018-08-30 | 2021-06-15 | 中国科学院地理科学与资源研究所 | 一种土壤容重数据空间格网化构建方法 |
CN110243347A (zh) * | 2019-06-25 | 2019-09-17 | 西华师范大学 | 一种基于近景摄影测量的侵蚀沟监测方法 |
CN110375668A (zh) * | 2019-07-08 | 2019-10-25 | 西北农林科技大学 | 基于点云数据的黄土坡面微地形曲面重构方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schmidt et al. | Comparison of DEM data capture and topographic wetness indices | |
CN107421501A (zh) | 一种结合倾斜摄影、rtk和bim技术的土石方测量计算方法 | |
CN101750015B (zh) | 基于数字图像技术的沙石坑土方量测量方法 | |
CN106643529A (zh) | 基于无人机影像的山区农作物生长高度快速测量的方法 | |
CN108332719A (zh) | 坡耕地土壤耕作侵蚀速率快速估算方法 | |
Mourato et al. | Improving a DSM obtained by unmanned aerial vehicles for flood modelling | |
Xu et al. | Evaluation of gully head retreat and fill rates based on high-resolution satellite images in the loess region of China | |
CN105445137A (zh) | 一种基于摄影的坡面小区尺度细沟侵蚀三维动态监测方法 | |
CN107392908A (zh) | 作物株高测量方法及系统、田间定标装置 | |
CN106815432A (zh) | 矿山排土场边坡土壤侵蚀速率估算方法 | |
CN102565303A (zh) | 一种冲沟沟头溯源侵蚀速率的快速监测方法 | |
CN105606049B (zh) | 一种相对水面的水田农业机具作业姿态测量方法 | |
CN110470275A (zh) | 一种基于uav航测地形数据测量干枯河床沙波形态参数的方法 | |
CN108180897A (zh) | 坡耕地土壤水力侵蚀速率快速估算方法 | |
CN116448080B (zh) | 一种基于无人机倾斜摄影辅助土方开挖施工的方法 | |
Haug et al. | Geodetic mass balance of the western Svartisen ice cap, Norway, in the periods 1968–1985 and 1985–2002 | |
Wang et al. | Comparison of changes in glacier area and thickness on the northern and southern slopes of Mt. Bogda, eastern Tianshan Mountains | |
Kienast et al. | Geodetic measurements on glaciers and rock glaciers in the Hohe Tauern National Park (Austria) | |
Massuel et al. | Deriving bathymetries from unmanned aerial vehicles: a case study of a small intermittent reservoir | |
Römkens et al. | Soil erosion by water and tillage | |
Uysal et al. | 3D Shoreline mapping using an unmanned aerial vehicle | |
Dedrick et al. | Land forming for irrigation | |
Liu et al. | Use of UAV Photogrammetry for Monitoring Topographic Changes in the Tottori Sand Dunes, Japan | |
Luppichini et al. | A METHOD BASED ON BEACH PROFILE ANALYSIS FOR SHORELINE IDENTIFICATION | |
CN111272145B (zh) | 基于无人机影像的高潜水位沉陷区作物绝产边界确定方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180727 |
|
RJ01 | Rejection of invention patent application after publication |