NASA Develops Early Warning System for Detecting Forest Disturbances

“Soon after ForWarn moved into production, it detected previously unnoticed hail damage that posed a threat to a watershed. We would not have been able to do this work as efficiently without MATLAB.”

Challenge

Develop a system that uses satellite imagery to quickly detect forest disturbance threats from insects, drought, storms, blights, wildfires, and other events

Solution

Use MATLAB to process multispectral satellite images, construct multidimensional time-series data baselines, and analyze terabytes of data to help detect regionally evident forest disturbances

Results

  • New ideas implemented and tested in hours
  • Years of development time saved
  • Reusable production software delivered to growing user community
U.S. Forest Change Assessment Viewer map showing damage to the Asheville, North Carolina watershed following a 2012 hail storm.

U.S. Forest Change Assessment Viewer map showing damage to the Asheville, North Carolina watershed following a 2012 hail storm. Image courtesy ForWarn.

Working with the U.S. Department of Agriculture Forest Service and the U.S. Geological Survey, NASA’s Stennis Space Center developed ForWarn, an environmental monitoring and assessment tool that detects and tracks changes in forest vegetation nationwide. ForWarn software analyzes multispectral satellite imagery collected by MODerate resolution Imaging Spectroradiometer (MODIS) instruments on NASA’s Terra and Aqua satellites. Every 8 days, the software identifies in near real time potential forest disturbances caused by insects, drought, storms, blights, wildfires, and other events across the nation. The Forest Change Assessment Viewer for ForWarn enables users to view these changes on a map.

NASA engineers used MATLAB® to help develop and deploy ForWarn while working with other government agency partners.

“With MATLAB we produced prototypes and built products much faster than would have been possible using a traditional programming language, such as C,” says Duane Armstrong, chief of the Advanced Technology & Technology Transfer branch at the NASA Stennis Space Center. “The quick turnaround enabled by MATLAB let us demonstrate new capabilities to our customers at the Forest Service very quickly.”

U.S. Forest Change Assessment Viewer map showing vegetation changes in Gatlinburg, Tennessee, following a 2016 fire.

U.S. Forest Change Assessment Viewer map showing vegetation changes in Gatlinburg, Tennessee, following a 2016 fire (fire perimeter shown in dark red). Image courtesy ForWarn.

Challenge

To develop a near-real-time system for detecting and monitoring vegetation changes across the 48 contiguous states, NASA engineers needed to analyze multispectral images captured by the Aqua and Terra satellites, build multi-timescale baselines of the gathered data, and then compare incoming data with the baseline to identify anomalies.

Each day, MODIS data is generated for 14 gridded tiles that span the 48 states, with each tile comprising 4800 x 4800 pixels. Twice a year the system processes 15 to 25 terabytes of data when updating the yearly baselines used to identify changes in vegetation.

NASA engineers needed to process and analyze large time-series data sets of images stored as mutidimensional arrays. They needed a development environment for interactively developing and rapidly evaluating new analysis algorithms. After demonstrating their algorithms, they wanted to deploy them into production without recoding.

Solution

NASA engineers used MATLAB to develop the two key components of ForWarn: the Time Series Product Tool (TSPT), which temporally processes MODIS data, and the Phenological Parameters Estimation Tool (PPET), which uses the processed data to calculate the greenness magnitude and the day of year for the peak of the growing season and for other phenological parameters of interest to the U.S. Forest Service.

For TSPT, the team at NASA wrote a MATLAB script that retrieves MODIS data stored in hierarchical data format (HDF) files in the NASA archive and imports the data into MATLAB. TSPT invokes Mapping Toolbox™ functions to convert imported latitude and longitude data to a projected map coordinate system.

The TSPT Band Processing module, also developed in MATLAB, generates the normalized difference vegetation index (NDVI) from the time-series data, as well as soil, moisture, water, and other indices.

Working in MATLAB, the team developed algorithms for TSPT to eliminate pixels distorted by clouds, shadows, the view zenith angle, and other effects.

After merging data from the Aqua and Terra satellites into a single time series, TSPT uses Signal Processing Toolbox™ functions to identify and remove spikes and other outlying data points.

Once the TSPT algorithms have removed the outliers they further filter and resample the time series using Optimization Toolbox™ and Image Processing Toolbox™. The TSPT algorithms apply a Savitzky-Golay filter from Signal Processing Toolbox to interpolate values for any missing pixels.

The engineers used MATLAB, Optimization Toolbox, and Image Processing Toolbox to develop PPET, which performs curve fitting on the time-series data to identify vegetative states related to annual cyclical growing seasons, such as green-up, maturity, senescence, and dormancy. They later enhanced PPET to detect forest disturbances by identifying differences between daily satellite data and time-series baselines.

Using MATLAB Compiler™, the team created standalone executable versions of their MATLAB based algorithms, which can be run by users who do not have MATLAB installed.

The ForWarn team has won multiple awards for their efforts, including the Interagency Partnership Award, which recognizes federal employees from at least two different agencies who have “collaboratively accomplished outstanding work in transferring a technology.” The ForWarn team consists of the U.S. Forest Service, NASA, DOE Oak Ridge National Lab, U.S. Geological Survey, and University of North Carolina at Asheville.

Results

  • New ideas implemented and tested in hours. “Signal and image processing algorithms are built into MATLAB, so I don’t have to write them from scratch,” says Jerry Gasser, a computer engineer from Lockheed Martin. “Instead of writing and debugging several hundred lines of C code to test a new algorithm idea, I use MATLAB commands to verify it interactively and then write a script to test multiple parameter values. I get results in minutes or hours versus days or weeks.”
  • Years of development time saved. “The way MATLAB works with multidimensional data is perfect for analyzing satellite imagery,” says Armstrong. “If we had used traditional programming languages, such as C, it would have taken us much longer.”
  • Reusable production software delivered to growing user community. “ForWarn is the first near-real-time forest threat early warning system in the continental U.S.,” notes Armstrong. “Its growing user base includes over 7000 forestry experts, students, and land resource managers. We are also reusing the algorithms and modules developed in MATLAB for other customers and purposes, working with partners around the world.”

Acknowledgements