INSCOHSTATS

Description

From a set of coregistered interferograms, INSCOHSTATS creates an output file consisting of the average, standard deviation, minimum, and maximum coherence. A time-ordered stack of coherence layers can optionally be appended to the output.

Parameters

Name	Type	Caption	Length	Value range
MFILE *	String	Input folder, or text file	1 - 192
DBIC	Integer	Selected Channel (Default channel 1)	0 -	Default: 1
STACK	String	Stack Data	0 - 192	yes\| no Default: no
FLSZ	Integer	Filter size (pixels, lines) ODD	0 -
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 list of input interferograms to be processed. Only the PIX file format (*pix) is supported. All files listed in the MFILE or within the input folder pointed to by MFILE must correspond to interferograms generated by INSRAW, INSTOPO, INSADJUST, INSMLOOK or INSMODGOLD. Unwrapped interferograms generated by INSUNWRAP are not 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 interferograms (one per line) with their full pathname, fore example: MFILE=r"C:/InSAR_project/RAW_neterograms_list.txt"

where the text file contains

C:\InSAR_project\RAW\ref20141028_dep20141215.pix C:\InSAR_project\RAW\ref20141121_dep20141215.pix C:\InSAR_project\RAW\ref20141215_dep20150108.pix C:\InSAR_project\RAW\ref20131009_dep20131126.pix C:\InSAR_project\RAW\ref20131126_dep20131220.pix ...

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

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

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

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

Note: For more information about using MFILE, including the requirements for a valid file, see Using an MFILE with a CATALYST Professional algorithm.

DBIC

The input raster layer from which to generate the coherence layer. Only complex channels of type C32R or C16S are supported.

If not specified, the first raster layer is used by default.
Only one channel can be specified.

This file must not already exist.

STACK

This option indicates whether to include the time-ordered coherence layers in the output PIX file. The layers are sorted in lexicographic order by reference date followed by dependent date.

If STACK="no", only the four statistics layers are returned (see the details section).

FLSZ

This parameter specifies the filter size to be applied to the data.

Blank will apply no filtering.
A single ODD integer will provide a boxcar filter of that size.
A pair of ODD integers will provide the boxcar filter size in the x and y directions respectively.

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

INSCOHSTATS computes coherence statistics for a given series of interferograms. The output file (filo) contains:

Coherence Average: the average of coherence at each pixel.
Coherence Standard Deviation: the standard deviation of each pixel.
Coherence Maximum: the maximum coherency at each pixel from all the interferogram.
Coherence Minimum: the minimum coherency at each pixel from all the interferogram.

Optionally, when STACK = "YES", the output file (FILO) will also contain the time series of coherence layers used to compute the statistics channels. The output coherence stack is time-ordered based upon the reference and dependent acquisition times as defined by the interferogram's metadata.

It is mandatory that all input interferograms have the same size (lines, columns) and channel type (C32R or C16S). The Interferograms must be of the same source and contain one the following metadata:

"InSARContent" = "RawInterferogram", (from INSRAW)
"InSARContent" = "DeformationInterferogram", (from INSTOPO)
"InSARContent" = "TopographicInterferogram", (from INSTOPO)
"InSARContent" = "ADJUSTED_DeformationInterferogram", (from INSADJUST)

INSCOHSTATS accepts wrapped interferograms generated by INSRAW, INSTOPO, INSADJUST, INSMODGOLD and INSMLOOK.

The NoData values (as defined by the "NO_DATA_VAULE" metadata tag) are excluded from the statistical computations.

If the center of the averaging window falls on a NODATA pixel, no statistics are calculated for that pixel. If an averaging window contains NODATA values, they are not used to calculate the statistics of that pixels.

Workflow

INSCOHSTASTS is generally used within the Persistent Scatterers Interferometry workflow (PSI) and requires all interferograms identified by MFILE to be co-registered to a common reference frame with the exact same number of lines and columns.

This is usually accomplished by first ingesting the data (SARINGESTAOI), co-registering the files (INSCOREG) using a single-reference file or the small baselines (SBAS) technique and then generating the interferograms (INSRAW) for co-registered pairs of reference and dependent data.

Examples

Generate the coherence statistics layers, with no filtering, using a list of RAW interferograms.


EASI>MFILE="mfile.txt"  				!List of raw interferograms files
EASI>DBIC=								!selecting channel, default 1st channel
EASI>STACK= "no"                 		!process default ""=no, only statistics
EASI>FLSZ=                       		!no filter
EASI>FILO="INSCOHSTATS_stats.pix"		!output containing the raster layers
EASI>run INSCOHSTATS

Generate a pix file with a 5x5 filter containing only the statistics channels.


EASI>MFILE="mfile.txt"  				!List of raw interferograms files
EASI>DBIC=								!selecting channel, default 1st channel
EASI>STACK= "no"                 		!process default "no", only statistics
EASI>FLSZ= 5                   			! 5 x 5 filter
EASI>FILO="INSCOHSTATS_stats.pix"		!output containing the raster layers
EASI>run INSCOHSTATS

Generate a pix file with a 5x7 filter containing only the statistics channels.

EASI>MFILE="mfile.txt"  				!List of raw interferograms files
EASI>DBIC=								!selecting channel, default 1st channel
EASI>STACK= "no"                 		!process default ""=no, only statistics
EASI>FLSZ=   5,7              			!rectangular 5x7 filter size
EASI>FILO="INSCOHSTATS_5x7_stats.pix" !output containing the raster layers
EASI>run INSCOHSTATS

Generate a pix file with a 5x5 filter containing the statistics channels and the time-ordered series of coherence layers.

EASI>MFILE="mfile.txt"  				!List of raw interferograms files
EASI>DBIC=								!selecting channel, default 1st channel
EASI>STACK= "yes"                 		!append temporally ordered coherence layers
EASI>FLSZ=   5,5              			!rectangular 5x5 filter size
EASI>FILO="INSCOHSTATS_5x5_stack.pix" !output containing the raster layers
EASI>run INSCOHSTATS

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								!selecting channel 2, could be any valid channel
EASI>STACK= "no"                 		!process default ""=no, only statistics
EASI>FLSZ=                       		!no filter
EASI>FILO="INSCOHSTATS_ch_2_def.pix"  !output containing the raster layers
EASI>run INSCOHSTATS

Environments	PYTHON :: EASI
Quick links	Description :: Parameters :: Parameter descriptions :: Details :: Examples :: Related