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

Jonathan Greenberg greenberg at ucdavis.edu
Fri Jul 10 03:33:40 CEST 2009


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


-- 

Jonathan A. Greenberg, PhD
Postdoctoral Scholar
Center for Spatial Technologies and Remote Sensing (CSTARS)
University of California, Davis
One Shields Avenue
The Barn, Room 250N
Davis, CA 95616
Cell: 415-794-5043
AIM: jgrn307, MSN: jgrn307 at hotmail.com, Gchat: jgrn307



More information about the R-sig-Geo mailing list