Radiofrequency ablation of liver tumors: quantitative assessment of tumor coverage through CT image processing
© Passera et al.; licensee BioMed Central Ltd. 2013
Received: 10 February 2011
Accepted: 4 January 2013
Published: 16 January 2013
Radiofrequency ablation (RFA) is one of the most promising non-surgical treatments for hepatic tumors. The assessment of the therapeutic efficacy of RFA is usually obtained by visual comparison of pre- and post-treatment CT images, but no numerical quantification is performed.
In this work, a novel method aiming at providing a more objective tool for the evaluation of RFA coverage is described. Image registration and segmentation techniques were applied to enable the visualization of the tumor and the corresponding post-RFA necrosis in the same framework. In addition, a set of numerical indexes describing tumor/necrosis overlap and their mutual position were computed.
After validation of segmentation step, the method was applied on a dataset composed by 10 tumors, suspected not to be completed treated. Numerical indexes showed that only two tumors were totally treated and the percentage of a residual tumor was in the range of 5.12%-35.92%.
This work represents a first attempt to obtain a quantitative tool aimed to assess the accuracy of RFA treatment. The possibility to visualize the tumor and the correspondent post-RFA necrosis in the same framework and the definition of some synthetic numerical indexes could help clinicians in ameliorating RFA treatment.
Image-guided radiofrequency ablation is a powerful locoregional technique for the treatment of the non resectable primary and metastatic hepatic malignancies and a curative alternative to surgery in circumscribded tumors (inferior to 3 cm in diameter) . Total necrotization of the target-tumor is the crucial condition of an effective treatment. However, there are many factors that could lead to an incomplete tumor ablation, such as insufficient visibility of tumor on US images, heat propagation prevented by blood vessels, and dimensions and morphology of the target tumor [2, 3]. Therefore, it would be fundamental to assess the overlap between the tumor and the RFA necrotized tissues after the treatment. In current clinical procedures, the two major indicators of tumor necrotization and treatment efficacy are parenchyma changes seen on post-treatment CT and MR images and tumor recurrence in the patient’s follow-up. However, some subjectivity and uncertainty affect medical judgment about tumor coverage. The obstacles are: i) small contrast between necrosis and residual tumor, ii) inflammation process due to coagulative effects, iii) possible blood effusions .
Objective of this paper was to describe a method for quantitative assessment of RFA tumor coverage. To this purpose, we introduced a set of numerical indexes measuring overlap between tumor and RFA induced necrosis and their reciprocal position after realignment and segmentation of pre- and post-RFA CT images. In this work, segmentation of regions of interests (ROIs) was performed through a Fuzzy-C-means approach [5, 6], while spatial correspondence between ROIs was obtained by realigning images through a non-linear B-splines-based algorithm [7, 8] able to compensate for liver deformations.
The method was tested in an experimental protocol developed at the Department of Images for Diagnosis and Therapy of the Fondazione IRCCS Istituto Nazionale dei Tumori (Milan, Italy) during 2007. 5 metastases and 5 hepatocellular carcinomas (HCCs) smaller than 20 cm3 (1-4 cm, tumor diameter range) were selected from 10 patients (5 men, 5 women; age range 40-86 years) treated with RFA as they were very complex cases not suitable for surgery. These tumors were selected among those suspected not to be completely treated. This was done in order to better evidence the potentiality of the proposed method.
All patients underwent intravenous multiphase dynamic CT before and after RFA. Pre-RFA CT scans were performed 5-10 days before RFA, while post-RFA scans were performed 15-20 days after the treatment.
Intravenous multiphase dynamic CT protocol
In this study, a Siemens SOMATOM Sensation 16 CT scanner (Erlangen, Germany) was used (gantry rotation speed of 420 ms, generator of 60 KW, pitch 1-1.5, 120 kVp, 250-300 mAs). 120–140 ml of nonionic Iopamiron 370 contrast agent were injected at a rate of 3 ml/s. Delay times after injection were 45, 80, and 180 seconds for the arterial, hepatic venous, and equilibrium phases, respectively. For HCC and hypervascular metastases arterial phase was considered, for hypovascular metastases portal or equilibrium phases were evaluated.
All CT scans were performed on inspiration. CT images were reconstructed on a 512x512 grid with a slice thickness of 5 mm and a pixel size of 0.74 mm (increment 0.6 mm, reconstruction 1-5 mm, medium smooth kernel (i30F), window/level 300/40).
Depending on device availability, operator preference and tumor location, RFA was performed using a hook-tip needle housing 4 or 7 retractable curve electrodes (RITA Medical System, Mountain View, CA) or a 19 gauge MIRAS RC electrode (INVATEC, Brescia, Italia). The electrode was connected to a 460-KHz RF generator (Model 500 L; RITA Medical System, Mountain View, CA), which supplied a maximum power output of 110 W. In both cases, RFA was performed with real-time US guidance using a 3.5 MHz convex-probe (HDI 5000, ATL Ultrasound, Bathell, WA) and a guide device incorporated into the US probe in case of percutaneous puncture or during intraoperative approach. For all patients, the duration of RF energy application was in the range of 12-20 min. For HCC nodules, RFA procedure was performed after the interruption of their arterial supply by the occlusion of the hepatic artery with a balloon (diameter 11.5, length 2 cm) at the tip of a 7.0-F catheter (Medi-tech/Boston Scientific, Watertown, Mass). The occlusion balloon in the hepatic artery was filled with a mixture of saline solution and contrast materials.
RFA was performed by experienced physicians (70 RFA treatments per year).
In this work, the tumor coverage of RFA treatment was assessed by the integration between pre- and post-RFA CT images. Image processing was divided into three steps: (i) pre-processing (ii) segmentation and (iii) registration.
To improve image quality and segmentation/registration performances the following pre-processing steps were applied:
image noise reduction. CT image quality is often degraded by artifacts resulting from excessive X-ray quantum noise. In order to reduce image granularity and then to improve ROIs extraction, the original CT images were pre-filtered by a 5 × 5 median filter .
liver volume pre-segmentation. Segmentation of tumors and necrosis on the liver edge was difficult because their typical gray levels were close to the intercostal tissues ones. In order to overcome this problem, a binary mask of the liver clearly separating liver from thorax pixels, was obtained [10, 11] (Figure 1). To segment the liver region we used the Live–Wire algorithm  implemented in the MevisLab  software. Next, subsequent processing was applied to liver pixels only.
image contrast enhancement. In order to obtain a sharp distinction among different clusters composing liver images, remapping of image dynamic range was applied. As liver intensity distribution is similar to the Gaussian distribution, it is possible to define the characteristic liver range analyzing the histogram of liver intensity into the volume pre-segmented by the Live–Wire technique. The lower and the higher gray levels that bounds the liver range were identified by two thresholds corresponding to the 2-3% of the liver histogram peak . The defined liver range was then spread on the whole image dynamic range.
tumor. Reference pixels were foremost used to define the range of tumor cluster indexes. Actually, the kind of tumors considered in this study (metastasis and small HCCs) was characterized by some geometric regularity. Tumors looked rather compact in the center, while they usually tended to vanish at the periphery. This means that, in the clusterized image, the tumor had a concentric ring structure. The mean of cluster indexes of reference pixels on all slices were used to identify the type of tumor: if it was inferior to 3.5 was a metastasis, while if superior, a HCC. Therefore, pixels were classified with increasing (for metastasis) or decreasing (for HCCs) cluster indexes from the center to the periphery of the tumor (Figure 4).
necrosis. The implemented algorithm differentiated if the necrosis was uniformly hypodense or showed hyperdense areas due to some coagulative effects. In the clusterized image, uniformly hypodense necrosis were characterized by pixels belonging to the first two clusters corresponding to the darkest gray levels. These pixels were easily detected using the neighborhood conditions on the reference pixels picked by the operator. Instead, necrosis with hyperdense areas were identifiable by the presence of reference pixels with a cluster index greater than two. In the latter case, reference pixels were used to detect all hyperdense pixels and to further classify them iteratively with lower indexes. This classification allowed to treat irregular necrosis detection in the same way as uniform necrosis.
ROI extraction was particularly challenging for HCCs located very close to arterial vessels and for metastasis and necrosis on the liver edge. In these particular cases, morphological operators were used in order to distinguish real ROIs from those artificially extended due to partial volume effects at the liver periphery or to contrast enhanced opacified vessels.
The segmentation method was quite fast, taking about ten minutes considering user interaction in the live-wire technique and pixels picking.
Pre- and post-RFA images were registered through a B-spline free-form deformation algorithm [7, 8], using the normalized mutual information (NMI) as similarity measure. Liver shifts and deformations were modeled by an affine global transformation and a local transformation based on a free-form deformation model of a regularly spaced control points grid. This registration technique had several advantages: it worked in the 3D space; it was efficient and robust, as the employed similarity measure (NMI) was not influenced by intensity changes in the processed images.
Registration step (fully automatic) took about 40 minutes. This time was evaluated on a PC (Intel Pentium III). The use of a more powerful or dedicated PC could improve this performance.
Numerical Indexes for the RFA evaluation
After having extracted and registered the pre-RFA tumor and the post-RFA necrosis, it was possible to measure their overlap and mutual position. As shown in Figure 7 there were two possible situations: a totally treated tumor or a not-totally treated tumor. To describe these situation four numerical indexes were defined:
The residual tumor size. All pixels belonging to the pre-RFA tumor and not included into necrosis area were labeled as not treated pixels. Then, the residual tumor size (given by the number of non treated pixels) was computed in volume and in percentage with respect to the pre-RFA tumor size.
(ii)The tumor free margin (T.F.M.). In the case of totally-treated tumor, it was possible to calculate the minimum T.F.M. by dilating pre-RFA tumor edge in a isotropic way until the tumor remains into the necrosis. The product between number of dilating iterations and pixel dimension gave the desired thickness. As the probability of tumor recidivism was higher for limited extra-tumor necrosis thickness , this index might have clinical relevance to indicate possible pathway for recidivism.
The inter-barycentric distance (|B T – B N |). This index computed the distance between tumor and necrosis barycenter (BT and BN , respectively) providing a measure of ROIs centering.
The orientation index (O.I.). For each slice, an ellipse was fitted to each ROI by making equal second order central moments of the ellipse to those of ROIs. Then, the angle (in degrees) between the major axis of the two ellipses was measured. This angle gave indication of the reciprocal orientation between tumor and necrosis ROIs and of potential favorite orientations for necrosis growth.
The inter-barycentric distance and the O.I. provided complementary elements for a more objective RFA evaluation. In fact, a not–totally treated tumor also resulted in a misalignment with the post-ablation necrosis. In particular, an inter-barycentric distance |BT − BN | with the same magnitude order of the target tumor was the natural consequence of an inaccurate tumor centering by the necrosis. Besides, an O.I. in the critical range 60° − 120° indicated an objective difficulty of the necrosis in following the tumor morphology and orientation.
As segmentation was crucial in the proposed method, a detailed validation of this step was performed, following a classical approach reported in literature [16, 17]. An experienced radiologist manually traced the contour of 10 tumors and 10 necrosis in our datasets. Manual contours were considered as Ground Truth (GT) and compared with the output ROIs obtained with the semi-automatic algorithm.
The numerical indexes used in the evaluation of segmentation performance
Segmentation validation indexes [%]
The proposed method was applied to 8 out of 10 available cases as for 2 patients it was not possible to obtain a satisfactory realignment due to contingent image features: in the first case, serous fluid effusion in the abdominal cavity (ascites) caused the presence of a hypodense strip in the post–RFA image only; in the second one, the patient underwent not only RFA-treatment but also hepatic resection. In both cases the registration algorithm was not able to correct differences between pre- and post-RFA images. Such cases were thus used for the segmentation validation only. In the remaining 8 cases (4 HCCs and 4 metastasis), the complete RFA assessment procedure was performed.
The implemented segmentation algorithm took these considerations into account. In fact, for tumors the whole cluster range was assessed starting from the operator reference pixels, while for necrosis only the first two clusters (corresponding to the darkest and the most compact ROI area) were extracted, limiting risk of error by excess in necrosis segmentation procedure.
In order to assess the opposite risk of error by defect SPEC and P- indexes were calculated. Figure 8(right) shows that also the SPEC and P- mean values are very close to the optimum point on the top-right angle (SPEC > 96% and P− > 92%). This demonstrates a good balance between errors by excess and by defect in the segmentation process.
Assessment of RFA tumor coverage
The pre-RFA and post-RFA image registration and segmentation operations enabled to visualize in the same framework the overlap between the pre-RFA tumor and the corresponding necrosis. In this way, information was automatically integrated and it was not required a physician’s effort in finding anatomical markers for visually locating the tumor after the RFA-treatment.
The quantitative analysis of the tumor and post-ablation necrosis overlap for 8 cases
Tumor slices (#)
Tumor volume (#)
Residual tumor* (cm3) [%]
T.F.M. (mm) [min max]
|BT-BN| (mm) [mean ± std]
O.I. (degree) [mean ± std]
Case 1 META
0 - 1.73
11.59 ± 1.52
57.46 ± 22.22
Case 2 META
0.74 - 2.23
4.96 ± 3.54
33.59 ± 20.82
Case 3 META
0 - 2.16
7.60 ± 0.45
74.72 ± 5.27
Case 4 META
0 - 6.69
16.98 ± 3.38
56.60 ± 18.54
Case 5 HCC
1.93 - 3.87
6.61 ± 1.52
63.94 ± 4.50
Case 6 HCC
4.07 ± 0.17
25.31 ± 11.03
Case 7 HCC
0 - 3.64
7.57 ± 3.30
19.04 ± 12.75
Case 8 HCC
2.75 ± 2.41
31.33 ± 20.08
Discussion and conclusion
In this paper, we presented a novel method to assess the hepatic tumor coverage of RFA. In the following, some comments about the method, the results and the study limitations were reported.
A previous study  dealt with the problem of finding an objective RFA assessment method based on CT images segmentation and registration. There were three main differences between that study and the method proposed in this work. First of all, the method proposed in  was designed for pre- and post-RFA images fusion in order to enable easier understanding of the relationship between the tumor and ablation zone and to help to judge whether an ablative margin was ensured or not, but did not provide any numeric index for quantitative analysis and assessment of the RFA treatment. Secondly, in  ROI segmentation was performed manually, while in this work was based on a semi-automatic segmentation algorithm, requiring user interaction only for the selection of a few reference pixels needed to outline liver edges by the Live–Wire technique. As compared with manual segmentation, the algorithm used in this work was more efficient, as it was based on the automatic pixels gray level analysis performed through a Fuzzy-C-Means approach, it was reproducible thus avoiding uncertainty and subjectivity of manual edges outlining.
Finally, in our method, the pre- and post-RFA CT image realignment was based on the B-splines free form deformation algorithm [7, 8] that allowed a non-linear hepatic volumes registration, instead of a simpler rigid registration algorithm with additional manual adjustments, which compensated only for translation and rotation shifts that occurred between pre- and post-RFA acquisitions. As the liver is a soft tissue undergoing non-linear deformations, mainly caused by respiration, heart pulsations and adjacent organ movements, a non-linear registration is required for correct volumes realignment [20, 21].
There were peculiar properties which made the proposed algorithm particularly appealing for the assessment of RFA treatment in liver tumors: the short interaction time required by users to outline liver edges by the Live–Wire technique and to acquire the reference pixels on the tumor and necrosis ROIs (a few minutes for each image volume), and the accuracy, proved by good results of segmentation validation.
Even if in this study the tumor coverage of RFA was assessed on a small dataset of suspected non-completely treated tumors, it is worth noting that most of the tumors were not totally-treated (6 out of 8) and the percentage of a residual tumor was in the range of 5.12%-35.92%. In addition, T.F.M. was always well far from 1 cm-thick as recommended in current RFA guidelines .
There are several reasons that may explain these results. Before the treatment, the surgeon might be often aware that a large T.F.M. was not possible. In particular, the target tumor was frequently close to vital anatomical structures which had to be prevented from heat injury. Besides, there were physical limitations to heat spread associated to the liver blood flow (heat loss due to convection), that were enhanced also by O.I. and inter-barycentric distance indexes. Finally, the type of RFA treatment might affect the heat spread, such as in the case of RFA-treatment executed after occlusion of tumor blood supply. It was demonstrated that a larger necrosis area can be created when RFA treatment is performed in HCC nodules after their arterial supply occlusion [22–24]. In fact, there is a different temperature distribution within and around the HCC nodule and this phenomenon seems related to the difference in vascularization between HCC and the surrounding cirrhotic hepatic tissue . The latter has a dual blood supply and is nourished mainly by the portal vein, which provides about two-thirds of the blood flow. HCC, however, is nourished mainly by the hepatic artery, with the portal vein providing a minor blood supply and the main venous drainage. Therefore acute occlusion of arterial flow is soon followed by a decrease in pressure within the HCC nodule, which continues to be perfused by means of reversed portal flow and, in some cases, by small collateral arteries.Thus, although the blood flow supplying the HCC nodule is substantially impaired, changing from high to sluggish flow, the blood flow supplying the surrounding hepatic tissue is only marginally modified. This results in an almost complete lack of heat loss due to convection within the HCC nodule, whereas intact and perhaps even increased portal blood flow in the surrounding tissue acts as an efficient heat sink that prevents heat diffusion outside the HCC nodule . Therefore, the resultant necrosis reproduces the shape of HCC nodules and spare surrounding non-tumor tissues, producing a really tiny T.F.M.
The present study had some limitations and some improvements could be performed. First of all, a larger data-set for the quantitative analysis of RFA-treatment could allow an adequately statistical analysis of the results. In addition, the evaluation of segmentation was performed on the same image dataset, which had also used to develop the segmentation algorithm. It is likely that the performance will be slightly lower when larger dataset of unseen data are considered. Secondly, the use of isotropic and small voxels (slice thickness < 5 mm) could improve the segmentation algorithm performance. In this case, in fact, it could be possible to perform a 3D segmentation and to include in the clustering approach additional spatial information (by using, for example, the hidden Markov random fields theory) [25, 26]. Concerning registration method, the use of a synthetic pattern replacing tumor and necrosis ROIs, that prevents the registration algorithm from modeling fictitious deformation, had the drawback not to taking potential local deformations due to the heat coagulative effects into account. In a recent study , it was found that RF and microwave ablation both cause significant contraction of normal bovine liver and lung tissue ex vivo. Ablation-induced contraction appears to be tissue type and ablation modality specific. This phenomenon should be studied in detail and modeled into the registration process.
Nevertheless, this work represents a first attempt to obtain a quantitative tool aimed to assess the accuracy of RFA treatment, and its major contribution is the definition of several numerical indexes that could be helpful to quantitatively describe RFA treatment, pointing out potential limitations.
Finally, from a clinical point of view it would have been also very interesting to follow up the cases with incomplete tumor necrosis. Unfortunately we have a limited dataset consisting of 10 tumors only (5 HCC and 5 metastasis) almost heterogeneous in nature. For this purpose we are planning a new study including a larger dataset grouped for type of tumors. In this new study, we think to investigate other valuable indexes such as the largest axial diameter of tumor and necrosis and difference of the HU-units in order to provide possible surrogates for the modified RECIST criteria.
- Buscarini E, Savoia A, Brambilla G, Menozzi F, Reduzzi L, Strobel D, Hansler J, Buscarini L, Gaiti L, Zambelli A: Radiofrequency thermal ablation of liver tumors. Eur Radiol. 2005, 15: 884-894. 10.1007/s00330-005-2652-x.View ArticlePubMedGoogle Scholar
- Padma S, Martinie JB, Iannitti DA: Liver tumor ablation: percutaneous and open approaches. J Surg Oncol. 2009, 100 (8): 619-634. 10.1002/jso.21364.View ArticlePubMedGoogle Scholar
- Lencioni R, Crocetti L: Radiofrequency ablation of liver cancer. Tech Vasc Interv Radiol. 2007, 10 (1): 38-46. 10.1053/j.tvir.2007.08.006.View ArticlePubMedGoogle Scholar
- Decadt B, Siriwardena AK: Radiofrequency ablation of liver tumours: systematic review. Lancet Oncol. 2004, 5 (9): 550-560. 10.1016/S1470-2045(04)01567-0.View ArticlePubMedGoogle Scholar
- Xu R, Wunsch DC: 2nd: Clustering algorithms in biomedical research: a review. IEEE Rev Biomed Eng. 2010, 3: 120-154.View ArticlePubMedGoogle Scholar
- Yang Z, Chung FL, Shitong W: Robust fuzzy clustering-based image segmentation. Appl Soft Comput. 2009, 9 (1): 80-84. 10.1016/j.asoc.2008.03.009.View ArticleGoogle Scholar
- Rueckert D, Hawkes D, Sonoda L, Hill D, Leach M, Hawkes D: Non-rigid registration using free-form deformations: application to breast MR images. IEEE Trans Med Imag. 1999, 18: 712-721. 10.1109/42.796284.View ArticleGoogle Scholar
- Schnabel J, Rueckert D, Quist M, Blackall J, Castellano-Smith A, Hartkens T, Penney G, Hall W, Liu H, Truwit C, Gerritsen F, Hill D, Hawkes D: A generic framework for nonrigid registration based on nonuniform multilevel free-form deformations. MICCAI. 2001, 2208: 573-581.Google Scholar
- Ciecholewski M, Ogiela M: Automatic segmentation of neoplastic hepatic disease symptoms in CT Images. Modeling Decisions for Artificial Intelligence. 2007, Springer-Verlag, Berlin Heidelberg, 414-421.View ArticleGoogle Scholar
- Barrett WA, Mortensen EN: Interactive live-wire boundary extraction. Med Image Anal. 1997, 1 (4): 331-341. 10.1016/S1361-8415(97)85005-0.View ArticlePubMedGoogle Scholar
- Campadelli P, Casiraghi E, Esposito A: Liver segmentation from computed tomography scans: a survey and a new algorithm. Artif Intell Med. 2009, 45 (2-3): 185-196. 10.1016/j.artmed.2008.07.020.View ArticlePubMedGoogle Scholar
- MeVis Research GmbH Germany. http://www.mevislab.de 2008
- Lim SJ, Jeong YY, Ho YS: Automatic liver segmentation for volume measurement in CT Images. J Vis Commun Image R. 2006, 17 (4): 860-875. 10.1016/j.jvcir.2005.07.001.View ArticleGoogle Scholar
- Yim P, Vora A, Raghavan D, Prasad R, McAullife M, Ohman-Strickland P, Nosher JL: Volumetric analysis of liver metastases in computed tomography with the fuzzy C-means algorithm. J Comput Assis Tomogr. 2006, 30: 212-220. 10.1097/00004728-200603000-00008.View ArticleGoogle Scholar
- Hong J, Kaneko T, Sekiguchi R, Park K: Automatic liver tumor detection from CT. IEICE T INF SYST. 2001, E84-D (6): 741-748.Google Scholar
- Udupa JK, LeBlanc VR, Zhuge Y, Imielinska CZ, Schmidt H, Currie LM, Hirsch BE, Woodburn J: A framework for evaluating image segmentation algorithms. Comput Med Imaging Graph. 2006, 30 (2): 75-87. 10.1016/j.compmedimag.2005.12.001.View ArticlePubMedGoogle Scholar
- Pham D, Xu C, Prince J: Current methods in medical image segmentation. Annu Rev Biomed Eng. 2000, 2: 315-337. 10.1146/annurev.bioeng.2.1.315.View ArticlePubMedGoogle Scholar
- Alonso F, Algorri E, Flores-Mangas F: Composite index for the quantitative evaluation of image segmentation results. Proc. Ann. Int. Conf. of the IEEE Engineering in Medicine and Biology Society: 1-4 September 2004. 2004, San Francisco, 1794-1797.View ArticleGoogle Scholar
- Fujioka C, Horiguchi J, Ishifuro M, Kakizawa H, Kiguchi M, Matsuura N, Hieda M, Tachikake T, Alam F, Furukawa T: A feasibility study: evaluation of radiofrequency ablation therapy to hepatocellular carcinoma using image registration of preoperative and postoperative CT. Acad Radiol. 2006, 13: 986-994. 10.1016/j.acra.2006.05.011.View ArticlePubMedGoogle Scholar
- Crum WR, Hartkens T, Hill DL: Non-rigid image registration: theory and practice. Br J Radiol. 2004, 77 (Spec 2): S140-S153.View ArticlePubMedGoogle Scholar
- Mainardi L, Passera K, Lucesoli A, Potepan P, Setti E, Musumeci R: A method for dynamic subtraction MR imaging of the liver. BMC Med Imaging. 2006, 6: 5-10.1186/1471-2342-6-5.View ArticlePubMedPubMed CentralGoogle Scholar
- Rossi S, Garbagnati F, De Francesco I, Accocella F, Leonardi L, Quaretti P, Zangrandi A, Paties C, Lencioni R: Relationship between the shape and size of radiofrequency induced thermal lesions and hepatic vascularization. Tumori. 1999, 85 (2): 128-132.PubMedGoogle Scholar
- Rossi S, Garbagnati F, Lencioni R, Allgaier H, Alfonso A, Fornari F, Quaretti P, Tolla GD, Ambrosi C, Mazzaferro V, Blum H, Bartolozzi C: Percutaneous radiofrequency thermal ablation of nonresectable hepatocellular carcinoma after occlusion of tumor blood supply. Radiology. 2000, 217: 119-126.View ArticlePubMedGoogle Scholar
- Rossi S, Di Stasi M, Buscarini E, Quaretti P, Garbagnati F, Squassante L, Paties CT, Silverman DE, Buscarini L: Percutaneous RF interstitial thermal ablation in the treatment of hepatic cancer. AJR Am J Roentgenol. 1996, 167 (3): 759-768. 10.2214/ajr.167.3.8751696.View ArticlePubMedGoogle Scholar
- Zhang Y, Brady M, Smith S: Segmentation of brain MR images through a hidden markov random field model and the expectation-maximization algorithm. IEEE Trans Med Imag. 2001, 20: 45-57.View ArticleGoogle Scholar
- Azmi R, Norozi N: A new markov random field segmentation method for breast lesion segmentation in MR images. J Med Signals Sens. 2011, 1 (3): 156-164.PubMedPubMed CentralGoogle Scholar
- Brace CL, Diaz TA, Hinshaw JL, Lee FT: Tissue contraction caused by radiofrequency and microwave ablation: a laboratory study in liver and lung. J Vasc Interv Radiol. 2010, 21 (8): 1280-1286. 10.1016/j.jvir.2010.02.038.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2342/13/3/prepub