Novel retrospective, respiratory-gating method enables 3D, high resolution, dynamic imaging of the upper airway during tidal breathing

A Novel Robotic Bronchoscope System for Navigation and Biopsy of Pulmonary Lesions

Transbronchial biopsy sampling, as a minimally invasive method with relatively low risk, has been proved to be a promising treatment in the field of respiratory surgery. Although several robotic bronchoscopes have been developed, it remains a great challenge to balance size and flexibility, while integrating multisensors to realize navigation during complex airway networks. This paper proposes a novel robotic bronchoscope system composed by end effector with relatively small size, relevant actuation unit, and navigation system with path planning and surgical guidance capability. The main part of the end effector is machined by bidirectional groove on a nickel-titanium tube, which can realize bending, rotation, and translation 3 degrees of freedom. A prototype of the proposed robotic bronchoscope system is designed and fabricated, and its performance is tested through several experiments to verify the stiffness, flexibility, and navigation performance. The results show that the proposed system is with good environment adaptiveness, and it can become a promising biopsy method through natural cavity of the human body.

Simulating the effect of individual upper airway anatomical features on drug deposition

This study aims to systematically isolate different anatomical features of the human pharynx with the goal to investigate their independent influence on airflow dynamics and particle deposition characteristics in a geometrically realistic human airway. Specifically, the effects of the uvula, epiglottis and soft palate on drug particle deposition are studied systematically, by carefully removing each of these anatomical features from reconstructed models based on MRI data and comparing them to a benchmark realistic airway model. Computational Fluid Dynamics using established turbulence models is employed to simulate the transport of mono-dispersed particles (3 µm) in the airway at two flow-rates. The simulations suggest three findings: 1) widening the space between the oral cavity and oropharynx and where the soft palate is situated leads to the most dramatic reduction in drug deposition in the upper airway; 2) exclusion of the uvula and epiglottis: a) affects flow dynamics in the airway; b) alters regional deposition behaviour; c) does not significantly affect the total number of particles deposited in the pharynx; and 3) the space adjacent to the soft palate is a key determinant for aerosol deposition in the extrathoracic region and is related to mechanisms of flow acceleration, diversion and recirculation.

Automatic upper airway segmentation in static and dynamic MRI via anatomy-guided convolutional neural networks

PURPOSE

Upper airway segmentation on MR images is a prerequisite step for quantitatively studying the anatomical structure and function of the upper airway and surrounding tissues. However, the complex variability of intensity and shape of anatomical structures and different modes of image acquisition commonly used in this application makes automatic upper airway segmentation challenging. In this paper, we develop and test a comprehensive deep learning-based segmentation system for use on MR images to address this problem.

MATERIALS AND METHODS

In our study, both static and dynamic MRI data sets are utilized, including 58 axial static 3D MRI studies, 22 mid-retropalatal dynamic 2D MRI studies, 21 mid-retroglossal dynamic 2D MRI studies, 36 mid-sagittal dynamic 2D MRI studies, and 23 isotropic dynamic 3D MRI studies, involving a total of 160 subjects and over 20 000 MRI slices. Samples of static and 2D dynamic MRI data sets were randomly divided into training, validation, and test sets by an approximate ratio of 5:2:3. Considering that the variability of annotation data among 3D dynamic MRIs was greater than for other MRI data sets, we increased the ratio of training data for these data to improve the robustness of the model. We designed a unified framework consisting of the following procedures. For static MRI, a generalized region-of-interest (GROI) strategy is applied to localize the partitions of nasal cavity and other portions of upper airway in axial data sets as two separate subobjects. Subsequently, the two subobjects are segmented by two separate 2D U-Nets. The two segmentation results are combined as the whole upper airway structure. The GROI strategy is also applied to other MRI modes. To minimize false-positive and false-negative rates in the segmentation results, we employed a novel loss function based explicitly on these rates to train the segmentation networks. An inter-reader study is conducted to test the performance of our system in comparison to human variability in ground truth (GT) segmentation of these challenging structures.

RESULTS

The proposed approach yielded mean Dice coefficients of 0.84±0.03, 0.89±0.13, 0.84±0.07, and 0.86±0.05 for static 3D MRI, mid-retropalatal/mid-retroglossal 2D dynamic MRI, mid-sagittal 2D dynamic MRI, and isotropic dynamic 3D MRI, respectively. The quantitative results show excellent agreement with manual delineation results. The inter-reader study results demonstrate that the segmentation performance of our approach is statistically indistinguishable from manual segmentations considering the inter-reader variability in GT.

CONCLUSIONS

The proposed method can be utilized for routine upper airway segmentation from static and dynamic MR images with high accuracy and efficiency. The proposed approach has the potential to be employed in other dynamic MRI-related applications, such as lung or heart segmentation.

A minimally interactive method for labeling respiratory phases in free-breathing thoracic dynamic MRI for constructing 4D images

OBJECTIVE

Determination of end-expiration (EE) and end-inspiration (EI) time points in the respiratory cycle in free-breathing slice image acquisitions of the thorax is one key step needed for 4D image construction via dynamic magnetic resonance imaging. The purpose of this paper is to realize the automation of the labeling process.

METHODS

The diaphragm is used as a surrogate for tracking respiratory motion and determining the state of breathing. Regions of interest (ROIs) containing the hemi-diaphragms are set by human interaction to compute the optical flow matrix between two adjacent 2D time slices. Subsequently, our approach examines the diaphragm speed and direction and by considering the change in the optical flow matrix, the EE or EI points are detected.

RESULTS AND CONCLUSION

The labeling accuracy for the lateral aspect of the left lung and the lateral aspect of the right lung (0.630.71) is significantly lower (P < 0.05) than the accuracy for other positions (0.420.44), but the error in almost all scenarios is less than 1 time point. By comparing between automatic and manual labeling in 12 scenarios, we found out that 9 scenarios showed no significant difference (P > 0.05) between two methods. Overall, our method is found to be highly agreeable with manual labeling and greatly shortens the labeling time, requiring less than 8 minutes/ study compared to 4 hours/ study for manual labeling.

SIGNIFICANCE

Our method achieves automatic labeling of EE and EI points without the need for use of patient internal or external markers.

Segmentation of 4D images via space-time neural networks

Medical imaging techniques currently produce 4D images that portray the dynamic behaviors and phenomena associated with internal structures. The segmentation of 4D images poses challenges different from those arising in segmenting 3D static images due to different patterns of variation of object shape and appearance in the space and time dimensions. In this paper, different network models are designed to learn the pattern of slice-to-slice change in the space and time dimensions independently. The two models then allow a gamut of strategies to actually segment the 4D image, such as segmentation following just the space or time dimension only, or following first the space dimension for one time instance and then following all time instances, or vice versa, etc. This paper investigates these strategies in the context of the obstructive sleep apnea (OSA) application and presents a unified deep learning framework to segment 4D images. Because of the sparse tubular nature of the upper airway and the surrounding low-contrast structures, inadequate contrast resolution obtainable in the magnetic resonance (MR) images leaves many challenges for effective segmentation of the dynamic airway in 4D MR images. Given that these upper airway structures are sparse, a Dice coefficient (DC) of ~0.88 for their segmentation based on our preferred strategy is similar to a DC of >0.95 for large non-sparse objects like liver, lungs, etc., constituting excellent accuracy.

Thoracic Quantitative Dynamic MRI to Understand Developmental Changes in Normal Ventilatory Dynamics

BACKGROUND

A database of normative quantitative measures of regional thoracic ventilatory dynamics, which is essential to understanding better thoracic growth and function in children, does not exist.

RESEARCH QUESTION

How to quantify changes in the components of ventilatory pump dynamics during childhood via thoracic quantitative dynamic MRI (QdMRI)?

STUDY DESIGN AND METHODS

Volumetric parameters were derived via 51 dynamic MRI scans for left and right lungs, hemidiaphragms, and hemichest walls during tidal breathing. Volume-based symmetry and functional coefficients were defined to compare left and right sides and to compare contributions of the hemidiaphragms and hemichest walls with tidal volumes (TVs). Statistical analyses were performed to compare volume components among four age-based groups.

RESULTS

Right thoracic components were significantly larger than left thoracic components, with average ratios of 1.56 (95% CI, 1.41-1.70) for lung TV, 1.81 (95% CI, 1.60-2.03) for hemidiaphragm excursion TV, and 1.34 (95% CI, 1.21-1.47) for hemichest wall excursion TV. Right and left lung volumes at end-expiration showed, respectively, a 44% and 48% increase from group 2 (8 ≤ age < 10) to group 3 (10 ≤ age < 12). These numbers from group 3 to group 4 (12 ≤ age ≤ 14) were 24% and 28%, respectively. Right and left hemichest wall TVs exhibited, respectively, 48% and 45% increases from group 3 to group 4.

INTERPRETATION

Normal right and left ventilatory volume components have considerable asymmetry in morphologic features and dynamics and change with age. Chest wall and diaphragm contributions vary in a likewise manner. Thoracic QdMRI can provide quantitative data to characterize the regional function and growth of the thorax as it relates to ventilation.

Quantitative Dynamic Thoracic MRI: Application to Thoracic Insufficiency Syndrome in Pediatric Patients

Background Available methods to quantify regional dynamic thoracic function in thoracic insufficiency syndrome (TIS) are limited. Purpose To evaluate the use of quantitative dynamic MRI to depict changes in regional dynamic thoracic function before and after surgical correction of TIS. Materials and Methods Images from free-breathing dynamic MRI in pediatric patients with TIS (July 2009-August 2015) were retrospectively evaluated before and after surgical correction by using vertical expandable prosthetic titanium rib (VEPTR). Eleven volumetric parameters were derived from lung, chest wall, and diaphragm segmentations, and parameter changes before versus after operation were correlated with changes in clinical parameters. Paired analysis from Student t test on MRI parameters and clinical parameters was performed to detect if changes (from preoperative to postoperative condition) were statistically significant. Results Left and right lung volumes at end inspiration and end expiration increased substantially after operation in pediatric patients with thoracic insufficiency syndrome, especially right lung volume with 22.9% and 26.3% volume increase at end expiration (P = .001) and end inspiration (P = .002), respectively. The average lung tidal volumes increased after operation for TIS; there was a 43.8% and 55.3% increase for left lung tidal volume and right lung tidal volume (P < .001 for both), respectively. However, clinical parameters did not show significant changes from pre- to posttreatment states. Thoracic and lumbar Cobb angle were poor predictors of MRI tidal volumes (chest wall, diaphragm, and left and right separately), but assisted ventilation rating and forced vital capacity showed moderate correlations with tidal volumes (chest wall, diaphragm, and left and right separately). Conclusion Vertical expandable prosthetic titanium rib operation was associated with postoperative increases in all components of tidal volume (left and right chest wall and diaphragm, and left and right lung tidal volumes) measured at MRI. Clinical parameters did not demonstrate improvements in postoperative tidal volumes. © RSNA, 2019 Online supplemental material is available for this article. See also the editorial by Paltiel in this issue.

Biomechanics of the soft-palate in sleep apnea patients with polycystic ovarian syndrome

Highly compliant tissue supporting the pharynx and low muscle tone enhance the possibility of upper airway occlusion in children with obstructive sleep apnea (OSA). The present study describes subject-specific computational modeling of flow-induced velopharyngeal narrowing in a female child with polycystic ovarian syndrome (PCOS) with OSA and a non-OSA control. Anatomically accurate three-dimensional geometries of the upper airway and soft-palate were reconstructed for both subjects using magnetic resonance (MR) images. A fluid-structure interaction (FSI) shape registration analysis was performed using subject-specific values of flow rate to iteratively compute the biomechanical properties of the soft-palate. The optimized shear modulus for the control was 38 percent higher than the corresponding value for the OSA patient. The proposed computational FSI model was then employed for planning surgical treatment for the apneic subject. A virtual surgery comprising of a combined adenoidectomy, palatoplasty and genioglossus advancement was performed to estimate the resulting post-operative patterns of airflow and tissue displacement. Maximum flow velocity and velopharyngeal resistance decreased by 80 percent and 66 percent respectively following surgery. Post-operative flow-induced forces on the anterior and posterior faces of the soft-palate were equilibrated and the resulting magnitude of tissue displacement was 63 percent lower compared to the pre-operative case. Results from this pilot study indicate that FSI computational modeling can be employed to characterize the mechanical properties of pharyngeal tissue and evaluate the effectiveness of various upper airway surgeries prior to their application.

Upper airway dynamic imaging during tidal breathing in awake and asleep subjects with obstructive sleep apnea and healthy controls

We used magnetic resonance imaging (MRI) to quantify change in upper airway dimension during tidal breathing in subjects with obstructive sleep apnea (OSA, N = 7) and BMI-matched healthy controls (N = 7) during both wakefulness and natural sleep. Dynamic MR images of the upper airway were obtained on a 1.5 T MR scanner in contiguous 7.5 mm-thick axial slices from the hard palate to the epiglottis along with synchronous MRI-compatible electroencephalogram and nasal/oral flow measurements. The physiologic data were retrospectively scored to identify sleep state, and synchronized with dynamic MR images. For each image, the upper airway was characterized by its area, and linear dimensions (lateral and anterior-posterior). The dynamic behavior of the upper airway was assessed by the maximum change in these parameters over the tidal breath. Mean upper airway caliber was obtained by averaging data over the tidal breath. There was no major difference in the upper airway structure between OSA and controls except for a narrower airway at the low-retropalatal/high-retroglossal level in OSA than in controls. Changes in upper airway size over the tidal breath ((maximum - minimum)/mean) were significantly larger in the OSA than in the control group in the low retropalatal/high retroglossal region during both wakefulness and sleep. In the four OSA subjects who experienced obstructive apneas during MR imaging, the site of airway collapse during sleep corresponded to the region of the upper airway where changes in caliber during awake tidal breathing were the greatest. These observations suggest a potential role for dynamic OSA imaging during wakefulness.

Quantitative dynamic MRI (QdMRI) Volumetric Analysis of Pediatric Patients with Thoracic Insufficiency Syndrome

The lack of standardizable objective diagnostic measurement techniques is a major hurdle in the assessment and treatment of pediatric patients with thoracic insufficiency syndrome (TIS). The aim of this paper is to explore quantitative dynamic MRI (QdMRI) volumetric parameters derived from thoracic dMRI in pediatric patients with TIS and the relationships between dMRI parameters and clinical measurements. 25 TIS patients treated with vertical expandable prosthetic titanium rib (VEPTR) surgery are included in this retrospective study. Left and right lungs at end-inspiration and end-expiration are segmented from constructed 4D dMRI images. Lung volumes and excursion (or tidal) volumes of the left/right chest wall and hemi-diaphragms are computed. Commonly used clinical parameters include thoracic and lumbar Cobb angles and respiratory measurements from pulmonary function testing (PFT). 200 3D lungs in total (left & right, pre-operative & post-operative, end-inspiration & end-expiration) are segmented for analysis. Our analysis indicates that change of resting breathing rate (RR) following surgery is negatively correlated with that of QdMRI parameters. Chest wall tidal volumes and hemi-diaphragm tidal volumes increase significantly following surgery. Clinical parameter RR reduced after surgical treatment with P values around 0.06 but no significant differences were found on other clinical parameters. The significant increase in post-operative tidal volumes suggests a treatment-related improvement in lung capacity. The reduction of RR following surgery shows that breathing function is improved. The QdMRI parameters may offer an objective marker set for studying TIS, which is currently lacking.

Upper Airway Areas, Volumes, and Linear Measurements Determined on Computed Tomography During Different Phases of Respiration Predict the Presence of Severe Obstructive Sleep Apnea

PURPOSE

The objective of this study was to analyze the potential of using low-dose volumetric computed tomography (CT) during different phases of respiration for identifying patients likely to have severe obstructive sleep apnea (OSA), defined as a respiratory disturbance index (RDI) higher than 30.

PATIENTS AND METHODS

A prospective study was undertaken at the Ramathibodi Hospital (Bangkok, Thailand). Patients with diagnosed OSA (N = 82) were recruited and separated into group 1 (RDI, ≤30; n = 36) and group 2 (RDI, >30; n = 46). The 2 groups were scanned by low-dose volumetric CT while they were 1) breathing quietly, 2) at the end of inspiration, and 3) at the end of expiration. Values for CT variables were obtained from linear measurements on lateral scout images during quiet breathing and from the upper airway area and volume measurements were obtained on axial cross-sections during different phases of respiration. All CT variables were compared between study groups. A logistic regression model was constructed to calculate a patient's likelihood of having an RDI higher than 30 and the predictive value of each variable and of the final model.

RESULTS

The minimum cross-sectional area (MCA) measured at the end of inspiration (cutoff point, ≤0.33 cm²) was the most predictive variable for the identification of patients likely to have an RDI higher than 30 (adjusted odds ratio [OR] = 5.50; 95% confidence interval [CI], 1.76-17.20; sensitivity, 74%; specificity, 72%,), followed by the MCA measured at the end of expiration (cutoff point, ≤0.21 cm²; adjusted OR = 3.28; 95% CI, 1.05-10.24; sensitivity, 70%; specificity, 68%).

CONCLUSION

CT scanning at the ends of inspiration and expiration helped identify patients with an RDI higher than 30 based on measurement of the MCA. Low-dose volumetric CT can be a useful tool to help the clinician rapidly identify patients with severe OSA and decide on the urgency to obtain a full-night polysomnographic study and to start treatment.

Characterization of upper airway obstruction using cine MRI in children with residual obstructive sleep apnea after adenotonsillectomy

OBJECTIVES/BACKGROUND

Tonsillectomy and adenoidectomy (T&A) lead to resolution of obstructive sleep apnea (OSA) in most children. However, OSA persists in about 25-40% of children. Cinematic magnetic resonance imaging (cine MRI) can aid the management of persistent OSA by localizing airway obstruction. We describe our experience in implementing and optimizing a cine MRI protocol by using a 3 Tesla MRI scanner, and the use of dexmedetomidine for sedation to improve reproducibility, safety, and diagnostic accuracy.

PATIENTS/METHODS

Patients aged 3-18 years who underwent cine MRI for the evaluation of persistent OSA after T&A and failed positive airway pressure (PAP) therapy were included. Clinical data and the apnea-hyponea index were compared with quantitative and qualitative estimates of airway obstruction from imaging sequences.

RESULTS

A total of 36 children were included with a mean age of 9.6 ± 4.6 (SD) years with 40% over 12 years of age. Two-thirds of them were boys. Seventeen out of 36 children (47%) had Down syndrome. Single site and multilevel obstruction were identified in 21 of 36 patients (58%) and in 12 of 36 patients (33%), respectively. All cine MRIs were performed without complications. Multiple regression analysis demonstrated that a combination of the minimum airway diameter and body mass index z-score best predicted OSA severity (P = 0.002).

CONCLUSIONS

Cine MRI is a sensitive, safe, and noninvasive modality for visualizing upper airway obstruction in children with persistent OSA after T&A. Accurate identification of obstruction can assist in surgical planning in children who fail PAP therapy.

Computed Tomography Characterization and Comparison With Polysomnography for Obstructive Sleep Apnea Evaluation

PURPOSE

We hypothesized that computed tomography (CT) combined with portable polysomnography (PSG) might better visualize anatomic data related to obstructive sleep apnea (OSA). The present study evaluated the CT findings during OSA and assessed their associations with the PSG data and patient characteristics.

PATIENTS AND METHODS

We designed a prospective cross-sectional study of patients with OSA. The patients underwent scanning during the awake state and apneic episodes. Associations of the predictor variables (ie, PSG data, respiratory disturbance index [RDI]), patient characteristics (body mass index [BMI], neck circumference [NC], and waist circumference [WC]), and outcome variables (ie, CT findings during apneic episodes) were assessed using logistic regression analysis. The CT findings during apneic episodes were categorized regarding the level of obstruction, single level (retropalatal [RP] or retroglossal [RG]) or multilevel (mixed RP and RG), degree of obstruction (partial or complete), and pattern of collapse (complete concentric collapse [CCC] or other patterns).

RESULTS

A total of 58 adult patients with OSA were scanned. The mean ± standard deviation for the RDI, BMI, NC, and WC were 41.6 ± 28.55, 27.80 ± 5.43 kg/m², 38.3 ± 4.3 cm, and 93.8 ± 13.6 cm, respectively. No variables distinguished between the presence of single- and multilevel airway obstruction in the present study. A high RDI (≥30) was associated with the presence of complete obstruction and CCC (odds ratio 6.33, 95% confidence interval 1.55 to 25.90; and odds ratio 3.77, 95% confidence interval 1.02 to 13.91, respectively) compared with those with a lesser RDI.

CONCLUSIONS

An increased RDI appears to be an important variable for predicting the presence of complete obstruction and CCC during OSA. Scanning during apneic episodes, using low-dose volumetric CT combined with portable PSG provided better anatomic and pathologic findings of OSA than did scans performed during the awake state.

Real-time multislice MRI during continuous positive airway pressure reveals upper airway response to pressure change

PURPOSE

To determine if a real-time magnetic resonance imaging (RT-MRI) method during continuous positive airway pressure (CPAP) can be used to measure neuromuscular reflex and/or passive collapsibility of the upper airway in individual obstructive sleep apnea (OSA) subjects.

MATERIALS AND METHODS

We conducted experiments on four adolescents with OSA and three healthy controls, during natural sleep and during wakefulness. Data were acquired on a clinical 3T scanner using simultaneous multislice (SMS) RT-MRI during CPAP. CPAP pressure level was alternated between therapeutic and subtherapeutic levels. Segmented airway area changes in response to rapid CPAP pressure drop and restoration were used to estimate 1) upper airway loop gain (UALG), and 2) anatomical risk factors, including fluctuation of airway area (FAA).

RESULTS

FAA significantly differed between OSA patients (2-4× larger) and healthy controls (Student's t-test, P < 0.05). UALG and FAA measurements indicate that neuromuscular reflex and passive collapsibility varied among the OSA patients, suggesting the presence of different OSA phenotypes. Measurements had high intrasubject reproducibility (intraclass correlation coefficient r > 0.7).

CONCLUSION

SMS RT-MRI during CPAP can reproducibly identify physiological traits and anatomical risk factors that are valuable in the assessment of OSA. This technique can potentially locate the most collapsible airway sites. Both UALG and FAA possess large variation among OSA patients.

LEVEL OF EVIDENCE

1 Technical Efficacy: Stage 1 J. Magn. Reson. Imaging 2017;46:1400-1408.

Retrospective 4D MR image construction from free-breathing slice Acquisitions: A novel graph-based approach

PURPOSE

Dynamic or 4D imaging of the thorax has many applications. Both prospective and retrospective respiratory gating and tracking techniques have been developed for 4D imaging via CT and MRI. For pediatric imaging, due to radiation concerns, MRI becomes the de facto modality of choice. In thoracic insufficiency syndrome (TIS), patients often suffer from extreme malformations of the chest wall, diaphragm, and/or spine with inability of the thorax to support normal respiration or lung growth (Campbell et al., 2003, Campbell and Smith, 2007), as such patient cooperation needed by some of the gating and tracking techniques are difficult to realize without causing patient discomfort and interference with the breathing mechanism itself. Therefore (ventilator-supported) free-breathing MRI acquisition is currently the best choice for imaging these patients. This, however, raises a question of how to create a consistent 4D image from such acquisitions. This paper presents a novel graph-based technique for compiling the best 4D image volume representing the thorax over one respiratory cycle from slice images acquired during unencumbered natural tidal-breathing of pediatric TIS patients.

METHODS

In our approach, for each coronal (or sagittal) slice position, images are acquired at a rate of about 200-300ms/slice over several natural breathing cycles which yields over 2000 slices. A weighted graph is formed where each acquired slice constitutes a node and the weight of the arc between two nodes defines the degree of contiguity in space and time of the two slices. For each respiratory phase, an optimal 3D spatial image is constructed by finding the best path in the graph in the spatial direction. The set of all such 3D images for a given respiratory cycle constitutes a 4D image. Subsequently, the best 4D image among all such constructed images is found over all imaged respiratory cycles. Two types of evaluation studies are carried out to understand the behavior of this algorithm and in comparison to a method called Random Stacking - a 4D phantom study and 10 4D MRI acquisitions from TIS patients and normal subjects. The 4D phantom was constructed by 3D printing the pleural spaces of an adult thorax, which were segmented in a breath-held MRI acquisition.

RESULTS

Qualitative visual inspection via cine display of the slices in space and time and in 3D rendered form showed smooth variation for all data sets constructed by the proposed method. Quantitative evaluation was carried out to measure spatial and temporal contiguity of the slices via segmented pleural spaces. The optimal method showed smooth variation of the pleural space as compared to Random Stacking whose behavior was erratic. The volumes of the pleural spaces at the respiratory phase corresponding to end inspiration and end expiration were compared to volumes obtained from breath-hold acquisitions at roughly the same phase. The mean difference was found to be roughly 3%.

CONCLUSIONS

The proposed method is purely image-based and post-hoc and does not need breath holding or external surrogates or instruments to record respiratory motion or tidal volume. This is important and practically warranted for pediatric patients. The constructed 4D images portray spatial and temporal smoothness that should be expected in a consistent 4D volume. We believe that the method can be routinely used for thoracic 4D imaging.

Computational fluid dynamics upper airway effective compliance, critical closing pressure, and obstructive sleep apnea severity in obese adolescent girls

Obstructive sleep apnea syndrome (OSAS) is associated with anatomical abnormalities restricting upper airway size and functional factors decreasing pharyngeal dilator activity in sleep. In this study we hypothesized that OSAS is also associated with altered pharyngeal mechanical compliance during wakefulness. Five OSAS and six control obese girls between 14 and 18 years of age were studied. All underwent polysomnography, critical closing pressure (P_crit) studies, and dynamic MRI of the upper airway during awake tidal breathing. Effective airway compliance was defined as the slope of cross-sectional area vs. average pressure between maximum inspiration and maximum expiration along the pharyngeal airway. Pharyngeal pressure fields were calculated by using image-based computational fluid dynamics and nasal resistance. Spearman correlations were calculated to test associations between apnea-hypopnea index (AHI), P_crit, and airway compliance. Effective compliances in the nasopharynx (C_NP) and velopharynx (C_VP) were lower and negative in OSAS compared with controls: -4.4 vs. 1.9 (mm²/cmH₂O, P = 0.012) and -2.1 vs. 3.9 (mm²/cmH₂O, P = 0.021), respectively, suggesting a strong phasic pharyngeal dilator activity during inspiration in OSAS compared with controls. For all subjects, C_NP and AHI correlated negatively (r_S = -0.69, P = 0.02), and passive P_crit correlated with C_NP (r_S = -0.76, P = 0.006) and with AHI (r_S = 0.86, P = 0.0006). Pharyngeal mechanics obtained during wakefulness could be used to characterize subjects with OSAS. Moreover, negative effective compliance during wakefulness and its correlation to AHI and P_crit suggest that phasic dilator activity of the upper pharynx compensates for negative pressure loads in these subjects.

Evaluation of upper airway collapsibility using real-time MRI

PURPOSE

To develop and demonstrate a real-time MRI method for assessing upper airway collapsibility in sleep apnea.

MATERIALS AND METHODS

Data were acquired on a clinical 3 Tesla scanner using a radial CAIPIRIHNA sequence with modified golden angle view ordering and reconstructed using parallel imaging and compressed sensing with temporal finite difference sparsity constraint. Segmented airway areas together with synchronized facemask pressure were used to calculate airway compliance and projected closing pressure, Pclose , at four axial locations along the upper airway. This technique was demonstrated in five adolescent obstructive sleep apnea (OSA) patients, three adult OSA patients and four healthy volunteers. Heart rate, oxygen saturation, facemask pressure, and abdominal/chest movements were monitored in real-time during the experiments to determine sleep/wakefulness.

RESULTS

Student's t-tests showed that both compliance and Pclose were significantly different between healthy controls and OSA patients (P < 0.001). The results also suggested that a narrower airway site does not always correspond to higher collapsibility.

CONCLUSION

With the proposed methods, both compliance and Pclose can be calculated and used to quantify airway collapsibility in OSA with an awake scan of 30 min total scan room time. J. Magn. Reson. Imaging 2016;44:158-167.

Supra-Epiglottic Upper Airway Volume in Elderly Patients with Obstructive Sleep Apnea Hypopnea Syndrome

Objective

Small upper airway measurements areas and high body mass index are recognized risk factors for obstructive sleep apnea syndrome (OSAS) in non-elderly populations; however, there is limited information regarding elderly patients. We evaluated whether upper airway volume is associated with OSAS and OSAS treated with continuous positive airway pressure (CPAP) treatment and whether BMI is correlated with upper airway volume and measurements in elderly subjects.

Methods

In 60 volunteers aged 75.58±0.9 years: 20 OSAS, 20 OSAS chronically treated with CPAP, and 20 controls, semi-automatic segmentation, retropalatal distance and transverse diameter of the supra-epiglottic upper airway were evaluated using 3DT1-weighted magnetic resonance imaging. Anteroposterior to transverse diameter ratio was defined as retropalatar diameter/transverse diameter.

Results

There were no significant differences in supra-epiglottic upper airway volume between OSAS, CPAP treated patients, and controls. There were significant differences in retropalatal distance and anteroposterior to transverse diameter ratio between OSAS, CPAP treated patients, and controls (P = 0.008 and P<0.0001 respectively). There was a significant correlation between body mass index and retropalatal distance (P<0.05) but not with supra-epiglottic upper airway volume.

Conclusion

In elderly subjects, OSAS and body mass index are not associated with changes in supra-epiglottic upper airway volume but are associated with modification of pharynx shape.

Minimally interactive segmentation of 4D dynamic upper airway MR images via fuzzy connectedness

PURPOSE

There are several disease conditions that lead to upper airway restrictive disorders. In the study of these conditions, it is important to take into account the dynamic nature of the upper airway. Currently, dynamic magnetic resonance imaging is the modality of choice for studying these diseases. Unfortunately, the contrast resolution obtainable in the images poses many challenges for an effective segmentation of the upper airway structures. No viable methods have been developed to date to solve this problem. In this paper, the authors demonstrate a practical solution by employing an iterative relative fuzzy connectedness delineation algorithm as a tool.

METHODS

3D dynamic images were collected at ten equally spaced instances over the respiratory cycle (i.e., 4D) in 20 female subjects with obstructive sleep apnea syndrome. The proposed segmentation approach consists of the following steps. First, image background nonuniformities are corrected which is then followed by a process to correct for the nonstandardness of MR image intensities. Next, standardized image intensity statistics are gathered for the nasopharynx and oropharynx portions of the upper airway as well as the surrounding soft tissue structures including air outside the body region, hard palate, soft palate, tongue, and other soft structures around the airway including tonsils (left and right) and adenoid. The affinity functions needed for fuzzy connectedness computation are derived based on these tissue intensity statistics. In the next step, seeds for fuzzy connectedness computation are specified for the airway and the background tissue components. Seed specification is needed in only the 3D image corresponding to the first time instance of the 4D volume; from this information, the 3D volume corresponding to the first time point is segmented. Seeds are automatically generated for the next time point from the segmentation of the 3D volume corresponding to the previous time point, and the process continues and runs without human interaction and completes in 10 s for segmenting the airway structure in the whole 4D volume.

RESULTS

Qualitative evaluations performed to examine smoothness and continuity of motions of the entire upper airway as well as its transverse sections at critical anatomic locations indicate that the segmentations are consistent. Quantitative evaluations of the separate 200 3D volumes and the 20 4D volumes yielded true positive and false positive volume fractions around 95% and 0.1%, respectively, and mean boundary placement errors under 0.5 mm. The method is robust to variations in the subjective action of seed specification. Compared with a segmentation approach based on a registration technique to propagate segmentations, the proposed method is more efficient, accurate, and less prone to error propagation from one respiratory time point to the next.

CONCLUSIONS

The proposed method is the first demonstration of a viable and practical approach for segmenting the upper airway structures in dynamic MR images. Compared to registration-based methods, it effectively reduces error propagation and consequently achieves not only more accurate segmentations but also more consistent motion representation in the segmentations. The method is practical, requiring minimal user interaction and computational time.

Minimum Field Strength Simulator for Proton Density Weighted MRI

OBJECTIVE

To develop and evaluate a framework for simulating low-field proton-density weighted MRI acquisitions based on high-field acquisitions, which could be used to predict the minimum B0 field strength requirements for MRI techniques. This framework would be particularly useful in the evaluation of de-noising and constrained reconstruction techniques.

MATERIALS AND METHODS

Given MRI raw data, lower field MRI acquisitions can be simulated based on the signal and noise scaling with field strength. Certain assumptions are imposed for the simulation and their validity is discussed. A validation experiment was performed using a standard resolution phantom imaged at 0.35 T, 1.5 T, 3 T, and 7 T. This framework was then applied to two sample proton-density weighted MRI applications that demonstrated estimation of minimum field strength requirements: real-time upper airway imaging and liver proton-density fat fraction measurement.

RESULTS

The phantom experiment showed good agreement between simulated and measured images. The SNR difference between simulated and measured was ≤ 8% for the 1.5T, 3T, and 7T cases which utilized scanners with the same geometry and from the same vendor. The measured SNR at 0.35T was 1.8- to 2.5-fold less than predicted likely due to unaccounted differences in the RF receive chain. The predicted minimum field strength requirements for the two sample applications were 0.2 T and 0.3 T, respectively.

CONCLUSIONS

Under certain assumptions, low-field MRI acquisitions can be simulated from high-field MRI data. This enables prediction of the minimum field strength requirements for a broad range of MRI techniques.

Dynamic 3-D MR Visualization and Detection of Upper Airway Obstruction During Sleep Using Region-Growing Segmentation

GOAL

We demonstrate a novel and robust approach for visualization of upper airway dynamics and detection of obstructive events from dynamic 3-D magnetic resonance imaging (MRI) scans of the pharyngeal airway.

METHODS

This approach uses 3-D region growing, where the operator selects a region of interest that includes the pharyngeal airway, places two seeds in the patent airway, and determines a threshold for the first frame.

RESULTS

This approach required 5 s/frame of CPU time compared to 10 min/frame of operator time for manual segmentation. It compared well with manual segmentation, resulting in Dice Coefficients of 0.84 to 0.94, whereas the Dice Coefficients for two manual segmentations by the same observer were 0.89 to 0.97. It was also able to automatically detect 83% of collapse events.

CONCLUSION

Use of this simple semiautomated segmentation approach improves the workflow of novel dynamic MRI studies of the pharyngeal airway and enables visualization and detection of obstructive events.

SIGNIFICANCE

Obstructive sleep apnea (OSA) is a significant public health issue affecting 4-9% of adults and 2% of children. Recently, 3-D dynamic MRI of the upper airway has been demonstrated during natural sleep, with sufficient spatiotemporal resolution to noninvasively study patterns of airway obstruction in young adults with OSA. This study makes it practical to analyze these long scans and visualize important factors in an MRI sleep study, such as the time, site, and extent of airway collapse.

Pediatric sleep-related breathing disorders: advances in imaging and computational modeling

We understand now that sleep of sufficient length and quality is required for good health. This is particularly true for infants and children, who have the added physiologic task of growth and development, as compared to their adult counterparts. Sleep-related breathing disorders (SRBDs) are common in childhood and if unrecognized and not treated can result in significant morbidity. For example, children with obstructive sleep apnea (OSA) can exhibit behavioral, mood, and learning difficulties. If left untreated, alterations in the function of the autonomic nervous system and a chronic inflammatory state result, contributing to the risk of heart disease, stroke, glucose intolerance, and hypertension in adulthood.

Seeing sleep: dynamic imaging of upper airway collapse and collapsibility in children

Sleep disordered breathing in children ranges from snoring, which has a prevalence of 12%, to obstructive sleep apnea (OSA) syndrome, which has a prevalence of 2?3% in the general population [1]. The underlying causes of pediatric OSA are extremely complex. There are bony structural influences, as seen in craniofacial abnormalities, and soft tissue abnormalities, such as a large tongue, redundant soft tissue, or compliance/collapsibility issues. In some groups, such as those with Down syndrome, a combination of these factors comes into play.

Body-wide hierarchical fuzzy modeling, recognition, and delineation of anatomy in medical images

To make Quantitative Radiology (QR) a reality in radiological practice, computerized body-wide Automatic Anatomy Recognition (AAR) becomes essential. With the goal of building a general AAR system that is not tied to any specific organ system, body region, or image modality, this paper presents an AAR methodology for localizing and delineating all major organs in different body regions based on fuzzy modeling ideas and a tight integration of fuzzy models with an Iterative Relative Fuzzy Connectedness (IRFC) delineation algorithm. The methodology consists of five main steps: (a) gathering image data for both building models and testing the AAR algorithms from patient image sets existing in our health system; (b) formulating precise definitions of each body region and organ and delineating them following these definitions; (c) building hierarchical fuzzy anatomy models of organs for each body region; (d) recognizing and locating organs in given images by employing the hierarchical models; and (e) delineating the organs following the hierarchy. In Step (c), we explicitly encode object size and positional relationships into the hierarchy and subsequently exploit this information in object recognition in Step (d) and delineation in Step (e). Modality-independent and dependent aspects are carefully separated in model encoding. At the model building stage, a learning process is carried out for rehearsing an optimal threshold-based object recognition method. The recognition process in Step (d) starts from large, well-defined objects and proceeds down the hierarchy in a global to local manner. A fuzzy model-based version of the IRFC algorithm is created by naturally integrating the fuzzy model constraints into the delineation algorithm. The AAR system is tested on three body regions - thorax (on CT), abdomen (on CT and MRI), and neck (on MRI and CT) - involving a total of over 35 organs and 130 data sets (the total used for model building and testing). The training and testing data sets are divided into equal size in all cases except for the neck. Overall the AAR method achieves a mean accuracy of about 2 voxels in localizing non-sparse blob-like objects and most sparse tubular objects. The delineation accuracy in terms of mean false positive and negative volume fractions is 2% and 8%, respectively, for non-sparse objects, and 5% and 15%, respectively, for sparse objects. The two object groups achieve mean boundary distance relative to ground truth of 0.9 and 1.5 voxels, respectively. Some sparse objects - venous system (in the thorax on CT), inferior vena cava (in the abdomen on CT), and mandible and naso-pharynx (in neck on MRI, but not on CT) - pose challenges at all levels, leading to poor recognition and/or delineation results. The AAR method fares quite favorably when compared with methods from the recent literature for liver, kidneys, and spleen on CT images. We conclude that separation of modality-independent from dependent aspects, organization of objects in a hierarchy, encoding of object relationship information explicitly into the hierarchy, optimal threshold-based recognition learning, and fuzzy model-based IRFC are effective concepts which allowed us to demonstrate the feasibility of a general AAR system that works in different body regions on a variety of organs and on different modalities.

Real-time 3D magnetic resonance imaging of the pharyngeal airway in sleep apnea

PURPOSE

To investigate the feasibility of real-time 3D magnetic resonance imaging (MRI) with simultaneous recording of physiological signals for identifying sites of airway obstruction during natural sleep in pediatric patients with sleep-disordered breathing.

METHODS

Experiments were performed using a three-dimensional Fourier transformation (3DFT) gradient echo sequence with prospective undersampling based on golden-angle radial spokes, and L1-norm regularized iterative self-consistent parallel imaging (L1-SPIRiT) reconstruction. This technique was demonstrated in three healthy adult volunteers and five pediatric patients with sleep-disordered breathing. External airway occlusion was used to induce partial collapse of the upper airway on inspiration and test the effectiveness of the proposed imaging method. Apneic events were identified using information available from synchronized recording of mask pressure and respiratory effort.

RESULTS

Acceptable image quality was obtained in seven of eight subjects. Temporary airway collapse induced via inspiratory loading was successfully imaged in all three volunteers, with average airway volume reductions of 63.3%, 52.5%, and 33.7%. Central apneic events and associated airway narrowing/closure were identified in two pediatric patients. During central apneic events, airway obstruction was observed in the retropalatal region in one pediatric patient.

CONCLUSION

Real-time 3D MRI of the pharyngeal airway with synchronized recording of physiological signals is feasible and may provide valuable information about the sites and nature of airway narrowing/collapse during natural sleep.