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: jgrn...@hotmail.com, Gchat: jgrn307
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
