Influence of the measurement method of features in ultrasound images of the thyroid in the diagnosis of Hashimoto’s disease
© Koprowski et al.; licensee BioMed Central Ltd. 2012
Received: 22 October 2012
Accepted: 21 November 2012
Published: 28 November 2012
This paper shows the influence of a measurement method of features in the diagnosis of Hashimoto’s disease. Sensitivity of the algorithm to changes in the parameters of the ROI, namely shift, resizing and rotation, has been presented. The obtained results were also compared to the methods known from the literature in which decision trees or average gray level thresholding are used.
In the study, 288 images obtained from patients with Hashimoto’s disease and 236 images from healthy subjects have been analyzed. For each person, an ultrasound examination of the left and right thyroid lobe in transverse and longitudinal sections has been performed.
With the use of the developed algorithm, a discriminant analysis has been conducted for the following five options: linear, diaglinear, quadratic, diagquadratic and mahalanobis. The left and right thyroid lobes have been analyzed both together and separately in transverse and longitudinal sections. In addition, the algorithm enabled to analyze specificity and sensitivity as well as the impact of sensitivity of ROI shift, repositioning and rotation on the measured features.
Results and summary
The analysis has shown that the highest accuracy was obtained for the longitudinal section (LD) with the method of linear, yielding sensitivity = 76%, specificity = 95% and accuracy ACC = 84%. The conducted sensitivity assessment confirms that changes in the position and size of the ROI have little effect on sensitivity and specificity. The analysis of all cases, that is, images of the left and right thyroid lobes in transverse and longitudinal sections, has shown specificity ranging from 60% to 95% and sensitivity from 62% to 89%. Additionally, it was shown that the value of ACC for the method using decision trees as a classifier is equal to 84% for the analyzed data. Thresholding of average brightness of the ROI gave ACC equal to 76%.
The measurement of thyroid echogenicity is currently one of the most common and standardly performed measurements in ultrasound diagnosis. Measurements of this type have evolved over the years in accordance with progress and increase in the quality of ultrasound equipment. In the beginning [1–4], qualitative evaluation methods related to the areas of analysis and methods of description were explained. At that time, it was proved that normal thyroid echogenicity is higher than that of sternocleidomastoid and subhyoid muscles. Later, this approach was extended and the salivary gland was included in the analysis . With advances in computer technology and capabilities of digital recording and analysis, first papers on quantitative measurements [6–9] appeared. Those measurements were related to the use of basic methods of image analysis and processing in the diagnosis of, for example, Hashimoto’s disease [10, 11]. Due to imperfections introduced by the measurement method (scanning ultrasound pictures), this methodology has not been adopted in clinical practice. Scanning as well as other processes of non-digital image analysis introduce a significant error of the method and are not repeatable. The next stage were the methods of digital images analysis which ensured repeatability of measurement. They are mainly presented in Mailloux’s papers from the years 1984 to 1986 [12–14]. Those papers concern the application of texture analysis in ultrasound images. Nowadays, there are modern methods of analysis of ultrasound images. Although they are virtually limitless, there is still no clear method of disk image analysis that would give reproducible and unambiguous results. Many authors now attempt to use morphological and statistical methods in the analysis of texture of the thyroid lobe. In those methods, both the analysis of histograms, which gives partially correct results, as well as more advanced methods of texture analysis are used. These are, for example, methods [15, 16] which are based on the analysis of the areas indicated by the operator. The areas are analyzed by Co-occurrence Matrices. Then, Haralick’s coefficients are determined. The analysis of the Radon Domain  or Fuzzifying the Local Binary Patterns [18, 19] are further examples of the afore-mentioned methods. In recent works, an approach based on Support Vector Machines [20–22] can also be found. The results obtained using the Bayes classifier  or Gaussian mixture model  are interesting as well. In the literature, there are also other approaches to texture analysis, such as neural networks [25, 26] other [27–33] or dissertation . The methods of image analysis presented in those works need to be profiled to a specific application every time they are used. However, valuable evidence related to the measurement method and rough interpretation of ultrasound images of the thyroid arise from those works. For example, it was found that it is best to set the instrument to 10 MHz to achieve accuracy of results; the cut-off point is -69dB for Hashimoto’s disease [35, 36]. The authors of papers [37, 38] showed advanced methods of texture analysis of thyroid lobe images. Those methods were shaped to the diagnosis of Hashimoto’s disease. In paper , it was proved that only three of the ten features measured in an image are enough for a correct assessment of Hashimoto’s disease. These three features will be the basis of analysis in this paper.
59 healthy subjects aged 18 to 60,
73 patients with Hashimoto’s disease
The images were obtained with GE Logiq P5 ultrasound machine. The frequency of the transmitter was set to 10 MHz, and harmonic imaging option was turned off. All the images were recorded in DICOM format. During the test, the patient remained in the supine position and the doctor applied ultrasound heads to the right and left side of the thyroid.
For each subject, four ultrasound images were taken. Those were images of the right and left lobe of the thyroid in both transverse and longitudinal section. Due to thick errors caused by improperly performed acquisition, 288 images from patients with Hashimoto’s disease and 236 images from healthy subjects were further analyzed. The examined group was divided in equal proportions into learning, validation and test groups. Each ultrasound image was analyzed in great detail and, then, an expert physician selected for analysis a rectangular region (ROI) which covered the thyroid lobe in individual sections. Each time, the ROI included the greatest possible and most representative area of the patient’s thyroid lobe.
Preliminary image analysis
L GRAY input images were obtained from GE ultrasound machine with a resolution of M G ×N G =614×816 pixels. The first stage of image preprocessing was filtration done with the use of a median filter whose mask size is M h ×N h =3×3. The filtered images L MED were further used in subsequent stages of image analysis and processing. In the images (taken in transverse and longitudinal sections of the right and left thyroid lobe), an expert physician selected a rectangular area of analysis. Papers  and  describe an automated way of selecting this area of the thyroid, but only in transverse sections. The basis for its operation is a clearly visible artery calibrating the recognition system. The manually marked area of the thyroid lobe L S with a resolution of M s ×N s was analyzed. The results of the analysis are shown below.
The measured image features
The analysis of the thyroid lobe as texture in paper  proved that only 3 out of 10 different features are reliable in the assessment of Hashimoto’s disease. These features are: smoothness- w(1), minimum brightness after removing noise- w(2) and the percentage number of areas 8×8 in the square-tree decomposition- w(3). The ways to calculate individual values of the features are discussed in detail below:
Smoothness defined by the formula (1) is relatively easy to interpret because it is a standardized measure based on a standard deviation of the mean.
On the basis of pre-tests and preliminary analyses, a noise threshold of 0.2 was set. The value of i* formulated in this way constitutes another feature, i.e. w(2).
w(3) – percentage of instances of areas 8×8 obtained for the 10% threshold as a result of a square-tree decomposition.
The features w(1) to w(3) are the basis for further analysis.
A qualitative assessment of the measurement of echogenicity and its impact on the results obtained in the classification of Hashimoto’s disease was conducted using a statistical approach [41, 42]. A discriminant analysis was used for the following five options:
linear- linear discriminant analysis,
diaglinear- linear discriminant analysis but with a diagonal covariance matrix estimate (naive Bayes classifiers),
quadratic- quadratic discriminant analysis,
diagquadratic- quadratic discriminant analysis but with a diagonal covariance matrix estimate (naive Bayes classifiers),
mahalanobis- using the distance Mahalanobis with stratified covariance estimates.
It was assumed that the discriminatory variables w(1), w(2), w(3) represent a three-dimensional normal distribution (although previous studies carried out with the use of multivariate discriminant functions confirm the correctness of the classification, even in violation of this assumption). Divisibility of the variables is retained. This divisibility is reflected in the systematic difference in mean values between groups. Also the equality of covariance matrices is preserved. Empirical studies show that the assumption of equal group covariance matrices can be omitted.
These specific types of discriminant analysis were used to classify patients from healthy subjects. Assuming the classification results in terms of the following results: TP- true positive, TN- true negative, FP- false positive, FN- false negative, sensitivity was defined as TPR = TP / (TP + FN) and specificity as SPC = TN / (FP + TN). In addition, the analysis was performed for the following groups of data:
LO– images of the left transverse section of the thyroid,
RO– images of the right transverse section of the thyroid,
LRO– images of the left and right transverse sections of the thyroid,
RD– images of the right longitudinal section of the thyroid,
LD– images of the left longitudinal section of the thyroid,
RLD– images of the right and left longitudinal section of the thyroid,
RLOD– images of the right and left transverse and longitudinal sections of the thyroid.
Table showing the dependence of the results of specificity from the types of classification and the analyzed areas
Table showing the dependence of the results of sensitivity from the types of classification and the analyzed areas
Sensitivity to the change of parameters
The measured area (ROI), image Ls, underwent affine transformations in order to determine the dependence between the analyzed features w(1), w(2), w(3) and the size of the analyzed area as well as its position and rotation. The sensitivity analysis of these changes will be considered in subsequent sections.
sensitivity of features w(1), w(2) and w(3) to affine transformations of the ROI,
sensitivity of classification results to affine transformations of the ROI.
The analysis of changes in the value of w(1), w(2) and w(3) is important in this case because it points to their direct link with affine transformations (rotation, resizing and repositioning of the ROI). A direct comparison enables to assess the correctness of the formulation of features and their sensitivity to, for example, image rotation. This, in turn, enables to indicate which feature (and to what extent) depends on the position of the ultrasound head. It is also a condition to modify the formulation of a given feature so that it is only slightly dependent on the rotation.
Regardless of these results, the quality of the classification results for affine transformations -derived on the basis of all the features w(1), w(2) and w(3) – was assessed. The results demonstrate sensitivity of the algorithm which is considered as a measurement (diagnostic) method.
organs immediately adjacent to the thyroid lobe,
image borders - moved or enlarged ROI may not exceed the limits of the image,
ROI cannot be smaller than 10×10 pixels - this limitation is recognized in the definition of the coefficients w(1), w(2) and w(3).
Therefore, rounded values of the changes in the ROI position in the range of ±20 pixels of the ROI and its size of 10×10 to 90×90 were adopted. These values do not result in a breach of any of the above restrictions on the ROI for any of the analyzed images.
The only correct position and size of the ROI are determined by a specialist physician. Results and their impact on the value of accuracy will be observed (calculated) during ROI shifting, resizing or rotating.
The algorithm sensitivity to the resize of the marked area
The graphs (Figure 7) indicate that the feature w(3) is most sensitive to the change in the size of Ls. A small change in the number of rows and columns of the covered range indicated by the doctor influences considerably the value of the feature w(3). These changes concern tens of percent for the decrease or increase of the area by a few pixels. Changes to the feature w(1) are much milder and amount to several percent. When the area Ls increases, the value of the feature w(1) slightly changes. Only a significant reduction in the area Ls increases the error for the measurement of the feature w(1) up to 100%. Sensitivity of the feature w(3) to a change of the size of the area Ls looks completely different. The value of the feature varies by only a few percent in the full measured range. Discrete changes, visible on the graph, are due to the definition of the feature w(3) which is based on a histogram. Changes in the value of w(3) result from a change in the shape of the histogram. Therefore, the number of pixels of a given brightness and its proportion in relation to different brightness have to change. For the image Ls, it means that the texture changes. For this reason, changes in the size of the area Ls only slightly affect changes in the value of w(3).
The algorithm sensitivity to the change of the marked area position
It can be observed that values of the features w(1), w(2) and w(3) behave differently when the position of the Ls in the axes of rows and columns is changed. For the extreme positions, i.e. Δm=20, Δn=20 of the feature w(1), Δm=−20, Δn=−20 of the feature w(2) and Δm=−20, Δn=20 of the feature w(3), maximum values are achievable. Thus, a significant shift (more than 20 pixels) of the area Ls affects the results to a considerable extent. Globally, the feature w(3) is least sensitive to shifts of the area Ls.
The algorithm sensitivity to rotation around its own axis
The graph in Figure 9 shows that sensitivity to the rotation of the analyzed area is the highest for the feature w(3). The value of the feature w(2) changes slightly whereas the value of the feature w(1) changes oscillating. These oscillations result from the modification (due to rotation) of the Ls image content into new areas which contribute significantly to the value of STD and, therefore, to the value of the feature w(1).
In summary, the presented algorithm is least sensitive to the rotation of the area Ls and the feature w(2) is least sensitive to affine transformations (rotation, repositioning and resizing).
Assessment of the classification method sensitivity to affine transformations of the ROI
Assessment of sensitivity presented in the previous sections is determined on the basis of the results obtained from the individual features w(1), w(2) and w(3). These results are meaningful when the features are considered separately. However, in the case of the presented algorithm for classification, they form a coherent whole equally influencing the decision function. Therefore, it becomes legitimate to analyze sensitivity of the classification method to the presented affine transformations – ROI shifting and resizing (Ls). Ls image rotation will not be analyzed because, as it has been proved in previous sections, its influence on the results is negligibly small.
ROI shift in the range of ±10 pixels in the row or column axis slightly affects the results of specificity and sensitivity (changes of less than 0.05),
for a shift to the left or to the top by 10 pixels, SPC increases by approximately 0.03,
an increase in the size of the ROI by 13 pixels in rows or by 7 to 8 pixels in rows and columns causes a significant increase in specificity and sensitivity by approximately 0.03.
In conclusion, the choice of the area conducted by the expert and the algorithm are very good. ROI shifts in the range of ±10 pixels in the row or column axis as well as a decrease or increase in the ROI do not significantly affect the results. Therefore, the algorithm is resistant to fluctuations of the ROI (of both position and size) and its rotation (as demonstrated in the previous section).
Comparison with other results
In the literature described in the introduction [1–39], authors present several original methods of ultrasound image analysis. These methods are very interesting from the point of view of an ultrasound operator as they increase the accuracy and efficiency of diagnosis. Verification of sensitivity of the presented algorithms to changes in parametres, such as position, size and rotation of the ROI, is also an important feature for operators. This sensitivity analysis is important from the point of view of medical practice and interindividual variation. These elements may significantly influence the obtained results which testify to the quality of the algorithm. It may be that the advantage of one approach over the other forces highly accurate and precise indication of the ROI.
the histogram analysis of our algorithm fulfills a similar function as a classical analysis of the histogram described in paper . However, in that paper only one feature is taken into account, namely w(2) which is the minimum brightness, but after the removal of noise. Noise is defined as pixels whose sum is less than 20% of the calculated maximum amount of pixels.
the analysis of the features of our algorithm is similar to the analysis of another set of features (entropy, sum variance and mean value) presented in paper . Accuracy obtained there reaches 93.6%. However, the example given does not apply to Hashimoto’s disease.
comparison of methods of Co-occurrence matrix with the Radon transform and Muzzolini’s spatial features is shown in paper . However, the results shown do not relate directly to Hashimoto’s disease and do not analyze the impact of changes in the position of ROI on the obtained results.
simple analysis of the areas associated with Hashimoto’s disease is shown in paper . The results were obtained depending on the analysis method; sensitivity in the range of 71% to 88% and specificity in the range of 67% to 91%. These results are comparable with the results obtained with our algorithm, i.e. sensitivity 76% and specificity 95%. It should be noted that in the quoted paper , ROI areas were carefully selected by experts and some of the artifacts were manually eliminated.
In addition, the results of sensitivity, specificity and accuracy obtained from this discriminant analysis were compared in detail with other known methods [10, 12, 39]. Calculations were performed for the same group of 73 patients with Hashimoto’s disease and 59 healthy subjects. The images concerned only the left thyroid lobe in cross section (LD).
The following results were obtained:
Method 1: a classification method based on thresholding of mean values of brightness levels – sensitivity 92% and specificity 40%. ACC=76%,
Method 2: a method that uses decision trees described in paper  - sensitivity 88% and specificity 76% ACC=84% for a pruned decision tree,
Method 3: a discriminatory classification method proposed in this paper - sensitivity 76% and specificity 95%. ACC=84%.
The exact differences between the three methods are described in detail below.
Therefore, in the assessment of Hashimoto’s disease, more than one feature needs to be taken into account. Moreover, DICOM files should be analyzed directly and one of the two of the compared classifiers should be used (discriminant analysis or decision trees- method 3 and 2). Not only the absolute values of ACC but also the dynamics of their changes for small ROI displacements should be taken into account when analyzing the changes in results caused by ROI displacements.
Method 1. The first method is based on thresholding of echogenicity mean value (described in detail in  ). When applied to these data, it enables to obtain a result of ACC equal to 76% for the gray level threshold set to 25% of luminance (Figure 11). The range of average gray levels in the ROI for the analyzed cases was between 10% and 39% of saturation. Therefore, the graph shown in Figure 11 was carried out for different values of the threshold changed in the range of 13% to 36% in increments of 2.9% (assuming a step which is the tenth part of the range of 39%-10%). It can be observed that for the threshold value of 25%, shift of the ROI in the range of ±20 pixels affects significantly the value of accuracy– ACC changes by 7%. For the other threshold settings, the value of changes of Δm remains at a similar level, not exceeding 10%. In no sequence, a maximum for the value of Δm=0 is visible. Changes in the accuracy for different Δm do not have a well-defined direction of growth. Thus, it can be ultimately assumed that in the method of thresholding of echogenicity average level, ROI repositioning affects the result of accuracy to the extent of less than 10%.
Method 2: Another method uses decision trees (described in detail in  ). When applied to the collected data, it enables to obtain accuracy at 84%. In this case, accuracy variation was evaluated as a function of changes in ROI size and shift (Δm, Δn, ΔM, ΔN). The results for the pruned decision tree are shown in Figure 12a). The best tree is the one that has a residual variance that is no more than one standard error above the minimum value along the cross-validation line. Figure 12a) shows that changes in accuracy for changes in the values Δm, Δn, ΔM and ΔN are similar to the ones observed for echogenicity thresholding method and change by about 10%. The results also show a range of changes in the value of accuracy for each shift or resize of the ROI. The greater changes in the size or shift of the ROI are, the higher accuracy rate of change becomes. For example, for ΔM revised from the value of −16 to −15 pixels, the change in ACC reaches 16% (95%-79%). For ΔM as well as Δm, Δn and ΔN close to zero, ACC changes are smaller and reach the values of 5, 10%.
Narrowing the analysis to observation of accuracy changes only as a function of Δm, the impact of pruning the decision tree on the results is shown in Figure 12b). The degree of cutting the decision tree is dependent on the level ranging from 0 to 6 where level = 0 means no tree pruning. Trees are pruned based on an optimal pruning scheme that first pruned branches give less improvement in error cost. It can be seen that accuracy values vary depending on the degree of cutting the decision tree. Δm changes affect the value of accuracy by 5% for the first level values. When the decision tree is pruned too much, it loses its ability of classification and the error of accuracy reaches 30%.
This paper presents the influence of a measurement method of echogenicity in the diagnosis of Hashimoto’s disease, with a particular reference to the assessment of the algorithm sensitivity to a change in the ROI position. Classification was performed using a discriminant analysis for the following five options: linear, diaglinear, quadratic, diagquadratic and mahalanobis. Transverse and longitudinal sections of the thyroid right and left sides were analyzed. The analysis showed that the highest accuracy was obtained for the longitudinal section (LD) with the linear method, obtaining sensitivity = 76%, specificity = 95% and ACC = 84%. The impact of changes in the location of the ROI on the results was shown in one example and, separately, for all the analyzed cases. A change in the ROI position has the greatest impact on the value of features w(1) and w(3). The feature w(3) showed the greatest dependence on both the ROI position and also change of its size in the measured range of ±20 pixels. The percentage changes in the feature w(3) in the measured range Δm=±20 pixels and Δn =±20 pixels exceed 100%, while the changes of the feature w(2) amount to 5, 10%. The change in the value of w(1) is between 50% and 60%. The analysis of the results (mainly in Figure 9), confirms low dependence (below 30%) of any feature w(1), w(2) or w(3) on the ROI rotation in the range of 0 to 180o. A significant variation in the features w(3), w(2) or w(1) is not meaningful in relation to changes in sensitivity and specificity for the analyzed group of patients. Sensitivity assessment studies confirm that changes in the ROI position and size have little effect on sensitivity and specificity. SPC changes from 60% to 74% and TPR from 75% to 83% in the analysis of all cases of RLOD. Comparing the obtained results with other methods (method 1,2) known from the literature is also interesting. In the case of the classification method which uses decision trees  - method 2, the dynamics of ACC changes was at 15% (from 80% to 95%) for the full ROI displacement by Δm=±20. In the case of the method of thresholding (method 1) of echogenicity average levels, ACC was 76% for Δm=0 and the variation range of ACC was from 71% to 80% for Δm=±20 pixels.
The study carried out under KBN N518 047 32/3670 grant “Automatic analysis of a thyroid image as a tool used to make changes in the thyroid more objective”.
- Marwaha RK, Tandon N, Kanwar R, Ganie MA, Bhattacharya V, Reddy DH, Gopalakrishnan S, Aggarwal R, Grewal K, Ganguly SK, Mani K: Evaluation of the role of ultrasonography in diagnosis of autoimmune thyroiditis in goitrous children. Indian Pediatric 2008, 45(4):279–284.Google Scholar
- Marcocci C, Vitti P, Cetani F, Catalano F, Concetti R, Pinchera A: Thyroid ultrasonography helps to identify patients with diffuse lymphocytic thyroiditis who are prone to develop hypothyroidism. Clin Endocrinol Metab 1991, 72: 209–213. 10.1210/jcem-72-1-209View ArticleGoogle Scholar
- Hayashi N, Tamaki N: Sonography of Hashimoto’s thyroiditis. J Clin Ultrasound 1986, 14: 123–126. 10.1002/jcu.1870140208View ArticleGoogle Scholar
- Gutekunst R, Hafermann W, Mansky T, Scriba PC: Ultrasonography related to clinical and laboratory findings in lymphocytic thyroiditis. Acta Endocrinol 1989, 121: 129–135.Google Scholar
- Raber W, Gessl A, Nowotny P, Vierhapper H: Thyroid ultrasound versus antithyroid peroxidase antibody determination: a cohort study of four hundred fifty-one subjects. Thyroid 2002, 12(8):725–731. 10.1089/105072502760258712View ArticleGoogle Scholar
- Loy M, Cianchetti ME, Cardia F, Melis A, Boi F, Mariotti S: Correlation of computerized gray-scale sonographic findings with thyroid function and thyroid autoimmune activity in patients with Hashimoto’s thyroiditis. J Clin Ultrasound 2004, 32(3):136–140. 10.1002/jcu.20008View ArticleGoogle Scholar
- Mailloux G, Bertrand M, Stampfler R, Ethier S: Computer analysis of echographic textures in Hashimoto disease of the thyroid. J Clin Ultrasound 1986, 14(7):521–527. 10.1002/jcu.1870140705View ArticleGoogle Scholar
- Baldini M, Orsatti A, Bonfanti MT, Castagnone D, Cantalamessa L: Relationship between the Sonographic Appearance of the Thyroid and the Clinical Course and Autoimmune Activity of Graves’ Disease. J Clin Ultrasound 2005, 33(8):381–385. 10.1002/jcu.20157View ArticleGoogle Scholar
- Schiemann U, Gellner R, Riemann B, Schierbaum G, Menzel J, Domschke W, Hengst K: Standardized grey scale ultrasonography in Graves’ disease: correlation toautoimmune activity. Eur J Endocrinol 1999, 141(4):332–336. 10.1530/eje.0.1410332View ArticleGoogle Scholar
- Mazzotii G, Sorvillo F, Iorio S, Carbone A, Romeo A, Piscopo M, Capuano S, Capuano E, Amato G, Carella C: Grey-scale analysis allows a quantitative evaluation of thyroid echogenicity in the patients with Hashimoto’s thyroiditis. Clin Endocrinol 2003, 59: 223–229. 10.1046/j.1365-2265.2003.01829.xView ArticleGoogle Scholar
- Uwe S, Avenhaus W, Konturek J, Gellner R, Hengst K, Gross M: Relationship of clinical features and laboratory parameters to thyroid echogenicity measured by standardized grey scale ultrasonography in patients with Hashimoto’s thyroiditis. Med Sci Monit 2003, 9(4):49–53.Google Scholar
- Mailloux G, Bertrand M, Stampfler R, Ethier S: Computer analysis of echographic textures in Hashimoto disease of the thyroid. JCU J Clin Ultrasound 1986, 14(7):521–527. 10.1002/jcu.1870140705View ArticleGoogle Scholar
- Mailloux GE, Bertrand M, Stampfler R, Ethier S: Texture analysis of ultrasound B-mode images by segmentation. Ultrason Imaging 1984, 6(3):262–277. 10.1016/0161-7346(84)90012-9View ArticleGoogle Scholar
- Mailloux GE, Bertrand M, Stampfler R, Ethier S: Local histogram information content of ultrasound B-mode echographic texture. Ultrasound Med Biol 1985, 11(5):743–750. 10.1016/0301-5629(85)90108-5View ArticleGoogle Scholar
- Smutek D, Sara R, Sucharda P, Tjahjadi T, Svec M: Image texture analysis of sonograms in chronic inflammations of thyroid gland. Ultrasound Med Biol 2003, 29(11):1531–1543. 10.1016/S0301-5629(03)01049-4View ArticleGoogle Scholar
- Bastanfard M, Jalaeian B, Jafari S: Analysis of Sonogram Images of Thyroid Gland, Based on Wavelet Transform. World Academy of Science, Engineering and Technology; 2007:27–39.Google Scholar
- Savelonas MA, Iakovidis DK, Dimitropoulos N, Maroulis D: Computational Characterization of Thyroid Tissue in the Radon Domain. IEEE Int Symp Comput-Based Med Syst 2007, 33: 189–192.Google Scholar
- Dimitris KI, Eystratios GK, Maroulis D: Fuzzy Local Binary Patterns for Ultrasound Texture Characterization. Image Anal and Recognit Lecture Notes in Computer Sci 2008, 5112: 750–759. 10.1007/978-3-540-69812-8_74View ArticleGoogle Scholar
- Dimitris KI, Eystratios GK, Dimitris M: Fusion of fuzzy statistical distributions for classification of thyroid ultrasound patterns. Artif Intell Med 2010, 50(1):33–41. 10.1016/j.artmed.2010.04.004View ArticleGoogle Scholar
- Keramidas EG, Maroulis D, Iakovidis DK: A Thyroid Nodule Detection System for Analysis of Ultrasound Images and Videos. Springer Science+Business Media, J Med Syst 2012, 36(3):1271–1281.View ArticleGoogle Scholar
- Kotropoulos C, Pitas I: Segmentation of ultrasonic images using Support Vector Machines. Pattern Recognit Lett 2003, 24(4–5):715–727.View ArticleMATHGoogle Scholar
- Tsantis S, Glotsos D, Spyridonos P, Kalatzis G, Dimitropoulos N, Nikiforidis G, Cavouras D: Improving diagnostic accuracy in the classification of thyroid cancer by combining quantitative information extracted from both ultrasound and cytological images, 1st International Conference “From Scientific Computing to Computational Engineering” 1st IC-SCCE 8–10. Â© IC-SCCE, Athens; 2004:087–088.Google Scholar
- Šára R, Smutek D, Sucharda P, Svacina S: Systematic construction of texture features for Hashimoto’s lymphocytic thyroiditis recognition from sonographic images. AIME 2001, 2101: 339–348.MATHGoogle Scholar
- Tesa L, Smutek D 6th International PhD Workshop on Systems and Control. In Ultrasonography diagnostics using gaussian mixture model. Jozef Stefan Institute, Izola, Slovenia; 2005.Google Scholar
- Sujana H, Swarnamani S, Suresh S: Application of artificial neural networks for the classification of liver lesions by image texture parameters. Ultrasound Med Biol 1996, 22(9):1177–1181. 10.1016/S0301-5629(96)00144-5View ArticleGoogle Scholar
- Cavouras D, et al.: Computer image analysis of ultrasound images for discriminating and grading liver parenchyma disease employing a hierarchical decision tree scheme and the multilayer perceptron neural network classifier. In Med Inform Eur 1997, 2: 522–526.Google Scholar
- Bidi ABH A thesis submitted to the Cardiff University in candidature for the degree of Doctor of Philosophy, Institute of Medical Engineering and Medical Physics. In Image texture analysis of transvaginal ultrasound in monitoring ovarian cancer. Cardiff School of Engineering, Cardiff University; 2011.Google Scholar
- Chan K: Adaptation of ultrasound image texture characterization parameters. In Proc of the 20th Ann Int Conf of the IEEE Eng in Med and Biol Soc 1998, 2: 804–807.Google Scholar
- Keramidas EG, Iakovidis D, Maroulis D, Karkanis SA: Efficient and effective ultrasound image analysis scheme for thyroid nodule detection. Lect Notes Comput Sci 2007, 4633: 1052–1060. 10.1007/978-3-540-74260-9_93View ArticleGoogle Scholar
- Hirning T, et al.: Quantification and classification of echographic findings in the thyroid gland by computerized B-mode texture analysis. Eur J Radiol 1989, 9(4):244–247.Google Scholar
- Svec M, Sara R, Smutek D Research Report CTU–CMP–2003–22, Center for Machine Perception. In Sensitivity Analysis for Reproducibility of Ultrasound Image Classification. Czech Technical University; 2003.Google Scholar
- Gul K, Dirikoc A, Kiyak G, Ersoy PE, Ugras NS, Ersoy R, Cakir B: The Association Between Thyroid Carcinoma and Hashimoto’s Thyroiditis, The Ultrasonographic and Histopathologic Characteristics of Malignant Nodules. Thyroid 2010, 20(8):873–878. 10.1089/thy.2009.0118View ArticleGoogle Scholar
- Korzyńska A, Hoppe A, Strojny W, Wertheim D: Investigation of a combined texture and contour method for segmentation of light microscopy cell images. Proceedings of The Second IASTED International Conference on, Biomedical Engineering; 2004:234–239.Google Scholar
- Toufik EA Master’s thesis, Faculty of Electrical Engineering. In Automatic classification of the thyroid gland diseases by a histogram. Czech Technical University, Prague, Czech Republic; 2001.Google Scholar
- Zieleźnik W, Witkowska A, Małyszek-Tumidajewicz J, Koprowski R, Stępień B, Owczarek A, Koprowski R, Wójcik W: The new method of analysis thyroid echogenicity in the patients with Hashimoto’s disease. sent to Thyroid 2012.Google Scholar
- Zieleźnik W, Małyszek‐Tumidajewicz J, Stęchły T, Stępień B, Wójcik W, Owczarek : The analysis of thyroid ultrasound echogenicity in the patients with hashimoto’s disease. Eur Thyroid J 2011, 5: 160.Google Scholar
- Koprowski R, Wróbel Z, Zieleźnik W: Automatic ultrasound image analysis in Hashimoto’s disease. Lect Notes in Comput Sci - Adv in Pattern Recognit 2010, 6256: 98–106. 10.1007/978-3-642-15992-3_11View ArticleGoogle Scholar
- Koprowski R, Wróbel Z, Zieleźnik W: Analysis of thyroid ultrasonogram in Hashimoto’s disease, The 2010 International Congress on Computer Applications and Computational Science. IRAST, Singapore; 2010:448–451.Google Scholar
- Koprowski R, Wróbel Z, Zieleźnik W, Małyszek J, Witkowska A, Wójcik W: Relevance of features derived from ultrasound images of the thyroid in the diagnosis of Hashimoto’s disease. Biomed Eng Online 2012, 11: 48. 10.1186/1475-925X-11-48View ArticleGoogle Scholar
- Gonzalez RC, Woods RE: Digital Image Processing Using Matlab. Printed in United States of America, Prentice Hall; 2008.Google Scholar
- McLachlan GJ: Discriminant Analysis and Statistical Pattern Recognition. Wiley; 2004.MATHGoogle Scholar
- Banaka I, Kaltsas G, Antoniou S, Kanakis G, Zilos A, Baltas CS, Thomas D: Prognostic value of vascularity index for the diagnosis of autoimmune thyroid disease. JBR–BTR 2011, 94: 185–190.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.