Change detection

From AWF-Wiki
(Difference between revisions)
Jump to: navigation, search
Line 9: Line 9:
 
[[File:multi_temporal.png|250px]]
 
[[File:multi_temporal.png|250px]]
 
==== Principal component analysis ====
 
==== Principal component analysis ====
# Add the raster layers of the years 1992 (''tm_920526_mul.tif''), 2005 (''etm_050623_mul.tif'') and 2010 (''etm_000515_mul.tif'') into a [[QGIS]] project. It should be available in the [[Course data|course data]].
+
# Add the raster layers of the years 1992 (''tm_920526_mul.tif''), 2000 (''etm_000515_mul.tif'') and 2005 (''etm_050623_mul.tif'') into a [[QGIS]] project. It should be available in the [[Course data|course data]].
 
# Install PCA plugin under {{mitem|text=Plugins --> Manage and Install Plugins...}}.
 
# Install PCA plugin under {{mitem|text=Plugins --> Manage and Install Plugins...}}.
 
# Open {{mitem|text=PCA plugin}} [[File:Pca.png|25px]].
 
# Open {{mitem|text=PCA plugin}} [[File:Pca.png|25px]].
Line 15: Line 15:
 
#* Set {{button|text=Number of output Principal Components}} to 1.
 
#* Set {{button|text=Number of output Principal Components}} to 1.
 
#* Save the {{button|text=Output Raster File}} as ''pca1992_1''.
 
#* Save the {{button|text=Output Raster File}} as ''pca1992_1''.
#* Repeat this procedure for the raster files of 2005 and 2010 and adapt the name of each {{button|text=Output Raster File}}.
+
#* Repeat this procedure for the raster files of 2000 and 2005 and adapt the name of each {{button|text=Output Raster File}}.
 
# Merge the three output raster files with {{mitem|text=Toolbox --> GDAL/OGR --> [GDAL] Miscellaneous --> Merge}} and save {{button|text=Merged}} output as ''pca_merge''.
 
# Merge the three output raster files with {{mitem|text=Toolbox --> GDAL/OGR --> [GDAL] Miscellaneous --> Merge}} and save {{button|text=Merged}} output as ''pca_merge''.
  
 
==== Amplitude ====
 
==== Amplitude ====
# Add the raster layers of the years 1992 (''tm_920526_mul.tif''), 2005 (''etm_050623_mul.tif'') and 2010 (''etm_000515_mul.tif'') into a [[QGIS]] project. It should be available in the [[Course data|course data]].
+
# Add the raster layers of the years 1992 (''tm_920526_mul.tif''), 2000 (''etm_000515_mul.tif'') and 2005 (''etm_050623_mul.tif'') into a [[QGIS]] project. It should be available in the [[Course data|course data]].
 
# Open {{mitem|text=Toolbox --> OTB --> Miscellaneous --> Band Math}}.
 
# Open {{mitem|text=Toolbox --> OTB --> Miscellaneous --> Band Math}}.
# Calculate the amplitude for each raster layer (1992, 2005, 2010) with the use of the bands 5-4-3.
+
# Calculate the amplitude for each raster layer (1992, 2000, 2005) with the use of the bands 5-4-3.
 
#* For the year 1992, set ''tm_920526_mul.tif'' as {{button|text=Input image list}}.
 
#* For the year 1992, set ''tm_920526_mul.tif'' as {{button|text=Input image list}}.
 
#* Type <tt>sqrt(im1b4^2 + im1b5^2 + im1b3^2)</tt> as {{button|text=Expression}}.
 
#* Type <tt>sqrt(im1b4^2 + im1b5^2 + im1b3^2)</tt> as {{button|text=Expression}}.
 
#* Save the{{button|text=Output image}} as ''amplitude1992''.
 
#* Save the{{button|text=Output image}} as ''amplitude1992''.
#* Repeat this procedure for the raster files of 2005 and 2010 and adapt the name of each {{button|text=Output Raster File}}.
+
#* Repeat this procedure for the raster files of 2000 and 2005 and adapt the name of each {{button|text=Output Raster File}}.
 
# Merge the three output images with {{mitem|text=Toolbox --> GDAL/OGR --> [GDAL] Miscellaneous --> Merge}} and save {{button|text=Merged}} output as ''amplitude_merge''.
 
# Merge the three output images with {{mitem|text=Toolbox --> GDAL/OGR --> [GDAL] Miscellaneous --> Merge}} and save {{button|text=Merged}} output as ''amplitude_merge''.
  

Revision as of 16:30, 1 February 2016

Contents

Prerequisites of multi-temporal image analysis

  1. Acquisition of remote sensing data by the same (or similar) sensor recorded using same spatial resolution, viewing geometry, spectral bands, time of day.
  2. Optical calibrationcorrecting uninteresting differences due to illumination, atmospheric conditions, topography, sensor calibration. DNs are converted to physical units (surface reflectance).
  3. Clear definitions and classification scheme
  4. Exact spatial coregistration

Change detection techniques

Multi-temporal color composites

Multi temporal.png

Principal component analysis

  1. Add the raster layers of the years 1992 (tm_920526_mul.tif), 2000 (etm_000515_mul.tif) and 2005 (etm_050623_mul.tif) into a QGIS project. It should be available in the course data.
  2. Install PCA plugin under Plugins --> Manage and Install Plugins....
  3. Open PCA plugin Pca.png.
    • Set tm_920526_mul.tif as Input Raster File.
    • Set Number of output Principal Components to 1.
    • Save the Output Raster File as pca1992_1.
    • Repeat this procedure for the raster files of 2000 and 2005 and adapt the name of each Output Raster File.
  4. Merge the three output raster files with Toolbox --> GDAL/OGR --> [GDAL] Miscellaneous --> Merge and save Merged output as pca_merge.

Amplitude

  1. Add the raster layers of the years 1992 (tm_920526_mul.tif), 2000 (etm_000515_mul.tif) and 2005 (etm_050623_mul.tif) into a QGIS project. It should be available in the course data.
  2. Open Toolbox --> OTB --> Miscellaneous --> Band Math.
  3. Calculate the amplitude for each raster layer (1992, 2000, 2005) with the use of the bands 5-4-3.
    • For the year 1992, set tm_920526_mul.tif as Input image list.
    • Type sqrt(im1b4^2 + im1b5^2 + im1b3^2) as Expression.
    • Save theOutput image as amplitude1992.
    • Repeat this procedure for the raster files of 2000 and 2005 and adapt the name of each Output Raster File.
  4. Merge the three output images with Toolbox --> GDAL/OGR --> [GDAL] Miscellaneous --> Merge and save Merged output as amplitude_merge.

NDVI ratio

  1. Add the raster layers of the years 1992 (tm_920526_mul.tif) and 2005 (etm_050623_mul.tif) into a QGIS project. It should be available in the course data.
  2. Open Toolbox --> OTB --> Feature Extraction --> Radiometric indices.
    • Set tm_920526_mul.tif as Input Image.
    • Set Red Channel to 4 and NIR Channel to 5.
    • Set Available Radiometric Indices to ndvi.
    • Save the Output Image as ndvi1992.
    • Repeat this procedure for the raster file of 2005 and adapt the name of the Output Image to ndvi_2005.
  3. Calculate the ratio of both raster images with the Band Math tool (Toolbox --> OTB --> Miscellaneous).
    • Load ndvi_1992.tif and ndvi_2005 to the Input image list.
    • Type im1b1/im2b1 as Expression.
    • Save theOutput image as ndvi_ratio.
Personal tools
Namespaces

Variants
Actions
Navigation
Development
Toolbox
Print/export