System and method for providing real-time motion correction by utilizing navigators

[0002] The present invention relates to a system and method for tracking and correcting, in real-time, a motion of an anatomical object. In particular, the system and method uses octant navigators to track and correct certain features of the anatomical object (e.g., a brain) during a three dimensional scan thereof.

[0003] Several approaches have been utilized to correct for a movement of a subject during scans. Such conventional methods require an off-line, post-processing of the images. Indeed, an image correction is not theoretically possible in some cases because certain areas of k-space may not be sampled, while others are oversampled. The earliest navigator-based methods for the motion correction utilized straight-line navigators to detect a linear motion. Such technique may be useful in chest examinations where the diaphragm and associated organs translate along a particular axis. However, these conventional methods do not quantify and determine the magnitudes or degrees of rotations of the objects being examined, or portions thereof.

[0004] Other conventional systems and methods use correlated volumes to calculate rotations and translations. While this technique allows for a correction of certain types of motions, it cannot be applied to three dimensional single image structural scans. Indeed, the techniques employed by these systems and methods are computationally intensive. Similar approaches have been used on-line and off-line for two dimensional sequences. In particular, slices can be correlated to calculate the transformation from one slice to the next, and correct for such changes. However, while this conventional method has certain benefits, the correction of the motion provided therein only applies to the object motion for an in-plane movement.

[0005] Another conventional approach for the motion correction utilizes a set of navigator pulses. For example, circular navigators have been proposed to gauge a rotation within sequences. The circular navigators are 30 ms in duration, and three of these circular navigators are required to characterize the motion of the object about three cardinal axes. An iterative approach is then taken to correct the motion of the object since the navigators are not usable for the out of plane rotations. Therefore, this conventional method provides an approximate correction. Another set of three circular navigators is then acquired for a final (fine) adjustment. While this approach fulfills the theoretical need to compensate for all three axes of motions, there are certain impracticalities associated with such conventional approach. Most importantly, the entire procedure is relatively time-consuming, i.e., adding this procedure to, e.g., a standard functional magnetic resonance imaging (“fMRI”) acquisition sequence is not practical.

[0006] It is therefore the object of the present invention to use navigator echoes during each radio frequency acquisition to track the motion of an object (e.g., the head of a patient) and to compensate for such motion by correcting the gradients during the scan thereof. Another object of the present invention is to utilize the navigators to track and correct for a motion of the object in real-time. Yet another object of the present invention is to provide an improvement in the basic signals to measure the brain physiologic changes that are being investigated. Another object of the present invention is to improve the quality of anatomical images obtained. Still another object of the present invention is to assess (and possibly measure and record) the rigid-body motion (e.g., rotations and translations in 3 dimensions) without necessarily correcting for it.

[0007] These and other objects of the present invention are achieved with a method and system for correcting for a motion of an object. With these system and method, the navigator data and map data are obtained for the object. Then, the navigator data is compared to the map data to generate comparison data. A translation and/or a rotation of the object is corrected, in real time, as a function of the comparison data. The navigator can be preferably an octant navigator.

[0008] In one exemplary embodiment of the present invention, a scanning sequence can be used to determine a position of the object. This scanning sequence may include a signal portion which includes at least one radio frequency signal, an octant navigator portion which includes at least one octant navigator, and a spoiler portion provided for reducing a signal magnitude of the scanning sequence. The octant navigator is provided for allowing a measurement of at least one of a rotation and a translation of the object. The navigator portion is advantageously provided between the signal portion and the spoiler portion.

[0009] In another advantageous embodiment of the present invention, the octant navigator is a navigator that traces the outline of an octant on the surface of a sphere in k-space. This enables a rotation about the three cardinal axes, and a translation in all three directions to be achieved in a single read after a single radio frequency pulse.

[0010] It is also advantageous, according to the present invention, to first obtain a pre-mapping of the region of k-space a small number of degrees in each direction from the initial octant navigator. This pre-mapping technique preferably eliminates the need for an iterative, approximate solution. By comparing the actual navigator with a local pre-mapped k-space map, it is possible to determine the true rotations and translations using a single subsequent octant navigator. It is preferable that the rotations are not so large as to shift the navigator outside the boundaries of the map. Since the gradients can be preferably corrected after each navigator is acquired, the accumulated rotation may exceed the limits of the map, provided that the inter-navigator rotation is constrained by the limits of the map.

[0011] In still another embodiment of the present invention, the pre-mapping technique may be combined with an optimal radius selection method using concentric octant navigators. This procedure selects the optimal radius for an individual object to be used for the octant navigators in the map, and along with the subsequent navigators. Also, the optimal offset angles relative to each of the three axes may be selected. The octant navigators can also be inserted between the phase steps of a three-dimensional sequence to obtain a single motion corrected structural scan. They can also be inserted within and/or between the volumes of a three dimensional echo planar image used for functional imaging or diffusion imaging. Further, a slab selection may be used to eliminate parts that do not move as a rigid body with the part of interest.

[0012] According to an exemplary embodiment of the present invention, the path of the octant navigant does not have to be the edge of an octant. In fact any path on a sphere in k-space may be traced. Indeed, any path that extends in all three dimensions (i.e., preferably not a path that lies completely within a single plane) can be utilized. It is within the scope of the present invention to extract the navigator information. While all paths do not have to be equal or yield the same SNR in the rotation/translation estimates, the use of the octant is preferable. According to the present invention, a preliminary exploration and analysis of the k-space representation of a volume can yield an optimal path. It is possible, with an appropriate spoiling, to implement this exemplary scheme in two dimensional sequences. In this case, more than one map may be collected to gauge the effect on the navigators of the excitation of the different slices due to the imaging part of the sequence.

[0013] One of the exemplary advantages of the present invention is that it utilizes a navigator that is capable of reading in a single path (which may use only one radio frequency (“RF”) pulse that may be shared with the imaging part of the sequence. It is also within the scope of the present invention not to require an additional RF pulse for the navigant which can fully qualify a rigid body motion (rotations/translations) for which the navigant should occupy 3 dimensions, i.e. preferably not contained within a plane. The navigator used by the present invention can be applicable in 2 dimensional and 3 dimensional sequences for motion correction.

[0014] For a more complete understanding of the present invention and its advantages, reference is now made to the following description, taken in conjunction with the accompanying drawings, in which:

[0015] FIG. 1 shows an exemplary embodiment of a system according to the present invention which tracks and corrects a motion of an anatomical object in real time;

[0016] FIG. 2 shows an exemplary embodiment of a method according to the present invention for tracking and correcting a motion of an anatomical object in real time;

[0017] FIG. 3 shows an illustration of an exemplary navigator path in k-space that circumscribes a particular section (e.g., an octant) on a spherical surface of k-space;

[0018] FIG. 4 shows an exemplary sequence structure provided to determine an optimum radius and offset angles for a navigator;

[0019] FIG. 5 shows an exemplary illustration of possible combinations of perturbations that are tested by the system and method according to the present invention;

[0020] FIG. 6 shows an exemplary octant navigator pulse sequence kernel according to the present invention;

[0021] FIG. 7 shows an enlarged illustration of exemplary gradient control signals as provided in the navigator pulse sequence kernel of FIG. 6;

[0022] FIG. 8 shows a detailed illustration of a portion of the method according to the present invention in which optimal radius and angles are determined;

[0023] FIG. 9 shows an exemplary map of a region of k-space surrounding edges of the octant navigator;

[0024] FIG. 10 shows an exemplary structure for the map illustrated in FIG. 9;

[0025] FIG. 11 shows an illustration of the navigators arranged in a predetermined number of partitions, each having a particular number of lines;

[0026] FIG. 12 shows an exemplary portion of the method according to the present invention in which the image of the object is transformed and corrected (in real-time) using the navigators; and

[0027] FIG. 13 shows another exemplary embodiment of the method according to the present invention for tracking and correcting a motion of the anatomical object in real time.

[0028] An exemplary embodiment of a system according to the present invention which tracks and corrects a motion of an anatomical object in real time is illustrated in FIG. 1. As shown, this system includes a processing arrangement 40 which is connected (either wirelessly or via a wired connection) to a sensor arrangement 50. The sensor arrangement 50 can be a conventional magnetic resonance (“MR”) sensor arrangement with which a subject 60 can be scanned. The exemplary conventional sensor arrangement is described in U.S. Pat. No. 6,118,273, the entire disclosure of which is incorporated herein by reference.

[0029] For example, this conventional sensor arrangement generally includes a gradient magnetic field coil which is designed to generate a gradient magnetic field, a radio-frequency (“RF”) coil which is designed to generate a high frequency magnetic field in the magnetic field region, and an RF probe which is designed to detect the MR signals generated from the subject 60. It is known to use gradient field coil units in the gradient magnetic field coil which are provided in three axial directions, crossing perpendicular to one another.

[0030] The sensor arrangement 50 may also include a magnetic field source which can be actuated so that the coil units generate the gradient magnetic fields in response to the signals of the coil. With the gradient field coil units, the magnetic fields in three axial directions (e.g., slicing, phase encoding and readout directions) can be applied to the space where the subject 60 is positioned. Using the sensor arrangement 50, the RF coil generates a high frequency magnetic field in the form of a pulse to respond to the signals generated by an RF transmitter which is also an element of the sensor arrangement 50. An RF probe, which is provided substantially near the subject 60, generates the RF signals which are detected by signal detector of the sensor arrangement 50. Further details of the conventional sensor arrangement 50 are provided in U.S. Pat. No. 6,118,273.

[0031] These RF signals (or data equivalent therefore) are received by the processing arrangement 40 which can be a general purpose computer (e.g., a Pentium® based computer), a specific purpose computer, or more than one computer or combination of computer(s) and special purpose hardware, especially designed to receive and process the data corresponding to the RF signals. After the processing arrangement 40 receives the RF signals, it determines if and to what extent at least a portion of the subject 60 that is being scanned has moved or rotated by tracking the subject 60. If the processing arrangement 40 determines that at least a portion of the subject 60 has moved (e.g., the patient's head), the processing arrangement 40 corrects this motion in real time using the data it obtained regarding such movement. This data can include a series of navigators, and more particularly, octant navigators.

[0032] The navigator is generally known as a scanning sequence which can be used for determining and acquiring positional data. The octant navigator is a navigator which measures a rotation of an object about the three cardinal axes, and the object's translation along each axis. In particular, the octant navigator is a navigator which, using a path through k-space (or a phase space), circumscribes an octant on a spherical surface which can be related to the object. The path can be understood as approximately similar to a trace of a quarter of a circle in each of the x-y, y-z and z-x planes.

[0033] FIG. 2 shows an exemplary embodiment of a top-level diagram of the method according to the present invention for tracking and correcting a motion of an anatomical object, in real time, using the octant navigators. It should be noted that this exemplary method can be implemented by the processing arrangement 40 of the present invention illustrated in FIG. 1. As shown in FIG. 2, the optimal radius and optimal angles for the octant navigators are determined in step 100. Then, in step 110, the map of the octant navigators is obtained, i.e., prior to the scanning procedure. Finally, the images with the interleaved octant navigators are established and collected (step 120). In addition, the image is preferably corrected using the map obtained in step 110, and the motion data corresponding to the motion and rotation of the object, which was determined using the octant navigators.

[0034] FIG. 3 illustrates an exemplary k-space path or trajectory 140 of the octant navigator shown on a surface of a sphere 130 which preferably relates to the object to be tracked. As shown in this drawing, the rotation about all three axes is sampled in a single exemplary sweep. For example, the axis R in the drawing corresponds to right-left axis, A to anterior-posterior axis, and S to superior-inferior axis (which anatomy correspond to the physical slice, phase and read-out directions according to the orientation of the subject in the magnetic field). This exemplary path 140 of the octant navigator is preferably used during step 100 of the method shown in FIG. 2 (i.e., the determination of the optimal radius ρ and angles Θ_Afor the octant navigators). A varying gray-scale bar 145 is provided near the sphere 130 and the path 140 to show the time (e.g., in tenths of milliseconds) that it takes certain portions of the path to be traced and/or sampled. The sampling time can be adjusted according to the desired signal to noise ratio (“SNR”), e.g., a time of less than 1.5 ms for the octant navigator portion is also possible. As shall be described in further detail below, the octant navigators are preferably formed using the gradients of control signals generated by the sensor arrangement 50 on the three axes. These gradients are used in a k-space readout for tracing the edge of the octant of the sphere 130, and to form the path 140 for the octant navigators. The quadrant of the circles (i.e., of the sphere 130) can be achieved by using, e.g., sinusoidal gradients.

[0035] An exemplary sequence pre-scan structure is shown in FIG. 4. This structure provides individual octant navigator pulse sequence kernels 150 (as shall be described in more detail below) for each perturbation (i.e., small angle or movement), for each offset or base angle and for each radius of the path 140. As shown in this drawing, the exemplary structure may include K-number of perturbations for each of L-number of offset/base angles (e.g., 6 base angles) and M-number of radii (e.g., 16 sampled radii). Thus, the sequence of kernels 150 should be tested for K×L×M times. The offset/base angles θ_Acan be a set of large angles, e.g., (π/4,−π/4, 0) (offsets in radians relative to the three axes), and the perturbations PR can be denoted by (Δα, Δβ, Δγ), where Δα, Δβ, and Δγ are small angles such as, e.g., 1°. An exemplary illustration of the perturbations scheme is illustrated in FIG. 5.

[0036] It should be understood that both the accuracy of the information provided by the sequence kernels 150 and the timing to obtain the optimal radius and angles for the octant navigators must be considered. While the accuracy of the information determined from the sequence kernels 150 is enhanced by sampling a large number of perturbations PR, a large number of the base angles θ_Aand a large number of radii ρ, the time to obtain such accurate information can be significant. However, while utilizing a small number of the perturbations PR, base angles θ_Aand radii ρ may reduce the processing time, the accuracy of the results could be reduced to a point of unusability.

[0037] Therefore, the number of the perturbations PR, base angles θ_Aand radii ρ should be set so that an acceptable accuracy is obtained, while the time it takes to determine the optimal radius and angles is also reasonable. In one exemplary embodiment of the present invention, the sequence is tested for 19 perturbations PR for each of the 6 offset/base angles Θ_A. In this exemplary embodiment, the set of 6 base angles θ_AIis provided for each of 16 radii ρ. The 19 perturbations is, e.g., the minimum number of small perturbations required to fit a quadratic function in 3 dimensions to the local error landscape, where the error is between the central and the perturbed navigator signal vectors. Thus, the number of tests performed on the sequence is 19×6×16. It should be understood that other combinations of perturbations, offset/base angle and radii numbers are conceivable, and are within the scope of the present invention.

[0038] FIG. 6 shows a graph of an embodiment of an exemplary octant navigator pulse sequence kernel 150, which includes a first section of gradient control signals 155, a second section of radio frequency (RF) pulses 160 and other (e.g., slab-selective) gradients which is provided immediately before the first portion (i.e., before the section of the gradient control signals), and a third portion for spoiler gradients 170 immediately following the first section 155. Thus, the first section, which includes the gradient control signals is provided between the second and third portions. The exemplary second portion includes three gradients 180, i.e., a phase gradient G_φ, a readout gradient G_roand a slice gradient G_S. An enlarged view of the exemplary enlarged illustration of the gradients which form the octant navigator 200 is shown in FIG. 7. The spoiler gradients 170 are provided immediately following the gradient control signals G_φ, G_ro, G_Sto spoil or reduce the magnitude of the transverse magnetization that might otherwise remain and contribute as noise in the subsequent read of the next image line and navigator. In case of the readout gradient G_ro, the spoiler gradient 170 associated therewith substantially suppresses further echoes in the subsequent blocks of the sequence. FIG. 6 also shows that an analog-to-digital converter (“ADC”) 190 provides a read-out signal, which occurs during a central part of the gradient waveforms. During this time, the octant navigator 200 and the transitions at the navigator's corners are traced out. In an exemplary embodiment of the present invention, the samples acquired during the two corner transitions are discarded, and the remaining three series of samples constitute the octant navigators 200 for the right-left or R axis, anterior-posterior or A axis and superior-inferior or S axis.

[0039] As an example, the flip angle of the RF pulse 160 should preferably match the angle of the sequence in which the octant navigator pulse sequence kernel 150 is embedded, and should match the flip angle of the navigator kernel in the map. If the octant navigator pulse sequence kernel 150 includes a slab-selective pulse, a slab selective gradient pulse can simultaneously be applied on the read out gradient G_ro(e.g., in the superior-inferior direction) to eliminate certain portions of the subject 60 from the navigator due to their inability to rigidly transform when a portion of the subject 60 moves. These unwanted portions may be the neck or lower jaw of a patient. If this is the case, the navigator may preferably be in a block of its own (i.e., with its own RF pulse) so that the image itself will not be also subject to the slab selection and the neck/jaw disappear from the final image volume.

[0040] FIG. 8 shows a flow diagram of an exemplary embodiment of the method according to the present invention to determine the optimal radius and the optimal angles for the octant navigator. In particular, K-number of combinations of perturbations PR are tested (i.e., for each offset/base angle and each radius) in step 300. Then, an error between the octant navigator at a current location and as initially positioned is calculated for each combination in step 310. Upon the completion of step 310, a Hessian matrix is calculated as a mean square difference using the matrix of the calculated errors (in step 320) as follows: $H\left(\Theta ,\rho \right)=\frac{{\partial }^{2}E}{\partial B^2},$

[0041] where

[0042] E is a matrix of errors (E_ijk), E_ijk=N_ijk−N_o,

[0043] N_ijkis the octant navigator at (i.Δα, j.Δβ, k.Δγ),

[0044] N_ois the octant navigator at the initial location, and

[0045] B=(Δα, Δβ, Δγ),

[0046] and i, j, k are all the combinations of i, j, kε{−1,0,1} such that |i|+|j|+|k|≦2, i.e. 19 combinations.

[0047] Then, to determine the optimal angle and radius of the octant navigator (step 330), the eigen values of H(θ,ρ) can be utilized to characterize the error surface for θ, ρ as follows. For example, a set of 3 eigen values can be found for H(θ,ρ). Then, if $D=\prod _{i=1}^3{\lambda }_i$

[0048] (where D is a product of eigen values), the θ, ρ combination for which D is maximal is determined as the optimal selection for the angle and radius of the octant navigator. It is preferable that the method according to the present invention ensures that all eigen values are positive for this D(θ,ρ). Thereafter, the optimal angle θ and the optimal radius ρ are passed to the next step in the process to obtain the navigator map (step 110 shown in FIG. 2).

[0049] Thereafter, the method and system according to the present invention preferably acquires the octant navigator map (step 110 of FIG. 2). The navigator map consists of a series of the navigator kernels 150 rotated at various angles. The angles are adjusted, e.g., in pairs so that lines in k-space parallel to each edge of the octant (on the surface of the sphere 130) is collected. In one exemplary implementation of the present invention, the entire navigator map can be acquired in approximately 3 seconds. An example of the navigator map consisting of 31 lines parallel to each edge of the octant, illustrated on the surface of the sphere 130, is shown in FIG. 9. As shown in this drawing, a map region 350 provides the borders surrounding three edges of the octant navigator. Thus within the map region 350, it is possible to track the motion and/or rotation of at least a portion of the subject to which corresponds the position, angle and direction of the respective octant navigators.

[0050] The navigator map shown in FIG. 9 can be acquired using an exemplary sequence map structure illustrated in FIG. 10 that is illustrated for the entire path 140 of the octant navigator. This exemplary sequence structure illustrates how the region 350 can be mapped by selecting n×n×n (e.g., 5×5×5) combinations of the perturbations PR of the octant navigator 150 (Δα, Δβ, Δγ). In the exemplary 5×5×5 sequence map, every combination of 5 angles close to zero can be acquired therein. Although this procedure measures the data directly, a large amount of the data is most likely not unique. It is preferable to acquire the angles which are provided at a few degrees to either side of each edge of the octant. In this manner the resolution is greater, and the space can be efficiently sampled. The details of this technique shall be described in further detail below.

[0051] For an exemplary n×n×n implementation of the sequence map structure, the structure can be defined as a rotation map (δ) which is the collection of the octant navigators for the following rotations:

[0052] {q(α,β,γ)|αεC∪βεC∪γεC}

[0053] where

[0054] C={−δ, −δ+2δ/(n−1), −δ+2.2δ/(n−1), . . . , δ} or a set of n angles appropriate to the expected range (±δ) of the motion of the subject, and

[0055] q is the quaternion by which the navigator is rotated, defined by:

[0056] q(α,β,γ)=q_α,RR·q_{β, RA}·q_γ,RSand

[0057] q(θ,φ,Φ)=q_αθ,RRA·q_βφ,RAS·q_γΦ,RRS

[0058] where

[0059] RR=[+1, 0, 0],

[0060] RA=[0, +1, 0],

[0061] RS=[0, 0, +1],

[0062] RRA=(1/2)·[+1, −1, 0],

[0063] RAS=(1/2)·[0, +1, −1], and

[0064] RRS=(1/2)·[−1, 0, +1],

[0065] where q_θ,Ris a quaternion representing a rotation of θ about the vector R in the RAS coordinate system.

[0066] In particular, it is possible to collect the octant navigators pre-turned about the optimal angle Θ and at the optimal radius ρ using the following:

[0067] Δα=−2Δα, −Δα, 0, Δα, 2Δα,

[0068] Δβ=−2Δβ, −Δβ, 0, Δβ, 2Δβ, and

[0069] Δγ=−2Δγ, −Δγ, 0, Δγ, 2Δγ,

[0070] thus providing 5×5×5=125 combinations for each of, e.g., three maps, with the optimal angles having the values of 10°, 3° and 1°.

[0071] It is also possible to select q+q+q (e.g., 31+31+31) combinations of the perturbations PR, where each angle is perturbated separately, as follows:

[0072] 31 steps of Δα, with Δβ=Δγ=0,

[0073] 31 steps of Δβ, with Δα=Δγ=0, and

[0074] 31 steps of Δγ, with Δα=Δβ=0.

[0075] Thus, the q+q+q rotation map is the collection of octant navigators along the path of the octant navigators, and the rotations can be defined as {q(θ,φ,Φ)|(θ,φ,Φ)ε{(C,0,0)∪(0,C,0)∪(0,0,C)}}, with the variables being defined above in a similar manner as for the n×n×n rotation sequence map. C can be defined as, e.g., a set of 31 angles spanning the extent of the desired map in each direction, and dependent on the expected maximum inter-navigator rotation.

[0076] Thereafter, a post-processing procedure can be implemented to interpolate the q+q+q rotation sequence map. By using this map, the data for the arbitrary combinations of angles is preferably calculated from the initial sequence map. An interpolation of scalar values on the surface of a sphere 130 should be used, which can be easily achieved by a bilinear or quadratic interpolation using quaternions to linearize the distances on the spherical surface.

[0077] As an alternative to interpolating the maps, the incoming navigators may be matched to the closest fit in the map(s). By fitting a quadratic error surface to this mean-square error match and its spatially closest neighbors in the map, and finding the minimum of this surface, a refined estimate of the matching angle may be calculated without the need for explicitly interpolating the navigator map(s). More specifically, the best fit rotation angle may be calculated as A_min=−inv(H)*G where H is the Hessian matrix calculated in the region of the best fit in a similar manner to that described above, and G is the gradient (e.g., first partial derivatives of the error with respect to the rotation angles). A_minand G can both be column vectors of 3 elements corresponding to the three angles of rotation. In addition, A_mincan be the refined offset from the best fit.

[0078] FIG. 11 shows an exemplary embodiment of a scan sequence of the navigators to be used for collecting and extrapolating the image as a function of the data obtained from comparison of the pre-mapped data and the actual navigator data (step 120 of FIG. 2). This drawing shows the scan sequence of the conventional imaging kernel with embedded navigators in each block, preferably between the image line readout and the spoilers, arranged in a predetermined number of partitions, each having a particular number of lines which are provided via Cartesian coordinates. The scan which utilizes these kernels is preferably the only scan where the first read-out of the ADC is activated. In an exemplary embodiment of the present invention, there may be two reads per block. The first read is provided for the line, and the second read is provided for the navigator. The first read not being activated for the map or prescan. Indeed, a line (at any one time) of p k-space samples is scanned, read or acquired using the sequence. Thus, the entire scan of n×m lines results in an n×m×p voxel k-space volume.

[0079] Then, the k-space volume (m×n×p) is transformed into an image space volume (m×n×p) using, e.g., the Fourier Transform. In addition, the octant navigator is acquired for each line, and this navigator is compared with the data of the navigator map to determine the magnitudes of the rotations of the object (Δα, Δβ, Δγ). The translations can then be determined by comparing the phase of the octant navigator in question with its best match in the map (Δx, Δy, Δz) using the following equation F [j(x+Δx, y+Δy, z+Δz)=F(u, v, w)e^{j2π(ux+vy+wz)}, with u, v, w being spatial frequencies in each of the three dimensions of k-space. The values are then provided back to the next block that is rotated by (Δα, Δβ, Δγ) before the acquisition when the next line is scanned. In particular, the gradients of each block are rotated according to the octant navigator from the previous block. The translations are corrected by adjusting the phase of the k-space volume data. This can be achieved during the scan or thereafter. In addition, it is possible for all the corrections (including the rotation) to be performed offline after the scan is completed.

[0080] FIG. 12 shows an exemplary embodiment of a portion of the method according to the present invention which utilizes the above-described scan and the scan sequence, and in which the image of the object is transformed and corrected (in real-time) using the navigator. In this exemplary method, the gradients associated with the first scan line of the volume to be scanned are set up in step 500, and the k-space is acquired for such scanned line (step 510). Then, in step 520, the octant navigator associated with such scanned line is acquired. The rotation of the navigator is calculated using the comparison between the actual navigator at the scanned line and the data from the map (step 530), and then the gradient is adjusted using the calculated rotation (step 540).

[0081] Thereafter, in step 550, it is determined whether all lines were scanned or acquired. If not, the gradients (possibly modified/rotated) are set up for the next line to be scanned or acquired (step 560), and the process is returned to step 510 to acquire the next line and begin the calculations for the octant navigator associated with the next line. Otherwise, the phase for the translations are corrected using the previously corrected gradients and the phase errors between the navigators and the map (step 570). Then, in step 580, the k-space volume of the object is transformed into the image, with the rotation and/or motion taken in consideration. It should be understood that the steps 500-580 are preferably performed by the processing arrangement 40 shown in FIG. 1. It is possible for the image to be displayed on the display device 10, stored in the database 20 or printed by the printer 30 (step 590) to be possibly used for further analysis and/or for comparison purposes. In addition, the trace of the motion (e.g., a table/graph of rotations/translations) may similarly be displayed/stored/printed whether the motion was corrected or simply measured.

[0082] Other embodiment of the system and method according to the present invention are conceivable. For example, FIG. 13 shows another exemplary embodiment of the method according to the present invention for tracking and correcting a motion of the anatomical object in real time. In particular, the optimal radius is selected in step 600. This radius selection has been described in great detail above, e.g. in a “pre-mapping radius selection” or in a “concentric navigator” scan. Then, in step 610, the navigator map is created. Two exemplary embodiments of the navigator maps (i.e., n×n×n and q+q+q) have been described above in detail. The octant navigator for the selected scanned line is acquired (step 620), and this octant navigator is matched and/or compared to the data of the acquired map (step 630). Thereafter, in step 640, the motion of the object can be predicted, and in step 650, the respective gradients are corrected using the predicted motion of the object and the results of the calculation of step 630. Furthermore, certain models (e.g., an ARMA model) can be used to forecast the position of the object so that further motion between the time of the navigator and the subsequent correction of the gradients can be accounted for, even if the motion is oscillatory. A feedback may also be achieved using a feedback control system with a state estimator such as a linear-quadratic-Gaussian controller with the state estimated by a Kalman filter.

[0083] In step 660, it is determined whether all image lines were scanned or obtained. If not, the next line is scanned or obtained (step 670), and the process is returned to step 620 to begin the calculations for the octant navigator associated with the next line. Otherwise, in step 680, the image and/or data corresponding to the motion of the object can be displayed on the display device 10. In addition, the image may be stored in the database 20 and/or printed by the printer 30 (step 690). Furthermore, the motion data (e.g., the rotations and/or translations) may be displayed, stored or printed.

[0084] The foregoing exemplary embodiments merely illustrate the principles of the present invention. Various modifications and alterations to the described embodiments will be apparent to those skilled in the art in view of the teachings herein without departing from the spirit and scope of the invention being defined by the appended claims. For example, the above exemplary embodiment have been described with particular reference to and use by the MRI imaging systems and techniques. However, the exemplary of the system and method according to the present invention can be used by other imaging systems and methods, such as CAT scanning systems and methods.