From b12328b8071017b61576f70ca77aefc6a2d69510 Mon Sep 17 00:00:00 2001
From: Francesco Sabatini <francesco.sabatini@idiv.de>
Date: Fri, 13 Mar 2020 14:58:55 +0100
Subject: [PATCH] Fixed 05_ExtractEnv to account for loc uncert

---
 code/05_ExtractEnvironment.Rmd | 12 ++++++++++--
 1 file changed, 10 insertions(+), 2 deletions(-)

diff --git a/code/05_ExtractEnvironment.Rmd b/code/05_ExtractEnvironment.Rmd
index b834e31..45db874 100644
--- a/code/05_ExtractEnvironment.Rmd
+++ b/code/05_ExtractEnvironment.Rmd
@@ -274,12 +274,21 @@ knitr::kable(tmp.sum,
                   full_width = F, position = "center")
 ```
 
+## 4 Elevation variability
+Reimport output (calculated in `04_buildHeader`)
+```{r, message=F}
+elevation.out <- read_csv("../_derived/elevatr/elevation.out.csv")
+```
+
 
-## 4 Create output and export
+## 5 Create output and export
 
 ```{r}
 soilclim <- header %>% 
   dplyr::select(PlotObservationID) %>% 
+  left_join(elevation.out %>% 
+              dplyr::select(PlotObservationID, Elevation_median, Elevation_q2.5, Elevation_q97.5, Elevation_DEM.res=DEM.res), 
+            by="PlotObservationID") %>% 
   left_join(header.shp@data %>% 
               dplyr::select(PlotObservationID) %>% 
               bind_cols(chelsa.out) %>%
@@ -288,7 +297,6 @@ soilclim <- header %>%
               bind_cols(isric.sd.out) %>% 
               distinct(), 
             by="PlotObservationID")
-
 ```
 
 ```{r, echo=F}
-- 
GitLab