ABSTRACT
Purpose: Analysis of the abnormal motion of thoraco-abdominal organs in respiratory disorders such as the Thoracic Insufficiency Syndrome (TIS) and scoliosis such as adolescent idiopathic scoliosis (AIS) or early onset scoliosis (EOS) can lead to better surgical plans. We can use healthy subjects to find out the normal architecture and motion of a rib cage and associated organs and attempt to modify the patient's deformed anatomy to match to it. Dynamic magnetic resonance imaging (dMRI) is a practical and preferred imaging modality for capturing dynamic images of healthy pediatric subjects. In this paper, we propose an auto-segmentation set-up for the lungs, kidneys, liver, spleen, and thoraco-abdominal skin in these dMRI images which have their own challenges such as poor contrast, image non-standardness, and similarity in texture amongst gas, bone, and connective tissue at several inter-object interfaces. Methods: The segmentation set-up has been implemented in two steps: recognition and delineation using two deep neural network (DL) architectures (say DL-R and DL-D) for the recognition step and delineation step, respectively. The encoder-decoder framework in DL-D utilizes features at four different resolution levels to counter the challenges involved in the segmentation. We have evaluated on dMRI sagittal acquisitions of 189 (near-)normal subjects. The spatial resolution in all dMRI acquisitions is 1.46 mm in a sagittal slice and 6.00 mm between sagittal slices. We utilized images of 89 (10) subjects at end inspiration for training (validation). For testing we experimented with three scenarios: utilizing (1) the images of 90 (=189-89-10) different (remaining) subjects at end inspiration for testing, (2) the images of the aforementioned 90 subjects at end expiration for testing, and (3) the images of the aforesaid 99 (=89+10) subjects but at end expiration for testing. In some situations, we can take advantage of already available ground truth (GT) of a subject at a particular respiratory phase to automatically segment the object in the image of the same subject at a different respiratory phase and then refining the segmentation to create the final GT. We anticipate that this process of creating GT would require minimal post hoc correction. In this spirit, we conducted separate experiments where we assume to have the ground truth of the test subjects at end expiration for scenario (1), end inspiration for (2), and end inspiration for (3). Results: Amongst these three scenarios of testing, for the DL-R, we achieve a best average location error (LE) of about 1 voxel for the lungs, kidneys, and spleen and 1.5 voxels for the liver and the thoraco- abdominal skin. The standard deviation (SD) of LE is about 1 or 2 voxels. For the delineation approach, we achieve an average Dice coefficient (DC) of about 0.92 to 0.94 for the lungs, 0.82 for the kidneys, 0.90 for the liver, 0.81 for the spleen, and 0.93 for the thoraco-abdominal skin. The SD of DC is lower for the lungs, liver, and the thoraco-abdominal skin, and slightly higher for the spleen and kidneys. Conclusions: Motivated by applications in surgical planning for disorders such as TIS, AIS, and EOS, we have shown an auto-segmentation system for thoraco-abdominal organs in dMRI acquisitions. This proposed setup copes with the challenges posed by low resolution, motion blur, inadequate contrast, and image intensity non-standardness quite well. We are in the process of testing its effectiveness on TIS patient dMRI data.
ABSTRACT
Purpose: Thoracic insufficiency syndrome (TIS) affects ventilatory function due to spinal and thoracic deformities limiting lung space and diaphragmatic motion. Corrective orthopedic surgery can be used to help normalize skeletal anatomy, restoring lung space and diaphragmatic motion. This study employs free-breathing dynamic MRI (dMRI) and quantifies the 3D motion of each hemi-diaphragm surface in normal and TIS patients, and evaluates effects of surgical intervention. Materials and Methods: In a retrospective study of 149 pediatric patients with TIS and 190 healthy children, we constructed 4D images from free-breathing dMRI and manually delineated the diaphragm at end-expiration (EE) and end-inspiration (EI) time points. We automatically selected 25 points uniformly on each hemi-diaphragm surface, calculated their relative velocities between EE and EI, and derived mean velocities in 13 homologous regions for each hemi-diaphragm to provide measures of regional 3D hemi-diaphragm motion. T-testing was used to compare velocity changes before and after surgery, and to velocities in healthy controls. Results: The posterior-central region of the right hemi-diaphragm exhibited the highest average velocity post-operatively. Posterior regions showed greater velocity changes after surgery in both right and left hemi-diaphragms. Surgical reduction of thoracic Cobb angle displayed a stronger correlation with changes in diaphragm velocity than reduction in lumbar Cobb angle. Following surgery, the anterior regions of the left hemi-diaphragm tended to approach a more normal state. Conclusion: Quantification of regional motion of the 3D diaphragm surface in normal subjects and TIS patients via free-breathing dMRI is feasible. Derived measurements can be assessed in comparison to normal subjects to study TIS and the effects of surgery.
