我正在对某些可以具有多个波段/通道的图像进行操作(参见多光谱和高光谱遥感图像)。
我们假设我们事先不知道给定图像中存在多少条带。所有频道都命名为"band1"、"band2"...等等。阅读我们使用的图像
library(sp)
library(rgdal)
image <- readGDAL("image.tif')
我们首先将其纳入data.frame
:
imdata <- image@data
在某些情况下,特定"像素"的所有波段的内容可能为零。因此,如果它是一个四波段图像,那么我们可以使用以下命令从data.frame
中删除这些行:
imData <- imData[!((imData$band1==0) & (imData$band2==0) & (imData$band3 == 0) & (imData$band4 == 0)),]
但假设我的影像有 256 个波段。如何从数据集中删除此类行?
尝试,
imData <- imData[!(rowSums(imData) == 0),]