首页 > 代码库 > 中国农田土壤重金属富集状况及其空间分布研究(张小敏)

中国农田土壤重金属富集状况及其空间分布研究(张小敏)

1、统计重金属数据于Excel表格中,命名为coordinate.xls,这里仅以Cr为例进行记录说明。

记录号经度纬度Cd(镉)(二声)
166甘肃省白银市 104.1836.5520.59
64甘肃省天水市天水果园种植地区105.7234.580.07
3广东省广州市 113.2723.130.32
4广东省佛山市 113.1223.020.32
5广东省江门市 113.0822.580.24
6广东省惠州市 114.4223.120.19
80广东省汕头市 116.6823.350.154
31广东省深圳市 114.0522.550.39
14湖南省长沙市 112.9328.236.9
15湖南省株洲市 113.1327.830.61
29湖南省湘潭市 112.9327.830.46
121湖南省怀化市 11027.570.67
147湖南省沅江市 112.2728.930.23
148湖南省衡阳市水口山铅锌矿112.626.573.55
153湖南省长沙市 112.9328.230.39
160湖南省长株潭 113.0327.750.36
151湖南省衡阳市某区112.5726.90.51
8重庆市  106.5529.570.492
16河南省洛阳市 112.4534.620.21
34河南省郑州市 113.6234.750.156
37河南省安阳市滑县114.5235.580.1079
371河南省安阳市安阳县114.3536.10.0931
372河南省安阳市武陟县113.3835.10.0627
373河南省安阳市新乡县113.835.20.0991
374河南省安阳市延津县114.235.150.0971
41河南省新乡市 113.935.348
58银川市兴庆区 106.2838.680.2078
60江西省赣州市赣县(王母渡镇东埠头钨矿区)11525.870.19
122辽宁省葫芦岛市 120.8340.7260.94
20辽宁省沈阳市 123.4341.80.26
45湖北省孝感市孝南区113.9230.921.47
57湖北省武汉市城郊114.330.70.26
149湖北省大冶市大冶矿区及其周边114.9730.14.05
21山东省青岛市 120.3836.07130
94山东省临沂市 118.3535.050.033
66山东省泰安市 117.0836.20.091
46山东省聊城市阳谷县(十五里园镇)115.7836.120.13
50香港  114.0822.21.3
99浙江省金华市富硒区蒋堂镇119.2528.70.46
123浙江省湖州市 120.0830.90.25
63安徽省铜陵市凤凰山117.8230.934.36
25安徽省合肥市 117.2531.830.23
61新疆乌鲁木齐市 87.6243.820.53
65天津市宁河宁河117.8239.330.21
651天津市静海静海116.9238.930.15
652天津市宝坻宝坻117.339.720.2
653天津市蓟县蓟县117.440.050.31
654天津市武青1武青1117.0339.380.215
655天津市东丽东丽117.339.080.63
656天津市北辰北辰117.1339.221.24
657天津市西青西青11739.130.29
658天津市津南津南117.3838.980.39
33山西省太原市 112.5537.870.14
73上海市黄浦江松江区黄埔大桥(饮用水源)121.2231.030.16
74上海市 宝山区121.4831.40.56
76上海市 宝山区121.4831.40.19
1421福建省龙岩上杭县 116.4225.050.301
1422福建省泉州晋江市 118.5824.820.213
1423福建省漳州龙文区 117.7224.520.141
1424福建省宁德霞浦县 12026.880.056
30福建省漳州市 117.6524.520.35
129河北省保定市清苑县115.4838.770.53527
143内蒙古鄂尔多斯市棋盘井106.8839.430.1169
1431内蒙古鄂尔多斯市薛家湾111.2339.870.5876
1432内蒙古鄂尔多斯市树林召110.0340.40.0661
1433内蒙古鄂尔多斯市东胜11039.820.061
11内蒙古包头市 109.8340.650.18
103江苏省徐州市 117.1834.270.14
104江苏省连云港市 119.2234.60.147
105江苏省宿迁市 118.2833.970.13
106江苏省盐城市 120.1533.350.123
107江苏省淮安市 119.0233.620.139
108江苏省扬州市 119.432.40.154
109江苏省泰州市 119.9232.450.144
110江苏省南通市 120.8831.980.139
111江苏省苏州市 120.5831.30.213
112江苏省无锡市 120.331.570.187
113江苏省常州市 119.9531.780.177
114江苏省镇江市 119.4532.20.185
115江苏省南京市 118.7832.070.174
155江苏省常熟市 120.7531.650.145
117江苏省张家港市 120.5531.870.183
72陕西省榆林市神木县110.538.831.18
92陕西省延安市洛川县(京兆乡)109.4335.770.12
125黑龙江省双城市 126.3245.370.092
126黑龙江省黑河市张地营子乡127.4850.250.077
150黑龙江省大庆市 125.0346.580.126
136黑龙江省哈尔滨市市郊126.5345.80.13
58海南省五指山平均值109.5218.780.745
138贵州省黔南布依族苗族自治州贵定县107.2326.580.156
13贵州省贵阳市 106.6326.650.28
98贵州省遵义市正安县107.4328.550.38
144四川省都江堰市蒲阳镇(金凤乡)103.62310.94505
95四川省广元市利州区105.8332.430.18
162北京市顺义区 116.6540.130.136
84北京市顺义区 116.6540.130.39
93吉林省龙井市龙池村129.4242.770.6
159江苏省张家港 120.5531.870.2

2、导入到arcgis中生成要素类或者.shp文件,不详细说明具体过程。

3、将数据复制,粘贴到GS软件,如图,这里分别点XCoord,YCoord,Z,自动弹出对话框,输入列名,这一步很重要,会自动检测数据。

4、数据是否符合正态分布检测,不符合的需要进行一定的数据变换。

5、进行半方差分析,开始计算,要记录模型参数,后续步骤进行Kriging插值时会用到。

6、进行Kriging插值,这里进行插值的参数一定要和上图的模型参数相同。

 

7、按照掩膜裁剪结果图,用中国底图作为掩膜来裁剪。

 

中国农田土壤重金属富集状况及其空间分布研究(张小敏)