Quantify the Impacts of a Fire Using MODIS and Landsat Remote Sensing Data in Python


Welcome to Week 9!

Welcome to week 9 of Earth Analytics! This week you will work with MODIS and Landsat data to calculate burn indices.

Download Cold Springs Fire Data (~150 MB)

Or use earthpy et.data.get_data('cold-springs-fire')

Materials to Review For This Week’s Assignment

Please be sure to review:

  1. Chapter 11 on Calculating Normalized Burn Ratio (NBR) in Section 5 of the Intermediate Earth Data Science Textbook and
  2. Data Story on Cold Springs Fire in Section 7 of the Intermediate Earth Data Science Textbook

1. Complete the Assignment Using the Template for this week.

Note that this assignment is worth more points than a usual weekly assignment.

Click here to view the the example GitHub Repo with the assignment template.

2. Suggested Fire Readings

Please read the articles below to prepare for next week’s class.

Downloading from https://ndownloader.figshare.com/files/10960211?private_link=18f892d9f3645344b2fe
Extracted output to /root/earth-analytics/data/cs-test-naip/.
Downloading from https://ndownloader.figshare.com/files/10960109
Extracted output to /root/earth-analytics/data/cold-springs-fire/.
Downloading from https://ndownloader.figshare.com/files/10960112
Extracted output to /root/earth-analytics/data/cold-springs-modis-h4/.
Downloading from https://ndownloader.figshare.com/files/21941085
Extracted output to /root/earth-analytics/data/earthpy-downloads/landsat-coldsprings-hw

Homework Figure 1 - Grid of 3 Color InfraRed (CIR) Post-Fire Plots: NAIP, Landsat and MODIS

CIR Composite images from NAIP, Landsat, and MODIS for the post-Cold Springs fire.
CIR Composite images from NAIP, Landsat, and MODIS for the post-Cold Springs fire.

Homework Figure 3 - Difference NBR (dNBR) Using Landsat & MODIS

<xarray.DataArray (y: 44, x: 113)>
array([[2., 2., 2., ..., 2., 2., 2.],
       [2., 2., 2., ..., 2., 2., 2.],
       [2., 2., 2., ..., 2., 2., 2.],
       ...,
       [2., 2., 2., ..., 2., 2., 2.],
       [2., 2., 2., ..., 2., 2., 2.],
       [2., 2., 2., ..., 2., 2., 2.]])
Coordinates:
  * x            (x) float64 4.577e+05 4.577e+05 4.577e+05 ... 4.61e+05 4.61e+05
  * y            (y) float64 4.426e+06 4.426e+06 ... 4.425e+06 4.425e+06
    band         int64 1
    spatial_ref  int64 0
NBR images calculated from Landsat for pre- and post-Cold Springs fire.
NBR images calculated from Landsat for pre- and post-Cold Springs fire.
<xarray.DataArray (y: 3, x: 8)>
array([[3., 2., 2., 4., 4., 3., 3., 3.],
       [3., 4., 4., 4., 3., 4., 3., 3.],
       [4., 4., 4., 4., 4., 4., 3., 3.]])
Coordinates:
  * y            (y) float64 4.446e+06 4.446e+06 4.445e+06
  * x            (x) float64 -8.988e+06 -8.988e+06 ... -8.986e+06 -8.985e+06
    band         int64 1
    spatial_ref  int64 0
NBR images calculated from MODIS for pre- and post-Cold Springs fire.
NBR images calculated from MODIS for pre- and post-Cold Springs fire.
Classified difference NBR (dNBR) images calculated from Landsat and MODIS for the Cold Springs fire.
Classified difference NBR (dNBR) images calculated from Landsat and MODIS for the Cold Springs fire.

Landsat vs MODIS Burned Area

Burned Landsat class 4:
Burned Landsat class 5:
Burned MODIS class 4:
Burned MODIS class 5:

BONUS PLOT - Difference NDVI (dNDVI) Using Landsat & MODIS

Landsat and MODIS NDVI Normalized Difference from before and after the Cold Springs fire.
Landsat and MODIS NDVI Normalized Difference from before and after the Cold Springs fire.