WO1996017261A1

WO1996017261A1 - Method of minimizing dead-periods in order to obtain hardware optimized gradient waveforms

Info

Publication number: WO1996017261A1
Application number: PCT/US1995/011167
Authority: WO
Inventors: Ergin Atalar; Elliot Ross Mcveigh
Original assignee: Johns Hopkins University
Current assignee: Johns Hopkins University
Priority date: 1994-11-25
Filing date: 1995-09-01
Publication date: 1996-06-06
Anticipated expiration: 1997-05-25
Also published as: US5512825A; AU3502495A; EP0809810A1; JPH10509893A; EP0809810A4

Abstract

A method of minimizing dead-periods in magnetic resonance imaging pulse sequences employs a specimen disposed within magnetic field (6), a source of RF signals (2), a receiver (8) for receiving signals from the specimen responsive to RF pulses and emitting respective output signals. A computer (12) is provided for receiving the output signals from the receiver and establishing image information which may be displayed (16). For the dead-period, the minimum and maximum phase encoding step, the scan plane gradient pulse for the slice, phase encoding and readout directions are determined and the moments contained within the dead-period waveform is determined. The values are transformed into gradient amplifier coordinates and the minimum dead-period based on the dead-period waveform is determined. The dead-period is employed in establishing a hardware optimized waveform which may be trapezoidal. The trapezoidal waveform is preferably established by for each phase encoding step determining the starting and ending gradient levels and the moments contained within the dead-period waveform using the calculated minimum dead-period.

Description

METHOD OF MINIMIZING DEAD-PERIODS IN ORDER TO OBTAIN HARDWARE OPTIMIZED GRADIENT WAVEFORMS.

1 Field of the Invention The present invention relates to an improved method for reducing the dead- periods in magnetic resonance imaging pulse sequences and to associated apparatus and. more specifically, it relates to a simplified means tor reducing dead-periods by employing a calculated waveform and. most specifically, is particularly advantageous in imaging oblique planes 2 Description of the Prior Art

The advantageous use ot non-invasive and non-destructive test procedures has long been known in both medicine and industrial applications In respect ot medical uses that limit a patient s exposure to potentially damaging x-ray radiation, it has been known to accomplish this objective through the use of other non-invasive imaging procedures, such as. for example, ultrasound imaging and magnetic resonance imaging

See. for example. United States Patents 4.766.381. 5.099.208 and 5.352.959

In a general sense, magnetic resonance imaging involves providing bursts of radio frequency energy to a specimen positioned m a main magnetic field in order to induce responsive emission of magnetic radiation from the hydrogen nuclei or other nuclei The emitted signal may be detected in such a manner as to provide information as to the intensity of the response and the spatial origin ot the nuclei emitting the responsive magnetic signal. In general, the imaging may be performed in a slice or plane, or multiple planes, or three-dimensional volume with information corresponding to the responsively emitted magnetic radiation being delivered to a computer which stores the information in the form of numbers corresponding to the intensity of the signal The computer establishes a pixel value as by employing Fourier Transformations which convert the signal amplitude as a function of time to signal amplitude as a function or frequency The signals may be stored in the computer and may be delivered with or without enhancement to a video screen display, such as a cathode-ray tube, tor example wherein the image created by the computer output w ill be presented through regions ot contrasting black and white which vary in intensity or color presentations which \ arv in hue and intensitv The challenges involved m producing high resolution, accurate imases in a timely fashion become more complex and difficult to solve when imaging oblique planes Such oblique plane imaging is frequently employed, tor example, in breath-hold magnetic resonance cardiac imaging. Efforts to employ trapezoidal waveforms which have been employed in non-oblique planes in oblique imaging have resulted in the problems ot maximum slew rate and gradient level limitation because the waveforms applied to the gradient amplifiers would not be identical to those applied in non-oblique imaging. The waveforms applied to the amplifiers are a function of the orientation of the slice being imaged. If these non- oblique imaging waveforms were applied in oblique imaging, the gradients would not be able to produce the waveforms and image artifacts and distoπions result.

It has been suggested to solve this problem by avoiding application ot t or more gradient waveforms simultaneously. See Bernstein et al.. "Pulse Sequence Generated Oblique Magnetic Resonance Imaging: Application to Cardiac Imaging. " Med Phys.. Vol. 13. pp. 648-657 (1986), Erratum: Med. Phys. 14(1): 145. ( 1987)

It has also been suggested to decrease the slew rates and maximum allowed gradient levels in designing trapezoidal waveforms by a factor of the 3 The problem with this approach is that it tends to cause an unnecessary increase in the minimum possible TE and TR. Another approach, which also employs trapezoidal waveforms. calculates the slew rates and maximum gradient levels depending on the orientation ot the slice and the number of active gradients See Bernstein et al . "Angle-Dependent Utilization of Gradient Hardware for Oblique MRI." J Mag. Reson. Imag.. Vol 4. pp 105-108 (Jan. 1994). This last method will be referred to herein as the scan plane optimization" method. Another approach would be to employ gradient amplifiers and coils that allow higher gradient levels and higher slew-rates. This, of course, would involve the expense of acquiring new hardware.

In spite of the foregoing known systems, there remains, a ver> real and substantial need for a method of magnetic resonance imaging which is capable ot producing faster, high quality magnetic resonance images on an oblique plane in reduced data acquisition time. SUMMARY OF THE INVENTION

The method of the present invention permits employing conventional magnetic resonance imaging equipment which has the computer controlling operation ot the gradient coils and amplifiers programmed, either by appropriate software or hardware embedded means. In order to effect a reduction in TE and TR. in a preferred practice of the method, the duration of each dead-period is reduced. As employed herein, the term "dead-period" in a gradient-echo pulse sequence is the time interval between RF pulses and readout. A dead-period is a time interval in an MR pulse sequences where the waveform shape is not important, but its zeroeth. and/or higher moments are important. The nth moment of a gradient waveform in a dead-period is defined as

m -f '' g(t)t "ι dt ( 1)

where t, and t_c are the staπing and ending time of the waveform g.

As employed herein, the staπing and ending values and zeroeth and/or higher moments of the slice, readout and phase encoding waveforms of a dead-period will be referred to as "the parameters of the dead-period " The zeroeth moment ill be the dead-period waveform area.

In a preferred practice, for the maximum and minimum phase encoding step, the parameters of the dead-period are determined. These values are transformed into gradient amplifier coordinates and the performance characteristics of the amplifier are used to compute an optimal waveform trom the dead-period.

In the preferred practice ot the invention, the method involves tor each phase encoding step for the slice, phase encoding and readout directions, the staπing and ending gradient levels and the moments contained within the waveform are determined and are transformed into gradient amplifier coordinates The minimum dead-period is then employed with this calculated information to design a hardware optimized waveform for the desired imaging sequence. This approach is paπicularly impoπant in respect ot oblique scans.

The apparatus of the present invention provides magnetic field generating means tor establishing a main magnetic field on the specimen, magnetic field gradient generating means for establishing gradients in said main magnetic field, and RF signal generating means tor emitting pulsed RF signals to at least a poπion ot the magnetic field passing through the specimen Dead-period determining means to determine the minimum and maximum phase encoding steps by delivering for the slice, phase encoding and readout directions the scan plane pulse staπing and ending gradient levels and the zeroeth moment or zeroeth and higher moments contained within the dead-period waveform and transforming these values into gradient amplifier coordinates as by rotation after which the reduced or minimum dead-periods are calculated based upon performance characteristics of the amplifiers. While the preferred practice of the invention emplov s the minimum dead-period, reductions in the dead-period which do not reach as low as the minimum may be advantageously employed.

The design waveform, which is employable in oblique and non-oblique imaging, is obtained by design waveform creating means which for each phase encoding step for the slice, phase encoding and readout directions calculates the staπing and ending gradient levels and the gradient moments contained therein. These v alues are transformed to gradient amplifier coordinates by rotating the same A hardware optimized waveform is then determined employing the calculated minimum dead-periods The apparatus also provides receiver means for receiving signals emitted trom the specimen responsive to the RF pulse being emitted and emitting responsive output signals which are delivered to computer means for establishing image information related thereto and visual display means tor displaying images from the image information received trom the computer means.

It is an object of the present invention to provide a method and related apparatus for providing taster magnetic resonance imaging through reduction in the TE and TR values, without reduction in the signal to noise ratio.

It is a fuπher object of the present invention to provide such a stem which may be employed effectively in respect of oblique imaging.

It is another object of the present invention to provide such a sv stem hich is adapted for use with existing apparatus and will avoid the undesired need to pro . ide larger capacity amplifiers and coils which would permit higher gradient levels and higher slew rates. It is a ruπher object of the present invention to facilitate determination of the minimum dead-period times and an effective hardware optimized waveform in gradient amplifier coordinates by focusing on the staπing amplitude, the ending amplitude and gradient moments of the dead-period. It is a fuπher object of the present invention to provide such a system which facilitates determination of the minimum TE and TR for a given set oϊ gradient amplifiers and coils.

It is another object of the present invention to provide improved image quality and more rapid imaging for any magnetic resonance pulse sequence wherein a dead-period exists.

It is a fuπher object of this invention to provide such a system which reduces and. in the preferred form, minimizes TR and TE w ithout altering the slice profile or increasing receiver bandwidth.

It is a fuπher object of the present invention to provide such a system which is employable in breath-hold cardiac imaging.

It is yet another object of the invention to provide such a system wherein the duration of dead-periods may be reduced without altering the receiver bandwidth or the slice profile.

These and other objects of the invention will be more fully understood from the following detailed description of the invention on reference to the illustrations appended hereto.

BRIEF DESCRIPTION OF THE DRAWINGS Figure 1 is a schematic illustration of a magnetic resonance imaging system. Figure 2 is an illustration of a gradient-echo pulse sequence optimized for imaging non-oblique planes of the prior an.

Figure 3 is a plot of the waveforms corresponding to Figure 2 as applied to gradient amplifiers/coordinates.

Figure 4 is a plot of scan gradient values employed in the present invention and Figure 5 is a representation of the waveforms of Figure 4 after transrormation to gradient amplifier coordinates. Figure 6 illustrates hardware optimized trapezoids tor the gradient-ech pulse sequence showing amplifier currents and Figure 7 illustrates a corresponding gradient-echo pulse sequence which has been optimized

Figures 8 through 10 are logic diagrams showing a preferred practice o the present invention

Figures 11 and 12 are contour plots of TR variation as a function ot plan orientation for. respectively, the practice of the present invention and prior an scan plan optimization

Figure 13 is a plot of scan gradient coordinate values employing gradien moments and Figure 14 is a plot ot the waveforms of Figure 13 after conversion t gradient amplifier coordinates

Figure 15 is an array ot eight types ot hardware optimized trapezoida pulses employable in the determination of the minimum dead-period

Figure 16 is an array of waveforms of hardware optimized trapezoida pulses employable in determining waveforms for the amplifiers to be played duπng th minimum dead-period

Figure 17 is a plot of an optimized flow-compensated gradient-echo puls sequence employing hardware optimized trapezoidal pulses in the amplifier coordinat system Figure 18 is a flo -compensated gradient-echo pulse sequence optimize tor a non-oblique plane

DESCRIPTION OF THE PREFERRED EMBODIMENTS As used herein, the term "specimen ' refers to any object placed in th mam magnetic field for imaging and shall expressly include, but not be limited t members of the animal kingdom, including humans, biological tissue samples, and tes specimens removed from such members of the animal kingdom It shall includ inanimate objects which may be imaged by magnetic resonance or which contain wate or sources of other sensitive nuclei

Figure 1 shows a schematic representation of the general concept o magnetic resonance imaging An RF source 2 provides a pulse radio trequencv energ to the specimen which, in the form shown, is a patient 4 in the main magnetic field whic is created by magnetic field generator 6 The specimen is generallv aligned w ith th main magnetic field and the RF pulses are imposed perpendicular thereto. Where oblique imaging is employed, the angle of impingement of the RF vector representing the spatial gradient in the magnetic field will be angularly offset from either the x. \ . or z directions. This arrangement results in excitation of the nuclei within the area or volume to be imaged and causes responsive emission of magnetic energy which is picked up bv receiver 8.

The receiver 8 may be a coil which has a voltage induced in it as a result of such responsive emissions of magnetic energy. As a practical matter, separate coils or identical coils may be employed as the RF source 2 and the receiver 8. The signal emerging from receiver 8 passes through analog-to-digital (A/D) convener 10 and enters computer 12. Within the computer 12 the Fourier Transformations of signals conveπ the plot of amplitude versus time to a map of the distribution of frequencies by plotting amplitude versus frequency. The Fourier Transformations are performed in order to establish the intensity values and locations of specific pixels. These values may be stored, enhanced or otherwise processed and emerge to be displayed on a suitable screen, such as a cathode-ray rube 16. for example.

Referring to the prior an representations of Figures 2 and 3. the gradient- echo pulse sequence includes the RF pulse sequence 24. the slice pulse 26. the phase encoding pulse sequence 28 and the readout pulse sequence 30. It will be appreciated by looking at the legends at the uppermost poπion of Figure 2 that the regions during which an RF pulse is being generated are labeled with the expression "RF" and the period for readout of data is indicated by the legend "READ. " The intermediate dead-periods are labeled, respectively. Dl and D2. The present invention targets reduction in the D l and D2 time periods which are the TE periods and the total elapsed time in the cycle of data acquisition which is the TR period. In this manner, the present invention permits the use of existing equipment at a paπicular facility with a customized designed waveform which facilitates reduction in the elapsed time periods Dl. D2. in which the waveform shape is not impoπant for image formation, but the zeroeth. first and higher moments of the waveform are. Such reduction of time is paπicularly critical in respect of oblique imaging and imaging of objects which are in motion such as cardiac imaging, for example. Refemng to Figure 3. there are shown the RF pulses 24' and waveforms of Figure 2 as convened to gradient amplifier coordinate systems for application to _the gradient amplifier. For purposes of Figures 2 and 3. the imaging plane has been ob_tained by respective rotations of 30°, 30° and 45° rotations about the X. Y" . Z" axis. An example of a suitable 30°. 30°. 45° rotation matrix is shown in Equation 2.

¹ 0.612 0.612 0.500

R = -0.789 0.436 0.422 (2)

0.047 -0.660 0.750

Figure 3 is illustrative of the problems created as the darkened regions ot the respective gradient waveforms, as shown in pulse sequences 26^* . 28' and 30' . are regions wherein the waveform exceeds the maximum allowable gradient level and slew rate. As a result, effoπs to employ the pulses of Figure 2 applied to the gradient amplifiers when used in imaging oblique planes, are inadequate because the waveform shown in Figure 2 are no longer applied to the gradient amplifier. The application of the waveforms of Figure 3 to the gradient amplifiers which will vary depending upon the orientation of the oblique slice, will produce undesired image aπifacts and distoπions In this approach the waveforms that are applied to the amplifiers are actually the weighted sum of the slice, readout and phase encoding waveforms wherein the weights are a function of slice orientations. It is this type of problem with oblique imaging that the present invention is directed toward primarily, although, the solutions for oblique imaging are such that the present invention^'s techniques and associated apparatus mav advantageously be employed in non-oblique imaging.

The present invention provides a system for calculating reduced and preferably the minimum duration dead-periods for scanning a given plane, including oblique planes. This is accomplished by focusing on the requirements of the gradient pulses for the paπicular scan plane during the dead-period based upon the starting gradient level, the ending gradient level and the moments contained within the waveform in the dead-period of the slice, phase encoding and readout directions. These values are transformed from the scan plane coordinates to amplifier coordinates with an appropriate ro_ta_tion by an oblique rotation matrix. Within the amplifier coordinate system, the hardware optimized pulse for each gradient amplifier for a given staπing amplitude, ending amplitude, and moments with a maximuη. slew rate waveform. This, in effect, is a waveform which has been optimized based upon the hardware considerations. Referring in greater detail to Figures 4 and 5. the dead-period Dl has waveforms for slice direction 40, phase encoding direction 42 and readout direction 44 By determining the staπing gradient level s_s, s_p. s_r. respectively, and the ending gradient levels e,, e_p, e_r, as well as the zeroeth moment, i.e.. area included within the waveform between the respective staπing gradient levels and ending gradient levels, the minimum dead-periods based upon trapezoidal waveforms may be determined. In determining the zeroeth moment, the region under the waveform within the dead-period defined by Equation (1 ) with n=0 It is impoπant to note that the present invention is useful as to all magnetic resonance imaging wherein the shape of the waveform within the dead- period is not of importance. This will include essentially all magnetic resonance imaging except spectroscopy not employing gradients. The invention uses of the three parameters, i.e. , the staπing gradient value, the ending gradient value and the moments, the waveform shape is not a required parameter. With respect to the moments, the zeroeth moment will be employed alone for some pulse sequences and with higher moments tor other pulse sequences. As the hardware limitations imposed upon the process and associated method are related to the slew rate and maximum current that an amplifier can achieve, in the preferred embodiment action is taken to conveπ the values shown in Figure 4 to amplifier current coordinates from the imaging coordinates. The manner in which this transformation is achieved will now be considered. The present invention employs a 3x3 oπhogonal rotation matrix R that translates from the slice encoding waveform, the readout encoding waveform and the phase encoding waveform to the corresponding physical gradient amplifier coordinates This transformation can be accomplished by multiplication of the parameter vectors w ith the rotation matrix. R, in accordance with Equation 3.

R is a matrix determined by the oblique angle being employed. In Equation 3 lower case "s" and "e." respectively, represent the starting an ending gradient values and lower case "a" represents- the zeroeth or higher moment Subscripts "r." "p" and "s" represent, respectively, the readout, phase encoding and slic directions and lower case "x. y. z" represent the amplifier current directions. In th way. for each amplifier current direction, the parameters of staπing and ending gradie values and moments can be determined.

The waveform of desired reduction or the shoπest duration that satisfi the three parameters for the amplifier current directions to produce the reduced shoπest waveform satisfying these parameters is referred to as the "hardware optimize waveform. " As used herein, the term shall be deemed to refer to trapezoidal waveform triangular waveforms with maximum slew rate ramps or other waveforms designed t utilize specific hardware optimally. The equation for the moment covered by t hardware optimized waveform with arbitrary staπing and ending levels is a first second order equation. The resultant amplifier coordinate waveforms are shown i Figure 5 wherein for the "x.y.z" directions, respectively, s, and e, are the staπing an ending gradient values and "a," is the zeroeth moment and for the ^"y" direction, the and e, beginning and ending gradient values with zeroeth moment a and for the "z direction. s₂ and e_z with the zeroeth moment being __τ.

It is desired to keep the dead-period the same for each of the amplifie and for each phase encoding step. The calculation time for the global minimum dea period may be reduced by considering only the maximum and minimum phase encodi steps. As the required area for each amplifier is a linear function of the phase encodi step, either the maximum phase encoding step or the minimum phase coding step wi give the longest of the minimum dead-periods. In this manner, by selecting the longe of the minimum dead-periods, which will become the global minimum, each amplifi will have the same dead-period and the amplifier^'s performance limits ill not exceeded. After determining the global minimum for each dead-period, a hardware optimized waveform pulse that satisfies the three parameters, i.e.. the staπing and ending current values and the moment and fits into this minimum time is designed.

While the preferred embodiment of the present invention involves seeking the minimum dead-periods, it will be appreciated that in some instances, one may obtain benefits of the invention by employing a reduced dead-period which may be greater than the absolute minimum. Such reduced dead-periods, when attained through the use of the methods or apparatus of the present invention, shall be deemed "minimuu " for purposes of the invention. In such an approach, the targeted or desired degree of dead-period reduction will be the "minimum" in the practice of the invention.

In order to achieve extremely shoπ TR and TE values, high duty cycle gradient amplifiers should be employed. In order to decrease power consumption wherever possible, the gradient level is decreased and the duration of the hardware optimized waveform is increased. As the waveforms that will be applied to the amplifiers are designed, the power required from the amplifiers can be evaluated by integrating the square of the waveform over time and TR can be increased if necessary. An example of a design is shown in Figure 6 which illustrates the hardware optimized trapezoids for the gradient-echo pulse sequence. These waveforms represent the amplifier currents for the respective x, y, z axes and show the corresponding RF pulses. These were obtained through rotation of the oblique matrix by 30°. 30° and 45° rotations about the x. y^' and z" axis under the influence of rotation matrix R. These waveforms can be applied directly to the gradient amplifiers. In the alternative, the hardware optimized trapezoid can be transformed back to the readout, phase encoding and slice coordinate system by multiplying the gradient waveform vector by the inverse of the oblique rotation matrix R ' . This waveform, after transformation, is shown in Figure 7. It will be noted that the sequences have two dead-periods. Dl and D2. The waveform of Figure 7 can be employed directly to apply the calculated waveform on the scanner. It is noted that although the waveforms applied to the amplifier, as shown in Figure 6. are simple trapezoids. the slice, readout and phase encoding waveforms, as shown in Figure 7. are somewhat complicated as optimizing in the scan plane coordinates is complex. In conventional designs, as shown in Figures 2 and 3. the reverse is true in that the slice, readout and phase encoding waveforms are simple and the waveforms that are applied to the gradient amplifiers are complex To enhance understanding of this embodiment ot the invention, an example will be considered. Employing hardware optimized wave_forms in the form ot trapezoidal pulses. TR values in the range of 7 milliseconds on oblique planes for the General Electric Signa Version 5.2 scanner that has isotropic gradients with a maximum gradient of 1 Gauss/cm and a slew rate of 0.167 T/sec/cm has been achieved. The minimum possible TR with this pulse sequence depends on the angle ot the oblique plane and the specific imaging parameters, such as slice thickness. FOV (field of view), receiver bandwidth and number of phase encoding steps A shoπ axis image of the heart with a double oblique view was acquired using segmented k-space cardiac tagging as disclosed in McVeigh "Cardiac Tagging with Breath Hold CINE MRI. " Magn

Reson. Med.. Vol. 28, pp. 318-327, ( 1992). This image was obtained employing the following parameters- 5 mm slice thickness. 250 Hz/pixel receiver band. 360 mm FOV matrix size 256 x 128: paπial echo and paπial excitation, i.e.. 1 sine side lobe before and no side lobe after the main lobe of the sine RF pulse. Referring to Figure 8. the logic diagram shows a preferred practice of the present invention wherein the initial stage 70 involves calculating the slice and readout gradient levels along with the phase encoding area. For each dead-period, this information is processed by calculating the minimum dead-period 72 after which the longest of the minimum dead-periods is employed to design the hardware optimized waveform 74 which is applied to the gradient amplifiers This process is employed tor dead-period 76 as indicated by the arrow returning to block 72.

With reference to Figure 9. greater details regarding the determination ot minimum dead-period will be considered. For the slice, phase encoding and readout directions 80. the maximum and minimum phase encoding step 82 is determined b> calculating the staπing and ending gradient levels and the required area 84 Focusing on the maximum and minimum phase encoding step saves time. The next determination is whether the first or higher order moments are impoπant 86 These would be impoπant. for example, in the event the system were employed for velocity-encoded or flow- compensated computations. While these two examples employ the zeroeth and first order moments, other sequence uses will use the zeroeth moment (area) and/ or second or higher moments, for example The required moments are calculated 88 and the results 84. 88 directed toward the transformation of the values to gradient amplifier coordinates 90 It the moments are not impoπant. the output of the response to 84 is delivered direc_tlv _to the transformation to gradient amplifier coordinates poπion 90. For each gradien_t amplifier 92. there is calculated the minimum dead-period based on i_ts hardware limitation 94. The process is repeated for the next phase encoding step 96 and is fuπher repeated for each direction 98 of gradient amplifier. The longest of the minimum dead- periods 100 is then determined. This will be employed in the design of the hardware optimized waveform as shown in Figure 10.

As shown in Figure 10. for each phase encoding step 110. for the slice. phase encoding and readout direction 112. the staπing and ending gradient levels, along with the required moments 114 are determined. If the first or higher order moments are impoπant 113. the moments are calculated 117 and delivered along with the output of 1 16 to a transformation block 118. Such moments are impoπant. for example, in respect ot velocity-encoding or flow-compensation systems. If the moments are not of consequence, the output of 116 is delivered directly to 118. wherein the values are transformed to gradient amplifier coordinates through multiplication by the rotation matrix. For each gradient amplifier coordinate 120. the hardware optimized waveforms employing the calculated dead-periods design hardware optimized waveforms 122. For the next direction of gradient amplifier 124. the process is repeated. For each phase encoding step 126. the process is repeated.

Within the parameters of this example, the minimum TR and TE for the imaging parameters were determined. The minimum TR and TE values were calculated for all possible oblique orientations and phase and frequency encoding directions with a resolution of 5 ° . The results are summarized in Table 1 with the results being expressed in milliseconds.

TABLE 1

Hardware Optimized Waveform Scan Plane Opt.

Min Max Mean Std. Dev, Min Max Mean Std. Dev

TR 6.516 6.964 6.769 0.112 6.964 8.554 7.945 0.463

TE 2.588 2.952 2.744 0.101 2.952 3.648 3.355 0. 189 The TR and TE values were also calculated for the prior an scan plane optimization method which was described in the Bernstein et al. 1994 anicle cited hereinbefore. The hardware optimized waveform approach of the present invention in respect of minimum TR and TE established a reduction of up to 30 and 40 percent. respectively, in TR and TE as compared with the prior an scan plane optimization.

The TR and TE values as a function of slice orientation are shown in the contour plots of Figures 11 and 12. In these figures, the plane orientation is given by azimuth and polar angles of the plane normal. The frequency in phase directions of the RF pulse width that gave the minimum TR were employed. In comparing the hardware optimized waveforms, which in this embodiment were trapezoids of the present invention, of Figure 11 with the prior an scan plane optimization as shown in Figure 12. the TR variation for the former is less than for the latter. With hardware optimized trapezoid pulses of Figure 11. the TR becomes smaller as the position of the slice becomes oblique By contrast the scan plane optimization, as shown in Figure 12. causes an increase in TR as the slice becomes oblique. This behavior also was observed for TE with the pulses.

The highest TE value of the present invention was observed for non-oblique planes, whereas the scan plane optimization method of the prior an had the lowest TE at the non- oblique planes. For non-oblique imaging, hardware optimized trapezoidal pulses give the same minimum TR and TE as that obtained with scan plane optimization. As disclosed hereinbefore, the moment parameter of the dead-period waveform may be the zeroeth moment (area), the zeroeth moment and higher moments depending upon the type of pulse sequences involved. For example, with flow- compensated and velocity -encoded gradient-echo pulse sequences and similar pulses. consideration is given to the first moments. The problem solved by the present invention is essentially the same as that previously discussed, wherein for some oblique slices being imaged, the waveforms applied to the gradient amplifiers may total on a single gradient amplifier in excess of the slew rate and maximum amplitude specifications for the amplifier, such as was shown in Figures 2 and 3. The hardware optimized waveform pulses employed to create amplifier current waveforms in accordance with the present invention, may be employed in connection with a wide variety of pulse sequences, such as for flow-compensated or velocity-encoded gradient-echo pulse sequences, for example As shown in Figure 13. in addition to considering as to each the readout, phase encoding and slice directions, the staπing and ending gradient levels in addition to the areas a_r, a_, a_v (zeroeth moment), the first moments of the gradient waveforms, which have been designated, respectively, m.. m.. m, are considered These parameters (starting and ending gradients, zeroeth moments and first moments) are then convened to gradient amplifier current coordinates, such as are shown in Figure 14 by multiplying the parameter vectors by the appropriate rotation matrix. In addition to employing the convened zeroeth moments a,. a_v. a., the first moments m₅. m_v, m. are employed The shape of the dead-period waveform is not impoπant to this determination. The amplifier current coordinates are shown with the moments being respectively, m,. m,. and m. The next step is to calculate the minimum possible dead-period from the waveform that has the highest demand on one of the amplifiers. The solution for the minimum dead-period is simplified by the fact that the waveforms use tor this pulse will be maximum slew rate ramps or maximum amplitude constants. In a preferred approach, in determining the minimum dead-period, eight hardware optimized trapezoidal waveform types, such as that shown in Figure 15. will be considered with one of them being selected tor a given set of S-. s,. s_z. e,. e,, e_z, a,, a_^, a., and m„, m , m. parameters Depending upon the waveform type, a second or fouπh order polynomial equation will be solved to select the specific waveform in a manner known to those skilled in the an.

Equal duration constant waveforms are added to each plateau in the bipolar pulses of the hardware optimized trapezoidal waveforms obtained in the determination ot

Figure 15 to extend the pulse duration to the required period

In desigmng the waveforms to be employed during this minimum dead- period for this embodiment, the eight forms of hardware optimized trapezoidal avetorms shown in Figure 16 will be considered with the waveform amplitudes and durations being calculated by solving a second order equation.

In this manner, in the example represented by Figures 13 through 16. the parameters of Figure 14 and the waveform arrays of Figures 15 and 16 are emploved in establishing the hardware optimized waveforms

An example for a flow-compensated gradient-echo pulse sequence is show n in Figure 17 This design is the hardware optimized design for an oblique plane that was obtained bv 45° rotations about the x. y' and z^" axes for a pulse sequence ot the general type shown in Figure 18. A rotation matrix suitable for 45° , 45° . 45° rotation is show in Equation 4.

0.500 0.500 0.707 j

R = 1 -0.854 0.146 0.5001 (4)

0.146 -0.854 0.5001

If the pulse sequence of Figure 18 were used for this specific oblique plane, the slew rat maximum amplitude limitations of the gradient amplifiers would be violated.

In order to confirm the operability of the invention in the context of th tlow-compensated or velocity -encoded embodiments, an example will be considered. Th minimum achievable TE for the imaging conditions used for cardiac imaging wer considered with a slice thickness of 5 cm FOV (field of view) of 280 mm: paπial ech panial excitation: a pixel band with a 250 Hz/pixel: RF pulse duration of 2.5 msec. I the optimization were carried out for a non-oblique plane, the minimum TE would b 4.63 msec which, as shown in Table 2. equals the minimum value for oblique imagin by the prior an scan plane optimization method.

TABLE 2

Hardware Optimized Trapezoidal Scan Plane Opt.

Mm Max Mean Dev. Min Max Mean Dev

3.93 4.78 4.46 . 0.21 4.63 6.48 6.02 0.39

An alternate approach would be to calculate the maximum gradient lev and slew rate reduction factor based on the scan plane. The statistics obtained employe using that method are shown in Table 2.

It will be appreciated, therefore, in respect of flow-compensation an velocity-encoding, the present invention offers a substantial reduction in TE by a no iterative method which does not require user interaction. Both anisotropic and isotropi gradient sets can be utilized optimally .

While it will be well known to those skilled in the an how to determin the staπing and ending gradient value and the moments for the waveform in the dea period, a specific example of computation of the same will be considered with a gradient- echo pulse sequence. A "sine" RF pulse with one left-side lobe and no right-side lobe is employed. Assuming the following parameters. Equations 5 through 7 will determine the slice gradient value, relationships 8 through 10 determine the readout direction gradient value, and Equation 11 through 13 determine the three values for phase encoding direction. For slice gradient:

1

O) γthick t_RF β

0 lb) m. = γthick

(7)

For readout direction:

^S _P=° (8ι

*_p = ° i (10) ^m _P γFOV_v

For phase encoding direction:

^Sr υ (11⁾

1

^Cr = (12) γFOV_A

XRES_L m. = (13) γFOV wherein.

T- Magnetogyric ratio thick Slice thickness t_RF Rf pulse duration r phase encoding step

_. YRES- l _ __ YRES-l

FOV, Field of view along the readout direction

FOV Field of view along the phase encoding direction t . Sampling rate of the MR signal

YRES. Resolution along the Y direction XRES: Resolution along the X direction XRES_L. Number of samples of the MR signal before echo

As the invention is adapted for use with conventional hardware, althoug if desired, it may be used employed with specialized hardware employable with magneti resonance imaging, disclosure of the same need not be provided herein in detail. Whil reference to visual display means has been made herein so as to suggest "real-time display, it will be appreciated that visual display can be accomplished by storing th image information and subsequently displaying the same on film and hard copies, to example. Such storage shall be deemed to be embraced with the reference to visu display means as an indirect means of getting information from the computer means t the visual display means.

The reduction in dead-period duration is reduced by the practices ot th invention without reduction in signal to noise ratio. It will be appreciate, theretore. that the present invention provides a improved method and associated apparatus for enhancing or maintaining image qualit of magnetic resonance images while reducing the length of the uat acquisition perio All of this is accomplished in a manner which focuses

three specific variabl within the dead-period and. in general, is independent ot the shape of the wavefor within the dead-period. The invention employs the staπing amplitude, the endin amplitude, and gradient moments of the dead-period to establish a reduced or minimum dead-period. The reduced or minimum dead-periods are transformed to gradient amplifier coordinates and are employed in combination with the staπing and ending amplitudes and the gradient moments to create a hardware optimized waveform, which may be a trapezoidal waveform. The hardware optimized waveform for both oblique and non-oblique imaging is provided. The invention is also employable with velocity-encoded or flow-compensated pulse sequences, for example. All of this is accomplished in a manner which facilitates employing the gradient amplifiers that exist at a paπicular facility to maximum advantage. Whereas paπicular embodiments of the present invention have been described above for purposes of illustration, it will be appreciated by those skilled in the an that numerous variations of the details may be made without depaning from the invention as described in the appended claims.

Claims

- 20 -

WE CLAIM-

1 A method ot minimizing the duration of dead-periods in magnetic resonance imaging pulse sequences comprising providing a specimen disposed within a main magnetic field, a source of pulsed RF signals, receiver means for receiving signals trom said specimen responsive to said RF pulses and emitting responsive output signals, computer means tor receiving said output signals from said receiver means and establishing image information related thereto and display means for displaying images from said imaging information. determining for said dead-periods for the scan plane gradient pulse for the slice, phase encoding and readout directions determining the staπing and ending gradient values and moments of the dead-period waveform, transferring said values and moments to gradient amplifier coordinates. determining minimum dead-periods from said gradient amplifier values, and employing said minimum dead-periods to establish a hardware optimized waveform based upon said minimum dead-period values

2 The method of claim 1 including employing said scan plane determination to determine the maximum and minimum phase encoding steps tor said dead-period

3 The method of claim 2 including determining the longest of said minimum dead-pertods. and employing the longest said minimum dead-period to establish said optimized waveform.

4 The method of claim 3 including effecting transformation of said scan plane minimum dead-period values to gradient amplifier coordinates by multiplication by an oblique rotation matrix

5 The method of claim 2 including establishing said waveform by for each phase encoding step determining the staπing and ending gradient levels and the zeroeth moments containe within the dead-period waveform, transforming said values into gradient amplifier coordinates, and 67

- 21 - employing said transformed gradient amplifier coordinates in determining said hardware optimized waveform.

6. The method of claim 5 including employing said method to image an oblique plane.

7. The method of claim 6 including effecting said minimum dead-period determination without regard to the shape of said dead-period waveform.

8. The method of claim 5 including effecting said hardware optimized waveform determination without regard to the shape of said dead-period waveform.

9. The method of claim 5 including employing each said hardware optimized waveform respectively on said gradient generating amplifiers.

10. The method of claim 5 including employing said method while remaining within the amplifier gradient level and slew rate limitations of said gradient amplifiers.

11. The method of claim 5 including employing said method in breath-hold cardiac imaging.

12. The method of claim 5 including employing in minimum dead-period determination in addition to said scan plane gradient pulse, staπing and ending gradient values and said zeroeth moment for said slice, phase encoding and readout directions, at least one higher moment of the dead-period gradient waveforms.

13. The method of claim 12 including employing a trapezoidal waveform as said hardware optimized waveform.

14. The method of claim 13 including in employing said minimum dead-periods to establish a trapezoidal waveform employing for each phase encoding step for slice, phase encoding and readout directions staπing and ending gradient values, and the zeroeth and first moments of the dead-period gradient waveforms.

15 The method of claim 14 including employing said method in

at least one pulse sequence selected trom the group consisting of velocity -encoded pulse sequences and flow- compensated pulse sequences

16 The method of claim 12 including in determining said dead-period after transformation ot said scan plane coordinates to gradient amplifier coordinates employing hardware limitations in determining said minimum dead-period.

17 The method of claim 6 including after determining said hardware optimized waveform transtorming said hardware optimized waveform to scan plane coordinates

18 The method of claim 17 including employing a trapezoidal waveform as said hardware optimized waveform

19 The method of claim 18 including employing the inverse of a rotation matrix in effecting said conversion to scan plane coordinates

20 The method of claim 6 including employing said method on a plurality ot dead-periods within a gradient-echo pulse sequence

21 Apparatus for minimizing dead-periods in magnetic resonance imaging pulse sequences comprising magnetic field generator means for establishing a main magnetic field on a specimen, magnetic field gradient generating means for establishing gradients in said magnetic field.

RF signal generating means for emitting pulsed RF signals to at least a pomon of said magnetic field passing through said specimen. dead-period determining means for determining phase encoding steps by determining for the slice, phase encoding and readout directions, the pulse staπing and ending gradient values and the moments contained within the dead-period waveform. means for transforming these scan plane dead-period values and moments to gradient amplifier coordinates. design waveform creating means for employing a dead-period determined by said dead-period determining means and establishing a hardware optimized waveform based thereon, receiver means for receiving signals emitted from said specimen responsive to said RF pulses and emitting responsive output signals. computer means for receiving said output signals from said receiver means and establishing image information related thereto, and visual display means for displaying images from said image information received from said computer means.

22. The apparatus of claim 21 including said dead-period determining means having means for determining the minimum and maximum said phase encoding steps.

23. The apparatus of claim 22 including said dead-period determining means having means for determining the longest said minimum dead-periods, and said design waveform creating means employing said longest minimum dead-period as said dead-period.

24. The apparatus of claim 23 including said means for transforming including multiplication by a rotating matrix.

25. The apparatus of claim 20 including said dead-period determining means having means for effecting said dead-period determination in oblique plane magnetic resonance imaging.

26. The apparatus of claim 24 including said rotating matrix being an oblique plane rotating matrix.

27. The apparatus of claim 22 including said design waveform creating means having means for establishing said design waveform for each phase encoding step for the slice, phase encoding and readout direction by determining the staπing and ending gradient values and the moments contained within said dead-period waveform and transforming said dead-period values and moments into gradient amplifier coordinates and creating said hardware optimize waveform therefrom.

28. The apparatus of claim 21 including said dead-period determining means and said design wavefor creating means having means for effecting said determinations with respect to an obliqu plane of said specimen.

29. The apparatus of claim 28 including said dead-period determining means effecting said determinatio without regard to the shape of said dead-period waveform.

30. The apparatus of claim 29 including said design waveform creating means effecting said hardwa optimized waveform determination without regard to the shape oi said dead-peri waveform.

31. The apparatus of claim 30 including said apparatus having means for employing each said hardwa optimized waveform respectively on each said gradient generating amplifiers.

32. The apparatus of claim 30 including said design waveform creating means creating said trapezoid waveform which will remain within the gradient levels and slew rates of said amplifier

33. The apparatus of claim 21 including said dead-period determining means having means for determini said phase encoding steps by determining a higher moment of said dead-period wavefor for said slice, phase encoding and readout directions in addition to determining sa staπing and ending gradient values and the zeroeth moment.

34. The apparatus of claim 27 including said design waveform means employing in establishing sa hardware optimized waveform for each phase encoding step the staπing and endi gradient values and moments for the slice, phase encoding and readout directions.

35. The apparatus of claim 18 including said apparatus being employable in breath-hold magnetic resonan cardiac imaging.

36. The apparatus of claim 33 including said apparams being employable with respect to at least one pulse sequence selected from the group consisting of velocity-encoded pulse sequences and flow-compensated pulse sequences.

37. The apparatus of claim 33 including said higher moment being a first moment.

38. The apparatus of claim 28 including said design waveform creating means having means for effecting said hardware optimized waveform as a trapezoidal waveform.

39. The apparatus of claim 28 including said design waveform means having means for transforming said hardware optimized waveform to scan plane coordinates.

40. The apparatus of claim 39 including said design waveform means employing the inverse of a rotation matrix in effecting said conversion to scan plane coordinates.