Dakai Jin

Quantitative Dual-Energy Computed Tomography Supports a Vascular Etiology of Smoking-induced Inflammatory Lung Disease

RATIONALE Endothelial dysfunction is of interest in relation to smoking-associated emphysema, a component of chronic obstructive pulmonary disease (COPD). We previously demonstrated that computed tomography (CT)-derived pulmonary blood flow (PBF) heterogeneity is greater in smokers with normal pulmonary function tests (PFTs) but who have visual evidence of centriacinar emphysema (CAE) on CT. OBJECTIVES We introduced dual-energy CT (DECT) perfused blood volume (PBV) as a PBF surrogate to evaluate whether the CAE-associated increased PBF heterogeneity is reversible with sildenafil. METHODS Seventeen PFT-normal current smokers were divided into CAE-susceptible (SS; n = 10) and nonsusceptible (NS; n = 7) smokers, based on the presence or absence of CT-detected CAE. DECT-PBV images were acquired before and 1 hour after administration of 20 mg oral sildenafil. Regional PBV and PBV coefficients of variation (CV), a measure of spatial blood flow heterogeneity, were determined, followed by quantitative assessment of the central arterial tree. MEASUREMENTS AND MAIN RESULTS After sildenafil administration, regional PBV-CV decreased in SS subjects but did not decrease in NS subjects (P < 0.05), after adjusting for age and pack-years. Quantitative evaluation of the central pulmonary arteries revealed higher arterial volume and greater cross-sectional area (CSA) in the lower lobes of SS smokers, which suggested arterial enlargement in response to increased peripheral resistance. After sildenafil, arterial CSA decreased in SS smokers but did not decrease in NS smokers (P < 0.01). CONCLUSIONS These results demonstrate that sildenafil restores peripheral perfusion and reduces central arterial enlargement in normal SS subjects with little effect in NS subjects, highlighting DECT-PBV as a biomarker of reversible endothelial dysfunction in smokers with CAE.

A Robust Algorithm for Thickness Computation at Low Resolution and Its Application to In Vivo Trabecular Bone CT Imaging

Adult bone diseases, especially osteoporosis, lead to increased risk of fracture which in turn is associated with substantial morbidity, mortality, and financial costs. Clinically, osteoporosis is defined by low bone mineral density; however, increasing evidence suggests that the microarchitectural quality of trabecular bone (TB) is an important determinant of bone strength and fracture risk. Accurate measures of TB thickness and marrow spacing is of significant interest for early diagnosis of osteoporosis or treatment effects. Here, we present a new robust algorithm for computing TB thickness and marrow spacing at a low resolution achievable in vivo. The method uses a star-line tracing technique that effectively deals with partial voluming effects of in vivo imaging with voxel size comparable to TB thickness. Also, the method avoids the problem of digitization associated with conventional algorithms based on sampling distance transform along skeletons. Accuracy of the method was examined using computer-generated phantom images, while the robustness of the method was evaluated on human ankle specimens in terms of stability across a wide range of voxel sizes, repeat scan reproducibility under in vivo conditions, and correlation between thickness values computed at ex vivo and in vivo imaging resolutions. Also, the sensitivity of the method was examined by evaluating its ability to predict the bone strength of cadaveric specimens. Finally, the method was evaluated in a human study involving 40 healthy young-adult volunteers (age: 19-21 years; 20 males and 20 females) and ten athletes (age: 19- 21 years; six males and four females). Across a wide range of voxel sizes, the new method is significantly more accurate and robust as compared to conventional methods. Both TB thickness and marrow spacing measures computed using the new method demonstrated strong associations (R2 ∈ [0.83, 0.87]) with bone strength. Also, the TB thickness and marrow spacing measures allowed discrimination between male and female volunteers (p ∈ [0.01, 0.04]) as well as between athletes and nonathletes (p ∈ [0.005, 0.03]).

A New Fuzzy Skeletonization Algorithm and Its Applications to Medical Imaging

Skeletonization provides a simple yet compact representation of an object and is widely used in medical imaging applications including volumetric, structural, and topological analyses, object representation, stenoses detection, path-finding etc. Literature of three-dimensional skeletonization is quite matured for binary digital objects. However, the challenges of skeletonization for fuzzy objects are mostly unanswered. Here, a framework and an algorithm for fuzzy surface skeletonization are developed using a notion of fuzzy grassfire propagation which will minimize binarization related data loss. Several concepts including fuzzy axial voxels, local and global significance factors are introduced. A skeletal noise pruning algorithm using global significance factors as significance measures of individual branches is developed. Results of application of the algorithm on several medical objects have been illustrated. A quantitative comparison with an ideal skeleton has demonstrated that the algorithm can achieve sub-voxel accuracies at various levels of noise and downsampling. The role of fuzzy skeletonization in thickness computation at relatively low resolution has been demonstrated.

A robust and efficient curve skeletonization algorithm for tree-like objects using minimum cost paths

Conventional curve skeletonization algorithms using the principle of Blums transform, often, produce unwanted spurious branches due to boundary irregularities, digital effects, and other artifacts. This paper presents a new robust and efficient curve skeletonization algorithm for three-dimensional (3-D) elongated fuzzy objects using a minimum cost path approach, which avoids spurious branches without requiring post-pruning. Starting from a root voxel, the method iteratively expands the skeleton by adding new branches in each iteration that connects the farthest quench voxel to the current skeleton using a minimum cost path. The path-cost function is formulated using a novel measure of local significance factor defined by the fuzzy distance transform field, which forces the path to stick to the centerline of an object. The algorithm terminates when dilated skeletal branches fill the entire object volume or the current farthest quench voxel fails to generate a meaningful skeletal branch. Accuracy of the algorithm has been evaluated using computer-generated phantoms with known skeletons. Performance of the method in terms of false and missing skeletal branches, as defined by human experts, has been examined using in vivo CT imaging of human intrathoracic airways. Results from both experiments have established the superiority of the new method as compared to the existing methods in terms of accuracy as well as robustness in detecting true and false skeletal branches. The new algorithm makes a significant reduction in computation complexity by enabling detection of multiple new skeletal branches in one iteration. Specifically, this algorithm reduces the number of iterations from the number of terminal tree branches to the worst case performance of tree depth. In fact, experimental results suggest that, on an average, the order of computation complexity is reduced to the logarithm of the number of terminal branches of a tree-like object.

Automated Assessment of Pulmonary Arterial Morphology in Multi-row Detector CT Imaging Using Correspondence with Anatomic Airway Branches

Multi-row detector CT (MDCT) provides high resolution structural and functional imaging that has been helpful in studying altered physiology, making early diagnosis, and evaluating treatments in pulmonary research. There is growing evidence suggesting that pulmonary vascular dysfunction plays a major role in progression of centrilobular emphysema, a component of chronic obstructive disease (COPD). Few studies have attempted to quantify central pulmonary vessel morphology and to compare these measurements across COPD groups. However, the scope of vascular structures examined in such studies has been limited, primarily, due to lack of an automated and standardized method of comparing matching vessel branches. In this paper, we present a fully automated method, using a novel arc skeletonization and a local correspondence analysis, to identify matching pulmonary arteries by linking those with anatomically defined specific airway branches. This method provides a standardized way of establishing correspondence between matched pulmonary arteries for intra- and inter-subject scans. The accuracy and repeatability of the method was examined on non-contrast MDCT scans of 10 normal subjects. It was observed that 83% of the arteries classified by our automated method agree with “true” arteries as labelled by an interactive manual artery-vein separation tool. Repeat scan intra-class correlation of arterial morphological measures over six anatomic airway branches was observed as 91%.

A semi-automatic framework of measuring pulmonary arterial metrics at anatomic airway locations using CT imaging

Pulmonary vascular dysfunction has been implicated in smoking-related susceptibility to emphysema. With the growing interest in characterizing arterial morphology for early evaluation of the vascular role in pulmonary diseases, there is an increasing need for the standardization of a framework for arterial morphological assessment at airway segmental levels. In this paper, we present an effective and robust semi-automatic framework to segment pulmonary arteries at different anatomic airway branches and measure their cross-sectional area (CSA). The method starts with user-specified endpoints of a target arterial segment through a custom-built graphical user interface. It then automatically detect the centerline joining the endpoints, determines the local structure orientation and computes the CSA along the centerline after filtering out the adjacent pulmonary structures, such as veins or airway walls. Several new techniques are presented, including collision-impact based cost function for centerline detection, radial sample-line based CSA computation, and outlier analysis of radial distance to subtract adjacent neighboring structures in the CSA measurement. The method was applied to repeat-scan pulmonary multirow detector CT (MDCT) images from ten healthy subjects (age: 21-48 Yrs, mean: 28.5 Yrs; 7 female) at functional residual capacity (FRC). The reproducibility of computed arterial CSA from four airway segmental regions in middle and lower lobes was analyzed. The overall repeat-scan intra-class correlation (ICC) of the computed CSA from all four airway regions in ten subjects was 96% with maximum ICC found at LB10 and RB4 regions.

Automated cortical bone segmentation for multirow-detector CT imaging with validation and application to human studies.

PURPOSE Cortical bone supports and protects human skeletal functions and plays an important role in determining bone strength and fracture risk. Cortical bone segmentation at a peripheral site using multirow-detector CT (MD-CT) imaging is useful for in vivo assessment of bone strength and fracture risk. Major challenges for the task emerge from limited spatial resolution, low signal-to-noise ratio, presence of cortical pores, and structural complexity over the transition between trabecular and cortical bones. An automated algorithm for cortical bone segmentation at the distal tibia from in vivo MD-CT imaging is presented and its performance and application are examined. METHODS The algorithm is completed in two major steps-(1) bone filling, alignment, and region-of-interest computation and (2) segmentation of cortical bone. After the first step, the following sequence of tasks is performed to accomplish cortical bone segmentation-(1) detection of marrow space and possible pores, (2) computation of cortical bone thickness, detection of recession points, and confirmation and filling of true pores, and (3) detection of endosteal boundary and delineation of cortical bone. Effective generalizations of several digital topologic and geometric techniques are introduced and a fully automated algorithm is presented for cortical bone segmentation. RESULTS An accuracy of 95.1% in terms of volume of agreement with manual outlining of cortical bone was observed in human MD-CT scans, while an accuracy of 88.5% was achieved when compared with manual outlining on postregistered high resolution micro-CT imaging. An intraclass correlation coefficient of 0.98 was obtained in cadaveric repeat scans. A pilot study was conducted to describe gender differences in cortical bone properties. This study involved 51 female and 46 male participants (age: 19-20 yr) from the Iowa Bone Development Study. Results from this pilot study suggest that, on average after adjustment for height and weight differences, males have thicker cortex (mean difference 0.33 mm and effect size 0.92 at the anterior region) with lower bone mineral density (mean difference -28.73 mg/cm(3) and effect size 1.35 at the posterior region) as compared to females. CONCLUSIONS The algorithm presented is suitable for fully automated segmentation of cortical bone in MD-CT imaging of the distal tibia with high accuracy and reproducibility. Analysis of data from a pilot study demonstrated that the cortical bone indices allow quantification of gender differences in cortical bone from MD-CT imaging. Application to larger population groups, including those with compromised bone, is needed.

A New Approach of Arc Skeletonization for Tree-like Objects Using Minimum Cost Path

Traditional arc skeletonization algorithms using the principle of Blums transform, often, produce unwanted spurious branches due to boundary irregularities and digital effects on objects and other artifacts. This paper presents a new robust approach of extracting arc skeletons for three-dimensional (3-D) elongated fuzzy objects, which avoids spurious branches without requiring post-pruning. Starting from a root voxel, the method iteratively expands the skeleton by adding a new branch in each iteration that connects the farthest voxel to the current skeleton using a minimum-cost geodesic path. The path-cost function is formulated using a novel measure of local significance factor defined by fuzzy distance transform field, which forces the path to stick to the centerline of the object. The algorithm terminates when dilated skeletal branches fill the entire object volume or the current farthest voxel fails to generate a meaningful branch. Accuracy of the algorithm has been evaluated using computer-generated blurred and noisy phantoms with known skeletons. Performance of the method in terms of false and missing skeletal branches, as defined by human expert, has been examined using in vivo CT imaging of human intrathoracic airways. Experimental results from both experiments have established the superiority of the new method as compared to a widely used conventional method in terms of accuracy of medialness as well as robustness of true and false skeletal branches.

A new algorithm for trabecular bone thickness computation at low resolution achieved under in vivo condition

Adult bone diseases, especially osteoporosis, lead to increased risk of fracture associated with substantial morbidity, mortality, and financial costs. Clinically, osteoporosis is defined by low bone mineral density (BMD); however, increasing evidence suggests that the micro-architectural quality of trabecular bone (TB) is an important determinant of bone strength and fracture risk. Accurate measurement of trabecular thickness and marrow spacing is of significant interest for early diagnosis of osteoporosis or treatment effects. Here, we present a new robust algorithm for computing TB thickness and marrow spacing at a low resolution achievable in vivo. The method uses a star-line tracing technique that effectively deals with partial voluming effects of in vivo imaging where voxel size is comparable to TB thickness. Experimental results on cadaveric ankle specimens have demonstrated the algorithms robustness (ICC > 0.98) under repeat scans of multi-row detector computed tomography (MD-CT) imaging. It has been observed in experimental results that TB thickness and marrow spacing measures as computed by the new algorithm have strong association (R2 ∈ {0.85, 0.87} ) with TBs experimental mechanical strength measures.

Filtering Non-Significant Quench Points Using Collision Impact in Grassfire Propagation

The skeleton of an object is defined as the set of quench points formed during Blum’s grassfire transformation. Due to high sensitivity of quench points with small changes in the object boundary and the membership function (for fuzzy objects), often, a large number of redundant quench points is formed. Many of these quench points are caused by peripheral protrusions and dents and do not associate themselves with core shape features of the object. Here, we present a significance measure of quench points using the collision impact of fire-fronts and explore its role in filtering noisy quench points. The performance of the method is examined on three-dimensional shapes at different levels of noise and fuzziness, and compared with previous methods. The results have demonstrated that collision impact together with appropriate filtering kernels eliminate most of the noisy quench voxels while preserving those associated with core shape features of the object