|
|
||||||||
Gastrointestinal Imaging |
1 From the Diagnostic Radiology Unit and Center for Anatomical, Functional and Molecular Imaging Research (L.H., I.L.A., L.A., B.E.V.B.) and Department of Surgery (J.L.), Université Catholique de Louvain, Saint-Luc University Hospital, Ave Hippocrate 10, B-1200 Brussels, Belgium; and Departments of Biomedical Engineering (Z.C.) and Electrical Engineering and Computer Science (B.M.D.), Vanderbilt University, Nashville, Tenn. Received November 7, 2003; revision requested January 28, 2004; final revision received April 13; accepted May 12. Supported in part by grant 3.4578.00 from the Fonds National de la Recherche Scientifique (Belgium) and by NIH grant 4R33CA09135203. Address correspondence to L.H. (e-mail: hermoye@rdgn.ucl.ac.be).
| ABSTRACT |
|---|
|
|
|---|
MATERIALS AND METHODS: The institutional review board approved this retrospective study and waived the requirement for informed consent. The semiautomatic segmentation algorithm is based on geometric deformable models and the level-set technique. It entails (a) placing initialization circle(s) on each image section, (b) running the algorithm, (c) inspecting and possibly manually modifying the contours obtained with the segmentation algorithm, and (d) placing lines to separate the liver segments. For 18 living donors (eight men and 10 women; mean age, 34 years; age range, 2546 years), two observers each performed two semiautomatic and two manual segmentations on contrast materialenhanced T1-weighted MR images. Each measurement was timed. Actual graft weight was measured during surgery. The time needed for manual and that needed for semiautomatic segmentation were compared. Accuracy and repeatability were evaluated with the Bland-Altman method.
RESULTS: Mean interaction time was reduced from 25 minutes with manual segmentation to 5 minutes with semiautomatic segmentation. The mean total time for the semiautomatic process was 7 minutes 20 seconds. Differences between the actual volume and the estimated volume ranged from 223 to +123 mL for manual segmentation and from 214 to +86 mL for semiautomatic segmentation. The 95% limits of agreement for the ratio of actual graft volume to estimated graft volume were 0.686 and 1.601 for semiautomatic segmentation and 0.651 and 1.957 for manual segmentation. Semiautomatic segmentation improved estimation in 15 of 18 cases. Inter- and intraobserver repeatability was higher with semiautomatic segmentation.
CONCLUSION: Use of the semiautomatic segmentation algorithm substantially reduces the time needed for volumetric measurement of liver segments while improving both accuracy and repeatability.
© RSNA, 2004
| INTRODUCTION |
|---|
|
|
|---|
In particular, the volume of the potential graft must be measured to ensure sufficient liver function after surgery. In the recipient, the ratio between the graft size and the recipients body weight ideally must be higher than 0.8%1% (5) or the ratio between the graft size and the standard liver volume (which is calculated according to surface area) (6) must be higher than 50% (7). In the donor, a remnant liver volume of 30% of the standard liver volume is usually considered to be the lower limit (8). Lower values are associatedwith higher morbidity and mortality in both the recipient and the donor owing to the small-for-size syndrome. In pediatric patients weighing less than 15 kg, a left lateral segment (Couinaud [9,10] segments II and III) graft is sufficient. In adults, a right lobe (Couinaud segments VVIII) graft is usually necessary because of volume constraints in the recipient.
Preoperative liver segmentation has proved useful for measuring the graft volume before living donor liver transplantations in previous studies (1114). In these studies, the liver segments were manually delineated on each image section. The delineated areas were multiplied by the section thickness to obtain volumes and summed to obtain the total volume of the liver segments. This process is tedious and time consuming. To compensate for this problem, automatic segmentation techniques have been proposed (15,16). These methods involve the use of sequences of thresholding, morphologic operations (ie, mathematic operations, such as image dilation, erosion, opening, and closing, that are based on shape), and model deformations. These techniques are complex and require long computation times. Because these algorithms were developed for CT, they are usually not directly usable with MR images, which are increasingly used for the preoperative screening of potential liver donors (3,13). Furthermore, these studies were focused on three-dimensional visualization. The actual organ volume was not used as the standard of reference.
The purpose of our study was to compare the accuracy and repeatability of a semiautomatic segmentation algorithm with those of manual segmentation for determining liver volume in living liver transplant donors at MR imaging.
| MATERIALS AND METHODS |
|---|
|
|
|---|
MR Imaging
Abdominal images were acquired during preoperative screening by using a 1.5-T whole-body MR imaging unit (Gyroscan NT Intera T15; Philips Medical Systems, Best, the Netherlands). To measure the liver volumes, a contrast materialenhanced T1-weighted fast-field-echo sequence was used. Acquisition parameters were as follows: field of view, 400 x 320 mm; acquisition matrix, 304 x 194; reconstruction matrix, 512 x 512; section thickness, 7 mm; section gap, 13 mm; repetition time msec/echo time msec, 135/2.9; flip angle, 60°; and acquisition time, 14 seconds per set of images. Twenty sections were acquired to cover the whole liver. The echo time was chosen to create a phase cancellation between fat and water signals. This creates a black outline at the border of the liver in which the voxels contain both liver tissue and peritoneal fat. This effect is beneficial for segmentation. A 0.1-mmol dose of gadodiamide (Omniscan; Amersham Health, Cork, Ireland) per kilogram of body weight was injected at a rate of 2.5 mL/sec with a power injector (Spectris; Medrad, Indianola, Pa). Four sets of images were acquired before administration of the contrast agent and during the arterial, portal venous, and delayed phases. Only the images acquired during the portal venous phase were used for segmentation.
Manual Segmentation
Two independent trained observers (L.H. and I.L.A., who both had 2 years of experience with liver MR imaging) each performed two manual segmentations retrospectively for each patient. To prevent bias, the observers were blinded to the results of their first measurement, the results of the other observer, the results of the other segmentation method, and the actual graft volume. The measurements were performed in random order, and successive measurements in the same patient were performed after a time interval of at least 2 weeks. Manual segmentation was performed at a GE Advantage Windows workstation (GE Medical Systems, Milwaukee, Wis).
The method that was described by Heymsfield et al (17) for measuring the volume of the whole liver and adapted by Kawasaki et al (12) for measuring the segmental volumes was used. On each section, the outlines of the caudate lobe (segment I), the left lateral segment (segments II and III), the left medial segment (segment IV), and the right lobe (segments VVIII) were manually drawn by using the computer mouse (Fig 1). The anatomic landmarks used to separate the liver segments were the gallbladder and the middle hepatic vein for the right lobeleft medial segment separation; the falciform ligament for the left medial segmentleft lateral segment separation; and the inferior vena cava, portal bifurcation, and ligamentum venosum for the delineation of the caudate lobe (18). For sections on which the anatomic landmarks were not visible, the separation lines were mentally interpolated. The volume corresponding to each outline was obtained by multiplying the area of the outline by the section thickness (including the gap). The total volumes of the liver segments were obtained by summing the volume on each section. For each measurement, the total segmentation time was recorded.
|
|
Specific features of the algorithm used herein include a speed function called the accumulative speed function (19), which permits stopping the progression of the contours at weak boundaries, automatic estimation of the parameters required by this algorithm on the basis of the initial circles, and a modification of the algorithm initially proposed by Pan and Dawant to prevent early contour stopping caused by the presence of the hepatic vessels that appear as very bright structures in our images. The overall segmentation process involves the following steps: (a) placing initialization circle(s) on each section, (b) running the algorithm, and (c) inspecting and possibly manually modifying the contours obtained with the segmentation algorithm.
With the segmentation algorithm, two lines must be drawn to separate the liver into segments. As in the manual segmentation process, the liver is separated into the caudate lobe, left lateral segment, left medial segment, and right lobe (Fig 1). Separation is performed by drawing lines on a few sections on which the landmarks are visible. These lines are then automatically interpolated to create surfaces that separate the liver into segments. The caudate lobe must be manually delineated. Each measurement was timed. The initialization time (ie, the time required to place the initialization circles), the computing time (ie, the time required to automatically delineate the liver contours), and the time for manual corrections and delineation of the liver segments were recorded separately. The total time was the sum of these three times. The interaction time was the sum of the initialization time and the time for manual corrections and delineation of the liver segments.
Surgical Procedure
The recipients received a left lateral segment (segments II and III), left lobe (segments IIIV), or right lobe (segments VVIII) graft (25). Intraoperative ultrasonography was used in the donors to delineate the intrahepatic vascular anatomy. The arterial, venous, and biliary elements of the corresponding graft were isolated. Parenchymal transection was performed without inflow occlusion. The liver graft was flushed successively with Hartmann (4°C) and University of Wisconsin solutions through a portal vein cannula. The corresponding hepatic artery was then ligated, and the corresponding hepatic vein was transected and sutured progressively. The graft was weighed by using a calibrated scale after the flushing and the preparation work. Graft weight was converted into graft volume by using an assumed unit density for the liver (26).
Effect of Section Thickness
The influence of the section thickness on the volumes measured with the semiautomatic segmentation algorithm was assessed in five patients with normal livers. In each patient, the contrast-enhanced T1-weighted fast-field-echo sequence used for preoperative screening of the living donors was performed three times with increasing section thicknesses (5, 7, and 10 mm). The section gap (1 mm) was not modified. For the images acquired with each sequence, one observer (L.H.) performed two semiautomatic segmentations as described above. The estimated volumes of the caudate lobe (segment I), the left lateral segment (segments II and III), the left medial segment (segment IV), the right lobe (segments VVIII), and the total liver, as averaged across the five patients, are reported for the three section thicknesses. The influences of section thickness on the number of sections necessary to cover the whole liver and on the interaction and total segmentation times were also studied.
Statistical Analysis
The total time for manual segmentation was compared with the total and interaction times for semiautomatic segmentation. For the semiautomatic segmentation method, the percentages of the total time devoted to initialization, computing, and manual corrections and delineation of the liver segments were calculated.
For accuracy comparisons, the mean graft volumes ± standard deviations measured during surgery and the mean volumes (ie, the mean of the four observations [two observations by each of two observers]) estimated with the manual and semiautomatic segmentations were averaged across all of the donors. In addition, for each donor, the difference between the actual graft volume (measured during surgery) and the estimated volume (the mean of the four observations) was calculated separately for the manual and semiautomatic segmentations. These differences were plotted against the averages of the actual and the estimated volumes. Scatterplots were drawn to illustrate the relationships between the actual and estimated graft volumes. The identity line was indicated on these plots. The number of cases in which use of the semiautomatic segmentation algorithm improved the accuracy of the estimation was counted.
The statistical analysis was performed with the method described by Bland and Altman (27,28). The 95% limits of agreement were calculated separately for the actual volume versus the volume determined with manual segmentation and the actual volume versus the volume determined with semiautomatic segmentation. Because the measurement error was proportional to the mean, a logarithmic transformation (natural logarithm) was used (29). Independence after the transformation was checked by using a Bland and Altman plot (27,28). Normality of the differences was checked with the Kolmogorov-Smirnov test. Because for each segmentation method the mean of four measurements was used, the variance of the differences was corrected for repeated measurements (28). The within-subject variance was calculated with one-way analysis of variance. The 95% limits of agreement were calculated as the mean difference ± 1.96 (adjusted standard deviation of the differences). The limits of agreement expressed in a logarithmic scale were transformed back into being expressed in a natural scale by using an exponential function. Because the calculations were performed with the log-transformed quantities, these limits of agreement must be interpreted in terms of ratios.
The Bland and Altman repeatability coefficient (28,30) was used as a measurement of inter- and intraobserver agreement for the manual and semiautomatic segmentations. The calculations were performed separately for the caudate lobe (segment I), the left lateral segment (segments II and III), the left medial segment (segment IV), the right lobe (segments VVIII), and the whole liver. Because for each of these tests the range of values was narrow, logarithmic transformation was not necessary. For the measurement of interobserver agreement, only the first measurement of each observer was taken into account (31). Intraobserver agreement was measured for both observers. The repeatability coefficient was calculated as 2.77 times the within-subject standard deviation. By definition, the measurement error is smaller than the repeatability coefficient for 95% of the observations. Therefore, a smaller repeatability coefficient indicates better repeatability of the method.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
Repeatability Measurements
The repeatability coefficients for both the manual and the semiautomatic segmentation methods are reported in Table 4. Lower repeatability coefficients (ie, better repeatability) were obtained with the semiautomatic segmentation method for all segments except the caudate lobe (segment I).
|
|
| DISCUSSION |
|---|
|
|
|---|
The main advantage of the semiautomated method is the time savings it confers. Initialization time might be further reduced by using three-dimensional segmentation. In this case, initialization on one or two sections would be sufficient. Automatic initialization could even be considered. However, three-dimensional segmentation would increase the computation time. Furthermore, its accuracy would have to be evaluated. Manual corrections and placing of separation lines are the most time-consuming steps. If the initialization circles are judiciously placed, manual corrections are only necessary at the borders, at which the contrast with the surrounding organs is very weak. This occurs especially in lean patients because in other patients, phase cancellation between the liver and the surrounding fat creates a black line surrounding the liver.
The improvement in accuracy can be explained by the more precise delineation of the liver contours. The accuracy of manual segmentation is indeed limited by the observers dexterity. If the automatic segmentation algorithm fails, manual corrections can quickly remove the mistakes. The positioning of lines on the anatomic landmarks that separate the segments and the automatic interpolation of these lines on the sections on which the landmarks are not visible are also favorable factors. Furthermore, the apparent border of the liver depends on the intensity window center and width. The manual adjustment of image contrast can thus influence the volume measured during manual segmentation (32), whereas automatic segmentation is independent of this setting. Nevertheless, even with semiautomated segmentation, a factor in the range of 0.6861.601 can still separate the estimated volume from the actual graft size. In adult grafts (whether from the right or the left lobe), an overestimation and an underestimation of graft size would lead to a higher risk of small-for-size syndrome in the recipient and in the donor, respectively.
Several reasons can explain the discrepancies observed between the estimated and the real volumes. The surgical context may impose the necessity of using slightly different anatomic landmarks to separate the liver segments. Furthermore, during the preparation of the graft, the blood vessels are flushed. Therefore, the blood vessels included in the segmented liver could lead to an overestimation of the graft size, on the order of 29% (33). Because this effect is more important for right- or left-lobe grafts, this may explain the systematic overestimation observed for these grafts. Partial volume effects related to the stomach and the heart are often present in upper sections (12). On these sections, the delineation of the liver must be based on both the image contrast and the contours on the adjacent sections.
Section thickness is also a limiting factor. With 7-mm sections and 3-mm gaps, the liver border is averaged over 10 mm. As shown in the results of our study, the left lateral segment, the right lobe, and the whole liver are the areas whose measurements are most affected by the section thickness. These liver parts have some oblique interfaces (relative to the imaging plane) with the surrounding organs (Fig 1) and are therefore sensitive to partial volume effects. In contrast, the borders of the caudate lobe and the left medial segment are more perpendicular to the imaging plane and are less sensitive to partial volume effects. Similar results were reported previously (34). It was shown that estimations of the volume of the bladder, whose shape varies along the z-axis, were higher with thinner sections, whereas the volume of the rectum, a cylindric organ that is quasiperpendicular to the imaging plane, was insensitive to the thickness of transverse sections.
In the present study, the actual and estimated graft volumes were similar to those reported in previous studies (12,35). The differences between the actual and the estimated graft volumes were in the normal range (11). Because most previous studies (11,12,35) involved the use of correlation coefficients and regression analysis to quantify the relationship between actual and estimated graft volumes, comparison of the 95% limits of agreement is difficult. Because the use of correlation coefficients and regression analysis is debated in the context of methods for comparing studies (27), their use was avoided. The main drawback of correlation coefficients and regression analysis relates to the fact that they measure the strength of the linear relationship between two variables, not the agreement between them. Furthermore, they depend on the scale and range of the measurements. Data that seem to be in poor agreement can produce quite high correlations (27). The statistical methods used in this study are similar to those used by Bakker et al (31) for renal volume measurements.
The use of automatic processes always increases repeatability (36). As expected, use of the semiautomatic segmentation method led to higher intra- and interobserver repeatability (except for the caudate lobe, for which the segmentation was manual in both methods). Manual determination of the liver contours, mental interpolation of the separations between the segments, and dependence on the display window settings (32) can explain the lower repeatability achieved with manual segmentation. However, even with semiautomatic segmentation, several variability sources exist. The manual placement of the initialization circles slightly modifies the final segmentation. The best results are obtained when two or three initialization circles are placed in representative areas on each section to take variations in the hepatic signal intensity into account. Manual corrections are left to the observers own judgment. Interobserver variability can result, especially in sections in which partial volume effects are present. The positioning of the separation lines between the liver segments can also introduce variability in patients in whom the anatomic landmarks are unclear.
The study had several limitations. Only 18 donors were included. This comparison could be performed with a larger number of donors within the framework of a multicenter study. The effects of the inclusion of blood vessels in the segmented liver and those of partial volume effects and section thickness have already been discussed. Moreover, in this study, the actual graft volume was inferred from its weight (11,12,35) while assuming a unit density for the liver (26). Water displacement (37) would allow direct measurement of the graft volume but would be cumbersome in clinical practice.
Finally, the accuracy of determining the segmental and subsegmental boundaries with straight lines has been questioned (18,38). The true anatomic limits of the liver subsegments are determined by the portal anatomy, which shows large variations (38). To take this variability into account, the use of liver segmentation algorithms based on the portal anatomy has been proposed (15,39). However, these methods are complex and are not widely used. In addition, most previous analyses of discrepancies between radiologic determinations of the segmental anatomy and the actual anatomy of liver segments and subsegments have concerned the limits between the subsegments of the right lobe (38,40,41). In contrast, the limits used in the present studynamely, the limits between the left lateral segment and the left medial segment and those between the left and the right lobeconform fairly to vertical planes (42).
In conclusion, in this validation study performed with living donors for liver transplantation, the use of semiautomatic liver segmentation saved a substantial amount of time while improving both accuracy and repeatability. Other potential applications of the technique include the determination of segmental volumes before partial hepatectomy for tumors (43) and the monitoring of liver regeneration (44).
| FOOTNOTES |
|---|
Author contributions: Guarantors of integrity of entire study, L.H., B.E.V.B.; study concepts, L.H., Z.C., J.L., B.M.D., B.E.V.B.; study design, L.H., I.L.A., B.M.D., B.E.V.B.; literature research, L.H., I.L.A., J.L., B.E.V.B.; clinical studies, L.H., I.L.A.; data acquisition, L.H., I.L.A., Z.C., J.L., L.A.; data analysis/interpretation, L.H., I.L.A., J.L., L.A., B.E.V.B.; statistical analysis, L.H.; manuscript preparation and editing, L.H.; manuscript definition of intellectual content, revision/review, and final version approval, all authors
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Y. Nakayama, Q. Li, S. Katsuragawa, R. Ikeda, Y. Hiai, K. Awai, S. Kusunoki, Y. Yamashita, H. Okajima, Y. Inomata, et al. Automated Hepatic Volumetry for Living Related Liver Transplantation At Multisection CT Radiology, September 1, 2006; 240(3): 743 - 748. [Abstract] [Full Text] [PDF] |
||||
![]() |
A.-J. Lemke, M. J. Brinkmann, T. Schott, S. M. Niehues, U. Settmacher, P. Neuhaus, and R. Felix Living Donor Right Liver Lobes: Preoperative CT Volumetric Measurement for Calculation of Intraoperative Weight and Volume Radiology, September 1, 2006; 240(3): 736 - 742. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| RADIOLOGY | RADIOGRAPHICS | RSNA JOURNALS ONLINE |