Unit 11
Unit 11
11.1 INTRODUCTION
Image filtering and band ratioing are image enhancement processes. In the previous unit, you have
learnt that image enhancement is a technique used to improve the readability and capacity of
remotely sensed images for the extraction of meaningful data by humans.One of the image
enhancement techniques, known as contrast enhancement (also called radiometric enhancement)
is used to increase contrast and brightness of an image for extracting useful information (mostly) by
human interpreters. This technique is discussed in detail in unit 10.
Image filtering and band ratioing are other important techniques of image processing. A satellite
image or picture may be enhanced or modified through the use of filtering. Filtering includes the
image processing steps of edge enhancement, sharpening, and smoothing.
The digital image filtering makes images suitable for land use studies and for extracting many
geological structures.
Block 2 Image Pre-processing Techniques
……………………………….……….............…………………………….……………………….................……
The band ratioing technique is used to enhance the spectral difference
between electromagnetic bands by dividing the measure of reflectance for a
pixel in one image band by its geometrically corresponding pixel from another
image band. It is largely used to minimise illumination impacts in
landscape/terrain variation or to highlight objects of interest. In this unit, we will
discuss digital image filtering approaches, basics of band ratioing and
significance of commonly used spectral indices.
This filter's response is just the average of the pixels in the vicinity of the filter
mask. Consider the following picture patch, which is depicted in Fig. 11.2a. If
the aforementioned kernel is applied to the image's central 3x3 block, we obtain
the low pass filtered output, which is depicted in Fig. 11.2b. 5x5, 7x7 and other
sizes are also possible for the kernel.
10 20 30 40 50 10 20 30 40 50
20 20 60 80 70 20 31 47 52 70
30 60 30 90 20 30 47 61 55 20
40 70 100 40 10 40 62 64 47 10
50 100 80 10 50 50 100 80 10 50
(a) (b)
Fig. 11.2: The logic of applying an average linear filter in an image: a) Input image
before filter; and b) Output image after using an average linear filter.
(a)
(b) (c)
Fig. 11.3: Images showing :a) Original image; b) Noise added input image; and c)
Averaging filtered output image.(Source:https://bhuvan-
app1.nrsc.gov.in/imagegallery/bhuvan.html)
10 20 30 40 50
20 20 60 80 70
30 60 30 90 20
40 70 100 40 10
50 100 80 10 50
10 20 30 40 50 10 20 30 40 50
20 20 60 80 70 20 20 60 80 70
30 255 30 90 20 30 40 30 90 20
40 70 100 40 10 40 70 100 40 10
50 100 80 10 50 50 100 80 10 50
a) B
(a) (b)
Fig. 11.5: Filtered image obtained after applying a 3x3 median filter on image
shown in Fig. 11.4: a) Neighbourhood of shaded pixel and b) Median is
replaced as 40.
(a) (b)
(c)
Fig. 11.6: Filtering of an image: a) Laplacian filter kernel; b) Input image; and c)
sharpening filter output image (Source: https://bhuvan-app1.nrsc.gov.in/
image gallery/bhuvan.html for a)
Fig. 11.7: Schematic showing filtering in frequency domain. (Source: modified after
Gonzales and Woods, 2018)
(d) (e)
Fig. 11.8: Schematics showing: a) Perspective view of frequency domain low-
pass filter; b) 2-D image; c) Low pass filter transfer function radial
cross-section; d) Noise added input image; and e) Low pass filtered
output image in frequency domain. (Source: Gonzales and Woods, 2018)
(a) (b)
Fig. 11.9: Satellite images showing: a) Noise-added input image; and b) Output of
frequency domain high pass filter . (Source:https://bhuvan-
app1.nrsc.gov.in/imagegallery/bhuvan.html)
SAQ I
a) Describe various types of image filters.
b) What will happen if you make the low pass kernel larger, from 3x3 to 7x7?
(a) (b)
(c) (d)
Fig. 11.10: True colour and band ratio images. a) Landsat TM true colour
composite (RGB); b) Thematic mapper band (TM) 3/TM4: This ratio
(d)
has defined barren lands and urban area uniquely but it could not
define water body, forests and croplands ;c) TM2/TM3: this ratio has
distinguished croplands, barren lands sharply but it has not
separated croplands, forests and water body. Both forests and water
body have appeared as lighter tone and barren land appeared as dark
tone and it did not enhance urban area. Since chlorophyll has strong
reflectance in the band 2 (0.52 -0.60 um) region and strong absorption
in the band 3 (0.63 -0.69 um) region, vegetation has appeared as
brighter tone; and d) TM3/TM5: This ratio enhanced barren lands,
highways, street patterns within the urban areas and urban built-up or
cemented areas. It could not enhance clear water but it enhanced
turbid water. This ratio is useful for observing differences in water
turbidity. Barren lands, highways, urban and built-up areas have
appeared as lighter tone and forests, water body and croplands
appeared as dark tone.
Ratio also gives unique information which is not present in any single band and
is useful for differentiating among various features such soils and vegetation.
Mathematically, the ratio function can be expressed as:
BVi,j,r = BVi,j,k/BVi,j,l
Where BVi,j,r is the value of output ratio for the pixel at row i and column j
Contributor: Prof. B. Sathya Bama 229
Block 2 Image Pre-processing Techniques
……………………………….……….............…………………………….……………………….................……
BVi,j,kand BVi,j,l represent the value of brightness at the same location in band k
and l.
If BVi,j,l is equal to zero, then it would become difficult to compute the ratio. In
this case there are some available alternatives, which used by assigning some
value to BVi,j,l between 0 and 1.
Band ratioing technique is generally applied on a multispectral image. Ratio is
an effective technique for selectively enhancing spectral features (Fig. 11.10).
Ratio images derived from different band pairs are often used to generate ratio
colour composites in a Red, Blue and Green (RGB) display. Many indices, such
as the normalised difference vegetation index (NDVI) have been developed
based on both differencing and ratio operations.
Band ratio is an effective technique for suppressing topographic shadows. For a
given incident angle of solar radiation, the radiation energy received by a land
surface depends on the angle between the land surface and incident radiation
(Fig. 11.11). Therefore, solar illumination on a land surface varies with terrain
slope and aspect, which results in topographic shadows. In a remotely sensed
image, the spectral information is often occluded by sharp variations of
topographic shadowing. DNs in different spectral bands of a multi-spectral
image are proportional to the solar radiation received by the land surface and
its spectral reflectance.
Fig. 11.11: Ratio of pixel values in NIR region to the corresponding pixel value in
the visible red region of the spectrum. The ratios for the illuminated and
shaded slopes are very similar, although pixel value differs by a factor
of more than two. Hence, an image made up of NIR:R ratio values at
pixel positions will exhibit a much reduced shadow or topographic
effect. (Source: redrawn after Liu and Mason, 2009)
Fig. 11.12: Healthy vegetation absorbs most of the visible light that hits it and
reflects a large portion of NIR light. Unhealthy or sparse vegetation
reflects more visible light and less NIR light
The most popular and used satellite instrument for collecting NDVI is US
National Oceanic and Atmospheric Administration’s Advanced Very High
Resolution Radiometer (AVHRR) satellite. It is sensitive to wavelengths from
0.55 - 0.7 µm and 0.73 - 1.0 µm, both of which are idealised in NDVI
calculation. AVHRR’s detectors measure intensity of light being reflected from
the different bands. Landsat TM is also used to calculate NDVI but because its
band wavelengths differ (uses bands 3 and 4), it is most often used to create
images with greater detail covering less area. Moderate Resolution Imaging
Spectroradiometer (MODIS) sensor of NASA also has an NDVI standard
product.
11.4.3 Soil Adjusted Vegetation Index
Similar to the NDVI, the Soil Adjusted Vegetation Index (SAVI) is employed in
regions with low (40%) vegetative cover. In order to calculate SAVI, the
following is used:
SAVI = (NIR – Red) / (NIR + Red) x (1 + L)
Where, L is a correction factor, which varies from 0 for very high vegetation
cover to 1 for very low vegetation cover. For intermediate vegetation cover, a L
value of 0.5 is commonly used; when L is equal to zero, the SAVI becomes the
same equation as the NDVI.
11.4.4 Atmospherically Resistant Vegetation Index
An enhanced index called the Atmospherically Resistant Vegetation Index
(ARVI) is used to account for the influence of the atmosphere. It works best in
places with a lot of atmospheric aerosol, like tropical areas with soot pollution.
The benefit of this index is that it uses the blue band in addition to the red and
NIR bands. A self-correction method for the atmospheric influence on the red
band gives the ARVI greater tolerance to atmospheric effects than other
indices. The difference in radiance between the blue and red bands is
performed in order to adjust the radiance in the red channel (Kaufman and
Contributor: Prof. B. Sathya Bama 233
Block 2 Image Pre-processing Techniques
……………………………….……….............…………………………….……………………….................……
Tanre, 1992).Green vegetation typically has ARVI values between 0.20 and
0.80, with a range of -1 to 1. The formula below is used to compute ARVI:
ARVI = (NIR – RB) / (NIR + RB)
RB = Red – γ(Blue band – Red band)
Where RB = red band and γ = 1
11.4.5 Normalised Difference Snow Index
Normalised Difference Snow Index (NDSI) is used to identify snow cover and is
often used for hydrological forecasting. It is calculated by using the following:
NDSI = (Visible – SWIR)/(Visible + SWIR)
Where SWIR = Shortwave infrared
Snow cover reflects strongly and appears as brilliant as clouds in the visible
spectrum. Thus, it is visually challenging to discern from cloud cover. Snow
does not reflect much light in the shortwave infrared region of the spectrum,
allowing clouds and snow cover to be distinguished from one another. While the
threshold may change based on the imagery and kind of snow or ice, higher
NDSI values (> 0.4) often indicate the presence of snow.
Fig. 11.14: NDVI of the satellite image given in Fig. 11.13 highlighting vegetation
rich areas in green.
It may be noted that when a considerable portion of the soil surface is exposed,
the soil reflectance can manipulate the NDVI values (Fig.11.15). Light reflected
from the soil can have a significant effect on NDVI values and can the values by
up to 20%. Compared to other indices, the ARVI resists the atmospheric effects
and helpful to identify vegetation at different levels (Fig. 11.16).
Fig.11.15: SAVI of a satellite image highlighting the regions of high soil content in
brown.(Source:https://bhuvan-app1.nrsc.gov.in/imagegallery/bhuvan.html)
Fig. 11.16: ARVI of the satellite image shown in Fig. 11.19 highlighting vegetation
rich areas in green.
SAQ II
a) List out the any four spectral Indices.
b) What is the use of NDVI?
11.6 ACTIVITY
Write the formulae of the following:
i) BVi,j,r =
ii) ARVI =
iii) NDSI =
11.7 SUMMARY
Image filtering and band ratioing are image enhancement processes.
Filtering could be a procedure for modifying or enhancing a satellite image
or picture.
The digital image filtering makes images suitable for land use studies and
for extracting many geological structures.
Spatial and frequency filtering are two widely used approaches to digitally
filter images.
The band ratioing technique is used to enhance the spectral difference
between electromagnetic bands by dividing the measure of reflectance for a
pixel in one image band by its geometrically corresponding pixel from
another image band.
It is largely used to minimise illumination impacts in landscape/terrain
variation or to maximise phenomenon in objects of interest.
It is an effective technique for suppressing topographic shadows.
11.9 REFERENCES
Jensen J R (2018) Digital Image Processing: A Remote Sensing
Perspective, 4th edition, Pearson India Education Services Pvt Ltd, Uttar
Pradesh, India.
Kaufman YJ and Tanre D (1992) Atmospherically Resistant Vegetation
Index (ARVI) for EOS-MODIS. IEEE Transactions on Geoscience and
Remote Sensing, 30, 261-270.
Liu JG and Mason P J (2009) Essential Image Processing and GIS for
Remote Sensing, Wiley-Blackwell, U.K.
https://support.esri.com/en-us/gis-dictionary/band-ratio
11.10 ANSWERS
SAQ I
a) Spatial and frequency domain filters.
b) If we increase the size of the low pass kernel from 3x3 to 7x7, the image get
further smoothed and you will get a blurred appearance.
SAQII
a) SR, NDVI, SAVI, ARVI
Contributor: Prof. B. Sathya Bama 237
Block 2 Image Pre-processing Techniques
……………………………….……….............…………………………….……………………….................……
b) Normalized Difference Vegetation Index (NDVI) enhances the vegetation
and more specifically the healthy vegetation.
Terminal Questions
1. Please refer to section 11.2.
2. Please refer to section 11.3.
3. \Please refer to subsections 11.4.1 and 11.4.2.