# Analysis of Soil Fertility Based on FUMF Algorithm

Hang Chen, Guifen Chen, Yating Hu, Liying Cao, Lixia Cai, Sisi Yang

► To cite this version:

Hang Chen, Guifen Chen, Yating Hu, Liying Cao, Lixia Cai, et al.. Analysis of Soil Fertility Based on FUMF Algorithm. 9th International Conference on Computer and Computing Technologies in Agriculture (CCTA), Sep 2015, Beijing, China. pp.564-573, 10.1007/978-3-319-48357-3_53 . hal-01557799

HAL Id: hal-01557799

https://inria.hal.science/hal-01557799

Submitted on 6 Jul 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Analysis of Soil Fertility Based on FUMF Algorithm

Hang Chen<sup>1,2,a</sup>, Guifen Chen<sup>2,b,*</sup>, Yating Hu<sup>2,c</sup>, Liying Cao<sup>2,d</sup>, Lixia Cai<sup>2,e</sup>, Sisi Yang<sup>1,f</sup>

<sup>1</sup>Institute of Scientific and Technical Information of Jilin, Changchun 130033, China

<sup>2</sup>College of Information Technology, Jilin Agricultural University, Changchun 130118, China

<sup>a</sup>chenhang0811@163.com, <sup>b</sup>guifchen@163.com, <sup>*</sup>Corresponding author <sup>c</sup> 65447539@qq.com, <sup>d</sup>45515189@qq.com, <sup>e</sup>419513823@qq.com, <sup>f</sup>1046235989@qq.com

**Abstract.** The soil nutrition is an important indicator of soil fertility. The method K-means and FCM are always used to evaluating the soil fertility, but the cluster number need to be set, and the outlier couldn't be eliminated accurately, and there is the deviation between the real result and the soil fertility. So the paper applied the FUMF to analysis the soil nutrient data of Nong An county for eight years, 2005-2012. The result show that the low fertility soils gradually decreased from 2005 to 2012 by precision fertilization, and the moderate and high fertility soil was rising, the overall soil fertility of Nong An had improved significantly. The analysis result was consistent with the actual situation, The FUMF algorithm is proved that was an effective evaluate method of the soil fertility evaluation. It has the practical significance to analyze the large number of soil fertility of high complexity and interactive, it also provided the technical support for precision fertilization decision-making.

**Keywords:** fuzzy clustering, soil fertility, fertility analysis, precision fertilization, FUMF

## 1 Introduction

Soil nutrient content is an important symbol of fertility and productivity of arable land, also it is an important indicator of soil fertility evaluation. With the arrival of precision agriculture era, spatial variability and correlations of wide variety agricultural data which have complex links relationship are more significantly. The attendant massive, diverse and dynamic changes, incomplete, uncertain and a series of characteristics.

Since the 1990s, Data Mining and geographic information systems technology in the agricultural sector has been increasingly widely used. DM and GIS technology can effectively statistics and analysis of massive, complex data. DM Clustering algorithms can dig out the knowledge of soil fertility evaluation from soil nutrient data analysis. Li Lianghou etc. put forward the application of clustering analysis which is in site classification and soil fertility evaluation[3]. Zheng Chao etc. improved rough K- means algorithm, and put forward the rough K- means clustering algorithm based on density weighted[5]. Chen Guifen etc. put forward a weighted spaces fuzzy dynamic clustering algorithm, and proved the validity of method in evaluation of soil fertility[6]. But conventional K-means, FCM and other clustering algorithms have some limitations on soil fertility evaluation. Such as K- means is hard clustering algorithm that can only get a hard divide. Although FCM can get fuzzy clustering divide, both algorithms require artificially set the number of clusters. So it can not eliminate outlier accurately or
solve the problem of soil fertility data including complex, dynamic, and interactive fuzzy. Whatever,
the clustering results presence of a certain error with the real fertility. For this reason the paper use
FUMF algorithm to analyze and evaluate soil fertility.

National measuring territories precise fertilizer projects in Jilin Province for over 10 years. During this period a large number of soil samples were collected and sample of soil nutrients were determined and analyzed. All of this could lay the foundation for soil fertility status by using DM and GIS technology. Thus, this paper use large amounts of data by successive years of soil testing precise fertilizer projects that from Nong An county in Jilin Province. Then, we use GIS and Matlab technical conducted a rapid unsupervised multiscale fuzzy clustering for soil nutrient data from 2005 to 2012. The results show that FUMF algorithm is an effective method for soil fertility evaluation and has practical significance when analyze large amounts of high complexity, strong interaction soil fertility factors. So, it is can provide a technical support for the precise fertilization decision.

2 Experiments and Methods

2.1 The situation of research area

Nong'an is located in Songliao Plain, Changchun, Jilin. specific in northwest of Changchun city away
from 60 kilometers, north latitude 43 ° 54'-44 ° 56 ', longitude 124 ° 32'-125 ° 45'; The zone is in the
temperate semi-humid continental monsoon climate. So, monsoon features obviously, four distinct
seasons, abundant sunshine, less rainfall and the annual average temperature of 4.6 ℃, annual average
sunshine hours 2590h, the average annual rainfall 507mm; On the one hand, there are diverse
landforms such as high mesa, mesa, two terraces, a terrace, floodplain, sand dunes, depressions, gullies
and so on. Thus, most soil is chernozem, meadow soil and black soil; On the other hand, they grow
corn sorghum, wheat, millet and soybean and other crops production as the mainstay. It is arguably one
of the country's important commodity grain production bases and its total grain production ranked first
in the major grain-producing counties.

## 2.2 Collection and analysis of sample data

On the basis of field research, we are cooperation with cropland capacity survey quality evaluation
office and considering soil types, land use, topography, cropping patterns, management measures and
production level and other factors according second national soil survey. After that, we determine the
sampling point through DGPS and RS systems. Then we can comprehensive analyze the survey plots
of soil testing precision fertilizing work from 2005 to 2012. It collected 23,976 samples, sampling map
of soil nutrients in Figure 1, Figure 2.

<figure><img src="image_2.png" /></figure>
<caption>**Fig.1.** Nong an soil fertility sampling map</caption>

Samples were collected depth from 0 to 20cm. Random multi-point sampling within the same plots. Whatever, after mixing the soil by quartering, we take 1.5 or 2kg bagging spare. Then, take it back to the laboratory for spare through dry naturally pulverized and sieving. Ultimately, index measuring soil nitrogen, phosphorus and potassium and other nutrients in which a total of 26 kinds of soil types. So, in this paper we analyze the collected 23976 data and calculate maximum, minimum, average value of nitrogen (N), phosphorus (P), potassium (K) according to the different soil types. The data of 2010 is shown in Table 1:

<table><caption>**Fig.2.** Part of the grid sampling map from Nong an</caption><thead><tr><th>town name</th><th>country name</th><th>plot name</th><th>Alkalystic N (mg/kg<sup>-1</sup>)</th><th>Available P (mg/kg<sup>-1</sup>)</th><th>Available K (mg/kg<sup>-1</sup>)</th><th>latitude</th><th>longitude</th></tr></thead><tbody><tr><td>Nong An</td><td>Xi Haolai</td><td>Jia Bei</td><td>154</td><td>28</td><td>208</td><td>44.58417</td><td>125.28982</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Jia Bei</td><td>136</td><td>31.3</td><td>217</td><td>44.49895</td><td>125.25127</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Jia Bei</td><td>132</td><td>12.8</td><td>217</td><td>44.4988</td><td>125.2509</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>132</td><td>16.3</td><td>198</td><td>44.49926</td><td>125.25072</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>121</td><td>21.4</td><td>208</td><td>44.5004</td><td>125.24821</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Gong Lubei</td><td>118</td><td>31.8</td><td>227</td><td>44.50205</td><td>125.24968</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Lu Nan</td><td>114</td><td>37.9</td><td>227</td><td>44.50206</td><td>125.24738</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>114</td><td>28</td><td>237</td><td>44.5024</td><td>125.24905</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Yi Jiedi</td><td>114</td><td>23.3</td><td>217</td><td>44.5054</td><td>125.24378</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>114</td><td>21.3</td><td>208</td><td>44.50712</td><td>125.24472</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>114</td><td>19.1</td><td>169</td><td>44.50855</td><td>125.24402</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>121</td><td>24.2</td><td>160</td><td>44.50305</td><td>125.23638</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>118</td><td>24.8</td><td>140</td><td>44.50648</td><td>125.24862</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>118</td><td>31.3</td><td>179</td><td>44.50616</td><td>125.24878</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>118</td><td>43.6</td><td>140</td><td>44.50523</td><td>125.24845</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>118</td><td>32.8</td><td>160</td><td>44.50502</td><td>125.24878</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>110</td><td>39.3</td><td>169</td><td>44.51012</td><td>125.2503</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>132</td><td>32.8</td><td>188</td><td>44.50512</td><td>125.25013</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>129</td><td>37.8</td><td>179</td><td>44.50522</td><td>125.25045</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>121</td><td>23.3</td><td>179</td><td>44.50755</td><td>125.25172</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>114</td><td>34.1</td><td>227</td><td>44.51202</td><td>125.2518</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>114</td><td>36.8</td><td>237</td><td>44.5058</td><td>125.24893</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>114</td><td>22.7</td><td>246</td><td>44.50618</td><td>125.24475</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Gong Lubei</td><td>114</td><td>30</td><td>217</td><td>44.50412</td><td>125.23805</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>118</td><td>35.8</td><td>208</td><td>44.50195</td><td>125.2379</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Gong Lubei</td><td>129</td><td>21.4</td><td>237</td><td>44.50885</td><td>125.24368</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Gong Lubei</td><td>140</td><td>30.8</td><td>227</td><td>44.51057</td><td>125.24563</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Lu Bei</td><td>114</td><td>36.3</td><td>217</td><td>44.51185</td><td>125.24018</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>118</td><td>30.8</td><td>208</td><td>44.51535</td><td>125.23695</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>118</td><td>36.3</td><td>198</td><td>44.50502</td><td>125.2354</td></tr></tbody></table>

<caption>**Table 1.** Soil fertility data</caption>
<table><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>118</td><td>42.8</td><td>227</td><td>44.50268</td><td>125.23857</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>103</td><td>32.1</td><td>198</td><td>44.51077</td><td>125.25425</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>San Jiedi</td><td>132</td><td>15.8</td><td>227</td><td>44.50702</td><td>125.2542</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Er Jiedi</td><td>165</td><td>10.8</td><td>237</td><td>44.51245</td><td>125.25438</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Gong Lubei</td><td>147</td><td>30.8</td><td>217</td><td>44.51385</td><td>125.2544</td></tr><tr><td>Nong An</td><td>Xi Haolai</td><td>Bei Dapian</td><td>143</td><td>21.3</td><td>198</td><td>44.5057</td><td>125.25578</td></tr></table>

According to soil grading standards of second soil survey, soil nutrients are divided into six levels, such as shown in Table 2.

<table><caption>**Table 2. Soil nutrient grading standards**</caption><thead><tr><th>project</th><th>Level 1</th><th>Level 2</th><th>Level 3</th><th>Level 4</th><th>Level 5</th><th>Level 6</th></tr></thead><tbody><tr><td>Alkalystic N(mg/kg)</td><td>&gt;150</td><td>120-150</td><td>90-120</td><td>60-90</td><td>30-60</td><td>&le;30</td></tr><tr><td>Available P(mg/kg)</td><td>&gt;40</td><td>20-40</td><td>10-20</td><td>5-10</td><td>3-5</td><td>&le;3</td></tr><tr><td>Available K(mg/kg)</td><td>&gt;200</td><td>150-200</td><td>100-150</td><td>50-100</td><td>30-50</td><td>&le;30</td></tr></tbody></table>

According to preliminary results of the analysis, we begin to accurate classification of soil fertility through data mining.

# 3 Results and Discussion

## 3.1 Fast unsupervised multiscale fuzzy clustering (FUMF) algorithm

First of all, we clustering the N, P, K three indicators of 23976 data by FUMF algorithm, the purpose is to eliminate the isolated samples point of each index; Then these three indicators data were normalized; Finally, set the parameters of weighted dimensional data for clustering analysis by FUMF.

We can accelerate UMF algorithm through nearest neighbor criterion and get FUMF. Well, FUMF method is divided into two stages:

The first stage: re-expression data by using the nearest neighbor criterion, the data is divided into $\bar{n}$ disjoint subsets $S_j$, Each subset's data represented by its representative point $c_j$ which is as a whole.

The second stage, implementation of weighted UMF algorithm.

FUMF algorithm is as follows:

Step1. The re-expression data, initialize the $m=1, c_m = \{x_i\}$, $i=2$ to $N$: $d(x_i, c_k) = \min_{1 \le j \le m} d(x_i, c_j)$, If $d(x_i, c_k) > Q$ and $m < q$ then $m=m+1, c_m = \{x_i\}$ Else $c_k = c_k \cup \{x_i\}$.

Step 2. Clustering UMF, set $j=1$, set a threshold $e > 0$ and $v^{(0)} = c_j$, then using the updated formula 3:
$$v^{(l+1)} = \frac{\sum_{k=1}^{\bar{n}} n_k \cdot c_k \cdot d(v^{(l)}, c_k)}{\sum_{k=1}^{\bar{n}} n_k \cdot d(v^{(l)}, c_k)} \quad (3)$$

Calculate convergence point of $c_j$, denote as $p_j$. If $j < \bar{n}$, then $j = 1 + j$, repeat step 2.

Step 3. If $\|p_a - p_b\| e$, The $S_a$ and $S_b$ of the data points into a class; otherwise, divided into different classes.

## 3.2 Soil nutrient content analysis

Through statistical analyze 23,976 samples of soil nutrient content, we summarizes the changes of soil nutrients from early, metaphase and anaphase data. As shown in Table 3:

<table><caption>**Table 3.** Nong an Part of the township of soil nutrient content in different years descriptive statistics</caption><thead><tr><th>YEAR</th><th>Index</th><th>Alkalystic N<br>(mg kg-1)</th><th>Available P<br>(mg kg-1)</th><th>Available K<br>(mg kg-1)</th></tr></thead><tbody><tr><td rowspan="5">2005</td><td>Range</td><td>9.0-193.0</td><td>2.0-145.0</td><td>8.0-960.0</td></tr><tr><td>Mean</td><td>105.79</td><td>18.7</td><td>138.1</td></tr><tr><td>CV(%)</td><td>21.83</td><td>67.15</td><td>51.15</td></tr><tr><td>Number of samples</td><td>2297</td><td>2297</td><td>2297</td></tr><tr><td>Range</td><td>19.0-211.0</td><td>2.1-96.9</td><td>40.0-382.0</td></tr><tr><td rowspan="5">2009</td><td>Mean</td><td>131.84</td><td>22.25</td><td>184.22</td></tr><tr><td>CV(%)</td><td>19.35</td><td>54.74</td><td>25.87</td></tr><tr><td>Number of samples</td><td>5115</td><td>5115</td><td>5115</td></tr><tr><td>Range</td><td>88.0-213.0</td><td>5.5-82.6</td><td>15.0-413.0</td></tr><tr><td>Mean</td><td>135.98</td><td>24.68</td><td>168.47</td></tr><tr><td rowspan="5">2012</td><td>CV(%)</td><td>13.77</td><td>41.97</td><td>28.19</td></tr><tr><td>Number of samples</td><td>6329</td><td>6329</td><td>6329</td></tr></tbody></table>

## 3.3 FUMF Analysis

Taking into account the soil sampling N, P, K three indicators' observed values are different. Data will inevitably be contaminated during sampling that resulting in some isolated points. Therefore before cluster analysis of soil nutrient, we need to pre-processing the data set. Pretreatment divided into the following steps:

(1) Executing clustering algorithm for N, P, K three indicators respectively. If it contains a small number of data points when clustering, indicating this category may be constituted by isolated point. In the experiment, we analyze categories which data points lower than 20 and delete those isolated points which beyond the normal range of values.

(2) Because of three indicators of N, P, K have differences in dimension as raw data. Therefore, each of
these three indicators were normalized so that the mean of each index is 0 and variance is 1.
After process the raw data, each sample as a data point for clustering. Due to Evaluation of soil fertility mainly depends on the content of P indexes, and P indexes are generally lower than the value of N, K.
Thus, we should weighted N, P, K as 1:10: 1 before performing clustering algorithm.

(3) Parameter settings: the convergence of the scale parameter is 0.15; convergence precision is $10^{-5}$;
maximum number of iterations is 100; fuzzy factor is m = 2; data reduction parameters is 0.8;
convergence scale parameter 0.14 multiplied mean value; After performing clustering algorithm to
pretreatment and weighted data, using inverse transform to get clustering results.

### 3.3.1 The initial precision Fertilization clustering results

In this paper, we collected 2297 samples from 27 towns in 2005 to establish the experimental data set (remove isolated points of 38 when prepossessing), all of the data come from Bajilei, Bangchai, Binghe, Fuquanlong, Gaojiadian, Halahai and so on. Then we clustering by FUMF. The clustering results shown in Table 4 and Figure 3(a 2005).

<table><caption>Table 4. Clustering results in 2005</caption><thead><tr><td>Category</td><td>Available N (mg kg-1)</td><td>Available P (mg kg-1)</td><td>Available K (mg kg-1)</td><td>The amount of data</td><td>Fertility</td></tr></thead><tbody><tr><td>1</td><td>105.30</td><td>11.30</td><td>119.05</td><td>824</td><td>low</td></tr><tr><td>2</td><td>102.83</td><td>13.87</td><td>122.90</td><td>623</td><td>medium</td></tr><tr><td>3</td><td>105.93</td><td>19.50</td><td>134.85</td><td>427</td><td>medium</td></tr><tr><td>4</td><td>112.61</td><td>64.58</td><td>140.52</td><td>21</td><td>high</td></tr><tr><td>5</td><td>105.14</td><td>54.18</td><td>448.18</td><td>2</td><td>high</td></tr><tr><td>6</td><td>102.83</td><td>32.72</td><td>135.60</td><td>95</td><td>high</td></tr><tr><td>7</td><td>126.20</td><td>47.54</td><td>399.26</td><td>6</td><td>high</td></tr><tr><td>8</td><td>108.14</td><td>38.36</td><td>158.07</td><td>26</td><td>high</td></tr><tr><td>9</td><td>111.16</td><td>40.29</td><td>144.33</td><td>38</td><td>high</td></tr><tr><td>10</td><td>107.81</td><td>44.89</td><td>138.15</td><td>51</td><td>high</td></tr><tr><td>11</td><td>132.25</td><td>60.23</td><td>125.78</td><td>4</td><td>high</td></tr><tr><td>12</td><td>101.39</td><td>29.55</td><td>131.55</td><td>132</td><td>medium</td></tr><tr><td>13</td><td>103.13</td><td>27.59</td><td>129.29</td><td>2</td><td>medium</td></tr><tr><td>14</td><td>107.73</td><td>54.61</td><td>107.13</td><td>5</td><td>high</td></tr><tr><td>15</td><td>132.78</td><td>60.18</td><td>264.45</td><td>3</td><td>high</td></tr></tbody></table>

### 3.3.2 The middle precision Fertilization clustering results

Experimental data sets with 5115 samples from 23 towns in 2009 (remove isolated points of 24 when prepossessing). The data come from Bajielei, Dehui, Gaojiadian, Halahai, Helong and so on. And clustering results table 5 and Figure 3(b 2009).

<caption>Table 5. Clustering results in 2009</caption>
<table><thead><tr><th>Category</th><th>Available N<br>(mg kg-1)</th><th>Available P<br>(mg kg-1)</th><th>Available K<br>(mg kg-1)</th><th>The<br>amount of<br>data</th><th>Fertility</th></tr></thead><tbody><tr><td>1</td><td>123.68</td><td>43.40</td><td>179.97</td><td>102</td><td>high</td></tr><tr><td>2</td><td>131.79</td><td>13.81</td><td>187.19</td><td>1499</td><td>low</td></tr><tr><td>3</td><td>133.58</td><td>22.55</td><td>185.53</td><td>1156</td><td>medium</td></tr><tr><td>4</td><td>135.92</td><td>17.73</td><td>189.62</td><td>957</td><td>medium</td></tr><tr><td>5</td><td>132.67</td><td>15.66</td><td>187.31</td><td>147</td><td>medium</td></tr><tr><td>6</td><td>134.86</td><td>19.69</td><td>186.40</td><td>441</td><td>medium</td></tr><tr><td>7</td><td>133.14</td><td>36.41</td><td>186.40</td><td>272</td><td>high</td></tr><tr><td>8</td><td>134.98</td><td>30.84</td><td>185.22</td><td>398</td><td>high</td></tr><tr><td>9</td><td>132.40</td><td>62.92</td><td>178.31</td><td>6</td><td>high</td></tr><tr><td>10</td><td>138.88</td><td>76.07</td><td>152.41</td><td>6</td><td>high</td></tr><tr><td>11</td><td>148.39</td><td>67.57</td><td>224.91</td><td>10</td><td>high</td></tr><tr><td>12</td><td>111.89</td><td>81.76</td><td>163.05</td><td>16</td><td>high</td></tr><tr><td>13</td><td>109.95</td><td>72.23</td><td>136.23</td><td>3</td><td>high</td></tr><tr><td>14</td><td>116.30</td><td>53.27</td><td>144.66</td><td>50</td><td>high</td></tr><tr><td>15</td><td>112.51</td><td>57.82</td><td>140.77</td><td>14</td><td>high</td></tr><tr><td>16</td><td>123.68</td><td>43.40</td><td>179.97</td><td>11</td><td>high</td></tr><tr><td>17</td><td>131.79</td><td>13.81</td><td>187.19</td><td>3</td><td>high</td></tr></tbody></table>

### 3.3.3 The late Precise fertilization clustering results

Experimental data sets with 6329 samples from 17 towns in 2012 (remove isolated points of 17 when prepossessing). All of the data come from Helong, Qiangang, Bajilei, Fuquanlong, Gaojiadian, Halahai, Huajia and so on. Then we clustering by FUMF. The clustering results shown in Table 6 and Figure 3(c 2012).

<table><caption>Table 6. Clustering results in 2012</caption><thead><tr><th>Category</th><th>Available N<br>(mg kg-1)</th><th>Available P<br>(mg kg-1)</th><th>Available K<br>(mg kg-1)</th><th>The<br>amount<br>of data</th><th>Fertility</th></tr></thead><tbody><tr><td>1</td><td>131.79</td><td>20.12</td><td>170.56</td><td>742</td><td>medium</td></tr><tr><td>2</td><td>134.85</td><td>37.22</td><td>160.72</td><td>561</td><td>high</td></tr><tr><td>3</td><td>133.05</td><td>21.75</td><td>166.18</td><td>399</td><td>medium</td></tr><tr><td>4</td><td>137.10</td><td>31.71</td><td>162.83</td><td>666</td><td>high</td></tr><tr><td>5</td><td>133.89</td><td>25.76</td><td>166.14</td><td>699</td><td>medium</td></tr><tr><td>6</td><td>133.15</td><td>28.58</td><td>166.18</td><td>593</td><td>medium</td></tr><tr><td>7</td><td>134.18</td><td>14.20</td><td>162.77</td><td>1312</td><td>low</td></tr><tr><td>8</td><td>133.13</td><td>23.58</td><td>164.96</td><td>260</td><td>medium</td></tr><tr><td>9</td><td>133.80</td><td>18.54</td><td>169.35</td><td>460</td><td>medium</td></tr><tr><td>10</td><td>135.91</td><td>15.77</td><td>165.41</td><td>316</td><td>medium</td></tr></tbody></table>
<table><tr><td>11</td><td>135.36</td><td>30.17</td><td>163.07</td><td>58</td><td>medium</td></tr><tr><td>12</td><td>134.90</td><td>47.52</td><td>161.72</td><td>204</td><td>high</td></tr><tr><td>13</td><td>140.12</td><td>63.81</td><td>133.31</td><td>25</td><td>high</td></tr><tr><td>14</td><td>134.71</td><td>45.89</td><td>170.48</td><td>4</td><td>high</td></tr><tr><td>15</td><td>136.68</td><td>35.46</td><td>160.15</td><td>4</td><td>high</td></tr><tr><td>16</td><td>131.79</td><td>20.12</td><td>170.56</td><td>9</td><td>high</td></tr></table>

<figure><img src="image_3.png" /><figcaption>Fig.3. Soil nutrient spatial clustering map</figcaption></figure>

### 3.4 Clustering analysis

When compared the clustering results from 2005,2009 to 2012, we can derive trend of soil fertility that soil fertility tend towards equilibrium and rise after precision fertilization.

<table><caption>Table 7. Clustering analysis table and soil fertility</caption><thead><tr><th rowspan="2">Year</th><th rowspan="2">The number of samples</th><th colspan="3">Classification Number</th><th colspan="3">The percentage of each category in the total sample</th></tr><tr><th>High</th><th>medium</th><th>low</th><th>high</th><th>medium</th><th>low</th></tr></thead><tbody><tr><td>2005</td><td>2259</td><td>251</td><td>1184</td><td>824</td><td>11.11%</td><td>52.42%</td><td>36.48%</td></tr><tr><td>2009</td><td>5091</td><td>891</td><td>2701</td><td>1499</td><td>17.50%</td><td>53.05%</td><td>29.44%</td></tr><tr><td>2012</td><td>6312</td><td>1473</td><td>3527</td><td>1312</td><td>23.34%</td><td>55.88%</td><td>20.79%</td></tr></tbody></table>

## 4 Conclusions

Through the clustering results we analyzed, soil fertility tend towards equilibrium and rise after precision fertilization, so it can reflect the trend of soil fertility better. The results can be seen from Table 7:

(1) We clustering according to the parameters which is set by clustering algorithm. The data is from 2005, 2009 and 2012 and the number of samples are 2259,5091 and 6312. Then we can derived its
cluster classification results are consistent with the actual high fertility, the fertility and low fertility referring soil grading standards. So, it could prove that FUMF algorithm is an effective method to soil fertility evaluation.

(2) The data were compared from 2005, 2009 and 2012, the clustering results show that high fertility soils were increased from 8.16% to 13.99% and 15.30%; The second soil fertility were increased from 53.64% to 62.80% and 65.80%; low soil fertility dropped 22.60% from 38.10% to 18.93% respectively. It is shown that from 2005 to 2012, the low soil fertility decreases and other soil fertility increase after precise fertilization. So, soil fertility has improved significantly.

(3) The analysis results are consistent with the actual situation, it is not only shows FUMF algorithm is an effective method for soil fertility evaluation, but also proved that after precise fertilization soil fertility has improved significantly in general. Therefore, we believe that the method is meaningful by using data mining to analyze fertility factors of large the high complexity and, strong interaction data. So it can provide technical support for precision fertilization decisions.

## Acknowledgment

This work was supported by the National "863" High-tech Project (2006AA10A309), Jilin province science and technology development projects(key science and technology research project): "The development and demonstration of corn production's monitoring and traceability system based on Internet of things technology"(20140204045NY), Jilin province agricultural committee projects: "The demonstration and generalization of Corn precision operation system based on Internet of things".

## References

1. Zeitouni K., 2002. A survey of spatial data mining methods databases and statistics point of views. In: Data Warehousing and Web Engineering, edited by Becker S. London: IRM Press, 229-242.

2. L K Sharma, O P Vyas, U S Tiwary, et al. A Novel Approach of Multilevel Positive and Negative Association Rule Mining for Spatial Databases. Machine Learning and Data Mining in Pattern Recognition. Springer Berlin/Heidelberg, 2005,3587:620-629.

3. Li Lianghou, Li Liyue. Application of Clustering Analysis in Site Classification and Soil Fertility Evaluation[J]. Proceedings of 2010 Third International Conference on Education Technology and Training 2010.

4. Guo Xiaojuan, Liu Xiaoxia, Li Xiaoling. Improvement and analysis of hierarchical clustering algorithm[J]. Computer Applications and Software, 25(6):243-244.2005.

5. Zheng Chao, Miao Duoqian, Wang Ruizhi. Improved K- means clustering algorithm based on rough density weighted[J]. Computer Science, 2009.

6. Guifen Cheng, LiYing Cao, Guowei Wang. Application of Weighted spatial fuzzy dynamic clustering algorithm in soil fertility evaluation in the black soil zone. Chinese Agricultural Science, 2009,42(10):3559-3563.

7. ZHANG JIANPING, LIU XIYU. Research and application of K-means algorithm Based on clustering Analysis[J]. Application of Computer, 2007,24(5):166-168.

8. Li Yuancheng, Yin Peipei. Speculative multithreading partitioning algorithm based on fuzzy clustering[J]. Journal of Computers, 2014.37(3): 580-592.

9. YANG YUE, GUO SHUXU. FCM image segmentation algorithm based on kernel function and spatial information[J]. Jilin University (Engineering Science), 2011.41(2): 283-287.
10.LIU FUHCENG, SHI XUEZHEG, YU DONGSHENG.Research of Mapping of soil properties in Taihu Based on geostatistics and GIS--In total soil nitrogen in cartography example[J].Journal of Soil,2004, 41(1): 20-27.

11.ZHAO YANFENG, SHI XUEZHEG, YU DONGSHENG.Discussion of mall-scale spatial variability of soil nutrients and its influencing factors-- In urban Areas in Wuxi City, Jiangsu Province[J].2006, 37(2): 214-219.

12.HUANG SHAOWEN, JIN JIYUN, YANG LIPING.Research Grain Crop spatial variability of soil nutrients and the partition management technology in Region of County[J].Journal of Soil, 2003, 40(1): 79-88.

13.YANG YANLI, SHI XUEZHENG, YU DONGCHENG.Research of region scale soil nutrients spatial variability and influencing factors[J].Geographic Sciences,2008, 28(6): 788-792.

14.GUO XUDONG, FU BOJIE, MA KEMING.Research on Spatial Variability of Soil Fertility Based on GIS and geostatistics---Case By Hebei Zunhua[J].Journal of Applied Ecology,2000, 11(4): 557-563.

15.SUN YUFEN, LU YANSHENG.High-dimensional data flow subspace clustering algorithm based on grid[J].Computer Science, 2007:199-203.

16.Shan Shiming, Yan Yan, Zhang Xianchao.Clustering algorithm based on K most similar clustering subspace[J].Computer Engineering, 2009: 4-6.

17.Qiu Baozhi, Zheng Zhijie,.Clustering algorithm based on the local density and dynamically generated mesh.[J].Computer Engineering,2010, 31(2):385-387..