|
|
||||||||
Dep. of Crop and Soil Sci., Univ. of Georgia, Athens, GA 30602 USA
fchen{at}arches.uga.edu
| ABSTRACT |
|---|
|
|
|---|
. The distribution of the surface SOC concentrations was predicted with two approaches. The first approach was to apply the relationship to individual pixels and then determine the distribution; the second approach was to classify the image and then apply the relationship to determine the class boundaries and means. Eight levels of surface SOC concentrations were classified in both approaches, and there was good agreement between the two approaches with a probability value near one using a paired t-test. The predicted and measured surface SOC concentrations, based on additional soil samples from 31 field locations, were compared using linear regression
. The surface SOC concentrations were correctly classified in 77.4 and 74.2% of cases for the two approaches. The procedures tested were accurate enough to be used for precision farming applications in agricultural fields.
Abbreviations: GPS, global positioning system SOC, soil organic C
| INTRODUCTION |
|---|
|
|
|---|
The dark color of soil is typically associated with high organic-matter concentration and high native fertility. Soils with thick, dark surface horizons are often separated from other soils at the highest categorical level in many soil classification systems, reflecting the differences in the genesis of soils as well as the importance of these soils as a medium for plant growth and independent natural bodies worthy of further study (Schulze et al., 1993). Research has been done concerning the relationships between soil color and soil organic matter. However, many of these studies were based on Munsell color notations for specific soils at specific locations (Alexander, 1969; Steinhardt and Franzmeier, 1979; Schulze et al., 1993) or for the purpose of designing spectral sensors (Pitts et al., 1983; Griffis, 1985; Smith et al., 1987).
There were attempts to quantify relationships between soil color and organic matter concentrations by Brown and O'Neal in the 1920's (Schulze et al., 1993). Later, color charts or tables that described the relationships between soil color and organic-matter concentration were developed by using visual color descriptions or Munsell soil color charts (Shields et al., 1968; Alexander, 1969; Steinhardt and Franzmeier, 1979). Shields et al. (1968) conducted a study of several Ap horizons in an attempt to distinguish between two soils, based on the soil color. Organic C was found to be correlated with soil color for both soils. Alexander (1969) developed a color chart for visually estimating the organic-matter concentration of Ap horizons from more than 300 Illinois soil samples. Steinhardt and Franzmeier (1979) correlated the organic-matter concentrations with the moist soil color for 262 samples of Ap horizons in Indiana. Both papers classified organic-matter concentrations into quantitative categories using the Munsell Color System as standards, and general relationships were developed for visually estimating organic-matter concentrations. Page (1974) used a color-difference meter to examine 96 soils from the Coastal Plain Region of South Carolina and found a curvilinear relationship between reflectance and percent organic matter in the 0 to 5% range. Research has shown that spectroscopic measurement of soil reflectance can give better accuracy in soil color measurement than visual matching (Schulze et al., 1993; Torrent and Barrón, 1993).
Reflectance in various spectral bands has been correlated with soil properties such as soil organic matter. Spectral sensors were designed to measure soil organic matter based on the relationship between light reflectance and soil organic matter (Pitts et al., 1983; Griffis, 1985; Smith et al., 1987; Shonk et al., 1991). Different algorithms were developed to transform the output reflectance to concentration of soil organic matter and soil moisture. Baumgardner et al. (1970) used 197 grid samples for a 25-ha field to correlate the soil organic-matter content to different wavelengths in 12 channels from the visual to infrared range, and a computer printout of soil pattern was generated. It was shown that the organic-matter content can be predicted from light reflectance with a linear or curvilinear relationship in the visual and infrared range (Baumgardner et al., 1970; Leger et al., 1979; Cihlar et al., 1987; Smith et al., 1987; Sudduth and Hummel, 1988; Shonk et al., 1991; Henderson et al., 1992). Research also showed that the relationship between soil organic matter and reflectance is poor if soil samples were collected from large geographic areas or different landscapes, such as soil samples from an entire state (Fernandez et al., 1988; Henderson et al., 1992; Schulze et al., 1993). The cause may be due to different types of parent materials (Henderson et al., 1992).
In previous research, there was no attempt to accurately determine the distribution of surface soil organic C (SOC) concentrations based on the reflected image intensity data for a field, which may be useful for precision farming. Relatively simple and inexpensive methods that would be both more accurate and less expensive than grid sampling are needed to develop maps of surface SOC concentrations. The method should employ only the minimum number of soil samples for organic C analysis to minimize the costs for creating maps. The objective of this study was to map the surface SOC concentrations for a field using an inexpensive remotely sensed image, a color slide, coupled with image processing and auto-classification technology and statistical approaches. A field located in Crisp County, Georgia, was selected for this study, in part because of its range and spatial distribution of surface SOC concentrations.
| Materials and methods |
|---|
|
|
|---|
|
|
To reduce the variance (noise) among the image pixels caused by micro-topography, film processing, and scanning, a low-pass filter was applied to the image with a mask in 5 by 5 cells before examining the relationships between image-intensity values and surface SOC concentrations. This is an average smoothing filter, as follows:
![]() | (1) |
Based on the locations for the 28 soil samples, the pixel values of these 28 locations were determined from the filtered image. The relationship between surface SOC concentrations and the pixel values for the 28 samples was developed by regression analysis. This relationship was applied to the original image, and then an image representing the distribution of surface SOC concentrations for the field was obtained. The result was called Pre_Result1. The filtered image was not used in this case because smoothing would remove real spatial variability in surface SOC concentrations.
An alternative approach was also used to perform a classification to the original image by a minimum-distance clustering algorithm (Jensen, 1986; Lillesand and Kiefer, 1987). This algorithm uses minimum spectral distance to assign a cluster (class) for each candidate pixel. The process begins with an arbitrary number of clusters (classes), and then it processes repetitively until meeting a certain stop condition (or conditions). The input parameters for the method include: (i) the maximum number of clusters to be considered (N); (ii) the convergence threshold (T), which is the maximum percentage of pixels whose class values are allowed to be unchanged between iterations; and (iii) the maximum number of iterations (M).
The process of this algorithm is as follows: (i) Arbitrarily initialize the mean for each of N clusters by simply dividing the image into N groups and then computing the mean for each group. (ii) For each pixel, compute the spectral distance between this pixel and each cluster mean, and assign the pixel to the cluster with the minimum distance between the cluster mean and the pixel. This process is repeated until the percentage of unchanged pixels is greater than or equals T, or the number of iterations is greater than or equals M. In each iteration, the mean of each cluster is recomputed, and these new means will be used for the next iteration. Initially, 20 classes were developed using this procedure.
The classified result was further processed to identify the surface SOC concentrations for each class based on the relationship between surface SOC concentrations and the pixel intensity values. The procedure was as follows: (i) compute the average image-intensity value and the histogram of image-intensity values for each class based on the original rectified image and the classified result. The image-intensity values of each class were extracted from the original image, whereas the boundary of each class was identified by the classified result; and (ii) determine the average surface SOC concentration and histogram of surface SOC concentrations for each class based on the relationship between surface SOC concentrations and the image-intensity values. The result was called Pre_Result2.
Based on the histogram of each class from Pre_Result2, Pre_Result2 was reclassified, and eight classes were derived from the reclassification; the result was referred to as Result2. Then according to the class range of Result2, Pre_Result1 was also classified into eight classes; and the result was referred to as Result1.
Further processing of the results (Result1 and Result2) was necessary because of two problems: The first was pixel values that were located outside the field; this area needed to be removed (for the reason of statistics and mapping). The second problem was the single-pixel classes in the results of the second classification. These single-pixel classes are mainly from spot noise in the original slide and scanning process.
The first problem was solved by measuring the field boundary with sub-meter accuracy GPS and discarding pixels outside the measured field boundary.
For the second problem, a majority algorithm was used to filter out single-pixel classes. This method sets a pixel value at location (i, j) to the pixel value that has the majority number in the filter mask. The process is as follows: (i) choose a suitable mask size and move the mask over the image. A mask with 3 by 3 cells was selected because it can effectively remove the single-pixel classes but keep all classes with five or more pixels; (ii) for each pixel at location (i, j), look for the pixel value Pm with the maximum number (majority) in the mask; and (iii) reassign the pixel value at location (i, j) to Pm. The final results were referred as Post_Result1 and Post_Result2. Comparisons between Post_Result1 and Post_Result2 were conducted by examining the area of each class and a histogram representing the degree of difference of uncommon class pixels.
The accuracy of the results obtained above was checked, based on the other 31 soil samples, which were different from those samples for model development. The check for each location was based on a point buffer (square buffer) with a buffer size of 5 by 5 pixels (10 by 10 m) to reduce the error caused from image rectification. For each location at (x, y), the buffer, with the center at (x, y), was overlaid on the result image. The average SOC concentrations within this buffer was computed, as follows:
![]() | (2) |
The measured data and the average value within the buffer were compared to check the accuracy of the final classification results, Post_Result1 and Post_Result2. Two approaches were used to check the accuracy: In the first approach, a relationship between measured and estimated values was developed by a linear regression. The r2 values were examined for the two methods. In the second approach, the measured and predicted values were classified into one of the eight classes based on the class scheme. For each location, the measured and the predicted surface SOC concentration values were examined to check if they were in the same class.
| Results and discussion |
|---|
|
|
|---|
|
![]() | (3) |
|
It could be found that a significant number of single-pixel classes existed from the above results. These single-pixel classes may show some details about the distribution of surface SOC concentrations. However, the field survey found that most of these details were not the true representation of the field distribution of surface SOC concentrations. In addition, these single-pixel classes caused too much image variance for analyzing and mapping the results. These single-pixel classes needed to be removed, which was done using the majority method. The results were then converted into vector format (classes are represented by polygons rather than by pixels), and a color scheme was applied to them for the output. Figures 5 and 6 show the results by using the first method and the second method respectively, with single-pixel classes removed.
|
|
|
|
0) were not significant at P < 0.15 for Post_Result1 and at P < 0.45 for Post_Result2 under a 0.95 confidence level, so they were not considered in the linear equations. From the scatter plots, we also noted that the prediction of the surface SOC concentrations <1.3% was better than that of the surface SOC concentrations >1.3%.
|
|
|
In addition, the methods developed in this research would have other advantages compared with grid sampling. At present, grid sampling for precision farming is labor-intensive and expensive both for soil sampling and for analysis. For example, 280 samples, based on a 0.405-ha (1-acre) grid size, would be taken in this field if the grid sampling method is used. For the method described in this paper, the number of samples (for developing the relationship between surface SOC concentrations and image-intensity values) was reduced to 28, which would be 10% of the number required to grid sample at a scale of 0.405 ha.
| Summary |
|---|
|
|
|---|
For further refinement, there are two things we may need to consider for use in other fields in the region. The first is the effect of noise from other soil properties, such as the soil Fe concentration. However, Fe concentration was as high as 1.2% in the original data and 1.1% in the test data and appeared to create no problems. Therefore, we concluded that the effect of Fe was not significant for predicting organic-C concentrations in this field. However, the effect of Fe concentrations may need to be considered in other fields because of its importance, either as variables in regression analysis or removing the effect of Fe concentration from the original image before the organic-C concentration analysis.
Another issue we may need to consider is the consistency of image-intensity values. In this study, we used a color slide as the source image. The lighting conditions at the time a photograph is taken and the development processing of slide film both may vary from one field to another and the scanning of color slides may change the true radiated values and may introduce noise. To avoid or minimize these sources of errors, a digital multispectral image could be used as the source image of a field. With digital imagery, the spectral radiance from surface objects is directly recorded and saved in image files and does not require a processing step that can vary from one image capture to another. Data from digital imagery could be more consistent and avoid the effects of variable lighting and processing associated with film development and scanning.
Received for publication October 29, 1998.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
F. Chen, D. E. Kissel, L. T. West, W. Adkins, D. Rickman, and J. C. Luvall Mapping Soil Organic Carbon Concentration for Multiple Fields with Image Similarity Analysis Soil Sci. Soc. Am. J., January 11, 2008; 72(1): 186 - 193. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Yadav and G. Malanson Progress in soil organic matter research: litter decomposition, modelling, monitoring and sequestration Progress in Physical Geography, April 1, 2007; 31(2): 131 - 154. [Abstract] [PDF] |
||||
![]() |
L. P. D'Acqui, C. A. Santi, and F. Maselli Use of Ecosystem Information to Improve Soil Organic Carbon Mapping of a Mediterranean Island J. Environ. Qual., January 9, 2007; 36(1): 262 - 271. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Persson Estimating Surface Soil Moisture from Soil Color Using Image Analysis Vadose Zone J., November 11, 2005; 4(4): 1119 - 1122. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. A. Fox and R. Metla Soil Property Analysis using Principal Components Analysis, Soil Line, and Regression Models Soil Sci. Soc. Am. J., September 29, 2005; 69(6): 1782 - 1788. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Fleming, D. F. Heermann, and D. G. Westfall Evaluating Soil Color with Farmer Input and Apparent Soil Electrical Conductivity for Management Zone Delineation Agron. J., November 1, 2004; 96(6): 1581 - 1587. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. R. Weaver, D. E. Kissel, F. Chen, L. T. West, W. Adkins, D. Rickman, and J. C. Luvall Mapping Soil pH Buffering Capacity of Selected Fields in the Coastal Plain Soil Sci. Soc. Am. J., March 1, 2004; 68(2): 662 - 668. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. B. Jaynes, T. C. Kaspar, T. S. Colvin, and D. E. James Cluster Analysis of Spatiotemporal Corn Yield Patterns in an Iowa Field Agron. J., May 1, 2003; 95(3): 574 - 586. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. G. Mueller and F. J. Pierce Soil Carbon Maps: Enhancing Spatial Estimates with Simple Terrain Attributes at Multiple Scales Soil Sci. Soc. Am. J., January 1, 2003; 67(1): 258 - 267. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. A. Fox and G. J. Sabbagh Estimation of Soil Organic Matter from Red and Near-Infrared Remotely Sensed Data Using a Soil Line Euclidean Distance Technique Soil Sci. Soc. Am. J., November 1, 2002; 66(6): 1922 - 1929. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. K. Johnson, J. W. Doran, H. R. Duke, B. J. Wienhold, K. M. Eskridge, and J. F. Shanahan Field-Scale Electrical Conductivity Mapping for Delineating Soil Condition Soil Sci. Soc. Am. J., November 1, 2001; 65(6): 1829 - 1837. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||