我想堆叠和写一些Landsat波段/tiff文件,以ENVI格式的BIP交织。然而,结果总是以BSQ的形式出现,即使我将带宽顺序更改为BIP。
下面是我的代码:
library(raster)
library(rgdal)
library(gdalUtils)
inbands <-list.files(pattern= "*.tif")
stk<-stack(inbands[2], inbands[3], inbands[4])
writeRaster(stk, "BIP_test", format="ENVI", bandorder='BIP') 这也不起作用
writeRaster(stk, "BIP_test", format="ENVI", options="INTERLEAVE=PIXEL", overwrite=TRUE) 如能提供任何协助,我们将不胜感激。
发布于 2015-01-09 20:55:31
我认为这应该是可行的:
writeRaster(stk, "BIP_test", format="ENVI", options="INTERLEAVE=BIP", overwrite=TRUE)根据ENVI格式上的GDAL格式页,"BIP“(而不是”象素“)是”交错“的参数。根据我对WriteRaster帮助文件的读取,bandorder='BIP'只适用于光栅包的本机文件格式。
https://stackoverflow.com/questions/27826557
复制相似问题