Research on Coal Exploration Technology Based on Satellite Remote Sensing

Coal is the main source of energy. In China and Vietnam, coal resources are very rich, but the exploration level is relatively low. This is mainly caused by the complicated geological structure, the low efficiency, the related damage, and other bad situations. To this end, we need to make use of some advanced technologies to guarantee the resource exploration is implemented smoothly and orderly. Numerous studies show that remote sensing technology is an effective way in coal exploration and measurement. In this paper, we try to measure the distribution and reserves of open-air coal area through satellite imagery. The satellite picture of openair coal mining region in Quang Ninh Province of Vietnam was collected as the experimental data. Firstly, the ENVI software is used to eliminate satellite imagery spectral interference. Then, the image classification model is established by the improved ELM algorithm. Finally, the effectiveness of the improved ELM algorithm is verified by using MATLAB simulations. The results show that the accuracies of the testing set reach 96.5%. And it reaches 83% of the image discernment precision compared with the same image from Google.


Introduction
Remote sensing is a method of acquisition, transmission, processing, and extraction of geographic information without contacting with the surface directly.Compared with the conventional methods, remote sensing technology has features such as "multipoint," "multiband," and "temporal".At the same time, it can collect the information of wavelengths that fall well outside the visible spectrum, which enlarges the range of observation.The evolution trajectory of an objective phenomenon in the time dimension can be provided by repeated observations through remote sensing.Over the past 20 years, remote sensing technology has been used successfully in the forestry, agriculture, geological survey, and marine ecology fields.Interestingly, remote sensing technology is also an indispensable research technique in resources and environment studies.For the purposes of exploring coal mines, remote sensing technology can access a wide range of coal mining area information despite of the limits in the geographical conditions.This technology has been found applicable in coal exploration and measurement.For instance, Guan [1] proposed that airborne remote sensing data can be applied to coal forecast by the experimental study of space remote sensing in Taiyuan Coal.Zou [2,3] studied the relationship between annular image, gravity, and magnetic anomalies on the satellite images of Hunan Xuefeng Coal.They found coal in Hunan Cenozoic cap by the use of the annular image.Later, based on multiple remote sensing messages for Landsat data and SPOT-5 data and airborne remote sensing data and other information, Mularz [4] explored the Belchatow lignite opencast mining area comprehensively and extracted change data of mining region by fusion of SPOT-5 panchromatic in different time and TM multispectral image eventually.Also based on hyperspectral remote sensing data, interfered radar data, and GPS positioning measurement techniques, Cloutis [5] detected the effect of environment caused by exploitation coal mine underground in Ruhrgebiet area of German.Liu et al. [6] found that the geological phenomena under the cover layer influence the soil layer or surface feature (such as vegetation) directly or indirectly via the study of TM image of Huaibei Coal, which can be used to find buried geological coal.Recently, Tan et al. [7] presented a new multiband, multipolarization, and multiangle method to analyze remote sensing data and explore the geology of deep-level coal information.So far, a large number of studies have shown that remote sensing technology is an effective way to search for coal resource.
Pankiewicz technique is often used to categorize variety of images in the area of remote sensing satellite imagery classification.Generally speaking, these categorized methods can be divided into unsupervised classification, supervised classification, and neural networks.The nearest neighbor algorithm, histogram method, and clustering method [8][9][10] are the most common unsupervised classification methods.The supervised classification methods mainly include the nearest neighbor algorithm, maximum likelihood estimation, and support vector machine [11][12][13].Artificial neural network based satellite image classification methods can be roughly divided into the PNN (Probability Neural Network), BP (Backpropagation) neural networks, and SOM (Self-Organized Feature Map) [14,15].
As is known, the variety of remote sensors is increasing and the image resolution of remote sensing is enhancing.However, the appropriate treatment, development, and application of the analytical tools are lagging behind.The traditional classification algorithm for the satellite image classification is likely to cause large scale and local minimum problems.Consequently, the speed and classification accuracy is far away from the demand.Therefore, to develop an accurate and fast automatic classification method in satellite imaging has always been a hot topic in the field of remote sensing.
To this end, this paper introduces a new single hidden layer feedforward extreme learning machine (ELM) and applies the algorithm to construct classifier for remote sensing satellite image classification.This novel algorithm is able to approach any complicated continuous functions and to learn new things with fast training speed and high accuracy.On the other hand, remote sensing technology has been widely applied to mining areas since 1970 [16][17][18][19][20][21].However, most of these methods are subjective and timeconsuming because the spectral characteristics are unknown.In order to improve the accuracy of classification, the same experiment has to be conducted repeatedly.Instead the spectral characteristics may provide more comprehensive information.Based on the satellite image-spectral characteristics, we propose a method to search and measure coal.First, we use the ENVI software to eliminate interference of spectral image data.Then, we propose an improved ELM algorithm to establish coal image classification model.Finally, we use coal image classification model to classify satellite images of the Quang Ninh coal region in Vietnam.Simulation results show that the proposed algorithm has merits over many traditional neural networks.

Remote Sensing Reflectance Spectrum Data Preprocessing
First, we downloaded the remote sensing images of coal mining of Vietnam Quang Ninh Province from USGS remote sensing image database.Landsat 5 TM is selected as the satellite species which is the fifth in the US Landsat series and is an optical Earth observation satellite.Its payloads are thematic mapper (TM) and multispectral imager (MSS).The image acquired by Landsat-5 satellite has been widely used in remote sensing of resources.The TM is divided into seven bands.The parameters of each band are shown in Table 1.
Since the sixth band is in hot infrared wavelength range, all the other six band (1st-5th and 7th) spectrums are used in this paper.
The original satellite imagery is vulnerable to the interference of atmosphere and surface reflectance of light, leading to distortion of the reflectance spectrum and geometry.Thus reflectance spectrum data must be verified before usage.The ENVI software is powerful in processing remote sensing image developed by American Exelis Visual Information Solutions, including input/output image data, scaling, image enhancement, correction, orthorectification, mosaic, information extraction, and image classification.ENVI can obtain the necessary information from remote sensing images quickly, easily, and accurately.Therefore we use ENVI to eliminate the influence of the surface reflections of the atmosphere and light.We can obtain a more realistic surface reflectance and reflectivity, surface temperature, and other physical model parameters, which can approach the actual spectral characteristics of the ground substance much better.ENVI spectral calibration consists of two steps, radiometric calibration and atmospheric correction.We open the calibration band and start the Landsat calibration tools.As shown in Figure 1, relevant parameters are set.A file must be created after six bands have been calibrated if only one band can be calibrated at a time.

Atmospheric Correction.
ENVI contains a lot of atmospheric correction models.In this paper, atmosphere is corrected by MODTRAN model of Landsat FLAASH.Since TM data is common in BSQ format during atmospheric correction data by FLAASH in BIL format, the data conversion is needed before FLAASH atmospheric correction.After the data conversion, we start FLAASH atmospheric correction.FLAASH parameter settings are shown in Figure 2. The effects before and after atmospheric correction are shown in Figures 3 and 4.
From reflectance spectrum of one point of the coal region in Figures 3 and 4, it can be seen that the coal of spectrum curve after correction is much closer to the real coal of spectrum curve.

Spectral Data Extraction.
After the spectral calibration of satellite images, we extract the spectral data samples required for the experiment using the tie point method on image.The satellite images of coal mining area in Quang Ninh Province of Vietnam are chosen for spectral sampling.As shown in Figure 4, we make a right click at one point on mine image, select "Z Profile (spectrum)," click "File," and select "Save Plot As" to output the sample data in txt file format.
Due to the fact that the spectral data used in MATLAB cannot be directly extracted from satellite images, we need to read the spectral data of satellite images by ERDAS IMAGINE

Extreme Learning Machines
The model of ELM proposed by Huang et al. [22] is constituted by input layer, single hidden layer, and output layer.The aim of ELM is to solve the above issues related to gradientbased algorithms.In ELM, input weights and single hidden layer biases are arbitrarily chosen without iterative adjust, and the only parameters to be learned in the training are the output weights which can be calculated by solving a single linear system [23].Therefore, ELM is an efficient learning algorithm that has been widely applied in regression and multiclass classification.
We first define  training samples, {, } = {  ,   }  =1 , where   ∈ R  and   ∈ R  are the th input and target vectors, and the parameters  and  are the dimensions of input and target vector.To seek a regression function from the input to the target, one popular form is the standard single hidden layer feedforward network (SLFN) [24], where  ℎ single hidden nodes fully connect the  input nodes to the  output nodes, which can be mathematically modeled as where   ∈ R  is the output vector of the th training sample,   ∈ R  is the input weight vector connecting the input nodes to the th hidden node,   is the bias of the th hidden node, and (⋅) denotes hidden nodes nonlinear piecewise continuous activation functions.The above  equations can be written compactly as where the matrix T is target matrix, The matrix H is the hidden layer output matrix, which can be randomly generated independent of the training data. = [ 1 ,  2 , . . .,   ℎ ]  (  ∈ R  ) is the output weight matrix between the hidden nodes and the output nodes.Thus, training SLFNs simply amounts to getting the solution of a linear system (2) of the output weights  [25,26].
According to the theoretical results in [27,28], ELM aims to reach the smallest training error: A simple representation of the solution of ( 4) is given explicitly by Huang et al. [24] as where H + is the Moore-Penrose generalized inverse of the hidden layer output matrix H.   hidden layer nodes of ELM neural network is 300.We change the number of ELM model, repeat the experiment 100 times, and average the results as the final result, which is shown in Table 3.

Classification Model of Coal Satellite Image Based on Improved ELM Algorithm
As shown in Table 3 and Figure 6, the accuracy of the test set by the integrated ELM model improves to more than 94%; the higher the number of models, the higher the accuracy of the model, but the longer the training time; more precisely, the training time of the 11-model is about equal to the 5model, but the accuracy of the test is about 1.4% higher than that of the latter; the accuracy of test set of the 11-model is almost equal to the 101-model while the former is much faster than the latter.Relatively speaking, the 11-model may be the most reasonable parameter.

Model Based on Stochastic Optimization Integrated ELM.
Because ELM weights and biases are random assigned, we further improve the integrated model.The design idea is summarized as follows: (i) Step one, train the individual model 100 times, compare the 100 results, and then save the  7, respectively, show the comparison results among the models of stochastic optimization integration ELM algorithm, basal ELM, integration ELM and traditional Vector Machine Support (SVM), and Backpropagation (BP).
From the numerical results, we conclude that stochastic optimization integrated XS-ELM has improved the accuracy of the test up to 96.55% which is 4% higher than the basic ELM and 1% than integrated ELM model.The accuracy of the proposed XS-ELM algorithm is better than SVM and BP neural network method and even reaches 96.5%.Its only drawback is that the training time is slightly longer.The average accuracy of training The average accuracy of testing

Application of Improved ELM Algorithm in Remote Sensing for Coal
The coal resource in Vietnam is rich in variety and good in quality.Vietnam is one of the most important countries for coal producing; the coal mines are mainly located in the northern region.Quang Ninh Province is known as the base of Vietnamese coal production, with a large open anthracite and a 100-year mining history till now.The area of coal reaches about 220 square kilometers with the volume about 3000 million tons in Quang Ninh Province known as one of the best anthracites in the world.We preprocess the remote sensing image data by ENVI.Then, we use the coal identification model to identify the remote sensing image.The identification results are shown in Figures 8(a), 9(a), 10(a), 11(a), and 12(a).Finally, we compare the coal mine image and the image of Google map which is acquired from Google HD map in accurate latitude and longitude.Then we circle out the areas of coal mining and make a careful revision.At last we conclude that the accuracy of the image of coal mine area is 100%.We compare the two images, select Google map image as the bottom layer, and   5. Figures 8-12 show the images of the coal mining area by the use of the basic ELM, integrated ELM, XS-ELM, SVM, and BP, respectively.We compare these images with Google maps by using MATLAB software.The accuracy of each method is shown in Table 5 (the basic ELM method is 78.5679%; the integrated ELM method is 81.3777%; the XS-ELM method is 83.0831%; the SVM method is 81.9066%;        the BP method is 80.7854%).Comparison results verify that the wide applicability of the proposed algorithm.The 83% accuracy of satellite image identification shows the effectiveness of this method.At the same time, the accuracy is higher than that of SVM and BP.The experimental results also show that the improved ELM algorithm can improve the accuracy of satellite image identification from 78.5% to 83%.

Conclusion
Remote sensing technology is an important method to find coal and also plays an important role in the investigation of all kinds of resources and mineral exploration.It has a great number of advantages, wide field of view, strong macro, low cost, and high efficiency.We illustrate the practical examples of coal by using spectral information of remote sensing images.The method can be applied to the exploration, measure the distribution, and reserve of open-air coal mining area.In this paper, we use artificial intelligence system and remote sensing spectral characteristics to study coal resource exploration.Simulation results verify the merits of the proposed improved ELM algorithm over some other conventional methods.With the development of remote sensing technology, it is expected that the proposed method could be applied to other fields in future, such as geology, geography, science, and coal sedimentology.

2. 1 .
Radiometric Calibration.Radiometric calibration is a process of quantifying the value of voltage recorded by the sensor or converting digital (DN) into absolute radiance value (reflectivity).ENVI provides tools specifically for radiometric calibration of Landsat satellite images, which can convert the DN value of TM into the value of radiation brightness or apparent reflectivity of atmospheric.

Figure 5 :
Figure 5: Effect of accuracy of training of different hidden layer nodes.

Figure 7 :
Figure 7: The accuracy of the five methods.

Figure 8 :
Figure 8: Images of identification of basic ELM and Google map.
(a) Identification image of integrated ELM (b) Google map image of coal of Quang Ninh, Vietnam

Figure 9 :
Figure 9: Images of identification of integrated ELM and Google map.
(a) Identification image of XS-ELM (b) Google map image of coal of Quang Ninh, Vietnam

Figure 10 :
Figure 10: Images of identification of XS-ELM and Google map.
(a) Identification image of SVM (b) Google map image of coal of Quang Ninh, Vietnam

Figure 11 :
Figure 11: Images of identification of SVM and Google map.
(a) Identification image of BP (b) Google map image of coal of Quang Ninh, Vietnam

Figure 12 :
Figure 12: Images of identification of BP and Google map.

Table 2 :
Effect of accuracy of ELM neural network of different hidden layer nodes.
Simulation results are shown in Table2and Figure5.It is seen that the accuracy of the training set and the test set is about 70% in the case that the nodes of hidden layer are relatively short; the accuracy is high and stable while the nodes of hidden layer are in 300 range and it can reach 92.9%; if the number of nodes in the hidden layer is improved again, the accuracy of the training set is improved, but the accuracy of the test set is decreased; and the fact that training time is short while nodes in the hidden layer are relatively small.
4.1.Model of Basic ELM.follows: if the weights and biases of the ELM neural network are randomly given, the results are different from each output.We integrate several ELM models and save every test set of each model.Then a final vote is made to determine the final result.This procedure can be named as "integrated ELM classification algorithm."In the experiment, the number of

Table 3 :
The effect of accuracy based on different number of ELM model.
Step two, repeat Step one 11 times and find out the weights and biases of 11 best ELM neural network; (iii) Step three, we integrate the 11 individual ELM neural network models to vote.The model has two identifications.They are labeled as Class 1 and Class 2. If there are six or more results in Class 1 in the 11 results, we consider the result as Class 1.It is similar to Class 2. The method can be called stochastic optimization integrated ELM (XS-ELM) based on its features.Table4and Figure

Table 4 :
The accuracy of the five methods.

Table 5 :
The accuracy comparing the map and the identification of different ELM methods.