"Madrid" = list(tiles = c("30TUK","30TUL","30TVK","30TVL")),
"Malaga" = list(tiles = c("30STF","30SUF","30SUG","30SVF")),
"Murcia" = list(tiles = c("30SWG","30SWH","30SXG","30SXH")),
"Navarra" = list(tiles = c("30TWM","30TWN","30TWP","30TXM","30TXN","30TXP")),
"Ourense" = list(tiles = c("29TNG","29TNH","29TPG","29TPH")),
"Palencia" = list(tiles = c("30TUM","30TUN","30TVM","30TVN")),
"Pontevedra" = list(tiles = c("29TMG","29TMH","29TNG","29TNH")),
"Salamanca" = list(tiles = c("29TPE","29TPF","29TQE","30TTK","30TTL","30TUL")),
"Segovia" = list(tiles = c("30TUL","30TVL","30TVM")),
"Sevilla" = list(tiles = c("29SQA","29SQB","30STF","30STG","30STH","30SUF","30SUG")),
"Soria" = list(tiles = c("30TVL","30TVM","30TWL","30TWM")),
"Tarragona" = list(tiles = c("30TYK","30TYL","31TCE","31TCF","31TCG")),
"Teruel" = list(tiles = c("30TWK","30TXK","30TXL","30TYK","30TYL")),
"Toledo" = list(tiles = c("30SUJ","30SVJ","30SWJ","30TTK","30TUK","30TVK")),
"Valencia" = list(tiles = c("30SXH","30SXJ","30SYH","30SYJ","30TXK","30TYK")),
"Valladolid" = list(tiles = c("30TTM","30TUL","30TUM","30TVL","30TVM")),
"Bizkaia" = list(tiles = c("30TVN","30TVP","30TWN","30TWP")),
"Zamora" = list(tiles = c("29TPG","30TTL","30TTM","30TUL","30TUM")),
"Zaragoza" = list(tiles = c("30TWL","30TWM","30TXL","30TXM","30TXN","30TYL","30TYM")),
"Otros_Territorios" = list(tiles = c("30STE","30SUE","30SVE","30SWE")),
"Ceuta" = list(tiles = c("30STE")),
"Melilla" = list(tiles = c("30SVE","30SWE")),
"SCT" = list(tiles = c("27RYL","28RBS","28RCS")),
"LPGC" = list(tiles = c("28RDR","28RDS","28RES","28RFS","28RFT"))
)
provincias <- c("Ceuta")
provincia <- "Ceuta"
print(provincia)
teselas <- paste0(diccionario_provincias[[provincia]]$tiles)
proj <- "+init=EPSG:25830"
descarga <- "Agua_y_Arroz/Teselas"
#descarga <- "Arroz"
anios <- c("2019","2020","2021")  #2018, 2019, 2020
trimestres <- c("01","04","07","10") #01, 04, 07, 10
#Capa de provincias.
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
#RUTAS
path <- paste0("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/",descarga)
dir_salida <- paste0(path,"/outputs/")
x <- list.files(path, pattern = ".tif", full.names = TRUE)
anio <- "2019"
for(trimestre in trimestres){
lista = list()
for(tesela in teselas){
for(i in x){
if(grepl(paste0(tesela,"_",anio,trimestre),i)){
lista = append(lista,i)
}
}
}
#Mosaicamos
rasters <- lapply(lista,raster)
r_mosaic <- do.call(mosaic, c(rasters, fun='max'))
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
}
lista
rasters <- raster(lista)
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_201901.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_201904.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
trimestre<-"04"
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
trimestre<-"07"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_201907.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
trimestre<-"10"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_201910.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
anio <- "2020"
trimestre <- "01"
lista = list()
for(tesela in teselas){
for(i in x){
if(grepl(paste0(tesela,"_",anio,trimestre),i)){
lista = append(lista,i)
}
}
}
lista
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202001.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
trimestre <- "04"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202004.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
trimestre <- "07"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202007.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
trimestre <- "10"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202010.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
anio <- "2021"
trimestre <- "01"
lista = list()
for(tesela in teselas){
for(i in x){
if(grepl(paste0(tesela,"_",anio,trimestre),i)){
lista = append(lista,i)
}
}
}
lista
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202101.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
trimestre <- "04"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202104.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
trimestre <- "07"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202107.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
trimestre <- "10"
r_mosaic <- raster("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/Agua_y_Arroz/Teselas/30STE_202110.tif")
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.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(
"Alava" = list(tiles = c("30TVN","30TWN")),
"Albacete" = list(tiles = c("30SWH","30SWJ","30SXH","30SXJ")),
"Alicante" = list(tiles = c("30SXG", "30SXH", "30SXJ","30SYH","30SYJ")),
"Almeria" = list(tiles = c("30SVE","30SVF","30SVG","30SWF","30SWG","30SWH","30SXG")),
"Asturias" = list(tiles = c("29TPH","29TPJ","29TQH","29TQJ","30TUN","30TUP")),
"Avila" = list(tiles = c("30TTK","30TTL","30TUK","30TUL")),
"Badajoz" = list(tiles = c("29SPC","29SPD","29SQC","29SQD","30STH","30STJ","30SUH","30SUJ")),
"Baleares" = list(tiles = c("31SCC","31SCD","31SDD","31SED","31TDE","31TEE","31TFE")),
"Barcelona" = list(tiles = c("31TCF","31TCG","31TDF","31TDG")),
"Burgos" = list(tiles = c("30TUM","30TUN","30TVL","30TVM","30TVN","30TWM","30TWN")),
"Caceres" = list(tiles = c("29SPD","29SQD","29TPE","29TQE","30STJ","30SUJ","30TTK","30TUK")),
"Cadiz" = list(tiles = c("29SQA","30STE","30STF","30STG","30SUF","30SUG")),
"Cantabria" = list(tiles = c("30TUN","30TUP","30TVN","30TVP")),
"Castellon" = list(tiles = c("30TXK","30TYK","30TYL","31TBE")),
"CiudadReal" = list(tiles = c("30SUH","30SUJ","30SVH","30SVJ","30SWH","30SWJ")),
"Cordoba" = list(tiles = c("30STG","30STH","30SUG","30SUH","30SVG")),
"Cuenca" = list(tiles = c("30SVJ","30SWJ","30SXJ","30TVK","30TWK","30TWL","30TXK")),
"Girona" = list(tiles = c("31TCG","31TCH","31TDG","31TDH","31TEG","31TEH")),
"Gipuzkoa" = list(tiles = c("30TWN","30TWP")),
"Granada" = list(tiles = c("30SUF","30SUG","30SVF","30SVG","30SWG","30SWH")),
"Guadalajara" = list(tiles = c("30TVK","30TVL","30TWK","30TWL","30TXK","30TXL")),
"Huelva" = list(tiles = c("29SPB","29SPC","29SQA","29SQB","29SQC","30STF","30STG","30STH")),
"Huesca" = list(tiles = c("30TXM","30TXN","30TYL","30TYM","30TYN","31TCG","31TCH")),
"Jaen" = list(tiles = c("30SUG","30SUH","30SVG","30SVH","30SWG","30SWH")),
"Coruna" = list(tiles = c("29TMH","29TNH","29TNJ","29TPJ")),
"LaRioja" = list(tiles = c("30TVM","30TVN","30TWM","30TWN","30TXM")),
"Leon" = list(tiles = c("29TPG","29TPH","29TQH","30TTM","30TUM","30TUN")),
"Lleida" = list(tiles = c("30TYL","30TYM","30TYN","31TCF","31TCG","31TCH","31TDG")),
"Lugo" = list(tiles = c("29TNG","29TNH","29TNJ","29TPG","29TPH","29TPJ")),
"Madrid" = list(tiles = c("30TUK","30TUL","30TVK","30TVL")),
"Malaga" = list(tiles = c("30STF","30SUF","30SUG","30SVF")),
"Murcia" = list(tiles = c("30SWG","30SWH","30SXG","30SXH")),
"Navarra" = list(tiles = c("30TWM","30TWN","30TWP","30TXM","30TXN","30TXP")),
"Ourense" = list(tiles = c("29TNG","29TNH","29TPG","29TPH")),
"Palencia" = list(tiles = c("30TUM","30TUN","30TVM","30TVN")),
"Pontevedra" = list(tiles = c("29TMG","29TMH","29TNG","29TNH")),
"Salamanca" = list(tiles = c("29TPE","29TPF","29TQE","30TTK","30TTL","30TUL")),
"Segovia" = list(tiles = c("30TUL","30TVL","30TVM")),
"Sevilla" = list(tiles = c("29SQA","29SQB","30STF","30STG","30STH","30SUF","30SUG")),
"Soria" = list(tiles = c("30TVL","30TVM","30TWL","30TWM")),
"Tarragona" = list(tiles = c("30TYK","30TYL","31TCE","31TCF","31TCG")),
"Teruel" = list(tiles = c("30TWK","30TXK","30TXL","30TYK","30TYL")),
"Toledo" = list(tiles = c("30SUJ","30SVJ","30SWJ","30TTK","30TUK","30TVK")),
"Valencia" = list(tiles = c("30SXH","30SXJ","30SYH","30SYJ","30TXK","30TYK")),
"Valladolid" = list(tiles = c("30TTM","30TUL","30TUM","30TVL","30TVM")),
"Bizkaia" = list(tiles = c("30TVN","30TVP","30TWN","30TWP")),
"Zamora" = list(tiles = c("29TPG","30TTL","30TTM","30TUL","30TUM")),
"Zaragoza" = list(tiles = c("30TWL","30TWM","30TXL","30TXM","30TXN","30TYL","30TYM")),
"Otros_Territorios" = list(tiles = c("30STE","30SUE","30SVE","30SWE")),
"Ceuta" = list(tiles = c("30STE")),
"Melilla" = list(tiles = c("30SVE","30SWE")),
"SCT" = list(tiles = c("27RYL","28RBS","28RCS")),
"LPGC" = list(tiles = c("28RDR","28RDS","28RES","28RFS","28RFT"))
)
diccionario_provincias
nombre <- "Mario"
nombre
lista <- list("Moi",10,"hola")
lista
class(nombre)
class(lista)
#Capa de provincias.
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
plot(provincias)
provincias <- c("SCT")
#Capa de provincias.
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
View(provincias)
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
provincias <- c("SCT")
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/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/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
provincia <- "SCT"
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
plot(capa_prov)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
x
tesela
i
teselas
i
x
#LIBRERIAS
require(rgdal)
require(raster)
library(doParallel)
#################################################################################################################
UseCores <- detectCores() -1
#Register CoreCluster
cl<- makeCluster(UseCores)
start_time <- Sys.time()
#VARIABLES
diccionario_provincias <- list(
"Alava" = list(tiles = c("30TVN","30TWN")),
"Albacete" = list(tiles = c("30SWH","30SWJ","30SXH","30SXJ")),
"Alicante" = list(tiles = c("30SXG", "30SXH", "30SXJ","30SYH","30SYJ")),
"Almeria" = list(tiles = c("30SVE","30SVF","30SVG","30SWF","30SWG","30SWH","30SXG")),
"Asturias" = list(tiles = c("29TPH","29TPJ","29TQH","29TQJ","30TUN","30TUP")),
"Avila" = list(tiles = c("30TTK","30TTL","30TUK","30TUL")),
"Badajoz" = list(tiles = c("29SPC","29SPD","29SQC","29SQD","30STH","30STJ","30SUH","30SUJ")),
"Baleares" = list(tiles = c("31SCC","31SCD","31SDD","31SED","31TDE","31TEE","31TFE")),
"Barcelona" = list(tiles = c("31TCF","31TCG","31TDF","31TDG")),
"Burgos" = list(tiles = c("30TUM","30TUN","30TVL","30TVM","30TVN","30TWM","30TWN")),
"Caceres" = list(tiles = c("29SPD","29SQD","29TPE","29TQE","30STJ","30SUJ","30TTK","30TUK")),
"Cadiz" = list(tiles = c("29SQA","30STE","30STF","30STG","30SUF","30SUG")),
"Cantabria" = list(tiles = c("30TUN","30TUP","30TVN","30TVP")),
"Castellon" = list(tiles = c("30TXK","30TYK","30TYL","31TBE")),
"CiudadReal" = list(tiles = c("30SUH","30SUJ","30SVH","30SVJ","30SWH","30SWJ")),
"Cordoba" = list(tiles = c("30STG","30STH","30SUG","30SUH","30SVG")),
"Cuenca" = list(tiles = c("30SVJ","30SWJ","30SXJ","30TVK","30TWK","30TWL","30TXK")),
"Girona" = list(tiles = c("31TCG","31TCH","31TDG","31TDH","31TEG","31TEH")),
"Gipuzkoa" = list(tiles = c("30TWN","30TWP")),
"Granada" = list(tiles = c("30SUF","30SUG","30SVF","30SVG","30SWG","30SWH")),
"Guadalajara" = list(tiles = c("30TVK","30TVL","30TWK","30TWL","30TXK","30TXL")),
"Huelva" = list(tiles = c("29SPB","29SPC","29SQA","29SQB","29SQC","30STF","30STG","30STH")),
"Huesca" = list(tiles = c("30TXM","30TXN","30TYL","30TYM","30TYN","31TCG","31TCH")),
"Jaen" = list(tiles = c("30SUG","30SUH","30SVG","30SVH","30SWG","30SWH")),
"Coruna" = list(tiles = c("29TMH","29TNH","29TNJ","29TPJ")),
"LaRioja" = list(tiles = c("30TVM","30TVN","30TWM","30TWN","30TXM")),
"Leon" = list(tiles = c("29TPG","29TPH","29TQH","30TTM","30TUM","30TUN")),
"Lleida" = list(tiles = c("30TYL","30TYM","30TYN","31TCF","31TCG","31TCH","31TDG")),
"Lugo" = list(tiles = c("29TNG","29TNH","29TNJ","29TPG","29TPH","29TPJ")),
"Madrid" = list(tiles = c("30TUK","30TUL","30TVK","30TVL")),
"Malaga" = list(tiles = c("30STF","30SUF","30SUG","30SVF")),
"Murcia" = list(tiles = c("30SWG","30SWH","30SXG","30SXH")),
"Navarra" = list(tiles = c("30TWM","30TWN","30TWP","30TXM","30TXN","30TXP")),
"Ourense" = list(tiles = c("29TNG","29TNH","29TPG","29TPH")),
"Palencia" = list(tiles = c("30TUM","30TUN","30TVM","30TVN")),
"Pontevedra" = list(tiles = c("29TMG","29TMH","29TNG","29TNH")),
"Salamanca" = list(tiles = c("29TPE","29TPF","29TQE","30TTK","30TTL","30TUL")),
"Segovia" = list(tiles = c("30TUL","30TVL","30TVM")),
"Sevilla" = list(tiles = c("29SQA","29SQB","30STF","30STG","30STH","30SUF","30SUG")),
"Soria" = list(tiles = c("30TVL","30TVM","30TWL","30TWM")),
"Tarragona" = list(tiles = c("30TYK","30TYL","31TCE","31TCF","31TCG")),
"Teruel" = list(tiles = c("30TWK","30TXK","30TXL","30TYK","30TYL")),
"Toledo" = list(tiles = c("30SUJ","30SVJ","30SWJ","30TTK","30TUK","30TVK")),
"Valencia" = list(tiles = c("30SXH","30SXJ","30SYH","30SYJ","30TXK","30TYK")),
"Valladolid" = list(tiles = c("30TTM","30TUL","30TUM","30TVL","30TVM")),
"Bizkaia" = list(tiles = c("30TVN","30TVP","30TWN","30TWP")),
"Zamora" = list(tiles = c("29TPG","30TTL","30TTM","30TUL","30TUM")),
"Zaragoza" = list(tiles = c("30TWL","30TWM","30TXL","30TXM","30TXN","30TYL","30TYM")),
"Otros_Territorios" = list(tiles = c("30STE","30SUE","30SVE","30SWE")),
"Ceuta" = list(tiles = c("30STE")),
"Melilla" = list(tiles = c("30SVE","30SWE")),
"SCT" = list(tiles = c("27RYL","28RBS","28RCS")),
"LPGC" = list(tiles = c("28RDR","28RDS","28RES","28RFS","28RFT"))
)
provincias <- c("Badajoz","Avila","Almeria","Alicante","Albacete","Alava","Melilla","Otros_Territorios","Zaragoza","Zamora","Valladolid","Bizkaia","CiudadReal","Lleida")
provincias <- c("SCT")
for(provincia in provincias){
print(provincia)
teselas <- paste0(diccionario_provincias[[provincia]]$tiles)
proj <- "+init=EPSG:4326"
descarga <- "Agua_y_Arroz/Teselas"
#descarga <- "Arroz"
anios <- c("2029","2020","2021")  #2018, 2019, 2020
trimestres <- c("01","04","07","10") #01, 04, 07, 10
#Capa de provincias.
provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_CAN_WGS84.shp")
#provincias <- readOGR("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/2_Cob/recintos_provinciales_25830.shp")
capa_prov <- provincias[provincias$NAMEUNIT ==provincia,]
#################################################################################################################
#RUTAS
path <- paste0("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/04_Descargas/Descargas/",descarga)
dir_salida <- paste0(path,"/outputs/")
#################################################################################################################
x <- list.files(path, pattern = ".tif", full.names = TRUE)
#print(x)
#aqui metemos las imgs para mosaicarlas.
for(anio in anios){
for(trimestre in trimestres){
lista = list()
for(tesela in teselas){
for(i in x){
if(grepl(paste0(tesela,"_",anio,trimestre),i)){
lista = append(lista,i)
}
}
}
#Mosaicamos
rasters <- lapply(lista,raster)
r_mosaic <- do.call(mosaic, c(rasters, fun='max'))
#Recortamos
capa <- mask(r_mosaic, capa_prov)
capa_clip <- crop(capa,capa_prov)
#Exportamos
if(file.exists(paste0(dir_salida,provincia))){
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
} else{
dir.create(paste0(dir_salida,provincia))
writeRaster(capa_clip, filename = paste0(dir_salida,provincia,"/",provincia,"_",anio,trimestre,".tif"),format='GTiff')
}
}
}
end_time = Sys.time()
print(end_time - start_time)
}
x
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
require(rgdal)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/LULUCF/PRODUCCION_LULUCF21/03_Scripts/Mosaico_Recorte_Provincia.R", echo=TRUE)
require(sf)
salida <- "D:/PRUEBAS/apanio_pedro/"
#########################################################################################################
#ZNT
path <- "//Imgwebserver/proyectos/FEGA/AMS/2024_AMS/JEG/ZNT/LIMITES"
huso <- "/UTM29/"
f <- "/LIMITES_CAPTURA_ETRS89-UTM29.shp"
#Formo la ruta de la capa:
capa <- paste0(path,huso,f)
#Cargo la capa:
capa_2 <- st_read(capa)
id <- list(capa_2$Id)
#df vacio
df2 <- data.frame(id=numeric(),
area_km2=numeric(),
zona = character(),
xmin = numeric(),
ymin = numeric(),
xmax = numeric(),
ymax = numeric(),
stringsAsFactors=FALSE)
ids <- length(capa_2$Id)
for(i in 1:ids){
fila_vec <- capa_2[capa_2$Id == capa_2$Id[i],]
fila_ext <- st_bbox(fila_vec)
print(fila_ext)
fila_ext_df <- c(fila_ext[1],fila_ext[2],fila_ext[3],fila_ext[4])
fila_ext_df <- as.data.frame(fila_ext_df)
fila_ext_df <- t(fila_ext_df)
#fila_df <- as.data.frame(fila)
#fila_id <- select(fila, Id, AREA_KM2, ZONA)
fila_vec_df <- as.data.frame(fila_vec)
fila <- fila_vec_df[,-ncol(fila_vec_df)]
fila_f <- cbind(fila,fila_ext_df)
df2 <- rbind(df2,fila_f)
}
write.csv2(df2, paste0(salida,"LIMITES_CAPTURA_ETRS89-UTM29.csv"))
#library(gdalUtils)
#require(rgdal)
#library(stars)
library(sf)
library(raster)
#library(exactextractr)
#library(cluster)
library(tidyverse)
library(DBI)
library(RPostgres)
# 2.0 Parametros de conexion a PostgreSQL
dvr_ <- RPostgres::Postgres()
db_ <- 'eikos_trabajo_2022'  ##Nombre de la BBDD
host_db_ <- '172.17.12.4'
db_port_ <- '5433'
db_user_ <- 'teledeteccion'  ##Tu usuario
db_password_<- 'teledeteccion-2023-dev' ##Tu contrasena
schema <-"-c search_path=tel_anomalias"
# 3.0 Conexion
con <- dbConnect(dvr_, dbname = db_, host=host_db_, port=db_port_,
user=db_user_, password=db_password_, options= schema)
# 4.0 Listado de tablas de la Base de Datos
dbListTables(con)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/2_ORG_DESC.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/3_CONTROL_DESCARGAS.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/4_UNSTACK_WITHNAMES_v2.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/5_Grafico_Tablas_CRUCE_SUP_v2_COBNUEVA_CURSO.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/6_LAGUNA_SECA_v3_2023.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/2_ORG_DESC.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/3_CONTROL_DESCARGAS.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/2_ORG_DESC.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/3_CONTROL_DESCARGAS.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/4_UNSTACK_WITHNAMES_v2.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/5_Grafico_Tablas_CRUCE_SUP_v2_COBNUEVA_CURSO.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/6_LAGUNA_SECA_v3_2023.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/2_ORG_DESC.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/3_CONTROL_DESCARGAS.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/4_UNSTACK_WITHNAMES_v2.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/5_Grafico_Tablas_CRUCE_SUP_v2_COBNUEVA_CURSO.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/4_UNSTACK_WITHNAMES_v2.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/4_UNSTACK_WITHNAMES_v2.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/5_Grafico_Tablas_CRUCE_SUP_v2_COBNUEVA_CURSO.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/5_Grafico_Tablas_CRUCE_SUP_v2_COBNUEVA_CURSO.R", echo=TRUE)
source("//imgwebserver/PROYECTOS/MITECO/DG_AGUA/LAGOS/CURSO_2023/3_Scripts/1_Procesado/6_LAGUNA_SECA_v3_2023.R", echo=TRUE)
