- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 964918, 9 pages

http://dx.doi.org/10.1155/2013/964918

## A Local Region of Interest Imaging Method for Electrical Impedance Tomography with Internal Electrodes

^{1}Department of Computational Science and Engineering, Yonsei University, Seoul 120-749, Republic of Korea^{2}The School of Electrical and Information Engineering, The University of Sydney, Sydney NSW 2006, Australia^{3}Impedance Imaging Research Center and Department of Biomedical Engineering, Kyung Hee University, Yongin 446-701, Republic of Korea

Received 13 November 2012; Revised 1 June 2013; Accepted 10 June 2013

Academic Editor: Bill Lionheart

Copyright © 2013 Hyeuknam Kwon et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Electrical Impedance Tomography (EIT) is a very attractive functional imaging method despite the low sensitivity and resolution. The use of internal electrodes with the conventional reconstruction algorithms was not enough to enhance image resolution and accuracy in the region of interest (ROI). We propose a local ROI imaging method with internal electrodes developed from careful analysis of the sensitivity matrix that is designed to reduce the sensitivity of the voxels outside the local region and optimize the sensitivity of the voxel inside the local region. We perform numerical simulations and physical measurements to demonstrate the localized EIT imaging method. In preliminary results with multiple objects we show the benefits of using an internal electrode and the improved resolution due to the local ROI image reconstruction method. The sensitivity is further increased by allowing the surface electrodes to be unevenly spaced with a higher density of surface electrodes near the ROI. Also, we analyse how much the image quality is improved using several performance parameters for comparison. While these have not yet been studied in depth, it convincingly shows an improvement in local sensitivity in images obtained with an internal electrode in comparison to a standard reconstruction method.

#### 1. Introduction

Electrical Impedance Tomography (EIT) is an attractive electrical properties imaging technique for medical applications due to its speed, safety, relatively low cost, and ability to display unique tissue conductivity information. Conventionally, surface electrodes are used to apply currents and record voltages from the boundary of the object of interest such as the thorax, breast, or head. These measurements are used to reconstruct the internal tissue impedances or conductivities using various image reconstruction methods where sensitivity matrix-based approaches are commonly used. Unfortunately, these measurements are insensitive to local impedance changes away from the measuring positions, whereas they are very sensitive to the boundary geometry and impedance changes near the electrodes. To improve the sensitivity and distinguishability in some local internal regions, the use of internal electrodes has been suggested by several researchers [1, 2]. These previous works were concentrated in cardiac applications of EIT as catheters are routinely introduced into various locations during cardiac monitoring, electrophysiology (EP mapping), or cardiac radiofrequency ablation (RFA).

EIT with internal electrodes could also be applied to monitoring RFA of liver tumors. RFA is widely used for the treatment of liver tumors such as hepatocellular carcinoma (HCC) and metastatic tumors since many patients are not eligible for surgical resection due to advanced disease, unfavorable location, or impaired clinical condition [3]. Numerous studies have reported that RFA is the most minimally invasive treatment for liver tumors and evaluated it as a successful therapeutic modality, providing reliable outcomes even though its posttreatment recurrence rate is higher than that of cryoablation and resection [4]. Ultrasound and magnetic resonance imaging (MRI) commonly provide the guidance of RFA because they allow real-time visualization of probe placement and high contrast images for tumor and normal tissues [5–8]. Contrast-enhanced computed tomography (CT) and MRI are widely used methods to evaluate ablated lesions by comparing the differences of lesions before and after RF ablation. However, none of these are able to continuously monitor the temperature and changes in properties of cancerous tissue and normal tissue during RFA process.

EIT does have this potential due to its high speed, sensitivity to the conductivity change of tumor and normal tissue at various frequencies [9, 10], and its ability to monitor tissue temperature in real time [11, 12]. In both cardiac EIT and EIT for monitoring liver tumour RFA, there are opportunities to add an additional internal electrode close to the region of interest (ROI) to improve EIT sensitivity to internal tissues and particularly to the tissue within the ROI. In both applications, a local imaging reconstruction method with high sensitivity within the ROI is desired since there are significant effects from conductivity changes outside the ROI such as ventilation, perfusion, movement of the lungs, diaphragm, gastric activity, and blood flow in large vessels [13].

In this paper, we propose a local sensitivity matrix-based imaging method to increase the sensitivity in a local region of interest near the internal electrode and decrease the sensitivity outside the region of interest to progress the application of EIT to cardiac monitoring and liver RFA. To further increase sensitivity to the ROI, we propose a new arrangement of surface electrodes where the electrode spacing is not equal but more near the ROI. In order to evaluate the performance of the local ROI imaging method, we perform numerical simulations and phantom experiments.

#### 2. Method

##### 2.1. Conventional EIT Reconstruction Method

Let a bounded domain represent the subject to be imaged. Surface electrodes for are attached to the boundary , where is the total number of electrodes. Let denote the complex conductivity at time , angular frequency of , and position . When we inject a sinusoidal current at an angular frequency of between a chosen pair of electrodes, a voltage distribution is formed at the position . Here, is used for expressing a slow-time change in the complex conductivity distribution, and is used for the fast-time change to represent time-harmonic fields. The induced time-harmonic potential satisfies the following elliptic partial differential equation [14]: where is the outward unit normal vector on and is the corresponding Neumann data on due to the injection current.

Static imaging in EIT is difficult due to its fundamental limitations in handling boundary geometry and uncertainty in electrode position. Time-difference EIT (tdEIT) and frequency-difference EIT (fdEIT) use time and frequency difference data, respectively, so that the data subtraction can effectively cancel out common errors related to boundary geometry [14–18]. The difference imaging in EIT is based on linear approximations of the following identities: (i)time-difference EIT , (ii)frequency-difference EIT .

Let the angular frequency be fixed. Let denote the time-harmonic potential due to th injection current between the adjacent pair of electrodes and . The boundary voltage between and due to the th injection current can be approximated as where is the surface area of the electrode. The last identity in (2) comes from (1) and divergence theorem. We collect number of boundary voltage data for a sequence of time : where is the transpose. Here, any index number must be understood as a modulus of .

The time-difference data due to the time change of has the following relation: The linearized method is based on the following rough approximation: where is the potential of (1) corresponding to a reference conductivity .

Discretizing the domain of interest into pixels and assuming that is constant on each pixel , the time-difference EIT problem of (5) can be changed to solve the following linear system: The th column of the sensitivity matrix is Hence, the EIT problem is to find a best linear combination of column vectors which produces : The column vector of represents sensitivity of current-voltage data at the fixed pixel , whereas row vectors of represent sensitivity distribution for a fixed current-voltage data. We refer to this approach (combined with a regularization) to reconstruct images as the conventional method throughout this paper.

##### 2.2. Local-ROI Imaging Method for EIT

The local-ROI imaging method for EIT is to provide the image of the conductivity change in a local region of interest (ROI) instead of the image in the entire domain . Let a domain be the local region of interest to be imaged. Imagine that the measured data in (5) is divided into two parts where is the voltage change in response to the conductivity perturbation in the local ROI and is the voltage change in response to in . With proper arrangement, we may assume that the first column vectors of are sensitivity vectors to pixels in ROI and the other column vectors are sensitivity vectors to pixels in . Then the sensitivity matrix can be decomposed into where represents a proper size of zero matrix, If we could extract the data by filtering out , the global problem (6) can be changed into the local problem: where . Figure 1(c) shows the reconstructed image of using the localized linear system (12) via numerical simulation. Comparing this local image with the standard EIT image reconstruction shown in Figure 1(b), it would be desirable to filter out to enhance image resolution.

For the local imaging in the ROI , we aim to develop a method of extracting from the full data . In order to eliminate the unrelated data in the linear system (6), we need to find an optimal matrix such that Here, gives a matrix at which is minimized. If satisfies , then it eliminates the unrelated data and we get the localized linear system corresponding to (12):

Note that the quantity depends on the electrode configuration and mesh structure that determines the structure of column vectors of as shown in Figure 2. The may not be small when the sensitivity matrix is highly ill-conditioned. If the column vectors are orthogonal to , then by choosing whose rows consist of the column vectors . But, this is not possible with the standard EIT electrode configuration. We try to find an optimal which minimizes . Indeed, cross-correlation of column vector in is big if is correlated with column vectors in ROI; cross-correlation is defined by Figure 3 shows that cross-correlation values decrease by placing internal electrode.

For finding proper , we propose the following minimization where is a suitable parameter. We should note that each is designed to be close and to be parallel to while orthogonal to for each . The first term in (16), , is small when is orthogonal to . The second term in (16), , is small if is parallel to . We define a matrix whose columns are consisted of : We multiply to the linearized system (6): Now, we have the linear system (18) with the modified sensitivity matrix with the modified data . Here, is regarded as a rough approximation of .

#### 3. Numerical Simulations

In order to analyse the boundary electrode position and the benefit of using an internal electrode with the proposed local-ROI imaging method, we prepared three different kinds of electrode configuration and applied the conventional and local-ROI imaging methods explained in the previous sections. Figure 4 shows the cylindrical phantoms with three different electrode configurations. Data obtained from all three electrode configurations were processed by the conventional TSVD reconstruction method. The local-ROI imaging method uses an internal electrode so it was only applied to data obtained with the two internal electrode configurations (Figures 4(b) and 4(c), (Models 1 and 2)). We carry out a total of five numerical simulations, the standard method on Model 0, 1, 2 and the local-ROI method on models 1, 2.

To compare sensitivity and robustness to noise in the suggested five cases, we simulated an object of 2 S/m conductivity with 0.1428 diameter at (0.8, 0) in the ROI of 1 S/m saline tank. The radius of the ROI is defined by 5/6 of distance between internal electrode and closest boundary electrode. We used two performance indexes to assess the improvement when using an internal electrode and the proposed linear system (18). First, we computed the singular value threshold to produce the same conductivity contrast of image as shown in Figure 5(a). We repeatedly reconstruct images, updating the singular value threshold, until the same conductivity contrast is produced between the anomaly and background in the reconstructed image. Then we compared the singular value thresholds as a lower number of threshold is an indication of better noise robustness. When we used an internal electrode, the truncated singular value threshold was lower than that without using the internal electrode. Also, the performance of local-ROI imaging method was improved when it was used with unequally spaced surface electrodes with denser spacing near by the ROI. Second, we are concerned about the effect of high contrast anomalies outside the ROI. To investigate this, we examined the sensitivity values of and for each simulation case. The normalized sensitivity values of and from within the ROI and out of ROI are shown separately in Figure 5(b). There was an improved relative sensitivity to the ROI region with the proposed method . To show sensitivity values, we calculated the ratio of the matrix norm and for inside ROI and outside ROI respectively, where of Model 0, 1, 2, of Model 1, 2}, is a submatrix of corresponding to ROI, is a submatrix of corresponding to outside ROI, .Note that a matrix norm shows how much deforms . So, and show maximum influence of submatrices and for multiplication of .

In order to compare the performance of the imaging methods in the reconstructed images, we placed multiple objects which had the same conductivity of 0.029 S/m with 0.01 diameter at (0.02, 0.015), (0.02, −0.015), (0.05, −0.015), and (0.05, −0.015) within the ROI of the 0.0418 S/m saline tank as shown in Figure 6(a). The conventional reconstructed image without an internal electrode was produced by the TSVD algorithm in Model 0 with diagonal current injection (Figure 6(b)). For comparison of the imaging methods using an internal electrode, we positioned an internal electrode at (0.035, 0) in Model 2. Figures 6(b) and 6(c) show images using the conventional TSVD reconstruction algorithm and the local-ROI imaging method. The proposed method with an internal electrode has better sensitivity and detectability in the ROI.

#### 4. Experimental Results

To evaluate the performance of the imaging methods in a physical model, we prepared a cylindrical saline tank with the same geometry, ROI, and electrode positions as Figures 4(a) and 4(c) since the denser clustering of boundary electrodes near the ROI (Figure 4(c)) produced better results than the conventional equally spaced electrodes (Figure 4(b)) in the simulation study. The diameter and height of saline tank were 14 cm and 6 cm, respectively. The boundary and internal electrodes were located 3 cm from the bottom in the direction. The and positions of the internal electrode were (3.5, 0) cm relative to the origin (0, 0) at the center of cylindrical tank. The diameter of internal electrode was 0.25 cm and it was covered with insulated rubber except for the end piece of exposed metal. The conductivity of saline was 0.042 S/m. All data was measured by the KHU Mark2.5 EIT system operated at 10 kHz [19].

We evaluated the imaging methods using two different situations. First, we located a piece of radish inside of ROI (case 0, 1, 2, 3 (5.6, 0)) and a piece of potato outside of ROI (case 1 (−5.6, 0), case 2 (−3.5, −4), case 3 (0, −5.6)) to assess how the sensitivity of each method is influenced by a high conductivity object in the surrounding area of ROI. The diameter and height of both objects were the same at 1.2 cm and 7 cm as shown in Figure 7. The conductivity of radish (0.038 S/m at 10 kHz) was 9.5% lower than saline and the potato (0.029 S/m at 10 kHz) was a higher contrast than the radish with 31% lower conductivity than saline. All data was obtained in the new configuration with an internal electrode and unequal surface electrode spacing as Model 2 in Figure 4(c). All measured data was processed by the conventional TSVD method and the local-ROI imaging method for comparison. We present the norm values of the difference of and compared to those values of reference case 0 because it is difficult to analyze images directly from the artifacts due to the large anomaly outside the ROI. Here, we processed reference value by conventional method in case 0 which included only one anomaly (radish) in the ROI. Figure 8(a) shows the effect of the high conductivity contrasted object located outside of ROI. There is less effect with the new local-ROI imaging method when the object outside the ROI is in the opposite hemicircle to the ROI. As the object approaches the ROI (case 1 to case 3) both methods perform similarly.

Figure 8(b) shows the singular value threshold required to produce the same conductivity contrast for each case. The local-ROI imaging method had similar threshold values for all cases with less than 1.73% variation. However the conventional TSVD method showed larger singular values and large dependence on the high contrast anomaly position.

Secondly, we evaluated the methods from the reconstructed images in the same configuration as the numerical simulation. We put four carrot objects with 0.04 S/m conductivity at 10 kHz in the 0.06 S/m saline tank. The position of each object was (1.5, 2), (5.5, 2), (5.5, −2), and (1.5, −2) cm, respectively. All objects had the same diameter of 2 cm.

We obtained an image in configuration Model 0 during diagonal current injection and TSVD reconstruction. For algorithm comparison, we measured data in Model 2 and applied both methods separately. Figures 9(a)–9(c) show the reconstructed images. The proposed local-ROI imaging method with an internal electrode can distinguish inner objects better than conventional methods.

#### 5. Conclusions and Discussions

The new local-ROI imaging method has been shown to improve the sensitivity in the ROI region and robustness to noise by comparison of sensitivity matrix values. It provided the approximated linear system with optimized sensitivity matrix to emphasize the detection in ROI. The ROI region is close to an internal electrode and surface electrodes were unequally spaced being more dense close to the ROI. The new method and setup is also less affected by the high conductivity contrasted object outside of ROI. These objects are mimicking a large vessel and several regions of hepatic and metastatic cancers treated by RFA, other gastric organs, or the lungs in cardiac monitoring.

Some performance indexes proposed in this paper described that the sensitivity and detectability were obviously improved in simulation and experiment results. However, the reconstructed images did not show improvements as dramatic as we initially expected. One of the major reasons was the ill-posedness of imaging problem. When we designed the optimal matrix that satisfied condition (16), it could not eliminate the effect of unrelated data completely. Also, the quality of the experimental reconstructed image was highly dependent on the position of an internal electrode and boundary electrodes because we placed the electrodes with denser spacing near by the ROI. Encouragingly with the pilot results of the local-ROI imaging method using an internal electrode, it shows the feasibility and suggests a new approach to improve the resolution of internal local region.

An additional improvement which was introduced in this paper was the nonequidistant spacing of electrodes, with the electrodes more densely spaced near the ROI in Model 2 (Figure 4(c)). This showed improved results in simulation and so was used in the experimental setup; however, this may be more sensitive to electrode positions, the location of the ROI, or noise. Interestingly in Figure 5(a) we found that this setup of non-equidistant electrodes (Model 2) showed less sensitivity with the conventional EIT method. This may reflect the preference for symmetry with EIT performed in a circular object. The single ROI we investigated was off centre and once we use our local-ROI focusing algorithm there is an improvement found in the non-equidistant electrode setup (model 2).

The application of EIT to RF ablation monitoring has high potential; however, we need to study more how to separate the conductivity variations which arise simultaneously from temperature changes and tissue property changes. While the temperature coefficient of conductivity in electrolytes of per degree is well known for electrolytes, tissues also exhibit a conductivity dependence on temperature-induced fluid volume shifts which is of the same order of magnitude [20].

Lower frequency measurements could be used to discriminate ablated and nonablated tumor and normal tissues. Liver tumour tissue has a higher conductivity than normal liver tissue over 10 Hz to 1 MHz shown by four terminal impedance measurements in excised tissue. Following an ablation, both tissue types showed significantly increased conductivity over the same frequency range indicating that electrical impedance may be used to differentiate tumor tissue diagnostically, for ablation planning and postablation assessment [10].

The suggested method may be applied in 3D domain without any additional changes to the algorithm. Here we used a 2D domain analysis because applying the suggested method requires a lot of computations on the column vectors of sensitivity matrix corresponding to pixels in ROI [21]. We need to optimize the local-ROI imaging algorithm by focusing on the ROI and investigate ways to increase the coverage area without moving the electrode position. We may combine the frequency-difference EIT [22, 23] with the local-ROI imaging method to apply frequency-difference data.

#### Acknowledgment

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIP) (2010-0018275) and Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (2011-0022779).

#### References

- B. M. Eyiiboglu and T. C. Pilkington, “Intracavity electrical impedance tomography,” in
*Proceedings of the 12th Annual International Conference of the IEEE on Engineering in Medicine and Biology Society*, p. 130, 1990. View at Publisher · View at Google Scholar - J. N. Tehrani, C. Jin, and A. L. McEwan, “Modelling of an oesophageal electrode for cardiac function tomography,”
*Computational and Mathematical Methods in Medicine*, vol. 2012, Article ID 585786, 10 pages, 2012. View at Publisher · View at Google Scholar · View at Scopus - A. R. Gillams, “The use of radiofrequency in cancer,”
*British Journal of Cancer*, vol. 92, no. 10, pp. 1825–1829, 2005. View at Publisher · View at Google Scholar · View at Scopus - E. K. Abdalla, J.-N. Vauthey, L. M. Ellis et al., “Recurrence and outcomes following hepatic resection, radiofrequency ablation, and combined resection/ablation for colorectal liver metastases,”
*Annals of Surgery*, vol. 239, no. 6, pp. 818–827, 2004. View at Publisher · View at Google Scholar · View at Scopus - N. J. McDannold and F. A. Jolesz, “Magnetic resonance image-guided thermal ablations,”
*Topics in Magnetic Resonance Imaging*, vol. 11, no. 3, pp. 191–202, 2000. View at Publisher · View at Google Scholar · View at Scopus - S. Terraz, A. Cernicanu, M. Lepetit-Coiffé et al., “Radiofrequency ablation of small liver malignancies under magnetic resonance guidance: progress in targeting and preliminary observations with temperature monitoring,”
*European Radiology*, vol. 20, no. 4, pp. 886–897, 2010. View at Publisher · View at Google Scholar · View at Scopus - S. N. Goldberg, G. S. Gazelle, and P. R. Mueller, “Thermal ablation therapy for focal malignancy: a unified approach to underlying principles, techniques, and diagnostic imaging guidance,”
*American Journal of Roentgenology*, vol. 174, no. 2, pp. 323–331, 2000. View at Google Scholar · View at Scopus - S. N. Goldberg, C. J. Grassi, J. F. Cardella et al., “Image-guided tumor ablation: standardization of terminology and reporting criteria,”
*Radiology*, vol. 235, no. 3, pp. 728–739, 2005. View at Publisher · View at Google Scholar · View at Scopus - S. Laufer, A. Ivorra, V. E. Reuter, B. Rubinsky, and S. B. Solomon, “Electrical impedance characterization of normal and cancerous human hepatic tissue,”
*Physiological Measurement*, vol. 31, no. 7, pp. 995–1009, 2010. View at Publisher · View at Google Scholar · View at Scopus - D. Haemmerich, D. J. Schutt, A. W. Wright, J. G. Webster, and D. M. Mahvi, “Electrical conductivity measurement of excised human metastatic liver tumours before and after thermal ablation,”
*Physiological Measurement*, vol. 30, no. 5, pp. 459–466, 2009. View at Publisher · View at Google Scholar · View at Scopus - F. Ferraioli, A. Formisano, and R. Martone, “Effective exploitation of prior information in electrical impedance tomography for thermal monitoring of hyperthermia treatments,”
*IEEE Transactions on Magnetics*, vol. 45, no. 3, pp. 1554–1557, 2009. View at Publisher · View at Google Scholar · View at Scopus - H. Griffiths and A. Ahmed, “Applied potential tomography for non-invasive temperature mapping in hyperthermia,”
*Clinical Physics and Physiological Measurement*, vol. 8, no. 4, supplement A, article 019, pp. 147–153, 1987. View at Publisher · View at Google Scholar · View at Scopus - D. S. K. Lu, S. S. Raman, P. Limanond et al., “Influence of large peritumoral vessels on outcome of radiofrequency ablation of liver tumors,”
*Journal of Vascular and Interventional Radiology*, vol. 14, no. 10, pp. 1267–1274, 2003. View at Google Scholar · View at Scopus - D. Holder,
*Electrical Impedance Tomography: Methods, History and Applications*, IOP, Bristol, UK, 2005. - P. Metherall, D. C. Barber, R. H. Smallwood, and B. H. Brown, “Three-dimensional electrical impedance tomography,”
*Nature*, vol. 380, no. 6574, pp. 509–512, 1996. View at Publisher · View at Google Scholar · View at Scopus - M. Cheney, D. Isaacson, and J. C. Newell, “Electrical impedance tomography,”
*SIAM Review*, vol. 41, no. 1, pp. 85–101, 1999. View at Google Scholar · View at Scopus - J. K. Seo, J. Lee, S. W. Kim, H. Zribi, and E. J. Woo, “Frequency-difference electrical impedance tomography (fdEIT): algorithm development and feasibility study,”
*Physiological Measurement*, vol. 29, no. 8, pp. 929–944, 2008. View at Publisher · View at Google Scholar · View at Scopus - J. K. Seo and E. J. Woo,
*Nonlinear Inverse Problems in Imaging*, Wiley Press, New York, NY, USA, 2012. - H. Wi, H. Sohal, T. I. Oh, A. L. McEwan, and E. J. Woo, “Multi-frequency electrical impedance tomography system with automatic self-calibration for long-term monitoring: KHU Mark2. 5,”
*IEEE Transactions on Biomedical Circuits and Systems*. In press. - E. Gersing, “Monitoring temperature-induced changes in tissue during hyperthermia by impedance methods,”
*Annals of the New York Academy of Sciences*, vol. 873, pp. 13–20, 1999. View at Publisher · View at Google Scholar · View at Scopus - W. R. B. Lionheart, “Uniqueness, shape, and dimension in EIT,”
*Annals of the New York Academy of Sciences*, vol. 873, pp. 466–471, 1999. View at Publisher · View at Google Scholar · View at Scopus - S. C. Jun, J. Kuen, J. Lee, E. J. Woo, D. Holder, and J. K. Seo, “Frequency-difference EIT (fdEIT) using weighted difference and equivalent homogeneous admittivity: validation by simulation and tank experiment,”
*Physiological Measurement*, vol. 30, no. 10, pp. 1087–1099, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. Kuen, E. J. Woo, and J. K. Seo, “Multi-frequency time-difference complex conductivity imaging of canine and human lungs using the KHU Mark1 EIT system,”
*Physiological Measurement*, vol. 30, no. 6, pp. S149–S164, 2009. View at Publisher · View at Google Scholar · View at Scopus