[MINC-users] mincstats histogram count limit?

[Alex Zijdenbos]

Hi Vladimir,

print_all_labels gives the same (correct) result as mincstats -count
-binvalue label,label,.... (green values). It only seems to be the counts
in minstats' histogram that are off.

By the way, these counts are printed in the histogram using %12g - that
shouldn't cap them at 1.67772e+07 I would think, but either way they should
probably be printed as integers.

-- A


On Fri, Nov 15, 2013 at 11:44 PM, Vladimir S. FONOV <
vladimir.fonov at gmail.com> wrote:

> Hello,
>
>
> what about print_all_labels ?
>
> On 2013-11-15, at 11:17 PM, Alex Zijdenbos <zijdenbos at gmail.com> wrote:
>
> > 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
> > _______________________________________________
> > MINC-users at bic.mni.mcgill.ca
> > http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users
>
> ---
> Best regards,
>
> Vladimir S. FONOV ~ v.s.fonov <at> ilmarin.info
>
>
>
>
>
> _______________________________________________
> MINC-users at bic.mni.mcgill.ca
> http://www.bic.mni.mcgill.ca/mailman/listinfo/minc-users
>
>