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

data linked to google drive

No related merge requests found
Showing with 12 additions and 8 deletions
+12 -8
......@@ -44,6 +44,8 @@ For example purpose, two adjacent ALS tiles are provided, in two versions:
* one with altitude values, in the folder "../data/aba.model/ALS/tiles.laz", these will be used for terrain statistics computations,
* one with (normalized) height values, in the folder "../data/aba.model/ALS/tiles.norm.laz", these will be used for tree detection and point metrics computation.
Files are too large to be hosted on the gitlab repository, they can be downloaded as a [zip file](https://drive.google.com/u/0/uc?export=download&id=1ripo-PLZ8_IjE7rAQ2RECj-fjg1UpC5i) from Google drive, and should be extracted in the folder "data/aba.model/ALS/" before proceeding with the processing.
```{r loadALS, include = TRUE}
cata.height <- lidR::catalog("../data/aba.model/ALS/tiles.norm.laz/")
cata.altitude <- lidR::catalog("../data/aba.model/ALS/tiles.laz/")
......@@ -388,11 +390,12 @@ par(mfrow = c(2, 2))
# plot raster and vector
raster::plot(metrics.map$stratum, main = "Ownership", legend=FALSE, col=c("green", "blue"), xaxt='n', yaxt='n')
legend("bottom", legend=c("private","public"), fill=c("green","blue"))
limits <- range(c(raster::values(prediction.map.private), raster::values(prediction.map.public)))
plot(public.cropped, col = NA, add = TRUE)
raster::plot(prediction.map.mixed, main = "Stratified model", zlim = c(0, 70), xaxt='n', yaxt='n')
raster::plot(prediction.map.mixed, main = "Stratified model", zlim = limits, xaxt='n', yaxt='n')
plot(public.cropped, col = NA, add = TRUE)
raster::plot(prediction.map.private, main = "Private model", zlim = c(0, 70), xaxt='n', yaxt='n', legend=FALSE)
raster::plot(prediction.map.public, main = "Public model", zlim = c(0, 70), xaxt='n', yaxt='n', legend=FALSE)
raster::plot(prediction.map.private, main = "Private model", zlim = limits, xaxt='n', yaxt='n', legend=FALSE)
raster::plot(prediction.map.public, main = "Public model", zlim = limits, xaxt='n', yaxt='n', legend=FALSE)
```
### Forest mask and thresholds
......
No preview for this file type
No preview for this file type
This source diff could not be displayed because it is too large. You can view the blob instead.
No preview for this file type
Supports Markdown
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