The field of the invention is systems and methods for spatially encoding nuclear magnetic resonance signals, such as may be used in magnetic resonance imaging (“MRI”).
Interventional procedures such as the crossing of chronic total occlusions (“CTOs”) using a wire could benefit a great deal from new imaging methods that enable the visualization of the device and the vessel wall during the procedure. This would be an improvement over the current x-ray methods for which the vessel is not visible due to a lack of blood flow. MRI has been proposed as a solution, but has not yet become widely used. One problem with using MRI for these applications is that there is a lack of interventional devices that enable imaging at the high spatial resolution required to adequately image the vessel wall and lumen.
In conventional MRI systems, magnetic field gradients are established to spatially encode nuclear magnetic resonance signals emanating from the volume-of-interest being imaged. These magnetic field gradients are stationary relative to physiologic movements occurring within the volume-of-interest. Thus, compensation for this motion either must be performed prospectively by limiting the times at which data is acquired to those times when little to no motion is occurring, or must be performed retrospectively by correcting the acquired data for the effects of the motion that occurred during data acquisition. In cardiac MRI, compensation for motion typically requires data acquisition during short quiescent periods where the heart is approximately motionless.
These limitations and problems are removed if the magnetic field gradients used to spatially-encode the nuclear magnetic resonance signals are able to move with the physiologic motion, such that a local frame of reference is created. Thus, there is a need to provide a system and method for spatially-encoding nuclear magnetic resonance signals in a local frame of reference, such as one that is allowed to move separately from an MRI system.
The present invention overcomes the aforementioned drawbacks by providing a device for spatially encoding nuclear magnetic resonance signals that is capable of moving with local motion sources while producing spatial-encoding magnetic fields, thereby resulting in the spatial encoding of these signals in a local frame of reference free of the local motion sources. By way of example, such local motion sources may include physiological motion. The present invention is also useful for other magnetic resonance imaging and spectroscopy applications in which it is advantageous to spatially-encode and acquire data from a local frame of reference that moves along with the device making the measurements. An example of such other applications includes oil well logging applications.
It is an aspect of the invention to provide a device for spatially-encoding nuclear magnetic resonance signals. The device includes a plurality of spatial-encoding elements configured to produce a spatial-encoding magnetic field in the presence of an external magnetic field. The plurality of spatial-encoding elements include at least one paramagnetic spatial-encoding element and at least one diamagnetic spatial-encoding element. The device further includes a support that is coupled to the plurality of spatial-encoding elements. The support is configured to space the plurality of spatial-encoding elements in a fixed arrangement and is also configured to move the spatial-encoding elements relative to each other in order to produce a plurality of different spatial-encoding magnetic fields in the presence of an the external magnetic field.
It is another aspect of the invention that the spatial-encoding elements may be shaped as spheres, spherical frustums, cones, conical frustums, circular discs, and the like.
It is yet another aspect of the invention that the device may include a radio frequency receiver coil coupled to the support and configured to receive nuclear magnetic resonance signals from a volume-of-interest adjacent the support.
It is yet another aspect of the invention that the device may include a rotation tracking system to measure a rotation of the device about its longitudinal axis. Thus, the device may include marks that indicate a rotational orientation of the support with respect to rotation about the longitudinal axis. The device may include an optical device coupled to the support and configured to identify the rotational orientation of the support by analyzing the marks. The device may include an optical fiber that couples the optical device to the support. The device may also include a sheath disposed about the support and configured to hold the optical fiber in a fixed position with respect to the support.
It is yet another aspect of the invention that the device may include a tracking system coupled to the support. Such a tracking system is configured to measure a yaw angle and a pitch angle of the device with respect to the external magnetic field. An example of a tracking system includes a first tracking member having a first magnetic susceptibility and a second tracking member having a second magnetic susceptibility different than the first magnetic susceptibility. The first tracking member is configured to be moved relative to the second tracking member, thereby altering local magnetic fields produced by the tracking system when the tracking system is positioned in the external magnetic field.
It is also an aspect of the invention to provide a method for producing images of a volume-of-interest with a magnetic resonance imaging (“MRI”) system using a device that includes spatial-encoding elements that produce spatial-encoding magnetic fields when exposed to a main magnetic field of the MRI system. The method includes providing the device to the volume-of-interest within the MRI system and operating the device to adjust the device into a configuration that produces a spatial-encoding magnetic field in response to a main magnetic field of the MRI system. The MRI system is used to acquire signal data from spins adjacent the device, in which the signal data is spatially encoded by the spatial-encoding magnetic field produced by the device. The foregoing steps are repeated while adjusting the configuration of the device to produce different spatial-encoding magnetic fields. An image of the volume-of-interest is then reconstructed from the acquired signal data.
The foregoing and other aspects and advantages of the invention will appear from the following description. In the description, reference is made to the accompanying drawings which form a part hereof, and in which there is shown by way of illustration a preferred embodiment of the invention. Such embodiment does not necessarily represent the full scope of the invention, however, and reference is made therefore to the claims and herein for interpreting the scope of the invention.
When exposed to an external magnetic field, such as the main magnetic field, B0, of an MRI system, paramagnetic materials become magnetized such that their magnetization is oriented in the same direction as the external magnetic field. Conversely, diamagnetic materials become magnetized such that their magnetization is oriented opposite the direction of the external magnetic field. This magnetization of paramagnetic and diamagnetic materials results in magnetic field perturbations. As an example, for a spherical object composed of a material with magnetic susceptibility, x, placed in the main magnetic field, B0, of an MRI system will produce a magnetic field perturbation, ΔBz, as follows:
where a is the radius of the spherical object and χ is the volume susceptibility of the spherical object.
It is an aspect of the present invention that objects with different magnetic susceptibilities can be arranged in certain configurations such that specific magnetic field perturbations will be generated when the objects are positioned within an external magnetic field. The inhomogeneous magnetic fields resulting from these magnetic field perturbations may then be used to spatially encode nuclear magnetic resonance signals; thus, when objects of different magnetic susceptibilities are so arranged, they may be referred to as spatial-encoding elements. By way of example, these spatial-encoding elements are incorporated into a device to allow the spatial localization of nuclear magnetic resonance signals generated in a volume-of-interest adjacent the device. The generated field perturbations will depend on the angle of such a device with respect to the external magnetic field in which it is positioned. An example of the magnetic fields produced by such spatial-encoding elements is illustrated in
The aforementioned magnetic field perturbations will change as the device is rotated about a roll axis by a roll angle, θ. Thus, by way of example, the device may be rotated through 360 degrees in a series of angular steps. At each angular step, spins in the volume-of-interest adjacent the device may be excited and a resulting nuclear magnetic resonance signal acquired. The device is rotated again by the angular step and more signals are acquired. By rotating the device through 360 degrees, sufficient information can be acquired to reconstruct an image of a volume-of-interest adjacent the device.
The signal value at each time point is the summation of all small signals from a continuous medium:
s(θk,tl)=∫x∫y∫zC(θk,x,y,z)ρ(x,y,z)e−iω(θ
where ρ(x,y,z) is the spin density of the medium; C(θk,x,y,z) is a time-independent scaling coefficient that incorporates flip-angle variations that may result from the excitation profile of the RF excitation pulse, the coil sensitivity profile, and so on; and ω(θk,x,y,z) is the three-dimensional angular frequency distribution due to the magnetic field perturbations generated by the spatial-encoding elements.
The linear relationship between signal and magnetization provides the a reconstruction process that is equivalent to finding a decoding matrix, F, that produces a magnetization vector with elements {tilde over (ρ)}n, where n=1, 2, . . . N for N voxels in the three-dimensional reconstruction grid. The reconstruction process using the decoding matrix, F, for voxel {tilde over (ρ)}n is given by the following equation:
Because the spatial-encoding fields generated by the spatial-encoding elements are nonlinear and non-bijective, the for a reconstructed voxel may be contaminated by signals from neighboring voxels. The spatial distribution of the signal that contributes to a particular reconstructed voxel, {tilde over (ρ)}n, can be characterized by a spatial response function (“SRF”) of the voxel:
where exp(−iω(θk,x,y,z)tl) is the phase of the magnetization for roll angle θk and time point tl. Practically, the SRF may be viewed as a linear combination of these phase maps.
As an example, the reconstruction technique may be viewed as computing the rows of the decoding matrix, Fnkl, for n=1, . . . , N, such that the SRF for each voxel is well localized to that voxel. Each row of the decoding matrix may then be used to compute the intensity of one voxel of the final reconstructed image using Eqn. (3). The rows of the decoding matrix may be calculated in one of at least two ways.
One approach for calculating the coefficients of the decoding matrix involves generating an approximation of an encoding matrix, E, from discretized phase evolution maps. An example of such an encoding matrix, E, is as follows:
where rn is the locations of a group of N spatial points that are to be reconstructed, which can be expressed as the following:
The encoding matrix, E, provided in Eqn. (5) assumes that there are N discrete voxels in the reconstructed image, that signals were acquired for K different roll angles, and that signals were recorded at L discrete sample times ending at tL. Elements of the encoding matrix may be calculated using a discrete representation of magnetic field perturbations surrounding the spatial-encoding elements. By way of example, the elements of the encoding matrix may be calculated for a two-dimensional plane of voxels located two millimeters in front of the spatial-encoding elements. It will be appreciated by those skilled in the art that this approach can be extended to three spatial dimensions, and that it is important to have elements in the encoding matrix that correspond to all of the spatial locations that contribute to the signal.
If the N points, rN, adequately cover the volume that contributes signal to s(θk,tl), then a good numerical approximation, {tilde over (s)}(θk,tl), of this signal can be given by the following equation:
{tilde over (s)}kl=Elkn{tilde over (ρ)}n (7).
The signal values, {tilde over (s)}kl, at each time point, tl, for each roll angle, θk, are calculated by multiplying one row of the encoding matrix
E
kl
=[C(θk,r1)e−iω
by the magnetization density vectors, {tilde over (ρ)}. In reality, the signal values at each time point will be a summation of all small signals from a continuous medium, as shown in Eqn. (2); thus, Eqn. (7) is a discretization of the continuous signal equation.
Because spatial encoding is a linear process, the reconstructed image voxel values, arranged as a vector in {tilde over (ρ)}, can be computed by solving a system of over-determined linear equations:
{tilde over (s)}=E{tilde over (ρ)} (9);
or equivalently by computing the decoding matrix:
{tilde over (ρ)}=F{tilde over (s)} (10).
If the spatial-encoding magnetic fields are linear, Eqns. (9) and (10) could be solved by performing a Fourier transform on the acquired signals. However, because of the nonlinear nature of the spatial-encoding magnetic fields, there is no one-to-one relationship between the acquired signals and spatial frequencies. Therefore, a k-space representation of the acquired signals will not be applicable in this technique. The encoding matrix, E, of Eqn. (5) is usually ill-conditioned, so taking the pseudo-inverse of the encoding matrix to compute the decoding matrix, F, will typically result in a poor image reconstruction. To improve upon this result, a cascade regularization technique may be used to reduce the ill-posed nature of the encoding matrix. In addition, to improve the accuracy of the image reconstruction process, truncated singular value decomposition (“TSVD”) and Tikhonov regularization may be used.
The reconstruction technique described above is considered a forward-problem solution for the spatial-encoding technique provided by the present invention. As described earlier, there is signal contamination in a reconstructed voxel that originates from other spatial locations. In the forward-problem solution described above, each element of the encoding matrix is an average of the evolved phase over the corresponding voxel, and this assumption will fail when a sharp gradient in the phase exists within a voxel. This inhomogeneous voxel then becomes a source of signal contamination for the other voxels and is itself reconstructed with attenuated signal. Therefore, it is advantageous to develop other reconstruction methods that can yield well-localized SRFs for all voxels.
In another approach for calculating the decoding matrix coefficients, an optimization technique can be used to solve for the decoding matrix coefficients in one row of the decoding matrix. As an example, this optimization may be done to minimize the least-square-error between an ideal SRF and a calculated SRF for a voxel, {tilde over (ρ)}n. It will be appreciated by those skilled in the art that other optimization functions may be used. For example, the l1-norm could be minimized in a sparse domain. A least-squares minimization will have the following form:
If a good global minimum is found, the SRF from the resulting optimized reconstruction coefficients will be fairly localized and signal contaminations from other spatial locations will be minimized. Global minima are only found when the target (ideal) SRF is realizable, given the particular set of phase maps. Thus, it is important to choose the appropriate size, shape, and position of each target SRF, corresponding to each desired voxel.
Once all of the decoding matrix coefficients for all voxels in the reconstruction grid have been computed, Eqn. (3) can be used to reconstruct a magnetization density matrix of size N. An approach such as the one just described can be considered as an inverse-problem solution for the spatial-encoding technique provided by the present invention.
Having described the spatial-encoding and image reconstruction processes, a discussion of various configurations of the spatial-encoding elements and their incorporation into different devices is now provided.
With initial reference to
The spatial encoding elements 102 may be arranged such that they are all collinear, as illustrated in
Intravoxel phase dispersions may be caused by inhomogeneous magnetic fields generated by the spatial-encoding elements 102 in the out-of-plane direction. The magnetic field perturbations generated by the spatial-encoding elements 102 vary significantly in the out-of-plane direction, which can result in signal loss in voxels. Thus, it may be advantageous to define a sensitive region of the device 100. One example of the sensitive region of the device 100 is illustrated in
As a practical matter, when the spatial-encoding elements 102 are arranged such that there is an inherent symmetry of the magnetic field perturbations produced by the spatial-encoding elements with respect to the center of the arrangement, images reconstructed using such a device may include image artifacts that manifest as replication of magnetization densities. To mitigate these artifacts, the spatial-encoding elements may be arranged such that there is a disconnect in the inherent symmetry of the magnetic field perturbations generated by the spatial-encoding elements. One example of such an arrangement is the purely collinear arrangement illustrated in
Optionally, the device 100 may be surrounded by a bio-compatible layer 160 so that the magnetic materials, such as the spatial-encoding elements 102, do not come into direct contact with blood and tissue. In addition, the bio-compatible layer 160 may be configured to provide a gap between the spatial-encoding elements 102 and the exterior of the device 100 such that the local magnetic fields in the volume-of-interest are less steep than they would be without the gap.
By way of example, when the spatial-encoding elements 102 are arranged as illustrated in
With the spatial-encoding elements 102 arranged as illustrated in
Referring now to
To measure the roll angle of the device 100, the device 100 may be configured to incorporate a rotation tracking system, as illustrated in
The yaw and pitch of the device can be measured in a number of different ways. By way of example, and referring to
In other configurations of the device, such as the one illustrated in
As described above, the disclosed spatial-encoding elements 102 produce magnetic field perturbations when placed in an external magnetic field. These magnetic field perturbations alter the resonance frequency of spins that exists at spatial locations affected by the magnetic field perturbations. By way of example, the spatial-encoding elements 102 produce magnetic field perturbations that result in resonance frequency offsets between −2000 Hz and +2000 Hz for an external magnetic field with a strength of 1.5 T. Therefore, a free-induction decay (“FID”) occurring in a region affected by the magnetic field perturbations will contains resonance frequencies in the range defined by these resonance frequency offsets. By rotating the device 100 about its roll axis 110, the frequency content of those FIDs that occur in a volume-of-interest adjacent the device 100 will remain the same. However, spins that do not experience the magnetic field perturbations generated by the spatial-encoding elements 102 will precess at their Larmor frequency, ω0=γB0, and their phases will contain no spatial information. Thus, exciting these spins will degrade the accuracy of the reconstructed image. Off-resonance excitation can be used to excite only those spins that are experiencing magnetic field perturbations due to the device 100 in order to restrict the region from which signals are acquired. As mentioned earlier, another approach to restricting the signal region is to incorporate RF receiver coils for signal acquisitions into the device 100.
To perform off-resonance excitation, a radio frequency (“RF”) pulse having frequency content that matches the frequency offsets noted above is used to excite the desired spins. This RF excitation pulse may be, for example, a composite RF pulse with a binomial distribution designed using the Shinnar-LeRoux algorithm. Spins that are precessing with different Larmor frequencies will experience different flip angles due to the excitation profile of these off-resonance RF excitation pulses, which are frequency dependent. It is contemplated that this excitation profile will lead to encoding matrix rows that are modulated by the flip angle profile of the RF excitation pulse, which will improve the reconstruction process.
In conventional slice-selective RF excitation, the RF excitation profile will be similar to the Fourier transform of the RF pulse in presence of a slice-selective gradient field. With this new technique, the magnetic field perturbations will assume the functions of the slice-selective gradient; however, because of the nonlinear nature of these magnetic field perturbations, the RF excitation profile will not resemble the Fourier transform of the RF excitation pulse. Thus, the off-resonance RF excitation pulse described above will excite only those spins precessing at Larmor frequencies within the bandwidth of the proposed off-resonance RF excitation pulse.
As described above, the magnetic field perturbations generated by the spatial-encoding elements 102 are dependent on the orientation of the device 100 with respect to the externally applied magnetic field. Therefore, it is desirable to provide a method for effectively measuring the orientation of the device 100 with respect to the external magnetic field. The measured yaw, pitch, and rotation angles can be used when calculating the encoding or decoding matrices to correct for any deviation of the device 100 from the external magnetic field. It is noted, however, that image reconstruction will not be significantly affected when the device 100 deviates from its alignment with the external magnetic field by a yaw angle, pitch angle, or both of plus or minus fifteen degrees. However, if the device 100 undergoes significant motion (e.g., inside coronary arteries), it may be necessary to measure the angle between the device 100 and the external magnetic field.
As the device 100 rotates around its roll axis 110, a positive image contrast in the plane perpendicular to the roll axis 110 will be changed as a result of the rotation of the magnetic field perturbations with the device 100. In addition to using a rotation tracking system incorporated into the device 100 as discussed above, the roll angle of the device 100 can be determined from two projection images of the device 100 that are acquired with the device rotated to two different roll angles, θ1 and θ2. It is contemplated that there is a bijective relationship between the roll angle, θ, of the device 100 and the length of the magnetic field perturbations at this angle. The length of the magnetic field perturbations can be calculated from the size of a hyperintense region present in a projection image of the device 100. The roll angle, θ, of the device can be determined from these two projection images as follows:
where a is the length of the magnetic field perturbations at the first roll angle and b is the length of the magnetic field perturbations at the second roll angle. It is possible to generalize this concept to measure both the yaw and pitch angles of the device 100 as well.
Without loss of generality and by way of example, assume that the device 100 is initially oriented along the z-axis and is parallel to the external magnetic field, B, and that the device is rotated to a roll angle of θ1. If the device 100 is moved, the following equation can be used to establish a relationship between the initial position of the device 100 and the new position of the device 100 using rotational angles:
However, the initial roll angle, θ1, is unknown and Eqn. (13) cannot be solved for the yaw and the pitch angles, ψ and Ω, respectively. To resolve the roll angle, the 100 device may be rotate by a known angle, such as ninety degrees, about the roll axis and the new position of the 100 device may be computed.
By combining Eqn. (13) and Eqn. (14), the yaw, pitch, and roll angles can be computed. A total of six measurements (x1,y1,z1,z2,y2,z2) are required to establish the following sets of linear equations.
The angle measurement technique may be used to locate the spatial position and orientation of device 100 under MRI guidance; however such a technique may suffer from its long acquisition time. For example, the long acquisition time may be unacceptable in some applications, such as positioning an interventional device in coronary arteries where there is significant motion in a short period of time. To address this issue, one-dimensional positive contrast projections can be obtained along each of three principal axes and used both to estimate the coordinates of the device and to measure the length of the hyperintense region adjacent the device.
One-dimensional projections are commonly used for positioning devices under MRI guidance, and can significantly reduce the acquisition time needed for localization. By way of example,
An example of a pulse sequence that may be used to obtain one-dimensional projections of the device 100 is illustrated in
Using the foregoing pulse sequence, the angle of the device 100 can be measured with respect to each principal axis before each data acquisition step. As described above, knowing the orientation of the device with respect to external magnetic field can be advantageous to improve the accuracy of the image reconstruction techniques described herein. With this proposed technique, the angle of the device 100, and therefore the angle of the spatial-encoding elements 102, can be measured respect to the external magnetic field before each angle-step of the device 100 in order to perform corrections to the decoding matrix, F. It is noted that the roll angle may be more precisely measured using a rotation tracking system, such as the optical system described above. For practical reasons, it is contemplated that using such a system instead of measuring the roll angle with an MRI system will be preferable.
Referring particularly now to
The pulse sequence server 1110 functions in response to instructions downloaded from the workstation 1102 to operate a gradient system 1118 and a radiofrequency (“RF”) system 1120. Gradient waveforms necessary to perform the prescribed scan are produced and applied to the gradient system 1118, which excites gradient coils in an assembly 1122 to produce the magnetic field gradients Gx, Gy, and Gz used for position encoding MR signals. The gradient coil assembly 1122 forms part of a magnet assembly 1124 that includes a polarizing magnet 1126 and a whole-body RF coil 1128.
RF excitation waveforms are applied to the RF coil 1128, or a separate local coil (not shown in
The RF system 1120 also includes one or more RF receiver channels. Each RF receiver channel includes an RF amplifier that amplifies the MR signal received by the coil 1128 to which it is connected, and a detector that detects and digitizes the I and Q quadrature components of the received MR signal. The magnitude of the received MR signal may thus be determined at any sampled point by the square root of the sum of the squares of the I and Q components:
M=√{square root over (I2+Q2)} (18);
and the phase of the received MR signal may also be determined:
The pulse sequence server 1110 also optionally receives patient data from a physiological acquisition controller 1130. The controller 1130 receives signals from a number of different sensors connected to the patient, such as electrocardiograph (“ECG”) signals from electrodes, or respiratory signals from a bellows or other respiratory monitoring device. Such signals are typically used by the pulse sequence server 1110 to synchronize, or “gate,” the performance of the scan with the subject's heart beat or respiration.
The pulse sequence server 1110 also connects to a scan room interface circuit 1132 that receives signals from various sensors associated with the condition of the patient and the magnet system. It is also through the scan room interface circuit 1132 that a patient positioning system 1134 receives commands to move the patient to desired positions during the scan.
The digitized MR signal samples produced by the RF system 1120 are received by the data acquisition server 1112. The data acquisition server 1112 operates in response to instructions downloaded from the workstation 1102 to receive the real-time MR data and provide buffer storage, such that no data is lost by data overrun. In some scans, the data acquisition server 1112 does little more than pass the acquired MR data to the data processor server 1114. However, in scans that require information derived from acquired MR data to control the further performance of the scan, the data acquisition server 1112 is programmed to produce such information and convey it to the pulse sequence server 1110. For example, during prescans, MR data is acquired and used to calibrate the pulse sequence performed by the pulse sequence server 1110. Also, navigator signals may be acquired during a scan and used to adjust the operating parameters of the RF system 1120 or the gradient system 1118, or to control the view order in which k-space is sampled. The data acquisition server 1112 may also be employed to process MR signals used to detect the arrival of contrast agent in a magnetic resonance angiography (“MRA”) scan. In all these examples, the data acquisition server 1112 acquires MR data and processes it in real-time to produce information that is used to control the scan.
The data processing server 1114 receives MR data from the data acquisition server 1112 and processes it in accordance with instructions downloaded from the workstation 1102. Such processing may include, for example: Fourier transformation of raw k-space MR data to produce two or three-dimensional images; the application of filters to a reconstructed image; the performance of a backprojection image reconstruction of acquired MR data; the generation of functional MR images; and the calculation of motion or flow images.
Images reconstructed by the data processing server 1114 are conveyed back to the workstation 1102 where they are stored. Real-time images are stored in a data base memory cache (not shown in
Having described the general structure of the provided device, and various examples of configurations thereof, a description of a general operation of the device is now provided. Referring now to
Advantageously, the device is free to move in conjunction with physiological motion, such as respiration, cardiac rhythm, and pulsatile flow. As a result, the spatial-encoding magnetic fields move along with the device and in conjunction with physiological motions. Because of this result, the device effectively images the VOI in a frame of reference that moves along with present physiological motions. By imaging in this local frame of reference, motion-related errors in the acquired image data are mitigated, and substantially motion-artifact free images can be reconstructed from the image data.
A determination is made at decision block 1208 as to whether the desired amount of image data has been acquired. If not, then the device is positioned in a different position, as indicated at step 1210, and the MRI system operated again to acquire more image data. By way of example, a new roll, pitch, and yaw angle of the device may then be calculated in this new position. In each different position, the device provides a different spatial encoding to the acquired image data. By acquiring a plurality of image data sets, each with a different spatial encoding, one of the image reconstruction techniques presented above can be solved. Thus, after the desired amount of image data has been acquired, images of the VOI are reconstructed, as indicated at step 1212. By way of example, reconstruction coefficients based on the measured roll, pitch, and yaw of the device may be calculated and used in the reconstruction process. It is noted that the image data acquisition and image reconstruction processes can take place in real-time, so as to provide visual feedback to a user, such as a medical practitioner performing an interventional or endoscopic procedure. For real-time display applications, pre-computed reconstruction coefficients can be used and matched to the measured roll, pitch, and yaw of the device using a look-up table.
The present invention has been described in terms of one or more preferred embodiments, and it should be appreciated that many equivalents, alternatives, variations, and modifications, aside from those expressly stated, are possible and within the scope of the invention.
This application claims the benefit of U.S. Provisional Patent Application Ser. No. 61/478,122, filed on Apr. 22, 2011, and entitled “SYSTEM AND METHOD FOR MAGNETIC RESONANCE IMAGING.”
Number | Date | Country | |
---|---|---|---|
61478122 | Apr 2011 | US |