similar to: rasterize SpatialPolygon object using a RasterBrick object

Hi, I would like to do coloring of map regions based on the region values "weight". The approach I am taking is first to break regions into equal intervals, classIntervals(spdf$weight,4)$brks #4 intervals in this case and coloring all regions within the interval with the same color col = brewer.pal(4,"RdYlGn")) The max "weight" is as well the boundary of the

Dear List, I'm trying to use different R packages for my Teaching Assistantship classes. And I cam out to an (apparently) very simple problem. I would like to retrieve the vertices coordinate of a SpatialPolygon data. I know this is stored in the "coords" slot, but I can't get access to it! I tried to coerce the SpatialPolygon into a data.frame but it doesn't work.

Hi all, I have some concerns regarding the rasterize option in R and I would like to know if the fun=max in rasterize in R provides similar results to the one achieved by using "Polygon to Raster using maximum-combined-area" in ArcGIS? I'm trying to rasterize a habitat layer to a raster of 10m spatial resolution using the function 'max' (e.g. r <- rasterize(ht, r,

Dear listers, I am trying to create a RasterLayer of the values of a rasterbrick object. The rasterbrick object has, for example, 100cells library(raster) r <- raster(ncol=10, nrow=10) r[]=1:ncell(r) s <- brick(r,r,r) s <- s * 1:3 Each cell of the rasterfinal will have the AREA UNDER CURVE formed by the values in each cell of the three original rasters. For example: s[4] has the

Necesitaría ayuda con una cuestión. Pongamos por ejemplo estos valores aleatorios para generar un RasterBrick s de 3 capas. library(raster) r <- raster(ncol=10, nrow=10) r[]=1:ncell(r) s <- brick(r,r,r) s <- s * 1:3 Habría alguna forma sin hacer un bucle for, que me devolviese un vector con los valores medios de cada una de las 3 capas. mean_s<-NULL for (i in 1:3) {

i'm sure most of us are avid slashdot readers, but... http://slashdot.org/article.pl?sid=01/07/13/1558239&mode=nested have fun dongoodman %--------------------------------------------deg3@ra.msstate.edu| |bleu| |`I invented the term 'object oriented', pobox 2516| |----| |and I can tell you mississippi state ms| | () | /. |I did

Is there a Linux tool that will splice several gsm sound clips together into one clip? In my agi script, I would like to use 'get_data' with one clip instead of multiple 'stream_file' so the user doesn't have to listen to the entire spiel before responding. Thanks, -- Michael Welter Introspect Telephony Corp. Denver, Colorado +1 303 674 2575 mike@introspect.com

libswfdec/swfdec_movie.c | 18 - libswfdec/swfdec_sprite.c | 308 ++---------------------------- libswfdec/swfdec_sprite.h | 36 +-- libswfdec/swfdec_sprite_movie.c | 357 ++++++++++++++++++++++++++++------- libswfdec/swfdec_sprite_movie.h | 12 - libswfdec/swfdec_swf_decoder.c | 32 ++- libswfdec/swfdec_swf_decoder.h | 11 +

Hello. I am thinking of developing a tool for computer game makers, enabling them using an ogg file to hold a collection of very short context audio clips. I am looking for suggestion on usefulness of the tool once it is made, and its design. By "very short context audio" I mean the audio clips that usually is the response to a mouse-click or key press. e.g. a clicking sound, a

Hi R users, could someone help me to find out which formulas, for standard error calculation, are used in following example: a=data.frame(weights=rep(c(10,1),c(4,1)),fpc=rep(41,5),uk=rep(1,5)) srs<-svydesign(id=~1, weights=~weights, data=a) srs1<-svydesign(id=~1, weights=~weights,fpc=~fpc, data=a) svytotal(~uk,srs) total SE uk 41 9 svytotal(~uk,srs1) total SE uk 41

Hi, I have 2 huge datasets - May and Jun - a miniscule sample of one is given below. I am trying to do 2 things with these datasets. I need to verify if the weighted average of variable A for a Reason in Jun is same/different from the same for May. To do this I am first computing the weighted average for each SubReason using a function I wrote. Where I need help is applying the function on both

Hola Jorge, > My idea was to use 'stackApply': Para esto que necesitas es más sencillo usar calc. Si quieres usar stackApply, debes indicar rep(1, 3) como índices en lugar de sólo 1. Así, la función será (los ... son necesarios en stackApply): myFUN <- function(y, ...) { x <- seq_along(y) trapz(x, y) } Y aplicada con calc: FINAL <- calc(s, myFUN) o con

Hello everybody out there using R, When putting multiple plots with thousands of data points into a single PDF file, this file can get huge and take a long time to open. The following post describes exactly the same problem in Matplotlib, as well as a nice fix for it: Matplotlib: multipage PDF with rasterized plots Particularly nice about it is, that it only rasterizes the points without

Hi dear user, Anybody can tell me how extract the names of ID in SpatialPolygons object, I am try to link a data frame attributes with spatial polygons, but the row names of data frame is not the same that ID poligons. Thank for all Fernando [[alternative HTML version deleted]]

libswfdec/swfdec_codec_video.c | 2 libswfdec/swfdec_movie.c | 3 test/trace/Makefile.am | 193 ++++++++++++ test/trace/place-object-5.swf |binary test/trace/place-object-5.swf.trace | 4 test/trace/place-object-6.swf

Dear list, My favorite output format is usually pdf. I can include the graphics in pdflatex documents and benefit from the scalable nature of vector graphic formats. However, I recently had to generate high-res 2D levelplot graphics as in the example below, N <- 100 # N <- 1000 # slow to diplay xy <- expand.grid(x=seq(0, 10, length=N), y=seq(0, 10, length=N)) xy <- within(xy,

Dear R-users, I’m currently trying to extract the value of a slot (area) but can’t find out how to do that. >str(overlperc) List of 1 $ :Formal class 'SpatialPolygons' [package "sp"] with 4 slots .. ..@ polygons :List of 1 .. .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots .. .. .. .. ..@ Polygons :List of 1 .. .. .. .. .. ..$

Hi, I seem to be having a problem adding the following two raster objects together - one is a rasterLayer, the other is a rasterBrick. The extent, resolution, and origin are the same, so according to my understand it should work. The objects look like so: > obs.clim class : RasterLayer dimensions : 60, 200, 12000 (nrow, ncol, ncell) resolution : 0.5, 0.5 (x, y) extent : -70,

On Thu, 24 May 2018 at 18:41 Mark R Payne <markpayneatwork at gmail.com> wrote: Hi, I seem to be having a problem adding the following two raster objects together - one is a rasterLayer, the other is a rasterBrick. The extent, resolution, and origin are the same, so according to my understand it should work. The objects look like so: > obs.clim class : RasterLayer dimensions : 60, 200,

Dear all, In grDevices R functions cairo_pdf and cairo_ps it is mentioned that when transparency (alpha channels) are used in vector output, it will rasterize the PDF or postscript exported graph at a resolution of 72 dpi : https://stat.ethz.ch/R-manual/R-devel/library/grDevices/html/cairo.html You can see the problem if you try library(ggplot2) cairo_ps(file = "test.eps",onefile =