Advanced Investigation of Remote Sensing to Geological Mapping of Zefreh Region in Central Iran

This study has tried to prove the ability of remote sensing techniques to extract information necessary for preparation of geological mapping of the earth’s surface using multi-spectral satellite images which are rich sources of Earth’s surface information. In this study, the surface geological mappings of Zefreh region have been investigated through ASTER, OLI, and IRS-PAN remote sensing data. To prepare the geological map, preprocessing steps and reducing noises from data using MNF algorithm were firstly carried out. Then a set of processing algorithms and image classification methods are included; the band rationing, color composite and pixel classification based on maximum likelihood, spectral and sub-pixel classification methods of spectral angle mapper (SAM), spectral feature fitting (SFF), linear spectral differentiation (LSU), hill-shade images and automatic lineament extraction were used. Confusion matrix was formed for all classified images through control points were randomly selected from 1:25,000 map of the region to determine the accuracy of obtained results, which indicated the maximum accuracy (up to 90%) of output images. Comparing the results obtained from these methods with the map prepared by ground operations confirmed accuracy results. Finally, the surface geology and fault map of Zafreh region was produced by combining detected geological formations and tectonic lineaments.


Introduction
Zefreh region is located on the Uremia-Dokhtar volcanic belt.Geology of the re- gion includes consists of the igneous, sedimentary and metamorphic.Various types of igneous rocks including tonalite, gabbro and diorite have infiltrated to the region during cutting-bending development of Qum-Zefreh fault system.The adjacent sedimentary rocks with intrusive masses has caused metamorphism in the region [1] Figure 1.This article studying on an area with complex geology to show the ability of remote sensing Images in creating Geological maps with good overall accuracy in comparing with the result of 1:25,000 map created with filed operations lately.To reach this goal a collection of remote sensing data with different spectral and spatial resolutions was processed in an Image processing algorithm complex in the form of flowchart Figure 2. The traditional ways of preparing geological maps are very time consuming and costly.In the event that remote sensing provides a variety of satellite images with different spatial and spectral resolutions with fast and easy access which the matter has increased the ability of this technology to detect accidents and natural disasters and planning for large-scale actions.Phenomena selectivity absorbs and reflect electromagnetic energy differently based on their different surface molecular  compound.Texture, pattern, color tone and geological concept are among the most useful features for identifying lithological units in remote sensing images.All of these features lead to creation of different spectral absorption characteristics (the spectral signature curve) [2].The brightness of each pixel is expressed as a number and its position inside the image is identified based on its row and column number.A cell includes one or several dark and light phenomena based on the area which it covers.Only one digit is recorded for a pixel which is the brightness grade average of total of the phenomena.As a result, details inside the pixels are not expressed and pixels are considered pure and hence, any pixel belongs to a certain class.It is while that in sub-pixel classification methods, not only a pixel can be placed in more than one class but also there is no obligation for the number of classes to place combined pixels in them.The matter has led that a set of processing and classification methods of basic pixel, spectral and sub-pixel is used to extract maximum information about available images [3].
The study region is one of the rural district of Isfahan city which is located at 75 kilometers distance of North East of Isfahan with 52 degrees and 15 minutes of longitude and 33 degrees latitude in the upper left point and 52 degrees and 23 minutes and 30 seconds of longitude and 45 minutes of longitude in the lower right point.Zefreh village located in a mountainous area and its height from sea level is 2070 meters.

Materials and Methods
In present study, the images of ASTER, LANDSAT8-OLI sensors and IRS single panchromatic band sensor have been used for create a geological map of the study area.ASTER sensor covers a wide spectral range with high spatial resolution by 14 bands from visible wave length to thermal infrared.In present study, B1 data of ASTER sensor has been used which is Geo referenced, but has CROSSTALK error.The error is created due to the light scattering by detector of Band 4 which the impacts of this noise are removed from the image using Crosstalk Software (ccs).To remove atmospheric effects from the used data, FLAASH Atmospheric Correction was applied on the data using ENVI software [4].LANDSAT8 satellite is equipped with two sensors (OLI operational land imaging sensor and TIRS infrared thermal sensor).These two sensors gather data on 9 short-wave bands and 2 thermal wavelength bands.The bands of OLI sensor of this satellite were used in present study.Radiometric and atmospheric corrections were applied on the data using ENVI 5.02 Software.The used image was Geo referenced due to utilizing level T1 of LANDSAT8 products and there was no need to geometric correction [5].
The IRS-PAN image with spatial resolution of 5.7 meter was also used as high-resolution spatial data for creating pan sharp images and for linement extraction operation.
The general schema of operation set performed on the input data that includes the preprocessing steps, the image processing methods, the accuracy calculation of the classified images and merging output images for creating final geological map, are shown in the flowchart of Figure 2. The Images processing algorithms used in this study have been selected based on the spatial and spatial power of the input data.In continue, the performance of each of these algorithms has been discussed.

Image Processing Algorithms
Minimum Noise Fraction Algorithm (MNF algorithm) is based on two consecutive operations [6].The first step is based on estimating the level of noise in data and shown by a correlation matrix.The second step is related to calculations for measuring correlation in data.In this step, several data packages are created including weighted data on standard deviation in raw data.Since all bands are participated in weighting each of the produced data package, the algorithm keeps the general information of data.Most of data related to the difference in earth surface effects are saved in the first data packages and the last packages include noises [7].The performance of Principal component analysis algorithm (PCA algorithm) defines satellite data on a new axis based on their frequency.The main component of this analysis is data integration with a high variance in the first few components and the next components mostly include noises [8].
The images of band ratio are used to display spectral contrast of certain spectral aspects.Therefore, the method can be used to separate rock units from each other.Another advantage of images' separation is to remove or reduce adverse effects resulting from exposure changes such as topography effects [9] [10].
The best reference to create a False Color Composite is spectral signature curve of targets, but as much as the number of sensing bands can help.Since non-visible bands are also used in color combination, the natural color of objects are not observed in produced image, but it would be easier to detect the effects in this status of real scheme by knowing the exact arrangement of assignment bands and how behave the effects in these bands [11].The method of determining the optimum index factor is used to select the best Tri-band color combination.The coefficient of optimum index depends on to two parameters of variance and correlation coefficient of bands in the band combination.
The main advantage of maximum likelihood classification method is to take into account changes of each class in classifying the considered pixel using covariance matrix.The method relies on the assumption that data of each used band as input of algorithm have normal distribution and it is only need to select a number of pixels to correctly estimate the mean vector and covariance-variance matrix [12].The technique classifies each pixel of satellite image only in a single class based on its features [13].
A pure pixel is one which covers a region with the same material on a wider area of a sensing pixel's size and provides a spectral profile of a certain event (effect) the Earth surface [6].Pixel purity index (PPI) image is formed in a way that the value of each pixel is equal to the number of times it has been sharply recorded.In fact, the algorithm separates purer pixels from mixed ones.The PPI technique is applied on the output of MNF algorithm in which bands with noise have been removed [14].
In Spectral angle mapper (SAM) method, each spectrum is considered corresponding to the vector in a space with dimensions equal to the number of bands.This method is not sensitive to the brightness of the pixel and the algorithm of indicating spectral angle calculates the similarity between a spectrum of pure member (T) and a reference spectrum (R) for N bands of Arc sinus point index of two spectrums as like equation of Equation ( 1) [15].
Therefore its output for each pixel is the angular distance between two spectrums which expressed in radians and varies from 0 to Π/2.Therefore the spectrum with lowest angle to the pure member would be the most similar spectrum [16].
The narrowest spectral range including the total of absorption features are used to obtain most accurate results from SFF classification.It would be possible to obtain the best accordance (fitness) with the minimum error by separating an absorption feature.To conduct SFF, the reference spectrum must be firstly extracted form an image or a spectral library and then the statistical sequence of data must be removed [17].Therefore, the results of SFF include a scaling image and a RMS image or a combined fit image (scale/RMS) to each reference spectrum.In fact, the image is a sensing from the depth of absorption feature associated with the frequency of considered materials.
Linear spectral unmixing algorithm (LSU) considers the reflected energy from each pixel as a linear combination of each of the various pure coves' reflections on the surface of that pixel.In fact, the reflections created by each of the pure covers on the surface of each mixed pixel are participated in the linear combination with the weights which are their mixing ratio [18].The mixing occurs when materials with different spectral features presented in a pixel of image [19].The algorithm finds the relative or absolute frequency of a number of spectral combinations or end members.The pure spectral result is a series of gray scale images for each end member and the error image of Root Mean Square [18] [20].
Hill shade images are prepared using DEM or images with high spatial resolution and change in the angle of azimuth and altitude of the sun [21] [22] [23].The technique is carried out through elevation data and highlighting it and changes in the radiation angle of the sunlight to the Earth' surface which makes lineaments visible in the image.Combining images at the pixel level is the most common level of combining images used to improve image quality to better visual interpretation.In combining images, multi-spectrum images with medium spatial resolution are mostly combined with an image with higher spatial resolution and provide a proper combination to visual interpretation.The resulting color combination involving more spatial details compare to the initial combination of multi-spectrum image.

Results
The main purpose of present study was to prepare geological mapping through combine the results obtained from algorithms of images' processing and classification.Firstly, the noises of ASTER and LANDSAT8 images were calculated and removed using MNF algorithm.The output of this algorithm was an image with 9 components for ASTER image and 7 components for LANDSAT8 image.The amounts of eigenvalue were studied to select components with higher level of variance among data packages of desired images Table 1.Components 8 and 9 in ASTER-MNF image and components 6 and 7 in LANDSAT8 were removed due to having the lowest value of variance.Figure 3 Represents the spectral diagram of sensing bands of ASTER and LANDSAT8 before and after implementing MNF algorithm.
Principal component Analyze algorithm (PCA) was implemented to create data packages with the least amount of duplicate information.The components were used to produce false color composite to create a lithology separation among data during image analysis.
The images of band ratio were prepared through dividing sensing bands of ASTER using bands of the image with special spectral signature of considered rock units which had been prepared to create separation Table 2.The band ratios (B6 + B8)/B7) were used to separation of carbonate and dolomite units.The The false color composite technique was also applied on remote sensing data using band ratios to create lithology separation in the region Table 3.The OIF index was applied on the images of LANDSAT8 and ASTER to select the best falsecolor composite which its results have been presented in Table 4 and Table 5.The best false color composite for ASTER data was created using band composite of RGB = 731 Figure 4, which had the highest level of frequency and lowest amount of duplicate information among color composites can be prepared from ASTER image sensing.The false color composite of Abdolsalam at the ratios of 12/14, 3/1 and 4/1was applied on the ASTER image sensing to create lithology separation to identify Dolomitic units, Ferro silicate and quartz-rich regions of the west sedimentary units of under study region Figure 5. Figure 6 represents the false color composite (RGB = 761) which has been prepared from bands of LANDSAT8 image and utilizes from the highest level of OIF index.In the Vander Mir false color composite has been prepared from band ratios of RGB = (B6/B7, B4/B2, and B3) of LANDSAT8 data, the E3PY unit is observable       and the around met volcanic units with pink.In the image, the dolomitic areas are also observable with light yellow Figure 9.
The technique of pixel purity index was applied on sensing image of ASTER with PPI ITERATIONS: 5000 and THRESHOULD VALUE: 2.5,000 using ENVI5.02Software to identify pure pixels.Then ASTER-MNF data were used in the N-dimensional visualer to classify the material of pure pixels of the image.
The extracted pure pixels were used for sampling rock units to perform classification methods of SAM, SFF and LSU.These three methods were used to classify Sedimentary units located in west part of study region in 4 different classes to have a proper criterion of comprising classification accuracy (Figure 10).
LANDSAT8 image was used to perform maximum likelihood classification method.Firstly the training samples were prepared at ARC GIS10.2 environment to select classes from map of the study area and available satellite images of Zefreh region which were overlaid and geo-referenced.The samples were clus-    tered in 8 classes.Then the algorithm of maximum likelihood classification was applied on the image using ENVI 5.02 Software.Then, MAJOORITY filter was applied on the classified image to create more uniformity in produced classes Figure 11.THE output classified images produced by SAM, SFF, LSU and maximum likelihood classification methods, used for calculating Confusion to determine the accuracy of classified images.A number of control points were determined to each class which each of the classes of classified images was comprised with them.The obtained results indicate the maximum overall accuracy.The accuracy level of 95% and Kappa coefficient of 0.943 were obtained for maximum likelihood basic pixel classification.Also, the overall accuracy of 86.25% and Kappa coefficient of 8167 obtained to spectral feature fitting classification.
The overall accuracy of 90% and Kappa coefficient of 0.9667 obtained for classification of LSU basic pixel and the maximum overall accuracy level of 97.5% and Kappa coefficient of 0.9667 were obtained for classification of basic spectrum of spectral angle mapper.Table 6 presents the related statistics to the algorithm of determining accuracy on the classified images.Software, At the last step of flowchart of provided processing algorithms, the gathered data were combined with each other by putting layers prepared from different processing methods on each other using GIS10.02Environment and faults and surface geological maps of study region were drawn.Figure 12 represents the faults map of the region which has been overlayed on IRS-PAN image.The faults maps have been prepared through comprising with Lineaments map of the region, with satellite images and removing the roads, waterways and folds.The diagram below the map indicates faults orientation toward the North in Zefreh region.Figure 13 represents the final geological map of the ground's surface in Zefreh region has been prepared by combining results obtained from applying various types of used processing methods in present study on ASTER, LANDSAT8 and IRS-PAN images.

Conclusions
The geological map from the surface of under study area was prepared using remote sensing data with an acceptable accuracy to compare to the map with scale of 1:25,000.In addition the fault maps of the region which prepared automatically are in accordance with Earth fact.The matter indicates the capability of proposed flowchart of present study to surface geological mapping.
At the end of this article, considering input data, we examine the operation of  the algorithms used in this research in the processing of remote sensing images to extract useful information in the preparation of a geological map.The algorithm of principal component analysis (PCA) is a useful method to separate noises and duplicated data from spatial information in separated data packages, which has been also utilized in minimizing noise algorithm.It is also proper to create false color composites with the least amount of duplicated information.According to the fact that false color composites and band ratio are used to emphasize on certain effects and are not directly posed in numerical methods of information extraction, utilizing from the technique significantly helped to visual lithology separation of the region.As it is clear in the images and obtained results, most of the units of images have been separated from each other by this method.In this method to obtain the best results, it would be proper to identify desired effect at first and investigate it in the term of spectrum and then to determine its reflective and absorption bands and based on it, the bands related to reflection and absorption ranges of target effect used in proper sensors with higher spectral resolution in the range of absorption and reflection.In preparing PAN-SHARP images, the Gram-Schmidt Conversion was identified as the most proper method of combining images in present study.The results obtained from this conversion had lower noise and better visual quality compare to other methods of image conversion.
The results of applying maximum likelihood classification on data indicated the matter that LANDSAT8 data are more capable to pixel base classifying compare to ASTER data which is probably due to the high spectral resolution of ASTER data.In lithology with high percentage of material mixing, each pixel of image reflects different spectrums from itself.The matter led that the algorithm didn't classify most of the pixels in related class after introducing training samples.But implementing the technique with LANDSSAT8 image classifies units with mixed material more proper than basic spectrum classification methods.However, it is needed to prepare and exact training samples from all parts of geological units with differences in color tone to achieve optimum results.The validation procedure table of this method was prepared using OLI satellite image which indicates classification accuracy of 94% and kappa coefficient of 0.9430 and hence, explains the proper accuracy of this method in image classification.
To implement spectrum base algorithms, the images of ASTER were used due to having high spectral resolution.The image was also used to prepare spectral library, because utilizing from this image provides better results in output of classification methods compare to the standard libraries.The reason of this matter is that slight changes in the composition of the minerals cause different spectral signature.To select spectral diagram from pure pixels, the PPI index was used which comprising the obtained results with geological map of the region and 1:250,000 scales indicate the medium accuracy of this method in identifying pure pixels.The method is not able to identify all pure pixels of image, especially pixels related to the rock units with mixed materials.But it is more proper in identifying uniform units.
The advantage of SAM method compare to other spectral classification methods is ease and speed of its implementation.However, all of the basic spectrum algorithms have disadvantages in identifying and classifying mixed units due to the nature of these methods.The results obtained from LSU algorithm indicate the distinguishing feature of this method compare to other basic spectrum algorithms in the tem of mixed lithology separation which is relatively impossible in other basic spectrum algorithms.However, obtaining accurate results by this method needs complete understanding of the geology of under study region to be able to separate desired member from other areas correctly in post-processing step.Spectral fitness factor (SFF) is a spectral classification method in which the spectrums inserted to the algorithm get threshold to make it possible for each spectrum to get a special signature through slight differences in similar spectrums.The nature of this algorithm is similar to Sam, but its implementing is more complex compare to other basic spectrum classification methods.However, the method is more sensitive to spectral mixture.It means that, the method provides acceptable responses only in the case of using for pure and uniform units and provides no acceptable response in the case of impurity in considered unit.The results obtained from determining the accuracy of classified images with basic spectrum methods indicate overall accuracy of 86.25% and Kappa coefficient of 0.8167 to spectral fitness basic pixel classification, overall accuracy of 90% and Kappa coefficient of 0.8667 to LSU classification and overall accuracy of 97.5% and Kappa coefficient of 0.9667 to spectral angel mapping method.The matter indicates high accuracy of Sam method in image classification compare to other basic spectrum classification methods.However, it is possible that different results are obtained in different regions and other methods show higher accuracy.The method of automatically extracting lineaments in GEOMATICA Software significantly accelerates lineaments mapping.To input of this algorithm, Images such as Hill-shade and IR-PAN are used in which the edges are visible.The method extracts all linear effects of the image such as waterways, roads and faults, which its results must be modified, based on the desired purposes and other linear effects used in lithology mapping and preparing final lineaments map.
In the case of utilizing from data with higher spatial and spectral resolution compare to the data have been used in present study, more accurate results would be obtained.
How to cite this paper: Mohammadizad, R. and Arfania, R. (2017) Advanced Investigation of Remote Sensing to Geological Mapping of Zefreh Region in Central Iran.

Figure 2 .
Figure 2. The flowchart of processing and data integration for geological mapping.

Figure 3 .
Figure 3. (a) spectral diagram of ASTER data before implementing MNF; (b) spectral diagram of ASTER data after implementing MNF; (c) spectral diagram of LANDSAT8 data before implementing MNF; (d) spectral diagram of LANDSAT8 data after implementing MNF.

Figure 10 .
Figure 10.Rule image classification for calculating accuracy of spectral base classifier methods.

Figure 12 .
Figure 12.The Faults map extracted from Lineaments map.

Figure 13 .
Figure 13.Final geological map of the Earth's surface in Zafreh region prepared by using satellite data of ASTER, LANDSAT8 and IRS-PAN.

Table 1 .
The table related to values of STDV and EIGENVALUE for data of ASTER and LANDSAT8.

Table 2 .
The band ratios applied on bands of ASTER image.
ratio of 5/4 was used to make visible the Ferro silicate areas.The ratio of 12/14 (thermal bands of ASTER sensor) was used to extract areas rich in quartz.In addition, NDVI formula was used to prepare gray scale image which indicated regions covered with vegetation.

Table 3 .
The produced false color composite using data of ASTER and LANDSAT8 to create lithology separation.
from bright green to dark green ranges.In addition, some parts have been distinguished with dark blue and little parts with orange color.The E3BR and E3BA which indicate igneous lithology with ultra-basic and probably serpentinized

Table 4 .
The table of OIF index for ASTER data.

Table 5 .
The table of OIF index for LANDSAT8 data.
combinations are observable with red color.The units of Naiband formation are observed with pink color and composites of yellow color.Intrusive units and water bodies are observable with yellow and ink red colors, respectively Figure7.B7/B6 Gad and Kusky color composite created by band ratios of LANDSAT8 (4/2, 6/5 and 7/6) shows Serpentinized units with color of red-black and met volcanic units are also observable in the image with a color range of green to yellow.In the image, dolomitic areas have been distinguished with green-blue and limestone areas with blue and dark-orange Figure8.B6/B4 Gad and Kusky color composite with band ratios of LANDSAT8 (B6/B4, B6/B2 and B7/B6) was created which indicates Serpentinized rock units with dark color (blue-black)