#LIBRERIAS
require(rgdal)
require(raster)
library(doParallel)
provincia <- "Albacete"
#Capa de provincias.
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/2_Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
#Capa de provincias.
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/02-Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
#LIBRERIAS
require(rgdal)
require(raster)
library(doParallel)
#################################################################################################################
UseCores <- detectCores() -1
#Register CoreCluster
cl<- makeCluster(UseCores)
start_time <- Sys.time()
#VARIABLES
diccionario_provincias <- list(
"Albacete" = list(tiles = c("30SXH")),
"Alicante" = list(tiles = c("30SYJ")),
"Almeria" = list(tiles = c("30SWG","30SXG")),
"Badajoz" = list(tiles = c("29SPD","29SQD","30STH","30STJ")),
"Baleares" = list(tiles = c("31TDE","31TEE")),
"Caceres" = list(tiles = c("29TPE","30STJ","30TTK")),
"Cadiz" = list(tiles = c("30STF")),
"Castellon" = list(tiles = c("30TYK")),
"Cordoba" = list(tiles = c("30SUG")),
"Girona" = list(tiles = c("31TDG","31TEG")),
"Granada" = list(tiles = c("30SVG","30SWG")),
"Huelva" = list(tiles = c("29SPB","29SPC","29SQB","30STG")),
"Huesca" = list(tiles = c("30TXM","30TYM")),
"Jaen" = list(tiles = c("30SVH")),
"Lleida" = list(tiles = c("30TYM")),
"Malaga" = list(tiles = c("30SUF")),
"Murcia" = list(tiles = c("30SXH")),
"Navarra" = list(tiles = c("30TWM","30TXM")),
"Sevilla" = list(tiles = c("30STG")),
"Tarragona" = list(tiles = c("30TYL","31TCF")),
"Teruel" = list(tiles = c("30TYL")),
"Valencia" = list(tiles = c("30SYJ")),
"Zaragoza" = list(tiles = c("30TXM"))
)
provincias <- c("Albacete")
provincias <- c("Valencia")
for(provincia in provincias){
print(provincia)
teselas <- paste0(diccionario_provincias[[provincia]]$tiles)
proj <- "+init=EPSG:25830"
descarga <- "Arroz"
#descarga <- "Arroz"
anios <- c("2022","2023","2024")
anios <- "2021"
#Capa de provincias.
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/02-Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
#################################################################################################################
#RUTAS
path <- paste0("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/04-Descargas/",descarga)
dir_salida <- paste0(path,"/provincias/")
#################################################################################################################
x <- list.files(path, pattern = ".tif", full.names = TRUE)
#aqui metemos las imgs para mosaicarlas.
for(anio in anios){
lista_inun = list()
lista_culti = list()
lista_cos = list()
for(tesela in teselas){
for(i in x){
if(grepl(paste0(tesela,"_inundado_",anio),i)){
lista_inun = append(lista_inun,i)
}
elif(grepl(paste0(tesela,"_cultivado_",anio),i)){
lista_culti = append(lista_culti,i)
}
elif(grepl(paste0(tesela,"_cosechado_",anio),i)){
lista_cos = append(lista_cos,i)
}
}
}
if(length(lista_inun>1)){
#Mosaicamos
rasters_inun <- lapply(lista_inun,raster)
rasters_culti <- lapply(lista_culti,raster)
rasters_cos <- lapply(lista_cos,raster)
r_mosaic_inun <- do.call(mosaic, c(rasters_inun, fun='max'))
r_mosaic_culti <- do.call(mosaic, c(rasters_culti, fun='max'))
r_mosaic_cos <- do.call(mosaic, c(rasters_cos, fun='max'))
#Recortamos
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}else{
r_mosaic_inun <- raster(lista_inun[1])
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
r_mosaic_culti <- raster(lista_culti[1])
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
r_mosaic_cos <- raster(lista_cos[1])
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
}
}
end_time = Sys.time()
print(end_time - start_time)
}
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
#LIBRERIAS
require(rgdal)
require(raster)
library(doParallel)
#################################################################################################################
UseCores <- detectCores() -1
#Register CoreCluster
cl<- makeCluster(UseCores)
start_time <- Sys.time()
#VARIABLES
diccionario_provincias <- list(
"Albacete" = list(tiles = c("30SXH")),
"Alicante" = list(tiles = c("30SYJ")),
"Almeria" = list(tiles = c("30SWG","30SXG")),
"Badajoz" = list(tiles = c("29SPD","29SQD","30STH","30STJ")),
"Baleares" = list(tiles = c("31TDE","31TEE")),
"Caceres" = list(tiles = c("29TPE","30STJ","30TTK")),
"Cadiz" = list(tiles = c("30STF")),
"Castellon" = list(tiles = c("30TYK")),
"Cordoba" = list(tiles = c("30SUG")),
"Girona" = list(tiles = c("31TDG","31TEG")),
"Granada" = list(tiles = c("30SVG","30SWG")),
"Huelva" = list(tiles = c("29SPB","29SPC","29SQB","30STG")),
"Huesca" = list(tiles = c("30TXM","30TYM")),
"Jaen" = list(tiles = c("30SVH")),
"Lleida" = list(tiles = c("30TYM")),
"Malaga" = list(tiles = c("30SUF")),
"Murcia" = list(tiles = c("30SXH")),
"Navarra" = list(tiles = c("30TWM","30TXM")),
"Sevilla" = list(tiles = c("30STG")),
"Tarragona" = list(tiles = c("30TYL","31TCF")),
"Teruel" = list(tiles = c("30TYL")),
"Valencia" = list(tiles = c("30SYJ")),
"Zaragoza" = list(tiles = c("30TXM"))
)
provincias <- c("Valencia")
provincia <- "Valencia"
print(provincia)
teselas <- paste0(diccionario_provincias[[provincia]]$tiles)
proj <- "+init=EPSG:25830"
descarga <- "Arroz"
#descarga <- "Arroz"
anios <- c("2022","2023","2024")
anios <- "2021"
#Capa de provincias.
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/02-Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
#RUTAS
path <- paste0("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/04-Descargas/",descarga)
dir_salida <- paste0(path,"/provincias/")
x <- list.files(path, pattern = ".tif", full.names = TRUE)
x
anio <- "2021"
lista_inun = list()
lista_culti = list()
lista_cos = list()
for(tesela in teselas){
for(i in x){
ifelse(grepl(paste0(tesela,"_inundado_",anio),i),lista_inun= append(lista_inun,i),
ifelse(grepl(paste0(tesela,"_cultivado_",anio),i),lista_culti = append(lista_culti,i),lista_cos = append(lista_cos,i)))
}
if(length(lista_inun>1)){
#Mosaicamos
rasters_inun <- lapply(lista_inun,raster)
rasters_culti <- lapply(lista_culti,raster)
rasters_cos <- lapply(lista_cos,raster)
r_mosaic_inun <- do.call(mosaic, c(rasters_inun, fun='max'))
r_mosaic_culti <- do.call(mosaic, c(rasters_culti, fun='max'))
r_mosaic_cos <- do.call(mosaic, c(rasters_cos, fun='max'))
#Recortamos
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}else{
r_mosaic_inun <- raster(lista_inun[1])
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
r_mosaic_culti <- raster(lista_culti[1])
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
r_mosaic_cos <- raster(lista_cos[1])
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
}
end_time = Sys.time()
print(end_time - start_time)
}
i
tesela
grepl(paste0(tesela,"_inundado_",anio),i)
grepl(paste0(tesela,"_cosechado_",anio),i)
paste0(tesela,"_cosechado_",anio)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
x
i
grepl(paste0(tesela,"_inundado_",anio),i)
grepl(paste0(tesela,"_cultivado_",anio),i)
for(tesela in teselas){
for(i in x){
ifelse(grepl(paste0(tesela,"_inundado_",anio),i),lista_inun= append(lista_inun,i),
ifelse(grepl(paste0(tesela,"_cultivado_",anio),i),lista_culti = append(lista_culti,i),lista_cos = append(lista_cos,i)))
}
if(length(lista_inun>1)){
#Mosaicamos
rasters_inun <- lapply(lista_inun,raster)
rasters_culti <- lapply(lista_culti,raster)
rasters_cos <- lapply(lista_cos,raster)
r_mosaic_inun <- do.call(mosaic, c(rasters_inun, fun='max'))
r_mosaic_culti <- do.call(mosaic, c(rasters_culti, fun='max'))
r_mosaic_cos <- do.call(mosaic, c(rasters_cos, fun='max'))
#Recortamos
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}else{
r_mosaic_inun <- raster(lista_inun[1])
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
r_mosaic_culti <- raster(lista_culti[1])
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
r_mosaic_cos <- raster(lista_cos[1])
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
}
end_time = Sys.time()
print(end_time - start_time)
}
lista_inun <- list()
lista_culti <- list()
lista_cos <- list()
for(tesela in teselas){
for(i in x){
ifelse(grepl(paste0(tesela,"_inundado_",anio),i),lista_inun <- append(lista_inun,i),
ifelse(grepl(paste0(tesela,"_cultivado_",anio),i),lista_culti <- append(lista_culti,i),lista_cos <- append(lista_cos,i)))
}
if(length(lista_inun>1)){
#Mosaicamos
rasters_inun <- lapply(lista_inun,raster)
rasters_culti <- lapply(lista_culti,raster)
rasters_cos <- lapply(lista_cos,raster)
r_mosaic_inun <- do.call(mosaic, c(rasters_inun, fun='max'))
r_mosaic_culti <- do.call(mosaic, c(rasters_culti, fun='max'))
r_mosaic_cos <- do.call(mosaic, c(rasters_cos, fun='max'))
#Recortamos
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}else{
r_mosaic_inun <- raster(lista_inun[1])
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
r_mosaic_culti <- raster(lista_culti[1])
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
r_mosaic_cos <- raster(lista_cos[1])
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"cosechado.tif"),format='GTiff')
}
end_time = Sys.time()
print(end_time - start_time)
}
lista_cos[1]
raster(lista_cos[1])
raster(lista_cos[1][1])
lista_cos[1][1]
lista_cos
lista_cos{1}
lista_cos[[1]]
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
length(lista_inun>1)
length(lista_inun)>1
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
#LIBRERIAS
require(rgdal)
require(raster)
library(doParallel)
#################################################################################################################
UseCores <- detectCores() -1
#Register CoreCluster
cl<- makeCluster(UseCores)
start_time <- Sys.time()
#VARIABLES
diccionario_provincias <- list(
"Albacete" = list(tiles = c("30SXH")),
"Alicante" = list(tiles = c("30SYJ")),
"Almeria" = list(tiles = c("30SWG","30SXG")),
"Badajoz" = list(tiles = c("29SPD","29SQD","30STH","30STJ")),
"Baleares" = list(tiles = c("31TDE","31TEE")),
"Caceres" = list(tiles = c("29TPE","30STJ","30TTK")),
"Cadiz" = list(tiles = c("30STF")),
"Castellon" = list(tiles = c("30TYK")),
"Cordoba" = list(tiles = c("30SUG")),
"Girona" = list(tiles = c("31TDG","31TEG")),
"Granada" = list(tiles = c("30SVG","30SWG")),
"Huelva" = list(tiles = c("29SPB","29SPC","29SQB","30STG")),
"Huesca" = list(tiles = c("30TXM","30TYM")),
"Jaen" = list(tiles = c("30SVH")),
"Lleida" = list(tiles = c("30TYM")),
"Malaga" = list(tiles = c("30SUF")),
"Murcia" = list(tiles = c("30SXH")),
"Navarra" = list(tiles = c("30TWM","30TXM")),
"Sevilla" = list(tiles = c("30STG")),
"Tarragona" = list(tiles = c("30TYL","31TCF")),
"Teruel" = list(tiles = c("30TYL")),
"Valencia" = list(tiles = c("30SYJ")),
"Zaragoza" = list(tiles = c("30TXM"))
)
provincias <- c("Albacete","Alicante","Almeria","Badajoz","Baleares","Caceres","Cadiz","Castellon","Cordoba","Girona","Granada","Huelva","Huesca","Jaen","Lleida","Malaga","Murcia","Navarra","Sevilla","Tarragona","Teruel","Valencia","Zaragoza")
provincia <- "Albacete"
print(provincia)
teselas <- paste0(diccionario_provincias[[provincia]]$tiles)
proj <- "+init=EPSG:25830"
descarga <- "Arroz"
#descarga <- "Arroz"
anios <- c("2022","2023","2024")
#Capa de provincias.
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/02-Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
#RUTAS
path <- paste0("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/04-Descargas/",descarga)
dir_salida <- paste0(path,"/provincias/")
x <- list.files(path, pattern = ".tif", full.names = TRUE)
x
anio <- "2022"
lista_inun <- list()
lista_culti <- list()
lista_cos <- list()
for(tesela in teselas){
for(i in x){
ifelse(grepl(paste0(tesela,"_inundado_",anio),i),lista_inun <- append(lista_inun,i),
ifelse(grepl(paste0(tesela,"_cultivado_",anio),i),lista_culti <- append(lista_culti,i),lista_cos <- append(lista_cos,i)))
}
if(length(lista_inun)>1){
#Mosaicamos
rasters_inun <- lapply(lista_inun,raster)
rasters_culti <- lapply(lista_culti,raster)
rasters_cos <- lapply(lista_cos,raster)
r_mosaic_inun <- do.call(mosaic, c(rasters_inun, fun='max'))
r_mosaic_culti <- do.call(mosaic, c(rasters_culti, fun='max'))
r_mosaic_cos <- do.call(mosaic, c(rasters_cos, fun='max'))
#Recortamos
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}else{
r_mosaic_inun <- raster(lista_inun[[1]])
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
r_mosaic_culti <- raster(lista_culti[[1]])
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
r_mosaic_cos <- raster(lista_cos[[1]])
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cosechado.tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cosechado.tif"),format='GTiff')
}
end_time = Sys.time()
print(end_time - start_time)
}
i
lista_cos
i
grepl(paste0(tesela,"_inundado_",anio),i)
tesela
teselas
anio
lista_inun <- list()
lista_culti <- list()
lista_cos <- list()
for(tesela in teselas){
for(i in x){
ifelse(grepl(paste0(tesela,"_inundado_",anio),i),lista_inun <- append(lista_inun,i),
ifelse(grepl(paste0(tesela,"_cultivado_",anio),i),lista_culti <- append(lista_culti,i),
ifelse(grepl(paste0(tesela,"_cosechado_",anio),i),lista_cos <- append(lista_cos,i),next)))
}
if(length(lista_inun)>1){
#Mosaicamos
rasters_inun <- lapply(lista_inun,raster)
rasters_culti <- lapply(lista_culti,raster)
rasters_cos <- lapply(lista_cos,raster)
r_mosaic_inun <- do.call(mosaic, c(rasters_inun, fun='max'))
r_mosaic_culti <- do.call(mosaic, c(rasters_culti, fun='max'))
r_mosaic_cos <- do.call(mosaic, c(rasters_cos, fun='max'))
#Recortamos
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}else{
r_mosaic_inun <- raster(lista_inun[[1]])
capa_inun <- mask(r_mosaic_inun, capa_prov)
capa_clip_inun <- crop(capa_inun,capa_prov)
r_mosaic_culti <- raster(lista_culti[[1]])
capa_culti <- mask(r_mosaic_culti, capa_prov)
capa_clip_culti <- crop(capa_culti,capa_prov)
r_mosaic_cos <- raster(lista_cos[[1]])
capa_cos <- mask(r_mosaic_cos, capa_prov)
capa_clip_cos <- crop(capa_cos,capa_prov)
}
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cosechado.tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip_inun, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_inundado.tif"),format='GTiff')
writeRaster(capa_clip_culti, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cultivado.tif"),format='GTiff')
writeRaster(capa_clip_cos, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,"_cosechado.tif"),format='GTiff')
}
end_time = Sys.time()
print(end_time - start_time)
}
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/LULUCF24/03-Scripts/0_Auxiliares/Mosaico_Recorte_Provincia_arroz.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2024/3_Scripts/1_Procesado/5_LagsSalinas.R", encoding = 'UTF-8', echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2024/3_Scripts/1_Procesado/6_HIDROPERIODO_CURSO.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2024/3_Scripts/1_Procesado/7_Grafico_Tablas_CRUCE_SUP_v2_COBNUEVA_CURSO.R", encoding = 'UTF-8', echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2024/3_Scripts/1_Procesado/8_LAGUNA_SECA_v3_2024.R", encoding = 'UTF-8', echo=TRUE)
