Pulmonary emphysema is a destructive disease of the lungs that is currently diagnosed via visual assessment of lung Computed Tomography (CT) scans by a radiologist. Visual assessment can have poor inter-rater reliability, is time consuming, and requires access to trained assessors. Quantitative methods that reliably summarize the biologically relevant characteristics of an image are needed to improve the way lung diseases are characterized. The goal of this work was to show how spatial point process models can be used to create a set of radiologically derived quantitative lung biomarkers of emphysema. We formalized a general framework for applying spatial point processes to lung CT scans, and developed a Shot Noise Cox Process to quantify how radiologically based emphysematous tissue clusters into larger structures. Bayesian estimation of model parameters was done using spatial Birth-Death MCMC (BD-MCMC). In simulations, we showed the BD-MCMC estimation algorithm is able to accurately recover model parameters. In an application to real lung CT scans from the COPDGene cohort, we showed variability in the clustering characteristics of emphysematous tissue across disease subtypes that were based on visual assessments of the CT scans.
Keywords: Birth–death MCMC; COPD; Cox cluster process; Emphysema.