Image segmentation of trabecular spongiosa by visual inspection of the gradient magnitude

Recent advances in physical models of skeletal dosimetry utilize high-resolution 3-dimensional microscopic computed tomography images of trabecular spongiosa. These images are coupled to radiation transport codes to assess energy deposition within active bone marrow and trabecular endosteum. These transport codes rely primarily on the segmentation of the spongiosa images into bone and marrow voxels. Image thresholding has been the segmentation of choice for bone sample images because of its extreme simplicity. However, the ability of the segmentation to reproduce the physical boundary between bone and marrow depends on the selection of the threshold value. Statistical models, as well as visual inspection of the image, have been employed extensively to determine the correct threshold. Both techniques are affected by partial volume effect and can provide unexpected results if performed without care. In this study, we propose a new technique to threshold trabecular spongiosa images based on visual inspection of the image gradient magnitude. We first show that the gradient magnitude of the image reaches a maximum along a surface that remains almost independent of partial volume effect and that is a good representation of the physical boundary between bone and marrow. A computer program was then developed to allow a user to compare the position of the iso-surface produced by a threshold with the gradient magnitude. The threshold that produces the iso-surface that best coincides with the maximum gradient is chosen. The technique was finally tested with a set of images of a true bone sample with different resolutions, as well as with three images of a cube of Duocell aluminium foam of known be the case and the method described here should be applied cautiously. With this in mind, the technique is well suited for trabecular spongiosa segmentation and will be used for future studies of skeletal dosimetry using microCT scanning.