search for: rastervis

Dear useRs, I'd like to announce the release of version 0.10-3 of rasterVis. This package provides a set of methods for enhanced visualization and interaction with the Raster* objects from the raster package. You will find some examples at: http://rastervis.r-forge.r-project.org/. Best, Oscar. ------------- Oscar Perpi??n Lamigueiro Dpto. de Ingenier?a El?ctrica EUITI-U...

Dear useRs, I'd like to announce the release of version 0.10-3 of rasterVis. This package provides a set of methods for enhanced visualization and interaction with the Raster* objects from the raster package. You will find some examples at: http://rastervis.r-forge.r-project.org/. Best, Oscar. ------------- Oscar Perpi??n Lamigueiro Dpto. de Ingenier?a El?ctrica EUITI-U...

Saludos estimados usuarios ! Estoy recién iniciándome en R y no he encontrado ejemplos de cómo graficar los archivos netcdf que genera el modelo meterológico WRF. ¿Alguien tiene por ahí algunos scripts que me permitan ir interiorizándome en el uso de paletas de colores, mapas y graficas de variables de este modelo? Anteriormente lo hacía vía GrADS, convirtiendo el NetCDf a GRIB, pero con las

summary: I can display a lon-lat map on a lattice::layerplot, and I can display a Lambert conformal conic (LCC) map on a spam::image, but I can't display an LCC map on a lattice::layerplot. Example follows. What am I doing wrong? details: I've been using `lattice` (via `rasterVis`) successfully to display global atmospheric data, which works well enough (though I am definitely intrigued by ggplot/ggmap). Particularly, I am able to overlay my plots with maps, which is essential for the sort of work I'm doing. However I'm currently unable to use lattice::layerplot for...

...max(map.IOAPI$x), > min(map.IOAPI$y), > max(map.IOAPI$y)) > map.IOAPI.shp <- > maptools::map2SpatialLines(map.IOAPI, proj4string=CRS) > return(map.IOAPI.shp) > } # end project.M3.boundaries.for.CMAQ The maps are then overlaid on data and plotted using package=rasterVis http://cran.r-project.org/web/packages/rasterVis/ with code like https://bitbucket.org/tlroche/gfed-3.1_global_to_aqmeii-na/src/95484c5d63502ab146402cedc3612dcdaf629bd7/vis_regrid_vis.r?at=master > visualize.layers( > nc.fp=monthly_out_fp, > datavar.name=monthly_out_datavar_name, &g...

Hola miembros de la lista, Estoy utilizando la función levelplot del paquete rasterVis para graficar un raster y quiero controlar la altura de la referencia de color (colorkey) del mapa, pero cuando cambio valores en el argumento "height" no parece provocar cambios. Alguna idea de por qué no está funcionando? > cobertura class : RasterLayer dimensions : 780, 129...

You need to load the package 'rasterVis' > library(rasterVis) HTH, Eric On Thu, Mar 8, 2018 at 5:11 PM, lily li <chocold12 at gmail.com> wrote: > Hi all, > > I ran the code: > > s <- stack(replicate(2, raster(matrix(runif(100), 10)))) > > xy <- data.frame(coordinates(sampleRandom(s, 10, sp=TRU...

I have a 3D field of a scalar variable (x, y, z, value). Is there a way to generate a vector field from this data--gradient at defined points? I found the rasterVis package for 2D data, but as yet nothing for 3D data. Thanks, Scott Waichler Pacific Northwest National Laboratory Richland, WA scott.waichler at pnnl.gov

Hi all, I ran the code: > s <- stack(replicate(2, raster(matrix(runif(100), 10)))) > xy <- data.frame(coordinates(sampleRandom(s, 10, sp=TRUE)), + z1=runif(10), z2=runif(10)) > levelplot(s, margin=FALSE, at=seq(0, 1, 0.05)) + + layer(sp.points(xy, pch=ifelse(pts$z1 < 0.5, 2, 3), cex=2, col=1), columns=1) + + layer(sp.points(xy, pch=ifelse(pts$z2 < 0.5, 2,

Hi all, I'm trying to add single points with known coordinates to a level plot, but could not find the proper answer. I got to know that layer() function is good for this, but I don't know which package is related to this function. The source is here: https://stackoverflow.com/questions/28597149/add-xy-points-to-raster-map-generated-by-levelplot but my question is a little different as I

This code will read binary file and display it as a map. may problem is that this code is using a continuous colour scheme, even though I have discrete data (which is a classification scheme). How can I map numbers to colours with raster? Please require(raster) conne <- file("C:\\lai.bin", "rb") sd<- readBin(conne, integer(), size=1, n=360*720, signed=F)

...To fix this, you will need reverse dimension again. And also if your longitudes ranges from 0 to 360, you will see the only the east side of the plot on a world map. West of Greencwich will be empty. They were the problems that I faced last year when I tried to plot netcdf files using lattice and rasterVis packages. > > As I the said, I haven't done extensive testing on what Bert sent, but on a toy 3-dimensional example I have it appeared to do what I want. > > Thanks again, > > -Roy > >> On Jun 1, 2017, at 12:22 PM, David L Carlson <dcarlson at tamu.edu>...

...To fix this, you will need reverse dimension again. And also if your longitudes ranges from 0 to 360, you will see the only the east side of the plot on a world map. West of Greencwich will be empty. ?They were the problems that I faced last year when I tried to plot netcdf files using lattice and rasterVis packages.? As I the said, ?I haven't done extensive testing on what Bert sent, ?but on a toy 3-dimensional example I have it appeared to do what I want. Thanks again, -Roy On Jun 1, 2017, at 12:22 PM, David L Carlson <dcarlson at tamu.edu> wrote: My error. Clearly I did not do en...

...To fix this, you will need reverse dimension again. And also if your longitudes ranges from 0 to 360, you will see the only the east side of the plot on a world map. West of Greencwich will be empty. They were the problems that I faced last year when I tried to plot netcdf files using lattice and rasterVis packages. > > > > > As I the said, I haven't done extensive testing on what Bert sent, but on a toy 3-dimensional example I have it appeared to do what I want. > > Thanks again, > > -Roy > > > On Jun 1, 2017, at 12:22 PM, David L Carlson <dcarls...

Thanks to all for responses/. There was a question of exactly what was wanted. It is the generalization of the obvious example I gave, >>> junk1 <- junk[, rev(seq_len(10), ] so that junk[1,1,1 ] = junk1[1,10,1] junk[1,2,1] = junk1[1,9,1] etc. The genesis of this is the program is downloading data from a variety of sources on (time, altitude, lat, lon) coordinates, but all

My error. Clearly I did not do enough testing. z <- array(1:24,dim=2:4) > all.equal(f(z,1),f2(z,1)) [1] TRUE > all.equal(f(z,2),f2(z,2)) [1] TRUE > all.equal(f(z,3),f2(z,3)) [1] "Attributes: < Component ?dim?: Mean relative difference: 0.4444444 >" [2] "Mean relative difference: 0.6109091" # Your earlier example > z <- array(1:120, dim=2:5) >