ENVISAMImageDifferenceTask

This task performs a type of image change detection known as spectral angle difference. It determines the spectral similarity between the Time 1 image spectra and the Time 2 image spectra for every pixel. Each pixel spectrum is represented by a vector in space whose dimensionality is equal to the number of bands. The task calculates the angle between the two vectors. The smaller the angle, the more similar the two spectra.

Before running this task, ensure that the two input images have the same number of bands and the same spatial dimensions. Use the BANDS keyword to ENVISubsetRasterTask to define spectral subsets, if needed. Also use ENVIImageIntersectionTask to ensure that the two images have the same dimensions, projection, and pixel size.

Example

This example performs a SAM image difference change detection on two co-registered Sentinel-2 multispectral images. The images are before and after a flood event along the Mala River in Peru in 2017. They have the same number of bands, columns, and rows; no preprocessing is necessary.

The source files are available on the ENVI Resource DVD (in the change_detection folder) or from the ENVI Tutorials web page. Click the "ENVI Tutorial Data" hyperlink. Copy the images to your local drive, and update the file references in the example with the correct locations.

; Start the application

e = ENVI()

; Open input files

File1 = 'MalaRiver_2017-02-20.dat'

Raster1 = e.OpenRaster(File1)

File2 = 'MalaRiver_2017-03-12.dat'

Raster2 = e.OpenRaster(File2)

; Get the task from the catalog of ENVITasks

Task = ENVITask('SAMImageDifference')

; Define task inputs

Task.INPUT_RASTER1 = Raster1

Task.INPUT_RASTER2 = Raster2

; Run the task

Task.Execute

; Get the collection of data objects currently

; available in the Data Manager

DataColl = e.Data

; Add the output to the data manager

DataColl.Add, Task.OUTPUT_RASTER

; Display the result

View = e.GetView()

Layer = View.CreateLayer(Task.OUTPUT_RASTER)

See More Examples below.

Syntax

Result = ENVITask('SAMImageDifference')

Input properties (Set, Get): INPUT_RASTER1, INPUT_RASTER2, OUTPUT_RASTER_URI

Output properties (Get only): OUTPUT_RASTER

Properties marked as "Set" are those that you can set to specific values. You can also retrieve their current values any time. Properties marked as "Get" are those whose values you can retrieve but not set.

Methods

This task inherits the following methods from ENVITask:

Properties

This task inherits the following properties from ENVITask:

COMMUTE_ON_DOWNSAMPLE

INPUT_RASTER1 (required)

Specify a multiband raster representing the earliest time.

INPUT_RASTER2 (required)

Specify a multiband raster representing the latest time.

OUTPUT_RASTER

This is a reference to the output raster of filetype ENVI.

OUTPUT_RASTER_URI (optional)

Specify a string with the fully qualified filename and path of the associated OUTPUT_RASTER. If you do not specify this property, the associated OUTPUT_RASTER will not be created. To force the creation of a temporary file set the property to an exclamation symbol (!).

More Examples

This example performs a SAM image difference change detection between a Landsat TM image from 1985 and a Landsat 8 image from 2014. It shows the effects of urban sprawl in Las Vegas, Nevada over a 30-year period. Because the number of bands is different between the two images, the example uses ENVISubsetRasterTask to create a spectral subset of the Landsat 8 image. The example also uses ENVIImageIntersectionTask to ensure that the two images cover the same geographic extent and that they have the same spatial dimensions.

The source files are available on the ENVI Resource DVD (in the LandsatCaseStudies\TimeSeries folder) or from the ENVI Tutorials web page. Click the "ENVI Tutorial Data" hyperlink. Copy the images to your local drive, and update the file references in the example with the correct locations.

; Start the application

e = ENVI()

; Open a Landsat TM image

File1 = 'LasVegasMay1985.dat'

Raster1 = e.OpenRaster(File1)

; Open a Landsat 8 image

File2 = 'LasVegasMay2014.dat'

Raster2 = e.OpenRaster(File2)

; Create a spectral subset of the Landsat 8 image that

; omits the coastal aerosol band, so that the bands match

; those of the Landsat TM image (Raster1)

Time2SubsetTask = ENVITask('SubsetRaster')

Time2SubsetTask.INPUT_RASTER = Raster2

Time2SubsetTask.BANDS = [1,2,3,4,5,6]

Time2SubsetTask.Execute

; Run the image intersection task so that the spatial dimensions

; are the same between the two rasters

ImageIntersectionTask = ENVITask('ImageIntersection')

ImageIntersectionTask.INPUT_RASTER1 = Raster1

ImageIntersectionTask.INPUT_RASTER2 = Time2SubsetTask.OUTPUT_RASTER

ImageIntersectionTask.Execute

; Run the SAM image difference task

Task = ENVITask('SAMImageDifference')

Task.INPUT_RASTER1 = ImageIntersectionTask.OUTPUT_RASTER1

Task.INPUT_RASTER2 = ImageIntersectionTask.OUTPUT_RASTER2

Task.OUTPUT_RASTER_URI = e.GetTemporaryFilename()

Task.Execute

; Get the collection of data objects currently

; available in the Data Manager

DataColl = e.Data

; Add the output to the data manager

DataColl.Add, Task.OUTPUT_RASTER

; Display the result

View = e.GetView()

Layer = View.CreateLayer(Task.OUTPUT_RASTER)

View.Zoom, 0.4

Version History

ENVI 5.5

Introduced

API Version

3.6