150 likes | 310 Views
McIDAS achievements. Compiled & Installed McIDAS-X v2008 Tested BOWTIE correction Enhanced MOD4 and MOD6 server code to reflect changes in v005 algorithms. Added ADDE server for MODIS AEROSOL Properties (MOD4), and MODIS CLOUD Properties (MOD6).
E N D
McIDAS achievements • Compiled & Installed McIDAS-X v2008 • Tested BOWTIE correction • Enhanced MOD4 and MOD6 server code to reflect changes in v005 algorithms. • Added ADDE server for MODIS AEROSOL Properties (MOD4), and MODIS CLOUD Properties (MOD6). • Implemented and tested export to netCDF for Hyrax & THREDDS catalog server
ADDE server enhancements mcidas2008/src --------------------------------------------------------- imgbowtie.c - new XRD command gennavblk.for - needed for IMGBOWTIE makefile - MODIS ADDE server modifications modx.h - SSEC 2008 version with mods modradir.cp - SSEC 2008 version modraget.cp - SSEC 2008 version mod3aget.cp - SSEC 2008 XRD version modiadir.cp - SSEC 2008 XRD version with mods modiaget.cp - SSEC 2008 XRD version with mods modladir.cp - SSEC 2008 XRD version modlaget.cp - SSEC 2008 XRD version modxadir.cp - SSEC 2007 version with mods modxaget.cp - SSEC 2007 version with mods kbxmodx.dlm - SSEC 2007 version with mods mcidas2008/data --------------------------------------------------------- SATBAND - modifications for MODIS ADDE servers
Vxeos required enhancements • Bow-tie correction fails on MODIS Swaths with missing scan lines -> centre location cannot be resolved (lat: -999 lon:-999) Task: check with Dominic if this can be fixed. • Tertiary (L3B) MODIS products (NDVI,LST,CRL) are currently produced in non-IMAPP compatible HDF structure (see /eOsphere/bin for relevant scripts) Task: check with Dominic if this can be fixed.
What is the Bow-tie problem? • Hyrax OPeNDAP URL, subsets HDF MOD02 using keyword constraints: http://processing.mne-ngic.mn:8080/opendap/data/hdf4/all_modis/MXD02/MOD02HKM.A2008234.0245.005.2008234025641.hdf.ascii?Latitude[0:4:25][0:4:25] • Server returns the following dataset in ascii: MODIS_SWATH_Type_L1B.Geolocation Fields.Latitude[0], 66.9258, 66.936, MODIS_SWATH_Type_L1B.Geolocation Fields.Latitude[1], 66.8532, 66.8642, MODIS_SWATH_Type_L1B.Geolocation Fields.Latitude[2], 66.7807, 66.7927, MODIS_SWATH_Type_L1B.Geolocation Fields.Latitude[3], 66.802, 66.8129, MODIS_SWATH_Type_L1B.Geolocation Fields.Latitude[4], 66.7296, 66.741, MODIS_SWATH_Type_L1B.Geolocation Fields.Latitude[5], 66.7507, 66.7611, MODIS_SWATH_Type_L1B.Geolocation Fields.Latitude[6], 66.6783, 66.689, • We expect a sort of records with decreasing latitudes, but we get unexpected results (in red).
Bow-tie correction (cont) (left: uncorrected, right: corrected)
Bow-tie limitations • Daily image bow-tie correction on “raw” MOD02/MYD02 1km has no limitation. • However, the HKM and QKM swaths are much larger than the required geo-location fields. • The geo-location field (MXD03) dimensions restrict the image output dimensions as follows: OUTPUT RES = INPUT RES (lines -1 /elems -2) (MXD03)
How to automate in McIDAS? • McIDAS Scheduler • Scheduler calls so-called McBASIC scripts. McBASIC is the BASIC dialect used by McIDAS (Man Computer Interactive Data Access System), a sem-automatic environment for handling large amounts of visual and statistical data, particularly satellite imagery. • Or use ‘cron’ but no visual displays.
McIDAS ICC-NGIC customizations • Complete TMACOMC1.MCB – In creating daily composite of ‘raw’ aka full-swath TERRA/AQUA MODIS • Create McBASIC script for weekly/monthly composite based on daily composite. Using FORM=CMIN conditional minimum to avoid cloudy pixel contamination.
Mcidas customizations • Creating weekly composite of bowtie-corrected TERRA/MODIS and AQUA. We want to automate (because we are lazy) this process. Using IMGOPER command with keyword FORM=CMIN conditional minimum expression. • Create Monthly composite of bowtie-corrected TERRA/MODIS and AQUA, which we again want to automate. Again using FORM=CMIN conditional minimum to avoid cloud contaminated pixels.
CMIN in TOA Reflectance CMIN(0)=(0,-0.1, 0.3,0.1,5….)=0.1 CMINmin new conditional minimum; uses only the term in the equation with the lowest value greater than the specified minimum; if all terms are less than min, the new value is used (default=0 0)
Mcidas customizations • Creating McBASIC script for routinely generating weekly composite of bowtie-corrected NDVI TERRA/MODIS and AQUA. We want to automate this process (because we are lazy). Using IMGOPER command with keyword FORM=CMAX conditional minimum expression. • Same as above but also for Monthly composite from Weekly composite
CMAX in NDVI CMAX(0.7)=(0,1, 0.3,0.8,0.5….)=0.5 CMAXmax new conditional maximum; uses only the term in the equation with the highest value less than the specified maximum; if all terms are greater than max, the new value is used (default=255 255). Theoretical range: -1 to +1 Biological range: -1 to +0.7 Observed range: -1 to >+1
Action list • Learn McIDAS: IMGBOWTIE, IMGOPER (with CMAX and CMIN), McBASIC, McIDAS Scheduler, Backup and File management • Solve vxeos incompatibilities (produce IMAPP compatible L3B products) by either: • Adjust vxeos file HDF formatting of tertiary products (VI’s, Fire detection, etc.), or; • Run IMAPP in parallel to vxeos • Sun glint issues with the NDVI are not observed to happen with other Direct Readout MODIS station(s) in Mongolia. Suggested solution: follow 2b to have redundancy in processing chains and products quality. • Complete McBASIC script for routinely generating weekly composite of bowtie-corrected NDVI TERRA/MODIS and AQUA