Research Article  Open Access
Zhenxing Zhang, Feng Gao, Bin Ma, Zhiqiang Zhang, "Extraction of Earth Surface Texture Features from Multispectral Remote Sensing Data", Journal of Electrical and Computer Engineering, vol. 2018, Article ID 9684629, 9 pages, 2018. https://doi.org/10.1155/2018/9684629
Extraction of Earth Surface Texture Features from Multispectral Remote Sensing Data
Abstract
Earth surface texture features referring to as visual features of homogeneity in remote sensing images are very important to understand the relationship between surface information and surrounding environment. Remote sensing data contain rich information of earth surface texture features (image gray reflecting the spatial distribution information of texture features, for instance). Here, we propose an efficient and accurate approach to extract earth surface texture features from remote sensing data, called gray level difference frequency spatial (GLDFS). The gray level difference frequency spatial approach is designed to extract multiband remote sensing data, utilizing principle component analysis conversion to compress the multispectral information, and it establishes the gray level difference frequency spatial of principle components. In the end, the texture features are extracted using the gray level difference frequency spatial. To verify the effectiveness of this approach, several experiments are conducted and indicate that it could retain the coordination relationship among multispectral remote sensing data, and compared with the traditional singleband texture analysis method that is based on gray level cooccurrence matrix, the proposed approach has higher classification precision and efficiency.
1. Introduction
Remote sensing technology can extract high resolution regional marine environmental information in time, especially for the complex sea area. Multispectral remote sensing data reflects the interested target or regional radiation characteristics through the electromagnetic spectrum of multiband, and it has the advantages of wide range, multiphase, multiband, and high resolution. Remote sensing image could enrich the spectral characteristics of landmark and find out more detailed information, such as the structure, shape, and texture. However, in virtue of the fact that same objects possess different spectral and different objects share same spectral, the applications of remote sensing data would be serious restricted if only spectral information is taken into consideration. The earth surface texture is a good solution to the problem because of the stability characteristics [1].
The classical texture extraction and analytic approaches include gray level cooccurrence matrix method [2, 3], wavelet analysis method [4, 5], Gabor spectrum method [6], and so forth. While all these methods could only be applied to analyze the information of single band in remote sensing images, for multispectral remote sensing data, all the bands should be processed separately, which would decrease the extraction efficiency badly.
Because of the geometric characteristics of the surface object, it has a unique texture features on the remote sensing images. So, the different surface objects can be extracted through the texture features. This paper utilizes gray level difference frequency spatial to extract texture features of multiband remote sensing data. We firstly conduct principal component analysis (PCA) on the eight bands of WorldviewII multispectral images and compress these data on basis of guaranteeing against loss of spectral information. Make gray difference statistics on the compressed principle components and establish the gray level difference frequency spatial. In the experiments, the gray level difference frequency spatial is used to extract texture features, and a comparison with Gray Level Cooccurrence Matrix (GLCM) is made. The experimental results indicate that the gray level difference frequency spatial has higher classification accuracy and efficiency.
2. WorldviewII Multispectral Remote Sensing Data
WorldviewII is one of the highest resolution remote sensing satellites, and it has the highest spatial resolution (0.46 m in the panchromatic band and 1.84 m in the multispectral bands). It provides high resolution multispectral data with eight bands, which include four conventional bands (red, green, blue, and nearinfrared 1) and four characteristic bands (coastal, yellow, red edge, and nearinfrared 2). The data analyzed in this paper are WorldviewII multispectral remote sensing image of the Sea Islands; the texture features of eight bands are extracted. Firstly, we calibrate the data and get the radiance data. Secondly, atmospheric correction is conducted to eliminate the influence of atmosphere and illumination, and the actual reflectance of surface objects is obtained. Finally, we make orthorectification on the data through a few control points, thus eliminating the geometric distortion.
3. Compression of Multispectral Remote Sensing Data
Principal component analysis could project the high dimensional data onto a lowdimensional space. It takes the variance in size as the evaluation standard of information quantity; the greater the variance, the more information it provides [7, 8]. On the premise of keeping useful information of multispectral remote sensing data, principal component analysis could reduce the correlation and redundant information in order to compress multispectral remote sensing data. We transform WorldviewII multispectral remote sensing data into a column vector as follows:
Principal component analysis makes a combination of through linear transformation and guarantees that has the largest variance after transformation, as shown in the following equation:where is the mdimensional space to be determined and is the covariance matrix of , thereby, the variance of could be computed as follows:
Thereby, solving the maximum value of is equal to seeking the vector that makes the largest. The length is limited to unit length, and then the question is converted to
In last equation, the covariance matrix could be expressed as follows:where is the characteristic value of and is satisfied. , where is the eigenvector corresponding to the unit orthogonal eigenvectors. Let and multiplying Equation (5) with and on the left and right side separately, we get
Let , then ; Equation (6) satisfies
Equation (4) can be rewritten as follows:
If , then , which indicates that the maximum value of is at the point of under the condition of , thereby the first PCA principle component could be expressed as . The contribution rate reflects the information quantity contained in each principle component, and the contribution rate of the ith principal component could be computed as follows:
The cumulative contribution rate of the first principal components is as follows:
The contribution rate indicates the ability that principle components reflect . It determines the number of principal components after compression of multispectral remote sensing data.
4. Texture Features Extraction from Multispectral Remote Sensing Data
4.1. Gray Level CoOccurrence Matrix
Gray Level Cooccurrence Matrix is the most direct and simplest texture analysis approach, which considers the spatial structure of remote sensing images [9]. It describes the image texture through the twoorder combined conditional probability density among image pixels [10]. Assume the remote sensing image is of size ; the gray level is ; the distance between two pixels is ; the angle is ; the gray levels are separately and ; the times that these two pixels appear simultaneously is which could be expressed as follows:Where and represent the number of remote sensing image pixels in a row and a column; , , , , are the pixel coordinates in the image. Figure 1 shows the spatial sketch map of GLCM. If the remote sensing image has picture gray levels, the size of the gray level cooccurrence matrix is , represents the distance of two pixels in the remote sensing image, represents the angle between the connection line of the two pixels and horizontal direction, and it is usually set as , , , and . The element at the ath row and bth column in represents the appearance times of all the pixel couples that are apart from each other in the direction, with gray values i and j, respectively. is related to the image, the step is usually set as , and the central pixel to be operated and compared with the directly adjacent pixel.
4.2. Establishment of the Gray Level Difference Frequency Spatial
Gray level difference frequency spatial is proposed as a texture extraction approach for the multiprinciple component, which is based on the Gray Level Cooccurrence Matrix. Assume the gray level of primary gradient remote sensing image is ; firstly, iterate the principle components’ remote sensing image with a rectangular window which is of size . Assume and are separately the horizontal and vertical spatial domains of the window and is the gray level. is the windows’ pixel set with ranking sequence in the row and column, specifying as the conversion formula of pixels in to G:
In different windows of principle components’ image, the appearance probability of the pixel couples with distance , direction , and gray level difference composing the GLDFS and it is named , in which , are coordinates of two pixels in the principal component ; , are gray levels corresponding to the two pixels; and locates in the direction of with distance . The two pixel gray level’s difference is expressed as follows:where . Make statistics of the gray level difference in each principal component and map the result to kdimensional space according to the distance and direction . The coordinates of the spatial points are , and k is the number of the principle components’ remote sensing images. Meanwhile, make statistics for each principal component at the four directions , , , and , thereby generating four ndimensional spaces, and indicates the number of the elements in the space; the statistical methods are as follows:
4.3. Texture Features’ Description of the Gray Level Difference Frequency Spatial
Haralick proposed Gray Level Cooccurrence Matrix 14 properties to describe the image’s texture features, and Sor proposed Gray Level Cooccurrence Matrix 10 properties aiming at the SAR images [11]. But some properties are relevant, resulting in information redundancy; besides, part of properties is not suitable for analysis of remote sensing images’ texture features [12–14]. This paper makes analysis of seven texture features with GLDFS and could be computed as follows:(1)Energy:(2)Entropy:(3)Autocorrelation:(4)Correlation:(5)Highlight degree of clustering:(6)Dark degree of clustering:(7)The maximum similarity:
The energy, entropy, autocorrelation, and correlation are the extentions of Haralick's method in highdimensional space, and the maximum similarity is generalization of Soh's method, the highlight degree, and the dark degree of clustering, which are the simulation calculations of human perception.
5. Analysis of Experimental Results
The experiment data are based on the multispectral remote sensing data of 8 different bands, which were recorded by WorldviewII on 8th April, 2009. In order to ensure the same landform features, the reef selected by the experiment is uncovered with plants consisting of sedimentary and metamorphic rocks. Gray value variance has a regular pattern when counted on space and texture feature shows up with strong intensify as well, so this paper tests and verifies the gray level difference frequency spatial texture feature extracting method through the classification of landform and compares it with the Gray Level Cooccurrence Matrix. The gray level difference frequency spatial and Gray Level Cooccurrence Matrix choose worldviewII remote sensing data of the same region which cohere well, and thus the comparability of the experiment data is assured. Experimental classification uses SVM, which is based on structural risk minimization principle. SVM is a better solution to small samples and nonlinear problems.
5.1. The Experimental Data Compression
PCA is applied to compression of 8 different bands of WorldviewII multispectral remote sensing data. In this paper, the first three principal components are derived as the feature of the gray level difference frequency spatial texture feature. Figure 2 shows the contribution rate as well as the cumulative contribution rate of each main component of the WorldviewII 8band multispectral data after principal component analysis transform. Figure 3 shows the three principal component image data after principal component analysis transformation and compression. Figure 3(a) shows the first principal component, and the contribution rate is 74.12%. Figure 3(b) shows the second principal component, and the contribution rate is 22.16%. Figure 3(c) shows the third principal component, and the contribution rate is 2.73%.
(a)
(b)
(c)
5.2. The Classification and Analysis of Sample Data Experiment
Deriving the texture features of the three principal components after being compressed with the gray level difference frequency spatial and deriving the texture features of the four regular bands and four special bands which are compared separately with Gray Level Cooccurrence Matrix, the distance d of both methods is 1; merely comparing the grayscale changes of the adjacent pixel, the window size is . Features including energy, selfcorrelation, correlation, cluster dark, and maximum selfsimilarity are selected; the sample size is 1350 among which the reef, sea, and foam samples are 350 separately. The SVM pattern is applied for the purpose of dividing the data into three. The result of landmark feature classification with the gray level difference frequency spatial and Gray Level Cooccurrence Matrix are shown in Tables 1 and 2. The recognition rate is based on the typical testing samples randomly selected from the samples including reef, sea, and foam data. Calculating the percentage of correct identification number for each landmark class with its corresponding total samples, the average value is the right recognition number of the three types of landmarks with the gross sample number. It can be figured out that the gray level difference frequency spatial is capable of dealing with three types of landmark features at the same time and has better recognition efficiency and quality than traditional Gray Level Cooccurrence Matrix.


The classification using the gray level difference frequency spatial and the classification using Gray Level Cooccurrence Matrix in conventional bands and characteristic bands of processing results about principal component data can be seen from Figure 4. In Figure 4, black indicates the reef, dark gray indicates foam, and light gray indicates sea. Figures 4(a)–4(d) are the classification result of conventional bands of the Gray Level Cooccurrence Matrix. Figures 4(e)–4(h) are the classification result of characteristic bands of the Gray Level Cooccurrence Matrix. Figure 4(i) is the classification result of the principal component data of the gray level difference frequency spatial. This paper presents a method named the gray level difference frequency spatial for the principal components data processing and has a better classification results with three kinds of surface features. The method of Gray Level Cooccurrence Matrix has poor identification with reef when it is processing blue bands and coastal bands, and some island and reef are identified as foam mistakenly; when green bands and yellow bands is processing, the identification to seawater is poor, some seawater is identified as reef and island mistakenly; when red band and red edge band are processing, the identification to seawater and foam is poor, some seawater and foam is identified as reef and island mistakenly; when No. 1 and No. 2 bands of nearinfrared are processing, seawater is identified as foam. The method of Gray Level Cooccurrence Matrix has high recognition to some landmark, but for the multilandmarks of multispectral remote sensing data, it cannot be guaranteed to have higher recognition rate in classification. If a variety of surface landmarks is classified, the singleband data need to be processed separately. In the experiment, texture features are extracted from the conventional band and special band data using Gray Level Co‐occurrence Matrix method and from the principal component data using gray level difference frequency spatial method. The extracted results are processed and used as SVM input variables of multilandmark classification. Finally, the classification time consumption is compared and analyzed. The efficiency of texture feature extracted by two methods using SVM classification is as shown in Table 3, and the results show that the gray level difference frequency spatial to the main component of compressed data has better efficiency in ensuring a higher recognition rate, but takes long time. The main reason is that the conventional band and special bands contain four bands of data; Gray Level Cooccurrence Matrix requires processing individual band and makes the results of 4 bands as SVM input variables to multilandmark classification. Due to the increased dimension of the input variables, the efficiency of calculation of SVM is reduced and the classification time becomes longer.
(a)
(b)
(c)
(d)
(e)
(f)
(g)
(h)
(i)

6. Conclusions and Discussions
According to the WorldviewII multispectral remote sensing data, this paper proposed a texture feature method based on grayscale difference in spacefrequency. This method compresses the multispectral remote sensing data after having it disposed, carries out grayscale statistic on the main components after the compression process, and builds the grayscale frequency difference space model. The method is validated by taking the data collected by WorldviewII. From the results, it can be seen that the gray level difference frequency spatial can extract texture features and recognize and classify multiband, multilandmark, and multispectral remote sensing data at the same time and has an advantage over tradition Gray Level Cooccurrence Matrix on consideration of both recognition efficiency and quality.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was funded by the National Key Research and Development Project of China (2017YFC1404100, 2017YFC1404102).
References
 J. Feng and N. Shu, “A novel texture feature extraction of hyperspectral,” Journal of WuHan University of Technology, vol. 31, no. 3, pp. 11–17, 2009. View at: Google Scholar
 R. M. Haralick, K. Shanmugam, and I. H. Dinstein, “Textural features for image classification,” IEEE Transactions on Systems, Man, and Cybernetics, vol. 3, no. 6, pp. 610–621, 1973. View at: Publisher Site  Google Scholar
 G. N. Srinivasa and G. Shobba, “Statstical texture analysis,” Proceedings of World Academy of Science, Engineering and Technology, vol. 11, pp. 196–201, 2006. View at: Google Scholar
 C. S. Lu, P. C. Chung, and C. F. Chen, “Unsupervised texture segmentation via wavelet transform,” Pattern recognition, vol. 30, no. 5, pp. 729–742, 1997. View at: Publisher Site  Google Scholar
 M. K. Bashar, T. Matsumoto, and N. Ohnishi, “Wavelet transformbased locally orderless images for texture segmentation,” Pattern Recognition Letters, vol. 24, no. 15, pp. 2633–2650, 2003. View at: Publisher Site  Google Scholar
 A. K. Jain and F. Farshid, “Unsupervised texture segmentation using Gabor Filters,” Pattern recognition, vol. 24, no. 12, pp. 1167–1186, 1991. View at: Publisher Site  Google Scholar
 A. Kaarna, “Integer PCA and wavelet transforms for multispectral image compression,” in Proceedings of International Geoscience and Remote Sensing Symposium, pp. 1853–1855, Sydney, Australia, July 2001. View at: Google Scholar
 C. Jutten and J. Herault, “Independent component analysis versus PCA,” in Proceedings of European Signal Processing Conference, pp. 643–646, Grenoble, France, September 1988. View at: Google Scholar
 G. M. Xian, “An identification method of malignant and benign liver tumors from ultrasonography based on GLCM texture features and fuzzy SVM,” Expert Systems with Applications, vol. 37, no. 10, pp. 6737–6741, 2010. View at: Publisher Site  Google Scholar
 R. F. Walker, P. T. Jackeay, and I. D. Longstaff, “Recent developments in the use of the Cooccurrence matrix for texture recognition,” in Proceedings of International Conference on Digital Image Processing, pp. 463–471, Santorini, Greece, July 1997. View at: Google Scholar
 L. K. Soh, “Texture analysis of SAR sea ice imagery using gray level Cooccurrence matrices,” IEEE Transactions on Geoscience and Remote Sensing, vol. 37, no. 2, pp. 780–795, 1999. View at: Publisher Site  Google Scholar
 R. Jobanputra and D. A. Clausi, “Preserving boundaries for image texture segmentation using gray level cooccurring probabilities,” Pattern Recognition, vol. 39, no. 2, pp. 234–245, 2006. View at: Publisher Site  Google Scholar
 D. A. Clausi, “An analysis of cooccurrence statistics as a function of gray level quantization,” Canadian Journal of Remote Sensing, vol. 28, no. 1, pp. 45–62, 2002. View at: Publisher Site  Google Scholar
 P. Philippe, O. Richard, and C. Dubes, “Performance evaluation for four classes of textural features,” Pattern Recognition, vol. 25, no. 8, pp. 819–833, 1992. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Zhenxing Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.