[MINC-users] mincstats histogram count limit?

[Vladimir S. FONOV]

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