Prospective volumetric assessment of the liver on a personal computer by nonradiologists prior to partial hepatectomy.

BACKGROUND: A small remnant liver volume is an important risk factor for posthepatectomy liver failure. ImageJ and OsiriX® are both free, open-source image processing software packages. The aim of the present study was to compare ImageJ and OsiriX® in performing prospective computed tomography (CT) volumetric analysis of the liver on a personal computer (PC) in patients undergoing major liver resection.
METHODS: Patients scheduled for a right hemihepatectomy were eligible for inclusion. Two surgeons and one surgical trainee measured volumes of total liver, tumor, and future resection specimen prospectively with ImageJ and OsiriX®. A radiologist also measured these volumes with CT scanner-linked Aquarius iNtuition® software. Resection volumes were compared with the actual weights of the liver specimens removed during surgery, and differences between the measured liver volumes were analyzed.
RESULTS: A total of 15 patients (8 men, 7 women) with a median age of 63 years (48-79 years) were included. There was a significant correlation between the measured weights of resection specimens and the volumes calculated prospectively with ImageJ and OsiriX® (r=0.89; r=0.83, respectively). There was also a significant correlation between the volumes measured with radiological software iNtuition® and the volumes measured with ImageJ and OsiriX® (r=0.93; r=0.95, respectively).
CONCLUSIONS: There were no major differences in total liver volumes, resection volumes, or tumour volumes for these three software packages. Prospective hepatic CT volumetry with ImageJ or OsiriX® is reliable and can be accurately used on a PC by nonradiologists. ImageJ and OsiriX® yield results comparable to the radiological software iNtuition®.

Introduction

Hepatectomies are performed increasingly more often, mainly because indications for liver resection are continuously being extended. Whereas in the past a number of exclusion criteria were applied to the selection of patients for hepatectomy, they have now been largely abandoned, and postoperative remnant liver volume and function have become the main determinants of resectability [1-3].

Postoperative liver dysfunction may occur when the extent of tumor involvement requires major liver resection (three or more segments), leaving a small postoperative remnant liver [4-6]. A reliable volumetric assessment of the part of the liver to be resected and the future residual liver volume should be a critical part of preoperative evaluation, particularly for extended resections, to prevent postoperative liver failure after resection. The safety of liver resection may increase if an estimate of minimal remnant liver volume is obtained via computed tomography (CT) volumetry [4, 6–9]. This is especially important if the underlying liver parenchyma is abnormal or has impaired functional capacity. In this context, repeated liver resections, increasingly extensive resections, and resections after aggressive chemotherapy are likely to lead to more patients with small (functional) remnant livers and hence more risk of postoperative liver failure.

As we [4, 10–12] and others [7, 13, 14] have shown before, liver volumes can be calculated accurately from CT or magnetic resonance imaging (MRI) scans. The liberal use of existing professional image-processing software is often limited by costs, access, flexibility, lack of radiology staff, and specific hardware requirements. In addition, the intended operation and the potential surgical scenarios should be known to the investigator to predict the remnant liver volume accurately, which requires the expertise of a liver surgeon. The Digital Imaging and Communications in Medicine (DICOM) standard has enabled volumetry on a stand-alone or network-attached personal computer (PC) remote from the radiology hardware (CT or MRI scanner).

In 2007, our group conducted a retrospective study and proposed a method to perform liver volumetry with ImageJ [11], which is a free open-source image analysis software package developed at the National Institutes of Health (NIH) [8]. In 2010, our group assessed the accuracy of OsiriX® for CT volumetry of the liver [12]. OsiriX® is a more user-friendly, freely downloadable open source image analysis software package for the Apple Mac OS. The accuracy of CT volumetry of the liver resection specimen prospectively prior to surgery has not yet been assessed for either of these software packages. Therefore, the objective of the present study was to compare ImageJ and OsiriX® when performing prospective CT volumetric analysis of the liver on a PC in patients undergoing major liver resection. The study also planned to compare the prospective measurements of ImageJ and OsiriX® to those of CT volumetry performed by a radiologist using classic radiology software (iNtuition®) linked to a CT scanner system.

Methods

Patients

Patients with primary or secondary liver tumors in otherwise normal livers (mostly colorectal cancer liver metastases) planned for a right hemihepatectomy at Maastricht University Medical Center between January 2008 and April 2009 were eligible for inclusion in this study. All patients had a contrast enhanced CT scan in their routine preoperative assessment either in our hospital or in one of the surrounding university-affiliated district general teaching hospitals. Patients were admitted to the hospital 1 day preoperatively, and routine blood tests were performed.

Surgical procedure

Liver resection was performed as detailed elsewhere [15]. In each patient, laparotomy was performed by bilateral subcostal incision, followed by intraoperative ultrasonographic (US) assessment of the liver. Once resectability had been confirmed, appropriate mobilization of the liver was performed as a prelude to hepatic parenchymal transection, which was undertaken using a Cavitron Ultrasonic Surgical Aspirator (Force GSU System; Valleylab, Boulder, CO, USA). Argon beam coagulation (Erbe, Tübingen, Germany), clips, and sutures were used for hemostasis. Immediately after liver resection, the weights of the resected specimens were recorded in the operating theater with a Mettler Toledo scale.

CT volumetry

For volumetric analysis, four-phase CT scans were used that were provided on CD-ROM on one of two viewers: eFilm Lite (eFilm Medical, Toronto, Canada) or SIENET MagicView 300 VA42D (Siemens, Erlangen, Germany).

Calculation of the liver volume using ImageJ

ImageJ (version 1.33) was downloaded from http://rsb.info.nih.gov/ij/download.html (accession date 01-01-08). The portovenous phase of all contrast enhanced CT scans was used to facilitate optimal identification of liver segments and the anatomic resection plane. Details on how volumetry was performed with ImageJ have been described previously [11]. Briefly, the relevant source images were individually copied (Ctrl + V) from the CT viewer on the electronic patient care system and opened with ImageJ by dragging them to the ImageJ main window in the appropriate order. In cases of tertiary referred patients, relevant CT slices were searched in the original CT viewer on the CD-ROM. CT-slices on a CD-ROM were usually saved as compressed DICOM files, which cannot be opened or viewed via Windows Explorer. Therefore, the “Tudor DICOM Viewer Plug-in” was downloaded (free) from http://imagejdocu.tudor.lu/ and installed in the ImageJ menu. This plug-in allowed us to preview the relevant CT files and open them directly in ImageJ.

The total liver, tumor, and resected specimen were manually outlined. The respective region of interest (ROI) of each slice was added to the ROI manager. The area of each ROI was calculated and then multiplied by the slice thickness. Finally, all slice volumes were added up to calculate the volume of each three-dimensional structure (i.e., total liver, resected specimen, tumor).

Prospective volumetric analysis: ImageJ

Volumetric analysis using ImageJ was performed prospectively by two nonradiologists (S.A.W.G.D., J.H.M.B.S.) during the preoperative evaluation of the patient and was blinded to the operating surgeon. For definition of the anatomic segments, Couinaud’s classification was used. In all right hemihepatectomies, transection followed Cantlie’s line from the top of the gallbladder, paralleling the middle hepatic vein. straight to the suprahepatic inferior caval vein. In all patients, the middle hepatic vein remained in situ with the liver remnant. Consequently, during the venous phase, the middle hepatic vein on CT was used to perform the volumetric measurement.

To assess the accuracy of the volumetric measurements of the liver with ImageJ, the volumes of the resected specimens were compared with the actual weights of the specimens. (The weights of the specimens had been recorded in the operating theater immediately after resection.) The actual weights of the resection specimens remained blinded to the investigators conducting CT volumetry until all patients had been included and analyzed prospectively. The functional remnant liver volume (FRLV) percentage was calculated with the following formula. FRLV [%] = 100% * [total liver volume – resected volume] / [total liver volume – tumour liver volume]

Prospective volumetric analysis: OsiriX®

The 32-bit OsiriX® version 3.3 was downloaded from: http://www.osirix-viewer.com. A 2.8-GHz Intel Core 2 Duo 24′′ iMac (Apple, Cupertino, CA, USA) was used. Details on how volumetry was performed with OsiriX® have been described previously [12]. Volumes calculated with ImageJ were compared to the volumes obtained by CT volumetry using OsiriX® by one of the liver surgeons (R.M.vanD.).

Comparison of ImageJ and OsiriX® with iNtuition®

Volumes calculated with ImageJ and OsiriX® were compared to the volumes obtained by CT volumetry performed by a specialized liver radiologist (R.S.A.vanS.) using radiology image analysis software from TeraRecon linked to an Aquarius iNtuition® workstation. Volumetric analysis performed with ImageJ and OsiriX® was blinded to the radiologist. The volume of the liver resection specimen measured by the radiologist was used as a second standard of reference instead of the actual weight of the resected specimen.

Reliability: ImageJ

For analysis of interuser variability of ImageJ, volume measurements of total liver volume, functional liver volume, and tumor volume performed before surgery by two investigators (S.A.W.G.D., J.H.M.B.S.) were compared. The first investigator had 3 years’ experience and the second investigator had no experience with the applied software and received instructions and training before starting with the measurements. Bland-Altman plots were used (see below, under Statistical Analysis) to analyze the reliability of the ImageJ method.

Time consumption: ImageJ

The time consumed to perform one volumetric assessment of a liver with ImageJ including the total liver, the tumor, and the resected specimen was recorded for each of the 15 livers. To assess time consumption, workload, and influence of slice thickness, the time needed to outline these three ROIs (total liver, tumor, resected specimen) was recorded per slice in five livers.

As for any volumetric measurement technique based on CT imaging data sets, the accuracy of the volume measurements is dependent on the slice thickness used for the volume calculation. The slice thickness used for volume measurement is also crucial for the time needed to perform the analysis because of manual outlining. Therefore, the influence of slice thickness on CT-based liver volume measurements using ImageJ was also investigated. Total liver volume (TLV) was analyzed based on different slice thicknesses in these patients undergoing a right hemihepatectomy. The relative change in TLV between volumetry performed with the standard of reference (0.5 cm) and volumetric measurements using thicker slices (1.0, 1.5, and 2.0 cm) was calculated with the following formulas. (TLV1.0 − TLV0.5)/TLV0.5 × 100%; (TLV1.5 − TLV0.5)/TLV0.5 × 100%; (TLV2.0 − TLV0.5)/TLV0.5 × 100%; The relative gain of time between volumetric measurements using different slice thicknesses was recorded and calculated in percentages as well.

Statistical analysis

Pearson’s coefficient was used to quantify the association between the volume of the resected specimen and the actual weight of the resection specimen. A value of P < 0.05 was considered to indicate statistical significance. The resulting regression line was described as a linear equation, and the correlation coefficient (r) was calculated. Changes of calculated volumes were tested against zero with the use of a Wilcoxon signed rank test with a theoretical mean of zero. Data are presented as the mean ± standard error of the mean (SEM).

The mean and standard deviation (SD) were calculated from the two volumetric results for each user using ImageJ. Next, the limits of agreement as proposed by Bland and Altman [16] were calculated. This method is based on the mean and SD of the difference between two ratings of the same patient. The mean difference between two users indicates systematic error. The SD of the difference between two users indicates random error. The 95% limits of agreement as in Bland and Altman’s plots result in both quantification of random error and systematic error (1.96 * SD) of the measurement method. Statistics were performed using Prism 5.0 for Windows (GraphPad Software, San Diego, CA, USA).

Results

Fifteen patients scheduled for right hemihepatectomy were included in the study. The median age of the patients was 63 years (48–79 years). None of the patients had preoperative evidence of underlying liver disease, and routine liver function tests were uniformly normal (Table 1).

Table 1

Patients’ characteristics

Characteristic	Median (range)
Age (years)	63 (48–79)
Sex (M/F)	(8/7)
Height (cm)	170 (155–191)
Weight (kg)	76 (55–127)
BMI (kg/m2)	26 (19–43)
AST (IU/l)	22 (7–52)
ALT (IU/l)	30 (8–61)
LDH (IU/l)	372 (312–469)
γGGT (IU/l)	37 (29–258)
ALP (IU/l)	98 (57–168)
Bilirubin (μM)	12 (7–15)

BMI body mass index, AST aspartate aminotransferase, ALT alanine aminotransferase, LDH lactate dehydrogenase, γGGT γ-glutamyltranspeptidase, ALP alkaline phosphatase

Relation between prospective CT-measured liver volume and resected liver weight

The mean weight of the resection specimens measured in the operating theater was 788.8 ± 53.7 g. The mean resected liver volumes calculated with ImageJ, OsiriX®, and iNtuition® were 987.7 ± 64.0, 944.7 ± 64.1, and 905.0 ± 59.9 ml, respectively. The mean ratios between the measured weight and volume measured by ImageJ, OsiriX®, and iNtuition® were 0.78 ± 0.03, 0.81 ± 0.04, and 0.83 ± 0.04, respectively. The mean remnant liver volume percentages measured by ImageJ, OsiriX®, and iNtuition® were 41.4 ± 1.6%, 36.2 ± 1.5%, and 37.8 ± 1.6%, respectively (Table 2).

Table 2

Results for 15 CT volumetry and intraoperative measurements

Measurement	ImageJ user 1	ImageJ user 2	OsiriX®	iNtuition®
Total liver volume (ml)	1736.0 ± 104.0	1785.0 ± 114.9	1536.0 ± 98.6	1540.0 ± 99.9
Resected liver volume (ml)	987.7 ± 64.0	1050 ± 78.6	944.7 ± 64.1	905.0 ± 59.9
Tumor volume (ml)	14.0 ± 4.3	18.9 ± 4.4	12.0 ± 3.7	12.4 ± 4.0
Remnant liver volume (%)	41.4 ± 1.6	39.0 ± 1.4	36.2 ± 1.5	37.8 ± 1.6
Resected specimen weight/volume ratio (g/ml)	0.78 ± 0.03	0.77 ± 0.03	0.81 ± 0.04	0.83 ± 0.04
Resected specimen weight (g)	788.8 ± 53.7

CT computed tomography

In one patient two surgical procedures were possible, either a right posterior sectionectomy or a right hemihepatectomy. Both volumes were calculated prospectively, and finally the volume for right posterior sectionectomy was included in the analysis. There was a significant correlation between the resected liver volumes calculated with ImageJ, OsiriX®, and iNtuition® prospectively and the actual measured weights of the resection specimens: r = 0.89 and P < 0.001; r = 0.83 and P < 0.001; r = 0.85 and P < 0.001; respectively (Fig. 1).

Fig. 1

Correlation between volume measured prospectively with ImageJ, OsiriX®, and iNtuition® and the resected specimen’s weight measured in the operating theater

Interobserver variability: ImageJ

Liver volumetry measurements using ImageJ were compared using Bland-Altman plots. The mean difference of TLV measurements between both users was 49.4 ± 86.1 ml. The 95% limits of agreement were −218.2 to 119.4 ml (Fig. 2a). The mean difference of the liver resection specimen measurements between both users was 53.0 ml ± 89.8 ml. The 95% limits of agreement were −228.9 to 122.9 ml (Fig. 2b). The mean difference of the %FRLV measurements between both users was 2.3 ± 4.6%. The 95% limits of agreement were −6.7% to 11.4% (Fig. 2c).

Fig. 2

a Bland-Altman analysis for interobserver variability of outlining the total liver volume with ImageJ. Horizontal solid line indicates the mean average between the two users. Horizontal dashed lines indicate 95% limits of agreement (mean ± 1.96 SD). b Bland-Altman analysis for interobserver variability of outlining the resectional liver volume with ImageJ. c Bland-Altman analysis for interobserver variability of the %functional remnant liver volume (%FRLV) with ImageJ

Time to perform measurement and effect of slice thickness

The time required for a single volumetric assessment of the liver with ImageJ (including total liver, tumor, resected specimen) was, on average, 20 minutes (range 15–30 minutes). There was no statistical difference between the TLVs based on 0.5 vs. 1.0 cm, 0.5 vs. 1.5 cm, and 0.5 vs. 2.0 cm. The change in volume (calculated in percentages) was not significantly different for volumetric measurements with different slice thicknesses. For a slice thickness of 1.0 or 1.5 cm, 80% of the measurements were within the 5% limit of acceptance for the TLV change. When using the 2.0 cm slice thickness, only 60% of the measurements were within this limit. The mean gain of time (calculated in percentages) by increasing slice thickness was 50% (0.5 vs. 1.0 cm), 64% (0.5 vs. 1.5 cm), and 72% (0.5 vs. 2.0 cm).

ImageJ versus OsiriX® and iNtuition®

There was a highly significant correlation between the volumes calculated with ImageJ and the volumetric analysis using OsiriX® (r = 0.95, P < 0.001). There was also a strong correlation between the volumes calculated with ImageJ and the volumetric analysis performed by a radiologist using radiology image analysis software linked to radiology hardware (r = 0.93, P < 0.001).

Discussion

The present study shows that prospective CT volumetric analysis of the liver in patients undergoing right hemihepatectomy with ImageJ and OsiriX® on the surgeon’s PC is valid and reliable. In this study, we found a strong significant correlation between resection weight and resection volume measured prospectively with ImageJ and OsiriX®. Practical advantages of ImageJ and OsiriX® were described previously by our group [11, 12]. In our center, the radiologists use the commercial image analysis software package iNtuition® for liver CT volumetry. There was a significant correlation between the resection weight and the resection volume measured prospectively with iNtuition®. There were no significant differences in the mean volumes between the three software packages. There was a strong correlation between the volumes measured with radiology software iNtuition® and the volumes measured with ImageJ and OsiriX®. This indicates that ImageJ and OsiriX® even if operated by nonradiologists, are at least as appropriate for CT volumetry of the liver as professional radiology software. ImageJ and OsiriX® bring preoperative liver CT volumetry within the surgeon’s reach, independent of the support of a radiologist. Hepatopancreatobiliary (HPB) surgeons should routinely perform CT volumetry during the preoperative assessment of patients undergoing a major liver resection. This is especially true for patients with compromised liver function due to neoadjuvant chemotherapy, cirrhosis, or steatosis or for patients scheduled for an extended liver resection.

As we [11] and others [17] concluded before, the time required to perform CT volumetry of the liver can be reduced by increasing slice thickness, thereby diminishing the number of CT slices that have to be outlined. Data from the current study show that to make a rough estimate of the FRLV (patients with a percentage FRLV of ≥30%) a slice thicknesses of 1.0 cm provides an optimal balance between accuracy and time efficiency. However, in patients with a marginally functional liver remnant volume (20–30%) greater accuracy is needed. In these patients, it is better to use a smaller slice thickness to avoid over- or underestimating the volume. Time required to outline all ROIs is approximately 20 minutes with both ImageJ and OsiriX®. However, OsiriX® has a significant advantage because it has an integrated picture archiving and communication system (PACS system) in which the data (CT scans) of all examined patients are automatically stored. Additionally, OsiriX® is more user-friendly, and ROIs can be selected semiautomatically, which saves time and lessens the workload [12].

Interobserver variability for ImageJ between the two users for volumetric assessment of the total and resected liver volumes was low. Bland-Altman plots revealed that there was one outperformer in the measurements of the total liver volume and the resection specimen volume. This resulted in relatively wide confidence intervals. However, it is mainly the %FRLV that is crucial. It has been shown elsewhere that to perform liver resection safely %FRLV has to be at least 26.6% in patients with normal liver parenchyma [4]. In the present study, the mean difference of the %FRLV measurements between both ImageJ users of 2.3 ± 4.6% was clinically acceptable because ImageJ users 1 and 2 predicted %FRLVs of 41.4 ± 1.6% and 39.0 ± 1.4%, respectively. In other words, the resection volume differences between the two ImageJ users should have had no clinical consequences. In patients in whom the preoperative liver volume assessment shows that the reported FRLV% approaches 26.6%, extra care is recommended.

A possible limitation of estimating the future remnant liver volume with CT volumetry is the systematic overestimation of the liver volumes compared to the resected specimens’ weight. A physical density <1.0 g/ml and systematic overestimation of the liver volume in the present study was in accordance with those of previous studies [7, 10–12]. Inaccurate estimation of the resected specimen volume by CT volumetry can occur because the liver on a CT scan is in a perfused state and the postoperative resection specimen is in an unperfused state. This applies to every radiologic application, and therefore some studies have suggested the implementation of conversion factors to improve volumetric accuracy [7, 18, 19]. These factors are mostly based on the intraoperative resection volume determined with the use of Archimedes principle. Unfortunately, this approach does not determine the influence of the blood content of the resection specimen. It might be that these conversion factors can lead to a more precise volumetric analysis of the liver, but their exact role remains to be determined. Patients with steatosis have more fat in the liver, and the relative densitie of fat is 0.9 g/ml. Hence, the volume is larger and the weight less than that of a normal liver. Ultimately, this changes the volume/weight standard and can also contribute to systemic overestimation.

In all right hemihepatectomies in the present study, virtual transection followed Cantlie’s line, and thus the middle hepatic vein (MHV) was used as a guide for prospective volume measurements. During surgery, the demarcation line after selective inflow occlusion is used for transection and may leave liver parenchyma outside (i.e., to the right of) the MHV in some patients. As a consequence, preoperative estimation of the future resected specimen volume may be overestimated in some individuals. In line with this, Karlo et al. showed in a retrospective study that, despite a highly accurate definition of the transection line on preoperative CT or MRI scans using anatomic landmarks on postoperative MRI scans, the resection specimen volume was still overestimated [7]. Systematic overestimation, however, does not affect the functional remnant liver volume as a percentage of the total functional liver volume, which is one of the most important indicators of the risk of postoperative liver failure.

Conclusions

User-friendly, easily accessible, free, downloadable software packages such as ImageJ and OsiriX® can be used reliably for prospective CT volumetry of resected liver specimens. The performance of ImageJ and OsiriX® is comparable to that of radiology software iNtuition®. OsiriX® is the most efficient and practical software package, but it can be used only on a Macintosh PC. For a Windows-based PC, ImageJ is a good alternative.