PSSARTSA

Description

PSSARTSA computes multi-temporal metrics form a co-registered stack of data. Fourteen metrics based upon amplitude statistics (see details section) are computed for the stack. The input values are automatically converted to amplitude when the metadata indicates units of intensity or decibel. If the input units are not specified, no conversion occurs. The option to append amplitude layers in temporal order is also provided.

Parameters

Name	Type	Caption	Length	Value range
MFILE *	String	Input folder, or text file	1 - 192
DBIC	Integer	Input raster channel	0 - 1	Default: 1
MASK	Integer	Area mask (bitmap or vector)	0 - 1
MASKFILE	String	Mask file	0 - 192
STACK	String	Stack Data	0 - 192	YES|NO Default: NO
FLSZ	Integer	Filter size (pixels, lines) ODD	0 - 2	1 -
FILO *	String	Database output file name	1 - 192
MONITOR	String	Monitor mode	0 - 3	ON, OFF Default: ON

* Required parameter

Parameter descriptions

MFILE

The name of a folder, or a text file that contains the input images to process. Only the PIX file format (*.pix) is supported.

All input files must be co-registered and have the same number of lines and columns.

When MFILE is a text file, it must contain a list of files (one per line) with their full pathname, for example: MFILE="C:/change_detection/sarimge_list.txt"

where the text file contains

C:/change_detection/sarimage_20201028.pix C:/change_detection/sarimage_20201121.pix C:/change_detection/sarimage_20201215.pix C:/change_detection/sarimage_20201009.pix C:/change_detection/sarimage_20201126.pix ...

Alternatively, an input folder can be specified: MFILE = "C:\change_detection"

The use of wildcards (*) is recommended to refine the search within the specified folder, for example:

MFILE = "C:\InSAR_project\RAW\*.pix"

MFILE = "C:\InSAR_project\RAW\ref*.pix"

DBIC

All input channels must be of the same size (lines, columns) and type.

Complex data are converted to amplitude before calculating the selected metrics.

MASK

The areas of the dataset to be included in the metric calculations. You can define inclusion areas using either bitmaps or polygon vectors. When a bitmap segment is specified, all pixels having a bitmap value of one define the area to mask. The bitmap segment must match the projection of the input. When a vector mask is specified, the masked pixels are those contained within the polygons.

MASKFILE

The name of the file containing the mask. If you do not specify a value for the MASK and MASKFILE parameters, no masking is performed. You cannot use the default value of MASK when MASKFILE is explicitly specified. MASKFILE is expected to be in a slant range projection of input files.

STACK

When STACK = "YES", all output amplitudes are time ordered and stack in a single PIX file.

When STACK = "NO", only the selected multi-temporal metric layers are produced.

FLSZ

FILO

Specifies the name of the PIX file to be generated.

This file must not already exist.

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

SAR MULTI-TEMPORAL METRICS

The multi-temporal metric is based on first-order statistics were derived from SAR intensity data.

1. Sum of Absolute changes (ASC) between adjacent temporal layers.
2. Mean: Represents the average value for each pixel.
3. Standard Deviation: Represents the standard deviation for each pixel.
4. Coefficient of Variation: Represents the ratio of the standard deviation to mean.
5. Median: Represents the median value (the central value in the ordered series) for each pixel.
6. Maximum Value: Represents the maximum value extracted from for each pixel.
7. Minimum Value: Represents the minimum value extracted from for each pixel.
8. Span: Maximum Value - Minimum Value
9. Maximum change: represents the maximum absolute change between consecutive acquisition dates.
10. Maximum change contribution to the ASC as a percentage %.
11. Index of Date of the Maximum change between layers.
12. Second Largest Change: represents the second largest absolute change between consecutive acquisition dates.
13. Second Largest Change contribution to the ASC in percentage %.
14. Index of Date of the Second Largest Change between adjacent layers.
15. to 14+N (optionally) temporally ordered layers.

Coregistration of the data layers must be performed by the INSCOREG module.

Pixels set as NoData values are not used in the calculations.

Required METADATA

All input files must contain the following file metadata: "Acquisition DateTime"

Examples

Generate the metrics only, with no averagin, for a list of SAR images contained in the text file.

EASI>MFILE="mfile.txt"				!List of raw interferograms files
EASI>DBIC=								!select channel, default 1st channel
EASI>MASK=								!select mask channel (Left Blank)
EASI>MASKFILE=							!no maskfile specified
EASI>STACK= "no"							!generate metrics only
EASI>FLSZ=								!no filter
EASI>FILO="PSSARTSA_stats.pix"		!output containing the raster layers
EASI>run PSSARTSA
      

Generate the metrics only without averaging for the areas defined by the second vector segment of the MASKFILE using the files in the specified directory.

EASI>MFILE="mfile.txt"				!directory containing SAR data sets
EASI>DBIC=								!select channel, default 1st channel
EASI>MASK= 2								!select mask channel (Left Blank)
EASI>MASKFILE= "MASK.pix"				!no maskfile specified
EASI>STACK= "no"							!process default ""=no, only statistics
EASI>FLSZ=								!no filter
EASI>FILO="PSSARTSA_Mask_stats.pix"	!output containing the raster layers
EASI>run PSSARTSA
      

Apply a 5x5 averaging filter before calculation the metrics..

EASI>MFILE="mfile.txt"				!directory containing SAR data sets
EASI>DBIC=								!select channel, default 1st channel
EASI>MASK=								!select mask channel (Left Blank)
EASI>MASKFILE=							!no maskfile specified
EASI>STACK= "no"							!process default "no", only statistics
EASI>FLSZ= 5								! 5 x 5 filter
EASI>FILO="PSSARTSA_stats.pix"		!output containing the raster layers
EASI>run PSSARTSA
      

Using the MFILE wildcard, generate a pix file with a 5x7 filter containing the metrics channels and the time-ordered amplitude layers.

EASI>MFILE="input/*_SLC.pix"			!wildcard identifying SAR data sets
EASI>DBIC=								!select channel, default 1st channel
EASI>MASK=								!select mask channel (Left Blank)
EASI>MASKFILE=							!no maskfile specified
EASI>STACK= "yes"							!append temporally ordered layers
EASI>FLSZ=   5,7							!rectangular 5x7 filter size
EASI>FILO="PSSARTSA_5x7_stats.pix"	!output containing the raster layers
EASI>run PSSARTSA
      

Generate a pix file containing the statistics channels from channel 2, could be from any valid channel.

EASI>MFILE="mfile.txt"				!List of raw interferograms files
EASI>DBIC= 2								!select channel 2, could be any valid channel
EASI>MASK=								!select mask channel (Left Blank)
EASI>MASKFILE=							!no maskfile specified
EASI>STACK= "no"							!process default ""=no, only statistics
EASI>FLSZ=								!no filter
EASI>FILO="PSSARTSA_ch_2_def.pix"		!output containing the raster layers
EASI>run PSSARTSA