A 有一个包含 4 个波段的 RasterLayer:
>rx<-raster("/media/karimdion/Passport/Essais/po_3804017_bgrn_0000000 tif")
> str(rx)
Formal class 'RasterLayer' [package "raster"] with 12 slots
@ file :Formal class '.RasterFile' [package "raster"] with 12 slots
@ nbands : int 4
@ bandorder : chr "BIL"
@ data :Formal class '.SingleLayerData' [package "raster"] with 13 slots
@ min : num 0
@ max : num 65535
@ band : int 1
@ history : list()
@ title : chr(0)
@ extent :Formal class 'Extent' [package "raster"] with 4 slots
@ xmin: num 655781
@ xmax: num 666701
@ ymin: num 4071522
@ ymax: num 4084598
@ rotated : logi FALSE
@ rotation:Formal class '.Rotation' [package "raster"] with 2 slots
@ geotrans: num(0)
@ transfun:function ()
@ ncols : int 5460
@ nrows : int 6538
@ crs :Formal class 'CRS' [package "sp"] with 1 slots
@ projargs: chr "+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
@ z : list()
我想将每个波段作为 RasterLayer 放在单个变量中。因此,我会得到 4 个具有 1 个波段的 RasterLayers。我怎样才能用 R 做到这一点?
文件名 f (此处使用示例文件以实现可重复性)
f <- system.file("external/rlogo.grd", package="raster")
对于特定的乐队,你可以这样做
r <- raster(f, band=2)
或者对于所有乐队
b <- brick(f)
或者(效率较低):
s <- stack(f)
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)