Applying Built-Up and Bare-Soil Indices from Landsat 8 to Cities in Dry Climates

Arid and semi-arid regions have different spectral characteristics from other climatic regions. Therefore, appropriate remotely sensed indicators of land use and land cover types need to be defined for arid and semi-arid lands, as indices developed for other climatic regions may not give plausible results in arid and semi-arid regions. For instance, the normalized difference built-up index (NDBI) and normalized difference bareness index (NDBaI) are unable to distinguish between built-up areas and bare and dry soil that surrounds many cities in dry climates. This paper proposes the application of two newly developed indices, the dry built-up index (DBI) and dry bare-soil index (DBSI) to map built-up and bare areas in a dry climate from Landsat 8. The developed DBI and DBSI were applied to map urban areas and bare soil in the city of Erbil, Iraq. The results show an overall classification accuracy of 93% (κ = 0.86) and 92% (κ = 0.84) for DBI and DBSI, respectively. The results indicate the suitability of the proposed indices to discriminate between urban areas and bare soil in arid and semi-arid climates.


Introduction
In recent decades, different remote sensing methods such as aerial photography and satellite imaging have become widely available as a source of data for mapping and monitoring land use and land cover.Due to rapid urbanization, land cover around many of the world's urban areas changes more quickly than ever before.The data obtained from remote sensing platforms provide up-to-date information and a general view of landscape characteristics and changes in urban areas [1].
Bare land, whether seasonal or permanently bare, plays an important role in many ecosystems, for instance, as a source of Aeolian deposits during sandstorms.However, the differentiation of base land from urban areas in arid and semi-arid environments has often proven difficult.There are different Land 2018, 7, 81 2 of 13 methods and techniques that have been used to identify urban areas from other land cover types, for example, classification into urban, rural and other classes.Night-light satellite data can also be used to distinguish between urban and non-urban areas [2].The Normalized Difference Built-up Index (NDBI) was proposed in 2003 as a new method which automatically maps built-up areas.This method was first applied in practice to extract data for a built-up area of Nanjing city in eastern China [1].Although the NDBI index was able to distinguish between built-up and vegetated or green and wet surroundings for the city of Nanjing, it was not successful in distinguishing between built-up and other land cover such as bare and dry soil that surround the city, due to the overlapping spectral reflectance for these land cover types [3].Although several techniques to observe land use and land cover have been developed, an evaluation of the performance and accuracy of these techniques in dry climate regions is rarely available.We have conducted several tests using NDBI and NDBaI derived from Landsat TM, ETM+ and OLI for LULC of Erbil.The NDBI results produced a high value for bare soil; in contrast, NDBaI produced a high value for built-up.Generally, NDBI is more efficient in places where the NDVI value is greater than 0 [1].Therefore, this method is not suitable for cities located in dry climates (Table 1).Spectral indices are a good approach to distinguish land cover types [4].Chen [5] proposed a bare soil index (BI) primarily for bare land extraction from satellite data.The normalized difference soil index (NDSI) has been used to detect signature changes in un-mixing coastal marsh from satellite images [6].A normalized difference bareness index (NDBaI) was proposed by Chen et al. [4] to distinguish bare land from other land use classes using Landsat data.As the output of this method showed a higher NDBaI value of built-up areas than bare land, this index appears not to be appropriate for use in cities in semi-arid environments.As-syakur et al. [7] proposed an enhanced built-up and bareness index (EBBI) and provided a case study of Denpasar city in Bali, Indonesia, with 90.5% overall accuracy of bare land detection [8].A normalized difference bare land index (NBLI) and unsupervised classification was used by Li et al. [8] to automatically map bare land from Landsat images.Bare soil indices, in particular most commonly NDBaI, are frequently used in humid regions with high accuracy while bare soil in dry areas requires improved indices.
Impervious surface materials and bare soil are both parts of the land surface [9].Sub-pixel and layered classification was proposed by Ji and Jensen [10] for coastal and urban environments.This method was later modified by Flanagan and Civco [11] by using artificial neutral networks to retrieve the impervious surface fraction.Lu and Weng [12] developed a new method for urban land use classification based on the combined use of impervious surface and population density.The Urban Index (UI) [13] performs better in the identification of urban areas but is unable to distinguish bare land and built-up areas accurately [14].Wu and Murry [15] employed a linear mixture pixel analysis (LMPA) to extract impervious distribution.Similarly, Xu [16] proposed a normalized impervious surface index (NDISI) for the estimation of impervious surfaces.It was demonstrated that NDISI can extract impervious surfaces efficiently.Based on NDBI, a continuous built-up index (BU c ) was developed by He et al. [3].Stathakis et al. [17] suggested a vegetation index built-up index (VIBI) to segment urban areas efficiently.Deng and Wu [18] proposed a biophysical composition index (BCI).Land in arid and semi-arid regions has different spectral characteristics from other climatic regions.Therefore, appropriate remotely sensed indicators of land use and land cover types need to be defined In 2013, Landsat 8 was launched with two sensors, the Operational Land Imager (OLI) and the Thermal Infrared Sensor (TIRS).To monitor long-term land cover and land use change, Landsat data are crucial.To minimize atmospheric absorption features, the width of several OLI bands in Landsat 8 was enhanced [19].In 2014, Bhatti and Tripathi proposed a built-up area extraction method (BAEM), and similarly, Zhou et al. used a built-up and bare-land index (BBI) to extract built-up and bare soil from Landsat 8 [20,21].Piyoosh and Ghosh [22] developed a normalized ratio urban index (NRUI) and modified normalized difference soil index (MNDSI) to distinguish between urban areas and soil from Landsat 8.
The reviewed literature on spectral indices leads us to conclude that these indices have not entirely successfully addressed the confusion between impervious surface and bare soil [18,23].Hence, this study focusses on identifying appropriate remotely sensed indices for arid and semi-arid regions derived from Landsat 8 data.It addresses the difference between bare soil and built-up classes that are inherent in other indices for mapping land use land cover classes in cities in dry climates in anticipation that this will provide better information than that obtainable from previously published spectral indices.Erbil city in Iraq was selected for this research and a Landsat 8 image was obtained in the dry season for examining the effectiveness of the proposed techniques.The objective of the current study was to determine the accuracy of the newly proposed indices (dry built-up index (DBI) and dry bare-soil index (DBSI)) in mapping built-up and bare soil areas in dry climate from Landsat 8.

Study Area
Erbil is located 412 m above sea level within 36 • 08 N-36 • 14 N and 43 • 57 E-44 • 03 E [24].The study area is located in the central part of northern Iraq serving as the capital of Kurdistan region (Figure 1).The population of the Erbil Governorate was estimated to be around 1,500,000 inhabitants in 2015 [25].
Land 2018, 7, x FOR PEER REVIEW 3 of 13 are crucial.To minimize atmospheric absorption features, the width of several OLI bands in Landsat 8 was enhanced [19].In 2014, Bhatti and Tripathi proposed a built-up area extraction method (BAEM), and similarly, Zhou et al. used a built-up and bare-land index (BBI) to extract built-up and bare soil from Landsat 8 [20,21].Piyoosh and Ghosh [22] developed a normalized ratio urban index (NRUI) and modified normalized difference soil index (MNDSI) to distinguish between urban areas and soil from Landsat 8.The reviewed literature on spectral indices leads us to conclude that these indices have not entirely successfully addressed the confusion between impervious surface and bare soil [18,23].Hence, this study focusses on identifying appropriate remotely sensed indices for arid and semi-arid regions derived from Landsat 8 data.It addresses the difference between bare soil and built-up classes that are inherent in other indices for mapping land use land cover classes in cities in dry climates in anticipation that this will provide better information than that obtainable from previously published spectral indices.Erbil city in Iraq was selected for this research and a Landsat 8 image was obtained in the dry season for examining the effectiveness of the proposed techniques.The objective of the current study was to determine the accuracy of the newly proposed indices (dry built-up index (DBI) and dry bare-soil index (DBSI)) in mapping built-up and bare soil areas in dry climate from Landsat 8.

Study Area
Erbil is located 412 m above sea level within 36°08′ N-36°14′ N and 43°57′ E-44°03′ E [24].The study area is located in the central part of northern Iraq serving as the capital of Kurdistan region (Figure 1).The population of the Erbil Governorate was estimated to be around 1,500,000 inhabitants in 2015 [25].
The city has a semi-arid continental climate, a rainy cool winter and dry and warm summer.The seasonal distribution of the precipitation in the study area varies and falls mostly as rain in winter and autumn [26].The average annual precipitation is 380.26 ± 108.88 mm.The annual air temperature is 22 °C, and July and August are the hottest months of the year [24].The dominant land use type in the city is residential land use and concrete blocks are the main material for buildings.In general, winter grains are the most common form of land cultivation in the area and depend on rainfall; therefore, in the summer, the majority of croplands are dry [27,28].The city has a semi-arid continental climate, a rainy cool winter and dry and warm summer.The seasonal distribution of the precipitation in the study area varies and falls mostly as rain in winter and autumn [26].The average annual precipitation is 380.26 ± 108.88 mm.The annual air temperature is 22 • C, and July and August are the hottest months of the year [24].The dominant land use type in the city is residential land use and concrete blocks are the main material for buildings.In general, winter Land 2018, 7, 81 4 of 13 grains are the most common form of land cultivation in the area and depend on rainfall; therefore, in the summer, the majority of croplands are dry [27,28].

Data
A Landsat 8 image from 1 July 2013 on path 139 and row 35 was used in this research.The data were provided as level L1T data captured under clear atmospheric conditions (cloud coverage = 0.01%) in the dry season.The image is in UTM projection (zone 38 N) and was retrieved from the United States Geological Survey (USGS), Earth Explorer website [29].

Methodology
To overcome the limitations in the aforementioned indices we used Blue, TIR1, SWIR1 and Green bands to separate built-up and bare areas in the dry climate of Iraq from Landsat 8 and to assess the performance of these indices using accuracy assessment.

Dry Built-Up Index
In general, all the four classes in the study area (bare soil: black colour, built-up: red colour, water: blue colour and green area: green colour) have their highest digital number (DN) in band 10; bare land as the first, built-up class as the second highest DN in this band (Figure 2).In contrast, the lowest DN of all classes is located in band 6 (SWIR-1).In most bands, the DN of barren land is higher than that of built-up areas; only in bands 1 and 2 is the spectral signature of the built-up class higher than that of bare land.The DN values of water is lower than the DN of all other classes across all bands.
Spectral disparity of built-up and bareness classes is low in the Blue band; however, it is large in the TIR1 band (Figure 2).The urban areas had higher blue reflectance than bare soil as a result of the type of building materials, mainly concrete, used for roof surfaces and walls.After testing Landsat 8 bands, we found that these classes can be better distinguished in Blue and TIR1 bands in the study area, thus the proposed equation for built-up area in dry climate is: where ρ Blue = band 2, ρ TIR1 = band 10.
A mask was applied to water surface pixels using water index NDWI.
where NDVI is normalized difference vegetation index, ρ N IR is the surface reflectance of band 5 and ρ RED is the surface reflectance of band 4 in Landsat 8.The reason for subtracting the DBI layer from the NDVI values is that built-up areas can be identified better [3,20,30].Through this step, DBI values of green vegetation will decrease and the features with higher DBI values but lower NDVI values will be enhanced.
The DBI values can be between −2 to +2 and the higher numbers represent more built-up areas.An appropriate threshold of the built-up class from the output of the equation can be used for mapping built-up and non-built-up areas.A test was conducted on the number of built-up pixels by selecting values 0.72 and higher as built-up for the study area while lower than this threshold as non-built up class.Note that the threshold of DBI and DBSI to extract built-up and bare-soil may vary based on the study site.

Dry Bareness Index
The result of testing common bare soil index (NDBaI) indicated that bare soil has lower mean NDBaI than mean NDBaI of built-up areas, meaning that this index was not able to accurately detect bare soil in the study area.DN of bareness and built-up classes is low in band 6 and high in band 3. Inspection of the Landsat 8 bands suggested that differentiation of these classes could be done based on spectral values in the SWIR1 and Green bands.In these bands, generally the DN of bare land is slightly higher than the DN of the built-up class.Thus, the proposed equation for bareness area in dry climate is the inverse of the Modified Normalized Difference Water Index (MNDWI; [31]): where = the surface reflectance of band 6, = the surface reflectance of band 3 in Landsat 8.
The DBSI values can be between −2 to +2, and higher numbers represent more bare soil.An appropriate threshold for the bare soil class can be used for mapping bare soil and non-bare soil areas.Based on a test carried out with a sample of bare soil pixels, a DBSI value 0.26 and higher was delineated as bare soil for the study area, and areas with lower values were delineated as other classes.

Accuracy Assessment
In order to assess the accuracy of classified images based of the DBI and DBSI indices, 300 points were randomly selected for each class using ArcGIS (Figure 3).Each reference point represents one pixel on the classified image.It means overall, 300 pixels were randomly considered.For DBI, 150 points were selected for the built-up class and 150 points for the non-built-up class.A similar procedure was repeated for DBSI, with 150 points for bare soil and 150 points for other classes.Afterwards, these pixels were compared with high-resolution satellite images from the online ArcGIS 10.3 basemaps (Image source: Esri, DigitalGlobe, GeoEye, i-cubed, Earthstar Geographics, CNES/Airbus DS, USDA, USGS, AEX, Getmapping, Aerogird, IGN, IGP, swisstopo, and the GIS User Community).

Dry Bareness Index
The result of testing common bare soil index (NDBaI) indicated that bare soil has lower mean NDBaI than mean NDBaI of built-up areas, meaning that this index was not able to accurately detect bare soil in the study area.DN of bareness and built-up classes is low in band 6 and high in band 3. Inspection of the Landsat 8 bands suggested that differentiation of these classes could be done based on spectral values in the SWIR1 and Green bands.In these bands, generally the DN of bare land is slightly higher than the DN of the built-up class.Thus, the proposed equation for bareness area in dry climate is the inverse of the Modified Normalized Difference Water Index (MNDWI; [31]): where ρ SW IR1 = the surface reflectance of band 6, ρ Green = the surface reflectance of band 3 in Landsat 8.
The DBSI values can be between −2 to +2, and higher numbers represent more bare soil.An appropriate threshold for the bare soil class can be used for mapping bare soil and non-bare soil areas.Based on a test carried out with a sample of bare soil pixels, a DBSI value 0.26 and higher was delineated as bare soil for the study area, and areas with lower values were delineated as other classes.

Accuracy Assessment
In order to assess the accuracy of classified images based of the DBI and DBSI indices, 300 points were randomly selected for each class using ArcGIS (Figure 3).Each reference point represents one pixel on the classified image.It means overall, 300 pixels were randomly considered.For DBI, 150 points were selected for the built-up class and 150 points for the non-built-up class.A similar procedure was repeated for DBSI, with 150 points for bare soil and 150 points for other classes.Afterwards, these pixels were compared with high-resolution satellite images from the online ArcGIS 10.3 basemaps (Image source: Esri, DigitalGlobe, GeoEye, i-cubed, Earthstar Geographics, CNES/Airbus DS, USDA, USGS, AEX, Getmapping, Aerogird, IGN, IGP, swisstopo, and the GIS User Community).

Mapping Built-Up Areas Using the Dry Built-Up Index
The DBI (Figure 4) was calculated from Equation ( 1) for the study area.Figure 4 illustrates that low DBI is present outside the city border.In contrast, high DBI is inside the city border and within the 100 m Street.Average DBI within the city borders is 0.74 while average DBI outside the city borders is 0.65 (Table 1).Figures 5 and 6 illustrate the comparison between DBI and NDBI. Figure 5 showed that the difference between the built up and bare-soil is higher in DBI than for NDBI and UI.DBI built-up land has the highest value based on the purpose of highlighting class by the index while in NDBI, UI and EBBI, bare-soil has the highest value.

Mapping Built-Up Areas Using the Dry Built-Up Index
The DBI (Figure 4) was calculated from Equation ( 1) for the study area.Figure 4 illustrates that low DBI is present outside the city border.In contrast, high DBI is inside the city border and within the 100 m Street.Average DBI within the city borders is 0.74 while average DBI outside the city borders is 0.65 (Table 1).Figures 5 and 6 illustrate the comparison between DBI and NDBI. Figure 5 showed that the difference between the built up and bare-soil is higher in DBI than for NDBI and UI.DBI built-up land has the highest value based on the purpose of highlighting class by the index while in NDBI, UI and EBBI, bare-soil has the highest value.

Mapping Built-Up Areas Using the Dry Built-Up Index
The DBI (Figure 4) was calculated from Equation ( 1) for the study area.Figure 4 illustrates that low DBI is present outside the city border.In contrast, high DBI is inside the city border and within the 100 m Street.Average DBI within the city borders is 0.74 while average DBI outside the city borders is 0.65 (Table 1).Figures 5 and 6 illustrate the comparison between DBI and NDBI. Figure 5 showed that the difference between the built up and bare-soil is higher in DBI than for NDBI and UI.DBI built-up land has the highest value based on the purpose of highlighting class by the index while in NDBI, UI and EBBI, bare-soil has the highest value.Land 2018, 7, 81 10 of 13

Mapping Bareness Areas Using Dry Bare-Soil Index
The DBSI map (Figure 5) was produced from Equation (4) for Erbil city.Bare soil includes both bare areas and seasonally bare farmland.Figure 5 illustrates that high DBSI is found outside the city borders, for instance, in the north and north-west of the study area.In contrast, low DBSI is seen inside the city borders and inside the 100 m street.The average DBSI within the city borders is 0.14 and DBSI outside the city borders is 0.21 (Table 2).This result confirms that the proposed index is able to delineate bare soil with a higher value of DBSI compared to other classes such as built-up areas.Figures 6 and 7 show the comparison between DBSI and NDBaI.

Mapping Bareness Areas Using Dry Bare-Soil Index
The DBSI map (Figure 5) was produced from Equation (4) for Erbil city.Bare soil includes both bare areas and seasonally bare farmland.Figure 5 illustrates that high DBSI is found outside the city borders, for instance, in the north and north-west of the study area.In contrast, low DBSI is seen inside the city borders and inside the 100 m street.The average DBSI within the city borders is 0.14 and DBSI outside the city borders is 0.21 (Table 2).This result confirms that the proposed index is able to delineate bare soil with a higher value of DBSI compared to other classes such as built-up areas.Figures 6 and 7 show the comparison between DBSI and NDBaI.

Acurracy Assessment
The results of the accuracy assessment show that the generated maps have an overall accuracy of 93% (κ = 0.86) and 92% (κ = 0.84) for DBI and DBSI, respectively (Tables 3 and 4).A few pixels of bare rock in high and mountainous area are mixed with built-up areas due to their reflectance resembling the main building materials in this area.The results indicate that the proposed indices can be used reliably for differentiating constructed and bare land from other land use classes in arid and semi-arid climates.

Acurracy Assessment
The results of the accuracy assessment show that the generated maps have an overall accuracy of 93% (κ = 0.86) and 92% (κ = 0.84) for DBI and DBSI, respectively (Tables 3 and 4).A few pixels of bare rock in high and mountainous area are mixed with built-up areas due to their reflectance resembling the main building materials in this area.The results indicate that the proposed indices can be used reliably for differentiating constructed and bare land from other land use classes in arid and semi-arid climates.

Discussion
To distinguish urban areas and bare land using spectral indices is commonly associated with low accuracy due to the high degree of homogeneity [7].Due to low moisture content in green areas surrounding cities, built-up and bare land indices that were developed for humid regions can often not sufficiently differentiate between built-up and bare land in drylands (Tables 1 and 2; Figures 5 and 6).In past studies, NDBI identified built-up areas accurately in cites in humid climates, for instance, in Colombo city in Sri Lanka [32], Montreal in Canada [30], São José dos Campos city in Brazil [23], Beijing and Guangzhou city in China [33,34].In contrast, NDBI has performed poorly for mapping built-up areas in the semi-arid cities of Urumqi and Shihezi in western China [35,36].Zhou et al. [21] also reported a low accuracy (57.4%) of NDBI for Zhengzhou city in China from Landsat 8 OLI data.
Using Blue with TIR1 and SWIR1 with Green bands, the proposed DBI and DBSI techniques can map built-up and bare land in dry climate zones at an accuracy of 93% and 92%, respectively.The results from DBI and DBSI are highly similar to land use classes from high-resolution satellite images.This result confirms that the proposed indices are able to detect built-up land from high DBI values and bare-soil land from high DBSI values.Applying the proposed indices in dry climate zones is shown to be more suitable.Thus, using the proposed indices for different backgrounds of urban areas may produce better and more accurate classification.This could be replicated in many fields of research, for example, in urban expansion and sustainability of cities.
Built-up areas have higher spectral values than bare land in the Blue band (Band 2), and in contrast, in the TIR1 band they have lower spectral values.The built-up index is accurate for distinguishing between land use/land cover classes using the TIR channels that display high emissivity in urban areas [7].The main disadvantage of use thermal band is it has lower spatial resolution; however, because in Erbil and similar semi-arid areas, built-up areas generally exhibit higher values in the Blue and lower in the TIR band, it improves the accuracy of classification.Diminishing NDVI in the DBI and DBSI indices significantly improved accuracy of built-up mapping and addresses the ambiguity in areas with high NDVI (i.e., green areas) that are often confused with built-up areas [3,20,30].Generally, high values of DBI and DBSI represent intense built-up and bare-soil, respectively, but the threshold to extract built-up areas and bare soil may vary between study sites.

Conclusions
Appropriate remotely sensed indicators of land use and land cover need to be defined for arid and semi-arid regions because these areas have different reflectance characters.As a consequence, indices developed for temperate regions may not give adequate results in dry climate areas.Since indices as NDBI and NDBaI are unable to distinguish between built-up areas and bare land that often surrounds cities in dry climates, this study proposes the application of two new spectral indices, DBI and DBSI derived from Landsat 8 data for mapping built-up and bare soil areas in dry climates.
The developed DBI and DBSI were applied to map urban areas and bare land in the city of Erbil, Iraq.The accuracy assessment shows an overall accuracy of 93% (κ = 0.86) and 92% (κ = 0.84) for DBI and DBSI, respectively.The results indicate that the proposed indices can be reliably used to differentiate built-up and bare land from other land use classes in arid and semi-arid climate.
semi-arid lands, as indices developed for other climatic regions may not give plausible results in those regions.

Figure 1 .
Figure 1.Map of Iraq, the Kurdistan Region and true colour Landsat image of the study area.Figure 1. Map of Iraq, the Kurdistan Region and true colour Landsat image of the study area.

Figure 1 .
Figure 1.Map of Iraq, the Kurdistan Region and true colour Landsat image of the study area.Figure 1. Map of Iraq, the Kurdistan Region and true colour Landsat image of the study area.

Figure 2 .
Figure 2. Mean spectral profiles for four typical land covers in the study area.

Figure 2 .
Figure 2. Mean spectral profiles for four typical land covers in the study area.

Figure 5 . 13 Figure 5 .
Figure 5.Comparison of three class samples between DBI and NDBI, UI, EBBI of the study area.Figure 5. Comparison of three class samples between DBI and NDBI, UI, EBBI of the study area.

Figure 7 .
Figure 7.Comparison between DBSI, NDBaI, BI and NDSI for samples of four classes in the study area.

Table 2 .
Average and SD of DBSI, NDBaI, BI and NDSI of the city and 10 km buffer outside the city.Indices Mean (City) SD (City) Mean (10 km Buffer) SD (10 km Buffer) Difference (Buffer-City)

Figure 7 .
Figure 7.Comparison between DBSI, NDBaI, BI and NDSI for samples of four classes in the study area.

Table 1 .
Average and standard deviation of DBI, NDBI, UI and EBBI of the city at 10 km buffer outside the city.

Table 2 .
Average and SD of DBSI, NDBaI, BI and NDSI of the city and 10 km buffer outside the city.