Total intra-cranial volume

How do I calculate Total intra-cranial volume?

Answer 1

[Walter Roberson]

Generally speaking, you start by segmenting the volume of data to figure out what is inside the cranium or not, creating a volume of logical values (inside --> true) .

Once that is done, nnz(VoxelIsInside) * voxel_height * voxel_width * voxel_depth gives you the total intra-cranial volume.

Voxels are often cubical, but depending on the modality might not be. For example at least at one time it was common for MRI to use a Z spacing that is twice as long as the X or Y spacing.

You should be able to figure out the voxel dimensions by examining DICOM headers -- but the details of the method can depend upon modality.

[BUHARI IBRAHIM]

Thank you for the prompt response.

I want to generate txt file from matlab that can be used for global calculation. How can i do it?

[Walter Roberson]

See writelines writematrix writecell writetable fprintf dlmwrite