[R-sig-Geo] plot polygons - transparent holes/island

I'd like to plot a SPDF with two objects, one of them buffering the other island. How can I set the island transparent? Similarly, keeping the outer polygon only, how can the hole still transparent? Thanks in advance, Paulo E. Cardoso [[alternative HTML version deleted]]

Re: [R-sig-Geo] NDVI time-series filter

Dear Jonathan, Interesting topic. Please send us also your code example and if possible part of your dataset (I assume it is a SpatialGridDataFrame?). Operations on lists can be speed up by using lapply and by implementing your own functions, consider also running Rprof to check which

Re: [R-sig-Geo] kriging as fish swim, not as crows fly

Multidimensional Scaling (MDS) is certainly an option to be considered. Its suitability depends, of course, on the specific application. In my case, which is conceptually similar to Martin's, that option was discarded. And I think I can explain why whith this single image:

Re: [R-sig-Geo] Overlay spatial lines on raster

Thanks Robert, The raster package is a magician's chest! I get an error, though: class(xter) [1] SpatialLinesDataFrame attr(,package) [1] sp class(r) [1] RasterLayer attr(,package) [1] raster lr = linesToRaster(xter, r, field = ID_GRAFIC ) [1] A rows cols [1,]1 -1 [2,] -1 -1

Re: [R-sig-Geo] plot polygons - transparent holes/island

On Sat, 30 Jan 2010, Paulo Eduardo Cardoso wrote: I'd like to plot a SPDF with two objects, one of them buffering the other island. How can I set the island transparent? Similarly, keeping the outer polygon only, how can the hole still transparent? Polygons objects are plotted largest to

Re: [R-sig-Geo] Overlay spatial lines on raster

Hi Agus, thank you. You were right about the cause of the problem --- I patched it in version 0.9.9-5, now on R-Forge (and should be available for automatic install within 24 hours); the function could use some optimization for speed but at least it works now for SpatialLines* with a larger extent

Re: [R-sig-Geo] plot polygons - transparent holes/island

Roger, I really thank you for you comments. You are always available to explain everything, and that's a fortune for forum users. But I still with the same problem. I did,'t understand how to leave the inner hole or the inner polygon transparent, in order to see any other layer plotter at first.

[R-sig-Geo] generate simulation data for a theoretical spatial model

Hi all, In classical statistics, we always need to generate a theoretical model such as y=a+b1*x1+b2*x2+e to study some new estimation content. I am wondering how to generate the similar spatial dataset for a theoretical model. Say y is response variable, x1 and x2 are explanatory variables. 1.