Multisensor Image Fusion in Remote Sensing
Multisensor Image Fusion in Remote Sensing
International Journal of
Remote Sensing
Publication details, including instructions for
authors and subscription information:
[Link]
To cite this article: C. Pohl & J. L. Van Genderen (1998): Review article Multisensor
image fusion in remote sensing: Concepts, methods and applications, International
Journal of Remote Sensing, 19:5, 823-854
This article may be used for research, teaching, and private study
purposes. Any substantial or systematic reproduction, redistribution,
reselling, loan, sub-licensing, systematic supply, or distribution in any form
to anyone is expressly forbidden.
The publisher does not give any warranty express or implied or make any
representation that the contents will be complete or accurate or up to
date. The accuracy of any instructions, formulae, and drug doses should
be independently verified with primary sources. The publisher shall not
be liable for any loss, actions, claims, proceedings, demand, or costs or
damages whatsoever or howsoever caused arising directly or indirectly in
connection with or arising out of the use of this material.
int. j. remote sensing, 1998 , vol. 19 , no. 5 , 823± 854
Review article
C. POHL
Western European Union± Satellite Centre (WEUSC), Research Division,
P.O. Box 511, 28850 Torrejo n de Ardoz (Madrid ), Spain
Downloaded by [North Carolina State University] at 01:55 14 January 2013
1. Introduction
Earth observation satellites provide data covering di erent portions of the electro-
magnetic spectrum at di erent spatial, temporal and spectral resolutions. For the
full exploitation of increasingly sophisticated multisource data, advanced analytical
or numerical data fusion techniques are being developed (Shen 1990 ). Fused images
may provide increased interpretation capabilities and more reliable results since data
with di erent characteristics are combined. The images vary in spectral, spatial and
temporal resolution and therefore give a more complete view of the observed objects.
It is the aim of image fusion to integrate di erent data in order to obtain
more information than can be derived from each of the single sensor data alone
(`1+ 1 = 3’). A good example is the fusion of images acquired by sensors sensitive to
visible/ infrared ( VIR) with data from active synthetic aperture radar (SAR). The
information contained in VIR imagery depends on the multispectral re¯ ectivity of
the target illuminated by sun light. SAR image intensities depend on the charac-
teristics of the illuminated surface target as well as on the signal itself. The fusion of
these disparate data contribute to the understanding of the objects observed.
Image fusion has many aspects to be looked at. Before being able to implement
and use an image fusion approach some of its questions that need to be answered
by the user include:
Ð What is the objective /application of the user?
Ð Which types of data are the most useful for meeting these needs?
0143± 1161/98 $12.00 Ñ 1998 Taylor & Francis Ltd
824 C. Pohl and J. L . van Genderen
Ð Which is the `best’ technique of fusing these data types for that particular
application?
Ð What are the necessary pre-processing steps involved?
Ð Which combination of the data is the most successful?
Ð platform;
Ð imaging geometry of optical and radar satellites;
Ð spectral, spatial and temporal resolution ( Pohl and Genderen 1993).
The availability of speci® c data plays an important role too. It depends on the
satellite coverage, operational aspects of the space agency running the satellite,
atmospheric constraints such as cloud cover, ® nancial issues, etc. ( Pohl 1996).
The next step is the choice of a suitable fusion level. The pre-processing steps
are depending on this. In case of pixel based image fusion the geocoding is of vital
importance. Related to the geometric correction of the data, details such as the
Ð geometric model
Ð ground control points (number, distribution, accuracy)
Ð digital elevation model
Ð resampling method etc.
fused data. Finally, bene® ts and limitations of image fusion are summarized in the
concluding section.
common interpretation and resolve di erences and furnish a better understanding
of the observed objects (Shen 1990).
This literature review mainly tackles pixel based image fusion. Techniques used
to fuse data at pixel level are described in § 4. In the literature a large number of
terms can be found. Some terms de® ne di erent approaches to image fusion, whilst
others can be used equivalently. Keys et al. ( 1990 ) and Franklin and Blodgett ( 1993 )
describe the fusion as the computation of three new values for a pixel based on the
known relation between the input data for the location in the image. The input
images are di erent in terms of spatial, spectral or temporal characteristics. Other
expressions accepted in this context are image merging (Carper et al. 1990), image
integration ( Welch and Ehlers 1988) and multi-sensor data f usion ( Franklin and
Blodgett 1993 ). A broader view of the term is presented by Mangolini ( 1994). He
Downloaded by [North Carolina State University] at 01:55 14 January 2013
describes data fusion as group of methods and approaches using multisource data
of di erent nature to increase the quality of information contained in the data. Often,
in the case of data fusion, not only remote sensing images are fused but further
ancillary data (e.g., topographic maps, GPS coordinates, geophysical information,
etc.) contribute to the resulting image ( Harris and Murray 1989 ). Data integration
comprises algorithms for image fusion too ( Nandhakumar 1990). This term is often
used in the context of geographical information systems (GIS) ( Ehlers 1993). Other
scientists evaluate multisensor image data in the context of combined ( Lichtenegger
1991 ), coincident (Crist 1984), complementary ( Koopmans and Forero 1993),
composited ( Daily et al. 1979 ) or co-registered ( Rebillard and Nguyen 1982 ) data
analysis. In these cases not always an alteration of the digital numbers amongst the
di erent image channels is involved. A simple overlay of multi-source data in a Red-
Green-Blue ( RGB) colour space integrates the data set. The replacement of one of
the three channels or parts of it with an image from another data source is called
substitution (Suits et al. 1988 ). Referring to a di erent level of image processing are
the words information f usion (Shufelt and McKeown 1990 ). Here the images are
already interpreted to reach the information or knowledge level before being fused.
A very wide ® eld of applications and approaches in image fusion are summarised by
synergy (Genderen et al. 1994 ) or synergism ( Harris and Graham 1976) of remote
sensing data. It requires the input of data that provide complementary rather than
redundant information.
Apart from the three levels at which fusion can be performed image fusion can
be applied to various types of data sets:
Ð single sensorÐ temporal ( Weydahl 1993 ), e.g., SAR multitemporal for change
detection
Ð multi-sensorÐ temporal ( Pohl and Genderen 1995 ), e.g., VIR/ SAR image
mapping
Ð single sensorÐ spatial (Cliche et al. 1985 ), e.g., high/ low resolution
panchromatic/ multi-spectral SPOT
Ð multi-sensorÐ spatial (Chavez et al. 1991), e.g., high/ low resolution
SPOT/ Landsat
Ð single dataÐ multi-sensor (Guyenne 1995 ), e.g., ERS-1/ ERS-2
Ð remote sensing data with ancillary data ( Janssen et al. 1990), e.g., image with
topographic maps, (the references given are not exhaustive but meant as
example).
Having introduced the concept of image fusion, the following section outlines some
goals of image fusion.
Multisensor image f usion in remote sensing 827
conventional image to map registration are the di erences in appearance of the
features used as control or tie points. The positioning of points is facilitated if they
are localised in similar `views’ (remote sensing data). Of advantage is an integrative
recti® cation approach which iteratively improves the registration accuracy ( Ehlers
1991 ). Multisensor stereo mapping requires accurately co-registered input data and
shows improved results in comparison with image-to-map registration ( Welch et al.
1985, 1990 ).
used. Some constraints exist depending on radiometric di erences of the images used
as stereo pair ( Domik et al. 1986, Bloom et al. 1988, Welch et al. 1990, Buchroithner
et al. 1992, Raggam et al. 1994, Toutin and Rivard 1995 ).
might have occurred in the area observed. T emporal image f usion is applicable to
images from the same sensor as well as to multiple sensor data. It is nearly impossible
to acquire multisensor data simultaneously so that fusion of data from di erent
sensors mostly includes a temporal factor. An overview on image enhancement
techniques for change detection from multisensor data including image di erencing,
ratioing and PCA is given in Mouat et al. ( 1993). The techniques are described in
§ 4. There is a necessity of correcting the input imagery for radiometric distortions
which occur especially in VIR data, i.e., atmospheric and illumination correction in
order to create compatible data sets. The images to be fused should be acquired at
similar seasons to account for seasonal changes which might in¯ uence the change
detection capabilities of this approach.
Multitemporal SAR data provides good potential for change detection analysis
Downloaded by [North Carolina State University] at 01:55 14 January 2013
because of its all-weather capability. The evaluation of single SAR imagery is di cult
due to the ambiguities contained in SAR images ( De Groof et al. 1992, Weydahl
1992, Kattenborn et al. 1993, Kohl et al. 1994 ).
Figure 2. Processing ¯ ow chart for pixel based image fusion ( Pohl 1996 ).
has to be pointed out that this is a major element of pixel level image fusion. These
fusion techniques are very sensitive to misregistration.
Subsequently, the data can be fused using one of the fusion techniques described
in this review. In some cases, especially if image data of very di erent spatial
resolution is involved, the resampling of the low resolution data to the pixel size of
the high resolution image might causes a blocky appearance of the data. Therefore
a smoothing ® lter can be applied before actually fusing the images (Chavez 1987).
The resulting image map can be further evaluated and interpreted related to the
desired application.
In general, the techniques can be grouped into two classes: ( 1 ) Colour related
techniques, and ( 2 ) Statistical/ numerical methods. The ® rst comprises the colour
composition of three image channels in the RGB colour space as well as more
sophisticated colour transformations, e.g., IHS and HSV. Statistical approaches are
developed on the basis of channel statistics including correlation and ® lters.
Techniques like PCA and regression belong to this group. The numerical methods
follow arithmetic operations such as image di erencing and ratios but also adding
of a channel to other image bands. A sophisticated numerical approach uses wavelets
in a multiresolution environment. The next sections describe these techniques in
more detail.
between any two of three bands. Other considerations related to OIF can be found
in Chavez et al. ( 1984 ), Keys et al. ( 1990 ) and Ye sou et al. ( 1993 a).
Another approach is to select the bands which are the most suitable to a certain
application. This requires a priori knowledge by the user (She eld 1985 ). Kaufmann
and Buchroithner ( 1994) suggest to select the three bands with the highest variance.
The principal component analysis is another solution to reduce the number of
channels containing the majority of the variance for the purpose of image fusion
( Ye sou et al. 1993 a). Others provide a certain band combination based on an
algorithm which takes into account the statistics of the scene including correlations
between image channels (She eld 1985 ).
blue. Together they form a colour composite that can be displayed with conventional
media, e.g., cathode ray tube (CRT). The colour composite facilitates the interpreta-
tion of multichannel image data due to the variations in colours based on the values
in the single channels. The grey scale value, used to select a set of red, green, and
blue brightness’, is stored in a look up table ( LUT ) that are the voltages sent to the
display tube. Operations on the LUT and the histogram of the image data can
enhance the colour composite for visual interpretation ( Russ 1995 ).
The possibilities of varying the composite are manifold. Depending on the selec-
tion of the input image channels the fused data will show di erent features. Very
important for the colour composite is the distribution of the available 0± 255 grey
values to the range of the data. It might be of advantage to invert input channels
before combining them in the RGB display with other data depending on the objects
Downloaded by [North Carolina State University] at 01:55 14 January 2013
A B
Jupp 1990 ).
1 1 1
AB AB
Ó 3 Ó 3 Ó 3
I R
1 1 2
n1 = Õ G (a)
Ó 6 Ó 6 Ó 6
n2 B ( 2)
1 1
Õ 0
Ó 2 Ó 2
AB
n2
H = tan Õ S = Ó n1 + n2
1 2 2
(b) (c)
n1
There are two ways of applying the IHS technique in image fusion: direct and
substitutional. The ® rst refers to the transformation of three image channels assigned
Multisensor image f usion in remote sensing 833
to I , H and S ( Rast et al. 1991). The second transforms three channels of the data
set representing RGB into the IHS colour space which separates the colour aspects
in its average brightness (intensity). This corresponds to the surface roughness, its
dominant wavelength contribution ( hue) and its purity (saturation) (Gillespie et al.
1986, Carper et al. 1990). Both the hue and the saturation in this case are related to
the surface re¯ ectivity or composition (Grasso 1993). Then, one of the components
is replaced by a fourth image channel which is to be integrated. In many published
studies the channel that replaced one of the IHS components is contrast stretched
to match the latter. A reverse transformation from IHS to RGB as presented in
equation ( 3 ) converts the data into its original image space to obtain the fused image
( Hinse and Proulx 1995 b).
The IHS technique has become a standard procedure in image analysis. It serves
Downloaded by [North Carolina State University] at 01:55 14 January 2013
colour enhancement of highly correlated data (Gillespie et al. 1986), feature enhance-
ment (Daily 1983), the improvement of spatial resolution (Welch and Ehlers 1987,
A B
Carper et al. 1990 ) and the fusion of disparate data sets ( Harris et al. 1990, Ehlers 1991).
1 1 1
AB AB
Ó 3 Ó 6 Ó 2
R I
1 1 1
G = Õ n1 ( 3)
Ó 3 Ó 6 Ó 2
B n2
1 2
Õ 0
Ó 3 Ó 6
The use of IHS technique in image fusion is manifold, but based on one principle:
the replacement of one of the three components (I , H or S ) of one data set with
another image. Most commonly the intensity channel is substituted. Replacing the
intensityÐ the sum of the bandsÐ by a higher spatial resolution value and reversing
the IHS tranformation leads to composite bands (Chavez et al. 1991 ). These are
linear combinations of the original (resampled ) multispectral bands and the higher
resolution panchromatic band (Campbell 1993).
A variation of the IHS fusion method applies a stretch to the hue saturation
components before they are combined and transformed back to RGB ( Zobrist et al.
1979 ). This is called colour contrast stretching (Gillepsie et al. 1986). The IHS
transformation can be performed either in one or in two steps. The two step approach
includes the possibility of contrast stretching the individual I , H and S channels. It
has the advantage of resulting in colour enhanced fused imagery ( Ehlers 1991). More
results using IHS image fusion are reported by Rast et al. ( 1991 ), Jutz and Chorowicz
( 1993 ), Koopmans and Richetti ( 1993), Oprescu et al. ( 1994), Smara et al. ( 1996 ),
and Yildimi et al. ( 1996 ). A closely related colour system to IHS (sometimes also
called HSI) is the HSV: hue, saturation and value ( Russ 1995 ). An example of HSV
image fusion was presented by Chiesa and Tyler ( 1990).
4.2.3. L uminance-chrominance
Another colour encoding system called YIQ has a straightforward transformation
from RGB with no loss of information. Y , the luminance , is just the brightness of a
panchromatic monochrome image. It combines the red, green, and blue signals in
proportion to the human eye’s sensitivity to them. The I and Q components of the
colour are chosen for compatibility with the hardware used. The I is essentially red
834 C. Pohl and J. L . van Genderen
minus cyan, while Q is magenta minus green. The relation between YIQ and RGB
is shown in equation ( 4 a , b ):
AB A BA B
Y 0´299 0´587 0´114 R
AB A BA B
( 4)
R 1´000 0´956 0´621 Y
Since the Y , I and Q components are less correlated than the RGB ones, this
transformation o ers better possibilities for enhancing an image. This was shown in
an example by Guo and Pinliang ( 1989).
cases the resulting di erence image contains negative values. Therefore a constant
has to be added to produce positive digital numbers. However, di erences do not
always refer to changes since other factors, like di erences in illumination, atmo-
spheric conditions, sensor calibration, ground moisture conditions and registration
of the two images, can lead to di erences in radiance values. In ratioing, two images
from di erent dates are divided, band by band if the image data have more than
one channel. If the intensity of the re¯ ected energy is nearly the same in each image
then the ratio image pixel is one, it indicates no change. The critical part of this
method is selecting appropriate threshold values in the lower and upper tails of the
distribution representing change pixel values. In this respect the normalisation of
the data is of advantage as indicated in equation ( 6 ), (Gri ths 1988):
Downloaded by [North Carolina State University] at 01:55 14 January 2013
XS3 Õ XS 2 TM 4 Õ TM 3
Õ +C (C = 128 for positive values) ( 6)
XS3 + XS2 TM 4 + TM 3
A ratio for spatial enhancement is summarised by equation ( 7). The aim of this
method is to maintain the radiometric integrity of the data while increasing the
spatial resolution.
DNXS(i)
DNHybridXS(j) = DNPAN ( 7)
DNSynPAN
where DNHybridXS(i) - i th band of the fused high resolution image, DNPAN - correspond-
ing pixel in high resolution input PAN image, DNXS(i) - super-pixel in i th band of
input low resolution XS image, DNSynPAN - corresponding pixel in low resolution
synthetic PAN image, created from low resolution multispectral bands that overlap
the spectral response of the input high resolution PAN, (Munechika et al. 1993 ).
That is also the aim of the Brovey Transform, named after its author. It is a
formula that normalises multispectral bands used for a RGB display, and multiplies
the result by any other desired data to add the intensity or brightness component
to the image. The algorithm is shown in equation ( 8) where D N fused means the DN
of the resulting fused image produced from the input data in n multispectral bands
multiplied by the high resolution image DNhighres .
DNb1
DNfused = DNhighres ( 8)
DNb1 + DNb2 + DNbn
4.3.2. Principal component analysis
The PCA is useful for image encoding, image data compression, image enhance-
ment, digital change detection, multitemporal dimensionality and image fusion. It is
a statistical technique that transforms a multivariate data set of intercorrelated
variables into a data set of new un-correlated linear combinations of the original
variables. It generates a new set of axes which are orthogonal.
The approach for the computation of the principal components ( PCs) comprises
the calculation of:
1. Covariance (unstandardised PCA) or correlation (standardised PCA) matrix
2. Eigenvalues, -vectors
3. PCs
An inverse PCA transforms the combined data back to the original image space.
The use of the correlation matrix implies a scaling of the axes so that the features
836 C. Pohl and J. L . van Genderen
receive a unit variance. It prevents certain features from dominating the image
because of their large digital numbers. The signal-to-noise ratio ( SNR) is signi® cantly
improved applying the standardised PCA (Singh and Harrison 1985, Shettigara
1992 ). Better results are obtained if the statistics are derived from the whole study
area rather than from a subset area ( Fung and LeDrew 1987). The PCA technique
can also be found under the expression Karhunen L oeve approach ( Zobrist et al.
1979 ).
Two types of PCA can be performed: selective or standard. The latter uses all
available bands of the input image, e.g., TM 1± 7, the selective PCA uses only a
selection of bands which are chosen based on a priori knowledge or application
purposes ( Ye sou et al. 1993 b). In case of TM the ® rst three PCs contain 98± 99 per
cent of the variance and therefore are su cient to represent the information.
Downloaded by [North Carolina State University] at 01:55 14 January 2013
with band addition. The high spatial resolution image is ® ltered with a small high-
pass-® lter resulting in the high frequency part of the data which is related to the
spatial information. This is pixel wise added to the low resolution bands ( Tauch and
KaÈhler 1988, Shettigara 1992, Jutz and Chorowicz 1993). Nevertheless, the HPF
method has limitations in passing on important textural information from the high
resolution band to the low resolution data (Shettigara 1992).
of an image channel that can be replaced by another image channel. This method
is called regression variable substitution ( RVS).
To achieve the e ect of fusion the replacement vector should account for a
signi® cant amount of variance or information in the original multi-variate data set.
The method can be applied to spatially enhance data or for change detection with
the assumption that pixels acquired at time one are a linear function of another set
of pixels received at time two. Using the predicted value obtained from the least-
square regression, the di erence image is the regression value-pixel of time one (Singh
1989, Shettigara 1992 ).
4.3.7. W avelets
A mathematical tool developed originally in the ® eld of signal processing can
also be applied to fuse image data following the concept of the multiresolution analysis
(MRA) (Mallat 1989). Another application is the automatic geometric registration
of images, one of the pre-requisites to pixel based image fusion ( Djamdji et al. 1993).
The wavelet transform creates a summation of elementary functions (= wavelets)
from arbitrary functions of ® nite energy. The weights assigned to the wavelets are
the wavelet coe cients which play an important role in the determination of structure
characteristics at a certain scale in a certain location. The interpretation of structures
838 C. Pohl and J. L . van Genderen
with SAR. The results of the triple sensor image fusion combine VIR/ SAR and high
spatial resolution in one image ( Pohl 1996).
4.4.4. Mosaic / other techniques
In order to solve the cloud cover problem e ectively, the mosaic approach o ers
a wide variety of possibilities in connection with other fusion techniques. SAR data
is introduced to areas of no information on the optical data, i.e., clouds and their
shadows. Likewise, SAR data from di erent sensors or orbits can reduce the regions
of foreshortening, layover and shadow. The idea to input optical data often fails
because the mountainous areas which are causing these geometric distortions in the
SAR are also the reason for cloud coverage. Therefore, the contribution of other
SAR imagery represents the more operational solution. Once the cloud / shadow
Downloaded by [North Carolina State University] at 01:55 14 January 2013
( VIR) and layover/ foreshortening / shadow mask (SAR) has been produced, it is
possible to introduce all types of data in the various elements of the mask ( Pohl 1996 ).
classi® cation results ( Brisco et al. 1983). Multitemporal SAR is a valuable data
source in countries with frequent cloud cover and successfully used in crop
monitoring. Especially, for Developing Countries the fusion of SAR data with VIR
is a cost e ective approach which enables continuous monitoring ( Nezry et al. 1993,
Mangolini and Arino 1996 a, 1996 b). Optical and microwave image fusion is also
well known for the purpose of identifying and mapping forest cover and other types.
The combined optical and microwave data provide a unique combination that allows
more accurate identi® cation, as compared to the results obtained with the individual
sensors ( Leckie 1990, Lozano-Garcia and Ho er 1993, Kachwalha 1993, Hinse
and Coulombe 1994, Aschbacher et al. 1994, Wilkinson et al. 1995). With the
implementation of fusion techniques using multisensor optical data the accuracy of
urban area classi® cation is improved mainly due to the integration of multispectral
Downloaded by [North Carolina State University] at 01:55 14 January 2013
with high spatial resolution (Gri ths 1988, Haack and Slonecker 1991, Ranchin
et al. 1996 ).
5.5. Geology
Multisensor image fusion is well implemented in the ® eld of geology and a widely
applied technique for geological mapping. It is a well known fact that the use of
multisensor data improves the interpretation capabilities of the images. Geological
Multisensor image f usion in remote sensing 841
features which are not visible in the single data alone are detected from integrated
imagery. In most cases VIR is combined with SAR based on the fact that the data
sets complement each other. They introduce information on soil geochemistry,
vegetation and land use ( VIR) as well as soil moisture, topography and surface
roughness (SAR) ( Daily et al . 1979, Blom and Daily 1982, Rebillard and Nguyen
1982, Reimchen 1982, Haydn et al . 1982, Aarnisalo 1984, Evans 1988, Baker and
Henderson 1988, Hopkins et al . 1988, Paradella et al . 1988 and Taranik 1988, Guo
and Pinliang 1989, Harris and Murray 1989, Koopmans et al . 1989, Harris et al .
1990, Grasso 1993, Jutz and Chorowicz 1993, Koopmans and Forero 1993,
Koopmans and Richetti 1993, Ye sou et al . 1993 a, 1993 b, 1994, Ray et al . 1995).
6.1.1. Pre-processing
All sensor-speci® c corrections and enhancements of image data have to be applied
prior to image fusion since the techniques refer to sensor-speci® c e ects. After image
fusion the contribution of each sensor cannot be distinguished or quanti® ed in order
to be treated accordingly. A general rule of thumb is to ® rst produce the best single-
sensor geometry and radiometry (geocoding, ® lter, line and edge detection, etc.) and
then fuse the images. Any spatial enhancement performed prior to image fusion will
bene® t the resulting fused image. An advantage is the possibility of ® ltering and
enhancing the data during the geocoding process to avoid multiple resampling. The
data has to be resampled at the pixel spacing required for the desired image fusion
( Pohl 1996 ).
The importance of geometric accuracy to avoid artefacts and misinterpretation
in pixel based image fusion should not be underestimated. Pixels registered to each
other should refer to the same object on the ground. This implies that the data
should be geocoded with sub-pixel accuracy. The DEM therefore plays an important
842 C. Pohl and J. L . van Genderen
role in this process. The need for DEMs of high quality and appropriate grid spacing
is evident.
[Link]. RGB
Digital numbers in single images in¯ uence the colours in the RGB composite.
This implies that the following factors need consideration:
The technique is simple and does not require CPU time-intensive computations.
RGB overlay protects the contribution from optical imagery from being greatly
a ected by speckle from SAR ( Pohl 1996).
[Link]. Brovey
The spectral content of the VIR data is preserved while introducing the texture
from SAR. The resulting image is not quite as sharp as the one produced from
multiplication only. The water-land boundaries are well de® ned in the fused images;
it allows colours to be assigned to the water currents (e.g., tidal inlets) ( Pohl 1996).
[Link]. PCA
Radiometric pre-processing plays an important role in relation to the spectral
content of the fused image. The appearance of SAR signi® cantly in¯ uences the feature
visibility in the fused VIR/ SAR image. As a consequence, features that are detectable
on SAR data can be introduced to the VIR data by image fusion to complement
the data (e.g., soil moisture, urban area, oceanographic objects) ( Pohl 1996). The
PCA technique enables the integration of more than three types of data ( Zobrist
et al . 1979 ).
Multisensor image f usion in remote sensing 843
Principal component SAR images show potential for topographic mapping. This
is valid, in particular, for the three-dimensional impression of topography and change
detection. The possibilities have not yet been fully explored, e.g., combination of
principal components with optical data ( Pohl 1996 ).
[Link]. IHS
The IHS o ers a controlled visual presentation of the data. The informative
aspects are presented in IHS using readily identi® able and quanti® able colour
attributes that can be distinctly perceived. Numerical variations can be uniformly
represented in an easily perceived range of colours ( Harris et al . 1990 ). Other
techniques produce images which are di cult to interpret quantitatively and
qualitatively because the statistical properties have been manipulated and the original
Downloaded by [North Carolina State University] at 01:55 14 January 2013
integrity of the data is disturbed. Related to the IHS technique the hue has to be
carefully controlled since it associates meaningful colour with well de® ned charac-
teristics of the input. An advantage of the IHS versus the RGB is the ability of
integrating four instead of three channels in the fused image ( Rothery and Francis
1987 ). IHS fused imagery have the capability of allocating data from the SAR to
cloud covered areas without having to identify the clouds at an earlier stage. The
speckle is preserved from the SAR data in the fused image. It shows similarities with
the Brovey transformation in terms of spectral content of the imagery. A disadvantage
is the reduced spatial detail compared to original optical data ( Pohl 1996).
[Link]. Mosaic
The mosaic has an important position amongst the image fusion techniques as
far as cloud removal from VIR and the replacement of radiometrically distorted
SAR data is concerned. The result depends very much on the quality of the mask
designed for the mosaic. This is a critical point for the optical imagery. The
identi® cation of foreshortening, layover and shadow areas in the SAR is based on
DEM calculations and pure geometry. These products are often delivered with the
SAR image itself. It is essential to match the histograms of the various input data
to each other. It can be used in combination with any other image fusion technique
( Pohl 1996 ).
the sensor. In addition, pixel based methods deal with very large data volumes and
long computation times are involved. Frequently, the di culty of comparing pixels
from heterogeneous data is mentioned in the literature. Feature level image fusion
uses corresponding structures which makes the geometry a less critical issue. The
pixel based methods have the advantage against feature or information level fusion
of using the most original data. It avoids a loss of information which occurs during
a feature extraction process (Mangolini 1994 ).
Acknowledgments
The authors thank Dr H. Ye sou for his valuable comments and suggestions
on the ® rst draft of the manuscript. ITC’s research on image fusion has been
externally funded by the CEC (CHRX-CT93-0310 ), EOSAT, ESA, EURIMAGE,
EUROSENSE, MCS, MPA, NLR and NASDA. Their support to this research is
gratefully acknowledged by the authors.
References
A arnisalo, J ., 1984 , Image processing and integration of geophysical, MSS and other data
as a tool for numerical exploration in glaciated Precambrian terrain. Proceedings of
the International Symposium on Remote Sensing of Environment, Proceedings of T hird
T hematic Conference, Remote Sensing for Exploration Geology, Colorado Springs
(Colorado Springs: ISRS), pp. 107± 128.
A hern, F . J ., G oodenough, D . G ., G rey, A . L ., and R yerson, R . A ., 1978, Simultaneous
microwave and optical wavelength observations of agricultural targets. Canadian
Journal of Remote Sensing, 4, 127± 142.
A lbertz, J ., and T auch, R ., 1994, Mapping from space Ð Cartographic applications of satellite
image data. GeoJournal , 32, 29± 37.
A rmour, B ., E hrismann, J ., C hen, F ., and B owman, G ., 1994, An integrated package for
the processing and analysis of SAR imagery, and the fusion of radar and passive
microwave data. In Proceedings ISPRS Symposium, 6± 10 June 1994, Ottawa, Canada
Downloaded by [North Carolina State University] at 01:55 14 January 2013
C hiuderi, A ., and F ini, S ., 1993, Data fusion techniques for remotely sensed images. T ime-
varying Image Processing and Moving Object Recognition, Proceedings of 4th W orkshop,
Florence, Italy, 10± 11 June 1993 , pp. 1± 8.
C liche, G ., B onn, F ., and T eillet, P ., 1985, Integration of the SPOT Pan channel into its
multispectral mode for image sharpness enhancement. Photogrammetric Engineering
and Remote Sensing, 51, 311± 316.
C omhaire, A . L ., O zer, A ., and J aspar, A ., 1994, The contribution of ERS-1 and SPOT
data in the morphologic and bathymetric study of the north sea and the Scheldt
Estuary. Proceedings of First ERS-1 Pilot Project W orkshop, T oledo, Spain, 22± 24
June 1994 , SP-365 (Paris: European Space Agency), pp. 297± 302.
C orves, C ., 1994, Assessment of multi-temporal ERS-1 SAR Landsat TM data for mapping
the Amazon river ¯ oodplain. Proceedings of First ERS-1 Pilot Project W orkshop,
T oledo, Spain, 22± 24 June 1994 , SP-365 (Paris: European Space Agency), pp. 129± 132.
C rist, E . P ., 1984, Comparison of coincident Landsat-4 MSS and TM data over an agricultural
region. T echnical Papers, Proceedings 50th Annual Meeting ASP-ACSM Symposium,
ASPRS, Washington DC, U.S.A. (Washington: ASPRS), pp. 508± 517.
D aily, M ., 1983, Hue-saturation-intensity split-spectrum processing of SEASAT radar imagery.
Photogrammetric Engineering and Remote Sensing, 49, 349± 355.
D aily, M . I ., F arr, T ., and E lachi, C ., 1979, Geologic interpretation from composited
radar and Landsat imagery. Photogrammetric Engineering and Remote Sensing, 45,
1009± 1116.
D allemand, J . F ., L ichtenegger, J ., K aufmann, V ., P audyal, D . R ., and R eichert, A .,
1992, Combined analysis of ERS-1 SAR and visible /infrared RS data for land
cover/ land use mapping in tropical zone: a case study in Guinea. Space at the Service
of our Environment, Proceedings First ERS-1 Symposium, Cannes, France, 6± 8
November 1992, SP-359 (Paris: European Space Agency), pp. 555± 561.
D e G roof, H ., S owter, A ., and S ieber, A . J ., 1992, The TREES project: Creating an
environment for the routine classi® cation and image processing of ERS-1 SAR imagery.
International Space Year: Space Remote Sensing, Proceedings 12th Annual I.E.E.E.
International Geoscience and Remote Sensing Symposium (IGARSS ’92), Houston,
U.S.A. , edited by R. Williamson with collaboration of T. Stein (New York: I.E.E.E.),
pp. 1343± 1345.
D esnos, Y .-L ., M ayer, T ., and S ardar, A . M ., 1996, Multitemporal ERS-1 SAR images of
the Brahmaputra ¯ ood plains in Northern Bangladesh. Earth Observation Quarterly,
51, 6± 10.
D jamdji, J .-P ., B ijaoui, A ., and M anieri, R ., 1993, Geometrical registration of images: The
multiresolution approach. Photogrammetric Engineering and Remote Sensing, 59,
645± 653.
D omik, G ., L eberl, F ., and C imino, J . B ., 1986, Multi incidence angle SIR-B experiment
over Argentina: Generation of secondary image products. I.E.E.E. T ransactions on
Geoscience and Remote Sensing, 24, 492± 497.
D uguay, G ., H older, G ., H owarth, P ., and L eD rew, E ., 1987, Integrating remotely sensed
data from di erent sensors for change detection. Proceedings of the I.E.E.E.
International Geoscience and Remote Sensing Symposium (IGARSS ’87), Ann Arbor,
U.S.A., 18± 21 May 1987 ( New York: I.E.E.E.), p. 333± 000.
848 C. Pohl and J. L . van Genderen
E hlers, M ., 1987 , Integrative Auswertung von digitalen Bildern aus der Satelliten-
photogrammetrie und-fernerkundung im Rahmen von geographischen Informations-
systemen. Wissenschaftliche Arbeiten der Fachrichtung Vermessungswesen der
UniversitaÈ t Hanover, Germany, Dissertation, No. 149.
E hlers, M ., 1991, Multisensor image fusion techniques in remote sensing, ISPRS Journal of
Photogrammetry and Remote Sensing, 46, 19± 30.
E hlers, M ., 1993, Integration of GIS, remote sensing, photogrammetry and cartography: the
geoinformatics approach. GIS GEO-Informations-Systeme , 6, 18± 23.
ESA, 1995 a, ERS-1 SAR images of ¯ ood in Northern Europe (January 1995). Earth
Observation Quarterly, 47, 13.
ESA, 1995 b, Multitemporal image Maastricht ( The Netherlands). Earth Observation Quarterly,
47, 14.
E ssadiki, M ., 1987, A combination of panchromatic and multispectral SPOT images for
topographic mapping. IT C Journal , 1987-1, 59± 66.
Downloaded by [North Carolina State University] at 01:55 14 January 2013
with geophysical data. Proceedings of the 12th Canadian Symposium on Remote Sensing,
V ancouver, Canada ( Vancouver: CSRS), pp. 923± 926.
H arris, J . R ., M urray, R ., and H irose, T ., 1990, IHS transform for the integration of radar
imagery with other remotely sensed data. Photogrammetric Engineering and Remote
Sensing, 56, 1631± 1641.
H arrison, B . A ., and J upp, D . L . B ., 1990, Introduction to image processing. MicroBRIAN
Resource Manual, Part 2 (Melbourne: CSRIO Publications).
H aydn, R ., D alke, G . W ., H enkel, J ., and B are, J . C ., 1982, Application of the IHS color
transform to the processing of multisensor data and image enhancement. Remote
Sensing of Arid and Semi-Arid L ands, Proceedings of International Symposium (Cairo:
ISRS), pp. 599± 607.
H inse, M ., and C oulombe, A ., 1994, Le radar ¼ seul ou avec d’autres?. Numerimage, Bulletin
d ’Information Quadrimestriel (Quebec) , 3, 2± 4.
H inse, M ., and P roulx, J ., 1995, La spatiocarte regionale. Numerimage, Bulletin d’Information
Downloaded by [North Carolina State University] at 01:55 14 January 2013
Quadrimestriel (Quebec) , 3, p. 2.
H opkins, H . R ., N avail, H ., B erger, Z ., M erembeck, B . F ., B rovey, R . L ., and
S chriver, J . S ., 1988, Structural analysis of the Jura mountains± Rhine Graben inter-
section for petroleum exploration using SPOT stereoscopic data. SPOT 1: Image
utilization, assessment, results, Proceedings CNES Conference, Paris, France, Cepadues
Editions (Paris: CNES), pp. 803± 810.
H ussin, Y . A ., and S haker, S . R ., 1996, Optical and radar satellite image fusion techniques
and their applications in monitoring natural resources and land use change. Synthetic
Aperture Radar, Proceedings European Conference (EUSAR ’96), KoÈnigswinter,
Germany, 26± 28 March 1996 (Berlin: VDE), pp. 451± 456.
J anssen, L . L . F ., J aarsma, M . N ., and L inden, E . T . M . van der , 1990, Integrating
topographic data with remote sensing for land-cover classi® cation. Photogrammetric
Engineering and Remote Sensing, 56, 1503± 1506.
J utz, S . L ., and C horowicz, J ., 1993, Geological mapping and detection of oblique extension
structures in the Kenyan Rift Valley with a SPOT/ Landsat-TM datamerge.
International Journal of Remote Sensing, 14, 1677± 1688.
K achhwalha, T . S ., 1993, Temporal and multisensor approach in forest /vegetation mapping
and corridor identi® cation for e ective management of Rajaji Nat. Park, Uttar
Pradesh, India. International Journal of Remote Sensing, 14, 3105± 3114.
K annen, A ., O berstadler, R ., and P rotman, F ., 1995, Flood mapping in Germany.
Application Achievements of ERS-1, New V iews of the Earth, Preliminary V ersion (Paris:
European Space Agency), pp. 124± 125.
K attenborn, G ., N ezry, E ., D e G randi, G ., and S ieber, A . J ., 1993, High resolution detection
and monitoring of changes using ERS-1 time series. Space at the Service of our
Environment, Proceedings of the Second ERS-1 Symposium, 11± 14 October, Hamburg,
Germany, I ( Paris: European Space Agency), pp. 635± 642.
K aufmann, K .-H ., and B uchroithner, M . F ., 1994, Herstellung und Anwendungs-
moÈ glichkeiten von Satellitenkarten durch digitale Kombination von Landsat-TM und
KWR-1000-Daten. Zeitschrift f uÈ r Photogrammetrie und Fernerkundung , 62, 133± 137.
K eys, L . D ., S chmidt, N . J ., and P hillips, B . E ., 1990, A prototype example of sensor fusion
used for a siting analysis. T echnical Papers 1990, ACSM-ASPRS Annual Convention,
Image Processing and Remote Sensing, 4, 238± 249.
K nipp, K .-H ., 1993, Speckle reduction applied to ERS-1 SAR images. International Mapping
f rom Space, Proceedings ISPRS Working Group IV / 2 W orkshop and Conference,
Hanover, Germany, 27 September± 1 October 1993 , 15, 177± 190.
K ohl, H .-G ., N ezry, E ., and D e G roof, H ., 1994, Crop acreage estimation with ERS-1 PRI
images. Earth Observation Quarterly, 46, 6± 9.
K oopmans, B . N ., and F orero, G . R ., 1993, Airborne SAR and Landsat MSS as com-
plementary information source for geological hazard mapping. ISPRS Journal of
Photogrammetry and Remote Sensing, 48, 28± 37.
K oopmans, B . N ., and R ichetti, E ., 1993, Optimal geological data extraction from SPOT-
Radar synergism with samples from Djebel Amour (Algeria), Red Sea Hills (Sudan),
Sirte Basin ( Libya) and Magdalena Valley (Colombia). From Optics to Radar, SPOT
and ERS Applications, Conference Proceedings, 10± 13 May 1993, Paris, France (Paris:
Cepad ), pp. 263± 274.
850 C. Pohl and J. L . van Genderen
O tten, M . P . G ., and P ersie, M . van, 1995, SAR image of ¯ oods using PHARS system.
Earth Observation Quarterly, 47, 15.
P aradella, W . R ., V itorello, I ., L iu, C . C ., M attos, J . T ., M eneses, P . R ., and D utea, L . V .,
1988, Spectral and spatial attribute evaluation of SPOT data in geological mapping
of Precambrian terrains in semi-arid environment of Brazil. SPOT 1: Image utilization,
assessment, results, Proceedings CNES Conference, Paris, France, Cepadue s Editions
( Paris: CNES), pp. 851± 860.
P ellemans, A . H . J . M ., J ordans, R . W . L ., and A llewijn, R ., 1993, Merging multispectral
and panchromatic SPOT images with respect to the radiometric properties of the
sensor. Photogrammetric Engineering and Remote Sensing, 59, 81± 87.
P erlant, F ., 1992, Example of SPOT/ ERS-1 complementary. Space at the Service of our
Environment, Proceedings First ERS-1 Symposium, 4± 6 November 1992, Cannes, France
( Paris: ESA), pp. 563± 568.
P erlant, F ., S empere, J .-P ., and G uerre, L .-F ., 1994, Production of SPOT/ ERSI image
maps. Proceedings of First ERS-1 Pilot Project Workshop, T oledo, Spain, 22± 24 June
1994 , SP-365 ( Paris: European Space Agency), pp. 331± 336.
P ohl, C ., 1995, SPOT/ ERS Image maps: Topographic map updating in Indonesia. SPOT
MAGAZINE, December 1995, 17± 19.
P ohl, C ., 1996, Geometric aspects of multisensor image fusion for topographic map updating
in the humid Tropics. IT C publication No. 39 ( Enschede: ITC), ISBN 90 6164 121 7.
P ohl, C ., and G enderen, J . L . van, 1993, Geometric integration of multi-image information.
Space at the Service of our Environment, Proceedings of the Second ERS-1 Symposium,
11± 14 October 1993, Hamburg, Germany, ESA SP-361 (Paris: European Space Agency),
pp. 1255± 1260.
P ohl, C ., and G enderen, J . L . van, 1995, Image fusion of microwave and optical remote
sensing data for map updating in the Tropics. Image and Signal Processing for Remote
Sensing, Proceedings EUROPT O ’95, Paris, France, 25± 29 September 1995 , SPIE
Vol. 2579, pp. 2± 10.
P ohl, C ., W ang, Y ., and K oopmans, B . N ., 1994. The 1995 ¯ ood in the Netherlands from
space. IT C Journal , 1994-4, 414± 415.
P ohl, C ., K oopmans, B . N ., and W ang, Y ., 1995, Provision of ¯ ood information in the
Netherlands. In Application Achievements of ERS-1, New V iews of the Earth,
Preliminary V ersion (Paris: European Space Agency), pp. 122± 123.
P olidori, L ., and M angolini, M ., 1996, Potentials and limitations of multisensor data fusion.
Fusion of Earth Data, Proceedings EARSeL Conference, Cannes, France, 6± 8 February
1996 , pp. 13± 19.
P radines, D ., 1986 , Improving SPOT image size and multispectral resolution. Earth Remote
Sensing using the L andsat T hematic Mapper and SPOT Systems, Proceedings SPIE
Conference, 660, pp. 78± 102.
P rice, J . C ., 1987, Combining panchromatic and multispectral imagery from dual resolution
satellite instruments. Remote Sensing of Environment , 21, 119± 128.
R aggam, J ., A lmer, A ., and S trobl, D ., 1994, A combination of SAR and optical line scanner
imagery for stereoscopic extraction of 3-D data. ISPRS Journal of Photogrammetry
and Remote Sensing, 49, 11± 21.
852 C. Pohl and J. L . van Genderen
SAR and optical sensor data (TM and IRS) for geological investigations.
Photonirvachak, Journal of the Indian Society of Remote Sensing, 23, 77± 86.
R ebillard, P ., and N guyen, P . T ., 1982, An exploration of co-registered SIR-A, SEASAT
and Landsat images. RS of Environment, RS for Exploration Geology, Proceedings
International Symposium, Second T hematic Conference, Forth W orth, U.S.A. (Forth
Worth: ISRS), pp. 109± 118.
R eimchen, T . H . F ., 1982, Location of economic gypsum deposits by integration of multi-
spectral, multitemporal, geophysical, -chemical, -botanical and geological data. RS of
Environment, RS for Exploration Geology, Proceedings International Symposium, Second
T hematic Conference, Forth Worth, U.S.A . (Forth Worth: ISRS), pp. 119± 125.
R ichards, J . A ., 1984, Thematic mapping from multitemporal image data using the principal
component transformation. Remote Sensing of Environment , 16, 35± 46.
R ogers, R . H ., and W ood, L ., 1990, The history and status of merging multiple sensor data:
an overview. T echnical Papers 1990, ACSM-ASPRS Annual Convention, Image
Processing and Remote Sensing, 4, 352± 360.
R othery, D . A ., and F rancis, P . W ., 1987, Synergistic use of MOMS-01 and Landsat TM
data. International Journal of Remote Sensing, 8, 501± 508.
R uss, J . C ., 1995, T he image processing handbook . Second edn (London: CRC Press).
S chistad-S olberg, A . H ., J ain, A . K ., and T axt, T ., 1994, Multisource classi® cation of
remotely sensed data: fusion of Landsat TM and SAR images. I.E.E.E. T ransactions
on Geoscience and Remote Sensing, 32, 768± 778.
S harpe, B ., and K err, A ., 1991, Multichannel fusion technique for resolution enhancement
during ground processing. Proceedings of the I.E.E.E. International Geoscience and
Remote Sensing Symposium (IGARSS ’91), Espoo, Finland, 3± 6 June 1991 (New York:
I.E.E.E.), pp. 1459± 1462.
S heffield, C ., 1985, Selecting band combinations from multispectral data. Photogrammetric
Engineering and Remote Sensing, 51, 681± 687.
S hen, S . S ., 1990, Summary of types of data fusion methods utilized in workshop papers.
Multisource Data Integration in Remote Sensing, Proceedings of W orkshop, Maryland,
U.S.A., 14± 15 June 1990 , NASA Conference Publication 3099 (Greenbelt, MD: NASA ),
pp. 145± 149.
S hettigara, V . K ., 1992, A generalized component substitution technique for spatial enhance-
ment of multispectral images using a higher resolution data set. Photogrammetric
Engineering and Remote Sensing, 58, 561± 567.
S hufelt, J ., and M c K eown, D . M ., 1990, Use of information fusion to improve the detection
of man-made structures in aerial imagery. Multisource Data Integration in Remote
Sensing, Proceedings of W orkshop, Maryland, U.S.A., 14± 15 June 1990 , NASA
Conference Publication 3099 (Greenbelt, MD: NASA), pp. 94± 110.
S imard, R ., 1982, Improved spatial and altimetric information from SPOT composite imagery.
Proceedings ISPRS Conference, Forth W orth, U.S.A., 6± 10 December 1982 (Forth
Worth: International Archive of Photogrammetry and Remote Sensing ), pp. 433± 440.
S ingh, A ., 1989, Digital change detection techniques using remotely-sensed data. International
Journal of Remote Sensing, 10, 989± 1003.
S ingh, A ., and H arrison, A ., 1985, Standardized principal components. International Journal
of Remote Sensing, 6, 883± 396.
Multisensor image f usion in remote sensing 853
S mara, Y ., B elhadj- A issa, A ., and S ansal, B ., 1996, Decorrelation and IHS color transforms
applied for multisource satellite images integration. Fusion of Earth Data, Proceedings
EARSeL Conference, Cannes, France, 6± 8 February 1996 ( Paris: European Space
Agency), pp. 107± 112.
S trobl, D ., R aggam, J ., and B uchroithner, M . F ., 1990, Terrain correction geocoding of a
multi-sensor image data set. Proceedings 10th EARSeL Symposium, T oulouse, France
( Paris: European Space Agency), pp. 98± 107.
S uits, G ., M alila, W ., and W eller, T ., 1988, Procedures for using signals from one sensor
as substitutes for signals of another. Remote Sensing of Environment , 25, 395± 408.
T anaka, S ., S ugimura, T ., and H igashi, M ., 1989, High resolution satellite image map from
the SPOT and Landsat TM data. Advanced Space Research, 9, 115± 120.
T aranik, J . V ., 1988, First results of international investigations of the applciations of SPOT-1
data to geologic problems, mineral and energy exploration. SPOT 1: Image utilization,
assessment, results, Proceedings CNES Conference, Paris, France, Cepadue s Editions
Downloaded by [North Carolina State University] at 01:55 14 January 2013
Y e’sou, H ., B esnus, Y ., and R olet, J ., 1993 a, Extraction of spectral information from Landsat
TM data and merger with SPOT panchromatic imageryÐ a contribution to the study
of geological structures. ISPRS Journal of Photogrammetry and Remote Sensing, 48,
23± 36.
Y e’sou, H ., B esnus, Y ., R olet, J ., and P ion, J . C ., 1993 b, Merging Seasat and SPOT imagery
for the study of geologic structures in a temperate agricultural region. Remote Sensing
of Environment , 43, 265± 280.
Y e’sou, H ., B esnus, Y ., and R olet, J ., 1994, Perception of a geological body using multiple
source remotely-sensed dataÐ relative in¯ uence of the spectral content and the spatial
resolution. International Journal of Remote Sensing, 15, 2495± 2510.
Y ildimi, H ., A lparstan, E ., and O zel, M . E ., 1996, Data fusion of sharpened di erence
image by IHS transform. Fusion of Earth Data, Proceedings EARSeL Conference,
Cannes, France, 6± 8 February 1996 ( Paris: European Space Agency), pp. 125± 130.
Z obrist, A . L ., B lackwell, R . J ., and S tromberg, W . D ., 1979, Integration of Landsat,
Downloaded by [North Carolina State University] at 01:55 14 January 2013