我正在对某些可能具有多个波段/通道的图像进行操作(参见多光谱和高光谱遥感图像)。
让我们假设我们事先不知道给定图像中存在多少个频带。所有的频道都被命名为"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个波段。如何从数据集中删除这些行?
发布于 2015-01-06 19:40:13
试一下,
imData <- imData[!(rowSums(imData) == 0),]https://stackoverflow.com/questions/27797692
复制相似问题