[MINC-users] mincstats histogram count limit?

[Alex Zijdenbos]

Hello all,

I seem to have stumbled across a limit to the histogram counts in
mincstats. I thought it may be a data type limitation, but going through
the source code I so far haven't been able to find a cause.

I have a classified, discrete volume only containing values 1, 2, and 3;
it's a 722x866x722 volume with 0.25mm isotropic voxels. To get the voxel
counts for each of the labels in the volume:

$ mincstats -clobber -quiet -none -integer_histogram -hist_range 1 3
-histogram hr_hist.txt hr.mnc
451431944
$ cat hr_hist.txt
# histogram for: hr.mnc
#  mask file:    (null)
#  domain:       0.5  3.5
#  entropy:      1.32723
# bin centres                 counts
  1                      1.67322e+07
  2                      1.67772e+07
  3                      1.67772e+07

However:

$ mincstats -clobber -none -count -binvalue 1,2,3 hr.mnc
Included Range:    0.5   1.5
# voxels:          16732212
Included Range:    1.5   2.5
# voxels:          42817611
Included Range:    2.5   3.5
# voxels:          35033439

The red and green results _should_ be the same; but in the histogram output
the seem to be capped at 1.67772e+07, which is clearly wrong.

Anybody able to figure this out?

Thanks,

-- A