## Create an MSImageSet object
spectra <- matrix(1:27, nrow=3)
mz <- 101:103
coord <- expand.grid(x=1:3, y=1:3)
msset <- MSImageSet(spectra=spectra, mz=mz, coord=coord)
# Find pixel indices
pixels(msset, x=2, y=2)
pixels(msset, coord=list(x=2, y=2))
pixels(msset, coord=list(x=c(2,3), y=c(2,3)))
# Find feature indices
features(msset, mz=102)
features(msset, mz=c(101,103))
features(msset, mz=c(102.2))
Run the code above in your browser using DataLab