R Notebook

PISCO MENSUAL

library(sp)
library(rgdal)
library(Rcpp) # for R 4.x versions, try also package “terra”
library(raster)
library(ncdf4)
# Definiendo un directorio de trabajo donde se encuentren los archivos input descargados
setwd("C:/Users/DELL/Desktop/UTEC 2022-2/Hidrología/VILCANOTA_PISCO/shapefile")

Warning: The working directory was changed to C:/Users/DELL/Desktop/UTEC 2022-2/Hidrología/VILCANOTA_PISCO/shapefile inside a notebook chunk. The working directory will be reset when the chunk is finished running. Use the knitr root.dir option in the setup chunk to change the working directory for notebook chunks.

# Leyendo el polígono CuencaVeronica.shp (layer: CuencaVeronica) desde la carpeta “files” con data source name (dsn)
cuenca.shape <- readOGR(dsn="veronica", layer="CuencaVeronica")

OGR data source with driver: ESRI Shapefile 
Source: "C:\Users\DELL\Desktop\UTEC 2022-2\Hidrología\VILCANOTA_PISCO\shapefile\veronica", layer: "CuencaVeronica"
with 1 features
It has 3 fields
Integer64 fields read as strings:  id 

# Conociendo el tipo de objeto después de la importación
class(cuenca.shape)

[1] "SpatialPolygonsDataFrame"
attr(,"package")
[1] "sp"

# Visualizando la cuenca con color cyan, se aprecia que se encuentra en coordenadas geograficas
plot(cuenca.shape, axes=T, col=c("green"))

# Conociendo el contenido del polígono importado
head(cuenca.shape@data)

# Transformando a UTM zona 18, WGS84 y visualizando la conversión
cuenca.utm <- spTransform(cuenca.shape, CRS("+proj=utm +zone=18 +ellps=WGS84 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"))
plot(cuenca.utm, axes=T, asp=1)

help("plot")
# Reproyectando a Geograficas WGS84 y visualizando la reconversión
cuenca.wgs <- spTransform(cuenca.utm, CRS("+proj=longlat +ellps=WGS84"))
plot(cuenca.wgs, axes=T, asp=1, cex.axis=1, las=1, col="dark green", main="Mapa de la cuenca Verónica")

# Extrayendo datos y visualizando el producto ráster PISCO
r <- stack("PISCOV3-MONTHLY.nc") # El archivo *.nc debe estar en la ruta directa setwd

# Visualizando la precipitación espacial del primer mes (enero) de 1981
plot(r[[1]])

# Ploteando la cuenca Verónica dentro del mapa de precipitaciones
plot(cuenca.wgs, add=T)

# Visualizando la precipitación espacial del primer mes (febrero) de 1981
plot(r[[2]])

# Ploteando la cuenca del Veronica dentro del mapa de precipitaciones
plot(cuenca.wgs, add=T)

# Delimitando el área de estudio al cuadrante que ocupa la Cuenca Veronica
r.basin <- crop(r, cuenca.wgs, snap="out")

# Ploteando el primer mes (enero) de 1981 en el cuadrante que ocupa la Cuenca Veronica
plot(r.basin[[1]])
plot(cuenca.wgs, add=T)

# Delimitando el área de estudio a la cuenca de la Cuenca Veronica
r.basin <- mask(r.basin, cuenca.wgs)

# Ploteando los meses de enero a diciembre de 1981
plot(r.basin[[1:12]])

  
# Extrayendo los datos y promediando todas las grillas de la cuenca de la Cuenca Veronica
pp.cuenca.mensual <- extract(r.basin, cuenca.wgs, fun=mean)
colnames(pp.cuenca.mensual) <- 1:ncol(pp.cuenca.mensual)

# Visualizando los 431 datos promediados
View(pp.cuenca.mensual)
range(pp.cuenca.mensual)

[1]  11.2841 253.9020

# Ploteando la serie de los 431 valores de precipitación mensual promedio areal
plot(pp.cuenca.mensual[1,], type="o", col="1", ylim=c(0,300), ylab="P (mm)", xlab = "Meses", main="Precipitacion promedio areal - Vilcanota (mm)")

Pasando datos mensaules a excel

# todos meses desde 1981
Pisco <- data.frame(pp.cuenca.mensual[1,])
write.csv(Pisco, "C:/Users/DELL/Desktop/UTEC 2022-2/Hidrología/Hidro R/Entregable 2/Pisco.csv", row.names=FALSE)