Re: [R-sig-Geo] Question raster package

2009-06-23 Thread Peter van Horssen
lt;- raster(fn, band=1) r2 <- raster(fn, band=2) r3 <- raster(fn, band=3) you can also first create a RasterStack of all bands s <- stack(fn, bands=-1) and then extract a RasterLayer whenever you want it r <- raster(s, layer=1) Robert 2009/6/22 Peter van Horssen : I have a ras

[R-sig-Geo] Question raster package

2009-06-22 Thread Peter van Horssen
I have a raster object with 3 bands: (How) Can I obtain the info for each different band ? Can I use a command like reclass on all three bands (with different classification definitions for each band) ? thanks in advance, Peter -- Peter van Horssen Bureau

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

2009-02-05 Thread Peter van Horssen
,raster.midpoint$raster.mid.x2, ms.arrows, angle= Data$rad, inches=0.5, adj=0.5, lwd=1.5, code = 2, col = 'blue'), lty="blank", add=TRUE) ##code Roger Bivand schreef: On Wed, 4 Feb 2009, Peter van Horssen wrote: Dear all, I have a plot problem which is related to t

[R-sig-Geo] symbols in spplot

2009-02-04 Thread Peter van Horssen
UE) plot(turbines, pch=18, add= TRUE, col = 'red') plot(metmast, pch = 19, col = 'blue', add= TRUE) plot(raster,pch=my.symbols(raster.midpoint$raster.mid.x1,raster.midpoint$raster.mid.x2, ms.arrows, angle= Data$rad, inches=0.5, adj=0.5, lwd=1.5, code = 2, col = &