[R-sig-Geo] writeRaster and RGB vs Gray

[MacQueen, Don]

Here's a little example script using functions from the raster package,

r1 <- r2 <- r3 <- raster(ncol=25,nrow=25,
     xmn=0, ymn=0, xmx=1, ymx=1)
r1[] <- round(runif(25^2,0,255))
r2[] <- round(runif(25^2,0,255))
r3[] <- round(runif(25^2,0,255))
s <- stack(r1,r2,r3)

plotRGB(s)

writeRaster(s,'s.tif', format='GTiff', overwrite=TRUE,
            bandorder='BSQ')


When I open the file s.tif in QGIS, it's fine; I see the same image as in
R.

But the Mac platform's default viewer, Preview.app, thinks the file's
color model is "Gray", rather than "RGB", and it does not render properly.
Preview does recognize, for example, a jpeg photo as having an RGB color
model.
(Adobe Reader (for Mac) wouldn't open it at all.)

Is there some better set of arguments to give writeRaster that will
somehow make the file more broadly recognizable as RGB?

-- 
Don MacQueen

Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062