Commit 05eb8e7c authored by Monnet Jean-Matthieu's avatar Monnet Jean-Matthieu
Browse files

Add licence of public forest GIS layer

parent fc3e659c
...@@ -3,8 +3,8 @@ title: "R workflow for ABA data preparation" ...@@ -3,8 +3,8 @@ title: "R workflow for ABA data preparation"
author: "Jean-Matthieu Monnet" author: "Jean-Matthieu Monnet"
date: "`r Sys.Date()`" date: "`r Sys.Date()`"
output: output:
pdf_document: default
html_document: default html_document: default
pdf_document: default
papersize: a4 papersize: a4
bibliography: "./bib/bibliography.bib" bibliography: "./bib/bibliography.bib"
--- ---
...@@ -491,7 +491,9 @@ hist(plots$D.mean.cm, main="Mean diameter", xlab="(cm)", ylab="Plot number") ...@@ -491,7 +491,9 @@ hist(plots$D.mean.cm, main="Mean diameter", xlab="(cm)", ylab="Plot number")
In this area, public forests are generally managed in a different way compared to private forests, resulting in different forest structures. Ownership is also linked to species composition. This classification could be used as a stratification when calibrating relationships between ALS metrics and forest parameters. In this area, public forests are generally managed in a different way compared to private forests, resulting in different forest structures. Ownership is also linked to species composition. This classification could be used as a stratification when calibrating relationships between ALS metrics and forest parameters.
Plots will be attributed to strata based on external GIS layers. The first layer is a vector map of public forests. The following map displays the plots colored based on ownership, after spatial extraction of this information from the polygons of public forests (black borders). Plots will be attributed to strata based on external GIS layers. The first layer is a vector map of public forests (Forêts publiques, ONF Paris, 2019.) which is available under the [Open Licence Version 2.0](https://www.etalab.gouv.fr/licence-ouverte-open-licence).
The following lines load the public forests layer and intersects it with the plot locations.
```{r stratumExtraction, include=TRUE, message=FALSE, warning=FALSE, fig.width = 4, fig.height = 6} ```{r stratumExtraction, include=TRUE, message=FALSE, warning=FALSE, fig.width = 4, fig.height = 6}
# load GIS layer of public forests # load GIS layer of public forests
...@@ -506,7 +508,7 @@ plots.sf$stratum <- factor(ifelse(is.na(plots.sf$FID), "private", "public")) ...@@ -506,7 +508,7 @@ plots.sf$stratum <- factor(ifelse(is.na(plots.sf$FID), "private", "public"))
plots.sf$FID <- NULL plots.sf$FID <- NULL
``` ```
The following lines maps the plots with the GIS layer defining the public stratum. The following map displays the plots colored based on ownership, after spatial extraction of this information from the polygons of public forests (black borders).
```{r stratumPlot, include=TRUE, message=FALSE, warning=FALSE, fig.width = 4, fig.height = 6} ```{r stratumPlot, include=TRUE, message=FALSE, warning=FALSE, fig.width = 4, fig.height = 6}
# project points into map system # project points into map system
......
...@@ -638,7 +638,8 @@ hist(plots$D.mean.cm, main = &quot;Mean diameter&quot;, xlab = &quot;(cm)&quot;, ...@@ -638,7 +638,8 @@ hist(plots$D.mean.cm, main = &quot;Mean diameter&quot;, xlab = &quot;(cm)&quot;,
<div id="add-stratum-information-from-external-data" class="section level3"> <div id="add-stratum-information-from-external-data" class="section level3">
<h3>Add stratum information from external data</h3> <h3>Add stratum information from external data</h3>
<p>In this area, public forests are generally managed in a different way compared to private forests, resulting in different forest structures. Ownership is also linked to species composition. This classification could be used as a stratification when calibrating relationships between ALS metrics and forest parameters.</p> <p>In this area, public forests are generally managed in a different way compared to private forests, resulting in different forest structures. Ownership is also linked to species composition. This classification could be used as a stratification when calibrating relationships between ALS metrics and forest parameters.</p>
<p>Plots will be attributed to strata based on external GIS layers. The first layer is a vector map of public forests. The following map displays the plots colored based on ownership, after spatial extraction of this information from the polygons of public forests (black borders).</p> <p>Plots will be attributed to strata based on external GIS layers. The first layer is a vector map of public forests (Forêts publiques, ONF Paris, 2019.) which is available under the <a href="https://www.etalab.gouv.fr/licence-ouverte-open-licence">Open Licence Version 2.0</a>.</p>
<p>The following lines load the public forests layer and intersects it with the plot locations.</p>
<pre class="r"><code># load GIS layer of public forests <pre class="r"><code># load GIS layer of public forests
public &lt;- sf::st_read(&quot;./data/aba.model/GIS/Public4Montagnes.shp&quot;, stringsAsFactors = TRUE, public &lt;- sf::st_read(&quot;./data/aba.model/GIS/Public4Montagnes.shp&quot;, stringsAsFactors = TRUE,
quiet = TRUE) quiet = TRUE)
...@@ -650,7 +651,7 @@ plots.sf &lt;- sf::st_join(plots.sf, public) ...@@ -650,7 +651,7 @@ plots.sf &lt;- sf::st_join(plots.sf, public)
plots.sf$stratum &lt;- factor(ifelse(is.na(plots.sf$FID), &quot;private&quot;, &quot;public&quot;)) plots.sf$stratum &lt;- factor(ifelse(is.na(plots.sf$FID), &quot;private&quot;, &quot;public&quot;))
# remove column # remove column
plots.sf$FID &lt;- NULL</code></pre> plots.sf$FID &lt;- NULL</code></pre>
<p>The following lines maps the plots with the GIS layer defining the public stratum.</p> <p>The following map displays the plots colored based on ownership, after spatial extraction of this information from the polygons of public forests (black borders).</p>
<pre class="r"><code># project points into map system <pre class="r"><code># project points into map system
plots.sf.transform &lt;- sf::st_transform(plots.sf, 4326) plots.sf.transform &lt;- sf::st_transform(plots.sf, 4326)
public.transform &lt;- sf::st_transform(public, 4326) public.transform &lt;- sf::st_transform(public, 4326)
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment