一种土壤铬含量与土地利用类型空间相关性的分析方法技术

本发明专利技术公开了一种土壤铬含量与土地利用类型空间相关性的分析方法，其步骤为：在分布图的研究区域上布置点位，监测每个点位上的土壤重金属铬含量数据；通过GIS空间数据插值得到研究区域土壤重金属铬含量的空间分布栅格图；以网格布点，获取网格节点；读取土地利用类型和土壤铬含量数值；计算土壤铬含量分布熵和土地利用类型分布熵以及联合分布熵；计算土壤铬含量和土地利用类型的相关系数；确定研究区域土壤铬含量和土地利用类型之间相关性。本发明专利技术在定性变量和定量变量之间建立了定量的空间相关性分析，分析结果可定量揭示土地利用类型对土壤铬含量分布的影响大小，可为区域重金属污染防治方案和治理决策提供参考，免于经验决策的盲目性。

【技术实现步骤摘要】
一种土壤铬含量与土地利用类型空间相关性的分析方法
本专利技术属于环境保护
，具体涉及一种土壤铬含量与土地利用类型空间相关性的分析方法。
技术介绍
随着矿产资源的不合理开发与利用、化肥农药的大量施用、工业化和城镇化的迅速发展，土壤污染日益严重。重金属是严重危害生态安全的土壤环境污染物之一，它们能够改变土壤化学组成，直接或间接破坏土壤的生态结构，其在土壤中的累积和迁移不仅影响动植物生长发育与农产品质量安全，而且可通过食物链直接危害人类健康。土壤中的铬(Cr)元素是典型的有毒有害重金属元素，其在土壤中的含量与自然条件和人类活动有着密切的联系。土壤本身就含有一定量的Cr，不同的母质、成土过程会影响土壤Cr元素的含量，这是土壤Cr元素的背景含量。而人类活动如工业生产过程会向环境中排放含Cr的废水、废气、废渣，农业生产过程中向农田施用含Cr农药和化肥，从而导致土壤中Cr含量的增加。关于土壤重金属污染的空间分布与评价是当今国内外学者及管理部门关注的重点，它能够对区域土壤重金属污染情况有总体掌握，为防治决策提供参考。进一步而言，土壤重金属污染的空间分布情况受什么因素影响呢？有研究发现城市区域土壤表层重金属污染的空间格局与人口密度特征、环境要素、交通条件、产业结构与布局、土地利用等内外部约束因素的空间分异特征密切相关。关于区域土壤中重金属元素的自然背景值已有定量化的调查结果，但人类活动是如何影响土壤重金属污染的，对于二者之间的关系很少有定量化的研究，其难点在于人类活动难以定量化表征。所以，需要一个新的技术方案来解决这个问题。
技术实现思路
专利技术目的：为了克服现有技术中存在的不足，提供一种土壤铬含量与土地利用类型空间相关性的分析方法，从定性上，研究区域内土壤重金属铬含量的空间分布；从定量上，得到土壤铬含量与土地利用类型的空间相关系数。技术方案：为实现上述目的，本专利技术提供一种土壤铬含量与土地利用类型空间相关性的分析方法，包括如下步骤：S1：获取用于研究区域土地利用类型的分布图；S2：在分布图的研究区域上布置p1个点位，监测每个点位上的土壤重金属铬含量数据；S3：基于步骤S2的测量数据，通过GIS空间数据插值得到研究区域土壤重金属铬含量的空间分布栅格图；S4：考虑到研究区域的面积及土地利用斑块的大小，以网格布点，共获取p2个网格节点，其中p2＞＞p1；S5：基于步骤S1中的分布图及步骤S3中得到的空间分布栅格图，读取步骤S4中p2个网格节点处对应的土地利用类型和土壤铬含量数值；S6：基于步骤S5读取的土地利用类型和土壤铬含量数值，计算土壤铬含量分布熵和土地利用类型分布熵；S7：确定土壤铬含量和土地利用类型的联合分布熵；S8：计算土壤铬含量和土地利用类型的相关系数K；S9：根据相关系数K确定研究区域土壤铬含量和土地利用类型之间相关性的大小。进一步的，所述步骤S3中空间分布栅格图的具体获取方式为：基于ArcGIS10.3平台空间分析模块IDW插值功能，将铬含量的点矢量图层插值为面栅格图层，插值过程利用被插值区域附近的12个点位数据，按照距离远近赋予权重，从而获得被插值区域的铬含量数据，经过空间插值处理，得到研究区域土壤铬含量空间分布栅格图。进一步的，所述的步骤S2中布置的p1个点位应较均匀的分布于研究区域内。进一步的，所述步骤S5中在网格中间部位取点，读取该点位的土壤铬含量数据以及对应的土地利用类型，得到p2个点位的数据系列。进一步的，所述步骤S6中土壤铬含量分布熵和土地利用类型分布熵的计算方法如下：其中X是随机变量(本专利技术中可指代土壤铬含量或土地利用类型二者中的任一变量)的整体，也称为信源，P(xi)为发生事件xi的概率，xi为变量X的子集，n为可能发生的事件(或状态)总数，式中对数的底a一般取2。进一步的，所述步骤S7中土壤铬含量和土地利用类型的联合分布熵的确定方法为：设定信源X代表土壤铬含量，信源Y代表土地利用类型，这两个信源X和Y为二维随机变量，则它们的联合分布概率表达为：p(xi,yj)(i＝1,2,...,n；j＝1,2,...,m)X，Y的联合信息熵H(X,Y)表达为：其中p(xi,yj)为两个信源X，Y的联合分布概率；xi，yi为信源X、Y的子集，即可能发生的事件，n，m为可能发生的事件(或状态)总数。进一步的，所述步骤S8中土壤铬含量和土地利用类型的相关系数K的计算方法为：因为H(X,Y)≤H(X)+H(Y)，因此K的取值范围为[0，1]。进一步的，所述步骤S9中相关性大小的具体判断方法为：当K＝0时，表示X，Y不相关；K＞0，表示X，Y具有相关性；K越大则说明二者之间的相关程度越高。有益效果：本专利技术与现有技术相比，可快速确定研究区域内土壤重金属铬含量的空间分布，并定量获得土壤铬含量与土地利用类型的空间相关系数，在定性变量和定量变量之间建立了定量的空间相关性分析，分析结果可定量揭示土地利用类型对土壤铬含量分布的影响大小，可为区域重金属污染防治方案和治理决策提供参考，免于经验决策的盲目性。附图说明图1为本专利技术分析方法的流程图；图2为土壤监测点位分布图；图3是土壤重金属Cr含量空间分布图。具体实施方式下面结合附图和具体实施例，进一步阐明本专利技术。本专利技术提供一种土壤Cr(铬)含量与土地利用类型空间相关性的分析方法，如图1所示，其包括如下步骤：S1：收集研究区域土地利用现状的分布图。S2：在研究区域较均匀布置52个点位，监测每个点位上的土壤重金属Cr含量数据，具体的布置图如图2所示。S3：基于步骤S2中的测量数据，通过数据插值得到研究区域土壤重金属Cr含量的空间分布栅格图，其具体方式为：基于ArcGIS10.3平台空间分析模块IDW插值功能，将铬含量的点矢量图层插值为面栅格图层，插值过程利用被插值区域附近的12个点位数据，按照距离远近赋予权重，从而获得被插值区域的铬含量数据，经过空间插值处理，得到研究区域土壤铬含量空间分布栅格图，具体如图3所示。S4：考虑到研究区域的面积及土地利用斑块的大小，以1.2km×1.2km大小的网格布点，共获取210个网格节点；S5：基于步骤S1中研究区域土地利用现状分布图及步骤3中得到的土壤Cr含量空间分布栅格图，读取步骤S4中210个网格节点处对应的土地利用类型和土壤Cr含量数值，建立概率矩阵，如表1所示：表1S6：计算得到Cr含量熵值H(X)为2.7，土地利用类型熵值H(Y)为2.63，计算得到二者之间的联合熵为4.62。S7：计算土壤Cr含量和土地利用类型的相关系数K＝0.13；S8：根据步骤S7中的相关系数确定研究区域土壤Cr含量和土地利用类型的空间相关系数K＝0.13，说明二者具有空间相关性，证明研究区域土壤Cr含量分布情况已受到本文档来自技高网...

【技术保护点】
1.一种土壤铬含量与土地利用类型空间相关性的分析方法，其特征在于：包括如下步骤：/nS1：获取用于研究区域土地利用类型的分布图；/nS2：在分布图的研究区域上布置p

【技术特征摘要】
1.一种土壤铬含量与土地利用类型空间相关性的分析方法，其特征在于：包括如下步骤：
S1：获取用于研究区域土地利用类型的分布图；
S2：在分布图的研究区域上布置p1个点位，监测每个点位上的土壤重金属铬含量数据；
S3：基于步骤S2的测量数据，通过GIS空间数据插值得到研究区域土壤重金属铬含量的空间分布栅格图；
S4：考虑到研究区域的面积及土地利用斑块的大小，以网格布点，共获取p2个网格节点，其中p2＞＞p1；
S5：基于步骤S1中的分布图及步骤S3中得到的空间分布栅格图，读取步骤S4中p2个网格节点处对应的土地利用类型和土壤铬含量数值；
S6：基于步骤S5读取的土地利用类型和土壤铬含量数值，计算土壤铬含量分布熵和土地利用类型分布熵；
S7：确定土壤铬含量和土地利用类型的联合分布熵；
S8：计算土壤铬含量和土地利用类型的相关系数K；
S9：根据相关系数K确定研究区域土壤铬含量和土地利用类型之间相关性的大小。


2.根据权利要求1所述的一种土壤铬含量与土地利用类型空间相关性的分析方法，其特征在于：所述步骤S3中空间分布栅格图的具体获取方式为：基于ArcGIS10.3平台空间分析模块IDW插值功能，将铬含量的点矢量图层插值为面栅格图层，插值过程利用被插值区域附近的12个点位数据，按照距离远近赋予权重，从而获得被插值区域的铬含量数据，经过空间插值处理，得到研究区域土壤铬含量空间分布栅格图。


3.根据权利要求1所述的一种土壤铬含量与土地利用类型空间相关性的分析方法，其特征在于：所述步骤S5中在网格中间部位取点，读取该点位的土壤铬含量数据以及对应的土地利用类型，得到p2个...

【专利技术属性】
技术研发人员：李秀霞，张子豪，胡逸宁，杨艺，权仕帆，
申请(专利权)人：南京信息职业技术学院，
类型：发明
国别省市：江苏;32