Skip to main content

OBJECT-BASED BUILDINGS DETECTION & CHANGE ANALYSIS USING MULTI TEMPORAL HIGH RESOLUTION REMOTE SENSING DATA

With urbanization and urban sprawl, land-use and land-cover change studies became an area of keen interest for researchers. Cause of rapid urbanization is peoples’ responses to economic opportunities, and government policies. New housing societies and business model adopted overall globe has provided opportunities for new land uses.
Image by piviso from Pixabay 
“Image Classification is the process of sorting pixels into a finite number of individual classes, categories of data based on their DN (Digital Number) values”. 
With advent of modern high resolution era of satellite images, PIXEL based, Object based and many other machine learning algorithms have been developed.

Currently the prospects of a new classification concept, Object-Based Classification, are being investigated. Recent studies have proven the accuracy of object-based classification over traditional classifiers. The object-based classification approach relies on basic principle to use other information such as shape, texture and 3D information from DEM (Digital Elevation Model), and ancillary GIS data.

Change detection is defined as detecting spatial and spectral change in an object over a time span probably using temporal data sets. With advent of 21st century and rapid urbanization, accurate and timely updated information regarding LULC (land use land cover land cover) change and its spatial variations is significant for urban planning and management.

Our research domain for this project is based on Object-Based Classification and post-classification change detection technique.  Data set used for the study was multi-temporal high-spatial resolution imagery acquired by IKONOS sensor. Area of Interest was Lafayette-west from 2005 to 2018. Object-Based Classification techniques including, Multi-resolution segmentation was done on three images, then an object-based fuzzy approach was applied to implement the land cover classification separately for each of the three years i.e. 2005, 2013 & 2018. The object based and fuzzy approach was selected because it provides accurate and better results for change detection by incorporating spatial information and user defined expert knowledge. After rule set definition buildings of AOI were extracted and refined for object based post-classification comparison. Results were graphically and statistically compared for year 2005, 2013, & 2018. E-cognition ArcMap by Esri, and Erdas Imagine by Intergraph were used for our research project.

1.1  Objective
To make use of high resolution remote sensing data by incorporating object-based fuzzy classification and post-classification change detection to detect buildings of Lafayette-west for year 2005, 2013 & 2018.

1.2  Study Area
Lafayette is a city in and the parish seat of Lafayette Parish, Louisiana located along the Vermilion River in the southwestern part of the state. The city of Lafayette is the fourth-largest in the state, with a population of 127,657 according to 2015 U.S. Census estimates, It is the principal city of the Lafayette, Louisiana Metropolitan Statistical Area, with a 2015 estimated population of 490,488. The larger trade area or Combined Statistical Area of Lafayette-Opelousas-Morgan City CSA was 627,146 in 2015. Its nickname is The Hub City.


Figure 1.2 Study Area Zone-1

1.3  Data set used

High Resolution Imagery provided by
·         IKONOS imagery of  2013  and 2018
3D Information extracted using
·         DEM and DSM of study area 2005
Reference Data
·         2005 building footprints
Ancillary Data
·         Address data
·         Zoning Data
·         Roads Data


1.4  Software used

·         Ecognition by Trimble
·         ArcGIS by ESRI
·         Erdas Imagine by Intergraph
·         MS Excel by Microsoft




Image Source: eCognition Image
Image Source: ESRI Image






Image Source: Microsoft 











Image Source: Intergraph 















METHODOLOGY

Data was acquired, based on data quality no digital image processing was required. Image was further treated on ecognition software. Object based Classification was done on 2013 and 2018 imagery. Use of DEM and DSM provided subjected further ease in building detection.
Approach was made to distinguish shadow and buildings using NIR band threshold, Water was classified from shadow class by applying area threshold. Vegetation was classified by applying NDVI threshold. Buildings were then classified on the basis of height threshold and sub classified into child classes by applying area threshold as we have different area ranges in building class. Building class was classified by selecting samples and using nearest neighbor classifier. Most of the area at first was classified into building. Roads were separated from the imagery and building class by applying length threshold, vegetation using the same NDVI approach and child classes of building class was classified from buildings by applying area and length to width ration threshold.

Detailed step by step approach:
1.  Open the ecognition software.
2.  Click on “create new project” option.
3.  Import the “Study area image” from the workspace.
4.  Assign weights to NIR band
5.  Click on the layers and rename each layer. Layer 1 as blue, layer 2 as green, and 3rd as red and 4th NIR.
6.  Go to “Process” option on the top toolbar and open “ process tree” All the operations in ecognition are executed through the “ process tree”
7.  Right click in “process tree” window and open the option “Append new”. Uncheck the “Automatic option” and type the name “Segmentation” and click OK.
8.  The parent class “segmentation” has been added. Now click the ‘segmentation” option and select the “Insert child” option to add sub classes.
9.  Click the first segmentation algorithm” Multi Resolution Segmentation” and select the execute option.
10. Execute  next “Multiresolution, Spectral difference
Note that Execute Multi Resolution segmentation Scale 70 value -average size of objects Shape (0.4) weight -Geometric form of the objects Compactness (0.6) Value- The higher the value, the more compact image objects maybe.
11.  Make classes in class hierarchy
12.  Create NDVI
13.  Apply threshold
14.  Define other rule set such as NIR, DN value. Length/width ratio etc….
15.  Export Classes
16.  Data Cleansing in ArcMap


RESULTS
High resolution imageries of 2013 and 2018 were classified by object based fuzzy classification and with the aid of building footprint, address point and zoning data and change detection and comparison was done. Results generated are


Buildings Count
Entities
2005
2013
2018
Total Building
131
203
234
Medical-Related MR
2
2
2
Planned Residential Development
107
180
212
Single Family Residential
6
5
4
Single_Two and Multi Family Residential
16
16
16
Table 3.1 Statistics of Building Detection


Change Detection Maps

Conclusion

Object based fuzzy approach proved to be an efficient accurate and time saving image classification approach. Due to Urbanization buildings are increasing with the passing years. Using fuzzy rules by giving appropriate thresholds the accuracy of the classification increases. Change detection analysis showed growth in planned residential development both in 2013 and 2018. Small growth was also observed in single and multifamily residential area. Availability of very high spatial and temporal resolution remote sensing data facilitates mapping highly complex and diverse urban environments. Project concluded that
·         Incorporating high resolution aerial digital images and elevation data, and it’s processing using object-based image analysis through various rule set for mapping urban land covers and quantifying buildings is no doubt an efficient and accurate technique.
·         The Project presents and demonstrates an approach for formulating an optimal land covers classification rule set over small representative training urban area image, and its subsequent transfer to the same, and different multi-sensor, multi-temporal images.
·         There also was very rigorous selection of the object features for all of the land cover classes and iterative classification and rules updating
·         Shadows missclassified as buildings are corrected with additional segmentation and using the shape and context features re-applied to the buildings class only
·         This process picks up most of the shadow correctly, but if some shadow segments over the roofs still remain, the building can be split into parts
·         Overall, the height information helped capture most of the buildings, which otherwise are very diverse in their spectral and shape characteristics


References
  •       K. Leukert, A. Darwish and W. Reinhardt, “Urban Land-cover Classification: An Object-based Perspective,” 2nd Joint Workshop on Remote Sensing and Data Fusion over Urban Areas, Berlin, May 2003.
  •       P. Mather, Computer processing of remotely-sensed images, Chichester: Wiley, 1999. [3] B. Jähne, Digitale Bildverarbeitung, Berlin: Springer, 1993.
  •       R. Haralick, L. Shapiro, Computer and robot vision, Vol. 1, Reading, Mass.: Addison-Wesley, 1992.
  •       Definiens imaging, eCognition user guide [Online]. Available: http://www.definiens-imaging.com/documents/index.htm, April 2003.
  •       NIMA, Geospatial Standards and Specifications [Online]. Available: http://164.214.2.59/publications/specs/index.html, April 2003.
  •       Lillesand, T., Kiefer, R., “Remote Sensing and Image Interpretation”, New York, USA: John Wiley & Sons, 2000.
  •       E. Mendoza, R. Dirzo Deforestation in Lacandonia (southeast Mexico): evidence for the declaration of the northernmost tropical hot-spot
  •       Biodiversity and Conservation, 8 (1999), pp. 1621-1641

Comments

Popular posts from this blog

GDAL & OGR (Importance & Installation)

GDAL & OGR libraries allow you to perform various GIS operations like, Flow Direction or working with DEM few raster & vector analysis outside GIS environment. Image by  Gerd Altmann  from  Pixabay   What Knowledge one should have :  Basic Programming, Python What software you need : Python IDLE or command line, GDAL libraries files (explained later) Installation : 1- Once you have Python installed. 2- Start python IDLE to check which version of python you are using available at top as shown        Python IDLE 3- Install two files available at link below based on your version of python. Download  Link :  http://www.gisinternals.com/query.html?content=filelist&file=release-1500-gdal-2-1-2-mapserver-7-0-2.zip Install first :  gdal-201-1500-core.msi  and then:   GDAL-2.1.2.win32-py2.7.msi * Note specify the path C:\Python27\ArcGIS10.3\Lib\site-packages\ for GDAL-2.1.2.win32-py3.2.msi as shown in figure for (ArcGIS lib installed in python) Pyth

“Rule Based Classification” using eCognition

Availability of high resolution satellite imagery leads to a growing number of change detection techniques and algorithms,  research is done using remote sensing and GIS. In this tutorial, we build on the classification category feature set by introducing neighborhood relationships and topological functions. Secondly, we use relative elevation values and fuzzy rules for the classification systems. This study demonstrates that the rule-based classifier is a significantly better approach than the classical per-pixel classifiers & object based NN method . Image Source: NASA  Methodology 1.    Open the eCognition software. 2.    Click on “create new project” option. 3.    Import the “Study area image” from the workspace. 4.    Assign weights to NIR band 5.    Click on the layers and rename each layer. Layer 1 as blue, layer 2 as green, and 3rd as red and 4th NIR. 6.   Go to “Process” option on the top toolbar and open “ process tree” All the operations in ec

Working With Raster Calculator In ArcMap

Raster Calculator extracts information from a Raster Image, based on user requirements. Image by James Osborne from Pixabay Working with Raster Calculator To explore the highest elevation areas in your DEM raster e.g. select Spatial Analyst Tools à Map Algebra   Raster Calculator Tool  Raster Calculator Double click on the layer e.g. in my case it is  DEM_Layer or any raster of your choice,   to enter this into the ‘ search ’ window. Click on the “ > ” symbol e.g. and select a number less than the maximum elevation . Query can change based on your criterion.  This arithmetic raster operation will select all cells with values above the defined threshold. In the example below a threshold of 500m was used.  A new layer appears on your map. The majority of the map has a 0 value representing false (values below the threshold), and the colored region has a value of 1 representing true ( eleva

Learning Object Based Classification Using eCognition

Object Based Classification Using eCognition tutorials by GIS HIVE . ImageSource:PixaBay   Lab 01 Performing Image Classification On High Resolution Data on Erdas Imagine Lab 02 Introduction to eCognition Basics Lab 03 Image Segmentation using eCognition Lab 04 Object Based Classification Using Standard NN eCognition Lab 05 Rule Based Classification using  eCognition GIS HIVE intends to promote Geo Spatial Sciences awareness, for that purpose GIS HIVE has launched Youtube Channel . Here you can explore various Geospatial tutorials.

What is Radiometric Resolution Of Satellite Image

Radiometric Resolution  is defined as magnitude of EM ( Electromagnetic Energy ) recorded by a sensor Units : Measured in Bits Formula: Each bit records an exponent Power of 2. 1 Bit= 2 2 Bit=4 4 Bit =16 8 Bit=255 digital values Note: 0 represents Black tone and maximum like 255 represents highest tone. image source  www.satimagingcorp.com 1 Bit will be Black & white 16 Bit will be colored 32 Bit will be more colored and will carry greater energy levels Higher Radiometric Resolution much variation in colors reflected. Image to left is of  higher radiometric resolution than that of Landsat below. image source: geoinfo.amu.edu.pl

Generating Contours From DEM In ArcMap

Contours are a useful way to visualize topography and to create a map that can be printed in black and white, or that mimics the standard USGS topo maps. Image by TheAndrasBarta from Pixabay What you need: DEM data of your Area. Extract your area from online freely  available   DEM and select  Spatial Analyst Tools -> Surface -> Contour.  Select the inputs as follows, with a 100m or 50m contour interval: since 10 m might take time. and click  OK. Contour tool in ArcMap Output may be some thing like depending on color scheme. Go to properties and label elevation values for a good visualization. Image source:  en.wikipedia.org

GeoDa Free GIS Software

What is GeoDa : GeoDa provides analysis of various data formats like  shapefile, GeoJSON, KML, SQLite and table format (CSV, XLS and DBF). Image by  Pexels  from  Pixabay   How it is related to GIS: It provides Geo-visualization, user can perform certain spatial analysis.  Good things about GeoDa:   GeoDa is an Open Source software,  Along with Spatial analysis, simpler & complex statistical data analysis performed in excel (charting, graphing ) can be performed in GeoDa.  Most suited for: It is most suitable for Students, Researchers. Downloadable link :  http://geoda.software.informer.com/1.8/  imagesource: GitHub

Climate Change Not Myth But Reality

Northern Areas of Pakistan are blessed with beauty, we often hear people quoting these beautiful areas as “Switzerland of Pakistan”, “Heaven on Earth”, but these words proved there reality when for the first time I visited Northern Pakistan. It is worth mentioning that Pakistan is dependent on frozen hydrological resources. Many rivers such as River Swat, Chitral, Gilgit, Hunza, Shigar, Shyok, Indus, Shingo, Astor, and Jhelum are fed by runoff generated by melting of SC of Northern areas. Image by Aqsa kamran from Pixabay Image by Abdullah Shakoor2 from Pixabay This runoff supports agriculture . Agriculture of Pakistan is largely dependent on the snow melting phenomenon. Anomalies witnessed in glaciers and their melting patterns will have massive impact on hydrological resources because the country’s 70% of fresh water depends upon these frozen hydrological resources in the high mountains of Himalaya and Hindu Ku

How To View Satellite Image In ArcMap

Viewing Satellite Map in Arc Map " Are you  interested  to study satellite image, go to google and search for freely available online satellite images"  Image by  Daryl Govan  from  Pixabay   Download them, i have added tmrect.bil   named   imagery  here  Bands Added In ArcMap  Right Click tmrect.bil in ArcCatalog Select preview section and examine tmrect.bil What is tmrect: It is Landsat TM image with 5 bands . Right Click on tmrect and select properties All information of satellite image is available in Raster section. Explore tmrect.bil by expanding it you will see 5 bands in drop down Now Launch Arc Map Insert a new Data Frame by clicking on insert as shown Spatial Reference Dialogue  Now to Add Data Right Click Task Data Frame and click Add Data Select tmrect.bil a window “Unknown Spatial Reference” is displayed. Click  OK Arc map will display Raster Image as