Homework Submission

Produce a Report

Create a new R markdown document. Name it: lastName-firstInitial-weeknumber.Rmd Within your .Rmd document, include the plots listed below. When you are done with your report, use knitr to convert it to html format. Submit both the .Rmd file and the .html file. Be sure to name your files as instructed above.

Use knitr Code Chunk Arguments

In your final report, use the following knitr code chunk arguments to hide messages and warnings and code as you see fit.

• message = FALSE, warning = FALSE Hide warnings and messages in a code chunk
• echo = FALSE Hide code and just show code output
• results = 'hide' Hide the verbose output from some functions like readOGR().

Answer the Questions Below in Your Report

1. What is the spatial resolution of NAIP, Landsat & MODIS data in meters?
   • 1b. How can differences in spatial resolution in the data that you are using impact analysis results?
2. Calculate the area of “high severity” and the area of “moderate severity” burn in meters using the dNBR for Landsat and MODIS respectively. State what the area in meters is for each data type (Landsat and MODIS) and each level of fire severity (high and moderate) in your answer.
   • 2b. Is the total area as derived from MODIS different from the area derived from Landsat? Why / why not? Use plots 4 and 5 to discuss any differences that you notice.
3. Write 1-3 paragraphs that describes the Cold Springs fire. Include:
   • 3b. Where and when the fire occurred.
   • 3c. What started the fire.
   • 3d. A brief discussion of how fire impacts the natural vegetation structure of a forest and also humans (i.e. are fires always bad?).
4. Describe how dNBR and NDVI can be used to study the impacts of a fire. In your answer, be sure to discuss which parts of the spectrum each index uses and why those wavelength values are used. (i.e. what about NIR light makes it useful for NDVI and what about SWIR light make it useful for a burn index).

Refer to your plots in your answer.

For all of your answers:

• Be sure to carefully proofread your report before handing it in.
• Be sure to cite atleast 2 of the assigned articles in your answers.
• Be sure to use proper citation format.

Include the Plots Below.

For all plots:

• Be sure to describe what each plot shows in your final report using a figure caption argument in your code chunks: fig.cap = "caption here.
• Add appropriate titles that tell someone reading your report what the map shows.
• Use clear legends as appropriate - especially for the classified data plots!
• Only include code that is directly related to creating your plot. Do not include intermediate plot outputs in your report - only the final plot (unless it helps you directly answer a question asked above).
• Use the overlay() function and your normalized_diff() function to computer both NDVI and NBR.

Plot 1 - Grid of Plots: NAIP, Landsat and MODIS

Use the plotRGB() function to create color infrared (also called false color) images using:

• NAIP data
• Landsat data
• MODIS data

in one figure collected pre-fire.

For each map be sure to:

• Overlay the fire boundary layer (vector_layers/fire-boundary-geomac/co_cold_springs_20160711_2200_dd83.shp).

• Use the band combination r = infrared band, g = red band, b = green band. You can use mfrow=c(rows, columns).

• Render the map to the extent of the fire boundary layer using either the ext = extent() plot argument. Or crop the data to the fire boundary extent.
• Be sure to label each plot with the data type (NAIP vs. Landsat vs. MODIS) and spatial resolution.

Use this figure to help answer question 1 above. An example of what this plot should look like (without all of the labels that you need to add), is here at the bottom of the page.

Plot 2 - Difference NBR (dNBR) Using Landsat Data

Create a map of the classified dNBR using Landsat data collected before and after the Cold Springs fire. Overlay the fire extent vector_layers/fire-boundary-geomac/co_cold_springs_20160711_2200_dd83.shp on top of the NBR map. Add a legend that clearly describes each burn “class” that you applied using the table provided in the lessons.

Be sure to use cloud free data.

Plot 3 - Difference NBR (dNBR) Using MODIS Data

Create a map of the classified dNBR using MODIS data collected before and after the Cold Springs fire. Overlay the fire extent vector_layers/fire-boundary-geomac/co_cold_springs_20160711_2200_dd83.shp on top of the NBR map. Add a legend that clearly describes each burn “class” that you applied using the table provided in the lessons.

Be sure to use cloud free data.

dNBR burn classes

Note that depending on how you scaled your data, you may need to scale the values below by a factor of 10.

Plot 4 - Difference NDVI Landsat Data

Create a map of the the difference in NDVI - pre vs post fire. Classify the change in NDVI as you see fit. Overlay the fire extent vector_layers/fire-boundary-geomac/co_cold_springs_20160711_2200_dd83.shp on top of the difference NDVI map. Add a legend that clearly describes each NDVI difference “class”.

Be sure to use cloud free NDVI data.

Plot 5 - Difference NDVI MODIS Data

Create a map of the the difference in NDVI - pre vs post fire. Classify the change in NDVI as you see fit. Overlay the fire extent vector_layers/fire-boundary-geomac/co_cold_springs_20160711_2200_dd83.shp on top of the difference NDVI map. Add a legend that clearly describes each NDVI difference “class”.

Be sure to use cloud free NDVI data.

Bonus Opportunity (2 points)

Find & download a cloud-free Landsat scene from the summer 2017 using Earth Explorer:

• Calculate dNBR using that scene compared to the pre-fire scene that you downloaded from Earth Explorer.
• Plot the data using the dNBR classification table.
• Calculate the total area of high vs. moderate severity burn for the fire scar.

Address the question: Has the fire scar area regenerated at all?

Homework Due: Monday November 6, 2017 @ 8AM.

Submit your report in both .Rmd and .html format to the D2l dropbox.