GCD

Description

Generates a difference image from two input images.

Parameters

Name	Type	Caption	Length	Value range
FILI_IN *	String	Input file (I1)	1 - 192
FILI_REF *	String	Input file (I2)	1 - 192
FILO_DIF *	String	Output file (GCD Image)	1 - 192
DBIC_REF *	Integer	Input reference channel/layer	1 - 1	1 -
DBIC_IN *	Integer	Input image channel/layer	1 - 1	1 -
DBOC_DIF *	Integer	Output difference channel/layer	1 - 1	1 -
MONITOR	String	Monitor mode	0 - 3	ON, OFF Default: ON

* Required parameter

Parameter descriptions

FILI_IN

Specifies the PCIDSK file that contains the input image data channel. The image must have the same dimensions and datatype (8-bit unsigned) as FILI_REF.

FILI_REF

Specifies the PCIDSK file that contains the image data to use as a reference. The reference image must have the same dimensions and datatype (8-bit unsigned) as FILI_IN.

FILO_DIF

Specifies the output PCIDSK file which will receive the GCD image.

The file must already exist and must have the same dimensions (pixels and lines) as FILI_IN, and FILI_REF.

DBIC_REF

Specifies the reference image channel or layer. This reference image channel must be the same size and datatype (8-bit unsigned) as DBIC_IN (Input).

DBIC_IN

Specifies the input image channel or layer. This image channel must be the same size and datatype (8-bit unsigned) as DBIC_REF (Input2).

DBOC_DIF

Specifies the image channel or layer to receive the GCD image data. This image channel must have the same dimensions as DBIC_IN (Input) and DBIC_REF (Input2), and must be 32-bit real (32R).

MONITOR

The program progress can be monitored by printing the percentage of processing completed. A system parameter, MONITOR, controls this activity.

Available options are:

• ON: Turns the monitor ON (default)
• OFF: Turns the monitor OFF (recommended if running in batch or background mode)

Details

GCD performs a Change Detection comparison of two input images, with the output going to a 32-bit real image of the same size. Input is restricted to 8-bit unsigned images, of the same X/Y (pixels and lines) dimensions.

GCD performs its change detection on a global scale, meaning it calculates difference on the whole image. The least squares minimization parameters, B0 and B1, are calculated once on the whole image (globally), then using the following formula:

FILI_DIF = Minimize || FILI_IN - (B1 * FILI_REF + B0) ||

The difference image is calculated.

Example

Determine the global change between channel 2 of the multispectral image database (ms.pix), and channel 1 of the panchromatic image database (pan.pix). Save the real GCD output image to the 1st channel (32R) of GCD.pix.

EASI>FILI_IN    = "ms.pix"
EASI>FILI_REF   = "pan.pix"
EASI>FILO_DIF   = "gcd.pix"
EASI>DBIC_IN    = 2
EASI>DBIC_REF   = 1
EASI>DBOC_DIF   = 1

EASI> run GCD

Algorithm

Global Change Detection obtains a difference image from two input images, I1(Input) and I2 (Input2), using the the following formula:

I1 -  (b1 * I2 + b0)            (1)

The coefficients b0 and b1 are calculated using a Least Square Minimization, namely:

Min || I1 -  (b1 * I2 + b0) ||  (2)

The coefficients b0 and b1 are calculated from the entire image databases, I1 (Input) and I2 (Input2), before the equation (1) is applied, to obtain the difference image.