[R-sig-Geo] Try to write GeoTiff with rgdal

[Roger Bivand]

On Fri, 23 May 2008, Agustin Lobo wrote:

> Many programs are not able to display tiff files with more
> than 3 bands. Test with qgis 0.10 or globalmapper. Or with matlab

Any other program using GDAL (like QGIS) will read files written by GDAL. 
ArcGIS is an exception, as it predictably does the wrong thing. My 
experience is that you have to read GeoTiff through the toolbox raster 
calculate statistics module, because ArcGIS guesses wrong about how the 
data are configured, and often ignores valid NA cells (hence the odd 
values). A three-band will typically be interpreted by ArcGIS as an RGB 
image. Doing it a band at a time, and going in through the toolbox (think 
Python script) is feasible. On the other hand, non-GDAL ENVI handles 
GDAL-written multiband GeoTiffs without any trouble that I've seen.

Welcome to the world of "interoperability"!

Roger

>
> Agus
>
>
>
> Lydia Gerharz wrote:
>>  Hi,
>>
>>  I have problems with writeGDAL. With:
>> 
>> > writeGDAL(mean2005, "mean.tif", drivername="GTiff", options="TFW=YES")
>>
>>  I export a SpatialGridDataFrame with 9 bands to a GeoTiff. But when I try
>>  to display this GeoTiff in ArcGIS, all 9 bands are not correctly displayed
>>  and all cells have the same value (1.17549e-038). I also tried another GIS
>>  (PcRaster) but it doesn't work. Can somebody tell me what I did wrong? In
>>  R all 9 bands are displayed correctly.
>>
>>  Thank you very much,
>>  Lydia
>>
>>
>>  ------------------------------------------------------------------------
>>
>>  _______________________________________________
>>  R-sig-Geo mailing list
>>  R-sig-Geo at stat.math.ethz.ch
>>  https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>

-- 
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: Roger.Bivand at nhh.no