Inundacion.R

library(mapview)
library(raster)

## Loading required package: sp

pal = mapviewPalette("mapviewTopoColors")

Profundidad <- raster("depth.tif", package = "mapview")
Prof <- stack(Profundidad)

mapview(Prof[[1]], col.regions = pal(100), at = seq(-0.2, 1, 0.2), legend = TRUE)

#la distribución de valores en el ráster
hist(Prof, main="profundidades de inundacion",
     col= "purple",
     maxpixels=22000000)

#Visualizar los datos

plot(Prof,
     main="Profundidad")

summary(Prof)

## Warning in .local(object, ...): summary is an estimate based on a sample of 1e+05 cells (62.81% of all cells)

##                Depth
## Min.    1.029968e-03
## 1st Qu. 9.464645e-01
## Median  1.673088e+00
## 3rd Qu. 2.313713e+00
## Max.    3.899666e+00
## NA's    5.984500e+04

#para filtrar los valores
image(Prof, zlim=c(0.5,0.9))

#Grafico de caja y bigote:
boxplot(Prof)

#Erosion

pal = mapviewPalette("mapviewTopoColors")

Erosion <- raster("Shear.tif", package = "mapview")
Erosion<- stack(Erosion)

mapview(Erosion[[1]], col.regions = pal(100), at = seq(-0.2, 1, 0.2), legend = TRUE)

# la distribución de valores en el ráster
hist(Erosion, main="profundidades de inundacion",
     col= "purple",
     maxpixels=22000000)

#Visualizar los datos

plot(Erosion,
     main="Erosion")

summary(Erosion)

## Warning in .local(object, ...): summary is an estimate based on a sample of 1e+05 cells (62.81% of all cells)

##                Shear
## Min.    6.090858e-06
## 1st Qu. 2.913535e-03
## Median  1.112604e-02
## 3rd Qu. 3.845156e-02
## Max.    1.776714e+01
## NA's    5.984500e+04

#para filtrar los valores
image(Erosion, zlim=c(0.5,0.9))

#Grafico de caja y bigote:
boxplot(Erosion)