I’m looking for a simple option to extract paraspinal muscle cross-sectional area (CSA) and its centroid location from MRI/CT scans (horizontal vertebral mid-planes). Is there already a feature/an extension available for this?
You can segment muscles using Segment Editor module and then use “Segment Cross-Section Area” module (provided by “Sandbox” extension in “Examples” category in Extension manager):
I have a follow-up question on this. I am writing a tutorial for this, and I am wondering if the requirement is that we should have isotropic data for correct estimation. Here is why I am asking this:
If I paint a 10mm diameter circle on the MRhead and calculate the area I get 243mm2, which is incorrect. If I resample the MRhead to be isotropic via CropVolume and paint the same area with the same diameter in a new segmentation, I get 309mm2 which is within measurement error. You get approximately 243mm2, if you scale the correct measurement by the slice spacing (1.3mm) which is different from XY spacing (1mm) in this dataset.
I attached the screenshots of the both segments (red from isotropic MRhead, and yellow is the original data).
I also suggest replacing row, column, and slice wording with the usual axial, coronal, sagittal.
Good catch. Isotropic spacing was not a requirement at all, this was just a bug that the automatic test did not reveal as it only checked computation along one axis. The issue is fixed now.
Measurement currently does not happen along anatomical axis but along an image IJK axis, so we cannot use anatomical terms. Row, column, slice can be mapped to any anatomical axis or oblique directions. We could change the axis names from row, column, slice to I, J, K but I don’t think it would be much better.
Another follow-up question just to make sure… Do the values in column B in the CSA output table (see below) represent the location of the centroid of the respective CSA?
You can hover the moues pointer over any column and wait for a few seconds without moving it then a tooltip appears that tells what that column contains. “Position” column contain slice center positions in RAS coordinate system.
You can get detailed statistics about each segment, such as centroid, oriented bounding box, etc. using Segment Statistics module.