AU783991B2 - Asymmetric radio frequency coils for magnetic resonance - Google Patents

Asymmetric radio frequency coils for magnetic resonance Download PDF

Info

Publication number
AU783991B2
AU783991B2 AU69966/01A AU6996601A AU783991B2 AU 783991 B2 AU783991 B2 AU 783991B2 AU 69966/01 A AU69966/01 A AU 69966/01A AU 6996601 A AU6996601 A AU 6996601A AU 783991 B2 AU783991 B2 AU 783991B2
Authority
AU
Australia
Prior art keywords
radio frequency
magnetic resonance
coil
field
current density
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired
Application number
AU69966/01A
Other versions
AU6996601A (en
Inventor
Stuart Crozier
David Michael Doddrell
Benjamin Lawrence
Kurt Luescher
Wolfgang Udo Roffmann
Desmond Yau
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NMR Holdings No 2 Pty Ltd
Original Assignee
NMR Holdings No 2 Pty Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from AUPR0059A external-priority patent/AUPR005900A0/en
Application filed by NMR Holdings No 2 Pty Ltd filed Critical NMR Holdings No 2 Pty Ltd
Priority to AU69966/01A priority Critical patent/AU783991B2/en
Publication of AU6996601A publication Critical patent/AU6996601A/en
Application granted granted Critical
Publication of AU783991B2 publication Critical patent/AU783991B2/en
Priority to AU2006201449A priority patent/AU2006201449B2/en
Anticipated expiration legal-status Critical
Expired legal-status Critical Current

Links

Landscapes

  • Magnetic Resonance Imaging Apparatus (AREA)

Description

AUSTRALIA
Patents Act 1990 COMPLETE SPECIFICATION FOR A STANDARD PATENT Name of Applicant: Actual Inventors: Address for Service: NMR Holdings No. 2 Pty Limited Stuart Crozier Benjamin Lawrence Desmond Yau Kurt Luescher Wolfgang Udo Roffmann David Michael Doddrell CULLEN CO., Patent Trade Mark Attorneys, 239 George Street, Brisbane, QId. 4000, Australia.
Asymmetric Radio Frequency Magnetic Resonance Invention Title: Coils for Details of Associated Provisional Applications: Australian Patent Application No.
PR0059 filed 11 September 2000 The following statement is a full description of this invention, including the best method of performing it known to us: 2 ASYMMETRIC RADIO FREQUENCY COILS FOR MAGNETIC RESONANCE FIELD OF THE INVENTION This invention relates to radio frequency coils for magnetic resonance (MR) applications. In particular, the invention is directed to asymmetric radio frequency coils for magnetic resonance imaging (MRI) machines.
The radio frequency coils of the invention may be used for transmitting a radio frequency field, receiving a magnetic resonance signal, or both transmitting a radio frequency field and receiving a magnetic resonance signal. When the radio frequency coil serves a transmitting function, it will normally be combined with a shield to reduce magnetic interference with external components of the magnetic resonance imaging system.
BACKGROUND OF THE INVENTION In magnetic resonance imaging (MRI) applications, a patient is placed in a strong and homogeneous static magnetic field, causing the otherwise randomly oriented magnetic moments of the protons, in water molecules within the body, to precess around the direction of the applied field. The part of the body in the homogeneous region of the magnet is then irradiated with radio-frequency (RF) energy, causing some 20 of the protons to change their spin orientation. The net magnetization of the spin ensemble is nutated away from the direction of the applied static magnetic field by the applied RF energy. The component of this net magnetization orthogonal to the direction of the applied static magnetic field acts to induce measurable signal in a receiver coil tuned to the frequency of precession. This is the magnetic resonance (MR) signal.
The useful RF components are those generated in at plane at 90 degrees to the direction of the static magnetic field. The same coil structure that generates the RF field can be used to receive the MR signal or a separate receiver coil placed close to the patient may be used. In either case the coils are tuned to the Larmor precessional frequency coo where o yB 0 and y is the gyromagnetic ratio for a specific nuclide and
B
0 is the applied static magnetic field.
A desirable property of radio frequency coils for use in MR is the generation of homogeneous RF fields over a prescribed region. Normally this region is central to the coil structure for transmission resonators. A well known example of transmission 3 resonators is the birdcage resonator, details of which are given by Hayes et. al. in The Journal of Magnetic Resonance, 63, 622 (1985) and U.S. Pat. No. 4,694,255.
In some circumstances it is desirable to generate a target field over an asymmetric region of the coil structure, a region that is asymmetric relative to the mid-length point of the longitudinal axis of the coil structure. This is potentially advantageous for patient access, conformation of the coil structure to the local anatomy of the patient and for use in asymmetric magnet systems.
One method that is known in the art for generating homogeneous fields over a volume that is asymmetric to the coil structure is to enclose one end of the cylindrical structure, a so-called 'end-cap' or dome structure (details of which are given by Meyer and Ballon in The Journal of Magnetic Resonance, 107, 19 (1995) and by Hayes in SMRM 5t annual meeting, Montreal, Book of Abstracts, 39 (1986)). These designs were applied to structures that surrounded only the head of a patient and, by their nature, prevent access to the top of the head. The limited access also makes these structures problematic for whole-body imaging as they substantially reduce access from one end of the magnet.
It is an aim of this invention to provide coil structures that generate desired RF fields within certain specific, and asymmetric portions of the overall coil structure, preferably without substantially limiting access from one end of the structure.
20 Asymmetric radio frequency coils can be used in conventional MR systems or in the newly developed asymmetric magnets of U.S. Patent No. 6,140,900.
SUMMARY OF THE INVENTION In a broad form, the invention provides a coil structure for a magnetic resonance device having a cylindrical space with open ends, the coil structure being adapted to generate a desired RF field within a specified portion of the cylindrical space. This portion is asymmetrically located relative to the mid-length point of the longitudinal axis of the cylindrical space.
A method for manufacturing a radio frequency coil structure for a MR device having a cylindrical space, preferably with open ends, may comprise the steps of selecting a target region over which a transverse RF magnetic field of a predetermined frequency is to be applied by the coil structure, the target region being 4 preferably asymmetrically located relative to the mid-length point of the longitudinal axis of the cylindrical space, calculating current density at the surface of the cylindrical space required to generate the target field at the predetermined frequency, synthesizing a design for the coil structure from the calculated current density in accordance with one of the methods discussed below, and forming a coil structure according to the synthesized design.
Preferably, the method for calculating the current density uses a time harmonic method that accounts for the frequency of operation of the RF coil structure and makes use of a complex current density.
The RF coils of the invention can be used as transmitter coils, receiver coils, or both transmitter and receiver coils. As discussed above, when the coil serves a transmitting function, it will normally be combined with a shield to reduce magnetic interference with external components of the magnetic resonance imaging system.
To avoid redundancy, the following summary of suitable methods for designing apparatus for transmitting a radio frequency field is in terms of a RF coil system which includes a main coil (corresponding to the "first complex current density") and a shielding coil (corresponding to the "second complex current density"), it being understood that these methods can be practiced with just a main coil.
o 20 For "mild" coil length to wavelength conditions, conditions in which the coil length is less than about one-fifth of the operating wavelength, a method for S-designing apparatus for transmitting a radio frequency field a field having a frequency of at least 20 Megahertz, preferably at least 80 Megahertz), receiving a 2 magnetic resonance signal, or both transmitting a radio frequency field and receiving a magnetic resonance signal, comprises: defining a target region a spherical or ellipsoidal region preferably Shaving a volume of at least 30 x10 3 cm 3 and asymmetrically located relative to the midpoint of the RF coil, the ratio of the distance between the midpoint of the •o.
target region and one end of the coil to the length of the coil is less than or equal to 0.4) in which the radial magnetic component Bx, By, or combinations of B, and By) of the radio frequency field is to have desired values substantially uniform magnitudes), said target region surrounding a longitudinal axis, the common longitudinal axis of the magnetic resonance imaging system and the RF coil; specifying desired values for said radial magnetic component of the radio frequency field at a preselected set of points within the target region a set of points distributed throughout the target region or a set of points distributed at the outer boundary of the target region); defining a target surface the shield coil surface) external to the apparatus on which the magnetic component of the radio frequency field is to have a desired value of zero at a preselected set of points on said target surface; determining a first complex current density function, having real and imaginary parts, on a first specified cylindrical surface the main coil surface) and a second complex current density, having real and imaginary parts, on a second specified cylindrical surface the shield coil surface), the radius of the second specified cylindrical surface being greater than the radius of the first specified cylindrical surface the radius of the second surface can be 10% greater than the radius of the first surface) by: defining each of the complex current density functions as a sum of a series of basis functions triangular and/or pulse or sinusoidal and/or cosinusoidal functions) multiplied by complex amplitude coefficients having real and imaginary parts; and determining values for the complex amplitude coefficients using an 20 iterative minimization technique a linear steepest descent or a conjugate gradient descent technique) applied to a first residue vector obtained by taking the difference between calculated field values obtained using the complex amplitude coefficients at the set of preselected points in the target region and the desired values at those points and a second residue vector equal to calculated field values obtained using the complex 25 amplitude coefficients at the preselected set of points on the target surface; and converting said first and second complex current density functions into sets of capacitive elements (as understood by persons skilled in the art, such capacitive elements will in general have some inductive and resistive properties) and sets of inductive elements (as understood by persons skilled in the art, such inductive elements will in general have some capacitive and resistive properties) located on the specified cylindrical surfaces by: 6 converting each of the first and second complex current density functions into a curl-free component Jcurw-fre and a divergence-free component Jdiv-free using the relationships: J,crI-fre V and J~dv-fee =VxS, where W and S are functions obtained from the respective first and second complex current density functions through the equations:
V
2 V. J
-V
2 S VxJ, and
-V
2 (n.S)=n.VxJ, where n is a vector normal to the respective first and second specified cylindrical surfaces and J is the respective first and second complex current density functions; (ii) calculating locations on the respective first and second cylindrical surfaces for the respective sets of capacitive elements by contouring the respective y functions; and (iii) calculating locations on the respective first and second cylindrical surfaces for the respective sets of inductive elements by contouring the respective functions neS (W and n*S are referred to herein and function as "streaming functions").
For "non-mild" coil length to wavelength conditions, conditions in which the coil length can be greater than about one-fifth of the operating wavelength, a method for designing apparatus for transmitting a radio frequency field or both transmitting a radio frequency field and receiving a magnetic resonance signal, :i comprises: defining a target region in which the radial magnetic component of the radio frequency field is to have desired values, said target region surrounding a longitudinal axis, the common longitudinal axis of the magnetic resonance imaging system and the RF coil; specifying desired values for said radial magnetic component of the radio frequency field at a preselected set of points within the target region; defining a target surface external to the apparatus on which the magnetic component of the radio frequency field is to have a desired value of zero; determining a first complex current density function, having real and imaginary parts, on a first specified cylindrical surface and a second complex current density, having real and imaginary parts, on a second specified cylindrical surface, the radius of the second specified cylindrical surface being greater than the radius of the first specified cylindrical surface: defining each of the complex current density functions as a sum of a series of basis functions sinusoidal and/or cosinusoidal functions) multiplied by complex amplitude coefficients having real and imaginary parts; and (ii) determining values for the complex amplitude coefficients by simultaneously solving matrix equations of the form: L BC (Eq.I)
BS
where A c AS, A C and A s are transformation matrices between current density space and magnetic field space whose components are based on time harmonic Green's functions, aC and aS are vectors of the unknown complex amplitude coefficients for the 2 first and second complex current density functions, respectively, BC is a vector of the go• 20 desired values for the radial magnetic field specified in step and B s is a vector S°whose values are zero, said equations being solved by: transforming the equations into functionals that can be solved using a preselected regularization technique, and solving the functionals using said regularization technique to obtain values for the complex amplitude coefficients; and converting said first and second complex current density functions into sets of capacitive elements and sets of inductive elements located on the specified cylindrical surfaces (preferably, this step is performed using the methods of step of S the first method aspect of the invention; also the methods of said step can be used S 30 independent of either the first or second method aspects of the invention to convert a complex current density function to sets of capacitive and inductive elements, to a manufacturable coil structure).
8 The regularization functional is preferably chosen so as to minimize the integral of the dot product of the first complex current density function with itself over the first specified cylindrical surface and to minimize the integral of the dot product of the second complex current density function with itself over the second specified cylindrical surface. Other regularization functionals that can be used include: the second derivative of the complex current density functions, and other functionals besides the dot product that are proportional to the square of the complex current density functions.
The complex amplitude coefficients are preferably chosen so that the first and second complex current density functions each has zero divergence.
The parenthetical statements set forth in connection with the summary of the first method also apply to the second method except where indicated. Both methods also preferably include the additional step(s) of displaying the locations of the sets of inductive elements on the first and second specified cylindrical surfaces and/or producing physical embodiments of those sets of elements.
In accordance with one aspect, the invention provides apparatus for use in a magnetic resonance system for transmitting a radio frequency field a field having a frequency greater than, for example, 20 megahertz, and preferably greater than megahertz), receiving a magnetic resonance signal, or transmitting a radio frequency io. 20 field and receiving a magnetic resonance signal, said apparatus and said magnetic o resonance imaging system having a common longitudinal axis, said apparatus comprising: a support member a tube composed of fiberglass, TEFLON, or other materials, which preferably can be used as a substrate for aphotolithography 25 procedure for printing conductive patterns and which does not substantially absorb RF energy) which defines a bore (preferably a cylindrical bore) having first and second open ends which are spaced from one another along the longitudinal axis by a distance L (the open ends are preferably both sized to receive a patient's body part which is to be imaged, the head, the upper torso, the lower torso, a limb, etc.); and a plurality of inductive elements copper or other metallic tracks or tubes) and a plurality of capacitive elements distributed or lumped elements) associated with the support member mounted on and/or mounted in and/or mounted to the support member); wherein if used as a transmitter, the apparatus has the following characteristics: the apparatus produces a radio frequency field which has a radial magnetic component which has a spatially-varying peak magnitude whose average value is Ar.avg; (ii) the apparatus has a homogenous volume within the bore over which the spatially-varying peak magnitude of the radio frequency field has a maximum deviation from Ar-avg which is less than or equal to 15% (preferably less than or equal to (iii) the homogeneous volume defines a midpoint M which is on the longitudinal axis, is closer to the first end than to the second end, and is spaced from the first end by a distance D such that the ratio D/L is less than or equal to 0.4 (preferably less than or equal to 0.25); and (iv) at least one of said inductive elements 1, 2, 3, or all) comprises a discrete conductor a copper or other metallic track or tube) which follows a sinuous path such that during use of the apparatus current flows through a first part of the conductor in a first direction that has both longitudinal and azimuthal components and through a second part of the conductor in a second direction that has both longitudinal and azimuthal components, said first and second directions being different.
example of such a sinuous path is shown in Figure 15 discussed below.
20 The terminology "if used as a transmitter" is used in defining both transmitting and receiving RF coils since by reciprocity, coils that transmit uniform radio frequency fields, receive radio frequency fields with a uniform weighting function.
In certain preferred embodiments, the homogenous volume is at least 30 x 10 3 cm 3 and/or L is at least 1.0 meter.
In other preferred embodiments, the first and second open ends have transverse cross-sectional areas A, and A 2 respectively, which satisfy the relationship:
IAI-A
2 1/AI 0.1.
*where A, is preferably at least 2 x 10 3 cm 2 •In connection with all of the aspects of the invention, the lengths of coils and of cylindrical spaces associated therewith are defined in terms of the inductive elements making up the coil. In particular, for a horizontally oriented coil, the length constitutes the distance along the longitudinal axis from the leftmost edge of the leftmost inductive element to the rightmost edge of the rightmost inductive element, corresponding definitions applying to other orientations of the coil.
BRIEF DESCRIPTION OF THE DRAWINGS The invention will be described by way of examples with reference to the drawings in which: FIGURE 1 illustrates a general layout of a cylindrical radio frequency coil system of length L and radius a. Changing the value of D changes the degree of symmetry of the 'target region' relative to the coil structure; FIGURE 2 shows a current distribution on the surface of a cylinder calculated using a quasi-static formulation, where the vertical axis is the z-axis; FIGURE 3 shows plots of J0 and J, along the length of a cylinder in the 0=0 and 0=7/2 planes respectively for a quasi-static formulation; FIGURES 4A and 4B show plots of J0 and J, along the length of a cylinder in the 4=0 and planes respectively for the real solutions of Table 4 (190MHz). In FIGURE 4A, z,,=20mm, while in FIGURE 4B, FIGURE 5 shows the current distribution on the surface of a cylinder calculated from the real solution of Table 4 (z,,=20mm; 190MHz); FIGURE 6 shows the current distribution on the surface ofa cylinder calculated from the real solution of Table 4 (zos=50mm; 190MHz); 20 FIGURE 7 shows a full conductor model of the contours of the symmetric current density distribution shown in FIG. 3.
FIGURE 8 shows a current density distribution for a 320mm length coil with radius 100mm having a homogeneous volume 50mm from the center of the coil *towards the coil's lower end.
FIGURE 9 shows contour lines for neS for the current density of FIG. 8.
FIGURE 10 shows contour lines for y for the current density of FIG. 8.
FIGURE 11 shows the contour lines for the stream function noS calculated using the second method.
FIGURE 12 shows a normalized magnetic field for a representative asymmetric coil in the transverse plane determined using the second method. The upper curve shows the distribution in the x direction and the lower curve shows the distribution in the y direction.
11 FIGURE 13 shows a transverse magnetic field along the z-axis for a representative asymmetric coil determined using the second method.
FIGURE 14 is a circuit-type layout for a representative asymmetric coil.
FIGURE 15 is a 3D perspective of an inductive geometry for a representative asymmetric coil. A partial shield is included for illustrative purposes.
FIGURE 16 is a flow chart illustrating the first method.
FIGURE 17 is a flow chart illustrating the second method.
FIGURE 18 is a flow chart illustrating a preferred method for calculating a coil geometry from a current density.
The foregoing drawings, which are incorporated in and constitute part of the specification, illustrate the preferred embodiments of the invention, and together with the description, serve to explain the principles of the invention. It is to be understood, of course, that both the drawings and the description are explanatory only and are not restrictive of the invention.
DETAILED DESCRIPTION OF EMBODIMENTS OF THE INVENTION As discussed above, the present invention relates to RF coils having prescribed properties and to methods for designing these and other types of RF coils. Figures 16- 18 illustrate the overall numerical procedures of the invention with reference to the 20 various equations presented below.
S: The methods described below are preferably practiced on a digital computer system configured by suitable programming to perform the various computational S: steps. The programming can be done in various programming languages known in the art. A preferred programming language is the C language which is particularly wellii i 25 suited to performing scientific calculations. Other languages which can be used include FORTRAN, BASIC, PASCAL, C, and the like. The program can be embodied as an article of manufacture comprising a computer usable medium, such as a magnetic disc, an optical disc, or the like, upon which the program is encoded.
The computer system can comprise a general purpose scientific computer and 30 its associated peripherals, such as the computers and peripherals currently being manufactured by DIGITAL EQUIPMENT CORPORATION, IBM, HEWLETT- PACKARD, SUN MICROSYSTEMS, SGI or the like. For example, the numerical procedures of the invention can be implemented in C-code and performed on a work 12 station. The system should include means for inputting data and means for outputting the results of the RF coil design both in electronic and visual form. The output can also be stored on a disk drive, tape drive, or the like for further analysis and/or subsequent display.
QUASI-STATIC ANALYSIS In outline, the object of the design process is to produce structures that generate target RF transverse fields (that is, Bx or By or combinations thereof) over a defined volume inside a cylinder. The value of Bx (for example) is specified at a number of chosen points inside this region and in one embodiment, inverse methods are used to compute the corresponding current distribution.
The geometry of the coil is shown in Figure 1. Surface current densities on the cylinder are defined as: J(r) J z) J, z)2 for p= a and -L/2 z L/2 (1) where r p zi L and a are the length and radius of the cylinder.
As is known in the art, in one formulation for the case of a symmetric coil structure, J(r) can be written as a summation of M sinusoidal basis functions: cos0 J sin(kz) m=l S 20 z)=sin J-C cos(kmz) ee 2) m (3) where Jcm are the current coefficients and I However, this expression is only applicable to problems where the magnetic 'field is symmetrical about z=0. Magnetic fields are expressed in terms of magnetic vector potentials defined as the convolution of surface current density with the corresponding Green's function in the spatial domain: 13 L/2 2r I I G(r d'dz' for pa (4) -LI2 0 where the dashed elements indicate source points and the undashed the field points.
The corresponding generalized Green's function in cartesian coordinates is: e-Jelr-r G(r I eJkr-r' ir r' r-r' where r and r' are the position vectors of field and source points respectively; ko is the wave-number given by ko=27n/o; and Xo is the free-space wavelength of the propagating field.
We consider first an initial quasi-static simplification in which the wavelength is assumed to be very long, or the frequency near zero. Equation then approaches: 1 G(r I =r r r r' (6) In cylindrical coordinates one solution of this is: G(r e p dkcos[k(z-z')]l (kp )Kp(kp') 0 15 (7) where r are again defined as the position vectors of field (source) respectively. Ip and Kp are the modified Bessel functions ofp h order. This solution is for the case S****where observation points are internal to the cylinder, p' p.
The vector potentials may then be expressed as: 1.I2 2l 20 a /2 0 14 L12 2;r Ar f f JG(r I cos( 0')d0'dz' 4r-L12 0 (9) L12 2rf A,(r Lopa f fJG(r I d0'dz' 41r -L12 0 Substituting equation into the above equations and carrying out the integration over the cylinder's surface, we obtain the expression for in the form of Fourier integral (series) or in terms of spectral variable k The relevant magnetic fields are calculated from the curl of the vector potentials to be: BP(r) Po a Z P fk Jcos kz JO k)IP (k p) K'P(ka)dk 7r 0 (1 BO(r) jia peP cszJ(~)~k)'(ad (12) Writing the currents as summation of and and using the following relationship: BP cos 0- BOsin 0 (13) we obtain B 9~a L Jcos/cz EJ. (ka)[kp cos' 01 (kp) cos 2q0 1 (kp)]dk 2 rp 0 m=1 (14) where sin(k k)L 2 sin(k k)L 2 (k k,)L2 (k +km)L/2 The result described in equations (14) and (15) is known in the art as a quasistatic solution, being only valid in applications where the free space wavelength is much longer than the length of the cylinder generating the RF fields (see for example, H. Fujita, L. S. Petropoulos, M. A. Morich, S. M. Shvartsman, and R. W. Brown 'A hybrid inverse approach applied to the design of lumped-element RF coils,' IEEE Trans. Biomedical Engineering, vol-46, pp. 353-361, Mar. 1999).
A target region is then defined as the volume within which a homogenous Bx is required. In particular, a target field value is defined as the fixed magnitude of the desired homogenous field Bx. The target region is, for example, specified as the volume inside a sphere of radius R centered at z=0. As an example, a target specification could be that R is set to be R 0.7a and max(Bx)-min(Bx) is specified at less than If we set Bx(rl)=B,(r 2 )=Bx(r 3 Bx(rN)= Bx~ where rl, r 2 r 3 rN are the target field points chosen inside the target region, a set of N linear equations is obtained from equation (14).
In particular, the following matrix equation is obtained: (16) where is a matrix containing the values of Bx at the defined observation points due Si to each of the current components; [Jo. is a column vector of the current coefficients Jci, Jc 2 Jc 3 Jcm and is a unity column vector.
If the number of linear equations equals the number of unknowns then 25 the numerical solution for [Jc] can be easily obtained by multiplying both sides of (16) by As an example of the quasi-static method, a coil of dimension: a=100 mm, L=320 mm was designed. The specified target field in spherical coordinates was: The current density coefficients up to M=3 were then calculated from the quasistatic expressions in equations (14) to be: Jcl Jc2 Jc3 43.834274291992190 -48.744728088378910 36.222351074218750 The corresponding resultant current distribution on the surface of the cylinder is shown in Figure 2 and illustrates the current paths required for the coil. In the graph of Figure 3, the current components J0 and J. are plotted along the z-direction for <=0 and 4=it/2 respectively. Contour plots (not shown) demonstrate that B, is homogeneous within the target region.
TIME HARMONIC ANALYSIS Instead of using the quasi-static approach, a time harmonic, full wave method is used to calculate the required current density. This is essential in cases where the operating wavelength and the dimension of the cylinder are of the same order, as the prior art quasi-static method described in the above section is no longer valid.
However, the mathematical derivation of a time-harmonic expression for Bx in the spectral domain or an equivalent of (14) is extremely difficult. A more direct approach is to define Bx directly in terms of the vector potentials of the generalized Green's function of equation rather than its spectral equivalent given in equation aA, (17) where A, are Ay are the vector potentials calculated from equation (Note that numerically, the calculation of equation (17) is carried out using a Gaussian-Quadrature integration routine over the prescribed 2-D cylindrical surface integral.) In the time-harmonic case, the wavelength or frequency of operation is taken into account and the wavenumber (ko) is non-zero. This leads in many instances to the 17 current density and its constituent coefficients being complex (that is, consisting of both real and imaginary components). The calculated current densities for the generation of specific transverse RF fields are such complex quantities.
To illustrate the use of complex current densities, current densities were calculated on the same cylinder as the quasi-static example discussed above at three different frequencies: 190 MHz, 300 MHz and 500 MHz. The same target field was used, the homogeneous volume was symmetrically located. Equation (16) was modified to include complex current densities and a complex version of the matrix [T] was obtained using equation The solution for then follows the same procedures for the quasi-static example except that complex quantities are now included. The calculated current coefficients are shown in the following table: Frequency (MHz) Jc Jc 2 Jc 3 190 38.6895+jl.9310 -42.7225-j2.1026 35.0854+j 1.6609 300 33.4992+j5.4433 -36.6628-j5.7392 35.1309+j4.9677 500 23.8764+j 13.0168 -25.3327-jl2.4308 39.07675+j 14.0778 Table 1. Computed current coefficients at different frequencies.
From the above table, it is clear that as the frequency increases, the imaginary part of the current densities becomes larger.
BASIS FUNCTIONS FOR ASYMMETRIC SYSTEMS In a preferred embodiment of the invention, open-ended asymmetric structures are designed for a time harmonic system. For the design of asymmetric coil structures, surface current densities (Amp/meter) on a cylinder are defined in the same manner as for symmetric structures: J(r) J J, z) for p=a and -L/2 z L/2 (18) The relationship between Jz and J under mild coil length/wavelength ratios, may be approximated as: J z) -a d z do (19) Or alternately Jz and J, may be decoupled and treated individually for the purposes of optimization. Currents on the coil are descretized along z and are represented by subdomain basis functions.
For asymmetric structures a new set of basis functions must be specified, and in one embodiment J, and J are written as a summation of triangular Tz() and pulse functions respectively as follows:
M
J z) cos0 JcmT, (z z) m=l a M J, sin -J{n lz n=1 (21) where the Jcm are the current coefficients and elsewhere (22) 1 z-z±1<1,/2 T(z-z m 0 elsewhere (23) 1 L and z ml. z (m and 2 2 2 19 Magnetic fields are expressed in terms of magnetic vector potentials defined as the convolution of surface current density with the corresponding Green's function in the spatial domain: LI12 2Z a f f G(r d 'dz' forp=a 4-L/2 0 (24) As previously, the corresponding generalized Green's function in cartesian coordinates is: e-Akol'-'l G(r r-r Ir-r' where r and r' are once again the position vectors of field and source points respectively; ko is the wave-number given by ko=2n/,o0; and X 0 and is the free-space wavelength of the propagating field.
The transverse field is written as: Bx aA a y az 15 (26) S where Az are A, are the vector potentials calculated by equation (24).
The homogeneous field is preferably asymmetrical about the z=0 plane. The target region is defined as a spherical volume of radius R centered at z=os,. A first i:i method of obtaining a set of current coefficients that produces the required 20 homogeneity is as follows.
Let (rl, r 2 r 3 rN) be the N field positions chosen inside the target region and Bx(r,) be the magnetic field at these points. The field at each target point is expressed in terms of M number of current coefficients as follows:
M
jC,,J m=l (27) where represents the field contribution at due to the basis function at z z, The residue is defined to be the difference between Bx(r,) and the specified target field Bx,, i.e., R,(J J, 2 JN) B(r, n=1..N (28) A functionf is defined as:
N
f(J,,J 2 JN)= R, (JI,J2 JN) (29) n=l In order to suppress large current variations in the solution, it is sometimes useful to also include the change of current amplitude as a constrained condition. In this
N-I
case a term w m, is added to equation where m, is the difference between the 1=1 adjacent current elements given by J Jil w is the weight used for bringing m, to the same order of magnitude as R,.
The task is to find a set of J 2 J such thatfis minimum, or fmin min{f(JI, J 2 J There are various ways in which this function may be minimized. In one embodiment, a linear steepest descent method is used for each variable. An advantage of using this minimization technique is that it allows either a o real or a complex solution to be readily obtained.
EXAMPLES OF THE FIRST METHOD Examples of the invention, designed using the time harmonic, full wave method of the first method are now detailed.
i i A coil of dimension: a=100 mm, L=320 mm was designed at different frequencies: 0 MHz, 190 MHz and 300 MHz. A set of target points was defined in spherical coordinates with respect to zos, where in terms of the parameters of Figure 1, :*Izosl IDI L/2. The radial positions of these points were r,=53.3mm; their angular 25 positions are given in the following table: N 0 1 00 00 2 450 00 3 900 00 4 450 300 00 450 6 450 450 7 450 600 8 00 900 9 900 900 450 1200 11 450 1350 12 450 1500 13 00 1800 Table 2. Angular positions of target points in spherical coordinates centered at zos.
The specified target field was B,=23.5 [iT. The requirement for homogeneity was specified as: max(Bx)-min(Bx)< 5% inside a spherical region of radius R z 0.7a for mm and R 0.55a for zos =50 mm. Eight current coefficients were used for each computation. Tables 3 to 5 show the calculated current coefficients for the different zos values at the three operating frequencies. For the non-static cases, the real (Re) and imaginary (Im) components of the complex currents are shown separately.
10 These coefficients are for use with the triangular and pulse functions discussed above.
Zos Jc Jc2 Jc3 J4 220.68 -188.62 -7.60 -42.38 1763.41 -959.25 -163.87 -52.156 Table 3. Computed current coefficients for quasi-static case.
Zos(mm) Re/Im J Jc2 Jc3 Jc4 Re 226.65 -195.19 17.93 -42.11 Im -159.73 +71.66 -5.96 -4.30 Re 815.32 -280.51 -68.98 0.22 Im -505.05 167.50 5.94 -2.68 Zos(mm) Re/Im J Jc6 Jc7 Jc8 Re -24.42 -33.19 -20.73 -55.24 Im +0.65 -1.77 -0.79 -1.81 Re -33.27 -30.49 -9.94 -63.67 Im 0.21 0.98 0.47 0.44 Table 4. Computed complex current coefficients at 190 MHz.
Zos Re/Im Jci Jc2 Jc3 Jc4 Re 194.73 -14.29 -22.91 Im 117.21 -21.10 -1.15 Re 634.66 -48.97 23.80 Im 0.92 3.09 9.78 4.87 Zos (mm) Re/Im Jcs J6 Jc7 J8 20 Re -23.01 -29.57 -10.18 -51.95 Im -1.32 -4.98 0.61 -10.29 Re -36.71 -27.10 -8.53 -53.72 Im 4.19 3.67 2.50 2.32 Table 5. Computed complex current coefficients at 300 MHz.
The transverse magnetic field (Bx) inside the cylinder was examined for each of the above cases. Computed results show that the resultant Bx values meet the required homogeneity for all the above cases and they show similar characteristics for a given Zos. For the 190 MHz case, the real components of the current densities J and Jz are plotted along the z-direction in the j=0 and =n/2 planes on the graphs of Figure 4A and 4B which correspond to Zos=20mm and Zos=50mm respectively. In both cases, J.
23 and J, along z are of an oscillatory nature. It is seen that both components experience a surge in magnitude near the positive end of the coil the end closest to the end where the target region is located).
Figures 5 and 6 show current distributions on the surface of a cylinder from the real solutions of Table 4 for zos equal to 20mm and 50mm, respectively.
Plots of the B, field created by the currents of Figure 4A and Figure 4B over the volume p 80 mm, -160 mm z 5 160 mm (not shown) revealed a homogeneous region in both cases that was asymmetrical about the center of the coil and satisfied the criteria for homogeneity. It was also found that there was a larger volume of homogeneity in the Zo,=20mm case than in the Zos=50mm case.
EXAMPLES OF THE SECOND METHOD As a representative illustration of the second method, the following sets of basis functions are used for the current density components on the surface of a cylinder with radius po and length L J mn ,Cos nC cos kz+ (2p m) a, q=0 p=0 n=1 m=p+l e 0 0 1v cmsin +in sin kmz(2pm) a. (31) q=0 p=0 n=1 m=p+lwhere Skm (32) km L Sf and the complex coefficients amnp and cmn, are to be calculated.
25 This form of the basis functions was selected based on the following rationale.
If the specified field has only one vector component that does not vary with 0, then it can be reasoned that the term in the basis functions is unnecessary as the axis can 2 be set in any direction in the xy plane. However for simplicity, 4 0 is set to coincide with the x-axis and the term will be necessary for any specified field that has a 2 y-component.
For specified fields in MRI applications, it is expected that the current density J is anti-symmetrical at 180°, so coefficients of even harmonics of 0 will be zero and hence, n Both components of J have a z dependence described with k,z where p=0,l and m is an integer series commencing at p. The index 2 p=0 describes basis functions that have J. 0 at z ±L 2 whereas the index p=l describes basis functions with a J, 0 at z ±L 2.
If the specified field is symmetric with respect to the z 0 plane, then all coefficients of even m terms will equal zero and hence m can be assigned m However, all terms of m are required if the specified region is asymmetric with respect to the z 0 plane as in the example coil design presented herein.
The next step in the method is to calculate the coefficients of the basis functions in equations (30) and (31) such that a specified homogeneous field is generated.
Because the time harmonic Green's functions are of a complex variable, the coefficients of the current density J are complex as well. If the specified field is that of a circularly polarized field, then the real and imaginary components would be the same albeit spatially separated by 90*. So if J(4, z) produces Ka x in the DSV, J(4 7 2, z) 20 would produce Kay where K is just a constant. The real part of the current that would produce a circularly polarized field is then: Jin real{J((, imag{J(( (33) 25 and this is the same current that would produce a linearly polarized field. Hence only B~Bi is needed to be specified in the DSV to produce J( z) and (33) can be used to obtain the final current density with only a real component.
Transmitting RF coils must have a shield to prevent eddy currents in external conductors and to provide a suitable RF ground. A shield may be approximated by specifying B 0 on a cylinder with a radius slightly larger than the shield radius. A current is approximated on the shield by the same set of basis functions as those for the inner coil except that the shield current now has both div-free and curl-free components (see below).
If and (c c denote column vectors of the J, and J: coefficients of the coil current (see equation (as) and (c
S
denote column vectors of the J, and J_ coefficients of the shield current (see equation and denotes a matrix, then the field in the DSV due to the coil currents is found in matrix form, as:
=(B
c (34) where the summation of the vectors(Bc,) and results in the target field vector(Bx, )within the DSV. Similar matrix equations result when applying B. 0, B, 0 and B. 0 as the condition at specified points near the outside of the surface of the shield cylinder: .".The matrix equations are usually not square because the number of points where the magnetic field is specified does not usually equal the number of unknown Scoefficients. The rank of the matrices are usually less than the number of unknown S coefficients as well, and so a regularization method must be used for a solution.
S: 25 It should be noted that equations and (36) can be written in more compact form as Eq. above.
REGULARIZATION
The regularization method chosen as a non-limiting example was to minimize some functional in terms of the current density and impose the extra conditions onto the matrix equations such as (34) to Numerous functionals are available and it is 26 convenient to choose a functional such that the resulting current density is easier to implement. The functional: min J Jds (37) So results in a current density similar to one of minimum power and hence is a useful choice. When the current density equations (30) and (31) are substituted into the above equation (37) and differentiated to find the coefficients for a minimum: (J.J)ds=0 (38)
S
0 ijk f a (J.J)ds=O (39) so Uv the resulting matrix equation is just some constant multiplied by the identity matrix.
This is because the basis functions are orthogonal. For a current density with no divergence, the wire implementation depends on the stream-function neS, where Z dp, Mpcos sin kz+(2p q=0 p=0 n=1 m=p+l e where n is normal to the cylindrical surface po(pokmamn nc d, PoPokanpq mnpq) (41) 5: d VVn*2S)22ds (41) n po Hence minimizing the variation in S over the surface would minimize the variation in its contours that ultimately determine the wire positions: s 25 min J(VV(n S)ds (41A)
SO
where V V(n S) gives a quantitative measure of the stream-function's variation.
27 In this case, because the basis functions on the coil form a set with zero divergence, the dm, coefficients of S reduce to: dnpq= pq (42) k, Applying the condition (42) results in a matrix with diagonal terms: nL 2 (43) being the regularization matrix, and being the identity matrix.
MATRIX SOLUTION An iterative method was used as a non-limiting example to calculate the coil coefficients and shield coefficients separately so that the error in each could be adjusted more exactly by adjusting two scalar penalty values and 2 Firstly, the initial coil currents were calculated using the expression: (44) where is the combination of and of (34) using the relationship of equation 20 (42) and is the regularization matrix calculated by applying (37) to the coil currents.
S" The field very close to the shield due solely to the coil currents is: where is the combination of and of(36). Shield currents are calculated to this eld: *o negate this field: ([MP]T ?X 2
[R
2
[AS"]
T
(bI) (46) where (S) T ((as )T of (36) and [R 2 1 is the regularization matrix calculated by applying (37) to the shield currents. The field in the DSV due to the shield currents is: (bDSVIS) (47) where of The field in the DSV due to the coil currents is: (bI,)svlc) (48) The error in the DSV is then: -(bI)svlc) -(bDsvls) (49) The new coil coefficient estimate is: (ac 1) (as) ]T X, [R (b 20 The calculation then loops back to (45) and is terminated when the difference in the error between two consecutive iterations falls below a certain predefined limit 0.01). The error SDS;V here is defined as: _I (b) ED)SV (b, (51) and the shield error F- is defined as the error in the field produced outside the shield: C= Jsc+ (52) The penalty scalars X 1 and ,X 2 affect the current density solution as well as determining the condition of the matrices to be inverted in equations (46) and Typically, the matrix problem becomes unstable (that is, where the matrix is tending towards being singular) below a certain value of X.
EXAMPLES IN WHICH CURRENT DENSITIES ARE CONVERTED TO DISCRETE COIL STRUCTURES By whatever method the complex current density functions are obtained, the next step in the process of the overall design method is to synthesize the coil structure form the current density.
In a preferred embodiment, a stream function method is employed to calculate the positions and dimensions of conductors such that the stream function of the current distribution in the conductors closely approximates the stream function of the calculated current density.
QUASI-STATIC CASE It can be seen from equations and that the divergence of the quasi-static current density is zero: V.J=0 (53) In true electromagnetic systems, the divergence is not zero but: 20 V.J =-jwcop (54) a.
11 where p is the charge density. However, for this analysis, the charge distribution is assumed to be negligible, which is true for static and low frequency systems.
a.
Because of zero divergence, a function S can be calculated that represents the magnetic flux produced from the current density as: 25 VxS=J where from equations and S =-ncos Ccoskmz (56) m=1 k., where n is the radial unit vector and is normal to the cylindrical surface. Note that ifJ is a surface current density, S will only have a component perpendicular to the surface containing J, and hence only the magnitude of S is needed. The magnitude of S is a function in two dimensions.
The current distribution in the conductors should be such that it produces the same function S as that produced by the calculated current density J. This is done by finding the contours of n*S and generating appropriate currents in each of the wires placed in the positions of the contours. If the contour 'heights' are chosen at equal increments between minimum and maximum, the currents in each of the wires following the contours should be equal. If the contour 'heights' are not equally spaced, the currents in the wires have to generated accordingly in proportion to the contour spacing.
By way of an example, consider equation (56) with six real coefficients cm=47.3518, 212.791, 314.615, 687.005, -3946.81 and 5365.36.
The values of the stream function dictates the shape and position of the contours and hence, is decided iteratively by investigating the electromagnetic simulations each contour configuration produces. For example, the contours may be taken at equal fractions of the stream-function's maximum value, resulting in equal current in all the wires. Alternatively, if the maximum is normalized to one, contours could be taken at ±0.697, ±0.394 and ±0.091. Because the interval (0.091) between zero and neS=0.091 does not equal the interval (0.303) between neS=0.091 and neS=0.394, the current in the lowest 0.091 contour will have to be 0.3 0.091/0.303) of the current flowing in the 0.394 contour and 0.697 contour (these two contours should have equal current).
20 The next task is to derive a workable model from these contours.
Design decisions in this process include: 1. A magnetic symmetry plane at 4 =0 is used for the stream function in the range -L/2<z<L/2 and <R.
2. Each contour is implemented as a conductor. Two cases were examined: one where the conductors are implemented as circular wires (tubes) of specified diameter and the other where the conductors are implemented as metal strips of specified width.
The diameter or width ofa conductor determines its self-inductance and hence its complex impedance. It is possible to adjust the current levels in the conductors by adjusting the diameter or width of the conductors accordingly.
S 30 3. The contours are connected depending on the sign of the stream function. If two contours are in close proximity and the stream function value has the same sign for both contours, then the contours are connected in parallel. If the contours have a stream function value of differing sign, the contours are connected in series.
31 4. Capacitors are added at the connections of the contours connected in series and also at the mid-point of contours that are not connected to other conductors.
However, there is a freedom to add capacitors wherever the designer considers suitable. The values of the capacitors are determined from an iterative approach using repeated frequency sweeps in the device simulation.
The source points are at the extreme ends of the contour, that is, either at both ends or at one end: either at z=L/2 and z=-L/2 or only at z=L/2.
A full implementation of this method results in the conductor configuration shown in Figure 7 for the current density defined by the six real coefficients set forth above, cm=47.3518, 212.791,314.615, 687.005, -3946.81 and 5365.36.
In a preferred embodiment, asymmetric structures are synthesized.
When triangular basis functions are employed to approximate the current density, only the current-density values are provided at points within the domain of each triangular function. Thus the problem is to find the appropriate contours of a twodimensional grid-like data set that will best approximate the current density that the data set is meant to represent. In the quasi-static problem, the current density has zero divergence because of the nature of the basis functions and thus the electric field at the surface is parallel with the current density. In the non-quasi-static case however, the electric field can have a component perpendicular to the current density and so an 20 approach has to include measures for the approximation of the charge density. In accordance with the preferred embodiments of the invention, the curl-free and divergence-free components of the current density are found separately and each is approximated in turn.
CURRENT-DENSITY COMPONENTS The current density J can be split into its curl-free and divergence-free i components by writing: J =VxS+Vy (57) From Laplace's law, it follows that 30 V y/=V.J (58) and the curl-free component of J is: J,1-free V(9 (59) 32 The divergence-free component can be found by applying curl to both sides of equation (57) to give: VxVxS=VxJ If the divergence of S is made zero: V-S=0 then:
-V
2 S= VxJ The divergence-free component is: Jdv-, =VxS (61) Since Jdv-free is only on the surface, and to be free of singularities in the differentiation normal to a surface, it is best if nxS=0 Hence, (60) also becomes a Poisson equation over the surface of J: -V2 (n nVxJ (62) The task then is to solve the two Poisson equations (58) and (62) for y and n*S and thereby using the contours of W and the contours of neS to approximate the entire current density J.
Note that Jdiv-free is the current density that directly influences the magnetic field tangential to the surface of the conductor. Jcurl-free is a measure of the quantity of 20 rate of change of charge density and this quantity directly influences the electric field normal to the surface. In other words, Jdiv-free can be thought of as the inductive component of the circuit and Jcurl-free as the capacitive component of the circuit.
Once S and y are calculated, the contours and configuration of the coils can be designed. The contours of n*S dictate the current configuration. Once n*S is found, its contours and hence the coil pattern can be calculated in a similar manner to the quasistatic case.
As examples of the above procedures, Figures 9 and 10 show the results of determining inductive and capacitive elements for complex current density functions obtained using the first method aspect of the invention and shown in Figure 8. The coil had a homogeneous region 50mm from the center towards the lower end; the length of the coil was 320mm and its radius was 100mm. The operating frequency was 190MHz.
33 A contour diagram of neS is plotted in Figure 9 and demonstrates the paths for the conducting strips and/or wires. The contours of y will be demarcation curves of varying areas of capacitance.
One mechanism to increase capacitance is to make an RF shield a contoured surface rather than just a cylindrical sheet, such that at regions of higher capacitance, the shield is closer to the coil than at regions where the capacitance is lower. Another mechanism is to add dielectric material of high permittivity at regions of higher capacitance compared to other regions, varying according to Y.
Figure 10 shows a contour plot of y for the same current density as shown in Figure 8. The contours in this case are not current paths but are markers to indicate the regions of increasing charge density and hence capacitance. The region of highest capacitive effect is in the center of the contours. Note that J,,rl-free is actually perpendicular to these contour lines.
As further examples of the above procedures, Figures 11 through 15 show the results of determining inductive and capacitive elements for complex current density functions obtained using the second method aspect of the invention.
The design objective was to produce an RF coil 20cm in diameter, 25cm in length with a DSV diameter of 10cm, offset by 2.5cm from the Z 0 plane. (The Z 0 plane passes half-way along the length of the coil). The coil was designed at the 20 frequency of 190MHz such that it could be tested in an available MRI machine.
oo: .The second method of the invention was used with the constraint that J_ 0 at z +L This constraint is a factor on how the resulting coil is excited because no current distribution flows out from the edge of the cylinder. After obtaining coefficients of the current density, the function S was calculated. The contours of the stream 25 function noS for half the coil are shown in Figure 11. These are the preliminary patterns for conductor positions. The resulting current density was tested using a :commercial method of moments package the FEKO program distributed by EM Software Systems, South Africa). The coil current density is approximated by Hertzian dipoles while the shield current density is ignored. Instead, a metallic shield o• 0 30 approximated by triangles as per the method of moments is positioned where the shield should be. The normalized magnetic field of the coil in the transverse plane is shown in Figure 12. The field varies within 10% over a distance of 13cm in the x-direction (upper curve) and 12cm in the y-direction (lower curve). The field variation along the z-axis is shown in figure 13 which shows a 10% variation over a distance of 11 cm shifted along the z-axis by 2.5cm. These simulated results generated by Hertzian dipoles approximating coil currents and a metal cylinder for the shield agree with the original field specification and target volume specification in the inverse program.
The coil patterns were then converted to conductor patterns. For this case, 2.2mm diameter wires were used for the conducting paths. Lumped elements were added to the model and adjusted such that the current distribution at the resonant frequency of 190MHz approximated the original calculated current distribution. The resulting coil is shown in Figure 15 with the circuit diagram shown in Figure 14. The values for the capacitors are: CI=3.9pF, C2=4.7pF, C3=3.3pF, C4=2.7pF, C5=2.2pF, C6 1OpF and C7 and C8 are variable capacitors.
When the coil was constructed, it realized an unloaded Q of 139, measured from the 3dB down power points either side of the resonant frequency In all of the examples presented throughout this specification, a homogenous target field has been specified, but non-homogeneous target fields may be required in some circumstances and these are readily handled by the methodology disclosed herein.
The invention has been described with particular reference to asymmetric radio 20 frequency coils and a method for designing such coils, in which the section of interest (the DSV) can be placed at an arbitrary location within the coil. It will also be ••understood by those skilled in the art that various changes in form and detail may be made without departing from the scope of this invention.
Throughout this specification including the claims, where the context permits, the terms "comprises" and "comprising" are to be interpreted in an exclusive sense, as i:i: including the stated integer(s), without necessarily excluding others.
o• .l a S..i

Claims (13)

1. Apparatus for use in a magnetic resonance imaging system for transmitting a radio frequency field, receiving a magnetic resonance signal, or transmitting a radio frequency field and receiving a magnetic resonance signal, said apparatus and said magnetic resonance imaging system having a common longitudinal axis, said apparatus comprising: a support member which defines a bore having first and second open ends which are spaced from one another along the longitudinal axis by a distance L; and a plurality of inductive elements and a plurality of capacitive elements associated with the support member; wherein if used as a transmitter, the apparatus has the following characteristics: the apparatus produces a radio frequency field which has a radial magnetic component which has a spatially-varying peak magnitude whose average value is Ar-avg; (ii) the apparatus has a homogenous volume within the bore over which the spatially-varying peak magnitude of the radio frequency field has a maximum deviation from Ar-avg which is less than or equal to (iii) the homogeneous volume defines a midpoint M which is on the longitudinal axis, is closer to the first end than to the second end, and is spaced from the first end by a distance D such that the ratio D/L is less than or equal to 0.4; and So (iv) at least one of said inductive elements comprises a discrete conductor S. which follows a sinuous path such that during use of the apparatus current flows through a first part of the conductor in a first direction that has both longitudinal and azimuthal components and through a second part of the conductor in a second direction that has both longitudinal and azimuthal components, said first and second directions being different.
2. The apparatus of Claim 1 wherein said first and second open ends have transverse cross-sectional areas A, and A 2 respectively, which satisfy the relationship: JAI-A 2 1/AI 0.1.
3. The apparatus of Claim 1 or 2 wherein AI is at least 2 x 103 cm 2
4. The apparatus of any preceding Claim wherein the bore is cylindrical. 36 The apparatus of any preceding Claim, wherein the homogenous volume is at least 30 x 10 3 cm 3
6. The apparatus of Claim 5 wherein L is at least 1.0 meter.
7. The apparatus of any one of Claims 1 to 4 wherein the spatially-varying peak magnitude has a maximum deviation from Ar,,,avg which is less than or equal to
8. The apparatus of Claim 7 wherein the homogenous volume is at least x 10 3 cm 3
9. The apparatus of Claim 8 wherein L is at least 1.0 meter. The apparatus of any preceding Claim wherein the D/L ratio is less than or equal to 0.25.
11. The apparatus of any preceding Claim wherein the radio frequency field has a frequency of at least 20 megahertz.
12. The apparatus of any preceding Claim wherein the radio frequency field has a frequency of at least 80 megahertz.
13. The apparatus of any preceding Claim wherein each of the plurality of inductive elements comprises a discrete conductor which follows a sinuous path such that during use of the apparatus current flows through a first part of the conductor in a first direction that has both longitudinal and azimuthal components and through a second part of the conductor in a second direction that has both longitudinal and azimuthal components, said first and second directions being different.
14. Apparatus for use in a magnetic resonance system for transmitting a radio frequency field substantially as hereinbefore described with reference to the *accompanying drawings. Apparatus for use in a magnetic resonance system for receiving a C:.•magnetic resonance signal substantially as hereinbefore described with reference to the 0• °accompanying drawings.
16. Apparatus for use in a magnetic resonance system for transmitting a ooo. radio frequency field and receiving a magnetic resonance signal substantially as hereinbefore described with reference to the accompanying drawings.
AU69966/01A 2000-09-11 2001-09-10 Asymmetric radio frequency coils for magnetic resonance Expired AU783991B2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
AU69966/01A AU783991B2 (en) 2000-09-11 2001-09-10 Asymmetric radio frequency coils for magnetic resonance
AU2006201449A AU2006201449B2 (en) 2000-09-11 2006-04-06 Method of Designing Asymmetric Radio Frequency Coils

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
AUPR0059 2000-09-11
AUPR0059A AUPR005900A0 (en) 2000-09-11 2000-09-11 Asymmetric radio frequency coils for magnetic resonance
AU69966/01A AU783991B2 (en) 2000-09-11 2001-09-10 Asymmetric radio frequency coils for magnetic resonance

Related Child Applications (1)

Application Number Title Priority Date Filing Date
AU2006201449A Division AU2006201449B2 (en) 2000-09-11 2006-04-06 Method of Designing Asymmetric Radio Frequency Coils

Publications (2)

Publication Number Publication Date
AU6996601A AU6996601A (en) 2002-09-12
AU783991B2 true AU783991B2 (en) 2006-01-12

Family

ID=25636115

Family Applications (2)

Application Number Title Priority Date Filing Date
AU69966/01A Expired AU783991B2 (en) 2000-09-11 2001-09-10 Asymmetric radio frequency coils for magnetic resonance
AU2006201449A Ceased AU2006201449B2 (en) 2000-09-11 2006-04-06 Method of Designing Asymmetric Radio Frequency Coils

Family Applications After (1)

Application Number Title Priority Date Filing Date
AU2006201449A Ceased AU2006201449B2 (en) 2000-09-11 2006-04-06 Method of Designing Asymmetric Radio Frequency Coils

Country Status (1)

Country Link
AU (2) AU783991B2 (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0546622A1 (en) * 1991-12-11 1993-06-16 Koninklijke Philips Electronics N.V. Magnetic resonance apparatus comprising a bird-cage RF coil
EP0571211A1 (en) * 1992-05-22 1993-11-24 General Electric Company Asymmetric radio frequency coil for magnetic resonance imaging

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2295020B (en) * 1994-11-03 1999-05-19 Elscint Ltd Modular whole - body gradient coil
US5777474A (en) * 1996-11-08 1998-07-07 Advanced Imaging Research, Inc. Radio-frequency coil and method for resonance imaging/analysis
WO2000006245A2 (en) * 1998-07-31 2000-02-10 Case Western Reserve University Method for correction of mri motion artifacts and main field fluctuation
AUPQ198899A0 (en) * 1999-08-03 1999-08-26 University Of Queensland, The A method of magnet design and magnet configuration

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0546622A1 (en) * 1991-12-11 1993-06-16 Koninklijke Philips Electronics N.V. Magnetic resonance apparatus comprising a bird-cage RF coil
EP0571211A1 (en) * 1992-05-22 1993-11-24 General Electric Company Asymmetric radio frequency coil for magnetic resonance imaging

Also Published As

Publication number Publication date
AU6996601A (en) 2002-09-12
AU2006201449B2 (en) 2008-09-18
AU2006201449A1 (en) 2006-05-04

Similar Documents

Publication Publication Date Title
Turner Gradient coil design: a review of methods
Giovannetti et al. A fast and accurate simulator for the design of birdcage coils in MRI
Jin et al. On the SAR and field inhomogeneity of birdcage coils loaded with the human head
Bogdanov et al. Coupled microstrip line transverse electromagnetic resonator model for high‐field magnetic resonance imaging
US20080054904A1 (en) Radiofrequency magnetic field resonator and a method of designing the same
EP0803737A2 (en) Radio frequency coils
US7348780B2 (en) Asymmetric radio frequency coils for magnetic resonance
JPH05269099A (en) Transverse gradient coil for head imaging
Leifer Theory of the quadrature elliptic birdcage coil
Lawrence et al. A time-harmonic inverse methodology for the design of RF coils in MRI
Crozier et al. Optimized small-bore, high-pass resonator designs
Fujita et al. A hybrid inverse approach applied to the design of lumped-element RF coils
JP2002512098A (en) Gradient magnetic field projection
AU783991B2 (en) Asymmetric radio frequency coils for magnetic resonance
Tomasi et al. Fast optimization of a biplanar gradient coil set
US7259562B2 (en) RF coil for a highly uniform B1 amplitude for high field MRI
Ghaly et al. Design, Simulation, Modeling, and Implementation of a Square Helmholtz Coil in Contrast with a Circular Coil for MRI Applications.
Ghaly et al. Design and modeling of a radiofrequency coil derived from a helmholtz structure
US4899108A (en) High frequency coil
GB2407388A (en) Method of designing asymmetric radio frequency coils for magnetic resonance
Zhilichev Analytical solutions of eddy-current problems in a finite length cylinder
While et al. 3-D gradient coil design—Initial theoretical framework
Forbes et al. An analysis and optimization of elliptical RF probes used in magnetic resonance imaging
Zhou et al. Modeling the self-capacitance of an irregular solenoid
Maunder et al. Metamaterial liner for MRI excitation—Part 2: Design and performance at 4.7 T

Legal Events

Date Code Title Description
MK14 Patent ceased section 143(a) (annual fees not paid) or expired