MOD04 - L2 Product Info From Http://modis-Atmos - Gsfc.nasa - Gov

MOD04_L2 Product Info
From http://modis-atmos.gsfc.nasa.gov/
Product Description
The MODIS Aerosol Product monitors the ambient aerosol optical thickness over the
oceans globally and over a portion of the continents. Further, the aerosol size distribution
is derived over the oceans, and the aerosol type is derived over the continents. Daily
Level 2 data are produced at the spatial resolution of a 10x10 1-km (at nadir)-pixel array.
There are two MODIS Aerosol data product files: MOD04_L2, containing data collected
from the Terra platform; and MYD04_L2, containing data collected from the Aqua
platform.
Research and Application
Aerosols are one of the greatest sources of uncertainty in
climate modeling. Aerosols vary in time in space and can lead
to variations in cloud microphysics, which could impact cloud
radiative properties and climate. The MODIS aerosol product is
used to study aerosol climatology, sources and sinks of specific
aerosol types (e.g., sulfates and biomass-burning aerosol),
interaction of aerosols with clouds, and atmospheric corrections of remotely sensed
surface reflectance over the land.
Data Set Evolution
Prior to MODIS, satellite measurements were limited to reflectance measurements in one
(GOES, METEOSAT) or two (AVHRR) channels. There was no real attempt to retrieve
aerosol content over land on a global scale. Algorithms had been developed for use only
over dark vegetation. The blue channel on MODIS, not present on AVHRR, offers the
possibility to extend the derivation of optical thickness over land to additional surfaces.
The algorithms will use MODIS bands 1 through 7 and 20 and require prior cloud
screening using MODIS data. Over the land, the dynamic aerosol models will be derived
from ground-based sky measurements and used in the net retrieval process.
Over the ocean, three parameters that describe the aerosol loading and size distribution
will be retrieved. Pre-assumptions on the general structure of the size distribution are
required in the inversion of MODIS data, and the volume-size distribution will be
described with two log-normal modes: a single mode to describe the accumulation mode
particles (radius < 0.5 m) and a single coarse mode to describe dust and/or salt particles
(radius > 1.0 m). The aerosol parameters we therefore expect to retrieve are: the ratio
between the two modes, the spectral optical thickness, and the mean particle size.
The quality control of these products will be based on comparison with ground stations
and climatology.
File Format Basics

MOD04_L2 product files are stored in Hierarchical Data Format (HDF). HDF is a
multi-object file format for sharing scientific data in multi-platform distributed
environments. HDF files should only be accessed through HDF library subroutine and
function calls, which can be downloaded from the HDF web site*. Each of the 53 gridded
parameters listed below is stored as a Scientific Data Set (SDS) within the HDF file.
* NOTE: MATLAB provides an HDF toolbox which allows the user to read HDF files
content. There are other softwares that can be downloaded from the NASA web sites to
view and display this type of files.
MOD04_L2 Dimension List
1. Cell_Along_Swath = 203 (typical size)
2. Cell_Across_Swath = 135 (typical size)
3. Solution_1_Land = 2
6. Solution_Ocean = 2
7. MODIS_Band_Land = 7
8. MODIS_Band_Ocean = 7
9. QA_Byte_Land = 5
10. QA_Byte_Ocean = 5
11. Solution_Index = 9
12. Number_of_Instrument_Scans = 203 (typical size)
13. Maximum_Number_of_1km_Frames = 1354 (typical size)
MOD04_L2 Scientific Data Set (SDS) List

Geolocation and Time Parameters
1.
Longitude
Description: Geodetic Longitude
Dimensions: (Cell_Along_Swath, Cell_Across_Swath)
Valid Range: -180 to +180 degrees east
2.
Latitude
Description: Geodetic Latitude
Valid Range: -90 to +90 degrees north
3.
Scan_Start_Time
Description: International Atomic Time at Start of Scan replicated across the Swath
Valid Range: 0.0 to 3.1558E+9 seconds since 1 January 1993 00:00:00
Solar and Viewing Geometry Parameters

4.
Solar_Zenith
Description: Solar Zenith Angle, Cell to Sun
Valid Range: 0 to +180 degrees
5.
Solar_Azimuth
Description: Solar Azimuth Angle, Cell to Sun
Valid Range: -180 to +180 degrees
6.
Sensor_Zenith
Description: Sensor Zenith Angle, Cell to Sensor
Valid Range: 0 to 180 degrees
7.
Sensor_Azimuth
Description: Sensor Azimuth Angle, Cell to Sensor
Valid Range: -180 to 180 degrees
Science Parameters
Combined Land and Ocean
8.
Scattering_Angle
Description: Scattering Angle
Valid Range: 0 to 180 degrees
9.
Optical_Depth_Land_And_Ocean
Description: Aerosol Optical Thickness at 0.55 m for both Ocean (best) and Land
(corrected)
Valid Range: 0 to 3
Land Only
10.
Aerosol_Type_Land
Description: Aerosol Type
Valid Range: 0 to 4 (0=mixed, 1=dust, 2=sulfate, 3=smoke, 4=heavy absorbing smoke)
11.
Continental_Optical_Depth_Land
Description: Continental Optical Thickness at 0.47 and 0.66 m
Dimensions: (Solution_1_Land, Cell_Along_Swath, Cell_Across_Swath)
Valid Range: 0 to 3
12.
Corrected_Optical_Depth_Land
Description: Corrected Optical Thickness at 0.47, 0.55, and 0.66 m
Valid Range: 0 to 3
13.
Estimated_Uncertainty_Land
Description: Uncertainty of Optical Thickness at 0.47 and 0.66 m
Valid Range: 0 to 2
14.
Mass_Concentration_Land
Description: Mass Concentration over Land
Valid Range: 0 to 1000 x 1.0e-6g/cm^2
15.
Angstrom_Exponent_Land
Description: Angstrom Exponent at 0.47 and 0.67 m

Valid Range: -0.5 to 3
16.
Reflected_Flux_Land
Description: Normalized Reflected Flux at 0.47 and 0.66 m
Valid Range: 0 to 1
17.
Transmitted_Flux_Land
Description: Normalized Transmitted Flux at 0.47 and 0.66 m
Valid Range: 0 to 1
18.
Cloud_Fraction_Land
Description: Cloud Fraction (%)
Valid Range: 0 to 100%
19.
Dust_Weighting_Factor_Land
Description: Dust Aerosol Weighting Factor
Valid Range: 0 to 1
20.
Number_Pixels_Percentile_Land
Description: Number of Pixels in Desired Percentile
Valid Range: 0 to 400
21.
Mean_Reflectance_Land
Description: Mean Reflectance at 5 bands
Dimensions: (MODIS_Band_Land, Cell_Along_Swath, Cell_Across_Swath)
Valid Range: 0 to 1
22.
STD_Reflectance_Land
Description: Standard Deviation of Reflectance at 5 bands
Dimensions: (MODIS_Band_Land, Cell_Along_Swath, Cell_Across_Swath)
Valid Range: 0 to 2
Ocean Only
23.
Solution_Index_Ocean_Small
Description: Solution number index (1 through 4) for small aerosol particles (for best and
average solutions). Indices of ocean models 1 through 4 correspond to accumulation
(small) mode models with effective radii 0.10, 0.15, 0.20, 0.25 m, respectively.
Dimensions: (Solution_Ocean, Cell_Along_Swath, Cell_Across_Swath)
Valid Range: 1 to 5
24.
Solution_Index_Ocean_Large
Description: Solution number index (5 through 9) for large aerosol particles (for best and
average solutions). Indices of ocean models 5 through 7 correspond to coarse (large)
mode models of marine (sea salt) particles with effective radii 1.0, 1.5, 2.0 m,
respectively. Indices of ocean models 8 and 9 correspond to coarse (large) mode models
of mineral dust particles with effective radii 1.5 and 2.5 m, respectively.
25.
Effective_Optical_Depth_Best_Ocean
Description: Aerosol Optical Thickness for Best Solution at 7 bands
Dimensions: (MODIS_Band_Ocean, Cell_Along_Swath, Cell_Across_Swath)
Valid Range: 0 to 3
26.
Effective_Optical_Depth_Average_Ocean
Description: Aerosol Optical Thickness for Average Solution at 7 bands
Valid Range: 0 to 3
27.
Optical_Depth_Small_Best_Ocean
Description: Aerosol Optical Thickness for Small Mode of Best Solution at 7 bands
Valid Range: 0 to 3
28.
Optical_Depth_Small_Average_Ocean
Description: Aerosol Optical Thickness for Small Mode of Average Solution at 7 bands
Valid Range: 0 to 3
29.
Optical_Depth_Large_Best_Ocean
Description: Aerosol Optical Thickness for Large Mode of Best Solution at 7 bands
Valid Range: 0 to 3
30.
Optical_Depth_Large_Average_Ocean
Description: Aerosol Optical Thickness at 7 bands for large mode of average solution
Valid Range: 0 to 3
31.
Mass_Concentration_Ocean
Description: Mass Concentration for Best and Average Solutions
Valid Range: 0 to 1000 x 1.0e-6g/cm^2
32.
Effective_Radius_Ocean
Description: Effective Radius of Both Solutions at 0.55 m
Valid Range: 0 to 5 m
33.
Cloud_Condensation_Nuclei_Ocean
Description: Column Number of CCN of Both Solutions at 0.55 m
Valid Range: 0.0 to 10000000000.0 CCN/cm^2
34.
Asymmetry_Factor_Best_Ocean
Description: Asymmetry Factor for Best Solution at 7 bands
Valid Range: 0 to 3
35.
Asymmetry_Factor_Average_Ocean
Description: Asymmetry Factor for Average Solution at 7 bands
Valid Range: 0 to 3
36.
Backscattering_Ratio_Best_Ocean
Description: Backscattering Ratio of Best Solution at 7 bands
Valid Range: 0 to 3
37.
Backscattering_Ratio_Average_Ocean
Description: Backscattering Ratio of Average Solution at 7 bands
Valid Range: 0 to 3
38.
Angstrom_Exponent_1_Ocean
Description: Angstrom Exponent for 0.550 and 0.865 m
39.
Angstrom_Exponent_2_Ocean
Description: Angstrom Exponent for 0.865 and 2.130 m

40.
Reflected_Flux_Best_Ocean
Description: Normalized Reflected Flux of Best Solution at 7 bands
Valid Range: 0 to 1
41.
Reflected_Flux_Average_Ocean
Description: Normalized Reflected Flux of Average Solution at 7 bands
Valid Range: 0 to 1
42.
Transmitted_Flux_Best_Ocean
Description: Normalized Transmitted Flux of Best Solution at 7 bands
Valid Range: 0 to 1
43.
Transmitted_Flux_Average_Ocean
Description: Normalized Transmitted Flux of Average Solution at 7 bands
Valid Range: 0 to 1
44.
Least_Squares_Error_Ocean
Description: Least Square Error Estimate
Valid Range: 0 to 1
45.
Small_Mode_Weighting_Ocean
Description: Small Mode Weighting Factor
Valid Range: 0 to 1
46.
Optical_Depth_Ratio_Small_Ocean
Description: Ratio of Optical Depth of Small Mode vs Effective Optical Depth at .55 m
Valid Range: 0 to 1
47.
Cloud_Fraction_Ocean
Description: Cloud Fraction in Percentage
Valid Range: 0 to 100%
48.
Number_Pixels_Used_Ocean
Description: Number of Pixels used for 0.55 m solution
49.
Mean_Reflectance_Ocean
Description: Mean Reflectances at 7 bands
Valid Range: 0 to 1
50.
STD_Reflectance_Ocean
Description: Standard Deviation of Reflectances at 7 bands
Valid Range: 0 to 2
Cloud Mask and QA Parameters

51.
Cloud_Mask_QA
Description: Cloud Mask QA at 10x10 km resolution
Valid Range: (bit mask)
52.
Quality_Assurance_Land
Description: Run-Time QA Flags for Land
Dimensions: (Cell_Along_Swath, Cell_Across_Swath, QA_Byte_Land)
53.
Quality_Assurance_Ocean
Description: Run-Time QA Flags for Ocean
Dimensions: (Cell_Along_Swath, Cell_Across_Swath, QA_Byte_Ocean)
Product Grid
The MODIS Aerosol (MOD04_L2) product contains data that has a spatial resolution
(pixel size) of 10 x 10 kilometers (at nadir). Each MOD04_L2 product file covers a fiveminute time interval (based on the start time of each MODIS Level-1B granule), which
means the MOD04_L2 output grid is 135 10-km (at nadir) pixels in width and 203 10-km
(nadir) pixels in length for nine consecutive granules. Every tenth granule has an output
grid size of 135 by 204 pixels.
Granule Coverage
The image below depicts MODIS Level-2 granule coverage during a single simulated
orbit. It should be noted that a granule of Level-2 MODIS data is defined as a single
Level-2 MODIS product HDF file.
Computation of 10-km Geolocation in L2 Atmosphere Products

The geolocation in the 04_L2 MODIS product is at 10-km resolution. This geolocation is
generated from 10x10 1-km L1B input and is computed by averaging the 4 central (5,5),
(5,6), (6,5), (6.,6) 1-km L1B input pixels in each 10x10 km area.
Sample File
A sample (Collection 005) MYD04_L2 HDF file (1.8MB), a fully populated daytime
granule, is available for download. The image (inset) was created from the L1B radiance
data for the same granule (RGB=1:4:3) 3 April 2007 at 2140 UTC. The sample HDF data
file can be used to familiarize yourself with the data format and/or as input to MODIS
visualization and analysis software on your local platform.
MYD04_L2 HDF File
Ordering Data
All MODIS Atmosphere data products are available to
the public (at no charge) through the Level 1 and
Atmosphere Archive and Distribution System (LAADS).
This new and user-friendly MODIS data ordering system gives the user convenient
means to simultaneously order several MODIS Data Sets, including Geolocation. This
system also works well for single products.
Web site: http://ladsweb.nascom.nasa.gov/data/search.html
HDF Product File

MODIS Level-2 Hierarchical Data Format (HDF) product files have standardized
filenames. The prefix MOD is reserved for files containing data collected from the Terra
(AM overpass) platform and MYD is reserved for files containing data collected from the
Aqua (PM overpass) platform.
Terra (AM) Platform:
Aqua (PM) Platform:
Note that:
all times are UTC time, not local time
A number of tools for working with MODIS HDF data are available for
download from this web site. These tools fall into three major categories,
which are listed below. Under each main category, more detail on the
function of each tool is offered.
1. Granule Locator Tools
1.
Locating Level-2 Granules from a Generic Global Map
2.
Locating Level-2 Granules from a Level-3 Global Image
3.
Locating Level-1B Granules using Text Input
4.
Locating Level-1B Granules using Level-1B Metadata
5.
Locating Level-2 Granules using Text Input
2. Spatial and Dataset Subsetting

1.
Subsetting Single Resolution HDF Files
2.
Subsetting Double Resolution HDF Files
3. Visualization & Analysis

1.
Visualizing HDF Data
2.
Extracting and Visualizing Bit Flags from Byte Data
3.
Creating the Atm. Standard Color Scale & Bar
4.
Creating Browse Images from Level-3 HDF Data
5.
Analyzing HDF Data
To learn more about and acquire these tools, visit the MODIS-Atmosphere:
http://modis-atmos.gsfc.nasa.gov/tools.html
Important Tables
TABLE 1: CHARACTERISTICS OF MODIS CHANNELS USED IN THE AEROSOL RETRIEVAL
Band Number
Bandwidth (m)
Weighted Central Wavelength

(m)
Resolution (m)
0.620 - 0.670
0.646
250
0.841 - 0.876
0.855
250
0.459 - 0.479
0.466
500
0.545 - 0.565
0.553
500
1.230 1.250
1.243
500
1.628 1.652
1.632
500
2.105 2.155
2.119
500
TABLE 2: CONTENTS OF MODIS C005 AEROSOL LEVEL 2 FILE (MOD04/MYD04): OCEAN

PRODUCTS
Name of Product (SDS)
Dimesions: 3rd Dimension
Type of product
Effective_Optical_Depth_Average_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Retrieved Primary
Effective_Optical_Depth_Best_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Retrieved Primary
Optical_Depth_Ratio_Small_Ocean_0_55micron
X,Y,2: average, best
Retrieved Primary
Solution_Index_Ocean_Small
Retrieved Primary
Solution_Index_Ocean_Large
Retrieved Primary
Least_Squares_Error_Ocean
Retrieved Diagnostic
Effective_Radius_Ocean
Derived
Optical_Depth_Small_Best_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Optical_Depth_Small_Average_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Optical_Depth_Large_Best_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Optical_Depth_Large_Average_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Mass_Concentration_Ocean
Derived
Derived
Asymmetry_Factor_Best_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Asymmetry_Factor_Average_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Backscattering_Ratio_Best_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Backscattering_Ratio_Average_Ocean
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Derived
Angstrom_Exponent_1_Ocean(0.55/0.86 micron)
Derived
Angstrom_Exponent_2_Ocean(0.86/2.1 micron)
Derived
Derived
Optical_Depth_by_models_ocean
X,Y,9: 9 models
Derived
Cloud_Fraction_Ocean
X,Y:
Diagnostic
Number_Pixels_Used_Ocean
X,Y:
Diagnostic
Mean_Reflectance_Ocean
X,Y:
Diagnostic
STD_Reflectance_Ocean
X,Y:
Diagnostic
Aerosol_Cldmask_Byproducts_Ocean
X,Y:
Diagnostic
Quality_Assurance_Ocean
X,Y,5 bytes
Diagnostic
X,Y: 0.55m
Joint (QAC0) **
Image_Optical_Depth_Land_And_Ocean
X,Y: 0.55m
Joint (QAC0)
Optical_Depth_Ratio_Small_Land_And_Ocean
X,Y: 0.55m
Joint (QAC0)
X = 135; Y = 203. If there is a 3rd dimension of the SDS, then the indices of it are given. The Retrieved
parameters are the solution to the inversion,whereas Derived parameters follow from the choice of
solution. Diagnostic parameters aid in understanding of the directly Retrieved or Derived products.
Experimental products are unrelated to the inversion but may have future applications. Joint products
are the combined land and ocean products, with associated QAC constraint (for over ocean) in parentheses.
**Based on evaluation of operational C005-O data, the QAC for quantitative studies should be limited to
QAC1 only.
Some of the ocean products are combined with products from land (discussed in the next section) as the
Joint products. For AOD, two joint products are reported, the Optical_Depth_Land_And_Ocean and the
Image_Optical_Depth_Land_And_Ocean. The first product is supposed to have more quantitative
meaning, so is constrained by QAC. In practice, however, the two joint products are identical over ocean,
meaning that the values of Effective_Optical_Depth_Average_Ocean (at 0.55m) are written into both
SDSs, regardless of QAC. The Optical_Depth_Ratio_Small_Ocean product is copied (regardless of QAC)
into Optical_Depth_Ratio_Small_Land_And_Ocean. Quality Assurance confidence (QAC) value ranges
from 0 (bad quality) to 3 (good quality).
A Retrieved parameter is one that is a solution to the AOD retrieval procedure. Derived parameters are
computed based on products directly retrieved. For example, the ngstrom Exponent is derived based on
the spectral AOD that characterizes the retrieved solution. Products that are Diagnostic include QA
parameters and those parameters that were calculated during intermediate steps. These diagnostic
parameters can be used to understand how the retrieval worked. Products denoted Experimental are
superfluous to the main retrieval procedure, but are useful in other applications.
Aerosol models are: Continental, Moderately Absorbing/Developing World, Non-absorbing/UrbanIndustrial, Absorbing/Smoke, Spheroid and Dust.
TABLE 3: CONTENTS OF MODIS C005 AEROSOL LEVEL 2 FILE (MOD04/MYD04): LAND
PRODUCTS
Name of Product (SDS)
Dimesions: 3rd Dimension
Type of product
X,Y,3: 0.47, 0.55, 0.66 m
Retrieved Primary
Corrected_Optical_Depth_Land_wav2p1
X,Y,1: 2.12 m
Retrieved Primary
Optical_Depth_Ratio_Small_Land
X,Y: (for 0.55 m)
Retrieved Primary
Surface_Reflectance_Land
X,Y,3: 0.47, 0.66, 2.12 m
Retrieved Primary
Fitting_Error_Land
X,Y: (at 0.66 m)
Retrieved By-Product
X,Y,5: 5 bytes
Diagnostic
Aerosol_Type_Land
X,Y:
Diagnostic
X,Y: (for 0.66/0.47 m)
Derived
Mass_Concentration_Land
X,Y:
Derived
Optical_Depth_Small_Land
X,Y,4: 0.47,0.55,0.66,2.12 m
Derived
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Diagnostic
STD_Reflectance_Land
X,Y,7: 0.47,0.55,0.66,0.86,1.2,1.6,2.12m
Diagnostic
Cloud_Fraction_Land
X,Y:
Diagnostic
Number_Pixels_Used_Land
X,Y:
Diagnostic
Path_Radiance_Land
X,Y,2: 0.47, 0.66 m
Experimental
Error_Path_Radiance_Land
X,Y,2: 0.47, 0.66 m
Experimental
Critical_Reflectance_Land
X,Y,2: 0.47, 0.66 m
Experimental
Error_Critical_Reflectance_Land
X,Y,2: 0.47, 0.66 m
Experimental
Quality_Weight_Path_Radiance_Land
X,Y,2: 0.47, 0.66 m
Experimental
Quality_Weight_Crit_Reflectance_Land
X,Y,2: 0.47, 0.66 m
Experimental
X,Y: (for 0.55 m)
Joint (QAC0) **
Image_Optical_Depth_Land_And_Ocean
X,Y: (for 0.55 m)
Joint (QAC0)
Optical_Depth_Ratio_Small_Land_And_Ocean
X,Y: (for 0.55 m)
Joint (QAC0)
From MYD04_L2(MODIS AQUA)

Deep_Blue_Aerosol_Optical_Depth_550_Land
Deep_Blue_Aerosol_Optical_Depth_Land
Deep_Blue_Angstrom_Exponent_Land
Deep_Blue_Single_Scattering_Albedo_Land
Deep_Blue_Surface_Reflectance_Land
X,Y:(for 0.55 m)
X,Y,3: 0.412, 0.47, and 0.66m
X,Y:(0.412-0.47) m
X,Y,3: 0.412, 0.47, and 0.66m
X,Y,3: 0.412, 0.47, and 0.66m
Derived
Derived
Derived
Derived
Derived
How to read MODIS AOD in Matlab:

1-To read from the hdf file:
When reading one file:
fname='filename.hdf';
hdfread(fname,'Parameter_name');
When reading several files
fid=fopen('filename_list.txt','r');
for loop with the size of filename_list
fname=fscanf(fid,'%s',1);
hdfread(fname,'Parameter_name');
2-Parameters from the hdf file that should be read:
TABLE 4: Used SDS from MOD04/MYD04
Parameter
Longitude
Latitude
Aerosol_Type_Land
Sensor_Zenith
Sensor_Azimuth
Solar_Zenith
Solar_Azimuth
Scattering_Angle
Surface_Reflectance_Land
Optical_Depth_Small_Land
Critical_Reflectance_Land
Path_Radiance_Land
Corrected_Optical_Depth_Land_wav2p1
Cloud_Mask_QA
Dimensions
Fill Value
Scale Factor
203x135
203x135
203x135
203x135
203x135
203x135
203x135
203x135
3x203x135
7x203x135
3x203x135
4x203x135
2x203x135
2x203x135
203x135
203x135
203x135x5 bytes
203x135x5
-999
-999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
-9999
none
none
1
1
1
0.01
0.01
0.01
0.01
0.01
0.001
0.0001
0.001
0.001
0.0001
0.0001
0.001
0.001
1
1
Description of the QA flags:

TABLE 5: Cloud_mask_QA flags
QA Flag Name
Cloud Mask
Number of Bits
1
Cloud Mask Quality Flag
Day/Night flag
Sun glint flag
Snow/Ice flag
Land/Water flag
Bit Value
0
1
0
1
2
3
0
1
0
1
0
1
0
1
Description
Undetermined
Determined
0-25% Cloudy pixels
25-50% cloudy pixels
50-75% cloudy pixels
75-100%cloudy pixels
Night
Day
Yes
No
Yes
No
Water (ocean)
Coastal
2
Desert
3
Land
TABLE 6: QUALITY ASSURANCE LAND FLAGS (5 bytes = 40 bits )
Flag Name
Number of Bits
Bit Value
Description
0.47 m Aerosol Optical
1
0
Not useful
Thickness Usefulness Flag
1
Useful
3
0
No Confidence(or Fill)
Thickness Confidence Flag
1
Marginal
2
Good
3
Very good
1
0
Not useful
Thickness Usefulness Flag
1
Useful
3
0
Thickness Confidence Flag
1
Marginal
2
Good
3
Very good
processing path flags
Dark Target Criteria
3
0
not met (Fill Value)
used in retrieval
1
0.01 < Ref (2.1 m) 0.05
2
0.05 < Ref (2.1 m) 0.10
3
0.10 < Ref (2.1 m) 0.15
4
0.15 < Ref (2.1 m) 0.25
5
0.25 < Ref (2.1 m) 0.40
Error Code
3
0
No error
(when Fill Values are
1
So Solar and illumination angles out-of-bounds in look-up
assigned)
table
2
Apparent reflectance out-of-bounds in look-up table
3
Number of cloud and water free pixels not met
4
Thresholds of 2.1 m not met
5
Thresholds of 3.8 m not met
6
Thin cirrus detection not met
High Solar Zenith Angle
1
0
No
( > 72 )
1
Yes
Increased Spatial Resolution
1
0
No
( 5x5 km )
1
Yes
Aerosol Type
2
0
Mixed
(over Land only)
1
Dust
2
Sulfate
3
Smoke
Thin Cirrus or
2
0
0 < (1.38 m) < 0.01; Correction is done
Stratospheric Aerosol Index
1
(1.38 m) < 0; No correction
2
(0.66 m) < 0.04; No correction
3
(1.38 m) > 0.01; No correction
input data resource flags
Total ozone
2
0
TOVS
1
TOMS
2
Climatology
3
GMAO
Total precipitable water
2
0
NCEP / GDAS
1
MOD05* NIR (MODIS Near-IR Water Vapor Retrieval)
2
Climatology
3
GMAO
Snow cover
2
0
MOD35* (MODIS Cloud Mask)
1
MOD10* (MODIS Eight-day Snow Cover)
Spares
6
TBD
Deep Blue Aerosol
1
0
Not useful
Usefulness Flag
1
Useful
Deep Blue Aerosol
2
0
Confidence Flag
1
Marginal
2
Good
3
Very good
Deep Blue Aerosol Type
2
0
Mixed
Note: Flags 2 and 3 are
1
Dust
reversed from the Aerosol
2
Smoke
Type (over land only) above
3
Sulfate
Deep Blue Aerosol
2
0
Optimal Retrieval Performed
Retrieving Condition
1
White Sand
2
Cloudy
3
Spare
(550 nm) > 5.0 (Out of Bounds)

TBD
Notes on the QA flags:

Please remember that when you read these flags, theyre in uint8 type. So what you need to do is the
following:
Youre gonna change the value to binary. Use the function dec2bin in Matlab to do the conversion.
Remember that you have to indicate whether you want to convert the number to 8 bits or 16 bits, etc.
For example:
Variable_binary=dec2bin(variable_decimal,16);
In the case of Cloud_Mask_QA, you have 8 bits.
cmqa=reshape(dec2bin(double(Cloud_Mask_QA),8),203,135,8);
cmqa is a variable of size 203x135x8 and if I want to see what is bit 0 I access (:,:,1) for example.
Also remember that bits are represented in Little-Endian.
In the case of Quality_Assurance_Land, you have 40 bits but theyre distributed in the following way:
The first 8 bits as indicated in the table by yellow correspond to (:,:,1)
Next 8 bits as indicated in the table by orange correspond to (:,:,2)
Again the ones indicated by green correspond to (:,:,3)
The ones indicated by blue correspond to (:,:,4) and finally the ones indicated by lavender correspond to
(:,:,5).
Qal1=reshape(dec2bin(double(Quality_Assurance_Land(:,:,1)),8),203,135,8)
;
3-For every file we need to calculate the minimum distance to find the closest pixel to the AERONET site.
First, we define latitude and longitude coordinates for the AERONET Site (the AERONET header
provides that information). Since we have the variables Longitude and Latitude from the hdf file, we can
calculate the minimum distance in the following way:
x1=longitude_AERONET
x=double(longitude_hdf)
y1=latitude_AERONET
y=double(latitude_hdf)
[R,C]= min[((x-x1)2+(y-y1)2)1/2]
[K,L]=min(R);
row_index=C(L);
column_index=L;
4-After the row and column indexes are calculated, do the following:
Assign a temporary variable for each parameter and take into account the following cases:
-Check if the [row index >203 or <1] or [column index >135 or <1], if this is the case then each
temporary variable should be set to NaN.
-For fill values, if any of the parameters at (row_index,column_index)= -9999, set them to
NaN.
-Else, temporary_variable=double(parameter(row_index,column_index))*scale_factor
5-For each parameter, create a variable that will get the temporary assigned value for each file.
For example:
AOD(1:3,file)=AODtemp;
MeanReflectance(1:7,file)=MeanReflectancetemp;
SensorZenith(file)=SensorZenithtemp;

MOD04 - L2 Product Info From Http://modis-Atmos - Gsfc.nasa - Gov

Uploaded by

Copyright:

Available Formats

MOD04 - L2 Product Info From Http://modis-Atmos - Gsfc.nasa - Gov

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

MOD04 - L2 Product Info From Http://modis-Atmos - Gsfc.nasa - Gov

Uploaded by

Copyright:

Available Formats

MOD04_L2 Product Info

File Format Basics

MOD04_L2 Scientific Data Set (SDS) List

Solar and Viewing Geometry Parameters

Description: Angstrom Exponent at 0.47 and 0.67 m

Description: Angstrom Exponent for 0.865 and 2.130 m

Cloud Mask and QA Parameters

Computation of 10-km Geolocation in L2 Atmosphere Products

HDF Product File

all times are UTC time, not local time

Locating Level-2 Granules from a Generic Global Map

Locating Level-2 Granules from a Level-3 Global Image

Locating Level-1B Granules using Text Input

Locating Level-1B Granules using Level-1B Metadata

Locating Level-2 Granules using Text Input

2. Spatial and Dataset Subsetting

Subsetting Single Resolution HDF Files

Subsetting Double Resolution HDF Files

3. Visualization & Analysis

Visualizing HDF Data

Extracting and Visualizing Bit Flags from Byte Data

Creating the Atm. Standard Color Scale & Bar

Creating Browse Images from Level-3 HDF Data

Analyzing HDF Data

Weighted Central Wavelength

TABLE 2: CONTENTS OF MODIS C005 AEROSOL LEVEL 2 FILE (MOD04/MYD04): OCEAN

Dimesions: 3rd Dimension

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,2: average, best

X,Y,3: 0.47, 0.55, 0.66 m

X,Y: (for 0.55 m)

X,Y,3: 0.47, 0.66, 2.12 m

X,Y: (at 0.66 m)

X,Y: (for 0.66/0.47 m)

X,Y,2: 0.47, 0.66 m

X,Y,2: 0.47, 0.66 m

X,Y,2: 0.47, 0.66 m

X,Y,2: 0.47, 0.66 m

X,Y,2: 0.47, 0.66 m

X,Y,2: 0.47, 0.66 m

X,Y: (for 0.55 m)

X,Y: (for 0.55 m)

X,Y: (for 0.55 m)

From MYD04_L2(MODIS AQUA)

How to read MODIS AOD in Matlab:

Description of the QA flags:

Cloud Mask Quality Flag

Sun glint flag

(550 nm) > 5.0 (Out of Bounds)

Notes on the QA flags:

You might also like