Using the latest version of EBImage and ver 8:6.6.9.7-5 of ImageMagick (debian package) several EBImage operations fail. In particular, cropping an image with: a<-readImage("xxxxxxxx") a[1:120,50:120] # fails with "incorrect number of dimensions"] I can get cropped images by doing: a[1:120,50:120,1] # display gives correct region with a red mask a[1:120,50:120,c(1,0,0)] # display gives correct region with a red mask a[1:120,50:120,c(1,0,1)] # display gives correct region with a yellow mask a[1:120,50:120,c(1,1,1)] # display gives correct region as grayscale a[1:120,50:120,c(0,1,1)] # display gives correct region with a yellow mask and so on. Has there been a change in EBImage that is not reflected in the documentation? I am trying to build histograms for each color channel to build classification models. Thanks, Michael