DESCRIPTION

Note: this command is no longer recommended for use. Please use i.pansharpen.

i.fusion.brovey performs a Brovey transformation using three multispectral and the panchromatic satellite image scene channels. Three new channels are calculated according to the formula:

                            DN_b1 
    DN_fused = ------------------------------- * DN_pan
                    DN_b1 + DN_b2 + DN_b3
The assignment of the channels depends on the satellite. The module's help text (above) suggests for LANDSAT-7, QuickBird and SPOT.

NOTES

The command changes temporarily to the high resolution of the panchromatic channels for creating the three output channels, then restores the previous region settings. The current region coordinates are respected. The three pan-sharpened output channels may be combined with d.rgb or r.composite.

EXAMPLE

Calculation of Brovey fusion map from North Carolina Landsat scene:
g.region raster=lsat7_2002_20 -p

# R, G, B composite at 28.5m 
d.rgb b=lsat7_2002_10 g=lsat7_2002_20 r=lsat7_2002_30

# Brovey fusion
i.fusion.brovey -l ms1=lsat7_2002_20 ms2=lsat7_2002_40 \
                   ms3=lsat7_2002_50 pan=lsat7_2002_80 \
                   outputprefix=brovey

# display at 14.25m
g.region raster=brovey.blue -p
d.rgb b=brovey.blue g=brovey.green r=brovey.red
Results:

 R, G, B composite of Landsat at 28.5m
R, G, B composite of Landsat at 28.5m
 R, G, B composite of Landsat Brovey fusion at 14.25m
R, G, B composite of Landsat Brovey fusion at 14.25m

Colors may be optionally optimized.

SEE ALSO

i.his.rgb, i.rgb.his, d.rgb, r.composite

REFERENCES

AUTHOR

Markus Neteler, ITC-irst, Italy

Last changed: $Date$