[R-sig-Geo] select a pixel with the right value respect the pixel around

2009-04-04 Thread gianni lavaredo
Dear R user, I am trying several methods to find the right solution. I have a raster DEM grid. Is it possible to select all pixel where the pixel value around (8 pixel or 24 or more) is >2/3 of value my pixel x? Gianni [[alternative HTML version deleted]] ___

Re: [R-sig-Geo] error in trmesh (alphahull package)

2009-04-04 Thread Murray Richardson
Ahhh - so simple, my mistake. Thanks! Murray José Manuel Blanco Moreno wrote: Hello, The problem is that your first three points CANNOT be collinear. This is a feature of trmesh in tripack. So: either you "unsort" your data (keeping the indices you can sort them again) or alter (e.g. jitter)

[R-sig-Geo] hscattergram between two different variables

2009-04-04 Thread facri...@libero.it
Dear all, I want to plot the cross hscattergram between two different variables. I used the following instruction: hscat(jura.pred$Cd~jura.pred$Zn, jura.pred, c(0, 0.1, 0.2, 0.3, 0.5, 1), alpha=c(0,45)) where jura.pred$Cd and jura.pred$Ni are the two variables, but the result obtained is th

[R-sig-Geo] extract value from grid and calculate weighted.mean of pixel

2009-04-04 Thread gianni lavaredo
Dear R user I ask some suggestio about this code that I wish to write. I have a raster grid a polygons shape files. I need to extract the pixel value inside the shape polygon and calculate the weight.mean and associate that value to polygon raster. Is there a package useful to do this that I could

Re: [R-sig-Geo] County-level map of the US

2009-04-04 Thread Corey Sparks
Hi Phil, I saw your post and have been doing this exact thing for the US recently. Do you have a county shapefile? You will need this to join your data (Y) to. I typically use ArcGIS to join my attribute data to my shapefiles, but I don't know if you have access to this software. This

Re: [R-sig-Geo] error in trmesh (alphahull package)

2009-04-04 Thread José Manuel Blanco Moreno
Hello, The problem is that your first three points CANNOT be collinear. This is a feature of trmesh in tripack. So: either you "unsort" your data (keeping the indices you can sort them again) or alter (e.g. jitter) the coordinates very slightly. Either solution will do the work. Cheers, José M

Re: [R-sig-Geo] rasterize shape file and output as .img

2009-04-04 Thread Roger Bivand
On Fri, 3 Apr 2009, milton ruser wrote: Hi Roger, I think I am making something wrong (see the error message, please). reser is the polygon object that I read from a shapefile. Thanks for your time. milton class(reser) [1] "SpatialPolygonsDataFrame" attr(,"package") [1] "sp" head(re...@da