Re: [R-sig-Geo] Clip a raster

2009-02-05 Thread Robert Hijmans
Etienne, I added a function clickBbox to the raster package for increased interactivity #install.packages(raster, repos=http://R-Forge.R-project.org;) require(raster) # read a raster file from disk filename - system.file(pictures/Rlogo.jpg, package=rgdal) r - rasterFromFile(filename, TRUE)

[R-sig-Geo] ArcView shape files and GeoBUGS

2009-02-05 Thread Tord Snäll
Dear all, I would like to read in ArcView (or ArcGIS) shape files from which I will pick information if fitting a geostatistical model using WinBUGS. In the GeoBUGS Used Manual it says GeoBUGS does not have an option for loading ArcView shape files directly. However, Ms Yue Cui at the

Re: [R-sig-Geo] ArcView shape files and GeoBUGS

2009-02-05 Thread Robert Hijmans
I think this is it: http://www.biostat.umn.edu/~brad/yuecui/index.html http://www.biostat.umn.edu/~brad/yuecui/convert.r #/09/23/03 The only different btw this convert.r and convert.s is that #in convert.s, we us AsciiToInt(character)-AsciiToInt(0) to change #character into numeric, in

Re: [R-sig-Geo] symbols in spplot

2009-02-05 Thread Peter van Horssen
Thanks Roger, your trick worked, reworking the polygon data to a raster and plotting the file with image() and then plotting the rest with plot(, add=TRUE) works great (see attached code). It shows that the famous KISS principle (keep is simple, stupid ) also applies to plotting maps in

Re: [R-sig-Geo] Summary: ArcView shape files and GeoBUGS

2009-02-05 Thread Roger Bivand
On Thu, 5 Feb 2009, Tord Snäll wrote: Thanks to you all for prompt replies! As evident from the replies pasted below there are different solutions to this problem. cheers, Tord Thanks for a very full summary, keeping the threaded information together. Just two more comments - the next

[R-sig-Geo] krige problem: IDW + irregular blocks

2009-02-05 Thread Mauricio Zambrano
Dear List, I need to carry out an IDW interpolation of precipitation values within a basin. I want to obtain the mean value in each one of the subcatchments, which are defined as polygons in a shapefile. I would like to know: How can I transform the subcatchments polygons (class