[R-sig-Geo] How do I set ENVI file interleave when using GDALTransientDataset?

[Roger Bivand]

On Fri, 10 Jul 2009, Roger Bivand wrote:

> On Thu, 9 Jul 2009, Jonathan Greenberg wrote:
>
>> I'm writing a function to mimic ENVI/IDL's "ENVI_SETUP_HEAD" using rgdal, 
>> and I was wondering how I set the interleave (BIP vs BSQ vs BIL)?
>
> You'll need to use the options= argument to pass through these values when 
> you instantiate the transient dataset. There is an example of its use in 
> ?writeOGR.

Sorry, ?writeGDAL, of course. The valid options are given on the GDAL 
supported formats pages for creation and/or copying - not (naturally) the 
reading options.

>
> Roger
>
>> The subsection of my code that I assume this belongs in is below:
>>
>>   driver="ENVI"
>>   d.drv <- new("GDALDriver", driver)
>>   tds.out <- new("GDALTransientDataset", driver = d.drv, rows = 
>> dims[2],cols = dims[1], bands = bands, type = type)
>>   gt <- c(offset[1] - 0.5 * cellsize[1], cellsize[1], 0.0,offset[2] + 
>> (dims[2] -0.5) * cellsize[2], 0.0, -cellsize[2])
>>   .Call("RGDAL_SetGeoTransform", tds.out, gt, PACKAGE = "rgdal")
>>   .Call("RGDAL_SetProject", tds.out, file_projection, PACKAGE = "rgdal")
>>   fn <- tempfile()
>>   saveDataset(tds.out, fn)
>> 
>> Thanks!
>> 
>> --j
>> 
>> 
>> 
>
>

-- 
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