Zonas aptas para el culivo de café en temperatura y precipitación
list.file
stack
sum, mean
# imporar libreria
require(raster)
## Loading required package: raster
## Loading required package: sp
library(readr)
# Temprattura promedio
img_temp=list.files("D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/", full.names = TRUE)
print(img_temp) # solo es una lista
## [1] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_01.tif"
## [2] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_02.tif"
## [3] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_03.tif"
## [4] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_04.tif"
## [5] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_05.tif"
## [6] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_06.tif"
## [7] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_07.tif"
## [8] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_08.tif"
## [9] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_09.tif"
## [10] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_10.tif"
## [11] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_11.tif"
## [12] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_temperatura/wc2.1_10m_tavg_12.tif"
temp_anual=stack(img_temp) # crea varios raster
names(temp_anual)=month.name # poner nombre mes de
plot(temp_anual)
promedio_temp = mean(temp_anual) # saca raster de la temp promedio
plot(promedio_temp)
# Precipitacion anual
img_precip=list.files("D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/", full.names = TRUE)
print(img_precip)
## [1] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_01.tif"
## [2] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_02.tif"
## [3] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_03.tif"
## [4] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_04.tif"
## [5] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_05.tif"
## [6] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_06.tif"
## [7] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_07.tif"
## [8] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_08.tif"
## [9] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_09.tif"
## [10] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_10.tif"
## [11] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_11.tif"
## [12] "D:/ESPECIALIZACION/SEMESTRE 1/1. Tratamiento de datos/Clase 5 Raster/Img_precipitacion/wc2.1_10m_prec_12.tif"
precipitacion=stack(img_precip)
precip_anual=sum(precipitacion)
plot(precip_anual)
# Consulta: muestre de temp entre 17°y 23°C
temp_cafe=promedio_temp>17&promedio_temp<23
plot(temp_cafe)
# Consulta: muestre las precipitaciones entre 1800 y 2000 mm
precip_cafe=precip_anual>1800&precip_anual<2000
plot(precip_cafe)
# Consulta: muestre zonas aptas cultivo cafe en temperatura y precipitacion
zonas_op_cafe=precip_cafe + temp_cafe
plot(zonas_op_cafe, main="zonas aptas de temp y precipt para el cultivo de café")
