Next Article in Journal
HoPE: Horizontal Plane Extractor for Cluttered 3D Scenes
Next Article in Special Issue
A Simulation Study Using Terrestrial LiDAR Point Cloud Data to Quantify Spectral Variability of a Broad-Leaved Forest Canopy
Previous Article in Journal
A Secure Multi-Tier Mobile Edge Computing Model for Data Processing Offloading Based on Degree of Trust
Previous Article in Special Issue
Transfer Learning for Soil Spectroscopy Based on Convolutional Neural Networks and Its Application in Soil Clay Content Mapping Using Hyperspectral Imagery
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Band Ratios Matrix Transformation (BRMT): A Sedimentary Lithology Mapping Approach Using ASTER Satellite Sensor

1
School of Earth Sciences, Damghan University, Damghan 3671641167, Iran
2
Korea Polar Research Institute (KOPRI), Songdomirae-ro, Yeonsu-gu, Incheon 21990, Korea
3
Centre for Advanced Modelling and Geospatial Information Systems (CAMGIS), Faculty of Engineering and Information Technology, University of Technology Sydney, Ultimo 2007, Australia
4
Department of Energy and Mineral Resources Engineering, Choongmu-gwan, Sejong University, 209 Neungdong-ro Gwangjin-gu, Seoul 05006, Korea
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(10), 3213; https://doi.org/10.3390/s18103213
Submission received: 18 August 2018 / Revised: 19 September 2018 / Accepted: 22 September 2018 / Published: 23 September 2018

Abstract

:
Remote sensing imagery has become an operative and applicable tool for the preparation of geological maps by reducing the costs and increasing the precision. In this study, ASTER satellite remote sensing data were used to extract lithological information of Deh-Molla sedimentary succession, which is located in the southwest of Shahrood city, Semnan Province, North Iran. A robust and effective approach named Band Ratio Matrix Transformation (BRMT) was developed to characterize and discriminate the boundary of sedimentary rock formations in Deh-Molla region. The analysis was based on the forward and continuous division of the visible-near infrared (VNIR) and the shortwave infrared (SWIR) spectral bands of ASTER with subsequent application of principal component analysis (PCA) for producing new transform datasets. The approach was implemented to ASTER spectral band ratios for mapping dominated mineral assemblages in the study area. Quartz, carbonate, and Al, Fe, Mg –OH bearing-altered minerals such as kaolinite, alunite, chlorite and mica were appropriately mapped using the BRMT approach. The results match well with geology map of the study area, fieldwork data and laboratory analysis. Accuracy assessment of the mapping result represents a reasonable kappa coefficient (0.70%) and appropriate overall accuracy (74.64%), which verified the robustness of the BRMT approach. This approach has great potential and capability for mapping sedimentary succession with diverse local–geological–physical characteristics around the world.

Graphical Abstract

1. Introduction

Preparation of geological maps and mineral exploration using conventional fieldwork investigations are time-consuming and expensive for geologists and mining companies [1,2,3]. Remote sensing plays a key role in geological mapping and mineral exploration, especially in complex, remote and inaccessible areas [4,5,6,7,8,9]. For the preparation of lithological maps using multi-spectral satellite imagery, numerous image processing methods have been developed and used by researchers [10,11]. Spectral indices and band ratios are the most common image processing methods for spectral enhancement to map lithological units [12,13,14,15,16,17].
Sedimentary rocks cover 75% of the Earth’s surface, but a few image processing methods have been proposed for mapping sedimentary lithological units [18,19,20,21]. The lack of specialized image processing method for mapping sedimentary rock units motivated us to develop a robust and easily applicable method to map the sedimentary strata. Some of the significant studies that used ASTER spectral bands for sedimentary rock mapping are summarized as follows. Ninomiya and Fu [18] proposed quartz (QI), carbonate (CI), and silica (SI) indices for the identification and separation of sedimentary rock units such as quartzite, dolomite and limestone, which mainly use the TIR bands of the ASTER. Öztan and Süzen [20] presented a method in order to identify the Ankara gypsum units. This method is based on a combination of decorrelation stretch, band ratio, feature-oriented principal component Analysis (FOPCA), and thermal indices such as CI, QI, and SI. Several Red-Green-Blue (RGB) composites of the band ratios 4/6, 4/9 and 8/6 along with 4/(6 + 8) and 9/(4 + 8) were used to identify gypsum units. Finally, the Sulfate Index (SI) = (10 × 12)/(11 × 11) was developed by the authors for mapping Ankara gypsum units.
Zadeh and Tangestani [21] mapped the sedimentary units (sandstone, tuff, conglomerate, siltstone and limestone) of the southeast of Iran in Kerman province using thermal bands of ASTER. Spectral Angel Mapper (SAM) method was used to prepare the geological map of the area. The Black Hatat area of eastern Oman was studied by Rajendran et al. [22] using Band ratio of 7 + 9/8 for limestone, 6 + 8/7 for dolomite, 2/1 for Fe3+ rich mafic rocks and 6/8 for rocks rich in Quartz (such as shale, schist, and sandstone). Pour et al. [8] developed several ASTER band ratios such as Fe-Minerals Index (Fe-MI) = (4/3) × (2/1), Al-OH-bearing alteration Minerals Index (Al–OH-MI) = (5 × 7)/(6 × 6) and Fe, Mg–OH-bearing alteration Minerals Index (Fe, Mg–OH-MI) = (7 × 9)/(8 × 8) for mapping sedimentary sequences in the Northern Victoria Land, Antarctica.
Accordingly, band ratios and PCA have been widely used for identifying rock and mineral units. Band ratios mostly rely on the absorption and reflection wavelengths of minerals, which have led to the establishment of several spectral indices for specific minerals or rocks. It is a threshold dependent technique. The principal component analysis (PCA) is a well-known method for lithological and alteration mapping in metallogenic provinces [14,15,23], which is a scene dependent method. Therefore, these widely used techniques for lithological–mineralogical mapping may have different outputs for a scene with different size and threshold. It results in low accuracy discrimination between lithological units, especially for sedimentary rock boundaries. Moreover, the application of band ratios for lithological mapping may contain some false results due to the spectral interference of the minerals with similar absorption and reflection properties. Hence, if only band ratio is used for producing a geological map of a study region, some of the spectral interference can result in lithological boundary discrimination with low accuracy.
To overcome these problems, Band Ratios Matrix Transformation (BRMT) was developed in this study. This algorithm depends on the type of rocks in the selected region of interest. It is because of consideration of all possible band ratios for extraction of the main components in the scene. A forward and non-repetitive Band Ratios (BR) matrix of VNIR and SWIR ASTER spectral bands was constructed and, subsequently, PCA transformation was implemented to the matrix. In this study, BRMT approach was established for sedimentary lithological mapping. The Deh-Molla sedimentary succession was selected for testing the BRMT method, which is located in the southwest of Shahrood city, Semnan Province, North Iran (Figure 1). This region has arid to semi-arid climate. Yearly precipitation is very low; hence, the surface is well-exposed due to sparse and non-existent vegetation cover. The Deh-Molla area hosts main stratigraphic units of the Alborz Range such as Fasham, Deh-Sufiyan and Deh-Molla. Moreover, stratigraphic rock units are well-exposed, and the region is easy to access through Damghan-Shahrood highway. Therefore, this area could be considered as a suitable case study area for remote sensing analysis especially for testing the BRMT technique that developed in this imagery. The objectives of this study were: (i) to develop the mathematics performance of BRMT algorithm and its effects on spectral characteristics of the input bands; (ii) to test and apply BRMT algorithm to the Deh-Molla sedimentary stratigraphic units from the imagery; and (iii) to compare the results of BRMT algorithm with well-established ASTER mineralogical/lithological indices and fieldwork and laboratory analysis.

2. Geological Setting and Stratigraphy of the Study Area

The Alborz Mountains of northern Iran extend about 2000 km along the southern margin of the Caspian Sea from Azerbaijan Province (west Iran) to western Afghanistan. The stratigraphic succession of the Alborz has been outlined to the central and eastern Alborz stratigraphic rock units. Their thickness is approximately 11–13 km and shows Pre-Cambrian to Phanerozoic age [24,25,26]. The study area is located near Deh-Molla Village, approximately 15 km southwest of the Shahrood city, Semnan Province, north Iran (Figure 1). The Deh-Molla stratigraphic section is part of the Central Alborz, which is located close to the boundary of eastern Alborz Range. According to geologic map of the Shahrood (scale 1:100000), sedimentary succession of this area started by the Bayandor Formation and reached to the Jurassic strata of the Shemshak Formation (Figure 2) [27].
The Bayandor Formation is the oldest exposed rock units of the Deh-Molla, consisting of 120 m of sandstones with intercalations of micaceous shale and dolomite. The Soltanieh Formation mainly consists of brown to gray dolomites with alternating shale members (latest Proterozoic to Early Cambrian) [26]. The Barut Formation composed of an alternation of dolomites and shale (Early Cambrian). The overlying Zagun Formation in the Deh-Molla section is not well exposed and consists of about 100 m of shales and sandstones. The Lalun Formation mainly consists of purple to red sandstones (Late Early Cambrian) (Figure 2). In the Deh-Molla section, it has about 400 m thickness.
The Fasham Formation is a newly introduced rock unit by Geyer et al. [25]. Type section of this rock unit is also located in the Deh-Molla section. Fasham Formation consists of 60–100 m of white quartz-arenitic sandstones [28]. This stratigraphic shift is a consequence of an unconformity recognized at the base of this quartzarenite unit. It is traditionally called Mila Formation sensu [26], which was subdivided into five members. Mila Group comprises Fasham, Deh-Sufyian, Deh-Molla, and Lashkarak Formations. The Deh-Sufyian Formation could be subdivided into two units that are in accordance with the informal “Member 1” and “Member 2” of the traditional Mila Formation. This formation predominantly composed of carbonate rocks. The basal unit mainly composed of dolostones, which is followed by limestones of the upper unit (Middle Cambrian) [25]. Lithologically the Deh-Molla Formation consists of two recognizable units that correspond to “Member 3” and “Member 4” of the Mila Formation. The lower unit dominated by siliciclastic sandstones with intercalations of carbonates interlayers. The upper unit consists of limestone, marlstone and shale [25].
Ordovician strata of the Alborz Mountain include Member 5 of the Mila, Lashkarak and Gheli Formations and unnamed Formation [25,26,29]. In this study, we preferred to use the term of the Lashkarak Formation for Lower Ordovician deposits of the Deh-Molla area. Lashkarak Formation is mainly composed of clastic deposits such as sandstone and shale. Some narrow igneous and limestone intercalations are also intermittently observed [30]. The Jeirud Formation overlay the Lashkarak deposits with a significant unconformable transition (Figure 2). The Jeirud at the base mainly composed of clastic deposits and towards to up-section contains fossiliferous limestones (Upper Devonian). Mesozoic rock units of the Deh-Molla area include Elika, Shemshak, Dalichay, and Lar Formations (see Figure 2). Carbonates of the Elika Formation (Early to Middle Triassic) are surrounded by Shemshak Formation clastic deposits. The Late Triassic Early Jurassic sandstones and coal-bearing shales of the Shemshak Formation with a thickness of more than 1100 m are one of the main stratigraphic rock units of the Deh-Molla area. These rock units with some coal-bearing strata are considerable for economic exploitation and have some coal mines. Quaternary alluvium deposits covered the river beds and foot-hill parts in the Deh-Molla area. Lithological compositions of the sedimentary formations reported in the Deh-Molla region are summarized in Table 1.

3. Materials and Methods

3.1. ASTER Data and Pre-Processing Analysis

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) sensor consists of three separate instrument subsystems with a total of 14 bands valuable for lithological–mineralogical mapping, including: (i) three bands in the visible-near infrared (VNIR) subsystem from 0.52 to 0.86 μm with a 15 m spatial resolution, which are suitable for iron oxide/hydroxide identification; (ii) six bands in the shortwave infrared (SWIR) subsystem from 1.60 to 2.43 μm with a 30 m spatial resolution, which are useful for clay mineral discrimination; and (iii) five bands in the thermal infrared (TIR) subsystem from 8.12 to 11.65 μm with a 90 m spatial resolution, which can be used for mapping silicate and carbonate minerals [31,32]. ASTER swath width is 60 km that each individual scene is cut to a 60 × 60 km2 area [33].
The ASTER Level 1 T (Precision Terrain Corrected Registered At-Sensor Radiance (AST_L1T)) data acquired on 5 August 2007 was obtained from NASA Earth data center (https://earthdata.nasa.gov/about/daacs/daac-lpdaac) for this study. AST_L1T contains calibrated at-sensor radiance, which corresponds with the ASTER Level 1B (AST_L1B), which has been geometrically corrected and rotated to a north-up UTM projection with WGS84 datum. The reflectance of VNIR and SWIR bands of the datasets has been corrected by incorporating the atmospheric condition through SCP plugin (QGIS 2.18) software [34]. Moreover, Crosstalk correction [35] was also applied to ASTER data. The VNIR reflectance data was resized to 30-m resolution of SWIR bands for further spatial and spectral image processing. It can help to maintain the maximum spectral properties of earth surface materials in the SWIR region, which conveys key spectral information for a wide variety of rock-forming minerals [36].

3.2. Fieldwork, Laboratory Data and Accuracy Assessment

For verification of the image processing analysis, the results were compared with the published geology map of the Shahrood (scale 1:100000) [27] prepared by geological survey of Iran, and fieldwork observation that was carried out on 15 August 2017 in the study area. Geological locations were measured by a GPS (Garmin eTrix 10, Nanjing Sifang Mapping Equipment Ltd., Jiangsu, China) survey. Ground photos were taken of the geomorphology and sedimentary rock units. Samples for laboratory studies were collected through a systematic sampling of fresh and surface-weathered sides of the representative sedimentary lithologies. The XRD analyses were implemented on bulk powder using an X-ray diffractometer, Advance-D8 XRD Bruker model, at Central Laboratory of Damghan University, Damghan, Iran. Moreover, several thin sections of the sedimentary rocks were provided for optical microscopy analyses of sedimentary sequences exposed in the study area.
The accuracy of image classification is most often reported as error matrix include of user’s accuracy (UA), producer’s accuracy (PA), and kappa coefficient (κ) [37]. UA is computed using the number of correctly classified pixels to the total number of pixels assigned to a particular class. PA defines the number of pixels correctly classified in a particular class as a percentage of the total number of pixels actually belonging to that class in the image. The Kappa Coefficient is generated from a statistical test to evaluate the accuracy of a classification [38]. Kappa essentially evaluates how well the classification performed as compared to just randomly assigning values. The Kappa Coefficient can range from −1 to 1. A value of 0 indicated that the classification is no better than a random classification. A negative number indicates the classification is significantly worse than random. A value close to 1 indicates that the classification is significantly better than random [39,40]. In this analysis, accuracy assessment was performed by comparing the map produced by remote sensing analysis and the reference geological map of the study area. For this purpose, the maps have been registered geometrically to each other using accurate GPS reading in the study area. Additionally, a similar classification scheme and spatial detail were used for this assessment.

3.3. Band Ratios Matrix Transformation (BRMT)

A forward and non-repetitive Band Ratios (BR) matrix of VNIR and SWIR ASTER spectral bands was constructed and. subsequently, PCA transformation was implemented to the matrix as a vector (see Equation (1)). Vectorization of the matrix is necessary for PCA implementation.
M m T = P C A [ n 1               n 2               n 3               n 4               n 5               n 6             n 7               n 8     n 9             n 10             n 11             n 12             n 13             n 14           n 15       n 16             n 17             n 18             n 19             n 20             n 21   n 22             n 23             n j               n 25             n 26     n 27             n 28             n 29             n 30     n 31             n 32             n 33     n 34             n 35     n 36     ] = P C A [ n 1 n 2 n 3 . . . . n 36 ] = [ P 1 P 2 P 3 . . . . P 36 ]
where M is the band ratios matrix (BRM), T is the transformation of the matrix using PCA method, and m indicates the number of input bands (assumed to be nine in Equation (1)), which is associated to a total number of VNIR and SWIR bands of ASTER datasets. P1–P36 are resultant principal components, and nj is the band ratio used for the construction of BRM in which the identity number (j) belongs to the set of J = {1, 2, 3, …, j, …, 36}.
Band ratio matrix (BRM) is a set of band ratios constructed by the forward and backward division of VNIR and SWIR bands, mutually (one by one). In the forward condition, the band identity number of numerators is smaller than denominators, i.e., the set of ratios b1/b2, b1/b3, b1/b4, etc.; b2/b3, b2/b4, b2/b5, etc.; b3/b4, b3/b5, b3/b6, etc.; and so on build the BRM. In the backward condition, the band identity number of numerators is larger than the denominator, i.e., the ratios b9/b8, b9/b7, b9/b6, etc.; b8/b7, b8/b6, b8/b5, etc.; b6/b5, b6/b4, b6/b3, etc.; and so on build the BRM. The following matrices in Equations (2) and (3) represent the forward and backward BRM, respectively. The matrix elements are the band ratios of ASTER SWIR and VNIR reflectance bands.
F o r w a r d = ( b 1 b 2         b 1 b 3         b 1 b 4         b 1 b 5         b 1 b 6         b 1 b 7       b 1 b 8         b 1 b 9 b 2 b 3         b 2 b 4         b 2 b 5         b 2 b 6         b 2 b 7         b 2 b 8       b 2 b 9 b 3 b 4         b 3 b 5         b 3 b 6         b 3 b 7         b 3 b 8         b 3 b 9 b 4 b 5         b 4 b 6         b 4 b 7         b 4 b 8         b 4 b 9 b 5 b 6         b 5 b 7         b 5 b 8         b 5 b 9 b 6 b 7         b 6 b 8         b 6 b 9 b 7 b 8         b 7 b 9 b 8 b 9 )
B a c k w a r d = ( b 9 b 1         b 9 b 2         b 9 b 3         b 9 b 4         b 9 b 5         b 9 b 6       b 9 b 7         b 9 b 8 b 8 b 1         b 8 b 2         b 8 b 3         b 8 b 4         b 8 b 5         b 8 b 6       b 8 b 7 b 7 b 1         b 7 b 2         b 7 b 3         b 7 b 4         b 7 b 5         b 7 b 6 b 6 b 1         b 6 b 2         b 6 b 3         b 6 b 4         b 6 b 5 b 5 b 1         b 5 b 2         b 5 b 3         b 5 b 4 b 4 b 1         b 4 b 2         b 4 b 3 b 3 b 1         b 3 b 2 b 2 b 1 )
BRM could be reconstructed for vide verity of multi-spectral remotely sensed datasets by manipulating the matrix elements. Equation (4) is proposed to calculate the total number of the matrix element, Sm, as follows:
S m = i = 1 m 1 [ r 1 = ( m 1 ) r 2 = ( m 2 ) r 3 = ( m 3 ) . . . r i = ( m i ) . . . r m 1 = m ( m 1 ) ] = m 2 m 2  
where m identifies the number of VNIR and SWIR bands, ri is the sum (total number) of elements in ith row of the matrix, and m-1 is the number of rows in the BRM. For example, in this study, for ASTER VNIR + SWIR bands, the Sm and ri were calculated using Equation (4) as follows:
  S m = i = 1 7 [ r 1 = ( 9 1 ) = 8 r 2 = ( 9 2 ) = 7 r 3 = ( 9 3 ) = 6 r 4 = ( 9 4 ) = 5 r 5 = ( 9 5 ) = 4 r 6 = ( 9 6 ) = 3 r 7 = ( 9 7 ) = 2 r 8 = ( 9 8 ) = 1 ] = 81 9 2 = 36  
The above matrices contain 36 band ratios and hold spectral properties of rock constituents in the image. This initial information is generally correlated and requires to be translated in new format. PCA is a transformation method for reducing intercorrelated variables into a few variables called principal components (PCs). This transformation is defined in such a way that the first principal component has the largest possible variance and each succeeding component has a lower variance [41,42]. Correspondingly, BRM forward and backward matrices were converted into new uncorrelated variables (components or bands herein as BT bands) using PCA analysis, which are called ft and bt for forward and backward conditions, respectively (Equations (5) and (6)).
  f t = [   P f 1       P f 2       P f 3       .       .       .       P f i       .       .       .       P f n   ]  
  b t = [ P b 1       P b 2       P b 3       .       .       .       P b i       .       .       .       P b n ]  
where pfi and pbi are the ith principal component of the forward and backward matrices, and n is equal to Sm.
Statistical result of PCA transformation (eigenvalues and eigenvector loadings) is effective to specify most variant components of ft and bt. Equation (7) was used to calculate the variance percent of the components as follow:
  % V i = e i i = 1 n e i × 100    
where ei indicates eigenvectors in ith component, i presents identity number of the components, and n is defined above. Most variant eigenvalues of the forward BRMT are reported in Table 1.
Magnitude and sign of the eigenvector loadings define the contribution of band ratios in each component [43,44]. Equation (8) shows the Pearson correlation coefficient between band ratios (of BRM) and the resulting components, i.e., ft and bt vectors.
  [ r ( p k ,   n j ) ] j , k = 1 s m = N l = 1 N p l n l [ l = 1 N p l × l = 1 N n l ] [ N l = 1 N p l 2 ( l = 1 N n ) 2 ] × [ N l = 1 N p l 2 ( l = 1 N n ) 2 ]  
where r is the set of correlations between jth band ratio (nj) and kth resulting principal components (pj), N indicates the total number of pixels in p (principal components) and n (band ratio) images, and l presents the pixel identity number.
Equation (9) is used to define the average correlation of band ratio vector (n1–n36) in each principal component. This equation calculates the negative and positive averages correlations separately.
  [ r ¯ k ] k = 1 s m = j = 1 m | r ( p k , n j ) | > 0.1 e  
where rk presents the absolute correlation coefficient of the kth principal component (pk,) and the associated band ratio vector (nj), “||” is the absolute symbol, e denotes the total number of band ratios candidates showing meaningful absolute correlation i.e., |r(pk, nj)| > 0.1, and ± symbol is used to switch between negative and positive correlations. The average correlation of band ratio candidates for p1p12 is presented in Table 2 for this analysis.
Equation (10) is used to define the average contribution of each band ratio in the construction of principal component vector. This conditional equation also calculates the negative and positive averages.
  [ r ¯ j ] j = 1 s m = ± k = 1 s m ( | r ( p k , n j ) | > 0.1 ) h  
where rj represents the absolute correlation coefficient of the jth band ratio (nj) and the associated principal component vector (pk), h denotes the total number of components showing meaningful absolute correlation i.e., |r(pk,nj)| > 0.1. The negative and positive correlations for n1–n36 band ratios are presented in Table 3.
To calculate the contribution percentage of the band ratios, Equation (11) is proposed, which is a type of normalization strategy.
  % C j = ± | r ( p k , n j ) | k = 1 s m ( | r ( p k , n j ) | × 100  
where %Cj indicates the contribution of jth band ratio in the construction of resulting principal components, i.e., ft or bt. Table 4 shows the contribution percentage of band ratios n1–n36 for building the resultant principal components.
Furthermore, Equation (8) is used to calculate the correlation coefficients among 36 band ratios, which the highly correlated band ratios (r > 0.90) are defined and listed in Table 5. Considering statistical results derived from Table 1, Table 2, Table 3, Table 4 and Table 5 for the selected sub-scene of the study area, it is evident that some of the BRMT bands show the similar manifestation of sedimentary formations due to high negative correlation, low contribution percentage and high correlation of the band ratios. However, a small number of the BRMT bands specifically map the sedimentary formations in the study area. Table 6 shows BRMT bands with specific manifestation (dark or bright pixels) for the sedimentary formations exposed in the Deh-Molla area. Subsequently, BRMT bands (abbreviated as BT bands) contain specific spectral information were selected as input datasets for post classification using rule image classifier tool. BT1, BT2, BT3, BT4, BT5, BT6, BT7, BT11, BT13, BT18 and BT32 were selected as rule images (input BRMT bands) (Table 7). The threshold value of 0.750 with maximum value option was applied for running the rule image classifier tool.

4. Results

4.1. Fieldwork and Laboratory Analysis Results

Fieldwork and laboratory analysis were carried out to identify the main petrographic and mineralogical characteristics of the sedimentary formations in the Deh-Molla area [27]. Sedimentological and stratigraphical fieldwork observations indicated that lithological units in the study area could be classified into three groups: (i) carbonates; (ii) clastics; and (iii) hybrids. Carbonate rocks are mainly composed of carbonate minerals such as calcite and dolomite. Clastic rocks include quartz, feldspars (potassium and plagioclase feldspars) and several accessory minerals such as clay minerals and mica. Hybrid rock units contain a mixture of carbonate and clastic rocks. Figure 3A–D shows some filed photographs of sedimentary formations exposed in the study area.
Laboratory analysis was performed for extracting petrographic information from thin sections provided for representative samples of sedimentary formations. Moreover, detailed mineralogical information was obtained for the samples using XRD analysis. Soltanieh, Bayandor, Deh-Sufiyan (Mila Group) and Elika Formations are carbonate dominated lithological units in the Deh-Molla area. Figure 4A shows the microphotograph of the Soltanieh representative sample. Dolomite is dominated mineral in the thin section. It seems that brown to red color zones between dolomite crystals are iron oxide minerals due to alteration processes. The XRD analysis for the Soltanieh sample shows dolomite as major minerals and quartz, calcite, orthoclase, mica and hematite as accessory minerals (Figure 5A). The XRD result for the Bayandor representative sample displays only dolomite mineral (Figure 5B). However, petrographic information for Deh-Sufiyan Formation shows calcite and clastic deposits such as fossil shells (echinoderms and brachiopods) as well as quartz and feldspars (Figure 4B). Detected minerals by XRD for Deh-Sufiyan are calcite (major mineral), quartz, orthoclase and microcline (Figure 5C). The XRD analysis reveals the domination of dolomite in Elika with a minor amount of hematite, mica and anorthoclase (Figure 5D).
Clastic rock units in the study area include Mila Group (Lalun, Fasham and Lashkarak) and Shemshak Formations (see Figure 3A–D). Microphotograph of the Lalun Formation displays that it composed of arkosic sandstone (Figure 4C). Quartz is the dominant mineral, and feldspar (mainly orthoclase, albite, and microcline) is also abundant (approximately 25%). In addition, small amounts of mica and clay minerals seem present (Figure 4C). Quartz, orthoclase, muscovite, albite, kaolinite and microcline were detected by XRD analysis (Figure 5E) in the Lalun representative sample. The Fasham Formation consist of mature quartzarenite sandstone (see Figure 4F), it has a white color feature that could be considered as the key bed during fieldwork. According to XRD analysis, it contains quartz as a major mineral with minor amounts of orthoclase and feldspar (Figure 5F). The Lashkarak Formation is an alternation of sandstones and shale (see Figure 3A). Quartz, calcite, chlorite, feldspars, kaolinite and muscovite were detected by XRD analysis in the Lashkarak representative sample (Figure 5G). The Shemshak Formation (see Figure 3B) is mainly composed of quartz, feldspars, mica, anorthoclase, kaolinite, microcline and chlorite based on XRD results (Figure 5H).
The Barut, Deh-Molla (Mila Group) and Jeirud Formations are considered as hybrid class in the study area. The Barut Formation consists of alternation of shale and dolomite beds (see Figure 3C). The XRD analysis results show the presence of quartz, chlorite, illite, kaolinite, montmorillonite and muscovite (Figure 5J). Microphotograph of a representative sample of the Deh-Molla Formation displays quartz, feldspar and few amounts of mica with carbonates cement (Figure 4E). According to XRD data, quartz, calcite, dolomite, feldspars, orthoclase and microcline are detected minerals for the Deh-Molla representative sample (Figure 5K). In the Jeirud Formation a mixture of carbonate and clastic deposits was identified (Figure 4D) in the microphotograph. In contrast with the Deh-Molla thin section (Figure 4E), calcite and dolomite are more dominated in Jeirud thin section (Figure 4D). Based on XRD analysis, calcite, dolomite, kaolinite, muscovite and quartz are dominated minerals, respectively (Figure 5L). Calcite and dolomite are major minerals, while kaolinite, muscovite and quartz are accessory minerals in the Jeirud representative sample.

4.2. BRMT Image Processing Results

Figure 6 shows resultant map derived from the post-classification (rule image classifier) of selected input BRMT bands for the Deh-Molla region. Sedimentary formations are discriminated based on specific spectral information extracted by BRMT approach. The resultant map shows six dominated classes (Figure 6). Color classes indicate different assemblages of minerals including quartz, calcite, dolomite, feldspar and mica. Red pixel class is generally governed in the study area especially in the southeastern and eastern parts. With reference to the geology map of the study area (see Figure 2), red pixel class (C1) is the Soltanieh and Bayandor Formations and Quaternary alluvium. These units mostly consist of carbonate rocks (dolomitic rocks) and weathered limestone. However, intercalations of sandstones and micaceous shale are also reported in the Soltanieh and Bayandor Formations, which are detected as pixels with different color classes in the red class background.
Green pixel class (C2) in the central and southwestern parts of the scene (Figure 6) is the Lalun Formation, which is mainly composed of sandstones. However, green pixels are also mapped in several parts of the study area especially in association with Quaternary alluvium in the southwestern part of the study area. The yellow pixel class (C3) as the third dominated class is governed the central north and northwestern parts of the study area, where the Shemshak Formation is located. It contains clastic deposits and coal-bearing shales. Magenta pixel class (C4) mostly detected with the Jeirud Formation, which is composed of clastic deposits and limestone. However, the magenta pixel is also distributed in the Mila Formation, which is consisting of carbonates and siliciclastic sandstones with intercalations of carbonates, limestone, marlstone and shale. Moreover, this pixel class could be seen with Quaternary alluvium in the southwestern part of the scene (Figure 6). Blue pixel class (C5) is mainly mapped with the Lashkarak Formation. Clastic deposits such as sandstone and shale were reported in the Lashkarak Formation. However, blue pixels could also be found in the Shemshak, Elika, Barut Formations and many other parts of the study area. Cyan pixel class (C6) is not observed specifically with any sedimentary formation in the study area. This class shows less abundance and generally distributed in the background of other classes. The admixture of the six detected classes is mapped with the Barut, Elika and Zagun formations (Figure 6), which are composed of an alternation of carbonates and shale and clastic deposits.
It is evident that similar spectral characteristics of sedimentary rocks influenced the results of BRMT in the classification map. Therefore, for finding the robustness of the BRMT, spectral signatures of the sedimentary formations in the study area compared with spectral properties extracted from BRMT (Figure 7A,B). High similarity in the spectral characteristics of sedimentary rocks is obvious in Figure 7A. Only the Soltanieh and Lalun Formations contain slightly different spectral signatures. It is probably due to their compositions with dominated dolomitic rocks (Soltanieh) and sandstones (Lalun). Spectral signatures of the sedimentary formations after applying BRMT display spectral differentiation between them especially from BT1 to BT3 (very strong).
Considering the statistical results, eigenvalues and %Vi of the forward BRMT for PC1–PC3 show the highest amount for them (see Table 1). As a result, they contain most spectral information in the image. These BRMT PCs (or BT bands) are assigned to RGB (Red, Green, and Blue) color composites for mapping sedimentary formations in the study area. Figure 8 shows the resultant image for Deh-Molla region. Comparing with the geology map of the study area, most of the sedimentary formations appear in recognizable tones. The Soltanieh, Bayandor, and Barut Formations and Quaternary alluvium are represented as yellow, greenish yellow, light orange and green colors, respectively (Figure 8). The Lalun and Zagun Formations manifested as cyan pixels with few blue and green pixels, which distributed in the boundary of these sedimentary formations. The Lashkarak, Jarud, Elika and Shemshak Formations appear as blue to purple color with some disseminations of cyan and blue pixels. The Mila Formations manifest mainly as mustard color. However, several mixing colors and variety of hues are observable in the background of sedimentary formations in the study area (Figure 8). With reference to the BRMT classification map of the Deh-Molla region, it is discernable that the BRMT PCs could not comprehensively identify the sedimentary formations compare to BRMT classification map (see Figure 6). Thus, it is required to consider other available statistical results for producing detailed color composites map of the sedimentary formations in the Deh-Molla region.
The deciphering of statistical results shows that Table 4 may contain some characteristics for finding the bands for high spectral differences. Table 4 shows the contribution percentage of 36 band ratios used for the construction of BRM. The BRs with high negative and positive contribution (C% > 3) are n1–n3, n9, n10, n21, n24–n28, n30, n31, and n33–n36. They contain visual information for considering that whether BRMT- based classification works properly or not. For producing RGB color composites, three BRs having strong contribution percentage were selected. Therefore, n1 (Positive C% = 5.408889 and Negative C% = 2.169569), n9 (Positive C% = 5.075753 and Negative C% = 1.716076), and n34 (Positive C% = 4.077037 and Negative C% = 3.132665) were used to assigning RGB color composites. Figure 9 shows a resultant map for the study area. Results indicate that most of the sedimentary formations were mapped more clearly compared to Figure 5; however, they match well with the BRMT classification map (see Figure 6). The resultant map was able to separate most of the sedimentary formations in the study area properly (Figure 9). The Lalun, Lashkarak, Jeirud, Elika, and Shemshak Formations are mapped and detected in differentiable tones. However, the Soltanieh, Bayandor, Barut, Zagun and Mila Formations and Quaternary alluvium appear in green to brown hues.
The BRMT classification map of the Deh-Molla region (Figure 6) does not include topographical information. However, Figure 7, Figure 8 and Figure 9 contain valuable topographical and texture information. The most obvious texture in the study area that could be seen in the image maps is two sectors separated by a NE–SW line. The SE sector seems to expose a large sedimentary region of recent alluvial fans and minor topography, giving a fine-grained texture. The NW sector is sedimentary hard rocks in a stratigraphical series. Some of the blue pixels in Figure 6 present slope shade effects. The SE sector, gathering young alluvial sediments coming from Lar formation, is composed of limestone. Effects of color compositing have greatly helped to enhance spectral differences in the imagery.
Considering the results achieved from the BRMT approach, it is evident that mapping sedimentary formation is a highly challenging task using ASTER spectral bands. However, this approach was able to extract appropriate information for mapping and discriminating sedimentary formations in the study region. Herein, the results of BRMT approach were compared with some well-established ASTER mineralogical indices such as OH-bearing-altered minerals Index (OHI) = (b7/b6 × b4/b6), Kaolinite Index (KLI) = (b4/b5 × b8/b6), Alunite Index (ALI) = (b7/b5 × b7/b8), Calcite Index (CI) = (b6/b8 × b9/b8), Dolomite Index (DI) = (b6 + b8)/b7) and Quartz Index (QI) = (b11/b10 × b11/b12) [12,14,18]. RGB color composite image map was produced by assigning Kaolinite Index (KLI), Calcite Index (CI) and Quartz Index (QI) for the study area (Figure 10). Moreover, OH-bearing-altered minerals Index (OHI), Alunite Index (ALI) and Dolomite Index (DI) were used to generate an RGB color composite image map (Figure 11). The resultant image maps indicate that the discrimination between the boundaries of sedimentary formations is very week using mineralogical indices. The pixels contain similar tones and colors almost are distributed in the background of all sedimentary formations in the study area (Figure 7 and Figure 8). The mineralogical–lithological indices were not capable of mapping and discriminating the sedimentary formations, including Soltanieh, Bayandor, Barut, Lalun, Mila Group, Jeirud, Shemshak and Elika, in the Deh-Molla region.

5. Discussion

Mapping sedimentary rock is a very challenging task by application of conventional image processing techniques to satellite remote sensing data in a variety of geological environments [8,9,20,45,46]. Lithological and mineralogical indices proposed by Ninomiya and Fu [18] and Ninomiya et al. [12] for mapping sedimentary strata are threshold dependent and could not perform successfully in different geological regions. Furthermore, they may contain some wrong outputs as a result of the spectral interference of the minerals with similar absorption and reflection characteristics. Consequently, if only these indices are used for investigation, the lithological boundary would be mapped with low accuracy. The lack of specialized image processing method for mapping sedimentary rock units encouraged us to develop BRMT to map the sedimentary strata.
In this study, BRMT was developed and established for mapping sedimentary lithology using VNIR and SWIR ASTER satellite sensor spectral bands. In this approach, all forward and non-repetitive possible band ratios of VNIR and SWIR spectral bands of ASTER were considered for extracting the image spectra. Subsequently, the PCA was applied to the selected band ratios for enhancing spectral dissimilarities of the sedimentary rock units. Several mathematical equations have been developed in this analysis. The deciphering of statistical results indicates high potential of the BRMT for increasing the spectral differences among sedimentary lithologies. The BRMT was capable of mapping sedimentary succession in Deh-Molla region, North Iran. Most of the lithostratigraphic units, including Soltanieh, Bayandor, Barut, Lalun, Mila group, Shemshak and Jeirud were mapped and discriminated in the study area. Fieldwork data, petrographic study and XRD analysis revealed three groups of sedimentary rocks such as carbonates, clastics and hybrids in the Deh-Molla region.
Petrographic study and XRD analysis emphasized the complexity of the sedimentary rock constituents and mineral assemblages in the study area. Dolomite, calcite and quartz with minor amounts of feldspars (mainly orthoclase, albite, and microcline) were the main constituents of the most of the sedimentary formation, which identified during the petrographic study. However, some alteration minerals (hematite, kaolinite and chlorite) were also in association with the main constituents of mineralogical content in representative samples of the sedimentary formations. Quartz, dolomite, calcite, orthoclase, microcline, anorthoclase, muscovite, albite, kaolinite, chlorite, montmorillonite and hematite were detected using XRD analysis. Therefore, the sedimentary formations exposed in the study area contain similar spectral characteristics with some small differences due to the variety of their mineralogical content. It indicates that the BRMT algorithm was able to reveal these small differences for geological mapping objectives.
Comparison of the results with the geological map of the study area, fieldwork and laboratory data indicate the robustness of the BRMT technique, despite similar constituents (mineralogical content) and spectral characteristics of sedimentary formations in the study area. Relative accuracy assessment was estimated for the BRMT classification map (Table 8). The overall accuracy of mapping results is 74.64% and the kappa coefficient is 0.70. It shows the robustness of the BRMT method, which was developed and established in this research. Figure 12 shows geological map that produced from the results of BRMT for the Deh-Molla region. It should be noted that Fasham Formation, introduced in this research, was not in the geological map of the study area, hence, corresponding accuracy assessment is not reported for Fasham Formation here.

6. Conclusions

The BRMT approach was developed in this study to characterize and discriminate the boundary of sedimentary rock formations in Deh-Molla region using VNIR and SWIR bands of ASTER remote sensing sensor. Sedimentary Formations exposed in the study region were discriminated and mapped based on specific spectral information extracted by BRMT approach. The results of fieldwork data and XRD analysis and petrographic information indicate three main sedimentary rock groups (carbonates, clastics and hybrids), which were successfully mapped by BRMT approach. Moreover, accuracy assessment of the mapping result indicates an appropriate kappa coefficient of 0.70 and appropriate overall accuracy of 74.64%. The results of this research demonstrated that the BRMT approach has great potential and capability for mapping sedimentary succession with diverse local–geological–physical characteristics in other regions around the world.

Author Contributions

G.A., writing—original draft preparation and developing the methodology and analysis; A.B.P., writing, re-constructing and editing the manuscript; B.P., editing, supervision, and funding acquisition; and M.S. and F.N., field data preparation, laboratory analysis and validation.

Funding

This research is funded by the UTS under grant numbers 321740.2232335 and 321740.2232357.

Acknowledgments

We acknowledge School of Earth Sciences, Damghan University for providing all facilities during this research. Thank to Korea Polar Research Institute (KOPRI) for great assistance during manuscript preparation. The support from the University of Technology Sydney (UTS) is also greatly acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bunte, M.K.; Williams, D.A.; Greeley, R. Geologic mapping of the Zal region of Io. Icarus 2008, 197, 354–367. [Google Scholar] [CrossRef]
  2. Carosi, R.; Montomoli, C.; Iaccarino, S. 20 years of geological mapping of the metamorphic core across Central and Eastern Himalayas. Earth-Sci. Rev. 2018, 177, 124–138. [Google Scholar] [CrossRef]
  3. White, O.L.; Moore, J.M.; McKinnon, W.B.; Spencer, J.R.; Howard, A.D.; Schenk, P.M.; Beyer, R.A.; Nimmo, F.; Singer, K.N.; Umurhan, O.M.; et al. Geological mapping of Sputnik Planitia on Pluto. Icarus 2017, 287, 261–286. [Google Scholar] [CrossRef] [Green Version]
  4. Ahmad, L.; Shah, M.T.; Khan, S.D. Reflectance spectroscopy and remote sensing data for finding sulfide-bearing alteration zones and mapping geology in Gilgit-Baltistan, Pakistan. Earth Sci. Inform. 2016, 9, 113–121. [Google Scholar] [CrossRef]
  5. Pour, A.B.; Hashim, M.; Park, Y.; Hong, J.K. Mapping alteration mineral zones and lithological units in Antarctic regions using spectral bands of ASTER remote sensing data. Geocarto Int. 2017. [Google Scholar] [CrossRef]
  6. Pour, A.B.; Hashim, M.; Hong, J.K.; Park, Y. Lithological and alteration mineral mapping in poorly exposed lithologies using Landsat-8 and ASTER satellite data: North-eastern Graham Land, Antarctic Peninsula. Ore Geol. Rev. 2017. [Google Scholar] [CrossRef]
  7. Pour, A.B.; Park, Y.; Park, T.S.; Hong, J.K.; Hashim, M.; Woo, J.; Ayoobi, I. Evaluation of ICA and CEM algorithms with Landsat-8/ASTER data for geological mapping in inaccessible regions. Geocarto Int. 2018. [Google Scholar] [CrossRef]
  8. Pour, A.B.; Park, Y.; Park, T.S.; Hong, J.K.; Hashim, M.; Woo, J.; Ayoobi, I. Regional geology mapping using satellite-based remote sensing approach in Northern Victoria Land, Antarctica. Polar Sci. 2018, 16, 23–46. [Google Scholar] [CrossRef]
  9. Pour, A.B.; Park, T.S.; Park, Y.; Hong, J.K.; Zoheir, B.; Pradhan, B.; Ayoobi, I.; Hashim, M. Application of Multi-Sensor Satellite Data for Exploration of Zn-Pb Sulfide Mineralization in the Franklinian Basin, North Greenland. Remote Sens. 2018, 10, 1186. [Google Scholar] [CrossRef]
  10. Joshi, N.; Baumann, M.; Ehammer, A.; Fensholt, R.; Grogan, K.; Hostert, P.; Jepsen, M.R.; Kuemmerle, T.; Meyfroidt, P.; Mitchard, E.T.A.; et al. A Review of the Application of Optical and Radar Remote Sensing Data Fusion to Land Use Mapping and Monitoring. Remote Sens. 2016, 8, 70. [Google Scholar] [CrossRef]
  11. Sabatakakis, N.; Nikolakopoulos, K.G.; Papatheodorou, G.; Kelasidis, G. A multisource approach for coastal mapping purposes: Limeni bay, Mani and surrounding area, southern Greece. Earth Sci. Inform. 2016, 9, 183–196. [Google Scholar] [CrossRef]
  12. Ninomiya, Y.; Fu, B.; Cudahy, T.J. Detecting lithology with Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) multispectral thermal infrared radiance-at-sensor data. Remote Sens. Environ. 2005, 99, 127–139. [Google Scholar] [CrossRef]
  13. Bhadra, B.K.; Pathak, S.; Karunakar, G.; Sharma, J.R. ASTER Data Analysis for Mineral Potential Mapping Around Sawar-Malpura Area, Central Rajasthan. J. Indian Soc. Remote Sens. 2013, 41, 391–404. [Google Scholar] [CrossRef]
  14. Pour, B.A.; Hashim, M. The application of ASTER remote sensing data to porphyry copper and epithermal gold deposits. Ore Geol. Rev. 2012, 44, 1–9. [Google Scholar] [CrossRef] [Green Version]
  15. Pour, A.B.; Hashim, M. Identifying areas of high economic-potential copper mineralization using ASTER data in the Urumieh–Dokhtar Volcanic Belt, Iran. Adv. Space Res. 2012, 49, 753–769. [Google Scholar] [CrossRef]
  16. Pour, A.B.; Hashim, M. Integrating PALSAR and ASTER data for mineral deposits exploration in tropical environments: A case study from Central Belt Peninsular Malaysia. Int. J. Image Data Fusion 2015, 6, 37–41. [Google Scholar] [CrossRef]
  17. Guha, A.; Kumar, V.K. New ASTER derived thermal indices to delineate mineralogy of different granitoids of an Archaean Craton and analysis of their potentials with reference to Ninomiya’s indices for delineating quartz and mafic minerals of granitoids—An analysis in Dharwar Craton, India. Ore Geol. Rev. 2016, 74, 76–87. [Google Scholar]
  18. Ninomiya, Y.; Fu, B. Spectral indices for lithologic mapping with ASTER thermal infrared data applying to a part of Beishan Mountains, Gansu, China. In Scanning the Present and Resolving the Future, Proceedings of the IEEE 2001 International Geoscience and Remote Sensing Symposium (Cat. No. 01CH37217), Sydney, NSW, Australia, 9–13 July 2001; IEEE: Piscataway, NJ, USA, 2001; pp. 2988–2990. [Google Scholar]
  19. Özkan, M.; Çelik, Ö.F.; Özyavaş, A. Lithological discrimination of accretionary complex (Sivas, northern Turkey) using novel hybrid color composites and field data. J. Afr. Earth Sci. 2018, 138, 75–85. [Google Scholar] [CrossRef]
  20. Öztan, S.N.; Süzen, L.M. Mapping evaporate minerals by ASTER. Int. J. Remote Sens. 2011, 32, 1651–1673. [Google Scholar] [CrossRef]
  21. Zadeh, M.H.; Tangestani, M.H. Comparison of ASTER thermal data sets in lithological mapping at a volcano-sedimentary basin: A case study from southeastern Iran. Int. J. Remote Sens. 2013, 34, 8393–8407. [Google Scholar] [CrossRef]
  22. Rajendran, S.; Hersi, O.S.; Al-harthy, A.; Al-wardi, M.; El-ghali, M.A.; Al-abri, A.H. Capability of advanced spaceborne thermal emission and reflection radiometer (ASTER) on discrimination of carbonates and associated rocks and mineral identification of eastern mountain region (Saih Hatat window) of Sultanate of Oman. Carbonates Evaporites 2011, 26, 351–364. [Google Scholar] [CrossRef]
  23. Pour, A.B.; Hashim, M. Identification of hydrothermal alteration minerals for exploring of porphyry copper deposit using ASTER data, SE Iran. J. Asian Earth Sci. 2011, 42, 1309–1323. [Google Scholar] [CrossRef]
  24. Assereto, R. The Paleozoic formations in central Elburz (Iran) (preliminary note). Rivista Italiana di Paleontologia e Stratigrafia 1963, 69, 503–543. [Google Scholar]
  25. Geyer, G.; Bayet-Goll, A.; Wilmsen, M.; Mahboubi, A.; Moussavi-Harami, R. Lithostratigraphic revision of the middle Cambrian (Series 3) and upper Cambrian (Furongian) in northern and central Iran. Newsl. Stratigr. 2014, 47, 21–59. [Google Scholar] [CrossRef]
  26. Stöcklin, J.; Ruttner, A.; Nabavi, M. New Data on the Lower Paleozoic and Pre-Cambrian of North Iran; Geological Survey of Iran: Tehran, Iran, 1964; Volume 1. [Google Scholar]
  27. Vaziri, S.H. Geological Map of the Shahrud, Geological Survey of Iran, Scale 1, 100,000; Geological Survey of Iran: Tehran, Iran, 2001. [Google Scholar]
  28. Lasemi, Y.; Amin-Rasouli, H. Interregional Unconformity on the Top Quartzite Unit of the Lalun Formation (Lower and Middle Cambrian Boundary): Evidance Indicating that the Top Quartzite is not the Base Quartzite of the Mila Formation. Geosciences 2006, 15, 128–139. [Google Scholar]
  29. Ghavidel-syooki, M. Palynostratigraphy and palaeogeography of the Cambro-Ordovician strata in southwest of Shahrud City (Kuh-e-Kharbash, near Deh-Molla), central Alborz Range, northern Iran. Rev. Palaeobot. Palynol. 2006, 139, 81–95. [Google Scholar] [CrossRef]
  30. Ghobadi-Pour, M.; Kebriaee-Zadeh, M.R.; Popov, L.E. Early Ordovician (Tremadocian) brachiopods from the Eastern Alborz Mountains, Iran Estonian. J. Earth Sci. 2011, 60, 65. [Google Scholar]
  31. Abrams, M.; Hook, S.; Ramachandran, B. ASTER User Handbook; version 2; Jet Propulsion Laboratory, California Institute of Technology: Pasadena, CA, USA, 2004. Available online: http://glcf.umd.edu/library/guide/aster_user_guide_v2.pdf (accessed on 23 September 2018).
  32. Yamaguchi, Y.I.; Fujisada, H.; Kahle, A.B.; Tsu, H.; Kato, M.; Watanabe, H.; Sato, I.; Kudoh, M. ASTER instrument performance, operation status, and application to Earth sciences. In Scanning the Present and Resolving the Future, Proceedings of the IEEE 2001 International Geoscience and Remote Sensing Symposium (Cat. No. 01CH37217), Sydney, Australia, 9–13 July 2001; IEEE: Piscataway, NJ, USA, 2001; pp. 1215–1216. [Google Scholar]
  33. Yamaguchi, Y.I.; Fujisada, H.; Kudoh, M.; Kawakami, T.; Tsu, H.; Kahle, A.B.; Pniel, M. ASTER instrument characterization and operation scenario. Adv. Space Res. 1999, 23, 1415–1424. [Google Scholar] [CrossRef]
  34. Congedo, L. Semi-Automatic Classification Plugin Documentation Release 5.3.2.1. Available online: https://media.readthedocs.org/pdf/semiautomaticclassificationmanual-v5/latest/semiautomaticclassificationmanual-v5.pdf (accessed on 23 September 2018).
  35. Iwasaki, A.; Tonooka, H. Validation of a crosstalk correction algorithm for ASTER/SWIR. IEEE Trans. Geosci. Remote Sens. 2005, 43, 2747–2751. [Google Scholar] [CrossRef]
  36. Hewson, R.D.; Cudahy, T.J.; Mizuhiko, S.; Ueda, K.; Mauger, A.J. Seamless geological map generation using ASTER in the Broken Hill-Curnamona province of Australia. Remote Sens. Environ. 2005, 99, 159–172. [Google Scholar] [CrossRef]
  37. King, T.V.V.; Clark, R.N. Verification of remotely sensed data. In Remote Sensing for Site Characterization; Kuehn, F., King, T., Hoerig, B., Pieters, D., Eds.; Springer: Berlin, Germany, 2000; pp. 59–61. [Google Scholar]
  38. Smeeton, N.C. Early History of the Kappa Statistic. Biometrics 1985, 41, 795. [Google Scholar]
  39. Villa, A.; Benediktsson, J.A.; Chanussot, J.; Jutten, J. Hyperspectral Image Classification with Independent Component Discriminant Analysis. IEEE Trans. Geosci. Remote Sens. 2011, 49, 4865–4876. [Google Scholar] [CrossRef]
  40. Askari, G.; Li, Y.; Moezzi Nasab, R. An adaptive polygonal centroidal voronoi tessellation algorithm for segmentation of noisy sar images. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2014, XL-2/W3, 65–68. [Google Scholar] [CrossRef]
  41. Cheng, Q.; Jing, L.; Panahi, A. Principal component analysis with optimum order sample correlation coefficient for image enhancement. Int. J. Remote Sens. 2006, 27, 3387–3401. [Google Scholar]
  42. Schowengerdt, R.A. Remote Sensing: Models and Methods for Image Processing, 3rd ed.; Academic Press, Elsevier: Burlington, MA, USA, 2007. [Google Scholar]
  43. Loughlin, W.P. Principal components analysis for alteration mapping. Photogramm. Eng. Remote Sens. 1991, 57, 1163–1169. [Google Scholar]
  44. Crosta, A.P.; Souza Filho, C.R.; Azevedo, F.; Brodie, C. Targeting key alteration minerals in epithermal deposits in Patagonia, Argentina, Using ASTER imagery and principal component analysis. Int. J. Remote Sens. 2003, 24, 4233–4240. [Google Scholar] [CrossRef]
  45. Yao, K.; Pradhan, B.; Idrees, M.O. Identification of rocks and their Quartz content in Gua Musang Goldfield using Advanced Spaceborne Thermal Emission and Reflection Radiometer Imagery. J. Sens. 2017, 2017, 6794095. [Google Scholar] [CrossRef]
  46. Testa, F.J.; Villanueva, C.; Cooke, D.R.; Zhang, L. Lithological and hydrothermal alteration mapping of epithermal, porphyry and tourmaline breccia districts in the Argentine Andes using ASTER imagery. Remote Sens. 2018, 10, 203. [Google Scholar] [CrossRef]
Figure 1. Geographical location of the Deh-Molla inside Shahrood geological map (scale 1:50000) and its location in North of Iran (red cube). Black rectangle shows the subset of the Deh-Molla sedimentary succession in the Shahrood geological map.
Figure 1. Geographical location of the Deh-Molla inside Shahrood geological map (scale 1:50000) and its location in North of Iran (red cube). Black rectangle shows the subset of the Deh-Molla sedimentary succession in the Shahrood geological map.
Sensors 18 03213 g001
Figure 2. Geological map of the Deh-Molla illustrates the sedimentary Formations in the area of study.
Figure 2. Geological map of the Deh-Molla illustrates the sedimentary Formations in the area of study.
Sensors 18 03213 g002
Figure 3. Outcrops of sedimentary lithological units in the Deh-Molla area. (A) Panoramic view of the Jeirud, Lashkarak, and Deh-Molla formations (Mila group); (B) carbonates of the Elika Formation showing fault contact with the Shemshak Formation; (C) dolomite and shales of the Barut Formation; and (D) panoramic view of the Lalun, Fasham, Deh-Sufyian and Deh-Molla (Mila group) Formations.
Figure 3. Outcrops of sedimentary lithological units in the Deh-Molla area. (A) Panoramic view of the Jeirud, Lashkarak, and Deh-Molla formations (Mila group); (B) carbonates of the Elika Formation showing fault contact with the Shemshak Formation; (C) dolomite and shales of the Barut Formation; and (D) panoramic view of the Lalun, Fasham, Deh-Sufyian and Deh-Molla (Mila group) Formations.
Sensors 18 03213 g003
Figure 4. Photomicrographs of lithological units (rock thin sections under a petrographic microscope). (A) Soltanieh Formation showing high abundance of dolomite; (B) Deh-Sufyian limestone composed of calcite minerals and fossil shells; (C) arkosic sandstone of Lalun Formation composed of feldspars (F) and quartz (Q); (D) Jeirud Formation includes hybrid limestone with sandy quartz grains; (E) Deh-Molla Formations composed of quartz (Q) and feldspar (F) and few amounts of mica (M) with carbonates cement; and (F) quartz arenite sandstone of Fasham Formation mainly composed of quartz (Q).
Figure 4. Photomicrographs of lithological units (rock thin sections under a petrographic microscope). (A) Soltanieh Formation showing high abundance of dolomite; (B) Deh-Sufyian limestone composed of calcite minerals and fossil shells; (C) arkosic sandstone of Lalun Formation composed of feldspars (F) and quartz (Q); (D) Jeirud Formation includes hybrid limestone with sandy quartz grains; (E) Deh-Molla Formations composed of quartz (Q) and feldspar (F) and few amounts of mica (M) with carbonates cement; and (F) quartz arenite sandstone of Fasham Formation mainly composed of quartz (Q).
Sensors 18 03213 g004
Figure 5. Results of XRD analysis shows minerals of representative samples collected from: (A) the Soltanieh Formation; (B) Bayandor Formation; (C) Deh-Sufiyan Formation (Mila group); (D) Elika Formation; (E) Lalun Formation (Mila group); (F) Fasham Formation (Mila group); (G) Lashkarak Formation (Mila group); (H) Shemshak Formation; (J) Barut Formation; (K) Deh-Molla Formation (Mila Group); and (L) Jeirud Formation.
Figure 5. Results of XRD analysis shows minerals of representative samples collected from: (A) the Soltanieh Formation; (B) Bayandor Formation; (C) Deh-Sufiyan Formation (Mila group); (D) Elika Formation; (E) Lalun Formation (Mila group); (F) Fasham Formation (Mila group); (G) Lashkarak Formation (Mila group); (H) Shemshak Formation; (J) Barut Formation; (K) Deh-Molla Formation (Mila Group); and (L) Jeirud Formation.
Sensors 18 03213 g005aSensors 18 03213 g005b
Figure 6. BRMT classification map of the Deh-Molla region.
Figure 6. BRMT classification map of the Deh-Molla region.
Sensors 18 03213 g006
Figure 7. (A) Spectral signatures of the sedimentary formations in the study area extracted form band ratios (BR) matrix; (B) spectral signatures of the sedimentary formations after applying BRMT.
Figure 7. (A) Spectral signatures of the sedimentary formations in the study area extracted form band ratios (BR) matrix; (B) spectral signatures of the sedimentary formations after applying BRMT.
Sensors 18 03213 g007
Figure 8. RGB Color composite of BRTM PC1, PC2, PC3 for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Figure 8. RGB Color composite of BRTM PC1, PC2, PC3 for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Sensors 18 03213 g008
Figure 9. RGB color composite image map of most effective BRs as n1 = Red, n9 = Green and n32 = Blue for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Figure 9. RGB color composite image map of most effective BRs as n1 = Red, n9 = Green and n32 = Blue for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Sensors 18 03213 g009
Figure 10. RGB color composite image map of KLI, CI, QI for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Figure 10. RGB color composite image map of KLI, CI, QI for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Sensors 18 03213 g010
Figure 11. RGB color composite image map of OHI, ALI, DI for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Figure 11. RGB color composite image map of OHI, ALI, DI for the Deh-Molla region. Qal, Quaternary alluviums; Js, Shemshak Formation; TRe, Elika Formation; Dj, Jeirud Formation; Ol, Lashkarak Formation; Єdm, Deh-Molla Formation; Єds, De-Sufian Formation; Єf, Fasham; Єl, Lalun; Єz, Zagun Formation; Єbt, Barut Formation; Єs, Soltanieh Formation; PЄbr, Bayandor Formation.
Sensors 18 03213 g011
Figure 12. Geological map the Deh-Molla region derived from the BRMT approach.
Figure 12. Geological map the Deh-Molla region derived from the BRMT approach.
Sensors 18 03213 g012
Table 1. Eigenvalues and %Vi of the forward BRMT for P1P12 for the selected subset of the study area.
Table 1. Eigenvalues and %Vi of the forward BRMT for P1P12 for the selected subset of the study area.
Pfi123456789101112
Eigenvalue0.44880.16370.03960.01920.00850.00710.00380.00270.00240.00230.00050.0004
%Vi64.1523.405.652.741.211.010.540.380.340.320.0730.06
Table 2. Positive and negative correlation averages of band ratio candidates for components 1–12 in selected subset of the study area.
Table 2. Positive and negative correlation averages of band ratio candidates for components 1–12 in selected subset of the study area.
pk123456789101112
+ r ¯ k 0.400.200.130.270.240.170.190.200.190.30--
r ¯ k 0.780.590.280.230.150.260.130.120.490.12-0.12
Table 3. Positive and negative correlation of the band ratios n1–n36 in the construction of resultant principal components for a selected subset of the study area.
Table 3. Positive and negative correlation of the band ratios n1–n36 in the construction of resultant principal components for a selected subset of the study area.
Band RatioPositiveNegativeBand RatioPositiveNegative
n10.36−0.40n190.00−0.47
n20.21−0.39n200.00−0.47
n30.24−0.54n210.00−0.64
n40.11−0.47n220.19−0.54
n50.14−0.47n230.21−0.55
n60.08−0.42n240.22−0.97
n70.10−0.41n250.20−0.97
n80.16−0.42n260.20−0.97
n90.33−0.32n270.59−0.12
n100.37−0.82n280.23−0.51
n110.04−0.30n290.17−0.92
n120.14−0.33n300.22−0.51
n130.05−0.49n310.23−0.39
n140.11−0.49n320.18−0.51
n150.05−0.40n330.20−0.50
n160.16−0.53n340.27−0.58
n170.12−0.42n350.31−0.50
n180.11−0.42n360.28−0.35
Table 4. Contribution percentage of band ratios n1–n36 in the construction of resultant principal components for a selected subset of the study area.
Table 4. Contribution percentage of band ratios n1–n36 in the construction of resultant principal components for a selected subset of the study area.
Band RatioPositive C%Negative C%Band RatioPositive C%Negative C%
n15.4088892.169569n1902.526033
n23.128042.091521n2002.521348
n33.6988082.903388n2103.465243
n41.6022832.554337n222.8513622.940772
n52.0625882.557068n233.2151672.987117
n61.1671842.274883n243.3069415.218326
n71.5825572.216687n252.9833185.225227
n82.493822.259561n263.1198335.23804
n95.0757531.716076n278.9979910.647781
n105.563764.453769n283.4388182.756674
n110.6002331.610137n292.5697074.984671
n122.1187821.789384n303.3856522.742368
n130.8359922.624651n313.541562.096501
n141.678632.619609n322.7859462.775086
n150.7806542.181139n333.0203962.717943
n162.5084532.862117n344.0770373.132665
n171.7540942.248454n354.6978842.710994
n181.7018262.247719n364.3201031.912448
Table 5. Band ratios with the correlation coefficient higher than 0.90 for a selected subset of the study area.
Table 5. Band ratios with the correlation coefficient higher than 0.90 for a selected subset of the study area.
Band Ratio NumberCorrelated Band Ratio Band Ratio NumberCorrelated Band Ratio
n6n7, n8n24n22
n11n12n25n22
n13n6, n7, n8n26n22, n23, n24, n25
n14n6,n7,n8n28n24, n25, n26
n15n7, n8, n14n29n24, n25, n26, n28
n17n13n30n24, n25, n26, n28, n29
n21n7, n8, n13, n14, n15n31n24, n25, n26, n28, n29, n30
n19n13, n14, n15, n17, n18n32n24, n25, n26, n28, n29, n30, n31
n20n13, n14, n15, n17, n19, n20n33n24, n25, n26, n28, n30, n31, n33
n23n22
Table 6. BR bands with a specific manifestation (dark or bright pixels) of the sedimentary formations in the study area.
Table 6. BR bands with a specific manifestation (dark or bright pixels) of the sedimentary formations in the study area.
Band RatioBright FormationDark Formation
n1 Lalun
n2Lashkarak, Barut
n3Lalun
n4, n5, n6, n7, n8Alluvials
n9 Deh-Molla, Jeirud
n10Lalun
n13, n14, n15Alluvial
n22Lalun
n23, n24, n25, n26Alluvial
n22, n23, n24, n25, n26Soltanieh, Deh-sufyianLalun
n27 Fasham
n28, n29, n30, n31, n32, n33
n35, n36Fasham
Table 7. Selected input BRMT bands with a specific manifestation (dark or bright pixels) of the sedimentary formations for post-classification.
Table 7. Selected input BRMT bands with a specific manifestation (dark or bright pixels) of the sedimentary formations for post-classification.
BT BandsBright FormationDark Formation
BT1LalunDeh-Sufyian, Deh-Molla, Alluvials
BT2 Lalun
BT3 Deh-Molla, Elika
BT4Soltanieh, Lalun, JeirudShemshak
BT5Deh-Sufyian
BT6Soltanieh
BT7Fasham
BT13LalunDeh-Sufyian, Soltanieh
BT18Barut
BT32Fasham
Table 8. Relative accuracy estimated for lithological mapping results based on BRMT method via Shahrood geology map.
Table 8. Relative accuracy estimated for lithological mapping results based on BRMT method via Shahrood geology map.
FormationsProd. Acc.User Acc.Prod. Acc.User Acc.
(Percent)(Percent)(Pixels)(Pixels)
Elika98.6360.543314/33603314/5474
Shemshak75.170.718304/110578304/11744
Jeirud59.8889.34723/78874723/5289
Lashkarak81.869.884388/53644388/6279
Deh-Molla5772.574092/71794092/5639
Deh-Sufyian3035.68628/2093628/1760
Fasham----
Lalun78.0581.868632/110598632/10545
Barut45.1272.321915/42441915/2648
Soltanieh77.1955.235949/77075949/10771
Bayandor66.4259.711242/18701242/2080
Alluvium85.5183.9226,057/3047226057/31048
Over. Acc.PercentPixel
74.6470811/94860
Kappa Coef.0.70

Share and Cite

MDPI and ACS Style

Askari, G.; Pour, A.B.; Pradhan, B.; Sarfi, M.; Nazemnejad, F. Band Ratios Matrix Transformation (BRMT): A Sedimentary Lithology Mapping Approach Using ASTER Satellite Sensor. Sensors 2018, 18, 3213. https://doi.org/10.3390/s18103213

AMA Style

Askari G, Pour AB, Pradhan B, Sarfi M, Nazemnejad F. Band Ratios Matrix Transformation (BRMT): A Sedimentary Lithology Mapping Approach Using ASTER Satellite Sensor. Sensors. 2018; 18(10):3213. https://doi.org/10.3390/s18103213

Chicago/Turabian Style

Askari, Ghasem, Amin Beiranvand Pour, Biswajeet Pradhan, Mehdi Sarfi, and Fatemeh Nazemnejad. 2018. "Band Ratios Matrix Transformation (BRMT): A Sedimentary Lithology Mapping Approach Using ASTER Satellite Sensor" Sensors 18, no. 10: 3213. https://doi.org/10.3390/s18103213

APA Style

Askari, G., Pour, A. B., Pradhan, B., Sarfi, M., & Nazemnejad, F. (2018). Band Ratios Matrix Transformation (BRMT): A Sedimentary Lithology Mapping Approach Using ASTER Satellite Sensor. Sensors, 18(10), 3213. https://doi.org/10.3390/s18103213

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop