FUSEPCT -- Data Fusion for Pseudocolour

Performs data fusion of a pseudocolour image with a black-and-white intensity image. The result is an output RGB colour image with the same resolution as the original B/W intensity image, but where the colour (hue and saturation) is derived from the resampled input pseudocolour image.

FUSEPCT is an EASI procedure which uses the REGPRO, PCE, IHS and RGB programs to perform data fusion.

Progress can be monitored. See MONITOR section.

PARAMETERS

FUSEPCT is controlled by the following global parameters:

Name     Prompt                                  Count     Type
FILI     Database Input File Name                1-64      Char
DBIC     Database Input Channel List             1         Int
DBPCT    Database Pseudocolour Segment           1         Int
FILO     Database Output File Name               1-64      Char
DBINT    Database Input Intensity Channel        1         Int
DBOC     Database Output Channel List            3         Int
RESAMPLE Resample Mode: NEAR/BILIN/CUBIC         4         Char
IHSMODEL IHS Model: CYLINDER/HEXCONE             8         Char

FILI

Specifies the name of the input PCIDSK file which contains the input pseudocolour image (DBIC) for data fusion.

 EASI>FILI="filename"

FILI and FILO can be the same PCIDSK file.

DBIC

Specifies one input channel on FILI for the pseudocolour image which defines the colour (hue and saturation components) for the fused output colour image.

 EASI>DBIC=channel

DBIC must be 8-bit.

DBPCT

Specifies the number of the pseudocolour table segment (type 121) on FILI used for encoding the input image (DBIC) into red, green and blue components.

 EASI>DBPCT=segment

FILO

Specifies the name of the output PCIDSK file which contains the input intensity image (DBINT) and output channels (DBOC) for data fusion.

 EASI>FILO="filename"

FILI and FILO can be the same PCIDSK file.

DBINT

Specifies the input black-and-white intensity image on FILO which defines the intensity component of the fused output colour image.

 EASI>DBINT=channel

Both DBINT and DBOC must be all 8-bit or all 32-bit real.
DBINT must not be the same as DBOC.

DBOC

Specifies three output RGB channels on FILO to store the results of data fusion. DBRGB on FILI defines the colour, and DBINT on FILO defines the intensity for the fused output colour image.

 EASI>DBOC=red,green,blue

Both DBINT and DBOC must be all 8-bit or all 32-bit real.
DBOC must not be the same as DBINT.
Duplicate output channels are NOT allowed.

RESAMPLE

Specifies the type of image resampling performed on the input RGB colour image (DBRGB). Three types of resampling are supported:

 EASI>RESAMPLE="NEAR"       | Nearest Neighbour
 EASI>RESAMPLE="BILIN"      | Bilinear Interpolation
 EASI>RESAMPLE="CUBIC"      | Cubic Convolution

Nearest Neighbour resampling is recommended if DBIC is a theme map or the result of classification.
Cubic Convolution resampling is recommended if DBIC is a elevation or other grey-scale image.
Resampling is not performed if FILI and FILO are the same.
Refer to the BACKGROUND section of the documentation for the REGPRO program for more information on resampling.

IHSMODEL

Specifies the type of IHS colour model to use. Refer to the IHS program documentation for more information on these two models.

 EASI>IHSMODEL="CYLINDER"    | Cylinder model
 EASI>IHSMODEL="HEXCONE"     | Hexcone model

DETAILS

FUSEPCT performs data fusion of an input pseudocolour image (DBIC and DBPCT) on an input file (FILI) with a black-and-white intensity image (DBINT) on an output file (FILO). The result is an output RGB colour image (DBOC) with the same resolution as the intensity image (DBINT). If the input and output files are different, the input pseudocolour image is resampled using the specified resampling method (RESAMPLE). Either the Hexcone or Cylinder IHS colour model (IHSMODEL) is used for data fusion.

FUSEPCT is an EASI procedure designed to STATUS and RUN like a PACE program. FUSEPCT performs data fusion by running the following programs:

REGPRO resamples and transfers the input pseudocolour image (DBIC) from FILI to the first output channel (DBOC) on FILO.
IIA transfers the pseudocolour table segment (DBPCT) from FILI to FILO.
PCE creates the red, green and blue components of a colour image from the resampled input image and pseudocolour table.
IHS performs an Intensity-Hue-Saturation transformation on the resampled colour image in DBOC.
RGB performs a Red-Green-Blue transformation (the inverse of IHS) using DBINT for the intensity channel, and the output of IHS for the hue and saturation channels.

The following EASI commands from the FUSEPCT procedure perform data fusion:

 rem
 rem Resample RGB image to output channels.
 rem
 ingeo=
 dbgc=
 outgeo=
 num0=dboc
 dboc(2)=
 order=1
 if (fili<>filo) run REGPRO
 rem
 rem Copy pseudocolour table to output database.
 rem
 dbsl=dbpct \ dbos=
 if (fili<>filo) run IIA
 rem
 rem Create RGB image from pseudocolour encoding.
 rem
 file=filo
 encoding="NORM"
 if (fili<>filo) dbic=dboc(1)
 if (fili<>filo) dbpct=lasc
 dboc=num0
 dbiw=
 run PCE
 rem
 rem Perform Intensity-Hue-Saturation transformation.
 rem
 dbic=dboc
 run IHS
 rem
 rem Perform Red-Green-Blue transformation.
 rem
 dbic=dbint(1),dboc(2),dboc(3)
 run RGB

MONITOR

Progress is monitored for four programs (REGPRO, PCE, IHS and RGB) run by the FUSEPCT procedure.

Program progress can be monitored by printing the percentage of completed processing in odometer fashion. A system parameter, MONITOR, controls this activity:

 EASI>MONITOR="ON"        | turn monitor ON  (default)
 EASI>MONITOR="OFF"       | turn monitor OFF (recommended if
                          | running in batch/background mode)

EXAMPLE

This example fuses the USGS Land Use/Land Cover theme map on channel 6 of irvine.pix with a 10 metre resolution panchromatic Spot image on eltoro.pix. The result is a 10 metre resolution Spot image with the land use classes overlaid in translucent colours.

Add three 8-bit channels to eltoro.pix for storing the result:

 EASI>FILE="eltoro.pix"
 EASI>PCIOP="ADD"
 EASI>PCIVAL=3
 EASI>RUN PCIMOD

Enhance the input intensity image on eltoro.pix. Otherwise, the pseudocolour image colours overpowers the intensity image.

 EASI>FILE="eltoro.pix"
 EASI>DBIC=1
 EASI>DBLUT=2
 EASI>DBOC=1
 EASI>MASK=
 EASI>RUN LUT

Use the FUSEPCT procedure to perform data fusion. Nearest neighbour resampling must be specified, since the input image is a theme map. In this case, the Hexcone model produces better results than the Cylinder model:

 EASI>FILI="irvine.pix"
 EASI>DBIC=6
 EASI>DBPCT=7
 EASI>FILO="eltoro.pix"
 EASI>DBINT=1
 EASI>DBOC=2,3,4
 EASI>RESAMPLE="NEAR"
 EASI>IHSMODEL="HEXCONE"
 EASI>RUN FUSEPCT

Transfer channels 2, 3 and 4 on eltoro.pix to the video display using ImageWorks or the IVI program.

 EASI>FILE="eltoro.pix"
 EASI>DBIC=2,3,4
 EASI>VDOC=1,2,3
 EASI>DBIW=
 EASI>VDOW=
 EASI>RUN IVI

About PCI Help Gateway