require(raster)
## Loading required package: raster
## Loading required package: sp
require(gdal)
## Loading required package: gdal
## Warning in library(package, lib.loc = lib.loc, character.only = TRUE,
## logical.return = TRUE, : there is no package called 'gdal'
require(sp)
##Cargar los rasters 1:12
mapas= list.files("~/Especializacion_geomatica/Tratamiento_de_datos_espaciales/clase6/", full.names = TRUE)
temp_blobal=stack(mapas)
plot(temp_blobal[[1]])
## Cargar shapefile
mundo=shapefile("~/Especializacion_geomatica/Tratamiento_de_datos_espaciales/shapes_clase6/shape_global/g2008_0.shp")
plot(mundo)
##hacer un filtro por Colombia
head(mundo@data)
## AREA PERIMETER G2008_0_ G2008_0_ID ADM0_CODE ADM0_NAME LAST_UPDAT
## 0 649.4221006 900.013634 2 1 98 Greenland 20050415
## 1 0.1422211 4.510463 5 4 98 Greenland 20050415
## 2 0.1156834 3.409226 15 14 98 Greenland 20050415
## 3 0.1173564 2.728860 16 15 98 Greenland 20050415
## 4 95.1420039 345.882703 18 17 46 Canada 20050415
## 5 0.2801303 6.484889 19 18 98 Greenland 20050415
## CONTINENT REGION STR_YEAR0 EXP_YEAR0
## 0 Americas Northern America 0 0
## 1 Americas Northern America 0 0
## 2 Americas Northern America 0 0
## 3 Americas Northern America 0 0
## 4 Americas Northern America 0 0
## 5 Americas Northern America 0 0
Colombia=which(mundo@data$ADM0_NAME=="Colombia")
Pais=mundo[Colombia,]
plot(Pais,add=TRUE)
#realizar un corte solo para temperaturas de colombia
recorte=crop(temp_blobal,Pais)
plot(recorte)
col_temp=mask(recorte,Pais)
plot(col_temp)
#ponemos los nombres de los meses
names(col_temp)=month.name
##Instalar rasterVIS para visualización de resultados y observamos para dos meses al azar
require(rasterVis)
## Loading required package: rasterVis
## Loading required package: lattice
## Loading required package: latticeExtra
levelplot(col_temp[[1]],par.settings=BuRdTheme)+layer(sp.polygons(Pais))
levelplot(col_temp[[c(1,10)]],par.settings=BuRdTheme)+layer(sp.polygons(Pais))
#Zonas de paramo según temperatura - aproximada
plot(col_temp>=2&col_temp<=14)
cond_temp=col_temp>=2&col_temp<=14
cond_temp
## class : RasterBrick
## dimensions : 100, 73, 7300, 12 (nrow, ncol, ncell, nlayers)
## resolution : 0.1666667, 0.1666667 (x, y)
## extent : -79, -66.83333, -4.166667, 12.5 (xmin, xmax, ymin, ymax)
## crs : +proj=longlat +datum=WGS84 +no_defs
## source : memory
## names : January, February, March, April, May, June, July, August, September, October, November, December
## min values : 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
## max values : 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1
#realizamos la suma para juntar los 12 resultados
sum(cond_temp)
## class : RasterLayer
## dimensions : 100, 73, 7300 (nrow, ncol, ncell)
## resolution : 0.1666667, 0.1666667 (x, y)
## extent : -79, -66.83333, -4.166667, 12.5 (xmin, xmax, ymin, ymax)
## crs : +proj=longlat +datum=WGS84 +no_defs
## source : memory
## names : layer
## values : 0, 12 (min, max)
Apti_temp=sum(cond_temp)
levelplot(Apti_temp,par.settings=RdBuTheme)+layer(sp.polygons(Pais))
##Conteo de pixeles para ver aptitud de condición
Apti_temp[]
## [1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [25] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 NA NA NA
## [49] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [73] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [97] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0
## [121] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [145] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [169] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0
## [193] 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [217] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [241] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0
## [265] 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [289] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [313] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0
## [337] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [361] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [385] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 NA
## [409] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [433] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [457] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 NA
## [481] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [505] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [529] NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0
## [553] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [577] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [601] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 NA
## [625] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [649] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [673] NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 12 12 12 12 0 0 0 0 0
## [697] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [721] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [745] NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 12 12 0 0 0 0 0
## [769] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [793] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [817] NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0
## [841] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [865] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [889] NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [913] 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [937] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [961] NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [985] 0 11 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1009] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1033] NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0
## [1057] 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1081] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1105] NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0
## [1129] 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1153] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1177] NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0
## [1201] 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1225] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1249] NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0
## [1273] 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1297] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1321] NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0
## [1345] 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1369] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1393] NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0
## [1417] 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1441] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1465] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0
## [1489] 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1513] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1537] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0
## [1561] 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1585] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1609] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0
## [1633] 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA
## [1657] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1681] NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0
## [1705] 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA
## [1729] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1753] NA NA NA NA NA NA NA NA NA NA 0 NA NA 0 0 0 0 0 0 0 0 0 0 0
## [1777] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA
## [1801] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1825] NA NA NA NA NA NA NA NA NA NA NA 0 0 NA 0 0 0 0 0 0 0 0 0 0
## [1849] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA
## [1873] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1897] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0
## [1921] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA
## [1945] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [1969] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 NA 0 0 0 0 0 0 0 0
## [1993] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 NA NA NA NA NA NA
## [2017] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2041] NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0
## [2065] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 2 8 0 0 NA NA NA NA NA
## [2089] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2113] NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0
## [2137] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 0 0 NA NA NA NA
## [2161] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2185] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0
## [2209] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 12 4 NA NA NA
## [2233] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2257] NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0
## [2281] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 12 12 0 0
## [2305] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2329] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0
## [2353] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 12 12 0
## [2377] 0 0 NA NA NA NA NA NA NA 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2401] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0
## [2425] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 12 12
## [2449] 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA
## [2473] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0
## [2497] 0 0 0 0 0 0 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 10
## [2521] 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA
## [2545] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0
## [2569] 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2593] 0 12 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA
## [2617] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0
## [2641] 0 0 0 0 0 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2665] 0 0 0 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA
## [2689] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0
## [2713] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2737] 0 12 12 1 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA
## [2761] NA NA NA NA NA NA NA NA 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2785] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2809] 0 6 12 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2833] 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA
## [2857] NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2881] 0 4 12 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2905] 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA
## [2929] NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 10
## [2953] 0 12 12 12 12 12 4 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [2977] 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA
## [3001] NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3025] 8 0 12 12 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3049] 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA
## [3073] NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 12 0 0 0 0 0 0 0
## [3097] 12 8 12 8 12 12 12 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3121] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA
## [3145] NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 12 12 0 0 0 0 0 0
## [3169] 9 12 12 12 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3193] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA
## [3217] NA NA NA NA NA NA 0 0 0 0 0 0 0 12 0 0 0 12 12 0 0 0 0 0
## [3241] 0 12 12 12 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3265] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA
## [3289] NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 12 12 0 0 0 0
## [3313] 0 4 12 12 12 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3337] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA
## [3361] NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 12 12 0 0 0
## [3385] 0 0 10 10 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3409] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA
## [3433] NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 12 12 0 0
## [3457] 0 0 0 12 12 12 12 9 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3481] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA
## [3505] NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 12 0 0 0
## [3529] 0 0 0 0 12 12 0 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3553] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA
## [3577] NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 12 5 0
## [3601] 0 0 0 0 0 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3625] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA
## [3649] NA NA NA NA NA NA NA NA NA NA NA 0 NA 0 0 0 0 0 0 0 0 12 12 0
## [3673] 0 0 0 0 0 11 12 12 6 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3697] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA
## [3721] NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 12 12 2
## [3745] 0 0 0 0 0 0 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3769] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA
## [3793] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 12 12 12
## [3817] 0 0 0 0 0 0 11 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3841] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3865] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 12 12
## [3889] 0 0 0 0 0 0 0 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3913] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3937] 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 12
## [3961] 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [3985] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4009] 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 3 0 0 0
## [4033] 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4057] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4081] 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0
## [4105] 0 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4129] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4153] 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0
## [4177] 0 12 12 12 12 0 0 0 0 12 11 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4201] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4225] 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0
## [4249] 0 0 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4273] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4297] 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0
## [4321] 0 0 12 12 0 0 0 0 0 0 4 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4345] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4369] 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0
## [4393] 0 0 0 12 12 4 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4417] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4441] 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA 0 0 0 0 0 0 0 0 0
## [4465] 0 0 0 12 12 12 0 0 0 0 12 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4489] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4513] 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA 0 0 0 0 0 0 0 0
## [4537] 0 0 0 0 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4561] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4585] 0 0 0 0 0 0 0 0 0 0 0 NA 0 NA NA NA NA NA 0 0 0 0 0 0
## [4609] 0 0 0 0 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4633] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4657] 0 0 0 0 0 0 0 NA 0 0 0 NA NA 0 0 NA 0 0 NA 0 0 0 0 0
## [4681] 0 0 0 0 0 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4705] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4729] NA 0 0 0 0 0 0 0 0 NA NA NA NA NA NA 0 NA 0 0 0 0 0 0 0
## [4753] 0 0 0 0 0 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4777] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4801] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 NA NA 0 0 0 0 0
## [4825] 0 0 0 0 0 0 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4849] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4873] 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 NA NA NA 0 0 0
## [4897] 0 0 0 0 0 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4921] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4945] 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 NA NA NA NA
## [4969] 0 0 0 12 12 9 12 12 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [4993] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5017] 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5041] NA NA 0 12 12 12 12 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5065] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5089] 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5113] NA NA NA NA NA NA 12 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5137] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5161] 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5185] NA NA NA NA NA NA NA NA 12 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5209] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5233] 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5257] NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 NA 0 0 0 0 0 0 0
## [5281] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5305] 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5329] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 NA NA 0 0 0 0 0
## [5353] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5377] 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5401] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0
## [5425] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5449] 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5473] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5497] NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5521] 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5545] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5569] NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5593] 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5617] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5641] NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5665] 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA
## [5689] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5713] NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5737] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA
## [5761] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5785] NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0
## [5809] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA NA NA
## [5833] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5857] NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0 0 0 0
## [5881] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA NA
## [5905] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [5929] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0 0 0 0
## [5953] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA NA NA
## [5977] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6001] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0 0 0
## [6025] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA NA
## [6049] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6073] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0 0
## [6097] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA
## [6121] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6145] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0 0
## [6169] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA NA
## [6193] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6217] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0 0
## [6241] 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA NA
## [6265] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6289] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6313] NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 NA
## [6337] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6361] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6385] NA NA 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
## [6409] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6433] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6457] NA NA NA 0 0 0 0 0 0 0 0 NA NA 0 0 0 NA 0 0 0 0 0 0 0
## [6481] 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6505] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6529] NA NA NA NA NA 0 0 NA 0 0 NA NA NA NA NA NA NA NA NA 0 0 0 0 0
## [6553] 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6577] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6601] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA 0
## [6625] 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6649] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6673] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6697] NA 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6721] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6745] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6769] NA 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6793] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6817] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6841] NA NA 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6865] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6889] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6913] NA NA 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6937] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6961] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [6985] NA NA 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7009] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7033] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7057] NA NA NA 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7081] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7105] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7129] NA NA NA 0 0 0 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7153] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7177] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7201] NA NA NA NA NA NA NA 0 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7225] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7249] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7273] NA NA NA NA NA NA NA NA NA 0 NA NA NA NA NA NA NA NA NA NA NA NA NA NA
## [7297] NA NA NA NA
table(Apti_temp[])
##
## 0 1 2 3 4 5 6 8 9 10 11 12
## 3143 5 4 1 7 3 2 4 5 4 4 157
#Seleccionamos desde el mes 9 al 12 que para el ejercicio pueden ser
#zonas donde se puedan encontrar páramos
table(Apti_temp[]>=9)
##
## FALSE TRUE
## 3169 170
#Se tiene como resultado 170 pixeles
#Se va a simular un aumento de 1.5 grados sobre los datos
#para ver como afecta el aumento de la temperatura las posibles
#zonas para encontrar paramos
col_temp_fut=col_temp+1.5
plot(col_temp_fut)
#realizamos la condición de temperatura
plot(col_temp_fut>=2&col_temp_fut<=14)
condi_fut=col_temp_fut>=2&col_temp_fut<=14
sum(condi_fut)
## class : RasterLayer
## dimensions : 100, 73, 7300 (nrow, ncol, ncell)
## resolution : 0.1666667, 0.1666667 (x, y)
## extent : -79, -66.83333, -4.166667, 12.5 (xmin, xmax, ymin, ymax)
## crs : +proj=longlat +datum=WGS84 +no_defs
## source : memory
## names : layer
## values : 0, 12 (min, max)
Apti_temp_fut=sum(condi_fut)
plot(Apti_temp_fut)
##Juntamos la capa del presente y futuro
Apti_all=stack(Apti_temp,Apti_temp_fut)
names(Apti_all)=c("Presente","Futuro +1.5°grados C")
levelplot(Apti_all,par.settings=RdBuTheme)+layer(sp.polygons(Pais))
#hacemos la resta de ambos mapas para juntarlos y ver el cambio
cate_apti=Apti_all>=9
Cambio_apti=cate_apti[[2]]-cate_apti[[1]]
levelplot(Cambio_apti,par.settings=RdBuTheme)+layer(sp.polygons(Pais))
table(Cambio_apti[])
##
## -1 0
## 48 3291
#se determina que el cambio se presentará en 48 pixeles//
#3221 pixeles siguen igual frente a la linea base// se podria decir
#que se presentará un cambio del 1.4% a futuro sobre areas que en
#los datos orginales existían y a futuro desaparecen.
48/3221*100
## [1] 1.49022