CN108629059A - 一种基于dac的土壤重金属含量估测方法 - Google Patents

一种基于dac的土壤重金属含量估测方法 Download PDF

Info

Publication number
CN108629059A
CN108629059A CN201710166961.3A CN201710166961A CN108629059A CN 108629059 A CN108629059 A CN 108629059A CN 201710166961 A CN201710166961 A CN 201710166961A CN 108629059 A CN108629059 A CN 108629059A
Authority
CN
China
Prior art keywords
data
heavy metal
content
dem
soil
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
Application number
CN201710166961.3A
Other languages
English (en)
Inventor
史舟
郑胜云
何康
徐烨
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201710166961.3A priority Critical patent/CN108629059A/zh
Publication of CN108629059A publication Critical patent/CN108629059A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Processing Of Solid Wastes (AREA)

Abstract

本发明公开了一种基于DAC的土壤重金属含量估测方法，属于空间预测方法领域。本发明提出了一种基于DAC算法，建立了Hg含量与陆表特征环境变量之间的模型关系，实现了Hg含量的空间分布特征及含量的预测，为后续利用少量采样点来估算整个研究区土壤重金属含量的空间特征提供理论和算法依据。

Description

一种基于DAC的土壤重金属含量估测方法
技术领域
本发明涉及一种基于DAC的土壤重金属含量估测方法，具体涉及到一种基于DAC的土壤重金属含量估测方法。
技术背景
土壤重金属污染具有在土壤中移动性差、滞留时间长、不能被微生物降解等特点，对环境和生物的潜在危害极大。在全球化、粮食安全和生态文明建设背景下，包括重金属污染在内的农产品安全问题(绿色壁垒)已成为农产品质量安全、农业产业化发展和发展现代农业的最大障碍。加强耕地重金属污染控制和治理，指导农产品产地环境分类分级管理、保障我省粮食质量安全，确保农产品消费安全，促进全省现代农业持续健康发展迫在眉睫。
受经济、地形地貌等因素的制约，土壤地面调查采样点布局难免不规则或者样点数量有限，因此，如何基于有限采样点的土壤重金属调查数据来建模预测整个研究区土壤重金属含量的空间分布，成为目前土壤重金属地面调查的研究热点、难点，以及重点。针对该科学问题，提出了一种基于DAC算法，建立了Hg含量与陆表特征环境变量之间的模型关系，并论证Hg含量的可预测性。
发明内容
针对采样点布局不规则或采样点数量有限问题，本发明提出一种基于DAC的土壤重金属含量估测方法。本发明中的一种基于DAC的土壤重金属含量估测方法，基于DAC算法，建立了Hg含量与陆表特征环境变量之间的模型关系，并验证其可预测性。
本发明的具体技术方案如下：
基于DAC的土壤重金属含量估测方法，包括以下步骤：
步骤1)数据获取：获取待测区域地面调查得到的土壤重金属采样点数据、MODIS卫星遥感影像数据、ASTERGDEM卫星遥感影像数据以及待测区域地表温度数据；所述的MODIS卫星遥感影像为MOD13A2数据产品；
步骤2)数据预处理：将步骤1)将ASTER GDEM卫星遥感影像数据进行聚合计算得到DEM数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算得到植被指数数据；将反演而来的地表温度数据聚合计算到空间分辨率为1km；将土壤重金属数据整理入库后提取出Hg含量数据，并进行异常值的剔除处理；
步骤3)根据陆表环境变量与Hg元素含量之间的相关关系，进行分区建模：将步骤2)处理后的Hg含量数据作为因变量，以植被指数数据、DEM数据、地表温度作为自变量，将该数据集进行分类建模预测各个采样点的Hg含量。
作为优选，所述的分区建模过程中所采用的多元回归模型形式为：
其中，N表示参数估算模型中自变量个数；an表示第n个环境变量的系数；a0表示模型参数的常数项系数；yn表示Hg含量预测值；xn表示第n个环境变量；
a0和an的计算公式如下：
其中:k代表土壤重金属采样点的个数；xin代表第i个土壤重金属采样点的第n个环境变量的值，yi代表的是第i个土壤重金属采样点的Hg含量，代表所有土壤重金属采样点中第n个环境变量因子的均值，代表所有土壤地面调查样点中Hg含量的均值。
作为优选，所述的DEM数据的空间分辨率为1km。
作为优选，所述的植被指数数据的空间分辨率为1km。
作为优选，所述的分区建模过程中所采用的多元回归模型具体如下：
(1)当dem≤25.02
YHg＝0.07-0.0036Xlstd-2.2e-005×Xdem-0.03×Xndvi+0.001Xlon
(2)当lat>30.13且dem≤25.02
YHg＝Xndvi Xdem＝-3.66-0.002553×Xdem+0.13×Xlat
(3)当lat≤30.13且dem≤25.02
YHg＝-3.43-0.004382×Xdem+0.089×Xlat+0.038Xlstd+0.013Xlstn+0.08Xndvi
其中YHg是土壤重金属Hg含量，Xdem代表的是空间分辨率为1km的数字高程模型的栅格值，Xndvi代表的是空间分辨率为1km的植被指数栅格值，Xlat代表的是采样点的纬度，Xlon代表的是采样点的经度，Xlstd代表的是白天的地面温度，Xlstn代表的是晚上的地面温度。
本发明提出的一种基于DAC的土壤重金属含量估测方法，分区建立1km土壤调查采样点Hg含量与空间分辨率为1km的植被指数数据、DEM数据、1km的白天地面温度和1km晚上地面温度之间多元回归模型，验证实验表明其实现了Hg含量的预测，为后续利用少量采样点来估算整个研究区土壤重金属含量的空间特征提供了理论和算法依据。
附图说明
图1为实施例中实际值与预测值的拟合结果图。
具体实施方式
下面结合附图和具体实施例对本发明进一步说明。
基于DAC的土壤重金属含量估测方法，包括以下步骤：
步骤1)数据获取：获取待测区域地面调查得到的土壤重金属采样点数据、MODIS卫星遥感影像数据、ASTERGDEM卫星遥感影像数据以及待测区域地表温度数据；所述的MODIS卫星遥感影像为MOD13A2数据产品；
步骤2)数据预处理：将步骤1)将ASTER GDEM卫星遥感影像数据进行聚合计算得到空间分辨率为1km的DEM数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算得到空间分辨率为1km的植被指数数据；将反演而来的地表温度数据聚合计算到空间分辨率为1km；将土壤重金属数据整理入库后提取出Hg含量数据，并进行异常值的剔除处理；
步骤3)根据陆表环境变量与Hg元素含量之间的相关关系，进行分区建模：将步骤2)处理后的Hg含量数据作为因变量，以空间分辨率为1km的植被指数数据、DEM数据、地表温度作为自变量，将该数据集进行分类建模预测各个采样点的Hg含量。
分区建模过程中所采用的多元回归模型形式为：
其中，N表示参数估算模型中自变量个数；an表示第n个环境变量的系数；a0表示模型参数的常数项系数；yn表示Hg含量预测值；xn表示第n个环境变量；
a0和an的计算公式如下：
其中:k代表土壤重金属采样点的个数；xin代表第i个土壤重金属采样点的第n个环境变量的值，yi代表的是第i个土壤重金属采样点的Hg含量，代表所有土壤重金属采样点中第n个环境变量因子的均值，代表所有土壤地面调查样点中Hg含量的均值。
实施例1
选取浙江省作为研究区域(边界以最新的行政区划边界为准)，对2013年土壤重金属Hg含量进行估测。
步骤1)数据获取：获取待测区域的土壤重金属采样点相关数据(包括Hg含量数据、经纬度坐标等)、MODIS卫星遥感影像数据、ASTERGDEM卫星遥感影像数据以及待测区域地表温度数据；其中MODIS卫星遥感影像为MOD13A2数据产品；
步骤2)数据预处理：将步骤1)将ASTER GDEM卫星遥感影像数据进行聚合计算得到空间分辨率为1km的DEM数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算得到空间分辨率为1km的植被指数数据；将反演而来的地表温度数据聚合计算到空间分辨率为1km；将纸质文档的土壤重金属数据整理入库，提取出Hg含量相关数据，并进行异常值的剔除处理。
步骤3)根据陆表环境变量与Hg元素含量之间的内在联系，进行分区建模：将步骤2)处理后的Hg含量数据作为因变量，空间分辨率为1km的植被指数数据、DEM数据、地表温度作为自变量，将该数据集进行分类建模预测各个采样点的Hg含量。
最终得到分区域多元回归模型具体如下：
(1)当dem≤25.02
YHg＝0.07-0.0036Xlstd-2.2e-005×Xdem-0.03×Xndvi+0.001Xlon
(2)当lat>30.13且dem≤25
YHg＝Xndvi Xdem＝-3.657-0.002553×Xdem+0.125×Xlat
(3)当lat≤30.13且dem≤25
YHg＝-3.425-0.004382×Xdem+0.089×Xlat+0.038Xlstd+0.013Xlstn+0.08Xndvi
其中YHg是土壤重金属Hg含量，Xdem代表的是空间分辨率为1km的数字高程模型的栅格值，Xndvi代表的是空间分辨率为1km的植被指数栅格值，Xlat代表的是采样点的纬度，Xlon代表的是采样点的经度，Xlstd代表的是白天的地面温度，Xlstn代表的是晚上的地面温度。
步骤4)土壤重金属Hg含量验证：由步骤3)得到的各个采样点Hg含量预测值与各个点的实测值进行精度对比分析。所选用评价因子为均方根误差、平均绝对误差以及相关系数。各指标的计算公式如下：
式中MAE代表的是平均绝对误差，RMSE代表的是均方根误差，R2代表的是回归相关系数，Yk是地面调查采样点k的含量值，Ok是在地面调查采样点k处的预测值，是所有土壤地面调查样点Hg含量的均值，是在所有土壤地面调查样点Hg含量模型预测值的平均值。
经对比分析评价如下：经对比，如图1所示，土壤重金属Hg含量真实值与预测值相关系数R2为0.88，所以Hg含量的空间分布特征及含量是可预测的，且其精度较高，可为后续利用少量采样点来估算整个研究区土壤重金属含量的空间特征提供理论和算法依据。

Claims (5)

1.一种基于DAC的土壤重金属含量估测方法，其特征在于，包括以下步骤：
步骤1)数据获取：获取待测区域地面调查得到的土壤重金属采样点数据、MODIS卫星遥感影像数据、ASTERGDEM卫星遥感影像数据以及待测区域地表温度数据；所述的MODIS卫星遥感影像为MOD13A2数据产品；
步骤2)数据预处理：将步骤1)将ASTER GDEM卫星遥感影像数据进行聚合计算得到DEM数据；从MOD13A2数据产品中提取植被指数参量，经过异常值剔除处理后，通过聚合计算得到植被指数数据；将反演而来的地表温度数据聚合计算到空间分辨率为1km；将土壤重金属数据整理入库后提取出Hg含量数据，并进行异常值的剔除处理；
步骤3)根据陆表环境变量与Hg元素含量之间的相关关系，进行分区建模：将步骤2)处理后的Hg含量数据作为因变量，以植被指数数据、DEM数据、地表温度作为自变量，将该数据集进行分类建模预测各个采样点的Hg含量。
2.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的分区建模过程中所采用的多元回归模型形式为：
其中，N表示参数估算模型中自变量个数；an表示第n个环境变量的系数；a0表示模型参数的常数项系数；yn表示Hg含量预测值；xn表示第n个环境变量；
a0和an的计算公式如下：
其中:k代表土壤重金属采样点的个数；xin代表第i个土壤重金属采样点的第n个环境变量的值，yi代表的是第i个土壤重金属采样点的Hg含量，代表所有土壤重金属采样点中第n个环境变量因子的均值，代表所有土壤地面调查样点中Hg含量的均值。
3.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的DEM数据的空间分辨率为1km。
4.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的植被指数数据的空间分辨率为1km。
5.如权利要求1所述的基于DAC的土壤重金属含量估测方法，其特征在于，所述的分区建模过程中所采用的多元回归模型具体如下：
(1)当dem≤25.02
YHg＝0.07-0.0036Xlstd-2.2e-005×Xdem-0.03×Xndvi+0.001Xlon
(2)当lat>30.13且dem≤25.02
YHg＝Xndvi Xdem＝-3.66-0.002553×Xdem+0.13×Xlat
(3)当lat≤30.13且dem≤25.02
YHg＝-3.43-0.004382×Xdem+0.089×Xlat+0.038Xlstd+0.013Xlstn+0.08Xndvi
其中YHg是土壤重金属Hg含量，Xdem代表的是空间分辨率为1km的数字高程模型的栅格值，Xndvi代表的是空间分辨率为1km的植被指数栅格值，Xlat代表的是采样点的纬度，Xlon代表的是采样点的经度，Xlstd代表的是白天的地面温度，Xlstn代表的是晚上的地面温度。
CN201710166961.3A 2017-03-20 2017-03-20 一种基于dac的土壤重金属含量估测方法 Pending CN108629059A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710166961.3A CN108629059A (zh) 2017-03-20 2017-03-20 一种基于dac的土壤重金属含量估测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710166961.3A CN108629059A (zh) 2017-03-20 2017-03-20 一种基于dac的土壤重金属含量估测方法

Publications (1)

Publication Number Publication Date
CN108629059A true CN108629059A (zh) 2018-10-09

Family

ID=63686438

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710166961.3A Pending CN108629059A (zh) 2017-03-20 2017-03-20 一种基于dac的土壤重金属含量估测方法

Country Status (1)

Country Link
CN (1) CN108629059A (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100040260A1 (en) * 2007-12-20 2010-02-18 Image Tree Corp. Remote sensing and probabilistic sampling based method for determining the carbon dioxide volume of a forest
CN103544550A (zh) * 2013-11-08 2014-01-29 湖南科技大学 一种金属矿区土-水界面重金属污染负荷的预测方法
CN106021872A (zh) * 2016-05-10 2016-10-12 浙江大学 一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100040260A1 (en) * 2007-12-20 2010-02-18 Image Tree Corp. Remote sensing and probabilistic sampling based method for determining the carbon dioxide volume of a forest
CN103544550A (zh) * 2013-11-08 2014-01-29 湖南科技大学 一种金属矿区土-水界面重金属污染负荷的预测方法
CN106021872A (zh) * 2016-05-10 2016-10-12 浙江大学 一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
BIFENG HU等: "Application of portable XRF and VNIR sensors for rapid assessment of soil heavy metal pollution", 《PLOS ONE》 *
李力生等: "高光谱遥感数据在农作物重金属胁迫调查中的应用", 《安徽农业科学》 *
杨灵玉等: "基于Hyperion影像植被光谱的土壤重金属含量空间分布反演—以青海省玉树县为例", 《应用生态学报》 *

Similar Documents

Publication Publication Date Title
Pande et al. Study of land use classification in an arid region using multispectral satellite images
Grunwald Multi-criteria characterization of recent digital soil mapping and modeling approaches
Kidd et al. Operational sampling challenges to digital soil mapping in Tasmania, Australia
CN106021868B (zh) 一种基于多规则算法的遥感数据降尺度方法
Holden et al. Estimating the carbon stock of a blanket peat region using a peat depth inference model
Mishra et al. Spatial representation of organic carbon and active-layer thickness of high latitude soils in CMIP5 earth system models
CN106021872A (zh) 一种基于低分辨率卫星遥感数据的环境变量动态筛选建模降尺度方法
Hollaway et al. The challenges of modelling phosphorus in a headwater catchment: Applying a ‘limits of acceptability’uncertainty framework to a water quality model
CN113176393B (zh) 一种基于hasm模型的土壤有机碳储量三维估算方法及系统
Kim et al. A heuristic gap filling method for daily precipitation series
Yang et al. Relationship between soil depth and terrain attributes in karst region in Southwest China
Sheng et al. Digital soil mapping to enable classification of the salt-affected soils in desert agro-ecological zones
CN116777079A (zh) 一种基于贝叶斯层间结构模型的沙漠化侵扰灾害预测方法
Jiang et al. Spatiotemporal dynamics of soil erosion risk for Anji County, China
Brito et al. Spatial patterns of magnetic susceptibility optimized by anisotropic correction in different Alisols in southern Amazonas, Brazil
Liou et al. Prediction of soil erosion and sediment transport in a mountainous basin of Taiwan
Zhu et al. Soil salinity dynamics in arid oases during irrigated and non‐irrigated seasons
Laosuwan et al. Urban heat island monitoring and analysis by using integration of satellite data and knowledge based method
Afrasiab et al. Assessing the risk of soil vulnerability to wind erosion through conditional simulation of soil water content in Sistan plain, Iran
CN107798357A (zh) 一种基于knn的村庄用地复垦规划模拟方法
Yu et al. Water eutrophication evaluation based on multidimensional trapezoidal cloud model
Bindajam et al. Evaluating the impact of urbanization on the urban heat islands through integrated radius and non-linear regression approach
CN108629059A (zh) 一种基于dac的土壤重金属含量估测方法
CN108563974A (zh) 一种土壤重金属Hg含量的空间预测方法
CN107727078A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20181009