Extraction of Earth Surface Texture Features from Multispectral Remote Sensing Data

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). -e 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 single-band texture analysis method that is based on gray level co-occurrence matrix, the proposed approach has higher classification precision and efficiency.


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.e earth surface texture is a good solution to the problem because of the stability characteristics [1].
e classical texture extraction and analytic approaches include gray level co-occurrence 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.is 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 Worldview-II 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 Co-occurrence Matrix (GLCM) is made.e experimental results indicate that the gray level difference frequency spatial has higher classification accuracy and efficiency.

Worldview-II Multispectral Remote Sensing Data
Worldview-II 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 near-infrared 1) and four characteristic bands (coastal, yellow, red edge, and near-infrared 2).e data analyzed in this paper are Worldview-II 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.

Compression of Multispectral Remote Sensing Data
Principal component analysis could project the high dimensional data onto a low-dimensional 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 Worldview-II multispectral remote sensing data into a column vector W as follows: Principal component analysis makes a combination of W through linear transformation and guarantees that P has the largest variance after transformation, as shown in the following equation: where is the m-dimensional space to be determined and R is the covariance matrix of W, thereby, the variance of P could be computed as follows: ereby, solving the maximum value of D(P) is equal to seeking the vector L that makes D(P) the largest.e length L is limited to unit length, and then the question is converted to max D(P) � L T RL, In last equation, the covariance matrix R could be expressed as follows: where where a 1 , a 2 , • • • , a n is the eigenvector corresponding to the unit orthogonal eigenvectors.Let a ∈ R n and multiplying Equation ( 5) with a T and a on the left and right side separately, we get Let k � Γ T a, then a T a � k T k; Equation (6) satisfies Equation ( 4) can be rewritten as follows: max If a � a 1 , then a T 1 Ra 1 � λ 1 , which indicates that the maximum value of L T RL is at the point of a 1 under the condition of L T L � 1, thereby the first PCA principle component could be expressed as p 1 � a T 1 W. e 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: e cumulative contribution rate of the first l principal components is as follows: e contribution rate indicates the ability that principle components reflect W. It determines the number of principal components after compression of multispectral remote sensing data.

Texture Features Extraction from
Multispectral Remote Sensing Data 4.1.Gray Level Co-Occurrence 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 two-order combined conditional probability density among image pixels [10].Assume the remote sensing image is of size M × N; the gray level is L; the distance between two pixels is d; the angle is θ; the gray levels are separately i and j; the times that these two pixels appear simultaneously is P(i, j, d, θ) which could be expressed as follows: 2

Journal of Electrical and Computer Engineering
Where m and n represent the number of remote sensing image pixels in a row and a column; i, j 1, 2, x 2 , y 2 are the pixel coordinates in the image.Figure 1 shows the spatial sketch map of GLCM.If the remote sensing image has L picture gray levels, the size of the gray level cooccurrence matrix is L × L, d 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 0 ∘ , 45 ∘ , 90 ∘ , and 135 ∘ .e element at the a-th row and b-th column in P(i, j, d, θ) 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.d is related to the image, the step is usually set as δ 1, and the central pixel to be operated and compared with the directly adjacent pixel.

Establishment of the Gray Level Di erence Frequency
Spatial.Gray level di erence frequency spatial is proposed as a texture extraction approach for the multiprinciple component, which is based on the Gray Level Co-occurrence Matrix.Assume the gray level of primary gradient remote sensing image is g; rstly, iterate the principle components' remote sensing image with a rectangular window which is of size } are separately the horizontal and vertical spatial domains of the window and G 0, 1, • • • , g − 1 is the gray level.L n × L m is the windows' pixel set with ranking sequence in the row and column, specifying I as the conversion formula of pixels in L n × L m to G: In di erent windows of principle components' image, the appearance probability of the pixel couples with distance d, direction θ, and gray level di erence Δ composing the GLDFS and it is named P, in which (m1 Pi , n1 Pi ), (m2 Pi , n2 Pi ) are coordinates of two pixels in the i − th principal component L n × L m ; I(m1 Pi , n1 Pi ), I(m2 Pi , n2 Pi ) are gray levels corresponding to the two pixels; and (m2 Pi , n2 Pi ) locates in the θ direction of (m1 Pi , n1 Pi ) with distance d. e two pixel gray level's di erence 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 d and direction θ. e 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 0 ∘ , 45 ∘ , 90 ∘ , and 135 ∘ , thereby generating four n-dimensional spaces, and # indicates the number of the elements in the space; the statistical methods are as follows: Journal of Electrical and Computer Engineering

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][13][14].is 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) e maximum similarity: e energy, entropy, autocorrelation, and correlation are the extentions of Haralick's method in high-dimensional 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.

Analysis of Experimental Results
e experiment data are based on the multispectral remote sensing data of 8 different bands, which were recorded by Worldview-II 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 Co-occurrence Matrix.
e gray level difference frequency spatial and Gray Level Co-occurrence Matrix choose worldview-II 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.

e Experimental Data Compression.
PCA is applied to compression of 8 different bands of Worldview-II 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 Worldview-II 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%.

5.2.
e 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 gray-scale changes of the adjacent pixel, the window size is 13 × 13.Features including energy, self-correlation, correlation, cluster dark, and maximum self-similarity are selected; the sample size is 1350 among which the reef, sea, and foam samples are 350 separately.e SVM pattern is applied for the purpose of dividing the data into three.e result of landmark feature classification with the gray level difference frequency spatial and Gray Level Co-occurrence Matrix are shown in Tables 1 and 2. e 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 4 Journal of Electrical and Computer Engineering identi cation 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 gured out that the gray level di erence frequency spatial is capable of dealing with three types of landmark features at the same time and has better      a better classification results with three kinds of surface features.e method of Gray Level Co-occurrence 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 near-infrared are processing, seawater is identified as foam.e method of Gray Level Co-occurrence 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 single-band 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.e extracted results are processed and used as SVM input variables of multilandmark classification.Finally, the classification time consumption is compared and analyzed.e 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.e 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.

Conclusions and Discussions
According to the Worldview-II multispectral remote sensing data, this paper proposed a texture feature method based on gray-scale difference in space-frequency.is method compresses the multispectral remote sensing data after having it disposed, carries out gray-scale statistic on the main components after the compression process, and builds the gray-scale frequency difference space model.e method is validated by taking the data collected by Worldview-II.From the results, it can be seen that the gray level difference frequency spatial can extract texture features and recognize and classify multiband, multi-landmark, and multi-spectral remote sensing data at the same time and has an advantage over tradition Gray Level Co-occurrence Matrix on consideration of both recognition efficiency and quality.

Figure 1 :
Figure 1: e space description of the GLCM.

Figure 2 :Figure 3 :
Figure 2: e contribution rate and the cumulative contribution rate of each principal component.

Figures 4 (
Figures 4(a)-4(d) are the classification result of conventional bands of the Gray Level Co-occurrence Matrix.Figures 4(e)-4(h) are the classification result of characteristic bands of the Gray Level Co-occurrence Matrix.Figure 4(i) is the classification result of the principal component data of the gray level difference frequency spatial. is paper presents a method named the gray level difference frequency spatial for the principal components data processing and has

Figure 4 :
Figure 4: e classification processing results of GLDFS in conventional bands and characteristic bands and GLCM in principal component data: (a) blue band, (b) green band, (c) red band, (d) near-infrared 1 band, (e) coastal band, (f ) yellow band, (g) infrared band, (h) nearinfrared 2 band, and (i) principal component data.

Table 1 :
e recognition rate of GLDFS and GLCM of conventional bands.

Table 2 :
e recognition rate of GLDFS and GLCM of characteristic bands.