ASPRS

PE&RS July 2001

VOLUME 67, NUMBER 7
PHOTOGRAMMETRIC ENGINEERING & REMOTE SENSING
JOURNAL OF THE AMERICAN SOCIETY FOR PHOTOGRAMMETRY AND REMOTE SENSING

Highlight Article

Maximizing Information Content of Landsat Imagery for Coastal Zone Applications

Introduction
Landsat Thematic Mapper (TM) and Landsat 7 Enhanced Thematic Mapper Plus (ETM+) sensors provide multiple spectral bands containing information on reflected energy ranging from the visible through the near to mid-infrared wavelengths. Landsat imagery has been well documented for its utility in terrestrial applications such as vegetation studies and land use/land cover mapping; applications which rely primarily on the infrared bands (bands 4, 5 and 7) to discriminate between cover types. Landsat imagery has not been used as extensively for coastal zone applications, primarily due to its spatial resolution (30 meters), but also due to the difficulty of imaging aquatic environments.

For aquatic applications, the blue, green, and red wavelengths provide most of the information, because longer wavelengths (near-infrared and beyond) are absorbed by water. The transmittance and reflectance of shorter wavelengths can be affected by many factors. The primary factor is the atmosphere itself, where transmittance and reflectance are decreased due to Rayleigh and Mie scattering (Schowengerdt, 1997). In aquatic environments turbidity, depth, sun angle, and surface winds also affect the transfer of energy back to the sensor.

Many of these effects can be minimized by acquiring imagery during certain times of the year. In the coastal zone of Florida, conditions are best in the winter when biological productivity is low, bio-organic turbidity is minimal, and the climate is drier. If timed just right, when there are clear skies and water, Landsat imagery may display very little atmospheric or environmental degradation, and the data contained in the blue, green, and red bands (bands 1, 2, and 3 respectively) may be of high quality.

Creating the Cover Image
At the Florida Fish & Wildlife Conservation Commission’s Florida Marine Research Institute (FMRI), scientists were involved in a comprehensive study of Charlotte Harbor, Florida. Charlotte Harbor is located along the Gulf Coast of Florida south of Tampa (Figure 1). The Peace and Myakka Rivers flow into the Harbor, creating a dynamic estuary where salinities fluctuate with river flows and tidal forces. The shoreline is fringed with mangroves and marsh along natural stretches, and seawalls or rip-rap in areas of urban development. The Harbor is fairly shallow ( 0.25 - 6 meters), but there are a series of deeper man-made and natural channels near the river mouths and passes (7 - 20 meters).


Figure 1.
Charlotte Harbor is located along the Gulf Coast of Florida south of Tampa.
The Peace and Myakka Rivers both drain into the Harbor from the North.
Inlets and passes shelter the Harbor from the Gulf of Mexico.

A number of data sets were available to the interdisciplinary science team, including data from the Marine Resources Geographic Information System (MRGIS): a collection of spatial data maintained at FMRI. To help the team visualize the study area, two Landsat TM scenes were extracted from the MRGIS; both images were acquired on January 26, 1997. The images were cloud free, and had no apparent haze or smoke. The two scenes were rectified to sub-pixel accuracy using ground control points obtained from digital ortho-photo quadrangles (DOQs). They were then joined to form an image mosaic (base image) covering the entire Charlotte Harbor area. Because both scenes were acquired during the same satellite pass, no color balancing was necessary.

Typically a Landsat image is portrayed as a color infrared composite, displaying bands 4, 3, and 2 with red, green, and blue (RGB) respectively. This is not an optimal band combination for aquatic environments. Rather, a natural color composite using bands 3, 2, and 1 (RGB) or a false color composite using bands 4, 2, and 1 provide more detail. To increase contrast, land areas can be removed from the image (masked), allowing the full dynamic range to be applied to the features of interest.

Although the Landsat TM mosaic contained some very good information after the land mask was applied, it was still difficult to discern differences in water depth. Although deeper channels within the estuary were identifiable by their darker tones and sinusoidal patterns, submerged aquatic vegetation also appeared dark, and could be confused with deeper waters. Narrow channels were difficult to distinguish due to their smaller size, and the 30 meter pixel size of Landsat TM data. Slopes, drainage patterns, and shallows were also difficult to perceive.

To help with visualizing the depth of the Harbor, bathymetry data were obtained primarily from the National Oceanic and Atmospheric Administration National Ocean Service’s (NOAA/NOS) National Hydrographic Survey data (National Geophysical Data Center, 1998). Additional bathymetry data were compiled as necessary from NOAA nautical charts and other databases available from FMRI’s MRGIS. The shoreline and islands were included defining the areas of zero depth. The combined data, consisting of points, lines and polygons, were converted to a triangulated irregular network (TIN). A TIN represents a surface using contiguous non-overlapping triangular facets. The resolution of a TIN can vary internally in relation to the surface complexity and data availability. The TIN can be used as a master source for the creation of products such as grids or contour maps. For this project, a 10 meter bathymetry grid was produced using a linear interpolation method.

The bathymetry grid was co-registered to the base image using shoreline and island features. Displaying the bathymetry data alongside the Landsat TM mosaic was helpful, but truly useful information could be obtained by integrating them into a single image. To add the bathymetry information to the Landsat TM mosaic, hill-shading was applied to the bathymetry grid using a solar azimuth of 225°, a solar elevation of 45°, and a scaling factor of 75x. The solar angles and scaling factor are user defined, and the values were chosen to yield the most visually pleasing view of the Harbor. The Landsat TM mosaic was combined with the hill-shaded view using an overlay function, specifying a natural color band combination of 3, 2, 1 (RGB). The cover image was created by displaying the base image in the background as a false color composite, and overlaying the hill-shaded view in natural color.

The enhanced Landsat image provided a wealth of information not readily apparent in the Landsat imagery or in the bathymetry data alone. By combining bathymetry data and the spectral information in the visible TM bands, a synergistic picture of Charlotte Harbor emerged: a picture which was truly worth a thousand words!

References
National Geophysical Data Center. 1998. GEOphysical Data System for Hydrographic Survey Data (GEODAS). Volume 1, version 4.0.
Schowengerdt, R. A. 1997. Remote Sensing: Models and Methods for Image Processing. Academic Press, Chestnut Hill, MA, 522 p.

Acknowledgments
This work was supported in part under funding from the Department of the Interior, U.S. Fish and Wildlife Service, Federal Aid for Sport Fish Restoration Grant Number F-66. The author wishes to thank Bill Sargent (FMRI) for compiling the bathymetry data used in this project.

Author
Paula Houhoulis Smit
Currently the Landsat 7 Science Data Specialist with Emergent Information Technologies, Incorporated at the U.S. Geological Survey’s (USGS) EROS Data Center; working under indirect contract to NASA.

Formerly with
Florida Fish and Wildlife Conservation Commission, Florida Marine Research Institute, 100 8th Ave. S.E., St. Petersburg, FL 33701-5095, smit@usgs.gov.
Top Home