| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Full paper |
1 Department of Obstetrics and Gynecology, 2 Department of Medical Physics, 3 Department of Radiology, School of Medicine, University of Patras, 265 00 Patras, Greece
Correspondence: Professor George Panayiotakis, Professor of Medical Physics, Department of Medical Physics, University of Patras, School of Medicine, Patras 265 00, Greece. E-mail: panayiot{at}upatras.gr
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Transvaginal sonography (TVS) represents the cornerstone of the available diagnostic procedures for excluding endometrial carcinoma, mainly due to its non-invasiveness and safe clinical profile. The incidence of suboptimally visualized endometria at TVS has diminished markedly over the past decades as sonographic equipment has evolved, particularly with the development of multifrequency endovaginal transducers and harmonic imaging [3]. While Doppler [4] and sonohysterography techniques [5] enhance endometrial assessment, in conventional grayscale TVS, measurement of double layer endometrial thickness remains the most widely accepted quantitative parameter. Despite the increased performance of TVS in excluding disease, it is considerably less accurate in establishing definite pathology in the presence of thickened endometrium [2].
Several authors have demonstrated the value of endometrial morphology in addition to the measurement of double layer endometrial thickness, particularly in the assessment of the 4–10 mm endometrial thickness "grey zone" [6, 7]. These observations are congruent with the perception that the combined consideration of morphological characteristics and endometrial thickness increases the specificity and negative predictive value of sonography [3].
With the objective of overcoming subjectivity induced by qualitative assessment, the implementation of automated techniques providing objective morphological characteristics, such as computerized texture analysis, would be beneficial, given the value of morphology in endometrial tissue interpretation.
In digital images, texture reflects tonal (intensities of image pixels) and structural (spatial distribution of pixel intensities) properties [8, 9]. Texture analysis refers to algorithms that quantify texture content that may, or may not, be perceived visually. As medical images capture various properties of biological structures, texture analysis of medical images can provide quantitative metrics relevant to the structure, morphology and status of biological tissues [10].
Texture-based classification schemes have been successfully implemented in a variety of ultrasound applications [11–15]. Computerized TVS assessment of endometrial morphology has been applied mainly in assisted reproduction techniques [16–19]; to the authors' knowledge, computerized texture analysis has not been implemented for diagnosing endometrial malignancies in grayscale TVS.
The aim of this study is to investigate the feasibility of computerized texture analysis in characterizing endometrial tissue as depicted in two-dimensional (2D) grayscale TVS images. Furthermore, the effect of a wavelet-based image processing technique in segmentation and subsequent characterization of endometrial tissue is investigated.
| Methods and materials |
|---|
|
|
|---|
Among the 82 recruited patients, 17 were excluded from the study for various reasons, including suboptimal visualization of the endometrium attributed to medioverted or retroverted uteri (n = 12), cancelled surgery (n = 3) and conservative surgery–myectomy (n = 2). Thus, the study population consisted of 65 women, with an average age of 57.2 years (range 41–80 years). 22 patients were perimenopausal (33.8%) and 43 patients (66.2%) were post-menopausal. Uterine perimenopausal or post-menopausal bleeding was the main indication for surgery (n = 49); other indications included uterine fibroids (n = 7), adnexal masses (n = 2), genital prolapse (n = 2), cervical intraepithelial lesions (n = 2) and miscellaneous (n = 3).
All 65 patients underwent a 2D transvaginal scan on a single ATL HDI 3500 ultrasonic imager (Advanced Technology Laboratories, Bothell, WA), with a multifrequency C 5–9 MHz transvaginal transducer, 48 h or less before the gynaecological operation (D&C or hysterectomy) that provided the histological diagnosis. For patients who underwent D&C prior to hysterectomy, only the TVS examination that was performed prior to D&C was considered. One experienced radiologist performed all examinations. For each patient, scan settings (dynamic range, gain, brightness, contrast and focus) were adjusted to obtain a clear depiction of the endometrium. Following an initial assessment of the entire endometrium in multiple planes, endometrial thickness was measured as the maximal double layer thickness in mid-sagittal section, at the thickest area of the endometrium near the fundus and from the outermost border of both sides of the endometrium. If a hypoechoic halo surrounded the endometrium, it was not included in endometrial thickness measurements. At least three measurements of endometrial thickness were performed, and a mean thickness was calculated. To assess possible alterations in endometrial morphology, additional evaluation of the uterus was undertaken. For subsequent analysis, ultrasonic images, in DICOM format, were digitally recorded with 8 bits pixel depth and 768x576 pixels size on a magneto-optical drive.
A collaborative reading by a radiologist and a gynaecologist, who were blinded to the histological diagnoses, was undertaken in order to select the most representative longitudinal image for each patient in the data set of images illustrating a distinctly visible endometrium from the internal cervical os to the fundus, and minimal shadowing artefacts.
Regarding surgical interventions, 24 patients underwent D&C, 40 abdominal hysterectomy and 1 vaginal hysterectomy. Histological diagnoses revealed malignant pathology in 15 women; in the remaining 50 women, benignity was documented. Concerning the malignant cases, 13 patients were diagnosed with endometrial adenocarcinomas, 1 with leiomyosarcoma and 1 with malignant mixed müllerian tumour – carcinosarcoma. No patient was diagnosed with atypical endometrial hyperplasia. Among benign cases, dysfunctional, proliferative or secretory endometria were detected in 8 patients (of whom 4 were perimenopausal), endometrial polyps in 17, simple or complex typical endometrial hyperplasia in 7 and endometrial atrophy in 18 patients. Fibroids were reported in 16 patients, while adenomyosis was diagnosed in 8 patients; 5 patients had additional benign cervical pathology.
The relationship between endometrial thickness and histological diagnoses is presented in Figure 1
. Of note is the high proportion of patients falling in the 5–10 mm "grey zone" group. As expected, the greater the endometrial thickness, the higher the incidence of endometrial cancer.
|
In this study, an image processing technique was employed to enhance the contrast of significant characteristics and to facilitate the segmentation procedure. Each image in the data set was processed with a locally adaptive wavelet-based technique [21, 22]. The method is based on local modification of multiscale gradient magnitude values provided by the redundant dyadic wavelet transform. Contrast enhancement is performed by applying a local linear mapping operator on multiscale gradient magnitude values. The first four frequency scales are used to enhance image characteristics of different size.
Figures 2 and 3
represent application examples of the locally adaptive wavelet technique on two images from the data set. Figure 2a
depicts an endometrial adenocarcinoma presenting with increased endometrial thickness, whereas Figure 3a
depicts derangement of the uterine contour by a leiomyosarcoma. Both corresponding processed images (
Figures 2b and 3b
) illustrate a significant amelioration of endometrial perception.
|
|
Endometrial tissue segmentation
A radiologist and a gynaecologist, in consensus, retraced manually two regions of interest (ROIs) on each original and processed image in the data set. They were both blinded to the histological diagnoses. The first region of interest (ROI1) corresponds to the endometrium without including the endometrial border, and the second (ROI2) comprises the endometrium along with the adjacent area of the myometrium (endometrium plus the layer containing the endometrial–myometrial interface). Specifically, ROI1 and ROI2 boundaries were delineated on either side of the endometrial–myometrial interface at approximately equal small distances. In order to investigate the textural properties of the layer that contains the endometrial–myometrial interface, a third region of interest (ROI3) was obtained by subtracting the two manually segmented ROIs (ROI2–ROI1).
Figure 4
depicts the manually segmented ROIs in the original (a) and the corresponding processed image (b). It is evident that the determined ROIs corresponding to the original and processed images are not identical.
|
were calculated. Normalization was performed by decreasing the grey level value of each pixel by the mean value µ and dividing it by the standard deviation
. The modified grey level values were rescaled in the range of [0, 255] (8 bits pixel depth) [25]. 32 textural features were extracted from each of the determined ROIs (ROI1, ROI2 and ROI3) on both original and processed images, employing first and second order statistics texture analysis algorithms [26–28]. The extracted features are capable of quantifying the inherent textural properties of the determined ROIs.
First order statistics textural features
First order statistics provide different statistical properties (first four statistical moments) of the intensity histogram of an image [26]. They depend only on individual pixel values and not on the interaction or co-occurrence of neighbouring pixel values. In this study, four first order textural features were calculated from each ROI: mean, standard deviation, skewness and kurtosis.
Second order statistics textural features
The Grey Level Co-occurrence Matrix (GLCM) is a well-established robust statistical tool for extracting second order texture information from images [27, 28]. The GLCM characterizes the spatial distribution of grey levels in the selected ROI. An element at location (i,j) of the GLCM signifies the joint probability density of the occurrence of grey levels i and j in a specified orientation
and specified distance d from each other. Thus, for different
and d values, different GLCMs result from each selected ROI.
In this study, four GLCMs, corresponding to four different directions (
= 0°, 45°, 90° and 135°) and one distance (d = 1 pixel), were computed for each selected ROI. 14 features were derived from each GLCM. Specifically, the features studied were angular second moment, contrast, correlation, variance, inverse difference moment, sum average, sum variance, sum entropy, entropy, difference entropy, information measure of correlation 1, information measure of correlation 2, shade and prominence. Four values were obtained for each feature corresponding to the four matrices. The mean and range of these 4 values were calculated, comprising a total of 28 second order textural features.
The 32 (first and second order statistics) textural features extracted from each ROI of the original and processed images were used as input for logistic regression analysis to determine the most discriminating subset of features for differentiation of benign from malignant endometrial tissue. The histological findings of benign or malignant diagnoses were considered as the truth table.
Segmentation reproducibility
In order to assess intraobserver reproducibility in the segmentation procedure, the 130 ROIs (65 ROI1 and 65 ROI2) were blindly delineated twice (after a 2 month period), collaboratively by a radiologist and a gynaecologist, on each original and processed image in the data set. Both radiologist and gynaecologist were blinded to the histological diagnoses. The area overlap criterion between the two segmented ROIs was used, defined as the ratio of intersection to the union of the two segmented areas [29]. Reproducibility was defined by the percentage of cases with an area overlap value exceeding 80%. To assess whether the differences between the two segmented areas were statistically significant, a non-parametric statistical test (two-tailed Wilcoxon signed ranks test [30]) with a significance level of 0.05 was applied.
Tissue characterization
Logistic regression analysis [31] was applied to determine the optimal subset of textural features (descriptors) capable of differentiating benign from malignant endometrial tissue. Specifically, the forward stepwise feature selection method of the logistic regression analysis was performed using a statistical software package (NCSS Statistical Software 2004, Kaysville, UT).
In order to assess the effect of the segmentation procedure (contribution of the various ROIs) in the endometrial tissue characterization task, separate logistic regression models were generated for the original images. Furthermore, aiming to investigate the effect of the processing technique in the characterization task, additional regression models were generated for the processed images.
For each data set (original and processed images), separate models were constructed for each ROI. Specifically, three regression models were constructed, with each one using as inputs the same 32 textural features extracted from ROI1, ROI2 and ROI3, respectively. An additional regression model was constructed combining the 32 textural features extracted from ROI1 and the 32 textural features extracted from ROI3. Hereafter, this model will be referred to as the one corresponding to (ROI1/ROI3).
The aforementioned models were constructed in the form of logit(
) equations (logarithm of the odds) as follows:
|
|
where Xs are the independent variables (features), v is the number of features (descriptors) in the model, Y is the binary dependent variable (pathology), which has two possible values, y1 (1: malignant) and y2 (0: benign),
0 is the intercept and
s are the logistic regression coefficients. From these logit(
) values, the estimated probability of malignancy (
) for tissue depicted in a ROI can be obtained from the inverse logistic transformation:
|
|
The performance of each logistic regression model was evaluated by means of ROC analysis using the area under the receiver operating characteristic (ROC) curve (Az) and the corresponding standard error (SE) [32]. Comparative evaluation was performed among models constructed of features extracted from: (a) original ROIs, (b) processed ROIs and (c) original and the corresponding processed ROIs. Furthermore, the sensitivity level of 93.3% (i.e. which correctly characterized 14 out of 15 malignant endometria) was selected to compare the specificity for each model.
| Results |
|---|
|
|
|---|
Tissue characterization
Table 1
provides the best logistic regression models, as chosen by forward stepwise selection of textural features extracted from ROI1, ROI2 and ROI3, as well as for the model combining the features from ROI1 and ROI3 (ROI1/ROI3), in terms of logit(
) equations for both original and processed images. All features entered in the regression models are statistically significant (Wald test for each regression coefficient, p<0.05). As observed in Table 1
, the features included in the regression models are differentiated with respect to the selected ROI, suggesting that each ROI is characterized by different inherent textural properties. Thus, different optimal subsets of textural features (descriptors) can be used to differentiate between malignant and benign tissue for both original and processed images. Furthermore, for both original and processed data sets, the regression model corresponding to the endometrium plus the layer containing the endometrial–myometrial interface (ROI2) is different from the regression model produced by combining features from the same two regions (ROI1/ROI3), suggesting that these two regions possess independent inherent textural properties. Additionally, the textural features included in regression models corresponding to each ROI are differentiated with respect to the original and processed images.
|
|
|
| Discussion |
|---|
|
|
|---|
The performance of the regression models with respect to tissue characterization was significantly improved when the wavelet-based processing technique was applied prior to the segmentation task. Specifically, the regression model produced by combining features extracted from ROI1/ROI3 of the processed images achieved the highest performance (Figure 5b
). This high performance may be attributed to: (a) more accurate delineation of ROIs obtained by endometrial contrast enhancement in the processed images; (b) the additionally revealed textural information provided by the processing technique yielding features with higher discriminating ability; (c) the hypothesis of two separate regions, corresponding to ROI1 and ROI3, capturing significant differences in textural information, resulting in improved endometrial characterization compared with merging the two regions into one (reflecting ROI2). As suggested by the proposed method, the endometrium and the endometrial–myometrial interface layer should be analysed separately as distinct regions in processed images. Combining textural features from the distinct regions provides improved performance in endometrial characterization.
It would be useful to compare the results obtained from this study with results reported aiming to increase TVS specificity. However, this might not be feasible because of inhomogeneities in case samples (size, patients' characteristics, etc.) used, the sonographic criteria assessed and variability in the statistical analyses performed [6, 7]. Moreover, there is no published study of computerized characterization of endometrial tissue, as depicted in TVS.
The proposed computerized approach for objective characterization of endometrial tissue has shown promising preliminary results. However, we are aware that our study design suffers from a number of shortcomings: both perimenopausal and post-menopausal women were recruited, and not all presented with vaginal bleeding. This inhomogeneity might reflect non-consistency in endometrial texture, yet the histological reports documented functional endometrium only in a few patients. Furthermore, surgical interventions were not uniform; D&C is undoubtedly inferior to hysterectomy in terms of accuracy of the histological diagnoses [34]. However, we have pooled heterogeneous data in accordance with other related studies [6, 35], due to the small size of case sample.
Additional work is required to improve the overall performance of the proposed method. Efforts will focus on assessing the performance of the proposed method in newer generation ultrasound scanners and larger data sets. The processing technique used in this work was found to be fundamental to the task of tissue characterization. To improve further the image processing technique, a speckle reduction [36] method should be considered prior to contrast enhancement. Finally, we aim to investigate the potential increase in the diagnostic accuracy of TVS by combining computerized textural features with shape features for capturing irregularity variability in the endometrial–myometrial interface.
| Acknowledgments |
|---|
Received for publication July 27, 2006. Revision received November 15, 2005. Accepted for publication November 23, 2006.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| BJR | DMFR | IMAGING | ALL BIR JOURNALS |