View source for Albedo Data and Surface Reflectance Data Visualizations in SNAP using MODIS Imagery


==Introduction==

===Software===

We are going to use '''SNAP 7.0''' to do visualization and data analysis.

===Albedo===
Albedo is the ratio of the amount of irradiance that is reflected by a surface to the total incoming irradiance. Albedo measured on a scale from 0 to 1. For example, a white-colored object reflects most of the incoming light and absorbs very little whereas, a black object will absorb most of the incoming light.

===MODIS===
MODIS is a Moderate Resolution Imaging Spectroradiometer which aboard the Terra and Aqua satellite. 

MODIS has 36 bands in total, but we are going to use band 1 to band 7 in this project. Band wavelength for all 1 to 7 is in nm.

[[File:Modis_band.PNG | 300px]]


According to the wavelength and colour spectrum, Band 1 represents '''Red''', Band 2 represents '''NIR''', Band 3 represents '''Blue''' and Band 4 represents '''Green'''. Band 6 represents '''SWIR'''.


====Pros====

Unlike Landsat or Sentinel product, MODIS provides daily satellite images. Daily images can help us to better understand what is occurring in the earth, or atmosphere. Also, MODIS is commonly using for detecting forest fires, floods or climate changes.

Saving more time. You don't have to do pre-processing like Sentinel 1 images.

====Cons====

Resolution is way lower than Sentinel images. Most MODIS products have 500m or 1km resolution. 

Cloud sensitive, which means MODIS cannot penetrate clouds. Fortunately, MODIS has built-in a cloud mask to remove clouds. Moreover. you can find MODIS image by cloud cover percentage through granule filters when you are downloading the product from the NASA Earth data search. 

[[File:Cloud cover.JPG | 500px]]

==Main Objective==

1. Access and download MODIS imagery

2. Visualize MODIS imagery in SNAP

3. Compute NDSI band from MODIS imagery

4. Know some basic information about MODIS products

5. '''Getting started with MODIS products.'''

==Data Acquisition==

MODIS albedo and surface reflectance products can be download at [https://search.earthdata.nasa.gov/search NASA Earth Data Search]. You have to register an [https://urs.earthdata.nasa.gov/oauth/authorize?response_type=code&client_id=OLpAZlE4HqIOMr0TYqg7UQ&redirect_uri=https%3A%2F%2Fd53njncz5taqi.cloudfront.net%2Furs_callback&state=https%3A%2F%2Fsearch.earthdata.nasa.gov%2Fsearch%2Fgranules%3Fp%3DC1597928934-NOAA_NCEI%26tl%3D1559595988!4!! Earthdata account] before downloading the product. 

MODIS related product are open to public and free to download.

MODIS provides daily satellite images from the year 2000 to the present.

[[File:Earthdata.JPG | 500px]]


In this tutorial, we are going to visualize one albedo product (MOD43A3) and one surface reflectance product (MOD09A1). Both MODIS products are all passed the '''Validation at Stage-3 General Accuracy Statement''', which means both products can be used in '''scientific study'''. '''Validation at stage-3''' means all values have a ± 5% of error and all data should have a zenith angles less than 70 degrees. 

You can view the full description in the product selection page.


[[File:MOD43A3 V006.JPG | 800px]]

[[File:MOD09A1 V006.JPG | 800px]]

Before download the data, you need to understand the products name. For example '''MCD43A3.A2019155.h17v00.006.2019165212344''', MCD43A3 is the project name, A2019155 is the day 155 of 2019 which is June 4, 2019, h17 is the flying orbit number. For arctic region, each days albedo product has four different images, h14 to h17. [[File:Albedo14to17.PNG | 700 px]]

Each images covers a part of arctic area. If you combined image from h14 to h17, you will get the albedo image for whole polar region. For example, [[File:WHOLEPOLAR.PNG | 500 px]]

==Visualization==

===Albedo Visualizations===

Open the albedo product in SNAP and select '''Albedo_WSA_shortwave''' under the band selection. The MODIS albedo product provides three different types of albedo value, bidirectional, black sky albedo and white sky albedo. In this tutorial, we are going to visualize WSA (white sky albedo).

[[file:Shortwavealbedoband.JPG | 500px]]

Go to Colour Manipulation, select "Basic" as editor type and choose a proper colour ramp. For albedo product, we usually choose light colour to represent small albedo area and dark red colour to represent high albedo area.

[[file:Colourselection.JPG | 500px]]

The output map with albedo value is shown below, and you can view each colour correspond to the relative albedo value under (Colour Manipulation - Table). 

[[File:Albedo snap.JPG | 500px]] [[File:Tablealbedo.JPG | 500px]]

====Mask out====

What if you want to know the total area with an albedo value greater than 0.2.

You can go to '''Mask Manager - Creates a new mask based on a value range ''' [[File:Creates a new mask based on a value range.JPG | 40px]] and input the value range.

[[File:Newrangemask.JPG | 500px]]

Then a new mask has been created and you can see the masked out area below. The right image is the original image, and the left image is the masked image. Red areas albedo are greater than 0.2.

[[File:Maskedoutarea.JPG | 500px]] [[File:Albedo snap.JPG | 500px]]


====Ice Albedo Visualization====

Snow and Ice cover land should have a higher albedo than the urban areas. We are going to use '''albedo value in visible band''' as an example to demonstrate snow and ice cover land albedo.

Download an area with snow or ice cover area. For this tutorial, we are using the Ward Hunt ice shelf (Northern coast of Ellesmere Island, Nunavut, Canada) and the Northern coast of Greenland as a study area.

Open the download product in SNAP, and open '''Albedo_WSA_Band1''', '''Albedo_WSA_Band3''', '''Albedo_WSA_Band4''', and '''Albedo_WSA_vis'''. Select '''Albedo_WSA_vis''' and go to '''Colour Manipulation''' [[File:Colourmaniii.PNG | 100 px]] - '''Basic''' - '''Colour ramp''' to change the albedo map from black&white to colours map. You should say some output map like below,

[[File:Whole.PNG | 600 px]]

There are one RGB bands composition to help better visualize albedo images. Select '''band 2 (NIR)''' as Red band, set '''band 1 (red)''' as Green band, and set '''band 4 (green)''' as Blue band [[File:Band214setup.PNG | 300 px]]. The output map is shown below,

[[File:BAND214.PNG | 800 px]]

The output image has a clear outline of the ice-covered area, land area, and snow-covered area.

====Ice Albedo Calculation====

In the MOD43A3 product, they are also provided albedo in each different band. For example, we have albedo from band 1 to band 7. The '''Albedo_WSA_vis''' are calculated by the formula [[File:VISalbedo.PNG | 300 px]]. Wheres, α1 is band 1, α2 is band 4, α3 is band 3 because red is band 1, green is band 4 and blue is band 3. 

You can view each pixels albedo value from pixel info, and use the pixel information to calculate visible bands albedo value. 

[[File:Vis.PNG | 1500 px]]

According to the pixel information above, we can calculate visible bands albedo [[File:Albedoresult.PNG | 400 px]]. For the studied pixel, the calculated visible band albedo is 0.8212 and the provided visible band albedo is 0.8240. According to the user guide, the '''average error''' for visible bands is 4%, so the difference between the calculated result and provided value is within an acceptable range. 

There are two more formulas for NIR bands albedo calculation and shortwave bands albedo calculation. For NIR bands albedo calculation, the formula is [[File:NIRCALCULATION.PNG | 400 PX]]. For shortwave bands albedo calculation, the formula is [[File:SHORTCALCULATION.PNG | 400 px]].

===Surface Reflectance Visualization===

There are lots of RGB band combinations to visualize surface reflectance product. 

Open on the MOD09A1 product in SNAP and select it. Go to '''Window - Open RGB image window''' to select RGB - image channels.

The band combination 1-4-3 creates a true colour image, which represents how the human eye would see the area.  

[[File:143selection.JPG | 500px]]

[[File:143.PNG | 500px]]

The band combination 1-2-1 is used to highlight vegetation. All bright green areas are all vegetations.

[[File:Band121.JPG | 500px]]


The band combination 3-6-7 highlights snow and ice.

In this band combination, the snow is shown as red, clouds are white, vegetation is green, ice is pinkish-red on the landscape.

[[File:367.JPG | 500px]]

The band combination 1-4-6 also highlights the snow-covered area. The snow-covered area shows yellow on the landscape.

[[File:Band146.JPG | 500px]]

There are more band combinations such as 7-2-1 or 2-1-1.

The band combination 7-2-1 is used to observe fire burn scars on the surface from fires. Burn Scars are shown as red to brown.

The band combination 2-1-1 creates a false-color image. Vegetation is green, clouds appear white.

==NDSI==

The remote sensing index can be calculated by the MODIS image. 

NDSI is the Normalized Difference Snow Index. NDSI is commonly used in cloud or snow differentiation. Cloud and snow are all shown as bright colour in a true colour image, however, snow is getting darker when the wavelength is greater than 1.6mm. According to the MODIS band wavelength above, for MODIS data the NDSI is calculated as:

                    NDSI = (Band 4 - Band 6) / (Band 4 + Band 6) 

The calculation can be done by '''Band Maths''' which under '''Raster'''.

[[File:Bandmath123.JPG | 800px]]

Areas with NDSI > 0.4 are snow area. In the output map, the white area is shown as the snow-covered surface.

[[File:NDSIGZ.JPG | 500px]]


===NDSI analysis===

If you want to know the NDSI for a specific region or you want to calculate the mean, of maximum NDSI, you can use '''Analysis - Statistics''' to see those values. 

In this tutorial, just create a polygon or square [[File:Rectangulardrawing.PNG | 100 px]] above the image by the '''rectangular drawing tool''' and go to '''Mask Manager''' to select the mask you just created.  If you want to know the average NDSI for the Ottawa area, download the Ottawa boundary shapefile and go to '''Vector - Import - ESRI Shapefile''' to import the vector mask.

[[File:NDSImaskmanager.JPG | 500px]]

Then open the '''Analysis - Statistics''' window to see the statistics value.

[[File:NDSISTATISTICSVALUE.JPG | 1300px]]

==Other Product==

===NDVI===

NDVI (normalized difference vegetation index) can also be computed by MODIS surface reflectance imagery. The formula is NDVI = (band 2 - band 1) / (band 2 + band 1). 

However, MODIS has three different resolution (250m, 500m, 1km) products contain NDVI value. The product names are '''MOD13Q1, MOD13A1, MOD13A2'''. You can view more information from [https://modis.gsfc.nasa.gov/data/dataprod/mod13.php MODIS Vegetation Index Products]

===Snow cover===

The product '''MOD10A1''' identifies snow-covered land and snow-covered ice on inland water.

===Surface Temperature===

The 31st and 32nd bands of MODIS are using for surface and cloud temperature monitoring, and the bandwidth for the 31st band is 10.780 - 11.280 µm, for the 32nd band is 11.770 - 12.270 µm.

Therefore, MODIS has several products with surface temperature value. For example, '''MOD11A1''' for daily surface temperature, '''MOD11A2''' for 8-day surface temperature

==Reference==

NASA Products list: https://modis.gsfc.nasa.gov/data/dataprod/

MOD09A1 Product description: https://lpdaac.usgs.gov/products/myd09a1v006/

Getting Started with MODIS V6 Surface Reflectance Data: https://lpdaac.usgs.gov/resources/e-learning/getting-started-modis-v6-surface-reflectance-data-part-2-interpreting-quality-information/

Validating MODIS land surface reflectance and albedo products: methods and preliminary results : http://cursosihlla.bdh.org.ar/ET_Local_Regional_2014/Clase_5/3_Material_Lectura/2002_Liang%20et%20al._ALBEDO.pdf