Home Ultrasound Breast Tumor Image Computer-Aided Diagnosis With Texture and Morphological Features
Post
Cancel

Ultrasound Breast Tumor Image Computer-Aided Diagnosis With Texture and Morphological Features

Rationale and Objectives

Computer-aided diagnosis (CAD) systems based on shape analysis have been proved to be highly accurate in evaluating breast tumors. However, it takes considerable time to train the classifier and diagnose breast tumors, because extracting morphologic features require a lot of computation. Hence, to develop a highly accurate and quick CAD system, we combined the texture and morphologic features of ultrasound breast tumor imaging to evaluate breast tumors in this study.

Materials and Methods

This study evaluated 210 ultrasound breast tumor images, including 120 benign tumors and 90 malignant tumors. The breast tumors were segmented automatically by the level set method. The autocovariance texture features and solidity morphologic feature were extracted, and a support vector machine was used to identify the tumor as benign or malignant.

Results

The accuracy of the proposed diagnostic system for classifying breast tumors was 92.86%, the sensitivity was 94.44%, the specificity was 91.67%, the positive predictive value was 89.47%, and the negative predictive value was 95.65%. In addition, the proposed system reduced the training time compared to systems based only on the morphologic analysis.

Conclusions

The CAD system based on texture and morphologic analysis can differentiate benign from malignant breast tumors with high accuracy and short training time. It is therefore clinically useful to reduce the number of biopsies of benign lesions and offer a second reading to assist inexperienced physicians in avoiding misdiagnosis.

Because of changing lifestyles and the polluted environment, the mortality rate for malignant tumors has had the highest rank of major causes of death in recent years. Breast cancer is the most common cancer in women. According to 2003 statistics, 211,300 women are expected to be diagnosed with this disease, and only lung cancer causes more death in women ( ). Because of such a high incidence, breast cancer must be studied.

To reduce the mortality rate and extend patients’ lives, early detection and prompt treatment for breast cancer are very important. Detection of breast cancer usually consists of physical examination, imaging, and biopsy ( ). Although biopsy is the best way to accurately determine whether a tumor is benign or malignant, it is invasive and costs much more than other detection methodologies. Moreover, most biopsies are avoidable because the probability of positive findings at biopsy for cancer is very low, between 10 and 31% ( ). To avoid unnecessary biopsy, many researchers have investigated computer-aided diagnosis (CAD) systems based on medical imaging ( ). The aim of these CAD systems is to offer more objective evidence and increase the physician’s diagnostic confidence.

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<

Materials and methods

Data Acquisition

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<

Feature Extraction

Get Radiology Tree app to read full this article<

Morphologic Features

Get Radiology Tree app to read full this article<

Solidity=Convex_Area−Area∑Ni=1Convex_Areai−Areai/N S

o

l

i

d

i

t

y

=

C

o

n

v

e

x

_

A

r

e

a

A

r

e

a

i

=

1

N

C

o

n

v

e

x

_

A

r

e

a

i

A

r

e

a

i

/

N

where Area is the area of the tumor, Convex Area is the area of the convex hull of a tumor, and N is the number of tumors in the database. Figure 2 presents the convex hull of a tumor.

Figure 1, Two ultrasound breast tumor images: ( a ) benign tumor and ( b ) malignant tumor.

Figure 2, The convex hull and contour of a tumor.

Get Radiology Tree app to read full this article<

Texture Features

Get Radiology Tree app to read full this article<

γ(Δm,Δn)=A(Δm,Δn)A(0,0). γ

(

Δ

m

,

Δ

n

)

=

A

(

Δ

m

,

Δ

n

)

A

(

0

,

0

)

.

We define Equation (3) as follows:

A(Δm,Δn)=1Count(Δm,Δn)⋅∑M−1−Δmx=0∑N−1−Δny=0[fin(x,y)−f¯in][fin(x+Δm,y+Δn)−f¯in] A

(

Δ

m

,

Δ

n

)

=

1

C

o

u

n

t

(

Δ

m

,

Δ

n

)

x

=

0

M

1

Δ

m

y

=

0

N

1

Δ

n

[

f

i

n

(

x

,

y

)

f

¯

in

]

[

f

i

n

(

x

+

Δ

m

,

y

+

Δ

n

)

f

¯

i

n

]

where fin(x,y) f

i

n

(

x

,

y

) and fin(x+Δm,y+Δn) f

i

n

(

x

+

Δ

m

,

y

+

Δ

n

) are the gray levels of two pixels inside a tumor, f¯in f

¯

i

n is the mean value of fin(x,y) f

i

n

(

x

,

y

) , and Count(Δm,Δn) C

o

u

n

t

(

Δ

m

,

Δ

n

) is the number of a pair of pixels that are both inside a tumor; distance along both the x and y axes between the two pixels is Δ m and Δ n , respectively.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Tumor Segmentation

Get Radiology Tree app to read full this article<

Anisotropic Diffusion Filtering

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Stick Method

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Automatic Threshold Method

Get Radiology Tree app to read full this article<

Level Set Method

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<

(x,t=0)=±d (

x

,

t

=

0

)

=

±

d

where d is the distance from x to γ(t=0) γ

(

t

=

0

) , and the sign is chosen if the point x is outside (plus) or inside (minus) the initial hypersurface γ(t=0) γ

(

t

=

0

) . Thus, we have an initial function (x,t=0):RN→R (

x

,

t

=

0

)

:

N

ℜ with the property that

γ(t=0)=(x|(x,t=0)=0). γ

(

t

=

0

)

=

(

x

|

(

x

,

t

=

0

)

=

0

)

.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

(x(t),t)=0. (

x

(

t

)

,

t

)

=

0.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

t+∇(x(t),t)⋅x’(t)=0. t

+

(

x

(

t

)

,

t

)

x

(

t

)

=

0.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

t+F|∇|=0 t

+

F

|

|

=

0

with a given value of (x,t=0) (

x

,

t

=

0

) . This is the level set equation introduced by Osher and Sethian ( ).

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 3, Level set formulation of equation of motion. ( a ) and ( b ) show the curve γ and the corresponding surface φ(x,y) φ(x,y) at t = 0, and ( c ) and ( d ) show the curve γ and the corresponding surface φ(x,y) φ(x,y) at time t .

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Figure 4, The tumor segmentation with the proposed method. ( a ) Original image. ( b ) Result of the anisotropic diffusion filter on image ( a ). ( c ) Result of the stick method on image ( b ). ( d ) Result of the automatic threshold method on image ( c ). ( e ) Result of combining images ( c ) and ( d ). ( f ) Result of the level set method on image ( e ).

Get Radiology Tree app to read full this article<

Classifying with the SVM

Get Radiology Tree app to read full this article<

Figure 5, Separating hyperplane ( dashed lines , the margin).

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

W(α)=∑Ni=1αi−12∑Ni,j=1αiαjyiyjΦ(xi)⋅Φ(xj) W

(

α

)

=

i

=

1

N

α

i

1

2

i

,

j

=

1

N

α

i

α

j

y

i

y

j

Φ

(x

i

)

Φ

(

x

j

)

where {(xi,yi),i=1,…,N} {

(

x

i

,

y

i

)

,

i

=

1

,

,

N

} is the training example set S , each example xi∈Rn x

i

R

n belongs to a class labeled by yi∈{−1,1} y

i

{

1

,

1

} , and Φ(x) Φ

(

x

) denotes a mapping function that maps x into a high-dimensional feature space. The corresponding training examples (xi,yi) (

x

i

,

y

i

) with nonzero coefficients α i are called support vectors.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

W(α)=∑Ni=1αi−12∑Ni,j=1αiαjyiyjK(xi,xj) W

(

α

)

=

i

=

1

N

α

i

1

2

i

,

j

=

1

N

α

i

α

j

y

i

y

j

K

(

x

i

,x

j

)

where K is called a kernel function and must satisfy Mercer’s theorem ( ). Finally, the decision function becomes

f(x)=sgn(∑Ni=1αiyiK(xi,x)+b) f

(

x

)

=

sgn

(

i

=

1

N

α

i

y

i

K

(

x

i

,

x

)

+

b

)

where b∈R b

R is a constant.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Table 1

Four Common Types of Kernels Used in SVMs

Kernel typeK(xi,xj) K

(

x

i

,

x

j

) Linearx⋅z x

z Polynomial(γ⋅x⋅z+coef)d (

γ

x

z

+

c

o

e

f

)

d Gaussian radial basisexp(−γ⋅|x−z|2) exp

(

γ

|

x

z

|

2

) Sigmoidal neural networktanh(γ⋅x⋅z+coef) tanh

(

γ

x

z

+

c

o

e

f

)

Get Radiology Tree app to read full this article<

Results

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Table 2

Classification of Breast Tumors by the Proposed Method

Sonographic classification Benign ⁎ Malignant ⁎ Benign TN 110 FN 5 Malignant FP 10 TP 85 Total 120 90

TP, true positive; TN, true negative; FP, false positive; FN, false negative.

Get Radiology Tree app to read full this article<

Figure 6, ROC analysis of the proposed method. FPF, false-positive fraction; TPF, true-positive fraction.

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Table 3

Summary of Performance for Different Features

Item (a) (b) (c) (d) Accuracy 92.86% 88.57% 90.95% 83.33% Sensitivity 94.44% 85.56% 88.89% 85.56% Specificity 91.67% 90.83% 92.50% 81.67% PPV 89.47% 87.50% 89.89% 77.78% NPV 95.65% 89.34% 91.47% 88.29%

(a) The proposed method (combining morphological and texture features), (b) one morphological feature (which we used in the study), (c) six morphological features, and (d) texture features.

TP, true positive; TN, true negative; FP, false positive; FN, false negative.

Accuracy = (TP + TN)/(TP + TN + FP + FN)

Sensitivity = TP/(TP + FN)

Specificity = TN/(TN + FP)

Positive predictive value = TP/(TP + FP)

Negative predictive value = TN/(TN + FN)

Get Radiology Tree app to read full this article<

Get Radiology Tree app to read full this article<

Table 4

Time Comparison for Calculating Features of All Breast Tumor Images With Different Methods

Proposed method Morphologic Feature (one feature) Morphologic Features (six features) Time (sec) 112 66 812

Get Radiology Tree app to read full this article<

Discussion

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<

References

  • 1. Jemal A., Murray T., Samuels A., et. al.: Cancer statistics. CA Cancer J Clin 2003; 53: pp. 5-26.

  • 2. Esserman L.J.: New approaches to the imaging, diagnosis, and biopsy of breast lesions. Cancer J 2002; 8: pp. S1-S14.

  • 3. Bassett L.W., Liu T.H., Giuliano A.E., et. al.: The prevalence of carcinoma in palpable vs impalpable, mammographically detected lesions. AJR Am J Roentgenol 1991; 157: pp. 21-24.

  • 4. Rubin M., Horiuchi K., Joy N., et. al.: Use of fine needle aspiration for solid breast lesions is accurate and cost-effective. Am J Surg 1997; 174: pp. 694-696.

  • 5. Gisvold J.J., Martin J.K.: Prebiopsy localization of nonpalpable breast lesions. AJR Am J Roentgenol 1984; 143: pp. 477-481.

  • 6. Palmedo H., Grunwald F., Bender H., et. al.: Scintimammography with technetium-99m methoxyisobutylisonitrile: Comparison with mammography and magnetic resonance imaging. Eur J Nucl Med 1996; 23: pp. 940-946.

  • 7. Chang R.F., Wu W.J., Moon W.K., et. al.: Automatic ultrasound segmentation and morphology based diagnosis of solid breast tumors. Breast Cancer Res Treat 2005; 89: pp. 179-185.

  • 8. Chen D.R., Chang R.F., Huang Y.L.: Computer-aided diagnosis applied to US of solid breast nodules by using neural networks. Radiology 1999; 213: pp. 407-412.

  • 9. Stavros A.T., Thickman D., Rapp C.L., et. al.: Solid breast nodules: Use of sonography to distinguish between benign and malignant lesions. Radiology 1995; 196: pp. 123-134.

  • 10. Chang R.F., Wu W.J., Moon W.K., et. al.: Support vector machines for diagnosis of breast tumors on US images. Acad Radiol 2003; 10: pp. 189-197.

  • 11. Chang R.F., Wu W.J., Moon W.K., et. al.: Improvement in breast tumor discrimination by support vector machines and speckle-emphasis texture analysis. Ultrasound Med Biol 2003; 29: pp. 679-686.

  • 12. Burges J.C.: A tutorial on support vector machines for pattern recognition. Data Mining Knowledge Discov 1998; 2: pp. 121-167.

  • 13. Pontil M., Verri A.: Support vector machines for 3D object recognition. IEEE Trans Pattern Analy Machine Intell 1998; 20: pp. 637-646.

  • 14. Vapnik V.N.: The Nature of Statistical Learning Theory.1995.Springer-VerlagNew York

  • 15. Chapelle O., Haffner P., Vapnik V.N.: Support vector machines for histogram-based image classification. IEEE Trans Neural Networks 1999; 10: pp. 1055-1064.

  • 16. Heiden V.D., Duin R.P.W., de Ridder D., et. al.: Classification, parameter estimation and state estimation: An engineering approach using MATLAB.2004.WileyOxford

  • 17. Gerig G., Kubler O., Kikinis R., Jolesz F.A.: Nonlinear anisotropic filtering of MRI data. IEEE Trans Med Imaging 1992; 11: pp. 221-232.

  • 18. Perona J., Malik J.: Scale-scape and edge-detection using anisotropic diffusion. IEEE Trans Pattern Analy Machine Intell 1990; 12: pp. 629-639.

  • 19. Czerwinski R.N., Jones D.L., O’Brien W.D.: Edge detection in ultrasound speckle noise. IEEE Int Conf Image Processing 1994; 3: pp. 304-308.

  • 20. Otsu N.: A threshold selection method form gray-level histograms. IEEE Trans Systems Man Cyber 1979; 9: pp. 62-66.

  • 21. Sethian J.A.: Level Set Methods: Evolving Interfaces in Geometry, Fluid Mechanics, Computer Vision, and Materials Science.1996.Cambridge University PressCambridge

  • 22. Osher S., Sethian J.: Fronts propagating with curvature-dependent speed: Algorithms based on Hamilton-Jacobi equations. J Comput Phys 1988; 79: pp. 12-49.

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