PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
 
Med Biol Eng Comput. Author manuscript; available in PMC 2010 June 1.
Published in final edited form as:
PMCID: PMC2828869
NIHMSID: NIHMS169617

Ultrasonically-Induced Lorentz Force Tomography

Abstract

Electrical conductivity can be measured using the ultrasonically-induced Lorentz force. An ultrasonic wave is passed through tissue in the presence of a magnetic field. Moving charges in a magnetic field are subject to the Lorentz force, which acts as the source of current and potential. This paper shows that ultrasonically-induced Lorentz force imaging can be formulated in a way that makes it similar to tomography: an image can be reconstructed using waves propagating in various directions. More specifically, measuring the dipole strength for a particular direction and wavelength is equivalent to measuring the Fourier transform of the conductivity distribution at one point in frequency space. Measurements at a variety of wavelengths and directions are equivalent to mapping the Fourier transform of the conductivity distribution. The conductivity can then be found by an inverse Fourier transform.

Keywords: Ultrasound, Conductivity, Lorentz Force, Tomography, Fourier Transform

Introduction

Wen et al. [7] and Montalibet et al. [3,4] developed a technique to image electrical conductivity using the ultrasonically-induced Lorentz force. An ultrasonic wave is passed through tissue in the presence of a magnetic field. The wave causes the tissue--and the charged ions within the tissue--to move. A moving charge in a magnetic field is subject to the Lorentz force [1]. This force causes positive and negative charge to separate, and acts as the source of current and potential. Montalibet et al. [3] were able to detect the potential produced by a 1.5 MPa, 500 kHz ultrasonic wave in a 0.35 T magnetic field. They used these measurements to produce a one-dimensional image of conductivity gradients in biological tissue.

The goal of this paper is to analyze and extend this method to two or three dimensions, and to perform a preliminary calculation to prove that the method works in principle. We show that ultrasonically-induced Lorentz force imaging can be formulated in a way that makes it similar to tomography: reconstructing an image using waves propagating in various directions. More specifically, measuring the dipole strength for a particular direction and wavelength is equivalent to measuring the Fourier transform of the conductivity distribution at one point in frequency space. Measurements at a variety of wavelengths and directions are equivalent to mapping the Fourier transform of the conductivity distribution throughout frequency space, allowing us to determine the conductivity as a function of position by an inverse Fourier transform.

Methods

Let an ultrasound planar wave pass through tissue in the direction x. At one instant of time, the velocity u(x) of the tissue in the direction of propagation is

u(x)=uoeikxi^,
1

where uo is the amplitude of the wave and k is the spatial frequency. Assume the tissue has a conductivity distribution given by

σ(x,y,z)=σo+σ(x,y,z),
2

where σo is a uniform background conductivity and σ′ (x, y, z) specifies a localized spatial variation about this background. A uniform, steady magnetic field is in the z-direction, B = Bo k^ . The movement of charge in the magnetic field will give rise to a Lorentz force, producing a current density J L , which by the right-hand-rule is directed in the y-direction

JL=uoBoeikx(σo+σ)j^.
3

If separation of charge occurs, an Ohmic current density J C is also present

JC=(σo+σ)V,
4

where V is the induced voltage. Continuity of total current, J L + J C , implies

·((σo+σ)V)=·JL.
5

In a completely homogeneous conductivity distribution (σ′ = 0), there would be currents J L inside the object having a strength that is proportional to the conductivity (as opposed to the gradient of the conductivity). These strictly solenoidal currents (zero divergence) have nonzero curl ([nabla] ×). In an inhomogeneous conductivity distribution, charge accumulates where the gradient of the conductivity is non-zero. The electric field produced by these charges is purely irrotational (has zero curl), but does have a divergence ([nabla] ·), leading to Eq. 5.

If the deviations of the conductivity distribution are small compared to the background, then we can ignore σ′ compared to σo on the left-hand-side of Eq. 5, obtaining Poisson's equation

2V=uoBoσoσyeikx.
6

Following a standard multipole expansion for the potential [1], the dipole term is given by

V=14πσop·r^r2,
7

where the dipole is

p(k)=uoBorσyeikxdxdydz.
8

If we integrate by parts in the y-direction and assume that σ′ is localized, so it goes to zero at large distances, then the y-component of the dipole strength becomes

py(k)=uoBoeikx[σdydz]dx.
9

In this paper, we consider only a 2-dimensional slab of tissue that is uniform and of thickness δ in the z direction, giving

py(k)=uoBoδeikx[σ(x,y)dy]dx.
10

The integral in x has the form of a Fourier transform.

The two-dimensional Fourier transform of the conductivity distribution is defined as

σ^(kx,ky)=σ(x,y)ei(kxx+kyy)dxdy.
11

If we let ky = 0, this becomes

σ^(kx,0)=eikxx[σ(x,y)dy]dx.
12

Thus, comparing Eqs. 10 and 12, we find that measurement of py(k) provides information about the two-dimensional Fourier transform of the conductivity distribution at the point in frequency space kx = k and ky = 0.

This analysis is very similar to the method of two-dimensional image reconstruction from projections by Fourier transform, often used to reconstruct computed tomography or magnetic resonance images [2]. By analogy to CT, there is nothing special about the direction that we chose for the ultrasonic wave. Applying the ultrasonic wave in another direction and measuring the dipole moment perpendicular to that direction will provide information about the conductivity distribution along another line in frequency space. By taking measurements along many directions and using different wavelengths, we can map out the entire frequency plane. Once this information is known, a two-dimensional inverse Fourier transform will yield an image of the conductivity distribution.

Results

We tested this algorithm by a four-step procedure: 1) calculate V(r,[var phi]) numerically using Eq. 6 for various values of k and θ, 2) estimate the dipole strength p(k,θ) from the calculated V in each case, 3) use a bilinear interpolation to obtain p(kx, ky) from p(k,θ), and 4) perform a 2D inverse Fourier transform of p(kx, ky) to obtain the conductivity σ′(x,y).

1. Calculate V(r,[var phi]) numerically

We assume the tissue is circular, with radius 20 mm (Fig. 1a). The outer boundary of the tissue is sealed. The potential is solved using a finite difference approximation of Eq. 6 in polar coordinates (r, [var phi]), and overrelaxation. The space step is Δr = 0.2 mm and Δ[var phi] = 2° (200 points in the r-direction, and 180 points in the [var phi]-direction).

Figure 1
a) Schematic diagram of the tissue. The arrow in the direction θ shows the direction of the ultrasonic wave. The circle has a radius of 20 mm, and the dots indicate where the potential difference is measured. b) The conductivity distribution given ...

The source term in Eq. 6 depends on the values k and θ, and on the conductivity distribution σ′(x,y). As an example, consider the case of a Gaussian distribution for σ′ (Fig. 1b),

σ(x,y)=Cexp[x2+(yyo)2a2],
13

where C = 1 S m−1, yo = 3 mm, and a = 3 mm. Figure 1c shows the real part of the source term for k =1.18 mm−1 (corresponding to a wavelength of 5.33 mm) and θ = 50°. This source term results in the potential distribution in Fig. 1d.

In order to provide a more challenging test of the algorithm, we consider a conductivity distribution σ′(x,y) that is the sum of five Gaussians

σ(x,y)=n=15Cnexp[(xxn)2+(yyn)2an2],
14

where the constants Cn, an, xn, and yn are given in Table 1. Figure 2a plots σ′(x,y). We calculate the potential for 16 spatial frequencies ki = (i−1) Δk, where Δk = 0.393 mm−1, in 18 directions θj = (j−1) Δθ, where Δθ = 10°.

Figure 2
a) The conductivity distribution given by Eq. 14, using the parameters in Table 1. b) The conductivity distribution calculated using the algorithm presented in this paper, with 16 values of k and 18 values of θ. In both cases, the gray scale extends ...
Table 1
Parameters used in the conductivity distribution (Eq. 14)

2. Estimate the Dipole Strength

The potential distribution in Fig. 1d has two contributions: a part that has a spatial distribution very similar to the source function, and a part more distant from the source that has approximately a dipole pattern. Our algorithm requires determination of the dipole strength. We obtain a measure of the relative strength of the dipole by measuring the potential difference between opposite sides of the circular tissue, perpendicular to the direction of propagation of the ultrasonic wave (the dots in Fig. 1a). These measurements yield dipole values p(k, θ) at each k and θ.

3. Bilinear Interpolation

In order to perform an inverse Fourier transform, we need the dipole strength for a regular array of kx and ky values. Therefore, we use a bilinear interpolation to obtain p(kx, ky) from p(k,θ), for 32 values (for both positive and negative spatial frequencies) using Δkx = Δky = 0.393 mm−1. Values of kx and ky corresponding tok=kx2+ky2greater than 5.89 mm−1 are outside the range measured and are therefore set to zero.

4. Inverse Fourier Transform

Once we have p(kx, ky), we take the inverse Fourier transform to obtain σ′(x,y) [5]. The result is shown in Fig. 2b.

Discussion

The qualitative agreement between Figs. 2a and 2b is proof that in principle our method can be used to image the conductivity distribution. The method would be subject to many problems that arise when using Fourier techniques: aliasing, ringing at sharp discontinuities of conductivity, limitations on the spatial resolution caused by an upper limit of the values of the spatial frequency, etc. These problems are not immediately apparent in Fig. 2, because the conductivity distribution we consider is smooth, and therefore does not contain a lot of high frequency information.

Proof in principle does not imply that the method will be practical. We foresee several limitations of this method, some evident in our calculation and some not.

  1. We did not calculate the conductivity distribution σ′, but instead a quantity proportional to it. In other words, we did not obtain the proportionality constant. This is not a fundamental limitation, but arises because we use the potential difference across the tissue (ΔV between the two dots in Fig. 1a) as a relative measure of the dipole strength. In principle, one could fit the measured potential to the dipole distribution and obtain the dipole strength, as is commonly done in electroencephalography and magnetoencephalography studies [2]. Once the dipole strength is known, one could divide by the magnetic field strength, the amplitude of the ultrasonic wave, and the tissue thickness (Eq. 10) and thereby obtain the actual conductivity distribution σ′, including the scale factor.
  2. We assumed that the conductivity distribution could be divided into two parts: a uniform background conductivity and a spatially varying part that is localized. The value of the background conductivity did not play any role in our calculation (although it would enter the calculation if determining the scale factor, because it arises in Eq. 7 relating the potential to the dipole moment). If the spatially varying part is not localized, near field contributions may contribute to the measured potential, making it difficult to relate the potential to the dipole strength. Being able to obtain the dipole moment from the potential is therefore crucial for the method. Interestingly, this method should be more accurate if the potential is measured farther from the source, ensuring a purer dipole field, than if measured nearer. Apparently a radius of 20 mm in our calculation was sufficient to image a conductivity distribution that was localized within a radius of about 8 mm.
  3. Our calculation was for a two-dimensional conductivity distribution. Equation 7 is appropriate for a three-dimensional geometry, so any attempt to use Eq. 7 to relate the dipole strength to the measured potential would require a full 3D analysis. One could obtain a 3D image either by using an ultrasonic wave that is focused and localized in the z direction (analogous to MRI or CT scans), or using plane waves in all three directions to obtain the 3D Fourier transform.
  4. We assume that the ultrasonic wave is a plane wave with no change in amplitude as it propagates. Clearly the mechanical properties of the tissue could reflect or scatter the wave: this is in fact the basis for ultrasonic imaging [2]. Moreover, as frequency increases attenuation of the ultrasonic wave becomes greater. Such modifications to the ultrasonic wave will need to be studied before ultrasonically-induced Lorentz force tomography could become practical. In many cases, the transmitted wave is not dramatically affected by reflections. For instance, the acoustic impedance of muscle is 1.23 times the acoustic impedance of fat, implying that 99% of the incident energy of a 1 MHz ultrasonic wave on a muscle-fat boundary is transmitted and only 1% is reflected [2]. (Other tissues, such as bone, will have a much larger effect.) One can speculate about ways around this issue. If an ultrasound image was obtained without a magnetic field present, the mechanical effects on the signal might be determined and somehow subtracted out of the conductivity image. One way to do this might be to use the measured ultrasonic wave amplitude for u(x) in Eq. 1 instead of a plane wave. Moreover, if attenuation is exponential along the path of the wave, its effect might be incorporated into the analysis in a way similar to how the exponential absorption of x-rays is accounted for in CT reconstruction.
  5. This method of image construction assumes that both the amplitude and the phase of the Fourier transform can be measured. Obtaining the phase information may be difficult, because it assumes that the phase of the ultrasonic wave is known at the time that the dipole is measured. At the least, this will require careful timing of the measurement with respect to the wave.
  6. The method uses continuous wave, rather than pulsed wave, ultrasound. It is essential for our method that the wave be present simultaneously throughout the tissue. This approach is different than that proposed by Wen et al. [7] and Montalibet et al. [3, 4], who both used pulsed ultrasound to record reflections at conductivity gradients in a manner more analogous to traditional ultrasound imaging.
  7. Tissue anisotropy can produce additional effects during Lorentz force imaging that are not included in our analysis [6].
  8. We did not include effects of noise in our simulation. However, Fourier methods for image reconstruction are used widely and have been studied extensively, particularly in magnetic resonance imaging [2]. Noise is a problem in these applications, but not an insurmountable one.

Despite these limitations, we believe that our analysis provides a valuable proof-of-principle that ultrasonically-induced Lorentz force tomography is possible. Wen et al. [7] and Montalibet et al. [3,4] showed that the method is feasible experimentally, and our computations (Fig. 2) show that under a somewhat idealized geometry the method works. Further research is necessary to determine if this new type of tomography will be useful in practice and what are the relative strengths and weaknesses of the method compared to other techniques to measure electrical conductivity.

Acknowledgments

This research was supported by the National Institutes of Health grant R01EB008421, and by the SMaRT program at Oakland University, a Research Experience for Undergraduates (REU) that was funded by NSF grant DMR-055 2779.

References

1. Griffiths DJ. Introduction to Electrodynamics. 3rd Edition. Vol. 1999. Upper Saddle River, NJ: Prentice Hall;
2. Hobbie RK, Roth BJ. Intermediate Physics for Medicine and Biology. 4th Edition. New York: Springer; 2007.
3. Montalibet A, Jossinet J, Matias A. Scanning electric conductivity gradients with ultrasonically-induced Lorentz force. Ultrasonic Imag. 2001;23:117–132. [PubMed]
4. Montalibet A, Jossinet J, Matias A, Cathignol D. Electric current generated by ultrasonically induced Lorentz force in biological media. Med Biol Eng Comput. 2001;39:15–20. [PubMed]
5. Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes in Fortran 77: The Art of Scientific Computing. 2nd Edition. New York: Cambridge Univ Press; 1992.
6. Tseng N, Roth BJ. The potential induced in anisotropic tissue by the ultrasonically-induced Lorentz force. Med Biol Eng Comput. 2008;46:195–197. [PubMed]
7. Wen H, Shah, Balaban RS. Hall effect imaging. IEEE Trans Biomed Eng. 1998;45:119–124. [PMC free article] [PubMed]