Home Novel Logistic Regression Model of Chest CT Attenuation Coefficient Distributions for the Automated Detection of Abnormal (Emphysema or ILD) Versus Normal Lung
Post
Cancel

Novel Logistic Regression Model of Chest CT Attenuation Coefficient Distributions for the Automated Detection of Abnormal (Emphysema or ILD) Versus Normal Lung

Rationale and Objectives

We evaluated the role of automated quantitative computed tomography (CT) scan interpretation algorithm in detecting interstitial lung disease (ILD) and/or emphysema in a sample of elderly subjects with mild lung disease. We hypothesized that the quantification and distributions of CT attenuation values on lung CT, over a subset of Hounsfield units (HUs) range (−1000 HU, 0 HU), can differentiate early or mild disease from normal lung.

Materials and Methods

We compared the results of quantitative spiral rapid end-exhalation (functional residual capacity, FRC) and end-inhalation (total lung capacity, TLC) CT scan analyses of 52 subjects with radiographic evidence of mild fibrotic lung disease to the results of 17 normal subjects. Several CT value distributions were explored, including (1) that from the peripheral lung taken at TLC (with peels at 15 or 65 mm), (2) the ratio of (1) to that from the core of lung, and (3) the ratio of (2) to its FRC counterpart. We developed a fused-lasso logistic regression model that can automatically identify sub-intervals of −1000 HU and 0 HU over which a CT value distribution provides optimal discrimination between abnormal and normal scans.

Results

The fused-lasso logistic regression model based on (2) with 15-mm peel identified the relative frequency of CT values of over −1000 HU and −900 and those over −450 HU and −200 HU as a means of discriminating abnormal versus normal lung, resulting in a zero out-sample false-positive rate, and 15% false-negative rate of that was lowered to 12% by pooling information.

Conclusions

We demonstrated the potential usefulness of this novel quantitative imaging analysis method in discriminating ILD and/or emphysema from normal lungs.

Interstitial lung disease (ILD) is increasing in importance in part because of the aging population but also detection and/or incidence appear to be increasing. Data from the National Center for Health Statistics indicate that the age-adjusted mortality rate from pulmonary fibrosis has increased by 28.4% in men and by 41.3% in women between 1992 and 2014 .

Inter- and intrareader variability in the interpretation of radiographs for pulmonary fibrosis and pneumoconiosis has been long recognized as a potential issue for screening programs, epidemiologic studies, and medico-legal evaluations . Studies have found variable degrees of agreement for both parenchymal and pleural fibroses dependent on the extent of abnormalities and the training and medical specialty of the chest X-ray (CXR) readers . The National Institute for Occupational Safety and Health recommends using multiple, International Labour Organization radiographic pneumoconiosis classification system (ILO system)-trained readers and median profusion scores as the preferred reconciliation protocol to increase accuracy and precision in PA film classification . Clinical studies suggest that computed tomography (CT) imaging technology may become a gold standard for the evaluation of obstructive airways disease , but the literature is relatively limited on CT’s use in quantifying and objectively characterizing patterns of subtle interstitial fibrosis.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Materials and Methods

Data

Get Radiology Tree app to read full this article<

Table 1

Demographic and Clinical Characteristics of Study Populations

ILD Normal_P_ Value_N_ 52 17 Sex 0.0002 Male 48 (92%) 8 (47%) Female 4 (8%) 9 (53%) Age, year 77.35 ± 8.41 40.35 ± 16.66 <0.0001 Race 0.16 African American 0 (0%) 1 (6%) Caucasian 50 (96%) 16 (94%) Hispanic 2 (4%) 0 (0%) BMI 28.81 ± 5.25 25.65 ± 3.69 0.025 Smoking history 0.00006 Never 20 (38%) 17 (100%) Former 25 (48%) 0 (0%) Current 7 (14%) 0 (0%) Pack-years 42.39 ± 34.47 0 ± 0 0.003 Plaques <0.0001 No 32 (62%) 17 (100%) Yes 20 (38%) 0 (0%) Emphysema <0.0001 No 6 (12%) 17 (100%) Yes 46 (88%) 0 (0%) Bronchiectasis <0.0001 No 9 (17%) 17 (100%) Yes 43 (83%) 0 (0%)

BMI, body mass index; ILD, interstitial lung disease.

Table 2

Pulmonary Function in Patients with ILD Compared to Normal Subjects

ILD Normal_P_ Value_N_ 52 17 TLC 96.41 ± 22.20 103.41 ± 8.52 0.21 RV 96.08 ± 31.26 100.18 ± 17.44 0.61 FVC 82.71 ± 24.11 109.76 ± 19.62 <0.0001 FEV 1 78.40 ± 24.40 112.82 ± 20.66 <0.0001 DLCO 62.25 ± 19.91 126.71 ± 15.81 <0.0001 FEV 1 /FVC 68.88 ± 14.14 82.94 ± 6.43 0.0002

DLCO, diffusion capacity for the lungs of carbon monoxide; FEV 1 , forced expiratory volume in 1 s; FVC, functional vital capacity; RV, residual volume; TLC, total lung capacity.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 1, The plots of the eight function summaries of the CT attenuation value distribution of lung parenchyma. The diagrams in the upper panel, from left to right, plot the realizations of Vi,i=1,2,3,4 Vi,i=1,2,3,4 , with those of the normal lungs drawn as gray solid lines and the abnormal counterparts as dashed lines. The diagrams in the lower panel are those of V′i,i=1,2,3,4 V′i,i=1,2,3,4 . These plots show that the V s and V ′s of the normal lungs fluctuated within narrow bands, whereas those of the abnormal lungs had greater fluctuations with higher skewness to the right. The preceding contrast is generally amplified in the plots of the V ′s.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Logistic Regression Model

Get Radiology Tree app to read full this article<

logit(p)=β0+∑Mj=1βjV(νj) logit

(

p

)

=

β

0

+

j

=

1

M

β

j

V

(

ν

j

)

where logit(p)=log{p/(1−p)} logit

(

p

)

=

log

{

p

/

(

1

p

)

} is the logistic transformation; p is the probability that the subject is abnormal, i.e. having lung disease; M is the total number of bins, here 100; ν j the center of the jth j

t

h bin; and β j = β ( ν j ) the coefficient. The sparse, piecewise constancy assumption on β entails that the β j ‘s are piecewise constant and sparse.

Get Radiology Tree app to read full this article<

The fused-lasso estimator

Get Radiology Tree app to read full this article<

l(b)=∑ni=1wi{yilog(pi)+(1−yi)log(1−pi)} l

(

b

)

=

i

=

1

n

w

i

{

y

i

log

(

p

i

)

+

(

1

y

i

)

log

(

1

p

i

)

}

where w i are known weights, y i is 1 if the ith i

t

h subject is abnormal and 0 otherwise, p is given by Equation (1) with V there replaced by that of the ith i

t

h subject, and n=17+52=69 n

=

17

+

52

=

69 , the total sample size. The weights are generally set to be identically 1, but to mitigate the unbalanced group sizes (17 normal scans vs. 52 abnormal scan), we set the weights to be proportionally 52 for each normal scan and 17 for each abnormal scan. The ML estimator is, however, generally neither sparse nor piecewise constant. So, a different approach of estimation is desirable.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

l(b)−λ1∑Mj=1|βj|−λ2∑Mj=2|βj−βj−1| l

(

b

)

λ

1

j

=

1

M

|

β

j

|

λ

2

j

=

2

M

|

β

j

β

j

1

|

where λ 1 and λ 2 are two non-negative tuning parameters to be determined by cross-validation; the penalized likelihood estimator is then obtained by maximization . The two tuning parameters effectively determine the degree of sparsity and piecewise constancy in the function estimate bˆ b

. For instance, when both tuning parameters are zero, the estimation becomes unconstrained ML estimation resulting in a generally nonsparse estimator that is not piecewise constant. On the other hand, for very large λ 1 ( λ 2 ), the β estimates will be mostly zero (approach a constant function with few jumps). Thus, the choice of the two tuning parameters is pivotal.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Misclassification Rates

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Information Pooling

Get Radiology Tree app to read full this article<

P(A|Vi​,i∈S)P(N|Vi​,i∈S)=∏i∈SP(A|Vi)P(N|Vi) P

(

A

|

V

i

,

i

S

)

P

(

N

|

V

i

,

i

S

)

=

i

S

P

(

A

|

V

i

)

P

(

N

|

V

i

)

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Results

Get Radiology Tree app to read full this article<

Figure 2, The plot of the coefficient estimates βˆ β^ as a function of HU; 15-mm peel depth and modification = T. The leftmost figure in the upper panel shows the fussed-lasso estimate of β j as a function of ν j in model (1) with V there being V 1 computed with data derived from lung images modified by radiologists, and peel depth was equal to 15 mm. Note that the function estimate is essentially a piecewise constant function that is non-zero over a small subset of (−1000 HU, 0 HU), as a result of the fussed-lasso penalty enforced in maximizing the objection function (3). Other figures in the upper panels, from left to right, are similar figures with V in model (1) replaced by Vi,i=2,3,4 Vi,i=2,3,4 . The figures in the lower panel are counterparts of V′i,i=1,2,3,4 V′i,i=1,2,3,4 . (Color version of figure is available online).

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Table 3

Misclassification Rates Based on Classification Using V′1 V

1

Depth, Modification In-sample Out-sample FPR FNR TFR FPR FNR TFR 15 mm, True 3/17 ≈ 0.18 15/52 ≈ 0.29 18/69 ≈ 0.26 4/17 ≈ 0.24 15/52 ≈ 0.29 19/69 ≈ 0.28 15 mm, False 3/17 ≈ 0.18 15/52 ≈ 0.29 18/69 ≈ 0.26 4/17 ≈ 0.24 15/52 ≈ 0.29 19/69 ≈ 0.28 65 mm, True 3/17 ≈ 0.18 15/52 ≈ 0.29 18/69 ≈ 0.26 4/17 ≈ 0.24 16/52 ≈ 0.31 20/69 ≈ 0.29 65 mm, False 4/17 ≈ 0.24 15/52 ≈ 0.29 19/69 ≈ 0.28 4/17 ≈ 0.24 16/52 ≈ 0.31 20/69 ≈ 0.29

FNR, false-negative rate; FPR, false-positive rate; TFR, total false rate.

Table 4

Misclassification Rates Based on Classification Using V′2 V

2

Depth, Modification In-sample Out-sample FPR FNR TFR FPR FNR TFR 15 mm, True 0/17 = 0 7/52 ≈ 0.13 7/69 ≈ 0.10 0/17 = 0 8/52 ≈ 0.15 8/69 ≈ 0.12 15 mm, False 0/17 = 0 6/52 ≈ 0.12 6/69 ≈ 0.09 0/17 = 0 7/52 ≈ 0.13 7/69 ≈ 0.10 65 mm, True 2/17 ≈ 0.12 9/52 ≈ 0.17 11/69 ≈ 0.16 2/17 ≈ 0.12 11/52 ≈ 0.21 13/69 ≈ 0.19 65 mm, False 2/17 ≈ 0.12 9/52 ≈ 0.17 11/69 ≈ 0.16 2/17 ≈ 0.12 14/52 ≈ 0.27 16/69 ≈ 0.23

FNR, false-negative rate; FPR, false-positive rate; TFR, total false rate.

Table 5

Misclassification Rates Based on Classification Using V′3 V

3

Depth, Modification In-sample Out-sample FPR FNR TFR FPR FNR TFR 15 mm, True 2/17 ≈ 0.12 6/52 ≈ 0.12 8/69 ≈ 0.12 3/17 ≈ 0.18 7/52 ≈ 0.13 10/69 ≈ 0.14 15 mm, False 2/17 ≈ 0.12 6/52 ≈ 0.12 8/69 ≈ 0.12 3/17 ≈ 0.18 7/52 ≈ 0.13 10/69 ≈ 0.14 65 mm, True 1/17 ≈ 0.06 8/52 ≈ 0.15 9/69 ≈ 0.13 1/17 ≈ 0.06 9/52 ≈ 0.17 10/69 ≈ 0.14 65 mm, False 1/17 ≈ 0.06 8/52 ≈ 0.15 9/69 ≈ 0.13 1/17 ≈ 0.06 8/52 ≈ 0.15 9/69 ≈ 0.13

FNR, false-negative rate; FPR, false-positive rate; TFR, total false rate.

Table 6

Misclassification rates based on classification using V′4 V

4

Depth, Modification In-sample Out-sample FPR FNR TFR FPR FNR TFR 15 mm, True 0/17 = 0 9/52 ≈ 0.17 9/69 ≈ 0.13 1/17 ≈ 0.06 10/52 ≈ 0.19 11/69 ≈ 0.16 15 mm, False 0/17 = 0 9/52 ≈ 0.17 9/69 ≈ 0.13 1/17 ≈ 0.06 10/52 ≈ 0.19 11/69 ≈ 0.16 65 mm, True 0/17 = 0 8/52 ≈ 0.15 8/69 ≈ 0.12 2/17 ≈ 0.12 8/52 ≈ 0.15 10/69 ≈ 0.14 65 mm, False 0/17 = 0 8/52 ≈ 0.15 8/69 ≈ 0.12 2/17 ≈ 0.12 8/52 ≈ 0.15 10/69 ≈ 0.14

FNR, false-negative rate; FPR, false-positive rate; TFR, total false rate.

Table 7

Misclassification Rates Based on Classification Using V′1​,V′2​,V′3 V

1

,

V

2

,

V

3 and V′4 V

4

Depth, Modification In-sample Out-sample FPR FNR TFR FPR FNR TFR 15 mm, True 0/17 = 0 3/52 ≈ 0.06 3/69 ≈ 0.04 0/17 = 0 4/52 ≈ 0.08 4/69 ≈ 0.06 15 mm, False 0/17 = 0 3/52 ≈ 0.06 3/69 ≈ 0.04 0/17 = 0 4/52 ≈ 0.08 4/69 ≈ 0.06 65 mm, True 0/17 = 0 5/52 ≈ 0.10 5/69 ≈ 0.07 1/17 ≈ 0.06 6/52 ≈ 0.12 7/69 ≈ 0.10 65 mm, False 0/17 = 0 5/52 ≈ 0.10 5/69 ≈ 0.07 1/17 ≈ 0.06 6/52 ≈ 0.12 7/69 ≈ 0.10

FNR, false-negative rate; FPR, false-positive rate; TFR, total false rate.

Table 8

Misclassification Rates Based on Classification Using V′2​,V′3 V

2

,

V

3 and V′4 V

4

Depth, Modification In-sample Out-sample FPR FNR TFR FPR FNR TFR 15 mm, True 0/17 = 0 3/52 ≈ 0.06 3/69 ≈ 0.04 0/17 = 0 3/52 ≈ 0.06 3/69 ≈ 0.04 15 mm, False 0/17 = 0 3/52 ≈ 0.06 3/69 ≈ 0.04 0/17 = 0 3/52 ≈ 0.06 3/69 ≈ 0.04 65 mm, True 0/17 = 0 4/52 ≈ 0.08 4/69 ≈ 0.06 1/17 ≈ 0.06 5/52 ≈ 0.10 6/69 ≈ 0.09 65 mm, False 0/17 = 0 4/52 ≈ 0.08 4/69 ≈ 0.06 1/17 ≈ 0.06 5/52 ≈ 0.10 6/69 ≈ 0.09

FNR, false-negative rate; FPR, false-positive rate; TFR, total false rate.

Table 9

Misclassification Rates Based on Classification Using V′2 V

2 , and V′3 V

3

Depth, Modification In-sample Out-sample FPR FNR TFR FPR FNR TFR 15 mm, True 0/17 = 0 4/52 ≈ 0.08 4/69 ≈ 0.06 0/17 = 0 6/52 ≈ 0.12 6/69 ≈ 0.09 15 mm, False 0/17 = 0 4/52 ≈ 0.08 4/69 ≈ 0.06 0/17 = 0 6/52 ≈ 0.12 6/69 ≈ 0.09 65 mm, True 2/17 ≈ 0.12 8/52 ≈ 0.15 10/69 ≈ 0.14 2/17 ≈ 0.12 9/52 ≈ 0.17 11/69 ≈ 0.16 65 mm, False 2/17 ≈ 0.12 8/52 ≈ 0.15 10/69 ≈ 0.14 2/17 ≈ 0.12 9/52 ≈ 0.17 11/69 ≈ 0.16

FNR, false-negative rate; FPR, false-positive rate; TFR, total false rate.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Discussion

Get Radiology Tree app to read full this article<

logit(p)≈constant+0.6×V′2(−1000,−900)+1.3×V′2(−400,−200). logit

(

p

)

constant

+

0.6

×

V

2

(

1000

,

900

)

+

1.3

×

V

2

(

400

,

200

)

.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

logit(p)≈constant−1.8×V′3(−1000,−950)+0.4×{V′3(−950,−900)−2×V′3(−900,−800)+V′3(−700,−600)} logit

(

p

)

constant

1.8

×

V

3

(

1000

,

950

)

+

0.4

×

{

V

3

(

950

,

900

)

2

×

V

3

(

900

,

800

)

+

V

3

(

700

,

600

)

}

where the term enclosed in curly brackets is equal to {V′3(−700,−600)−V′3(−900,−800)}−{V′3(−900,−800)−V′3(−950,−900)} {

V

3

(

700

,

600

)

V

3

(

900

,

800

)

}

{

V

3

(

900

,

800

)

V

3

(

950

,

900

)

} which can be loosely interpreted as curvature, i.e., the second derivative of the V′3 V

3 over the interval between −950 HU and −600 HU. Thus, the odds of abnormality then increases with reduced hyperaerated lung tissue from FRC to TLC, and it also increases with the second derivative of the V′3 V

3 over the interval between −950 HU and −600 HU. Indeed, Figure 1 shows that for normal scans, V′3 V

3 is generally a concave function between −950 HU and −600 HU, whereas it may become a convex function for abnormal scans. The physiological basis for the concavity in the normal population versus the convexity in some of the subjects with abnormal scans is an interesting future research problem. Figure 2 displays that the function estimate for the logistic regression model using V′4 V

4 implies the following model:

logit(p)≈constant−0.4×V′4(−1000,−900)+1.9×V′4(−800,−700)−V′4(−500,−450). logit

(

p

)

constant

0.4

×

V

4

(

1000

,

900

)

+

1.9

×

V

4

(

800

,

700

)

V

4

(

500

,

450

)

.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Acknowledgments

Get Radiology Tree app to read full this article<

Appendix

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

p(d|vi​,i∈S)×p(vi​,i∈S)=p(vi​,i∈S;d)=p(d)p(vi​,i∈S|d)=p(d)∏i∈Sp(vi|d)={p(d)}−|S|+1∏i∈Sp(vi​;d)={p(d)}−|S|+1∏i∈Sp(d|vi)p(vi) p

(

d

|

v

i

,

i

S

)

×

p

(

v

i

,

i

S

)

=

p

(

v

i

,

i

S

;

d

)

=

p

(

d

)

p

(

v

i

,

i

S

|

d

)

=

p

(

d

)

i

S

p

(

v

i

|

d

)

=

{

p

(

d

)

}

|

S

|

+

1

i

S

p

(

v

i

;

d

)

=

{

p

(

d

)

}

|

S

|

+

1

i

S

p

(

d

|

v

i

)

p

(

v

i

)

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

p(D=1|vi​,i∈S)p(D=0|vi​,i∈S)={p(D=0)}|S|−1{p(D=1)}|S|−1∏i∈Sp(D=1|vi)p(D=0|vi), p

(

D

=

1

|

v

i

,

i

S

)

p

(

D

=

0

|

v

i

,

i

S

)

=

{

p

(

D

=

0

)

}

|

S

|

1

{

p

(

D

=

1

)

}

|

S

|

1

i

S

p

(

D

=

1

|

v

i

)

p

(

D

=

0

|

v

i

)

,

which becomes ∏i∈Sp(D=1|vi)p(D=0|vi) ∏

i

S

p

(

D

=

1

|

v

i

)

p

(

D

=

0

|

v

i

) , the product of the posterior odds when the prior odds p(D=1)p(D=0)=1 p

(

D

=

1

)

p

(

D

=

0

)

=

1 .

Get Radiology Tree app to read full this article<

Appendix

Supplementary Material

Get Radiology Tree app to read full this article<

Figure S1

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure S2

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure S3

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure S4

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure S5

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

References

  • 1. Olson A.L., Swigris J.J., Lezotte D.C., et. al.: Mortality from pulmonary fibrosis increased in the United States from 1992 to 2003. Am J Respir Crit Care Med 2007; 176: pp. 277-284.

  • 2. Fletcher C., Oldham P.: Problem of consistent radiological diagnosis in coalminers’ pneumoconiosis: an experimental study. Br J Ind Med 1949; 6: pp. 168.

  • 3. Impivaara O., Zitting A.J., Kuusela T., et. al.: Observer variation in classifying chest radiographs for small lung opacities and pleural abnormalities in a population sample. Am J Ind Med 1998; 34: pp. 261-265.

  • 4. Mulloy K.B., Coultas D.B., Samet J.M.: Use of chest radiographs in epidemiological investigations of pneumoconioses. Br J Ind Med 1993; 50: pp. 273-275.

  • 5. Mikulski M.A., Hartley P.G., Sprince N.L., et. al.: Risk and significance of chest radiograph and pulmonary function abnormalities in an elderly cohort of former nuclear weapons workers. J Occup Environ Med 2011; 53: pp. 1046-1053.

  • 6. Muir D., Bernholz C.D., Morgan W., et. al.: Classification of chest radiographs for pneumoconiosis: a comparison of two methods of reading. Br J Ind Med 1992; 49: pp. 869.

  • 7. Musch D., Higgins I., Landis J.: Some factors influencing interobserver variation in classifying simple pneumoconiosis. Br J Ind Med 1985; 42: pp. 346-349.

  • 8. Yerushalmy J., Garland L., Harkness J., et. al.: An evaluation of the role of serial chest roentgenograms in estimating the progress of disease in patients with pulmonary tuberculosis. Am Rev Tuberc 1951; 64: pp. 225-248.

  • 9. International Labour Organization (ILO) : Guidelines for the use of the ILO international classification of radiographs of pneumoconioses.Occupational Safety and Health Series No. 222002.International Labour Office, ILOGeneva Rev 2000

  • 10. National Institute of Occupational Safety and Health. (NIOSH) : Issues in classification of chest radiographs. NIOSH.2010.Centers for Disease Control and Prevention (CDC)Atlanta, GA

  • 11. Busacker A., Newell J.D., Keefe T., et. al.: A multivariate analysis of risk factors for the air-trapping asthmatic phenotype as measured by quantitative CT analysis. Chest 2009; 135: pp. 48-56.

  • 12. Al Jarad N., Poulakis N., Pearson M., et. al.: Assessment of asbestos-induced pleural disease by computed tomography—correlation with chest radiograph and lung function. Respir Med 1991; 85: pp. 203-208.

  • 13. Bessis L., Callard P., Gotheil C., et. al.: High-resolution CT of parenchymal lung disease: precise correlation with histologic findings. Radiographics 1992; 12: pp. 45-58.

  • 14. Friedman A.C., Fiel S.B., Fisher M.S., et. al.: Asbestos-related pleural disease and asbestosis: a comparison of CT and chest radiography. AJR Am J Roentgenol 1988; 150: pp. 269-275.

  • 15. Harkin T.J., McGuinness G., Goldring R., et. al.: Differentiation of the ILO boundary chest roentgenograph (0/1 to 1/0) in asbestosis by high-resolution computed tomography scan, alveolitis, and respiratory impairment. J Occup Environ Med 1996; 38: pp. 46-52.

  • 16. Mathieson J., Mayo J., Staples C., et. al.: Chronic diffuse infiltrative lung disease: comparison of diagnostic accuracy of CT and chest radiography. Radiology 1989; 171: pp. 111-116.

  • 17. Ochsmann E., Carl T., Brand P., et. al.: Inter-reader variability in chest radiography and HRCT for the early detection of asbestos-related lung and pleural abnormalities in a cohort of 636 asbestos-exposed subjects. Int Arch Occup Environ Health 2010; 83: pp. 39-46.

  • 18. Flohr T.G., Schaller S., Stierstorfer K., et. al.: Multi–detector row CT systems and image-reconstruction techniques 1. Radiology 2005; 235: pp. 756-773.

  • 19. Kalra M.K., Maher M.M., Toth T.L., et. al.: Techniques and applications of automatic tube current modulation for CT 1. Radiology 2004; 233: pp. 649-657.

  • 20. Hoffman E.A., Lynch D.A., Barr R.G., et. al.: Pulmonary CT and MRI phenotypes that help explain chronic pulmonary obstruction disease pathophysiology and outcomes. J Magn Reson Imaging 2015; Epub ahead of print

  • 21. Guo J., Fuld M.K., Alford S.K., et. al.: Pulmonary Analysis Software Suite 9.0: Integrating quantitative measures of function with structural analyses. First International Workshop on Pulmonary Image Analysis2008.pp. 283-292.

  • 22. Hu S., Hoffman E.A., Reinhardt J.M.: Automatic lung segmentation for accurate quantitation of volumetric X-ray CT images. IEEE Trans Med Imaging 2001; 20: pp. 490-498.

  • 23. Tibshirani R., Saunders M., Rosset S., et. al.: Sparsity and smoothness via the fused lasso. J R Stat Soc Series B Stat Methodol 2005; 67: pp. 91-108.

  • 24. MATLAB and Statistics Toolbox Release 2013b, The MathWorks, Inc., Natick, Massachusetts, United States.

  • 25. Liu J., Ji S., Ye J.: SLEP: sparse learning with efficient projections. Arizona State University2009.

  • 26. Gattinoni L., Caironi P., Pelosi P., et. al.: What has computed tomography taught us about the acute respiratory distress syndrome?. Am J Respir Crit Care Med 2001; 164: pp. 1701-1711.

  • 27. Lynch D.A., Godwin J.D., Safrin S., et. al.: High-resolution computed tomography in idiopathic pulmonary fibrosis: diagnosis and prognosis. Am J Respir Crit Care Med 2005; 172: pp. 488-493.

  • 28. Hartley P.G., Galvin J.R., Hunninghake G.W., et. al.: High-resolution CT-derived measures of lung density are valid indexes of interstitial lung disease. J Appl Physiol 1994; 76: pp. 271-277.

  • 29. Yilmaz C., Watharkar S.S., de Leon A.D., et. al.: Quantification of regional interstitial lung disease from CT-derived fractional tissue volume: a lung tissue research consortium study. Acad Radiol 2011; 18: pp. 1014-1023.

This post is licensed under CC BY 4.0 by the author.