250 likes | 391 Views
Destriping VIIRS brightness temperatures for SST. Karlis Mikelsons Marouan Bouali Alex Ignatov Yury Kihai. STAR JPSS Annual Meeting College Park, MD May 14, 2014. Striping in VIIRS data. Low amplitude Unidirectional artifact Strongly affects gradients. Destriping method.
E N D
Destriping VIIRS brightnesstemperatures for SST KarlisMikelsons MarouanBouali Alex Ignatov YuryKihai STAR JPSS Annual Meeting College Park, MD May 14, 2014
Striping in VIIRS data • Low amplitude • Unidirectional artifact • Strongly affects gradients Destriping of brightness temperatures...
Destriping method • Start with striped image • Calculate gradients • Discard “y” gradients in striped, but otherwise smooth regions • Poisson reconstruction (with DCT using FFT) yields approximate destriped image • Split the original image into destripedand striped components original image destriped component residual striped component Algorithm: M. Bouali, A. Ignatov, J. Atmos. Oceanic Technol., 31, 150-163 (2014). Destriping of brightness temperatures...
Destriping method: iterative refinement • At each iteration, contribution to destriped image is extracted from residual striped component • Repeat until destriped component contains (nearly) all useful information and residual is (nearly) reduced to stripes residual striped component destriped component Algorithm: M. Bouali, A. Ignatov, J. Atmos. Oceanic Technol., 31, 150-163 (2014). Destriping of brightness temperatures...
Destriping method: nonlinear filter residual (striped) filtered residual filtering domain nonlinear filter last iteration destriped component final destriped image Destriping of brightness temperatures...
Results – VIIRS band M12 (3.7mm) Destriping of brightness temperatures...
Results – VIIRS band M12 (3.7mm) Destriping of brightness temperatures...
Results – VIIRS band M15 (10.75mm) Destriping of brightness temperatures...
Results – VIIRS band M15 (10.75mm) Destriping of brightness temperatures...
Results – VIIRS band M16 (12.01mm) Destriping of brightness temperatures...
Results – VIIRS band M16 (12.01mm) Destriping of brightness temperatures...
Results – VIIRS band M12 (3.7mm) – day (glint) Striping in glint region primarily due to different viewing angle for detectors Study: Q. Liu, C. Cao, F. Weng, J. Atmos. Oceanic Technol., 30, 2478-2487 (2013). Destriping of brightness temperatures...
Results – VIIRS band M12 (3.7mm) – day (glint) • Areas outside glint region and onset of glint region are destriped • High amplitude striping in center of glint region are not destriped Destriping of brightness temperatures...
Results – VIIRS band M15 (10.76mm) – day longer wavelength bands unaffected by glint Destriping of brightness temperatures...
Results – VIIRS band M15 (10.76mm) – day longer wavelength bands unaffected by glint Destriping of brightness temperatures...
Results – VIIRS band M12 (3.7mm) – effect on cloud mask Cloud cover identification affected by striping Destriping of brightness temperatures...
Results – VIIRS band M12 (3.7mm) – effect on cloud mask • Destriping avoids striping in cloud cover • Also helps to identify larger areas as cloud free Destriping of brightness temperatures...
Results – MODIS Aqua band 20 (3.75mm) Destriping of brightness temperatures...
Results – MODIS Aqua band 20 (3.75mm) Destriping of brightness temperatures...
Results – MODIS Aqua band 31 (11.0mm) Destriping of brightness temperatures...
Results – MODIS Aqua band 31 (11.0mm) Destriping of brightness temperatures...
Results – MODIS Aqua band 32 (12.0mm) Destriping of brightness temperatures...
Results – MODIS Aqua band 32 (12.0mm) Destriping of brightness temperatures...
Performance – IDL vs C • C code is about 6 times faster • I/O is a significant factor for C version: ≈25% time (VIIRS) and ≈40% time (MODIS) Destriping of brightness temperatures...
Summary • We have developed a fast, production ready destriping code • Newly developed C version is 6 times faster than IDL • Results show a significant visual improvement for brightness temperatures • Destriping helps to determine the cloud cover • Future • Immediate: • Incorporate destriping into ACSPO (experimental mode initially) • Destripe additional IR bands (M13, M14 VIIRS; 22, 23, 29 MODIS) • Near term: • Incorporate destriping into NOAA operations • Destripe visible bands for cloud cover (M6, M7 VIIRS; 6, 7 MODIS) • Address saw-like modulations in glint areas (short wavelength bands) • Optimize codes for reprocessing of historic data • VIIRS • MODIS (Terra + Aqua) Destriping of brightness temperatures...