search for: rastervis

Displaying 16 results from an estimated 16 matches for "rastervis".

Did you mean: rasterisk
2011 Jul 24
0
rasterVis 0.10-3
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...
2011 Jul 24
0
rasterVis 0.10-3
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...
2012 Apr 03
4
R+netcdf, resultados de WRF
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
2013 Feb 13
1
[lattice] display a projected map on a layerplot
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...
2013 Apr 26
1
[newbie] how to find and combine geographic maps with particular features?
...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...
2014 Jun 10
2
Como controlar la altura de "colorkey" en levelplot de RatserVis
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...
2018 Mar 08
1
add single points to a level plot
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...
2013 Mar 21
1
vector field from a 3D scalar field
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
2018 Mar 08
0
add single points to a level plot
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,
2018 Mar 08
3
add single points to a level plot
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
2012 Oct 29
3
How can I map numbers to colours with raster?
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)
2017 Jun 01
3
Reversing one dimension of an array, in a generalized case
...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>...
2017 Jun 01
0
Reversing one dimension of an array, in a generalized case
...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...
2017 Jun 01
1
Reversing one dimension of an array, in a generalized case
...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...
2017 Jun 01
0
Reversing one dimension of an array, in a generalized case
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
2017 Jun 01
2
Reversing one dimension of an array, in a generalized case
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) >