Use of an optical flow method for the analysis of serial CT lung images

Preliminary Assessment of an Optical Flow Method (OFM) for Nonrigid Registration and Temporal Subtraction (TS) of Serial CT Examinations to Facilitate Evaluation of Interval Change in Metastatic Lung Nodules

RATIONALE AND OBJECTIVES

Accurate assessment of size change of lung nodules on chest computed tomography (CT) is important for diagnosis and response assessment. However, manual methods are time-consuming and error-prone. We therefore assessed whether an optical flow method (OFM) with temporal subtraction (TS) can facilitate detection and quantification of lung nodule change on serial CT datasets.

MATERIALS AND METHODS

Serial chest CT examinations were selected from 12 patients with multiple pulmonary metastases. Lung nodules were evaluated for change in size using: (1) OFM with TS and (2) reference standard visual and manual assessment. Average time required to assess interval change using both methods was recorded and compared. Concordance of agreement between OFM with TS and reference standard assessment for nodule change was examined.

RESULTS

285 solid pulmonary nodules were evaluated. The average time per nodule to assess interval change in nodule size by OFM with TS (mean 1.15 + 0.5 minutes) was significantly less (P = 0.02) than that the reference standard approach (mean 1.56 + 0.5 minutes). Agreement between OFM with TS and reference standard occurred for 63.2% of nodules overall (kappa = 0.50, standard error 0.35, P< 0.00001), and significantly increased with larger nodule size (kappa = 0.48 for nodules <5 mm; kappa = 0.94 for nodules >20 mm, P < 0.0001).

CONCLUSIONS

This preliminary study demonstrates the feasibility of an OFM with TS to assess for interval change in metastatic lung nodules on serial CT examinations with significantly improved reading speed and moderate agreement relative to reference standard assessment. Agreement improved with larger nodule size.

Recent Advances in Computed Tomography Imaging in Chronic Obstructive Pulmonary Disease

Lung imaging is increasingly being used to diagnose, quantify, and phenotype chronic obstructive pulmonary disease (COPD). Although spirometry is the gold standard for the diagnosis of COPD and for severity staging, the role of computed tomography (CT) imaging has expanded in both clinical practice and research. COPD is a heterogeneous disease with considerable variability in clinical features, radiographic disease, progression, and outcomes. Recent studies have examined the utility of CT imaging in enhancing diagnostic certainty, improving phenotyping, predicting disease progression and prognostication, selecting patients for intervention, and also in furthering our understanding of the complex pathophysiology of this disease. Multiple CT metrics show promise for use as imaging biomarkers in COPD.

Evaluation of Deformable Image Registration for Three-Dimensional Temporal Subtraction of Chest Computed Tomography Images

Purpose

To perform lung image registration for reducing misregistration artifacts on three-dimensional (3D) temporal subtraction of chest computed tomography (CT) images, in order to enhance temporal changes in lung lesions and evaluate these changes after deformable image registration (DIR).

Methods

In 10 cases, mutual information (MI) lung mask affine mapping combined with cross-correlation (CC) lung diffeomorphic mapping was used to implement lung volume registration. With advanced normalization tools (ANTs), we used greedy symmetric normalization (greedy SyN) as a transformation model, which involved MI-CC-SyN implementation. The resulting displacement fields were applied to warp the previous (moving) image, which was subsequently subtracted from the current (fixed) image to obtain the lung subtraction image.

Results

The average minimum and maximum log-Jacobians were 0.31 and 3.74, respectively. When considering 3D landmark distance, the root-mean-square error changed from an average of 20.82 mm for P_fixed to P_moving to 0.5 mm for P_warped to P_fixed. Clear shadows were observed as enhanced lung nodules and lesions in subtraction images. The lesion shadows showed lesion shrinkage changes over time. Lesion tissue morphology was maintained after DIR.

Conclusions

DIR (greedy SyN) effectively and accurately enhanced temporal changes in chest CT images and decreased misregistration artifacts in temporal subtraction images.

Using imaging as a biomarker for asthma

There have been significant advancements in the various imaging techniques being used for the evaluation of asthmatic patients, both from a clinical and research perspective. Imaging characteristics can be used to identify specific asthmatic phenotypes and provide a more detailed understanding of endotypes contributing to the pathophysiology of the disease. Computed tomography, magnetic resonance imaging, and positron emission tomography can be used to assess pulmonary structure and function. It has been shown that specific airway and lung density measurements using computed tomography correlate with clinical parameters, including severity of disease and pathology, but also provide unique phenotypes. Hyperpolarized ¹²⁹Xe and ³He are gases used as contrast media for magnetic resonance imaging that provide measurement of distal lung ventilation reflecting small-airway disease. Positron emission tomography can be useful to identify and target lung inflammation in asthmatic patients. Furthermore, imaging techniques can serve as a potential biomarker and be used to assess response to therapies, including newer biological treatments and bronchial thermoplasty.

Comparison between multivolume CT-based surrogates of regional ventilation in healthy subjects

RATIONALE AND OBJECTIVES

The assessment of regional ventilation is of critical importance when investigating lung function during disease progression and planning of pulmonary interventions. Recently, different computed tomography (CT)-based parameters have been proposed as surrogates of lung ventilation. The aim of the present study was to compare these parameters, namely variations of density (ΔHU), specific volume (sVol), and specific gas volume (ΔSVg) between different lung volumes, in relation to their topographic distribution within the lung.

MATERIALS AND METHODS

Ten healthy volunteers were scanned via high-resolution CT at residual volume (RV) and total lung capacity (TLC); ΔHU, sVol, and ΔSVg were mapped voxel by voxel after registering TLC onto RV. Variations of the three parameters along the vertical and horizontal directions were analyzed.

RESULTS

Along the vertical direction (from ventral to dorsal regions), a strong dependence on gravity was found in ΔHU and sVol, with greater values in the dorsal regions of the lung (P < .001), whereas ΔSVg was more homogeneously distributed within the lung. Conversely, along the caudocranial direction (from lung bases to apexes) where no gravitational gradient is present, the three parameters behaved similarly, with lower values at the apices.

CONCLUSIONS

ΔHU, sVol, and ΔSVg behave differently along the gravity direction. As the greater amount of air delivered to the dependent portion of the lung supplies a larger number of alveoli, the amount of gas delivered to alveoli compared to the mass of tissue is not gravity dependent. The minimization of gravity dependence in the distribution of ventilation when using ΔSVg suggests that this parameter is more reliable to discriminate healthy from pathologic regions.

Assessment of regional lung function with multivolume (1)H MR imaging in health and obstructive lung disease: comparison with (3)He MR imaging

PURPOSE

To introduce a method based on multivolume proton (hydrogen [(1)H]) magnetic resonance (MR) imaging for the regional assessment of lung ventilatory function, investigating its use in healthy volunteers and patients with obstructive lung disease and comparing the outcome with the outcome of the research standard helium 3 ((3)He) MR imaging.

MATERIALS AND METHODS

The institutional review board approved the HIPAA-compliant protocol, and informed written consent was obtained from each subject. Twenty-six subjects, including healthy volunteers (n = 6) and patients with severe asthma (n = 11) and mild (n = 6) and severe (n = 3) emphysema, were imaged with a 1.5-T whole-body MR unit at four lung volumes (residual volume [ RV residual volume ], functional residual capacity [ FRC functional residual capacity ], 1 L above FRC functional residual capacity [ FRC+1 L 1 L above FRC ], total lung capacity [ TLC total lung capacity ]) with breath holds of 10-11 seconds, by using volumetric interpolated breath-hold examination. Each pair of volumes were registered, resulting in maps of (1)H signal change between the two lung volumes. (3)He MR imaging was performed at FRC+1 L 1 L above FRC by using a two-dimensional gradient-echo sequence. (1)H signal change and (3)He signal were measured and compared in corresponding regions of interest selected in ventral, intermediate, and dorsal areas.

RESULTS

In all volunteers and patients combined, proton signal difference between TLC total lung capacity and RV residual volume correlated positively with (3)He signal (correlation coefficient R(2) = 0.64, P < .001). Lower (P < .001) but positive correlation results from (1)H signal difference between FRC functional residual capacity and FRC+1 L 1 L above FRC (R(2) = 0.44, P < .001). In healthy volunteers, (1)H signal changes show a higher median and interquartile range compared with patients with obstructive disease and significant differences between nondependent and dependent regions.

CONCLUSION

Findings in this study demonstrate that multivolume (1)H MR imaging, without contrast material, can be used as a biomarker for regional ventilation, both in healthy volunteers and patients with obstructive lung disease.

Extracellular matrix deposition and scaffold biodegradation in an in vitro three-dimensional model of bone by X-ray computed microtomography

The development of an in vitro model of bone and the optimization of tools for determining the biological processes occurring during bone repair remains a major goal in the field of bone tissue engineering. Recently, a model based on a three-dimensional co-culture of osteoblasts and osteoclast precursors in Skelite(TM) scaffolds was developed. Although induction of osteoblast and osteoclast differentiation was observed, a complete evaluation of bone deposition and biodegradation processes was missing due to technical limitations. In the current study, both X-ray computed microtomography and histological analysis were used to monitor these two key biological processes in the same in vitro model. Either osteoblasts or a combination of osteoblasts and osteoclasts were seeded on Skelite(TM) scaffolds. Scaffold biodegradation and increased bone deposition together with a more organized extracellular matrix were observed in the co-cultures, highlighting the role of osteoclasts in the determination and regulation of bone deposition. Results confirmed the potential and relevance of co-culturing osteoblasts and osteoclasts to resemble native tissue. The combination of X-ray computed microtomography and histology presented in this study could be useful in future studies for the validation and development of new in vitro culture systems for bone tissue engineering.

Computational analysis of thoracic multidetector row HRCT for segmentation and quantification of small airway air trapping and emphysema in obstructive pulmonary disease

RATIONALE AND OBJECTIVES

Obstructive pulmonary disease phenotypes are related to variable combinations of emphysema and small-airway disease, the latter manifested as air trapping (AT) on imaging. The investigators propose a method to extract AT information quantitatively from thoracic multi-detector row high-resolution computed tomography (HRCT), validated by pulmonary function testing (PFT) correlation.

MATERIALS AND METHODS

Seventeen patients with obstructive pulmonary disease who underwent HRCT and PFT within a 3-day interval were retrospectively identified. Thin-section volumetric HRCT in inspiration and expiration was registered and analyzed using custom-made software. Nonaerated regions of lung were segmented through exclusion of voxels > -50 Hounsfield units (HU); emphysematous areas were segmented as voxels < -950 HU on inspiratory images. Small-airway AT volume (ATV) was segmented as regions of lung voxels whose attenuation values increased by less than a specified change threshold (set from 5 to 300 HU in 25-HU increments) between inspiration and expiration. Inspiratory and expiratory total segmented lung volumes, emphysema volume (EV), and ATV for each threshold were subsequently calculated and correlated with PFT parameters.

RESULTS

A strong positive correlation was obtained between total segmented lung volume in inspiration and total lung capacity (r = 0.83). A strong negative correlation (r = -0.80) was obtained between EV and the ratio between forced expiratory volume in 1 second and forced vital capacity. Stronger negative correlation with forced expiratory volume in 1 second/forced vital capacity (r = -0.85) was demonstrated when ATV (threshold, 50 HU) was added to EV, indicating improved quantification of total AT to predict obstructive disease severity. A moderately strong positive correlation between ATV and residual volume was observed, with a maximum r value of 0.72 (threshold, 25 HU), greater than that between EV and residual volume (r = 0.58). The benefit of ATV quantification was greater in a subgroup of patients with negligible emphysema compared to patients with moderate to severe emphysema.

CONCLUSIONS

Small-airway AT segmentation in conjunction with emphysema segmentation through computer-assisted methodologies may provide better correlations with key PFT parameters, suggesting that the quantification of emphysema-related and small airway-related components of AT from thoracic HRCT has great potential to elucidate phenotypic differences in patients with chronic obstructive pulmonary disease.

Pulmonary kinematics from image data: a review

The effects of certain lung pathologies include alterations in lung physiology negatively affecting pulmonary compliance. Current approaches to diagnosis and treatment assessment of lung disease commonly rely on pulmonary function testing. Such testing is limited to global measures of lung function, neglecting regional measurements, which are critical for early diagnosis and localization of disease. Increased accessibility to medical image acquisition strategies with high spatiotemporal resolution coupled with the development of sophisticated intensity-based and geometric registration techniques has resulted in the recent exploration of modeling pulmonary motion for calculating local measures of deformation. In this review, the authors provide a broad overview of such research efforts for the estimation of pulmonary deformation. This includes discussion of various techniques, current trends in validation approaches, and the public availability of software and data resources.

Deformation analysis of 3D tagged cardiac images using an optical flow method

BACKGROUND

This study proposes and validates a method of measuring 3D strain in myocardium using a 3D Cardiovascular Magnetic Resonance (CMR) tissue-tagging sequence and a 3D optical flow method (OFM).

METHODS

Initially, a 3D tag MR sequence was developed and the parameters of the sequence and 3D OFM were optimized using phantom images with simulated deformation. This method then was validated in-vivo and utilized to quantify normal sheep left ventricular functions.

RESULTS

Optimizing imaging and OFM parameters in the phantom study produced sub-pixel root-mean square error (RMS) between the estimated and known displacements in the x (RMSx = 0.62 pixels (0.43 mm)), y (RMSy = 0.64 pixels (0.45 mm)) and z (RMSz = 0.68 pixels (1 mm)) direction, respectively. In-vivo validation demonstrated excellent correlation between the displacement measured by manually tracking tag intersections and that generated by 3D OFM (R >or= 0.98). Technique performance was maintained even with 20% Gaussian noise added to the phantom images. Furthermore, 3D tracking of 3D cardiac motions resulted in a 51% decrease in in-plane tracking error as compared to 2D tracking. The in-vivo function studies showed that maximum wall thickening was greatest in the lateral wall, and increased from both apex and base towards the mid-ventricular region. Regional deformation patterns are in agreement with previous studies on LV function.

CONCLUSION

A novel method was developed to measure 3D LV wall deformation rapidly with high in-plane and through-plane resolution from one 3D cine acquisition.

Determination of Respiratory Motion for Distal Esophagus Cancer Using Four-Dimensional Computed Tomography

PURPOSE

To investigate the motion characteristics of distal esophagus cancer primary tumors using four-dimensional computed tomography (4D CT).

METHODS AND MATERIALS

Thirty-one consecutive patients treated for esophagus cancer who received respiratory-gated 4D CT imaging for treatment planning were selected. Deformable image registration was used to map the full expiratory motion gross tumor volume (GTV) to the full-inspiratory CT image, allowing quantitative assessment of each voxel's displacement. These displacements were correlated with patient tumor and respiratory characteristics.

RESULTS

The mean (SE) tidal volume was 608 (73) mL. The mean GTV volume was 64.3 (10.7) mL on expiration and 64.1 (10.7) mL on inspiration (no significant difference). The mean tumor motion in the x-direction was 0.13 (0.006) cm (average of absolute values), in the y-direction 0.23 (0.01) cm (anteriorly), and in the z-direction 0.71 (0.02) cm (inferiorly). Tumor motion correlated with tidal volume. Comparison of tumor motion above vs. below the diaphragm was significant for the average net displacement (p = 0.014), motion below the diaphragm was greater than above. From the cumulative distribution 95% of the tumors moved less than 0.80 cm radially and 1.75 cm inferiorly.

CONCLUSIONS

Primary esophagus tumor motion was evaluated with 4D CT. According to the results of this study, when 4D CT is not available, a radial margin of 0.8 cm and axial margin of +/-1.8 cm would provide tumor motion coverage for 95% of the cases in our study population.