170 likes | 344 Views
DATA FUSION TO PRODUCE CLOUD-FREE TEMPORAL NDVI CROSS PLATFORM COMPOSITES USING TEMPORAL MAP ALGEBRA And SALDAS PROCESSING. Bijay Shrestha Charles O’ Hara Preeti Mali GeoResources Institute Mississippi State University. Undesired cloud cover in satellite images
E N D
DATA FUSION TO PRODUCE CLOUD-FREE TEMPORAL NDVI CROSS PLATFORM COMPOSITES USING TEMPORAL MAP ALGEBRAAnd SALDAS PROCESSING Bijay Shrestha Charles O’ Hara Preeti MaliGeoResources Institute Mississippi State University
Undesired cloud cover in satellite images Need for cloud-free satellite observation Making use of high temporal resolution of satellite data using cross-platform fusion Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua and MODIS Terra Rule based cross-platform fusion of MODIS Aqua and MODIS Terra Using Parallel Temporal Map Algebra Motivation
Vegetation Indices • Dimensionless, radiometric measures that function as indicators green vegetation. • Normalized Difference Vegetation Index (NDVI) • Healthy, chlorophyll-based vegetation strongly reflects near-infrared wavelengths and reflects relatively weakly in the visible red. • Range [-1, +1] • NDVI is highly useful in vegetation studies.
MODIS Datasets • 6 MODIS datasets used • MYD09GQK - 250m surface reflectance MODO9GQK data • MYDMGGAD - 1 Km Geo-location angles MODMGGAD data • MYD09GST - 1 Km surface reflectance MOD09GST quality data • MYD - Aqua • MOD - Terra
Pre-processing Surface reflectance day 1 NDVI day 1 Aqua and Terra NDVI cubes Surface reflectance day 2 NDVI day 2 … … Surface reflectance day N NDVI day N QMask day 1 Surf.Refl. Quality day 1 Input to Fusion Algorithm Surf. Refl. Quality day 2 QMask day 2 Aqua and TerraQuality Mask cube … … Surf. Refl. Quality day N QMask day N Geolocation Angles day 1 view zenith angle day 1 Geolocation Angles day 1 view zenith angle day 2 Aqua and Terra View zenith angle cube … … Geolocation Angles day N view zenith angle day N
TPV5 TPV5 TPV4 TPV4 Aqua MODIS Daily Data Products MOD09, Quality, and Scan Angle TPV2 TPV2 TPV3 TPV3 Cross-Platform Fused Daily Data Product TPV1 TPV1 DAILY DAILY FUSION Terra MODIS Daily Data Products MOD09, Quality, and Scan Angle Daily Metadata Product Fusion
Aug 20, 2005 Aug 21 2005 Aug 22 2005 Aug 23 2005 Aug 24 2005 Aug 25 2005 Aug 26 2005 Aug 27 2005 Aug 28 2005 Aug 29 2005 Aug 30, 2005 Aug 31, 2005 Sep 1, 2005 Sep 2, 2005 Sep 3, 2005 Sep 4, 2005 Sep 5, 2005 Sep 6, 2005 Sep 7, 2005 Sep 8, 2005 Sep 9, 2005 Sep 10, 2005 Sep 11, 2005 Sep 12, 2005 Sep 13, 2005 Sep 14, 2005 Sep 15, 2005 Sep 16, 2005 Sep 17, 2005 Sep 18, 2005 Sep 19, 2005 Sep 20, 2005 Results AQUA MODIS, Mississippi Delta during Hurricane Katrina
Aug 20, 2005 Aug 21 2005 Aug 22 2005 Aug 23 2005 Aug 24 2005 Aug 25 2005 Aug 26 2005 Aug 27 2005 Aug 28 2005 Aug 29 2005 Aug 30, 2005 Aug 31, 2005 Sep 1, 2005 Sep 2, 2005 Sep 3, 2005 Sep 4, 2005 Sep 5, 2005 Sep 6, 2005 Sep 7, 2005 Sep 8, 2005 Sep 9, 2005 Sep 10, 2005 Sep 11, 2005 Sep 12, 2005 Sep 13, 2005 Sep 14, 2005 Sep 15, 2005 Sep 16, 2005 Sep 17, 2005 Sep 18, 2005 Sep 19, 2005 Sep 20, 2005 Results (Contd.) TERRA MODIS, Mississippi Delta during Hurricane Katrina
Aug 20, 2005 Aug 21 2005 Aug 22 2005 Aug 23 2005 Aug 24 2005 Aug 25 2005 Aug 26 2005 Aug 27 2005 Aug 28 2005 Aug 29 2005 Aug 30, 2005 Aug 31, 2005 Sep 1, 2005 Sep 2, 2005 Sep 3, 2005 Sep 4, 2005 Sep 5, 2005 Sep 6, 2005 Sep 7, 2005 Sep 8, 2005 Sep 9, 2005 Sep 10, 2005 Sep 11, 2005 Sep 12, 2005 Sep 13, 2005 Sep 14, 2005 Sep 15, 2005 Sep 16, 2005 Sep 17, 2005 Sep 18, 2005 Sep 19, 2005 Sep 20, 2005 Results (Contd.) Cross Platform Fused NDVI
Results show 5 day look-back for Aqua Terra fusion Pixels that fulfill the quality criteria Temporal compositing based on sliding window can be performed. Temporal Weighting Quality Coefficient and Temporal Weighting Results (Contd.)
Day of Interest 1 1 1 1 1 1 3 1 3 1 1 1 3 1 3 3 3 Day n - 1 1 3 1 3 1 3 1 1 1 Day n 1 1 1 Days Day n + 1 Temporal Compositing NDVI comp = (NDVI-1 * 1+ NDVI0 * 3 + NDVI1 *1)/ (1+3+1) NDVI comp = (NDVI-1 * FQCC-1*1+ NDVI0 * FQCC0* 3 + NDVI1 *FQCC1*1)/ (FQCC-1*1 + FQCC0*3+FQCC1*1)
Zonal Processing • 9 random Zonal areas were identified in the AOI
Zonal Processing (Contd.) • Calculated Mean Pixel values for each zone for each day in the temporal period of interest Katrina
Area Of Interest • Argentina, South America • Time Period of Interest • August 2004 – April 2006
South American Land Data Assimilation Systems SALDAS Radiation Data 1/8 of a degree (about 15 km) NetCDF/Binary Files Data Sets Precipitation Solar Radiation Minimum Temperature Maximum Temperture Soil Moisture – 4 levels Temporal Range SALDAS Data Streams Eg. Solar Radiation DatSet
SALDAS Data Streams • Current processing • To ASCII using ENVI scripts • To Grids using AML scripts • Zonal delineation of study fields • Zonal Mean Processing of study fields