【发布时间】：2015-01-06 19:40:12
【问题描述】：
我正在处理可以具有多个波段/通道的某些图像（参见多光谱和高光谱遥感图像）。
假设我们事先不知道给定图像中存在多少条带。所有通道都被命名为“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[!Reduce('&' , imData==0),]可能还有!!rowSums(imData!=0) -
我猜第一个是
imData[!Reduce(&,as.data.frame(imData==0))),]