feedback.Rmd

---
title: "Feedback raster exercise"
author: "Loïc Dutrieux"
date: "Friday, January 09, 2015"
output: slidy_presentation
---

# Yesterday

* Real world example
    * Data from USGS as they come

## Step by step 

* Identify RED, NIR, and fmask layers using `list.file()`
    * `list.files(dir, pattern = glob2rx('*band3*.tif'), full.names = TRUE)`
    * `list.files(dir, pattern = glob2rx('*fmask*.tif'), full.names = TRUE)`
* (Calculate NDVI -> apply cloud mask) * 2
* Crop to intersecting extent
    * `intersect()`
* Substract the two layers


* No need to stack


---

# Visualization


<img src="figs/intersect.png" style="width: 1000px;"/>

---

# Visualization

<img src="figs/final.png" style="width: 1000px;"/>

---

# Visualization

<img src="figs/panels.png" style="width: 1000px;"/>


Google earth

---

# My solution:

[https://github.com/GeoScripting-WUR/ExerciseLesson4](https://github.com/GeoScripting-WUR/ExerciseLesson4)