FUSEPCT

Description

Parameters

fusepct(fili, dbic, dbpct, filo, dbint, dboc, resample, ihsmodel)

Parameter descriptions

FILI

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

FILI and FILO can be the same PCIDSK file.

DBIC

Specifies the input pseudocolor image channel that defines the color (hue and saturation components) for the fused output color image.

The channel must be 8-bit.

DBPCT

Specifies the pseudocolor table segment (type 121) used for encoding the input image into red, green, and blue components.

FILO

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

FILI and FILO can be the same PCIDSK file.

DBINT

Specifies the input black-and-white intensity image that defines the intensity component of the fused output color image.

DBINT (InputIntensity) and DBOC (Output) channels must be all 8-bit or all 32-bit real.

DBINT (InputIntensity) must not be the same as DBOC (Output).

DBOC

Specifies three output RGB channels to store the results of the data fusion.

DBIC (InputColor) defines the color, and DBINT (InputIntensity) defines the intensity for the fused output color image.

The channels specified by DBINT (InputIntensity) and DBOC (Output) must be either all 8-bit or all 32-bit real.

The channels specified by this parameter must not be the same as the channel specified by DBINT (InputIntensity).

Duplicate output channels are not allowed.

RESAMPLE

Specifies the type of resampling method used on the input RGB color image.

Resampling extracts and interpolates the gray levels from the original pixel locations to corrected locations.

Three resampling modes are supported.

• NEAR (NEAREST): Nearest-neighbor interpolation identifies the gray level of the pixel closest to the specified input coordinates and assigns that value to the output coordinates.
• BILIN (BILINEAR): Bilinear Interpolation determines the gray level from the weighted average of the four closest pixels to the specified input coordinates and assigns that value to the output coordinates.
• CUBIC: Cubic Convolution determines the gray level from the weighted average of the 16 closest pixels to the specified input coordinates and assigns that value to the output coordinates.

The default is NEAR (NEAREST).

Resampling is not performed if FILI and FILO are the same.

IHSMODEL

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

Details

FUSEPCT performs data fusion of an input pseudocolor image on an input file with a black-and-white intensity image. The result is an output RGB color image with the same resolution as the intensity image. If the input and output files are different, the input pseudocolor image is resampled using the specified resampling method. Either the Hexcone or Cylinder IHS color model is used for data fusion.

FUSEPCT performs data fusion by running the following functions:

• REGPRO resamples and transfers the input pseudocolor image specified by DBIC (InputColor) to the first output channel specified by DBOC (Output).
• IIA transfers the pseudocolor table segment specified by DBPCT (InputColorPCT) from the input file to the output file.
• PCE creates the red, green, and blue components of a color image from the resampled input image and pseudocolor table.
• IHS performs an Intensity-Hue-Saturation transformation on the resampled color image specified by DBOC (Output).
• RGB performs a Red-Green-Blue transformation (the inverse of IHS) using DBINT (InputIntensity) for the intensity channel, and the output of IHS for the hue and saturation channels.

 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 pseudocolor table to output database.
 rem
 dbsl=dbpct \\ dbos=
 if (fili<>filo) run IIA
 rem
 rem Create RGB image from pseudocolor 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

Examples

A USGS Land Use/Land Cover theme map on channel 6 of irvine.pix is fused with a 10-meter resolution panchromatic Spot image on eltoro.pix. The result is a 10-meter resolution Spot image with the land use classes overlaid in translucent colors.

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

from pci.pcimod import pcimod

file='eltoro.pix'
pciop='ADD'
pcival=[3]                     #Add three 8-bit channels

pcimod(file,pciop,pcival)
        

Step 2: Enhance the input intensity image in eltoro.pix. Otherwise, the pseudocolor image colors overpower the intensity image.

from pci.lut import lut


fili  = "eltoro.pix"
filo  = "out_eltoro.pix"
dbic  = [1]         #Use input image channel 1
dblut = [2]         #Using LUT segment 2
dboc  = [1]         #Overwrite input channel
mask  = []          #entire image

lut (fili=fili, dbic=dbic, dblut=dblut, dboc=dboc, filo=filo)
        

Step 3: Use FUSEPCT to fuse the data. Nearest neighbor 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.

from pci.fusepct import fusepct

fili='irvine.pix'
dbic=[6]                   #Use channel 6
dbpct = [7]                #Use PCT segment 7
filo='eltoro.pix'
dbint = [1]
dboc = [2,3,4]
resample = 'NEAR'
ihsmodel = 'HEXCONE'

fusepct(fili,dbic,dbpct,filo,dbint,dboc, resample,ihsmodel)