Quantitative CT Characteristics of Cluster Phenotypes in the Severe Asthma Research Program Cohorts.

Background Clustering key clinical characteristics of participants in the Severe Asthma Research Program (SARP), a large, multicenter prospective observational study of patients with asthma and healthy controls, has led to the identification of novel asthma phenotypes. Purpose To determine whether quantitative CT (qCT) could help distinguish between clinical asthma phenotypes. Materials and Methods A retrospective cross-sectional analysis was conducted with the use of qCT images (maximal bronchodilation at total lung capacity [TLC], or inspiration, and functional residual capacity [FRC], or expiration) from the cluster phenotypes of SARP participants (cluster 1: minimal disease; cluster 2: mild, reversible; cluster 3: obese asthma; cluster 4: severe, reversible; cluster 5: severe, irreversible) enrolled between September 2001 and December 2015. Airway morphometry was performed along standard paths (RB1, RB4, RB10, LB1, and LB10). Corresponding voxels from TLC and FRC images were mapped with use of deformable image registration to characterize disease probability maps (DPMs) of functional small airway disease (fSAD), voxel-level volume changes (Jacobian), and isotropy (anisotropic deformation index [ADI]). The association between cluster assignment and qCT measures was evaluated using linear mixed models. Results A total of 455 participants were evaluated with cluster assignments and CT (mean age ± SD, 42.1 years ± 14.7; 270 women). Airway morphometry had limited ability to help discern between clusters. DPM fSAD was highest in cluster 5 (cluster 1 in SARP III: 19.0% ± 20.6; cluster 2: 18.9% ± 13.3; cluster 3: 24.9% ± 13.1; cluster 4: 24.1% ± 8.4; cluster 5: 38.8% ± 14.4; P < .001). Lower whole-lung Jacobian and ADI values were associated with greater cluster severity. Compared to cluster 1, cluster 5 lung expansion was 31% smaller (Jacobian in SARP III cohort: 2.31 ± 0.6 vs 1.61 ± 0.3, respectively, P < .001) and 34% more isotropic (ADI in SARP III cohort: 0.40 ± 0.1 vs 0.61 ± 0.2, P < .001). Within-lung Jacobian and ADI SDs decreased as severity worsened (Jacobian SD in SARP III cohort: 0.90 ± 0.4 for cluster 1; 0.79 ± 0.3 for cluster 2; 0.62 ± 0.2 for cluster 3; 0.63 ± 0.2 for cluster 4; and 0.41 ± 0.2 for cluster 5; P < .001). Conclusion Quantitative CT assessments of the degree and intraindividual regional variability of lung expansion distinguished between well-established clinical phenotypes among participants with asthma from the Severe Asthma Research Program study. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Verschakelen in this issue.

Association of Computed Tomography Densitometry with Disease Severity, Functional Decline, and Survival in Systemic Sclerosis-associated Interstitial Lung Disease.

Rationale: Measuring disease extent and progression of systemic sclerosis-associated interstitial lung disease (SSc-ILD) is challenging, with recent studies suggesting potential utility of quantitative measurements from computed tomography (CT) scans.Objectives: To determine the associations of quantitative computed tomography (qCT) density-based measures with physiological parameters, visual CT scores, and survival in patients with SSc-ILD.Methods: Patients with SSc-ILD and volumetric high-resolution CT images with ≤1.25-mm slice thickness were retrospectively identified. Cardiothoracic radiologists produced visual CT scores of ground glass, reticulation, and honeycombing, with visual fibrosis score equaling the sum of reticulation and honeycombing. qCT measurements included high-attenuation areas (HAA), skewness, kurtosis, and mean lung attenuation (MLA). Associations of qCT measures with pulmonary physiology, visual CT scores, and mortality were analyzed using Spearman's rank correlation and Cox regression.Results: A total of 503 CT scans from 170 patients with SSc-ILD were included. qCT HAA, skewness, kurtosis, and MLA were associated with lung function and visual fibrosis scores, independent of age, sex, and pack-years, using both baseline and change data. Baseline and changes in qCT measures (except ∆skewness) were associated with mortality on unadjusted analysis. Changes in all qCT variables remained associated with survival after adjustment for baseline age, sex, pack-years, and lung function, but not when adjusting for changes in lung function. ∆HAA and ∆MLA were associated with survival after adjustment for age, sex, pack-years, and change in visual CT scores.Conclusions: CT density measurements correlate with physiologic impairment and visual CT scores in patients with SSc-ILD; however, they were not associated with survival independent of changes in pulmonary physiology. The clinical utility of more sophisticated qCT measures should be explored.

Dependence of subject-specific parameters for a fast helical CT respiratory motion model on breathing rate: an animal study.

To determine if the parameters relating lung tissue displacement to a breathing surrogate signal in a previously published respiratory motion model vary with the rate of breathing during image acquisition. An anesthetized pig was imaged using multiple fast helical scans to sample the breathing cycle with simultaneous surrogate monitoring. Three datasets were collected while the animal was mechanically ventilated with different respiratory rates: 12 bpm (breaths per minute), 17 bpm, and 24 bpm. Three sets of motion model parameters describing the correspondences between surrogate signals and tissue displacements were determined. The model error was calculated individually for each dataset, as well asfor pairs of parameters and surrogate signals from different experiments. The values of one model parameter, a vector field denoted [Formula: see text] which related tissue displacement to surrogate amplitude, determined for each experiment were compared. The mean model error of the three datasets was 1.00 ± 0.36 mm with a 95th percentile value of 1.69 mm. The mean error computed from all combinations of parameters and surrogate signals from different datasets was 1.14 ± 0.42 mm with a 95th percentile of 1.95 mm. The mean difference in [Formula: see text] over all pairs of experiments was 4.7% ± 5.4%, and the 95th percentile was 16.8%. The mean angle between pairs of [Formula: see text] was 5.0 ± 4.0 degrees, with a 95th percentile of 13.2 mm. The motion model parameters were largely unaffected by changes in the breathing rate during image acquisition. The mean error associated with mismatched sets of parameters and surrogate signals was 0.14 mm greater than the error achieved when using parameters and surrogate signals acquired with the same breathing rate, while maximum respiratory motion was 23.23 mm on average.

Comparison of low- and ultralow-dose computed tomography protocols for quantitative lung and airway assessment.

PURPOSE: Quantitative computed tomography (CT) measures are increasingly being developed and used to characterize lung disease. With recent advances in CT technologies, we sought to evaluate the quantitative accuracy of lung imaging at low- and ultralow-radiation doses with the use of iterative reconstruction (IR), tube current modulation (TCM), and spectral shaping. METHODS: We investigated the effect of five independent CT protocols reconstructed with IR on quantitative airway measures and global lung measures using an in vivo large animal model as a human subject surrogate. A control protocol was chosen (NIH-SPIROMICS + TCM) and five independent protocols investigating TCM, low- and ultralow-radiation dose, and spectral shaping. For all scans, quantitative global parenchymal measurements (mean, median and standard deviation of the parenchymal HU, along with measures of emphysema) and global airway measurements (number of segmented airways and pi10) were generated. In addition, selected individual airway measurements (minor and major inner diameter, wall thickness, inner and outer area, inner and outer perimeter, wall area fraction, and inner equivalent circle diameter) were evaluated. Comparisons were made between control and target protocols using difference and repeatability measures. RESULTS: Estimated CT volume dose index (CTDIvol) across all protocols ranged from 7.32 mGy to 0.32 mGy. Low- and ultralow-dose protocols required more manual editing and resolved fewer airway branches; yet, comparable pi10 whole lung measures were observed across all protocols. Similar trends in acquired parenchymal and airway measurements were observed across all protocols, with increased measurement differences using the ultralow-dose protocols. However, for small airways (1.9 ± 0.2 mm) and medium airways (5.7 ± 0.4 mm), the measurement differences across all protocols were comparable to the control protocol repeatability across breath holds. Diameters, wall thickness, wall area fraction, and equivalent diameter had smaller measurement differences than area and perimeter measurements. CONCLUSIONS: In conclusion, the use of IR with low- and ultralow-dose CT protocols with CT volume dose indices down to 0.32 mGy maintains selected quantitative parenchymal and airway measurements relevant to pulmonary disease characterization.

Standardizing CT lung density measure across scanner manufacturers.

PURPOSE: Computed Tomography (CT) imaging of the lung, reported in Hounsfield Units (HU), can be parameterized as a quantitative image biomarker for the diagnosis and monitoring of lung density changes due to emphysema, a type of chronic obstructive pulmonary disease (COPD). CT lung density metrics are global measurements based on lung CT number histograms, and are typically a quantity specifying either the percentage of voxels with CT numbers below a threshold, or a single CT number below which a fixed relative lung volume, nth percentile, falls. To reduce variability in the density metrics specified by CT attenuation, the Quantitative Imaging Biomarkers Alliance (QIBA) Lung Density Committee has organized efforts to conduct phantom studies in a variety of scanner models to establish a baseline for assessing the variations in patient studies that can be attributed to scanner calibration and measurement uncertainty. METHODS: Data were obtained from a phantom study on CT scanners from four manufacturers with several protocols at various tube potential voltage (kVp) and exposure settings. Free from biological variation, these phantom studies provide an assessment of the accuracy and precision of the density metrics across platforms solely due to machine calibration and uncertainty of the reference materials. The phantom used in this study has three foam density references in the lung density region, which, after calibration against a suite of Standard Reference Materials (SRM) foams with certified physical density, establishes a HU-electron density relationship for each machine-protocol. We devised a 5-step calibration procedure combined with a simplified physical model that enabled the standardization of the CT numbers reported across a total of 22 scanner-protocol settings to a single energy (chosen at 80 keV). A standard deviation was calculated for overall CT numbers for each density, as well as by scanner and other variables, as a measure of the variability, before and after the standardization. In addition, a linear mixed-effects model was used to assess the heterogeneity across scanners, and the 95% confidence interval of the mean CT number was evaluated before and after the standardization. RESULTS: We show that after applying the standardization procedures to the phantom data, the instrumental reproducibility of the CT density measurement of the reference foams improved by more than 65%, as measured by the standard deviation of the overall mean CT number. Using the lung foam that did not participate in the calibration as a test case, a mixed effects model analysis shows that the 95% confidence intervals are [-862.0 HU, -851.3 HU] before standardization, and [-859.0 HU, -853.7 HU] after standardization to 80 keV. This is in general agreement with the expected CT number value at 80 keV of -855.9 HU with 95% CI of [-857.4 HU, -854.5 HU] based on the calibration and the uncertainty in the SRM certified density. CONCLUSIONS: This study provides a quantitative assessment of the variations expected in CT lung density measures attributed to non-biological sources such as scanner calibration and scanner x-ray spectrum and filtration. By removing scanner-protocol dependence from the measured CT numbers, higher accuracy and reproducibility of quantitative CT measures were attainable. The standardization procedures developed in study may be explored for possible application in CT lung density clinical data.

SPIROMICS Protocol for Multicenter Quantitative Computed Tomography to Phenotype the Lungs.

Multidetector row computed tomography (MDCT) is increasingly taking a central role in identifying subphenotypes within chronic obstructive pulmonary disease (COPD), asthma, and other lung-related disease populations, allowing for the quantification of the amount and distribution of altered parenchyma along with the characterization of airway and vascular anatomy. The embedding of quantitative CT (QCT) into a multicenter trial with a variety of scanner makes and models along with the variety of pressures within a clinical radiology setting has proven challenging, especially in the context of a longitudinal study. SPIROMICS (Subpopulations and Intermediate Outcome Measures in COPD Study), sponsored by the National Institutes of Health, has established a QCT lung assessment system (QCT-LAS), which includes scanner-specific imaging protocols for lung assessment at total lung capacity and residual volume. Also included are monthly scanning of a standardized test object and web-based tools for subject registration, protocol assignment, and data transmission coupled with automated image interrogation to assure protocol adherence. The SPIROMICS QCT-LAS has been adopted and contributed to by a growing number of other multicenter studies in which imaging is embedded. The key components of the SPIROMICS QCT-LAS along with evidence of implementation success are described herein. While imaging technologies continue to evolve, the required components of a QCT-LAS provide the framework for future studies, and the QCT results emanating from SPIROMICS and the growing number of other studies using the SPIROMICS QCT-LAS will provide a shared resource of image-derived pulmonary metrics.

A controlled statistical study to assess measurement variability as a function of test object position and configuration for automated surveillance in a multicenter longitudinal COPD study (SPIROMICS).

PURPOSE: A test object (phantom) is an important tool to evaluate comparability and stability of CT scanners used in multicenter and longitudinal studies. However, there are many sources of error that can interfere with the test object-derived quantitative measurements. Here the authors investigated three major possible sources of operator error in the use of a test object employed to assess pulmonary density-related as well as airway-related metrics. METHODS: Two kinds of experiments were carried out to assess measurement variability caused by imperfect scanning status. The first one consisted of three experiments. A COPDGene test object was scanned using a dual source multidetector computed tomographic scanner (Siemens Somatom Flash) with the Subpopulations and Intermediate Outcome Measures in COPD Study (SPIROMICS) inspiration protocol (120 kV, 110 mAs, pitch = 1, slice thickness = 0.75 mm, slice spacing = 0.5 mm) to evaluate the effects of tilt angle, water bottle offset, and air bubble size. After analysis of these results, a guideline was reached in order to achieve more reliable results for this test object. Next the authors applied the above findings to 2272 test object scans collected over 4 years as part of the SPIROMICS study. The authors compared changes of the data consistency before and after excluding the scans that failed to pass the guideline. RESULTS: This study established the following limits for the test object: tilt index ≤0.3, water bottle offset limits of [-6.6 mm, 7.4 mm], and no air bubble within the water bottle, where tilt index is a measure incorporating two tilt angles around x- and y-axis. With 95% confidence, the density measurement variation for all five interested materials in the test object (acrylic, water, lung, inside air, and outside air) resulting from all three error sources can be limited to ±0.9 HU (summed in quadrature), when all the requirements are satisfied. The authors applied these criteria to 2272 SPIROMICS scans and demonstrated a significant reduction in measurement variation associated with the test object. CONCLUSIONS: Three operator errors were identified which significantly affected the usability of the acquired scan images of the test object used for monitoring scanner stability in a multicenter study. The authors' results demonstrated that at the time of test object scan receipt at a radiology core laboratory, quality control procedures should include an assessment of tilt index, water bottle offset, and air bubble size within the water bottle. Application of this methodology to 2272 SPIROMICS scans indicated that their findings were not limited to the scanner make and model used for the initial test but was generalizable to both Siemens and GE scanners which comprise the scanner types used within the SPIROMICS study.

Acute administration of ivacaftor to people with cystic fibrosis and a G551D-CFTR mutation reveals smooth muscle abnormalities.

BACKGROUND: Airflow obstruction is common in cystic fibrosis (CF), yet the underlying pathogenesis remains incompletely understood. People with CF often exhibit airway hyperresponsiveness, CF transmembrane conductance regulator (CFTR) is present in airway smooth muscle (ASM), and ASM from newborn CF pigs has increased contractile tone, suggesting that loss of CFTR causes a primary defect in ASM function. We hypothesized that restoring CFTR activity would decrease smooth muscle tone in people with CF. METHODS: To increase or potentiate CFTR function, we administered ivacaftor to 12 adults with CF with the G551D-CFTR mutation; ivacaftor stimulates G551D-CFTR function. We studied people before and immediately after initiation of ivacaftor (48 hours) to minimize secondary consequences of CFTR restoration. We tested smooth muscle function by investigating spirometry, airway distensibility, and vascular tone. RESULTS: Ivacaftor rapidly restored CFTR function, indicated by reduced sweat chloride concentration. Airflow obstruction and air trapping also improved. Airway distensibility increased in airways less than 4.5 mm but not in larger-sized airways. To assess smooth muscle function in a tissue outside the lung, we measured vascular pulse wave velocity (PWV) and augmentation index, which both decreased following CFTR potentiation. Finally, change in distensibility of <4.5-mm airways correlated with changes in PWV. CONCLUSIONS: Acute CFTR potentiation provided a unique opportunity to investigate CFTR-dependent mechanisms of CF pathogenesis. The rapid effects of ivacaftor on airway distensibility and vascular tone suggest that CFTR dysfunction may directly cause increased smooth muscle tone in people with CF and that ivacaftor may relax smooth muscle. FUNDING: This work was funded in part from an unrestricted grant from the Vertex Investigator-Initiated Studies Program.

Computed Tomography and Magnetic Resonance Imaging for Longitudinal Characterization of Lung Structure Changes in a Yucatan Miniature Pig Silicosis Model.

Medical imaging is a rapidly advancing field enabling the repeated, noninvasive assessment of physiological structure and function. These beneficial characteristics can supplement studies in swine by mirroring the clinical functions of detection, diagnosis, and monitoring in humans. In addition, swine may serve as a human surrogate, facilitating the development and comparison of new imaging protocols for translation to humans. This study presents methods for pulmonary imaging developed for monitoring pulmonary disease initiation and progression in a pig exposure model with computed tomography and magnetic resonance imaging. In particular, a focus was placed on systematic processes, including positioning, image acquisition, and structured reporting to monitor longitudinal change. The image-based monitoring procedure was applied to 6 Yucatan miniature pigs. A subset of animals (n= 3) were injected with crystalline silica into the apical bronchial tree to induce silicosis. The methodology provided longitudinal monitoring and evidence of progressive lung disease while simultaneously allowing for a cross-modality comparative study highlighting the practical application of medical image data collection in swine. The integration of multimodality imaging with structured reporting allows for cross comparison of modalities, refinement of CT and MRI protocols, and consistently monitors potential areas of interest for guided biopsy and/or necropsy.

Comparison of breathing gated CT images generated using a 5DCT technique and a commercial clinical protocol in a porcine model.

PURPOSE: To demonstrate that a "5DCT" technique which utilizes fast helical acquisition yields the same respiratory-gated images as a commercial technique for regular, mechanically produced breathing cycles. METHODS: Respiratory-gated images of an anesthetized, mechanically ventilated pig were generated using a Siemens low-pitch helical protocol and 5DCT for a range of breathing rates and amplitudes and with standard and low dose imaging protocols. 5DCT reconstructions were independently evaluated by measuring the distances between tissue positions predicted by a 5D motion model and those measured using deformable registration, as well by reconstructing the originally acquired scans. Discrepancies between the 5DCT and commercial reconstructions were measured using landmark correspondences. RESULTS: The mean distance between model predicted tissue positions and deformably registered tissue positions over the nine datasets was 0.65 ± 0.28 mm. Reconstructions of the original scans were on average accurate to 0.78 ± 0.57 mm. Mean landmark displacement between the commercial and 5DCT images was 1.76 ± 1.25 mm while the maximum lung tissue motion over the breathing cycle had a mean value of 27.2 ± 4.6 mm. An image composed of the average of 30 deformably registered images acquired with a low dose protocol had 6 HU image noise (single standard deviation) in the heart versus 31 HU for the commercial images. CONCLUSIONS: An end to end evaluation of the 5DCT technique was conducted through landmark based comparison to breathing gated images acquired with a commercial protocol under highly regular ventilation. The techniques were found to agree to within 2 mm for most respiratory phases and most points in the lung.

Very low-dose (0.15 mGy) chest CT protocols using the COPDGene 2 test object and a third-generation dual-source CT scanner with corresponding third-generation iterative reconstruction software.

OBJECTIVES: The purpose of this study was to evaluate the impact of ultralow radiation dose single-energy computed tomographic (CT) acquisitions with Sn prefiltration and third-generation iterative reconstruction on density-based quantitative measures of growing interest in phenotyping pulmonary disease. MATERIALS AND METHODS: The effects of both decreasing dose and different body habitus on the accuracy of the mean CT attenuation measurements and the level of image noise (SD) were evaluated using the COPDGene 2 test object, containing 8 different materials of interest ranging from air to acrylic and including various density foams. A third-generation dual-source multidetector CT scanner (Siemens SOMATOM FORCE; Siemens Healthcare AG, Erlangen, Germany) running advanced modeled iterative reconstruction (ADMIRE) software (Siemens Healthcare AG) was used.We used normal and very large body habitus rings at dose levels varying from 1.5 to 0.15 mGy using a spectral-shaped (0.6-mm Sn) tube output of 100 kV(p). Three CT scans were obtained at each dose level using both rings. Regions of interest for each material in the test object scans were automatically extracted. The Hounsfield unit values of each material using weighted filtered back projection (WFBP) at 1.5 mGy was used as the reference value to evaluate shifts in CT attenuation at lower dose levels using either WFBP or ADMIRE. Statistical analysis included basic statistics, Welch t tests, multivariable covariant model using the F test to assess the significance of the explanatory (independent) variables on the response (dependent) variable, and CT mean attenuation, in the multivariable covariant model including reconstruction method. RESULTS: Multivariable regression analysis of the mean CT attenuation values showed a significant difference with decreasing dose between ADMIRE and WFBP. The ADMIRE has reduced noise and more stable CT attenuation compared with WFBP. There was a strong effect on the mean CT attenuation values of the scanned materials for ring size (P < 0.0001) and dose level (P < 0.0001). The number of voxels in the region of interest for the particular material studied did not demonstrate a significant effect (P > 0.05). The SD was lower with ADMIRE compared with WFBP at all dose levels and ring sizes (P < 0.05). CONCLUSIONS: The third-generation dual-source CT scanners using third-generation iterative reconstruction methods can acquire accurate quantitative CT images with acceptable image noise at very low-dose levels (0.15 mGy). This opens up new diagnostic and research opportunities in CT phenotyping of the lung for developing new treatments and increased understanding of pulmonary disease.

Scatter correction associated with dedicated dual-source CT hardware improves accuracy of lung air measures.

RATIONALE AND OBJECTIVES: Accurate assessment of air density used to quantitatively characterize amount and distribution of emphysema in chronic obstructive pulmonary disease (COPD) subjects has remained challenging. Hounsfield units (HU) within tracheal air can be considerably less negative than -1000 HU. This study has sought to characterize the effects of improved scatter correction used in dual-source pulmonary computed tomography (CT). MATERIALS AND METHODS: Dual-source dual-energy (DSDE) and single-source (SS) scans taken at multiple energy levels and scan settings were acquired for quantitative comparison using anesthetized ovine (n = 6), swine (n = 13), and a lung phantom. Data were evaluated for the lung, inferior vena cava, and tracheal segments. To minimize the effect of cross-scatter, the phantom scans in the DSDE mode were obtained by reducing the current of one of the tubes to near zero. RESULTS: A significant shift in mean HU values in the tracheal regions of animals and the phantom is observed, with values consistently closer to -1000 HU in DSDE mode. HU values associated with SS mode demonstrated a positive shift of up to 32 HU. In vivo tracheal air measurements demonstrated considerable variability with SS scanning, whereas these values were more consistent with DSDE imaging. Scatter effects in the lung parenchyma differed from adjacent tracheal measures. CONCLUSION: Data suggest that the scatter correction introduced into the dual-energy mode of imaging has served to provide more accurate CT lung density measures sought to quantitatively assess the presence and distribution of emphysema in COPD subjects. Data further suggest that CT images, acquired without adequate scatter correction, cannot be corrected by linear algorithms given the variability in tracheal air HU values and the independent scatter effects on lung parenchyma.

Effect of lung inflation level on hyperpolarized 3He apparent diffusion coefficient measurements in never-smokers.

PURPOSE: To evaluate the effects of lung volume differences on apparent diffusion coefficient (ADC) measurements on a regional basis, with breath holds at volumes adjusted for differences in lung size across individuals according to the subject's vital capacity (VC). MATERIALS AND METHODS: This study was approved by the local institutional review board and was compliant with HIPAA. Informed consent was obtained from all subjects. Imaging was performed under a physician's Investigational New Drug application from the Food and Drug Administration. ADC changes as a function of inflation levels were evaluated in 24 healthy never-smokers across three lung volumes (20%, 60%, and 100% VC) on the basis of the spirometric data collected from each subject. Response variables based on lung volume and anatomic position were assessed with multifactorial analysis of variance followed by posthoc pair-wise testing. Imaging was performed with a 1.5-T magnetic resonance (MR) unit with use of a two-dimensional gradient-echo fast low-angle shot sequence. RESULTS: Significant differences in ADCs between lung volumes were observed for all inflation levels (20%, 60%, and 100% VC; P < .001), along with significant dependent-nondependent vertical gradients at 20% VC (P < .0001) and 60% VC (P < .0001, left lung only). In addition, significant differences between mean values in the left and right lungs with respect to those in the whole lung were observed at the lower lung inflation levels (20% and 60% VC, P < .01), reaching more uniform expansion at 100% VC. CONCLUSION: The results confirm known anatomic differences in patterns of regional inflation and ventilation with corresponding lung volume changes, emphasizing the need for tight control over lung volume when performing hyperpolarized helium 3 ((3)He) lung studies if (3)He MR imaging is to be used to follow up small longitudinal changes in lung abnormalities.

A process model for direct correlation between computed tomography and histopathology application in lung cancer.

RATIONALE AND OBJECTIVES: Multimodal imaging techniques for capturing normal and diseased human anatomy and physiology are being developed to benefit patient clinical care, research, and education. In the past, the incorporation of histopathology into these multimodal datasets has been complicated by the large differences in image quality, content, and spatial association. MATERIALS AND METHODS: We have developed a novel system, the large-scale image microtome array (LIMA), to bridge the gap between nonstructurally destructive and destructive imaging such that reliable registration between radiological data and histopathology can be achieved. Registration algorithms have been designed to align the multimodal datasets, which include computed tomography, computed micro-tomography, LIMA, and histopathology data to a common coordinate system. RESULTS: The resulting volumetric dataset provides an abundance of valuable information relating to the tissue sample including density, anatomical structure, color, texture, and cellular information in three dimensions. An image processing pipeline has been established to register all the multimodal data to a common coordinate system. CONCLUSION: In this study, we have chosen to use human lung cancer nodules as an example; however, the flexibility of the image acquisition and subsequent processing algorithms makes it applicable to any soft organ tissue. A novel process model has been established to generate cross registered multimodal datasets for the investigation of human lung cancer nodule content and associated image-based representation.