[R-sig-Geo] writeGDAL: pseudocolor image?

[Agustin Lobo]

Hi!
Is it possible to write a pseudocolor image as a  geotif with embedded
color table using
writeGDAL() or any other R command?
I mean an image such as:

$ gdalinfo pseudocolorv2.tif
Driver: GTiff/GeoTIFF
Files: pseudocolorv2.tif
       pseudocolorv2.tif.aux.xml
Size is 217, 232
Coordinate System is `'
Metadata:
  TIFFTAG_IMAGEDESCRIPTION=ImageJ=1.43l

Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  232.0)
Upper Right (  217.0,    0.0)
Lower Right (  217.0,  232.0)
Center      (  108.5,  116.0)
Band 1 Block=217x37 Type=Byte, ColorInterp=Palette
  Min=0.000 Max=255.000
  Minimum=0.000, Maximum=255.000, Mean=116.873, StdDev=70.651
  Metadata:
    STATISTICS_MAXIMUM=255
    STATISTICS_MEAN=116.87251708247
    STATISTICS_MINIMUM=0
    STATISTICS_STDDEV=70.651436875893
  Color Table (RGB with 256 entries)
    0: 0,6,21,255
    1: 0,13,98,255
...
255: 67,142,93,255

Thanks,

Agus