Hi,
Yes - see the final lines of siena_diff.cc :
cout << "AREA " << count*voxel_area << " mm^2" << endl;
cout << "VOLC " << total*voxel_volume << " mm^3" << endl;
cout << "RATIO " << (total*voxel_volume) / (count*voxel_area) << "
mm" << endl; /* mean perpendicular edge motion; l in the equations */
cout << "PBVC " << (calib*30*total*voxel_volume) /
(count*voxel_area) << " %" << endl;
So you should be able to use avwstats++ to get the count of nonzero
voxels and the sum of nonzero values ("total").
Note that voxel_volume is the vol in mm of a voxel and voxel_area is
vol to the power of 2/3.
Calib is the calibration factor that siena passes into siena_diff
(see the .siena logfile).
Cheers, Steve.
On 24 Oct 2006, at 23:27, Saemann Philipp wrote:
> Hello Steve & FSL team,
>
> on some flow maps of the siena program we'd like do some manual
> corrections -
> how can the pbvc be calculated after masking out some non-zero voxels?
>
> Can I just take the mean of all non-zero voxels in the flow image?
>
> thank you for any hints,
> Philipp
>
>
> Max-Planck-Institute of Psychiatry
> NMR Research Group
> Kraepelinsr. 2-10
> 80804 Munich
> Mail: [log in to unmask]
> Phone: 0049-89-30622-413
------------------------------------------------------------------------
---
Stephen M. Smith, Professor of Biomedical Engineering
Associate Director, Oxford University FMRIB Centre
FMRIB, JR Hospital, Headington, Oxford OX3 9DU, UK
+44 (0) 1865 222726 (fax 222717)
[log in to unmask] http://www.fmrib.ox.ac.uk/~steve
------------------------------------------------------------------------
---
|