Tengo dos rásteres con la misma resolución (10 m), pero la diferencia es que creé el primer ráster a partir del archivo vectorial (usando la función rasterize) y el otro ráster es el resultado de mi clasificación. Quiero aplicar la máscara de ráster obtenida del vector al resultado de la . . . Read more
Estoy intentando calcular un promedio del parámetro sobre una región específica (longitud=0:30, latitud=40:70). Mis datos se pueden encontrar aquí en formato .nc: https://downloads.psl.noaa.gov/Datasets/ncep.reanalysis.dailyavgs/pressure/hgt.2006.nc Lo que intenté hacer: biblioteca(ncdf4) biblioteca(raster) biblioteca(rgdal) r_brick
He estado ejecutando el siguiente código en RStudio para probar la funcionalidad de vectorplot en R: library(raster) library(rasterVis) library(lattice) <p>proj <- CRS(‘+proj=longlat +datum=WGS84’)</p> <p>df <- expand.grid(x=seq(-2, 2, .01), y=seq(-2, 2, .01)) df$z <- with(df, (3<em>x^2 + y)</em>exp(-x^2-y^2)) r1 <- rasterFromXYZ(df, crs=proj) df$z <- with(df, x<em>exp(-x^2-y^2)) r2 <- rasterFromXYZ(df, crs=proj) df$z . . . Read more