DIEST

Description

DIEST combines multiple layers of spatial data to derive a favorability model for predicting a geological event (mineral potential, landslide, etc.). Available algorithms include probability estimation, certainty factor estimation, and fuzzy membership estimation. DIGRP must be run before DIEST.

Parameters

Name	Type	Caption	Length	Value range
FILE *	String	Input file name	1 - 192
UNICON *	Integer	Uni-condition region IDs channel	1 - 1
DATFIL *	String	ASCII data file name	1 - 192
ALGO *	String	Estimation algorithm	1 - 16	See Description Default: CF
DBOC *	Integer	Output estimated favorability map channel	1 - 1
REPORT	String	Report mode	0 - 192	Quick links
MONITOR	String	Monitor mode	0 - 3	ON, OFF Default: ON

* Required parameter

Parameter descriptions

FILE

Specifies the name of the PCIDSK file containing the uni-condition regions channel.

UNICON

Specifies the channel containing the uni-condition region ID(s), produced by DIGRP.

DATFIL

Specifies the name of the ASCII data file that contains the data layer statistics, as produced by DIGRP.

ALGO

Specifies the estimation algorithm to be used.

DBOC

Specifies the output channel to receive the estimated favorability map (which is scaled to [0, 255]).

REPORT

Specifies where to direct the generated report.

Available options are:

• TERM: generates a report on the terminal (default)
• DISK: generates a report on file "IMPRPT.LST"
• OFF: turns off report generation
• <filename>: appends a report to the specified file

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

DIEST provides various algorithms for estimating the favorability of occurrence, by combining evidence from multiple input layers. All the statistics of input layers are gathered and saved in an ASCII data file by DIGRP in a preprocessing step. DIEST does not need to read data from the original input channels.

The favorability measure is estimated for each unique condition region rather than for each pixel. Because the number of unique condition regions is much smaller than the number of pixels in the image, this process significantly reduces computation.

The computed favorabilities for all the uni-condition regions are saved in a file specified by REPORT.

In addition, a raster image of favorabilities is generated where each pixel contains the favorability measure mapped to the range of 0-255 (the lowest favorability is mapped to 0 and the highest is mapped to 255). The mapping function is similar to histogram equalization; you may expect approximately the same number of pixels at each gray level. For example, to obtain the top 10% most favorable areas, threshold the output image at 255*90%. The pixels above this threshold should fall in the top 10% category.

Example

The following example assumes that the first 7 channels of land.pix are the data layers to be integrated. Segment 2 contains a bitmap of the known occurrence areas. Segment 3 contains a bitmap of the surveyed areas. First, by running DIGRP, the data layers are grouped into a uni-condition image and saved in channel 17. Some statistics of uni-condition regions are saved in the ASCII file land.dat. Then, by running DIEST, a favorability map is generated using the Certainty Factor algorithm (CF) and saved in channel 10.

EASI>FILE="land.pix"
EASI>DBIC=1,2,3,4,5,6,7
EASI>KOSN=2
EASI>SASN=3
EASI>UNICON=17
EASI>DATFIL="land.dat"

EASI>run digrp

EASI>ALGO="CF"
EASI>DBOC=10

EASI>run diest

Algorithm

For the Fuzzy Logic Gamma operator, the Gamma value must be 0 < Gamma < 1. Otherwise, a default Gamma value of 0.95 is used.