length(find(Y)) works correctly (splitting zeros from non-zeros in this
context) for any values in Y (not only binary) which allows very easy
modification when need to add thresholding of the probability map
length(fing(Y>threshold))
very convenient for working with c1*, c2*, c3* images from SPM-Segment
ATB,
Roman
|