The present invention relates to generating a vegetation index image of a geographic region based upon satellite measurements and, more particularly, to methods and systems for reducing error in such images.
Photosynthesizing plants absorb visible light, especially light in the blue and red wavelength bands, but reflect and scatter light in the near-infrared (NIR) wavelength band. Thus, vegetation, such as green leaves, appears “dark” to red and blue wavelength band reflectance sensors and appears “bright” to NIR wavelength band sensors. Soil and water absorb and reflect red and NIR wavelength bands approximately equally. By contrast, clouds and snow reflect visible light and absorb NIR. Thus, clouds appear “bright” to red and blue wavelength band sensors and “dark” to NIR wavelength band sensors.
Vegetation indexes (VI's) are used to identify and study vegetation. Two common indexes are the Normalized Difference Vegetation Index (NDVI) and the Enhanced Vegetation Index (EVI). The NDVI is determined using the equation NDVI=(NIR−Red)/(NIR+Red), where Red and NIR are the spectral reflectance measurements acquired in thus, take on values between 0.0 and 1.0. Hence, the NDVI varies between −1.0 (clouds and snow) and +1.0 (dense forests). The EVI is determined by the equation EVI=[G*(NIR−Red)]/[NIR+(c*Red)−(c*Blue)+L], where G is a gain factor, c is an absorption factor, L is a canopy correction factor, and Red, Blue and NIR are the spectral reflectance measurements acquired in the red, blue and near-infrared wavelength bands, respectively. The NDVI and EVI complement each other in global vegetation studies and improve the detection of vegetation changes.
VI's are typically determined by using a ratio of reflected light wavelength bands from the planetary surface measured by satellite remote sensors. It is possible for locations on the planetary surface to be obstructed from view of the satellite remote sensors by clouds, fog, smoke, haze or poor satellite viewing angles, for example. These obstructions interfere with the detection of the wavelength bands reflected from the planetary surface, thereby reducing the accuracy of the determined VI's and images based upon the VI's.
Presently, composite VI images from 8, 16 or 30 day periods are typically used in an attempt to compensate for the effect of obstructed locations. In theory, increasing the number of days improves accuracy because it becomes less likely for the location to be obstructed for the entire period. Gathering measurements over several days to make one composite image decreases the temporal resolution of the VI images. Additionally, the effect of some obstructions may remain in these composite images (e.g., thin or shallow clouds).
The present invention is directed to methods and systems for generating vegetation index images. The vegetation index images are generated by identifying one or more obscured locations within a geographic region, determining a vegetation index value for each location within the geographic region, replacing the determined vegetation index value for each obscured location with a predetermined replacement value for that location, and generating a vegetation index image based on the predetermined replacement values for the obscured locations and the determined vegetation index values for the non-obscured locations.
In exemplary embodiments of the present invention, obscured locations may be identified by obtaining a land surface temperature measurement for each location within the geographic region, determining a climate difference between the land surface temperature measurements and a corresponding climatology land surface temperature for each location within the geographic region, and identifying each location having a determined climate difference greater than a climatology tolerance value for that location as an obscured location. In other embodiments of the present invention, obscured locations may be identified by obtaining a land surface temperature measurement for each location within the geographic region, determining a previous difference between the land surface temperature measurements and a corresponding previous land surface temperature value for each of the locations within the geographic region, and identifying each location having a determined difference greater than a predefined difference value as an obscured location. Yet other embodiments of the present invention may use both methods simultaneously.
The predetermined replacement value for each of the identified obscured locations may be NO DATA, an estimate of a present vegetation index value for each obscured location based upon at least one previously determined vegetation index for that location, or an estimate of a present vegetation index value for each obscured location based upon one or more current vegetation index values for locations surrounding that obscured location.
In another exemplary embodiment of the present invention, a computer readable medium includes software configured to control a computer to implement a method for generating a vegetation index image. The vegetation index image is generated by identifying one or more obscured locations within the geographic region, determining a vegetation index value for each location within the geographic region, replacing the determined vegetation index value for each obscured location with a predetermined replacement value for that location, and generating a vegetation index image based on the predetermined replacement values for the obscured locations and the determined vegetation index values for the non-obscured locations.
Yet another exemplary embodiment of the present invention is a system for generating a vegetation index image of a geographic region. The system includes means for identifying one or more obscured locations within the geographic region, means for determining a vegetation index value for each location within the geographic region, means for replacing the determined vegetation index value for each obscured location with a predetermined replacement value for that location, and means for generating a vegetation index image based on the predetermined replacement values for the obscured locations and the determined vegetation index values for the non-obscured locations.
The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
The invention is best understood from the following detailed description when read in connection with the accompanying drawings, with like elements having the same reference numerals. When a plurality of similar elements are present, a single reference numeral may be assigned to the plurality of similar elements with a small letter designation referring to specific elements. When referring to the elements collectively or to a non-specific one or more of the elements, the small letter designation may be dropped. The letter “n” may represent a non-specific number of elements. This emphasizes that according to common practice, the various features of the drawings are not drawn to scale. On the contrary, the dimensions of the various features are arbitrarily expanded or reduced for clarity. Included in the drawings are the following figures:
The present invention may be used to generate a vegetation index (VI) image having fewer inaccuracies due to obstructions such as clouds, fog, smoke, and haze than VI images generated using composite data compiled over as long as 30-day time periods. Accordingly, as will be described below, more accurate VI images can be generated while maintaining temporal resolution.
In step 404, the difference between a climatology LST value and the LST measurement for each of the plurality of locations is determined. The climatology LST value is an expected LST value based upon historical climatological trends for the geographic region. In an exemplary embodiment, during step 404, processor 210 calculates the difference between the obtained LST measurement stored in memory 202 and the climatology LST value stored in memory 202 for each location 114 within the geographic region 112.
In step 406, the difference between the climatology LST and the measured LST is compared to a climatology tolerance value (e.g. 5 degrees Celsius). In an exemplary embodiment, if the difference between the climatology LST and the measured LST is greater than the climatology tolerance value, it may be suspected that the location is obscured and processing proceeds to step 408. Conversely, if the difference is less than the climatology tolerance value, it may be determined that the location is not obscured and processing proceeds to step 410 with that location identified as non-obscured. In an exemplary embodiment, during step 406, processor 210 compares the difference to a climatology tolerance value stored in memory 202 and identifies as obscured each location where the difference between the measured LST and the climatology LST is greater than the climatology tolerance value.
In step 408, the difference between the measured LST and a previously measured LST is determined for each location having a climatology difference greater than the climatology tolerance value. The previously measured LST is based upon the LST value from a previous day or composite period for that location. In an exemplary embodiment, processor 210 calculates the difference between the obtained LST measurement and the previously measured LST value stored in memory 202.
In step 412, the difference between the measured LST and the previously measured LST is compared to a predefined difference value (e.g., 10 degrees Celsius). In an exemplary embodiment, if the difference between the measured LST and the previously measured LST is greater than the predefined difference value, it may be determined that the location is obscured and processing proceeds to step 414 with that location identified as obscured. Conversely, if the difference is less than the predefined difference value, it may be determined that the location is not obscured and processing proceeds to step 410 with that location identified as non-obscured. In an exemplary embodiment, processor 210 compares the difference to the predefined difference value stored in memory 202 and identifies as obscured each location where the difference is greater than the predefined difference value.
Referring back to
In step 306, the spectral reflectance measurements are used to determine VI values for each location within the geographic region. For example, if NDVI values are to be determined, the equation NDVI=(NIR−Red)/(NIR+Red) is used, where Red and NIR stand for the spectral reflectance measurements acquired in the red and near-infrared wavelength bands, respectively. Similarly, if EVI values are to be determined, the equation EVI=[G*(NIR−Red)]/[NIR+(c*Red)−(c*Blue)+L] is used, where G is a gain factor, c is an absorption factor, L is a canopy correction factor, and Red, Blue and NIR stand for the spectral reflectance measurements acquired in the red, blue and near-infrared wavelength bands, respectively. In an exemplary embodiment, during step 306 processor 210 calculates the VI for each location 114a-n within the geographic region 112.
In step 308, VI values for obscured locations are replaced with a replacement value, e.g., by processor 210. In an exemplary embodiment, the replacement value for each obscured location may be a unique value representing a lack of data (e.g., “NO DATA”). In an alternative exemplary embodiment, the replacement value for each obscured location may be an estimate based upon at least one previously determined vegetation index for that location. In yet another exemplary embodiment, the replacement value for each obscured location may be an estimate based upon one or more current vegetation index values for locations surrounding that obscured location. Other suitable methods for determining replacement values will be understood by one of skill in the art from the description herein.
In an exemplary embodiment, processor 210 may calculate each replacement value as an estimate based upon at least one previously determined vegetation index stored in memory 202 for that obscured location. Alternatively, processor 210 may calculate each replacement value as an estimate based upon one or more current vegetation index values stored in memory 202 for locations surrounding that obscured location. Processor 210 may store calculated replacement values to memory 202 for use in subsequent replacement of VI values for obscured locations.
In step 310, a VI image is generated based on the replacement values for the obscured locations and the determined vegetation index values for the non-obscured locations. In an exemplary embodiment, the VI image is generated on a display 204 (
It is contemplated that the methods previously described may be carried out within a computer instructed to perform these functions by means of a computer-readable medium. Such computer-readable media include integrated circuits, magnetic and optical storage media, as well as audio-frequency, radio frequency, and optical carrier waves.
Although the invention is illustrated and described herein with reference to specific embodiments, the invention is not intended to be limited to the details shown. Rather, various modifications may be made in the details within the scope and range of equivalents of the claims and without departing from the invention.
This application is related to and claims the benefit of U.S. Provisional Application Ser. No. 60/961,132 entitled VEGETATION INDEX IMAGE GENERATION METHODS AND SYSTEMS filed on Jul. 19, 2007, the contents of which are incorporated herein by reference.
The U.S. Government has a paid-up license in the present invention and the right in limited circumstances to require the patent owner to license others on reasonable terms as provided for by contract as awarded by the National Aeronautics and Space Administration under funding number NASA Space Grant (NNG05G092H).
Number | Date | Country | |
---|---|---|---|
60961132 | Jul 2007 | US |