一篇文章带你搞定单因子污染指数和卡梅罗污染指数（11年土壤重金属污染示例讲解含代码）

一、题目分析

给出8种主要重金属元素在该城区的空间分布，并分析该城区内不同区域重金属的污染程度

在问题一中，根据三次方程内插法，模拟得到该城区的地势等高图及重金属元素污染的空间分布图。然后采用单因子污染指数和改进版的卡梅罗污染指数评价法对该城区的污染程度进行评价。

二、重金属的元素空间分布的代码

首先，我们由所给采样点的地理位置坐标，建立该城区的地势等高图和功能区分布图，==由于所给数据的离散性和随机性，我们得到的是一个离散的点图，这样的图无法很好地呈现该城区重金属分布的总体状况。因此我们先将城区网格化，采用Cubic （双三次插值），根据已知数据进行插值，对图进行修补。

注意这里插值，可选cubic或者v4插值，cubic真正的数值会有不少Nan值

graddata 插值使用时，注意坐标范围的选择，这里的选取表中 x 值和 y 值的最大值作为限定值

%不同功能区的分布散点图
clear,clc;
%绘图
A = xlsread(