differences for PR #494

Author: actions-user

04-raster-calculations-in-r.md

@@ -33,6 +33,35 @@ See the [lesson homepage](.) for detailed information about the software,
 data, and other prerequisites you will need to work through the examples in 
 this episode.
 
+### Load the Data
+
+For this episode, we will use the DTM and DSM from the NEON Harvard Forest 
+Field site and San Joaquin Experimental Range. If you don't still have
+them loaded, do so now and turn them into dataframes:
+
+# DSM (tree top) data for Harvard Forest
+DSM_HARV <- 
+  rast("data/NEON-DS-Airborne-Remote-Sensing/HARV/DSM/HARV_dsmCrop.tif")
+
+DSM_HARV_df <- as.data.frame(DSM_HARV, xy = TRUE)
+
+# DTM (bare earth) data for Harvard Forest
+DTM_HARV <- 
+  rast("data/NEON-DS-Airborne-Remote-Sensing/HARV/DTM/HARV_dtmCrop.tif")
+
+DTM_HARV_df <- as.data.frame(DTM_HARV, xy = TRUE)
+
+# DSM data for SJER
+DSM_SJER <- 
+  rast("data/NEON-DS-Airborne-Remote-Sensing/SJER/DSM/SJER_dsmCrop.tif")
+
+DSM_SJER_df <- as.data.frame(DSM_SJER, xy = TRUE)
+
+# DTM data for SJER
+DTM_SJER <- 
+  rast("data/NEON-DS-Airborne-Remote-Sensing/SJER/DTM/SJER_dtmCrop.tif")
+
+DTM_SJER_df <- as.data.frame(DTM_SJER, xy = TRUE)
 
 ::::::::::::::::::::::::::::::::::::::::::::::::::
 
@@ -64,11 +93,7 @@ with the influence of ground elevation removed.
 
 ::::::::::::::::::::::::::::::::::::::::::::::::::
 
-### Load the Data
 
-For this episode, we will use the DTM and DSM from the NEON Harvard Forest 
-Field site and San Joaquin Experimental Range, which we already have loaded 
-from previous episodes.
 
 :::::::::::::::::::::::::::::::::::::::  challenge
 

05-raster-multi-band-in-r.md

@@ -1,7 +1,7 @@
 ---
 title: Work with Multi-Band Rasters
-teaching: 40
-exercises: 20
+teaching: 30
+exercises: 15
 source: Rmd
 ---
 
@@ -76,6 +76,8 @@ ggplot() +
 
 <img src="fig/05-raster-multi-band-in-r-rendered-harv-rgb-band1-1.png" style="display: block; margin: auto;" />
 
+To import the green band, we would use `lyrs = 2`.
+
 :::::::::::::::::::::::::::::::::::::::  challenge
 
 ## Challenge
@@ -136,56 +138,8 @@ represent pixels with less red in them (less red was reflected). To plot an RGB
 image, we mix red + green + blue values into one single color to create a full 
 color image - similar to the color image a digital camera creates.
 
-### Import A Specific Band
-
-We can use the `rast()` function to import specific bands in our raster object
-by specifying which band we want with `lyrs = N` (N represents the band number we
-want to work with). To import the green band, we would use `lyrs = 2`.
-
-
-``` r
-RGB_band2_HARV <-  
-  rast("data/NEON-DS-Airborne-Remote-Sensing/HARV/RGB_Imagery/HARV_RGB_Ortho.tif", 
-       lyrs = 2)
-```
-
-We can convert this data to a data frame and plot the same way we plotted the red band:
-
+ 
 
-``` r
-RGB_band2_HARV_df <- as.data.frame(RGB_band2_HARV, xy = TRUE)
-```
-
-
-``` r
-ggplot() +
-  geom_raster(data = RGB_band2_HARV_df,
-              aes(x = x, y = y, alpha = HARV_RGB_Ortho_2)) + 
-  coord_equal()
-```
-
-<img src="fig/05-raster-multi-band-in-r-rendered-rgb-harv-band2-1.png" style="display: block; margin: auto;" />
-
-:::::::::::::::::::::::::::::::::::::::  challenge
-
-## Challenge: Making Sense of Single Band Images
-
-Compare the plots of band 1 (red) and band 2 (green). Is the forested area 
-darker or lighter in band 2 (the green band) compared to band 1 (the red band)?
-
-:::::::::::::::  solution
-
-## Solution
-
-We'd expect a *brighter* value for the forest in band 2 (green) than in band 1 
-(red) because the leaves on trees of most often appear "green" - healthy leaves 
-reflect MORE green light than red light.
-
-
-
-:::::::::::::::::::::::::
-
-::::::::::::::::::::::::::::::::::::::::::::::::::
 
 ## Raster Stacks in R
 

data/.gitignore

@@ -0,0 +1 @@
+4.1.0

data/Global/Boundaries/ne_110m_graticules_all/ne_110m_graticules_1.README.html

@@ -0,0 +1 @@
+UTF-8

data/Global/Boundaries/ne_110m_graticules_all/ne_110m_graticules_1.VERSION.txt

@@ -0,0 +1 @@
+GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.017453292519943295]]