Patent Grant
11709212
The field of the currently claimed embodiments of this invention relates to magnetic particle imaging (MPI) devices and methods. MPI is an imaging modality that constructs images of magnetic nanoparticle tracers in a region of interest.
A pulsed magnetic particle imaging system according to an embodiment of the current invention includes a magnetic field generating system that includes at least one magnet, the magnetic field generating system providing a spatially structured magnetic field within an observation region of the magnetic particle imaging system such that the spatially structured magnetic field will have a field-free region (FFR) for an object under observation having a magnetic nanoparticle tracer distribution therein. The pulsed magnetic particle imaging system also includes a pulsed excitation system arranged proximate the observation region, the pulsed excitation system includes an electromagnet and a pulse sequence generator electrically connected to the electromagnet to provide an excitation waveform to the electromagnet, wherein the electromagnet when provided with the excitation waveform generates an excitation magnetic field within the observation region to induce an excitation signal therefrom by at least one of shifting a location or condition of the FFR. The pulsed magnetic particle imaging system further includes a detection system arranged proximate the observation region, the detection system being configured to detect the excitation signal to provide a detection signal. The excitation waveform includes a transient portion and a substantially constant portion.
A method of imaging an object using a magnetic nanoparticle tracer according to an embodiment of the current invention includes providing the object with the magnetic nanoparticle tracer; applying a spatially structured magnetic field that has an FFR such that the FFR and surrounding regions of the spatially structured magnetic field intercept the object under observation at a region containing at least a portion of the magnetic nanoparticle tracer; exciting a portion of the magnetic nanoparticle tracer by at least one of changing a property of the FFR or a position of the FFR; detecting changes in magnetization of the magnetic nanoparticle tracer resulting from the exciting while the property of the FFR and the position of the FFR are substantially constant to obtain a detection signal; repeating the exciting and detecting for a plurality of different locations of the FFR within the object to obtain a plurality of detection signals; and processing the plurality of detections signals to render an image of a region of the object.
A device for use with or as a part of a pulsed magnetic particle imaging system according to an embodiment of the current invention includes a pulsed excitation system arranged proximate a sample observation region, the pulsed excitation system includes an electromagnet and a pulse sequence generator electrically connected to the electromagnet to provide an excitation waveform to the electromagnet, wherein the electromagnet provides a magnetic field within the sample observation region to generate an excitation signal from a sample when held by the sample holder in the sample observation region. The device further includes a detection system arranged proximate the sample observation region, the detection system being configured to detect the excitation signal from the sample to provide a detection signal The excitation waveform includes a transient portion and a substantially constant portion.
Further objectives and advantages will become apparent from a consideration of the description, drawings, and examples.
Some embodiments of the current invention are discussed in detail below. In describing embodiments, specific terminology is employed for the sake of clarity. However, the invention is not intended to be limited to the specific terminology so selected. A person skilled in the relevant art will recognize that other equivalent components can be employed and other methods developed without departing from the broad concepts of the current invention. All references cited anywhere in this specification, including the Background and Detailed Description sections, are incorporated by reference as if each had been individually incorporated.
Some embodiments of the current invention are directed to new and/or improved MPI systems and methods. For example, some embodiments use an alternative paradigm for spatially homogeneous excitation or drive fields in which pulsed waveforms containing substantially flat or constant components are used during the MPI scanning process.
The term field-free region (FFR) is intended to refer to a portion of the structured magnetic field that is below the saturation field strength associated with a magnetic nanoparticle tracer when a distribution of the tracer is contained within an object of interest. Note that the FFR does not mean that the magnetic field is zero throughout the region. The FFR can be a localized region around a point, localized around a line or localized around a plane, for example.
The pMPI system 100 also includes a pulsed excitation system 108 arranged proximate the observation region 106. The pulsed excitation system 108 includes an electromagnet 110 and a pulse sequence generator 132 electrically connected to the electromagnet 110 to provide an excitation waveform to the electromagnet 110. However, one should note that the pulsed excitation system 108 is not limited to only one electromagnet 110. There could be one, two, three, four or more electromagnets without limitation as to the particular number of electromagnets. The electromagnet 110 generates an excitation magnetic field within the observation region 106 to induce an excitation signal from the object under observation by at least one of shifting a location or condition of the FFR. The phrase “condition of FFR” can be a size, specific shape, addition of a biasing field, and/or the magnetic field distribution within the FFR, for example. The pulsed excitation system 108 as well as various embodiments of waveforms generated therefrom will be described in more detail below with respect to some particular examples. However, the general concepts of the current invention are not limited to only those particular embodiments.
The pMPI system 100 further includes a detection system 112 arranged proximate the observation region 106. The detection system 112 is configured to detect the excitation signal from the object under observation to provide a detection signal 114. In some embodiments, the detection system 112 can include one or more receiver coils such as receiver coil 116. In some embodiments, there can be one receiver coil, or two receiver coils, or three receiver coils, etc. without limitation as to the particular number of receiver coils. The abbreviations Tx and Rx mean transmitter and receiver, respectively. Some embodiments of the detection system 112 will be described in more detail below. However, the general concepts of the current invention are not limited to only the particular examples described. The excitation waveform produced by the pulsed excitation system 108 includes a transient portion and a substantially constant portion. The waveform may include one pulse or a plurality of pulses. The waveform may also include one or more portions that are not pulses. The term “pulse” refers to a portion of the signal that includes at least one transient portion that is distinguishable from adjacent portions of the signal that are substantially constant. Substantially constant means a slow time variation relative to a time variation of the transient portion. For example, in some embodiments the transient portion can be changing with time twice as fast as changes in the substantially constant portion, three times as fast, ten times as fast, many orders of magnitude faster, or anything in between without limitation for the general concepts of the invention. For example, a pulse can be a simple form, such as, but not limited to, a spike. However, a pulse can have a more complex structure to include a plurality of transient portions and a plurality of substantially constant portions. A pulse can be short in duration, i.e., time, compared to the adjacent substantially constant portions of the signal. The adjacent substantially constant portions are not required to be equal. Short can be at least twice as short, ten times as short, many orders of magnitude shorter, etc. without limitation to the general concepts of the current invention.
In some embodiments, the pMPI system 100 further includes an electromagnetic shield 118 arranged to fully enclose the observation region 106 therein to electromagnetically isolate the observation region 106 from at least the magnetic field generating system 102 and an environment around the pulsed magnetic particle imaging system 100. The electromagnetic shield 118 can be a metal shield, such as, but not limited to, a shield formed from copper plate and/or an alloy of copper. The thickness of this shield may be prescribed by the expected or desired Tx and Rx bandwidths. For example, the thickness may be chosen such that the skin depth blocks all Rx bandwidth interfering signals, but passes certain signals such as slow shift magnetic fields without significant attenuation. For example, the shield can be substantially pure copper about 2-5 mm thick. However, the electromagnetic shield 118 is not limited to these examples and can be constructed from other metals in a variety of forms, thicknesses, etc. We note that the use of a fully enclosed shielding system is important for bi-directional EMI shielding. The sensitive Rx signal may be shielded from external or environmental interfering and noise sources as well as from interfering and noise sources originating in other parts of the system 100, including the magnet systems 104, 105, 130, 136, etc. and passive magnetic field focusing system 120 that produce and/or shift the FFR. Additionally, and importantly, these magnet systems 104, 105, 130, 136, and the passive magnetic field focusing system 120 are themselves shielded from the strong fields produced by the Tx coils 110, preventing magnetic saturation and other untoward interactions.
In some embodiments, as will be described in more detail below, at least a transmission portion (e.g., Tx coils 110) of the pulsed excitation system 108 and a reception portion (e.g., Rx coils 116) of the detection system 112 are enclosed within the electromagnetic shield 118.
In some embodiments, the magnetic field generating system 102 can also include a passive magnetic field focusing element 120 such as a soft magnetic material configured into a shape and arranged relative to the at least one magnet 104 so as to focus, i.e., concentrate, or otherwise guide the magnetic field lines to a desired path therein.
In some embodiments, the detection system 112 can be configured to detect the excitation signal substantially only during the substantially constant portion of the excitation waveform to avoid feedthrough interference during the transient portions of the excitation waveform to provide the detection signal. For example, electronic circuits can be provided to provide such functionality. The detection system 112 can be, or can include, hard-wired components and/or can include programmable structures. Either way, these are to be interpreted as structural components since a programmable structure can be replaced with hard-wired components such as, but not limited to, ASICs and/or FPGAs, for example. In some embodiments, the detection system 112 can be configured to detect the excitation signal during the transient portion of the excitation waveform, or during both the transient portion and the substantially constant portion of the excitation waveform.
In some embodiments, the substantially constant portion of the excitation waveform is at least 500 nanoseconds and less than 500 milliseconds. In some embodiments, the substantially constant portion of the excitation waveform is constant to within about 10% of a target amplitude of the excitation waveform. In some embodiments, the transient portion of the excitation waveform has a duration of at least 100 nanoseconds and less than 100 microseconds. In some embodiments, the excitation waveform includes a magnetization preparation portion and a readout portion such that the magnetization preparation portion includes at least a fraction of the transient portion and the readout portion includes at least a fraction of the constant portion. In some embodiments, the magnetization preparation portion dynamically configures a state of tracer magnetization that is in a vicinity of the FFR based on magnetic relaxation properties of the tracer prior to the readout portion. In some embodiments, the magnetization preparation selectively nulls a signal from tracer associated with a specified relaxation state or physical location during the readout portion. In some embodiments, the magnetization preparation encodes the tracer magnetization in such a way as to improve spatial resolution in reconstructed images. In some embodiments, the substantially constant portion is longer than a relaxation time for the magnetic nanoparticle tracer in the FFR. In some embodiments, the readout portion is longer than a relaxation time for the magnetic nanoparticle tracer in the FFR. In some embodiments, the substantially constant portion is long enough to establish a steady-state magnetization in the magnetic nanoparticle tracer in said FFR. In some embodiments, the readout portion is long enough to establish a steady-state magnetization in the magnetic nanoparticle tracer in said FFR.
In some embodiments, the excitation waveform includes a plurality of pulses, each pulse of the plurality of pulses including at least one transient portion. In some embodiments, the excitation waveform further includes a plurality of constant portions such that the excitation waveform includes a plurality of magnetization preparation portions and a plurality of readout portions such that each magnetization preparation portion of the plurality of magnetization preparation portions includes at least one transient portion of at least one of the plurality of pulses, and such that each readout portion of the plurality of readout portions includes at least a fraction of at least one constant portion of the plurality of constant portions. In some embodiments, the excitation waveform is a square wave. In some embodiments, the excitation waveform is trapezoidal.
In some embodiments, the FFR is a field-free “line” defining a longitudinal, or “long”, direction and the excitation waveform is at least partially applied in the longitudinal direction of the field-free line so as to change a condition of the field-free line by displacing a corresponding field-free line structure in a magnetic field space but maintaining shape and location of the field-free line structure. In some embodiments, the FFR is a field-free line defining a longitudinal direction and the excitation waveform is at least partially applied in the plane orthogonal to the longitudinal direction of the field-free line so as to change the position of the field-free line. Here, the phrase “field-free line” does not require a perfectly one-dimensional FFR with zero width in the mathematical sense. The “field-free line” has a width that is smaller than the length of the line. For example, the width can be smaller than the length of the line by at least a factor of 2, or a factor of 10, or a factor of 100, or a factor of 1000, without limitation to these particular examples.
In some embodiments, the pulse sequence generator is configured to provide the plurality of pulses such that each have at least one of preselected shapes, magnitudes, widths or inter-pulse periods to provide a particular pulse sequence encoding. In some embodiments, the pulse sequence generator is configured to provide an excitation waveform that includes a plurality of pulses separated by inter-pulse portions, and each pulse of the plurality of pulses includes a transient portion.
In some embodiments, the pulse sequence generator is configured to provide an excitation waveform that includes a plurality of pulses separated by constant inter-pulse portions, and each pulse of the plurality of pulses has a constant portion between transient portions. In some embodiments, the excitation waveform can approximate a finite-duration square wave with finite transitioning slew rates. In some embodiments, the excitation waveform is trapezoidal. In some embodiments, each constant portion between successive pulses of the plurality of pulses is greater than a preceding constant portion for at least a portion of the pulsed waveform. However, the general concepts of the current invention are not limited to these particular examples.
In some embodiments, the pulsed excitation system includes an LR circuit powered by a linear amplifier in a non-resonant filter chain. In some embodiments, the LR circuit has an inductance between 1 microHenry and 50 microHenries. In some embodiments, the LR circuit has an inductance between 1 microHenry and 30 microHenries.
In some embodiments, the pulsed excitation system includes a resonant switcher circuit. In some embodiments, the detection system includes gain control circuitry to amplify portions of the detection signal that have less feedthrough contamination relative to portions of the detection signal that have more feedthrough contamination.
In some embodiments, the pMPI system 100 further includes a signal processor 122 configured to communicate with the detection system 112 to receive detection signals 114 therefrom. The signal processor 122 is further configured to generate an imaging signal for rendering an image corresponding to regions of the object under observation traversed by the FFR. The signal processor 122 can be, for example, a computer in the general sense, which can be, but is not limited to a desktop computer, a laptop computer, a tablet device, a smart phone, a microcontroller or other system-on-a-chip device, or any networked combination of such devices or portions of such devices. The signal processor can be a programmable processor and/or a hard-wired processor, such as, but not limited to an ASIC and/or an FPGA. Regardless of whether the signal processor 122 is programmable or not, it is to be interpreted as a structural component of the pMPI system 100.
In some embodiments, the pMPI system 100 further includes an object holder 124 configured to be arranged within the observation region 106 of the pMPI system 100. In some embodiments, the pMPI system 100 further includes a mechanical assembly 126 operatively connected to at least one of the object holder 124 or the magnetic field generating system 102 to at least one of translate or rotate a relative position of the FFR.
In some embodiments, the pMPI system 100 further includes a slow-shift electromagnetic system 128 that includes a slow-shift electromagnet 130 and a slow-shift waveform generator 132 disposed proximate the observation region 106 of the pMPI system 100. The slow-shift waveform generator 132 provides waveforms to the slow-shift electromagnet 130 to shift a position of the FFR on a time scale that is slow compared to a timescale of the excitation pulses. In some embodiments, the slow-shift electromagnetic system 128 can include a plurality of slow-shift electromagnets, such as, but not limited to, slow-shift electromagnets 130, 134, 136, 138 for example. In an alternative embodiment, the pMPI system 100 further includes a slow-shift waveform generator 132 that is configured to communicate with the magnetic field generating system 102 to use the same coils as the magnetic field generating system 102, such as, but not limited to, coils 104 and/or 105. In some embodiments, a slow-shift linear motor is configured to move the object being imaged in at least one dimension on a timescale that is slow compared to a timescale of the excitation pulses. In some embodiments, the pMPI system 100 is further configured with a field-free line and a rotating gantry and motor to provide relative rotation between the field-free line and imaging bore 106 containing the sample under study 124. In some embodiments, one or more of the following are mounted to a rotating gantry: main magnets and soft magnetic flux guide (120), slow-shift magnets (130, 134, 136, 138), Tx and Rx coil(s) (110, 116), Tx/Rx circuitry (108, 112), and electromagnet shield (118).
In some embodiments, the magnetic field generating system 102 is dynamically configurable. In addition to dynamic behavior pre-programmed into the sequence generator 132, by the user beforehand, feedback loops may inform an online algorithm or controller logic to provide real-time modifications to excitation waveforms to more fully conform to a desired ideal trajectory. For example, field and/or current sniffer elements may report realized waveforms to the controller for appropriate negative feedback-controlled modification of inputs to the power amplifiers (see, e.g.,
In some embodiments, the pMPI system 100 further includes a signal processing and image rendering system 122 configured to be in communication with the detection system 112 to receive the detection signal 114. The signal processing and image rendering system 122 can be configured to process the detection signal 114 and render an image corresponding to portions of the object under observation containing the magnetic nanoparticle tracer and that was addressed by the FFR.
In some embodiments, the signal processing and image rendering system 122 is configured to process the detection signal 114 and render the image with a spatial resolution of at least 1.5 mm. In some embodiments, the signal processing and image rendering system 122 is configured to process the detection signal 114 and render the image with a spatial resolution of at least 1000 μm to 100 μm. In some embodiments, the signal processing and image rendering system 122 is configured to process the detection signal 114 and render the image to represent at least one of density, mass, concentration or a derivative thereof of the tracer at corresponding image locations. In some embodiments, the signal processing and image rendering system 122 is configured to process the detection signal 114 and render the image to represent a local relaxation time of the tracer at corresponding image locations. In some embodiments, the signal processing and image rendering system 122 is configured to process the detection signal and render the image to represent at least one of a local viscosity, a pH, a magnetic nanoparticle binding event, a local oxidation state, a concentration of a biochemical analyte of interest, or a functionalized magnetic nanoparticle interaction at corresponding image locations. In some embodiments, the signal processing and image rendering system 122 is configured to process the detection signal and render the image to represent kinetic information of local biochemical processes interacting directly or indirectly with a possibly functionalized magnetic nanoparticle tracer.
Another embodiment of the current invention is directed to a method of imaging an object using a magnetic nanoparticle tracer. The method includes providing the object with the magnetic nanoparticle tracer; applying a spatially structured magnetic field that has an FFR such that the FFR and surrounding regions of the spatially structured magnetic field intercept the object under observation at a region containing at least a portion of the magnetic nanoparticle tracer; exciting a portion of the magnetic nanoparticle tracer by at least one of changing a property of the FFR or a position of the FFR; detecting changes in magnetization of the magnetic nanoparticle tracer resulting from the exciting while the property of the FFR and the position of the FFR are substantially constant to obtain a detection signal; repeating the exciting and detecting for a plurality of different locations of the FFR within the object to obtain a plurality of detection signals; and processing the plurality of detections signals to render an image of a region of the object.
In some embodiments, the detecting occurs subsequent to the exciting of the portion of the magnetic nanoparticle tracer. In some embodiments, the detecting occurs during the exciting of the portion of the magnetic nanoparticle tracer. In some embodiments, the detecting occurs both during and subsequent to the exciting of the portion of the magnetic nanoparticle tracer.
In some embodiments, the method further includes administering the magnetic nanoparticle tracer to the object under observation, wherein the magnetic nanoparticle tracer administered comprises magnetic nanoparticles that have an ensemble average diameter of at least 10 nm and less than 100 nm. In some embodiments, the magnetic nanoparticles of the magnetic nanoparticle tracer are at least 25 nm and less than 50 nm. In some embodiments, the magnetic nanoparticles of the magnetic nanoparticle tracer are uniform to within a variance of 5 nm.
In some embodiments, the property of the FFR and the position of the FFR are substantially constant for at least 500 nanoseconds and less than 500 milliseconds during the detecting. In some embodiments, the property of the FFR and the position of the FFR are substantially constant to within about 10% of an amount of change of the FFR and the position of the FFR during the exiting. In some embodiments, the changing the property of the FFR or the position of the FFR has a duration of at least 100 nanoseconds and less than 10 microseconds.
In some embodiments, the processing the plurality of detections signals to render the image of the region of the object renders at least one of a magnetic particle density image, or a magnetic relaxation dynamic parameter image. In some embodiments, the processing the plurality of detections signals to render the image of the region of the object renders a local viscosity. In some embodiments, the processing the plurality of detections signals to render the image of the region of the object renders a local binding state of the tracer.
Another embodiment of the current invention is directed to a device for use with or as a part of a pulsed magnetic particle imaging system. This embodiment can be constructed from components that are similar to or the same as corresponding components in
The device can further include a signal processor 122 configure to communicate with the detection system 112 to receive and process the detection signal 114. The signal processor 122 can be further configured to process the detection signal to determine a magnetization relaxation time for magnetic particles within the sample. The device according to this embodiment can further include a sample holder defining the sample observation region.
The following describes some particular embodiments of the current invention in more detail; however, the general concepts of the current invention are not to be limited to these particular examples.
Pulsed Encoding Techniques
Pulsed waveforms can be used to encode the MPI signal in ways generally not attainable when using sinusoidal continuous waveforms. For example, it is possible to temporally separate transmit feedthrough interference from tracer signal. Pulsed waveforms can also allow for measuring the steady-state magnetic physics of the nanoparticle tracer distribution, potentially bypassing blur effects seen in continuous sinusoidal MPI due to finite magnetic response times associated with the tracer. At the same time, pulsed waveforms can provide methods and encoding regimes that allow quantification of these magnetic response times, giving rise to enhanced image contrast and the possibility for measuring additional information. For example, relaxation images that quantify a measure of tracer relaxation physics and report in a spatially resolved manner may also be produced, in addition to tracer density images. In other embodiments, fully four-dimensional imaging datasets may be provided. This data may be obtained simultaneously from the same excitation pulse sequences or via the sequential application of different excitation pulse sequences.
Pulsed waveform sequences may also lead to a second form of spatial information encoding, distinct from the Langevin saturation typically used in MPI, via magnetic relaxation dynamics. More generally, it is possible to use this new encoding to shape tracer magnetization prior to or during encoding periods to affect resolution, SNR, and image contrast in useful ways. For example, in current MPI approaches magnetic relaxation generally acts to degrade resolution and image quality, while in pulsed MPI techniques we can leverage magnetic relaxation to improve resolution.
Requirements for Pulsed Encoding
As shown in
Magnetic Saturation
MPI images magnetic tracers that can magnetically saturate. The relationship between the applied magnetic field and the nanoparticle magnetization can be described by a magnetization curve, known as a M-H curve, that gives a nanoparticle ensemble's magnetization response for an applied magnetic field (See
Establishment of an FFR with Gradient Fields
A typical MPI scanning system contains one or more active or passive magnetic field sources that are used to generate a particular spatial field pattern, referred to herein as a FFR. A FFR is defined as the spatial region where the applied magnetic field is below a magnetic saturation value. This can lie in contrast to nearby spatial regions where the applied field is above a magnetic saturation value. FFRs often, but not necessarily, contain linear magnetic field gradients so the magnetic field transitions smoothly between the FFR and nearby saturating magnetic field regions. The FFR may be a field-free point, field-free line, or of a more general shape as depicted in
Homogeneous Fast Excitation Fields
In a typical MPI system, the tracer signal is generated by applying a spatially homogeneous but time-varying magnetic field in superposition with the gradient fields that create the FFR. Homogeneous fields applied in this manner will generally act to shift or translate the FFR in space and in time while maintaining the shape of the FFR. In some embodiments, such as in the case of homogeneous excitation with an applied field coaxial with the line in a FFL, the addition of a homogeneous field does not shift the mean location of the FFR and instead adds a magnetic bias to transform the FFR into a low-field region (LFR), for example.
When using coils that produce homogeneous fields, homogeneity requirements may be on the order of <10% or better across the Field of View (FOV), defined as the amount of magnetic field deviation from nominal, although coils with poorer homogeneity may be used. In MPI, it can be instructive to distinguish between homogeneous excitation fields and homogeneous shift fields even though both types of fields spatially are mathematically equivalent. Excitation fields have high slew rates and provide for fast translation of the FFR that can stimulate the particles to induce voltage signals received by inductive receiver coils, while shift fields can be orders of magnitude stronger and orders of magnitude slower and are used to slowly translate the mean location of the FFR structure across a large FOV inaccessible to the excitation field alone. Pulsed MPI encoding is largely concerned with the faster excitation fields, although shift fields are also still used and enable encoding of larger FOVs, just as in the case of standard MPI.
Substantially Constant Pulsed Excitation Waveform Components
Whereas a typical MPI trajectory includes a sinusoidally-varying and continuously transmitting excitation field (e.g., sinusoid at 45 kHz), some embodiments of the present invention provide for excitation waveforms with non-sinusoidal scanning trajectories. These waveforms are referred to herein as pulsed waveforms. A key feature of these pulsed excitation waveforms is the inclusion of some periodic or non-periodic components in which a magnitude of the field is held substantially constant for a period of time. By substantially constant, it is meant that the magnitude of the field, after some initial transient or ramping period and/or before a subsequent transient or ramping period, does not deviate from some desired or ideal field magnitude by more than a defined error. Example magnitudes of error include but are not limited to 10% of the desired value, 5% of the desired value, and 1% of the desired value. The period of time over which the field is maintained at this substantially constant value may range from 500 nanoseconds to 500 milliseconds; however, this range is not limiting.
Time Varying Pulsed Excitation Waveform Components
Substantially constant values in the excitation waveform may be separated by pulsed waveform components of various shape primitives. In general, these primitives may be chained together with precise timing to create more complex pulsed waveforms. One important type of such a pulsed waveform primitive is one that rapidly transitions through rising and/or falling edges. During these rapid field transitions, the time varying field waveform may be of any shape. For example, the waveform may slew linearly, exponentially, sinusoidally, or superpositions thereof as illustrated in
Magnetic Tracer
The magnetic tracers used in MPI are characterized by both their steady state response, frequently depicted as a M-H curve (
The time-varying response of magnetic nanoparticles to an applied magnetic field is governed by magnetic relaxation. Magnetic relaxation occurs because magnetic nanoparticles are not able to instantaneously respond to an applied magnetic field. In traditional sinusoidal approaches to MPI, when a characteristic time associated with a particle reorienting its magnetic moment is at a rate on the order of the period of the excitation sinusoid, the magnetic relaxation can cause a significant image blur. We can approximate this magnetic relaxation as a low pass filter applied to the time domain signal, which in turn manifests as an image domain blur. The result is an increase in the width of the point spread function from what would be predicted by the steady state Langevin M-H curve.
Larger particles can theoretically have improved steady-state resolution, depending on how the signal is acquired. It is well known that the Langevin equation describing tracer magnetization improves cubically with magnetic core diameter. This means that tracers with larger core sizes, such as larger than 25 nm, have steeper magnetization curves than tracers with smaller core sizes, e.g., smaller than 25 nm. However, with the typical amplitudes and frequencies used in traditional sinusoidal excitation MPI, magnetic relaxation in larger tracers, e.g., diameter greater than approximately 25 nm, can cause significant image blur that can preclude the use of the larger particles for high resolution imaging. Some embodiments of the current invention provide for methods and devices that use pulsed excitation waveforms that allow sampling of the steady-state magnetization for larger particles with relatively long magnetic response times without loss of resolution.
Components of a Pulse Sequence
Waveform primitives, such as transient pulses and periods that are substantially constant may be combined to create more complex pulsed waveforms and MPI excitation pulse sequences. Two basic components of a MPI pulse sequence may include one or more distinct magnetization preparation primitives and one or more distinct signal readout periods as illustrated conceptually in
Substantially Constant Periods
One or more periods where the magnetic field is kept substantially constant may be used both during magnetization preparation and during signal readout periods. During periods where the magnetic field is kept substantially constant, the magnetization of the tracer distribution evolves toward a new magnetization steady-state. This steady-state can be spatially varying because, for example, of the presence of an FFR and the interaction between the spatial field pattern it creates and filed-dependent tracer relaxation properties, and will evolve according to the magnetic relaxation physics governing the tracer. In general, the evolution of the magnetization is approximately step response-like, and so the inductively received signal records an impulse-like response.
Temporal Relaxation Encoding by a Gradient Field
If a gradient field or FFR structure is present during a substantially constant period of a pulsed waveform, this leads to a second encoding of spatial information related to the tracer distribution due to spatially-variant relaxation times.
As
This second temporal-spatial encoding phenomenon can be exploited both in a magnetization preparation aspect of a pulse sequence and during image reconstruction.
Magnetization Preparation Periods
One or more pulse components may prepare, setup, or otherwise dynamically shape the magnetization state of a tracer distribution prior to, or possibly simultaneous with, a readout period.
Readout During Magnetization Evolution
A general feature of pulsed waveforms is the possibility to define distinct signal readout periods, for example, subsequent to magnetization preparation periods. In some embodiments, a readout period will be comprised of a substantially constant period in which the evolution of the magnetization of a tracer distribution moving toward a steady-state position is observed. In some embodiments, this readout may immediately follow a rapidly transitioning rising and/or falling edge. In some embodiments, the temporal duration of a substantially constant readout period will be purposely designed to ensure the achievement of magnetic steady-state of the tracer distribution in a sample of interest. In other embodiments, a substantially constant readout period will be purposely designed not to achieve steady-state conditions. In some embodiments, a readout period may also include non-constant components.
Pulse Types in MPI
Pulsed waveforms may be used to drive homogeneous excitation coils in MPI. These pulsed waveforms will affect the state and/or location of the FFR and generate signal in the presence of a magnetic tracer distribution.
Projective Excitation by Translation of an FFR.
In some embodiments, pulsed waveforms will be applied by homogenous coils such that the effect of pulses is to translate the location of the FFR structure. During substantially constant periods of these waveforms, the FFR structure is stationary in the FOV.
If two substantially constant components bracket an arbitrary dynamic trajectory, and the duration of both substantially constant periods is sufficient to establish steady-state magnetization conditions for the relative geometry between the FFR and tracer distribution associated with each substantially constant component (and there is no change in the tracer distribution in the intervening time period), then an initial magnetization state and final magnetization state are sampled without dependence on the particular intermediate trajectory. The instantaneous time-domain received signal, which measures the time derivative of the magnetization evolution, will depend on the specific trajectory chosen, but the integral of this signal may not. In this case, where the duration of the bracketing substantially constant components exceeds the time required for all magnetic relaxation phenomena to complete, then the signal may be viewed as encoding a full local projection in space along or near the line connecting the initial and final FFR locations. If steady-state is not established at one or another of the bracketing substantially constant periods, then this result will not apply: The difference in magnetization between the two components will be weighted by aspects of the specific trajectory.
Excitation in the Direction of an FFL
In some embodiments, excitation pulses may be applied in the direction along the line of a field-free line (FFL). As illustrated in
Biasing of an FFL
In some embodiments, a homogeneous bias field in the direction of an FFL may be applied generally with distinct excitation pulses applied separately, for example, in the orthogonal plane, along the line, or both. Such a bias field establishes an LFL in place of an FFL.
Pulsed Excitation Without an FFR
In some embodiments, homogeneous pulsed excitation waveforms will be applied to a sample without the presence of gradient fields or a spatially varying FFR structure. In this regime, a linear biasing field may be slowly ramped which will take the entire sample through a trajectory in the applied field magnitude space that is linearly proportional to the experience of a point source in a scanning system when linear gradient fields are used. Such a system can be used to rapidly test the effect of MPI pulse sequences and measure the one-dimensional (1D) PSFs associated with a tracer sample or sample object in aggregate. Note that such an approach excites all tracer, regardless of location, in the same manner. Thus the 1D PSF describes the aggregate or average behavior of all of the tracer in the observation volume. This type of scan can be accomplished with small tabletop relaxometer systems such as an arbitrary waveform relaxometer (AWR) or with a full imaging scanner if the magnetic sources that generate the FFR can be turned off.
Because sampling is only in one dimension—the applied magnetic field space—these spatially homogeneous scans are high throughput and, for example, may be completed in less than one second. A slew of scans testing different pulse sequence parameters or entirely different pulse sequences may easily be tested and vetted. Such scans can also quickly provide representative information about the entire tracer distribution in a sample, in terms of resolution, total tracer mass and SNR, and magnetic relaxation dynamics. This information may be used to guide choice of scan parameters such as the amplitude and period of substantially constant field values during a subsequent and generally longer imaging scan or series of scans. Moreover, this information may also be useful in choosing signal processing and reconstruction parameters such as the receive bandwidth and receive filter cutoffs. These values may be obtained through a standardized calibration procedure or on the fly before any scan using a fast gradient-less scan.
Pulse Sequences
The encoding machinery available in pulsed MPI, including various pulse types, conceptual pulse waveform components, and coil arrangements driven by the pulsed waveforms can be combined to create complete pulsed MPI pulse sequences to fully sample an imaging FOV. An exemplary pulsed MPI pulse sequence diagram is illustrated in
Some exemplary pulse sequences are described herein, including square wave-like pulse sequences, steady-state recovery sequences, and an illustrative more complex SNR efficient sequence. Neither the discussion of these specific pulse sequence embodiments nor the exemplary forms of these pulse sequences illustrated are to be taken as limiting.
Square Wave-Like Excitation Pulse Sequences
One pulsed MPI pulse sequence that illustrates many of the benefits of encoding the MPI signal with pulsed waveforms is the use of square wave-like excitation waveforms. These waveforms contain alternating periods of a substantially constant value of the magnetic field separated by rapid transitions through rising and falling edges. These components are depicted in
If the half-period of the square wave-like excitation waveform is greater than or equal to the magnetic relaxation or magnetic response time of the tracer distribution, then the tracer distribution will reach a steady-state in magnetization by the end of each square wave half-period. Here, the magnetic relaxation time or magnetic response time refers to the amount of time required for the magnetic response of the tracers to settle to a steady-state after each step-like transition of the square wave-like waveform. This time may be precisely defined in various ways. For example, if the magnetic impulse response of a tracer can be described adequately as one or more exponential time constants, then the response time may be defined as greater than or equal to some multiple of the largest or maximum time constant. For example: tchar*2, tchar*3, tchar*4, or tchar*5 where tchar is a mean, bounding, maximum, or characteristic exponential time constant. The experimental data of
Tracer magnetic response times are generally a function of the absolute magnetic field experienced by the tracer such that the zero field time constant is the longest, with a monotonic decrease in time constant with increasing applied field strength (and subject to limits at very high field strengths). The zero field time constant and/or low field time constants associated with a tracer, and also associated with the center of the FFR in an MPI system, will thus in general be associated with the longest relaxation times and will set the half-period requirements to ensure steady-state is achieved. This will also have ramifications on the receive signal bandwidth. The magnetic relaxation dynamics need not strictly follow a single or multiple exponential model. The ramifications of using a square wave-like (or more general) excitation waveform containing substantially constant periods pertain to any real magnetic physics as they are always characterized by finite temporal responses. In contrast, as long as the total applied field is continuously temporally varying, such as in the case of a sinusoidal excitation waveform or the non-substantially constant periods of a pulsed excitation waveform, the tracer distribution is subject to continuous excitation and the relaxation physics will be dynamically varying. In the special case of a substantially constant waveform, as in pulsed MPI, there is no continuous temporal stimulation. There is only a decaying response generated by magnetic relaxation of the tracer distribution evolving according to a fixed spatial pattern and the strong field dependence of magnetic relaxation physics. It is possible that other sources of time-variation in the factors that affect magnetic relaxation, e.g., physiologic sources such as viscosity, binding events, pH, or chemical reactions, may still lead to time-varying relaxation dynamics during a substantially constant field period, which may be separately visible, or encoded, by a pulsed MPI encoding scheme.
For a single point source, or a distributed sample excited without the presence of a gradient field/FFR, the magnetization at the end of each half-period will correspond to a single point on the steady-state M-H curve for the tracer at a given total field. For superparamagnetic tracers, as typically used in MPI, this corresponds to a point on the steady-state Langevin curve. For a distributed point source in the presence of a gradient field/FFR, the total magnetization at the end of each half-period corresponds to the spatial integral of the tracer distribution density mapped to corresponding points on the tracer's M-H curve according to the specific shape of the FFR in use. This relationship between steady-state inducing constantly held periods and sampling of the M-H curve is illustrated in
Steady-State Recovery Excitation Encoding
Steady-state recovery sequences can encode tracer response dynamics obscured in other encoding regimes or if multiple response dynamics that have greatly different response times contribute to the overall magnetization evolution and thus received signal.
In some embodiments, steady-state recovery pulse sequences are constructed with variable-duration substantially constant periods separated by short excitation pulses as illustrated in
The variation in signal associated with each pulse is directly related to the magnetic response time of the tracers when the FFR isocenter is located at the position associated with the interpulse period. In some embodiments, a small number of pulses will be associated with a given interpulse FFR isocenter location, such as in the range of 3 to 10. In other embodiments, 100 s of pulses may be applied per interpulse FFR isocenter location for a dense sampling of the relaxation dynamics. In some embodiments, a large range of interpulse times will be sampled, for example, in the range of 5 microseconds to 100 milliseconds to fully capture the entirety of all extant relaxation processes. Such broad ranges may necessitate sparse sampling of the interpulse periods or otherwise indicate a tradeoff; in some embodiments the interpulse period variation will follow a power pattern such as square, cubic, quartic and so on while in others exponential sampling may characterize the interpulse period sampling pattern. In some embodiments, much smaller ranges will be queried to increase the sampling density over a specific time range of interest.
In some embodiments, no gradient field or FFR structure will be present. The steady-state recovery sequence may then be applied to a sample with homogeneous excitation coils only to query the zero field relaxation dynamics of the sample in aggregate. In some embodiments, a bias waveform may be applied such that multiple steady-state recovery sequences are associated with mean total field values that are nonzero.
In other embodiments, the steady-state recovery sequence may be applied in the presence of an FFR, but along a low or zero-gradient field direction. For example, the steady-state recovery excitation may be applied along the line in an FFL scanning system. In this manner, it is possible to encode relaxation information as a function of total applied field but also spatially localized and associated with the mean location of an FFR or LFR isocenter.
Steady-state recovery approaches can have advantages over, for example, use of a square wave-like pulse sequence for encoding tracer relaxation dynamics. In general, the magnetic relaxation of a tracer may have both a significant Brownian and Neel component and/or significant relaxation components that occur very quickly, for example, within 10 microseconds, and other significant relaxation components that occur on the timescale of 10 s or 100 s of microseconds or even in the range of milliseconds. In a square wave-like encoding regime, the shortest relaxation components may completely dominate the time domain signal due to the time derivative nature of inductive signal reception. Encoding and observing longer and/or Brownian responses directly in the time-domain signal, which are of great interest for physiologic contrast and sensing, may be difficult in this regime because of time domain SNR, bandwidth restrictions, and/or noise floor limitations. A steady-state recovery approach, as described herein, can be a powerful tool for accurately and robustly measuring these longer and/or Brownian components of the magnetic relaxation signal. In general, multiple encoding regimes may be applied, such as including both square wave-like and steady-state recovery components in a single or sequential excitation waveforms. The time-domain data associated with each pulse in the steady-state recovery sequence may contain similar or equivalent information as a single square wave-like cycle. Using multiple encoding techniques simultaneously or in series can encode different relaxation information, for example in the case of tracers with heterogenous relaxation dynamics, in ways that lead to well-conditioned and robust reconstruction problems.
SNR Efficient Excitation Pulsed Waveforms
More complex encoding regimes and pulse sequences may be desirable, for example, to provide greater SNR efficiency.
Other System Sequencing
As illustrated in
Slow Shifting
To cover a large imaging field of view (FOV), e.g., larger than displacements associated with excitation waveforms, slowly varying homogeneous waveforms may be superposed with repeated excitation waveforms. This allows the mean location of the FFR isocenter to slowly translate over time as depicted in
Multiple slow shifting waveforms may be associated with each of the principal axes that are to be shifted in a given scan. For FFL scanning systems, two shifting components aligned with the plane orthogonal to the line of the FFL may be used in acquiring projection images and successive projection images may be taken after a rotation of the relative geometry between the FFL and sample about a fixed axis as depicted in
Gradient Waveforms
When using electromagnets as the source of the gradient fields that generate the FFR structure, gradient waveforms of different types may be desired depending on the application as illustrated in exemplary gradient waveforms shown in
It may also be desirable to dynamically modulate the strength of the gradient during an imaging scan as shown in the bottom of
Signal Processing and Reconstruction
Pulsed MPI encoding regimes admit multiple methods of forming images from received data.
After completion of digital signal processing steps, one or more resulting signals will be directly transformed, or combined and then transformed, into an image domain. A common pattern in pulsed MPI is to process the time domain signal such that a single value or a small number of values are associated with each of the repeated aspects of the excitation pulse sequence as illustrated in
It may also be desirable to use a system matrix method of reconstruction. In such a method, a pulsed MPI system matrix is constructed from many sets of pulsed MPI acquisition data. In particular, a point source may be placed at the center locations of final pixel/voxel locations. For each location, a pMPI acquisition, such as described in
One or more images constructed as described previously may undergo final postprocessing steps. For example, one or more images reconstructed at different resolutions may be combined into a single image, one or more images may be used to create one or more new images with colorized contrast, or one or more tracer density images may be combined with one or more relaxation images to create four-dimensional imaging data sets with a relaxation spectral information axis. In general, the production of distinct tracer density and tracer relaxation images is a feature of pulsed MPI.
Tracer Density Reconstruction
Images of the tracer density may be reconstructed using various gridding methods. The time domain data associated with substantially constant readout periods or more complex and repeated pulsed excitation waveforms may be integrated and the resulting value gridded to a pixel or voxel location. The gridding procedure may use the known x-space trajectory of the FFR isocenter to map time domain signals to image space locations. In some embodiments, the mean location of the FFR isocenter during a repeated pulsed excitation or the single location of the FFR isocenter associated with a substantially constant readout period will provide the rule for this mapping.
In some embodiments using an FFL-based or other projective FFR-based pulsed MPI system, projection reconstruction may first be applied to raw time domain data collected at different relative angles between the FFL/FFR and imaging volume prior to calculating image domain values and gridding or interpolating these values. In such a manner, time domain data may be associated with locations resolved in all three spatial dimensions prior to reducing the time domain data to image domain values.
Direct Integrated Gridding.
The integral of the time domain data associated with readout periods may be directly interpolated or gridded onto a regular pixel or voxel grid based on the location of the FFR during the readout. The grid location may also be partially determined by the locations of the FFR prior to readout, for example, during a magnetization preparation or a previous readout period.
Square Wave PSF Reconstructions with AWR
The experimental data points closest to the theoretical Langevin curve in
Integrated Gridding with Projection Reconstruction Techniques
In some embodiments, projection reconstruction techniques may be applied instead of direct integrated gridding. This may be a desirable approach to use larger pulsed excitation amplitudes without incurring a large projective blurring penalty. A pulsed excitation may be viewed as a local projection along the excitation direction. Quadrature excitation coils or other mechanisms of rotating the relative geometry between the sample and excitation direction of the pulses may be used to provide samplings of a local partial field-of-view (pFOV) at different excitation angles, enabling local projection reconstruction techniques as illustrated in
In other embodiments, many local pFOV projective excitations may be summed together or otherwise combined to obtain a full FOV projection. In this manner, a single full FOV projection reconstruction may be performed after many angles are sampled as illustrated in
These projection reconstruction approaches may provide improved SNR efficiency due to the SNR benefits inherent in projection reconstruction approaches. In general such a benefit will depend on the number of steps required per angle, the number of angles, and the number of projection reconstructions per slice. From an encoding and reconstruction perspective, when used in an excitation direction orthogonal to the line direction of an FFL, these approaches represent a way of emulating the use of a field-free plane or field-free slab FFR geometry. In general, inverse formulations, forward model formulations, and iterative reconstruction algorithms may be desirable to solve these projection reconstruction problems.
Magnetization Image Reconstruction
When a pulsed MPI excitation sequence is used in which the magnetization of the tracer distribution achieves a steady-state by the end of each successive readout period, it is possible to consider reconstructing images of the magnetization directly, either as an intermediate or final image.
In the case where steady-state is established during each readout or more generally between successively repeated pulsed excitation waveforms, the area under the curve of the signal associated with the received coil during each repeated excitation period is proportional to the finite difference between the absolute magnetization state established in one excitation cycle and the previous cycle. In some embodiments, this constitutes a recursive relationship between all successive excitations. With known initial conditions, these recursion equations can be solved to yield a signal proportional to the absolute magnetization associated with the end (steady-state) of each pulsed excitation. These values may be gridded to mean FFR locations, as described previously, to obtain a magnetization image.
In some embodiments, finite difference and algebraic methods such as inversion of a Toeplitz matrix may be employed to solve the recursion equations. Depending on the pulse sequence encoding, this method of reconstruction may suffer from poor conditioning and be ill-posed. A priori assumptions or specific encoding regimes may be leveraged to improve this situation. For example, non-Cartesian FFR trajectories in which a single FFR location is used as a reference or anchor between sampling of many successive FOV locations may be employed. In some embodiments a zero field excitation in which the excitation source and possibly shifting sources are switched off will establish an absolute reference and break the recursive relationship.
Such a magnetization image may be characterized by a PSF that lacks the finite support and even symmetry typically desired in medical imaging. In some embodiments, a well-posed magnetization image reconstruction will be followed by the application of finite difference or numerical derivative techniques to realize a final image with a Langevin PSF more typical in MPI. In some embodiments, both procedures will be formulated into one reconstruction algorithm, for example using a forward model description of the system. Smoothing or regularizing measures may be taken when performing such transformations to avoid significant amplification of image noise. These may leverage robust a priori assumptions, for example, assumptions about the spatial bandwidth of the expected tracer PSF. An image so constructed will not suffer from projective-blur as in the case of direct integrated gridding.
Relaxation-Weighted and Relaxation-Filtered Reconstruction
Some pulsed waveforms with magnetization preparation sequences will inherently weight the signal using the spatially variant magnetic relaxation conditions present in the context of a gradient or FFR as illustrated in
During reconstruction, selective integration, windowing, temporal decomposition, various temporal weightings, or filters can be applied to use this secondary spatial encoding to improve the resolution of tracer density images, possibly at the expense of SNR. This represents a coupling of the tracer density information with relaxation dynamics. One embodiment with resulting experimental data is shown in
In some embodiments, relaxation weighting may be performed by selective windowing or thresholding of the time domain data associated with readout periods and prior to signal integration and gridding steps. Windowing may be performed in various ways, and distinct images may be produced associated with one of many different windows applied to the data. For example, multiple non-overlapping windows may be applied to a data set as illustrated in
More general weighting procedures may be applied to the time domain data, where the windowing and thresholding discussed previously represents a binary or step-like weighting. Weighting procedures may be applied to individual or across multiple decomposed data sets prior to integration and gridding of the time domain data. For example, the temporal data associated with the entire time domain readout period, or with one or more temporal decompositions, may be weighted by a linear ramp or by more complex temporal weightings such as piecewise combinations of linear ramps and unity weights, polynomial functions, or exponential functions. These weightings may be used to accentuate the spatio-temporal encoding via the relaxation dynamics. Smoother or more continuous weightings may provide more flexibility in or more favorable conditions associated with the SNR-resolution tradeoff intrinsic in relaxation weighting as discussed herein. The ability to tradeoff SNR and resolution via temporal weighting is a second method of trading off SNR and resolution in pulsed MPI, independent, for example, of the SNR-resolution tradeoff possible in choosing excitation amplitudes when using direct integrated gridding reconstruction methods.
In some embodiments, filters may be applied to achieve relaxation-weighted improvements in resolution. In canonical, single-tone continuous wave systems, it is well established that higher harmonics of the fundamental excitation frequency are associated with higher spatial frequency information in the image domain, which leads to a fundamental resolution-receive bandwidth tradeoff in which increasing the overall received bandwidth is required to improve resolution and vice versa. However, due to the effect of relaxation and the use of quiescent half-periods long enough to establish steady-state magnetization, we can have a very different relationship in pulsed MPI where retaining only the data associated with a small subset of the lower frequency signal bandwidth can lead to improved resolution, as shown in
The signal from tracer farther from the FFR isocenter is characterized by signal that both occurs earlier in the readout period and which is also generally and relatedly characterized by temporally more rapidly varying signal components. From a Fourier domain perspective, this means higher frequency components. Thus, relaxation generally admits a spatio-frequency separability. Time-domain or frequency domain filters (after application of a Fourier transform) may also be used to improve resolution, similar to the application of temporal windows. Unlike many of the time-domain methods described previously, these frequencies may allow for improved resolution without loss of SNR. In some embodiments, it may be possible to improve resolution and SNR simultaneously. This is possible because, while filtering out data associated with higher harmonics does generally remove information and thus signal intensity, it also removes noise due to a reduced bandwidth. If the reduction in noise is greater than or equal to the loss in signal, we maintain or exceed the original SNR. In general, the exploitation of relaxation information in reconstruction may be provided by simple temporal windowing, filtering, or both.
Certain pulsed excitation waveforms may be viewed as simultaneous excitation of the tracer distribution with the fundamental frequency of the periodically repeating exciting waveform (e.g., square wave) along with higher harmonics per the Fourier decomposition of the particular excitation waveform. In this context, the received time domain response signal from the tracer will also be concentrated in small bands surrounding the harmonics of the fundamental frequency. The size of the bands around each harmonic is dictated by other features of the system such as the shift field slew rates. This leads to a similar receive bandwidth structure as in the canonical single-tone continuous wave systems. But as
In general, relaxation weighting, windowing, filtering, or more general decompositions of the time domain data may be used to generate many images of differing resolutions and SNR which may be later combined into one or fewer images using, for example, multi-resolution image processing techniques. The windows and decompositions used may be non-overlapping or may include overlapping. It is also possible and may be desirable to formulate system matrix approaches to relaxation-weighted imaging, similar to that described earlier for the general case. In the case of relaxation-weighted system matrix reconstruction, a transformation of the data prior to construction of the system matrix, such as using x-space gridded and windowed and/or filtered point source images as the matrix rows or columns will be required. In this manner, a single calibration dataset taken to construct a system matrix may be used to construct many different system matrices, including the cases with and without relaxation weighting/filtering. The raw data associated with a sample scan will need to be transformed accordingly prior to solving the inverse formulation.
Relaxation Image Reconstruction
In general, magnetic relaxation information is directly encoded in the received time domain signal in pulsed MPI. This magnetic relaxation information may be measured, fit, quantified, or otherwise characterized. This quantification can be done without the use of gradient fields in a non-imaging, sensor, or spectral format or may be performed in the context of gradient fields to produce relaxation images and/or 4D imaging data sets. When producing relaxation images, measures of magnetic relaxation phenomena are gridded to image space locations based on the known trajectory of the FFR isocenter to create images of relaxation with information distinct from and largely or completely orthogonal to tracer density images. With some encoding regimes, it is possible to associate a relaxation map, as a function of applied field, to each pixel or voxel in the image space. If no gradient fields or FFR structure is present, a relaxation map as a function of applied field that is representative of the entire sample volume may be rapidly produced. In other embodiments, a 4D dataset describing tracer density as a function of space and a relaxation-associated variable may be provided. In some embodiments, it is possible to use system matrix methods to construct relaxation images. Partial x-space gridding and fitting, measuring, or quantification of dynamic parameters may be desirable in the construction of a suitable system matrix.
In MPI, magnetic relaxation of the tracer may be characterized by multiple and possibly interacting phenomena, for example Neel, Brownian, or ferromagnetic physical relaxation processes. A distinction may also be made between faster magnetic relaxation processes and slower ones in heterogeneous scenarios, regardless of the specific nature of the relaxation physics. Pulsed MPI encoding admits multiple ways of encoding relaxation information, including with respect to these different components, such that reconstruction processes can observe, query, or accentuate different aspects of magnetic relaxation.
In general in pulsed MPI, and especially when using tracers with larger magnetic core sizes, magnetic relaxation phenomena may be sensitive to local microenvironmental processes and conditions. Variables such as viscosity, pH, reactive oxygen species concentration, biochemical state, tracer binding events, state or viability of the cell when a tracer is used to label cells, the kinetics or dynamics of various biochemical processes, and many others may all affect the magnetic relaxation properties of the tracer. This may be especially true when a tracer exhibits significant Brownian relaxation phenomena in which the tracer physically rotates, including both the magnetic core and outer shell, in the process of aligning with the local applied field. The temporal dynamics of a tracer undergoing Brownian relaxation are a strong function of both the applied field strength which acts as a magnetic torque and any sources of rotational drag. In pulsed MPI, when a substantially constant period of the excitation field is applied, the magnetization of the tracer distribution evolves toward a new steady-state based on the interplay between the magnetic torques imparted by the FFR structure and local microenvironmental conditions. When the same pulsed excitation is applied repeatedly at different mean locations of the FFR isocenter, then any differences in local microenvironmental variables of interest will provide contrast in terms of the observable magnetic relaxation. Many of these variables may be associated with physiologic processes of interest and may therefore directly or indirectly encode physiologic contrast into the MPI signal. When a relaxation image is reconstructed from these data, physiologic contrast not observable in the tracer density image may be provided.
Relaxation images so constructed may have various postprocessing steps applied to accentuate or identify imaging aspects of interest, for example, particular physiologic contrast or signals of interest. In some embodiments, thresholding of the relaxation image may be used to separate distinct phenomena of interest or to colorize the image. It is also possible to combine relaxation images and tracer density images, possibly constructed from data from the same scan, into one or more colorized images.
Direct Response Gridding
In some embodiments, the raw time domain signal associated with each readout period directly encodes a magnetization impulse response. For example, in the presence of a substantially constant period subsequent to a rapid transition, the magnetic tracer distribution will evolve to a new steady-state magnetization distribution in a step response-like manner. The time derivative applied by use of an inductive receive coil will realize an impulse-like response observable directly in the raw time-domain signal. In these circumstances, one or more characteristic relaxation parameters may be measured or fit for each readout period and gridded or interpolated onto a regular pixel or voxel grid to form a relaxation image. In some embodiments, the image domain grid point associated with each step is the known FFR isocenter location during that step.
Square Wave-Like Relaxation Map or PSF in the AWR
An AWR can be used to measure the aggregate relaxation time constant of a sample as a function of applied field strength, as shown in
Although in the embodiment shown in
In general, relaxation dynamics may be fit against phenomenological or physical models of magnetic tracer relaxation more complex than a simple linear first order process described by an exponential time constant. In some embodiments, relaxation parameters other than time constants may be fit, measured, or derived. For example, parameters such as peak delay time, FWHM, or finite support time as defined by some threshold such as an SNR threshold may be more instructive, better characterize the temporal dynamics, or accentuate contrast of interest in relaxation imaging.
Steady-State Recovery Image Reconstruction
In some embodiments, pulse sequences may query the time to achieve steady-state for tracers near a given mean location of the FFR using a steady-state recovery pulsed excitation waveform. Embodiments of this type of encoding and associated raw data are shown in
Many methods may be used to quantify the relaxation dynamics queried by a steady-state recovery excitation. In some embodiments, an envelope will be fit to the raw data as shown in
Excitation Along the Line
In some embodiments pulsed excitation will be directed along the line of an FFL. In this regime, the line is not translated in the orthogonal plane but instead transitions from an FFL to an LFL or from an LFL at one strength or polarization to an LFL at another strength or polarization. The mean location of the FFL or LFL in the orthogonal plane is constant regardless of the excitation pulses applied in the line direction. In some embodiments, it is possible to query relaxation of tracer along the line at different applied field strengths. Similar to relaxation maps constructed with pulsed encoding in an AWR or when no gradient is present, the relaxation dynamics of the tracer distribution may be probed against the total applied field strength with an excitation pulse sequence. However, this information is now spatially localized by the excitation line. This excitation may be applied at different mean locations of the FFR and/or FFL. Either in a projection mode or following full tomographic projection reconstruction, 1D relaxation maps or other non-scalar values may be associated with individual pixels and/or voxels. In general, these multi-value data sets may be reduced to scalar values of merit to yield scalar relaxation images as in other relaxation image embodiments. These data may also be reduced or transformed into 4D data sets describing tracer intensity as function of space and a temporal relaxation-associated variable.
Temporal Decomposition
As in the case of tracer density imaging, the spatio-temporal encoding inherent to the time-domain signal itself may be used in the construction of relaxation images. The same or similar weighting, windowing, or filtering procedures as applied in relaxation weighted tracer density imaging may be applied to the time domain signal prior to relaxation image formation. In some embodiments, a relaxation parameter of merit, such as an exponential time constant, may be fitted to each of a number of distinct windowed decompositions or distinct filtered versions of the time-domain signal. Because earlier aspects of the readout signal will be weighted more by fast relaxation processes, including those due to particles farther from the FFR isocenter and later aspects by slower relaxation processes, including those due to particles closer to the FFR isocenter, measures of relaxation parameters associated with different windows or filters provide a mechanism of filtering out or selecting for relaxation times associated with specific phenomena. In some embodiments, this procedure will improve the resolution of the relaxation images. In some embodiments, this procedure can be used to filter out or select for distinct relaxation states such as bound versus unbound, regardless of spatial weighting.
4D Imaging
Certain embodiments of forming distinct tracer density images and relaxation state images have been described. Mathematically, these can be described as yielding distinct images of the form rho(x, y, z) and r(x, y, z) where rho corresponds to a scalar tracer density image and r corresponds to a scalar relaxation image. In general, it may also be possible to combine the information more explicitly to form 4D image datasets of the form I(x, y, z, tau) in which a reconstructed dataset contains intensity information about tracer as a function of space and relaxation time. Such a reconstruction would yield a point-spread function (PSF) and images that more fully characterize the MPI process and provides spectral information about tracer dynamics. In some embodiments the scalar intensity variable I collapses to the tracer density rho if the 4D image is projected along the tau direction. In some embodiments, windowed decomposition techniques applied to the raw time domain data may be used to form such a 4D image dataset. In other embodiments, mathematical tools such as Laplace transforms may be leveraged to recover, deconstruct, or reconstruct relaxation time or relaxation spectral data. In other embodiments, the raw time domain data may be projected onto a suitable basis set, or transformed and then projected onto a suitable basis set. The basis set may, for example, be associated with known physical relaxation processes that characterize a tracer.
Postprocessing Approaches
Various postprocessing techniques may be leveraged in the context of pulsed MPI. In general one or more images may be reconstructed by the signal processor for each scan and one or more scans may be associated with a single sample. A set of images so obtained may be processed in various ways to produce a final output for a user. For example, single images may be processed, filtered, or thresholded. Multiple images may be similarly processed, filtered, or thresholded prior to being combined in some manner.
Multiresolution Image Combination
In some embodiments, a plurality of images may be obtained from encodings or reconstructions that lead to differing native resolutions and/or SNR levels. For example, the use of dynamically varying gradients, multiple scans taken at different fixed values of gradient magnitude, the use of dynamically varying pulsed excitation amplitudes, multiple scans taken at different fixed values of excitation amplitudes, and/or scans taken with different pulsed MPI encodings, e.g., through different magnetization preparation designs, may all be combined to produce a smaller set or single image. The final image or images may be designed to maximize one or more of resolution, SNR, and image contrast and may be dependent on the desired application.
Density and Relaxation Image Combination
In some embodiments, one or more distinct density and/or one or more distinct relaxation images may be combined into one or more new images. In some embodiments, relaxation information will be used to colorize a tracer density image to improve image contrast and map relaxation information onto the tracer density images. Relaxation information may be directly overlaid onto density images and/or anatomical reference images from another modality with a continuously varying color map. In some embodiments, continuously varying relaxation information may be thresholded or mapped onto distinct colors by other means. For example, forward model-based algorithms may be employed to unmix relaxation species of interest, such as MPI tracers in a bound or unbound state or in the separation of distinct tracers if multiple are applied to a sample volume. These results may be used to form standalone images or applied and overlaid onto density images as a color mapping.
Projection and Slicing of 4D Data
In some embodiments, signal processing and reconstruction will provide 4D MPI datasets where an MPI intensity value is mapped against three spatial coordinate variables and a temporal relaxation variable. Projections along various dimensions may allow seamless interrogation of the information-rich dataset. For example, projection along one or more spatial dimension will allow for visualization of relaxation spectral information against the remaining spatial dimensions. Projection along the relaxation dimension may recover a tracer density image. Slices through the 4D data set may provide easily visualized 3D data sets such as MPI intensities localized in space as one steps through the enumeration of the relaxation spectral variable.
Applications of Pulsed MPI Encoding
Enhanced Imaging through Use of Larger Core Tracers
Using pulsed encoding and reconstruction techniques, all current applications of MPI may be improved in terms of imaging parameters such as resolution, SNR, and image contrast. For example, use of tracers with mean core sizes larger than 25 nm with pulsed MPI techniques can improve resolution in current areas of active MPI development such as stem cell tracking, cancer diagnostic imaging, lung ventilation and perfusion imaging, cardiac and angiography applications, stroke diagnostic imaging, and perfusion imaging such as cerebral blood flow (CBF) and cerebral blood volume (CBV).
Examples of Applications of Pulsed MPI
The unique aspects of pulsed MPI such as quantification of relaxation information, the prospect of 4D imaging, and the ability for magnetization shaping through magnetization preparation may provide for new, previously unforeseen applications of MPI. With the ability to robustly quantify relaxation information, and especially that associated with environmentally-sensitive physics such as Brownian relaxation, these pulse sequences may be used across a broad range of new molecular imaging applications. Pulsed MPI pulse sequences can allow for all of the following, optionally with tracers with a mean magnetic core size of 25 nm or more: in vivo viscometric imaging, pH imaging, inflammation imaging, reactive oxygen species imaging, imaging of binding events using functionalized tracers, design and use of MPI magnetic switches, observation of biochemical kinetics via functionalized or specialized tracers, elucidating bleeding in the body, infection imaging, diagnosis of pulmonary embolism and perfusion disorders of the lung, sentinel node biopsy.
In some embodiments pulsed MPI methods and apparatuses are used to image and/or quantify relaxation properties for ex vivo samples of interest. These samples may have originated from whole organisms in which MPI tracers were introduced, e.g., by systemic delivery or local injection, or MPI tracers may have been introduced in a separate step after harvesting. These samples may include tissue and organ samples, biological fluids, and other biologically-derived samples. These samples may also include non-biological samples of interest such as tracer samples. In some embodiments, pulsed MPI methods and apparatuses are used to image and/or quantify relaxation properties for live animal subjects after MPI tracers of interest are introduced, such as by, but not limited to, systemic injection into the circulatory system or direct tissue injection. In some embodiments, pulsed MPI methods and apparatuses are used to image and/or quantify relaxation properties for human subjects after MPI tracers of interest are introduced, such as by, but not limited to, systemic injection into the circulatory system or direct tissue injection.
Pulsed MPI Hardware
Multiple hardware systems are required to realize pulsed MPI excitation in practice.
Generation of an FFR
A crux of all MPI systems is spatial localization of the signal through manifestation of an FFR structure of some type, as illustrated in
Moving, Exciting, and Manipulating the FFR-Sample Geometry
A typical MPI scan may be thought of as shifting the location of the FFR pattern relative to an imaging FOV of interest to sample the imaging FOV of interest as depicted in
In general, the same magnetic field sources, such as in the case of electromagnets, used to create the FFR pattern may be used to accomplish some or all components of the time-varying FFR trajectory. In other embodiments, distinct magnetic source components provide one or more of these shifting fields. Both of these cases are illustrated in
In general, fast (e.g., bandwidths >1 kHz) excitation or drive waveforms may be distinguished from slow shift waveforms (e.g., bandwidths <10 Hz) for conceptual clarity, differences in the way they interact with system components such as eddy-current shielding subsystems, and in consideration of the imaging signals. However, both types of waveforms act to shift or condition the FFR pattern over time, possibly differing substantially in field magnitude and power requirements in addition to speed or slew rates. When using electromagnets to shift linear gradient FFR designs, excitation and shift waveforms are typically applied to coils arranged to produce spatially homogenous, but time-varying, fields. Homogeneous fields will act to shift or translate many FFR implementations in time while maintaining the spatial integrity of the FFR pattern throughout. In some embodiments, such as excitation along the line of an FFL, the FFR is not translated but rather transformed as previously described. When using such homogeneous coils, homogeneity requirements may be on the order of 90% or better, in terms of deviation in coil sensitivity over the imaging FOV, although coils with poorer homogeneity may be used.
Pulsed Excitation Transmit Systems
In pulsed MPI, fast excitation or transmit hardware systems must be capable of realizing time-varying magnetic fields that are described by pulsed waveforms. Whereas a typical MPI trajectory includes a fast, temporally sinusoidally-varying excitation or drive field, some embodiments of the present invention contemplate the use of pulsed excitation waveforms with non-sinusoidal temporal variations that lead to non-sinusoidal scanning trajectories. These non-sinusoidal excitation waveforms may include periodic or non-periodic components in which a specific value of the field is held substantially constant for a period of time as well as rapid transitions through rising and/or falling edges, as well as more arbitrary transient waveform components.
To realize these pulsed waveforms, specific transmit electronics and coil designs are required. The ultimate goal is for the current through one or more excitation coils, and therefore the induced magnetic field, to follow a prescribed temporal trajectory defined by pulse sequences. In particular, the ability to support both substantially constant periods of various durations and rapid transients and/or excitation waveforms with larger bandwidths than the single tone of traditional sinusoidal transmit systems is required.
Linear Amplifier and Low Inductance LR Circuit
In some embodiments, the transmit, drive, or excitation electronics chain contains an LR circuit powered by a linear amplifier in a non-resonant filter chain that does not substantially alter this inductance as illustrated in
A pre-emphasized voltage waveform so constructed may be supplied to the power amplifier in real time during a scan to achieve maximal field slew rates with minimal ringing or other forms of distortion. Common mode rejection, filtering components, and shielding systems may be installed. For maximum flexibility in temporal resolution in applying pre-emphasis, care must be taken in choosing the sampling rate for the waveform controlling the power amplifier. This sampling rate may be approximately 1 MHz or may require higher sampling rates such as 10 MHz or higher. Similar results in terms of driving coils with pulsed waveforms may be obtained by operating a linear amplifier in controlled current mode in which an internal feedback control loop works to achieve pulsed waveform components such as rapid transitions through rising and/or falling edges given an ideal waveform target. To achieve prompt and rapid feedback adjustment, the linear amplifier may be modified to allow for high voltage slew rates, or in the case where high voltage slew rates are not possible, set up with a predictive computational algorithm or electronic circuit to avoid instability in the feedback loop and avoid ringing.
Resonant Switcher Circuit
In some embodiments, a switcher circuit may be used to realize pulsed excitation waveforms. In this approach, a transmit circuit system switches between various states according to precise timing diagrams as depicted in
Transmit Coil Designs
The pulsed waveforms may be applied to transmit, drive, or excitation coils. These excitation coils will generally be designed to provide a spatially homogeneous but time-varying field to the imaging FOV. To facilitate certain desirable signal encoding strategies and MPI trajectories in pulsed MPI, this disclosure contemplates multiple transmit coil arrangements. To support certain rotating or projective encoding strategies, it may be desirable to have excitation coils with principal axes that are orthogonal to each other. For example, in an FFL MPI system in which the direction of the line is taken as the y-axis, excitation coils may be placed concentric with the x- and z-axes as shown in
To achieve the rapid transitions through rising and/or falling edges required in some pulse sequences, high slew rates (of voltages, currents, or magnetic fields) are required and it is important to prevent or minimize sources of ringing and distortion. Excitation coils may be designed with a low inductance, for example, with an inductance value between 1 microHenry and 30 microHenries. This will facilitate faster current transitions (shorter rise times for each pulse) in embodiments where the coil is part of an LR circuit and driven with a linear voltage amplifier. In some embodiments, the transmit coil is a single wire loop that may be hollow for passage of cooling fluid, such as a liquid or gas. In some embodiments Litz wire will be used to construct excitation coils. These coils may be of solenoidal, saddle, Helmholtz-like, or other coil designs, depending on factors such as geometric relationship with the fixed imaging bore and the need for low inductance coils. In some embodiments, the transmit coil is made up of a single Helmholtz or Maxwell pair. In some embodiments, the transmit coil has a small diameter to reduce inductance.
In some applications in which an FFL-like scanning system is used, it may be desirable to install transmit coils concentric with the axis of the line of the FFL as illustrated in
Methods of Shifting the FFR-Sample Relative Geometry for Full FOV Coverage
As illustrated in
An important component of an imaging scanner is shielding between an inner region and an outer region. As depicted in
Cooling
Many aspects of an MPI scanning system may require cooling during operation. In some embodiments, the transmit or excitation coils will be actively cooled as well as shift magnets and FFR-generating electromagnets. In some embodiments, aspects of the eddy-current shielding system may be cooled as well. In some embodiments, water or alternative cooling fluids such as the fluoro-carbon based Fluorinert will be used to cool hollow electromagnet wire and coils. In some embodiments, elements will be cooled by placement in thermal fluid baths.
Safety Considerations.
Care must be taken in using such large slew rates in terms of specific absorption rate (SAR) considerations for the subject. In some cases, magnetostimulation effects may apply as well. The high slew periods in pulsed MPI are by definition short and thus may represent a small duty cycle inside of a single excitation period. In the context of the larger periodicity of the pulsed excitation waveform, pulsed MPI waveforms will be subject to the same or similar measures of SAR limitation as, for example, calculated in a root-mean-square (RMS) sense.
Pulsed MPI Receive Systems
In MPI, inductive receiver coils are used to detect the magnetization signal from magnetic nanoparticle tracers. Traditionally, a continuous sinusoidal excitation waveform is applied while the signal is received. Because the transmit and receive coils are strongly coupled, in general, this leads to significant excitation feedthrough into the receiver coil. Without mitigation, this feedthrough is orders of magnitude larger than the tracer signal. Geometric decoupling strategies such as use of receiver coils wound in a gradiometric fashion with zero net area with respect to the geometry defined by the transmit coil are typically used.
Pulsed MPI has direct feedthrough across a wide frequency band. In some embodiments inductive decoupling via finely-tuned gradiometers and/or active cancellation will be employed. In some embodiments, inductive decoupling is improved by spatial-shimming of the gradiometer coils relative to the transmit coil. In some embodiments, inductive decoupling is improved by a tunable voltage divider (to finely adjust amplitude of a portion of the receive coil) and/or a tunable capacitor (to finely adjust the phase differences between portions of the receive coil).
In addition to a gradiometric receiver coil design, a typical sinusoidal MPI receive chain may contain a notch filter at the fundamental excitation frequency as well as one or more other filters (low pass, high pass, or band pass) that, among other purposes, endeavor to fully remove all received signal at the fundamental frequency. It is typical to remove this component due to an inability to separate tracer signal at the fundamental frequency from excitation feedthrough and possibly the linear diamagnetic signal of imaged materials in the case of large samples. Loss of the tracer signal at the fundamental harmonic has major ramifications for signal linearity and shift invariance (LSI) in MPI and necessitates important steps in reconstruction to recover LSI. This requirement adds constraints on signal acquisition and encoding trajectories, adds complexity to the reconstruction process, and reduces overall signal-to-noise ratio (SNR). Further digital signal processing and baseline subtraction approaches may be leveraged for feedthrough mitigation and signal conditioning once the received signal is digitized and prior to image reconstruction.
A key aspect of pulsed MPI is the ability for temporal decoupling of the feedthrough from the nanoparticle signal. This is not possible in standard continuous sinusoidal wave MPI. When using pulsed waveforms, the circumstances of excitation feedthrough will be very different, in general. When using pulsed waveforms with substantially constant components over a period of time, the excitation feedthrough is not present after the substantially constant value is attained and an amount of time has passed in which the system response falls to zero. Furthermore, if the substantially constant value is arrived at after a rapidly transitioning rising and/or falling edge, and the system response time is very short as previously described, then there is only a very short period of time during each excitation cycle or period in which feedthrough is significantly present. This allows complete temporal separation of excitation feedthrough from tracer signal in the limit of short system response times. For finite system response times, the degree of feedthrough and signal separation is a function of the system response time and the magnetic tracer response time. In general, system response times much shorter than the tracer response time may be desired.
Although direct feedthrough can be limited to a very short window of time with pulsed encoding strategies, without mitigation, it may damage sensitive pre-amplifier systems set with gains ideal for the part of the signal in which the magnetic tracer signal alone is present. If the gain of the pre-amplifier system or systems is reduced to accommodate this peak, one or more of reduced dynamic range, reduced SNR, and increased discretization error will result. To avoid these adverse effects and maximize received signal dynamic range, multiple strategies may be leveraged. In some embodiments, a variable gain amplifier selectively amplifies the temporally delayed nanoparticle signal while not amplifying the initial direct feedthrough signal. In some embodiments, a blanking signal, precisely correlated with the excitation feedthrough aspects of each excitation cycle, period, or pulse, may be sent to a variable gain amplifier (VGA) system to switch between a low gain during the period of feedthrough and a high gain at other times. The system response time may be characterized and the blanking signal used to lower the gain of the receiver system beginning some finite time before the system response begins and ending some small finite time after the system response decays to a small value of choice. Subject to the response time of this VGA system itself, the excitation feedthrough signal can be greatly attenuated and the dynamic range and total amplification of desirable signal in readout can be maximized.
In some embodiments, crossed diodes circuits which reject voltages beyond a certain threshold value, as depicted in
In many pulsed MPI encoding strategies, especially those that leverage rapidly transitioning rising and/or falling edges, a significant part of the tracer signal may be contained in a short period of time following the rising and/or falling edge. It is important that the receiver system is configured with a bandwidth capable of fully capturing this signal in a Nyquist-limit sense. Sampling rates in the range of 1 MHz may suffice, but sampling rates up to 10 MHz or more may be required.
In some embodiments, detection of the raw received signal across a wide bandwidth is desired. In such embodiments, a receive bandwidth ranging from DC up to 10 MHz or more may be desired. In other embodiments, a reduced signal bandwidth may be desirable to modulate SNR and/or spatial resolution in image reconstruction. In some embodiments, a large signal bandwidth up to 10 MHz or more will be sampled by the data acquisition system and subsequent digital signal processing will reduce the bandwidth as desired during reconstruction. In other embodiments, the receiver systems may be configured with bandwidth-limiting filters for noise suppression. For example, a small bandwidth may be kept around the fundamental harmonic (this bandwidth determined by system parameters such as shift field slew rates and a fundamental frequency referring to the periodic excitation waveform), or a bandwidth around some small number N of harmonics of the fundamental periodic excitation waveform may be desired.
Modularity of Pulsed Hardware Components
Pulsed MPI requires major deviations from traditional MPI methods in the excitation and receive chains. Other aspects of the hardware may be maintained. For example, the magnetic field source components that generate the FFR and slowly (relative to the excitation periods) shift the FFR to cover a large imaging FOV may be constructed and used in the same or a similar manner as used in a sinusoidal MPI system. It may be desirable to design an MPI scanning system in which modular pulsed MPI excitation and receive systems are swappable or switchable with modular sinusoidal excitation and receive systems, as illustrated in
Pulsed MPI Signal and Imaging Equations in 1D
Here we provide a basic derivation of one dimensional pulsed MPI methods using ideal square wave excitation. This analysis is for a small number of idealized pulsed MPI methods, is intended to further illustrate some of the basic concepts, and is in no way limiting.
Steady-State Square Wave Imaging
The steady-state magnetization of a distribution of superparamagnetic magnetic particles in a field of view is described as:
MSS(t,x)=mρ(x)(kH(t,x))
where m is the magnetic moment of a single particle, rho is the particle density, L is the Langevin function, k is a tracer-specific term, and H is the applied field.
The 1D square wave excitation contribution to the location of the MPI FFR can then be expressed as:
where Delta x is the x-space displacement of the FFR and is equal to the square wave excitation amplitude in Tesla multiplied by the linear gradient strength of the MPI device in Tesla per meter. The voltage on a receiver coil, by inductive physics is:
Considering the response at each n-th step excitation, we can derive a signal equation:
where Delta Mn(x) is the difference between the steady-state magnetization states before and after the n-th square wave step excitation, hm(t, x) is the magnetic particle relaxation impulse response that is a function of position because of the strong field-dependence of tracer relaxation and the presence of strong gradient fields in MPI. This signal equation expresses that, in the ideal square wave case, the received voltage signal is a scaled form of the tracer relaxation impulse response integrated over space. This explicitly shows how tracer impulse responses, as well as steady-state information (via Delta Mn) are directly encoded in the raw time domain data.
We can further establish an indexed signal equation by integrating the signal under each half-period as previously described:
where we assume the half-period of the square wave excitation is long enough to fully establish steady-state conditions everywhere by the end of the half-period. In this case, the impulse response components integrate to unity and we are left with a finite difference between steady-state magnetization states before and after the step. We can express this result as a convolution relation as shown in the final step. 2Δx=2 AG−1 is the FFR excursion distance in a square wave step excitation where A is the square wave amplitude and G is the MPI gradient strength.
Finally we can establish an imaging equation by dividing out the known scalar quantities:
where our imaging PSF h(x) itself can be decomposed into a steady-state Langevin component, with no relaxation effects, and a rectangular blur which is controllable through our choice of excitation amplitude.
Relaxation Imaging
It is possible to directly quantify the magnetic relaxation process that evolves over each square wave half-period and construct relaxation images with this information. In one such embodiment, we can approximate our received signal for each half-period as a weighted exponential:
sn(t)≈B1ΔMnAn exp(−t/πn)|t
We can then fit the time constant τn and the amplitude An for each half-period and grid to the mean location of the FFR for the half-period. We may also enforce specific a priori knowledge or expectation, such as An=τn−1, if we believe a Debye exponential model adequately describes the relaxation process. This provides an image of relaxation dynamics over space rather than tracer mass/concentration.
Time Domain Relaxation Weighted Imaging
We can understand the effect of relaxation weighting to improve spatial resolution by exploiting the secondary spatial encoding of relaxation dynamics as follows. Unlike our steady-state square wave analysis, we do not integrate the entire signal associated with a square wave half-period, but only some latter part specified by a windowing threshold.
where srw(n) is our indexed signal equation for relaxation weighted reconstruction. Delayed windowing in this manner leads to the presence of a spatial weighting function wn(x). Following similar analysis as in the fully integrated case, we arrive at a modified imaging equation for relaxation weighted reconstruction:
where our relaxation weighted PSF, hw, is the fully integrated, steady-state square wave PSF (h(x)) point-wise multiplied by our weighting function w(x). In the case of an exponential relaxation dynamics, the relaxation function may be approximated as follows
Our relaxation weighted PSF can be understood as the relaxation-free steady-state square wave PSF point-wise multiplied by another narrow function, allowing us to improve resolution, potentially beyond the Langevin limit.
More complicated pulsed MPI excitation trajectories can be mathematically described in a similar manner. Arbitrary pulsed waveforms may most naturally be expressed as recurrence relations, for example. These few 1D mathematical descriptions of specific square wave approaches are intended for further illustration of these simple examples only and are in no way limiting.
The embodiments illustrated and discussed in this specification are intended only to teach those skilled in the art how to make and use the invention. In describing embodiments of the invention, specific terminology is employed for the sake of clarity. However, the invention is not intended to be limited to the specific terminology so selected. The above-described embodiments of the invention may be modified or varied, without departing from the invention, as appreciated by those skilled in the art in light of the above teachings. It is therefore to be understood that, within the scope of the claims and their equivalents, the invention may be practiced otherwise than as specifically described.
The application is a continuation of U.S. application Ser. No. 15/998,525, filed Aug. 16, 2018, which claims priority to U.S. Provisional Application No. 62/546,395 filed Aug. 16, 2017, the entire contents of which are hereby incorporated by reference.
This invention was made with government support under Grant Numbers EB019458 and MH106053 and awarded by the National Institutes of Health. The government has certain rights in the invention.
| Number | Name | Date | Kind |
|---|---|---|---|
| 5965214 | Crossfield et al. | Oct 1999 | A |
| 6054924 | Dames et al. | Apr 2000 | A |
| 6076007 | England et al. | Jun 2000 | A |
| 6144300 | Dames | Nov 2000 | A |
| 6204766 | Crossfield et al. | Mar 2001 | B1 |
| 6230972 | Dames et al. | May 2001 | B1 |
| 6323769 | Dames | Nov 2001 | B1 |
| 6369965 | Dames et al. | Apr 2002 | B1 |
| 6486655 | Crossfield | Nov 2002 | B1 |
| 6577237 | Dames | Jun 2003 | B1 |
| 6595419 | Doyle et al. | Jul 2003 | B1 |
| 20060238194 | Gleich et al. | Oct 2006 | A1 |
| 20100066363 | Brazdeikis et al. | Mar 2010 | A1 |
| 20120119739 | Gleich | May 2012 | A1 |
| 20120126800 | Vernickel et al. | May 2012 | A1 |
| 20130241548 | Gleich et al. | Sep 2013 | A1 |
| 20150015247 | Goodwill et al. | Jan 2015 | A1 |
| 20170020407 | Weber | Jan 2017 | A1 |
| Number | Date | Country |
|---|---|---|
| 3143929 | Mar 2017 | EP |
| 2006525506 | Nov 2006 | JP |
| 2012525193 | Oct 2012 | JP |
| 101 623 116 | May 2016 | KR |
| 9631790 | Oct 1996 | WO |
| 9748990 | Dec 1997 | WO |
| 9813708 | Apr 1998 | WO |
| 9815851 | Apr 1998 | WO |
| 9909436 | Feb 1999 | WO |
| 9948044 | Sep 1999 | WO |
| 0010123 | Feb 2000 | WO |
| 2011116229 | Sep 2011 | WO |
| Entry |
|---|
| International Search Report dated Dec. 14, 2018 from corresponding International Application No. PCT/US2018/000231. |
| Written Opinion dated Dec. 14, 2018 from corresponding International Application No. PCT/US2018/000231. |
| Karsten et al., “The Use of Flying Null Technology in the Tracking of Labware in Laboratory Automation”, 4 pages. |
| Crossfield, “Have null, will fly”, IEE Review, (Jan. 2001), pp. 31-34. |
| Sparavigna, Manuscript entitled “Labels discover physics: the development of new labelling methods as of promising research field for applied physics”, 16 pages. |
| Extended European Search Report issued in corresponding European Application No. 18846668.4 dated Apr. 9, 2021. |
| English translation of as-filed Korean Application No. 20150090507, which was published as KR 101 623 116. |
| Notice of Reasons For Rejection issued in corresponding Japanese Application No. 2020-508019 dated Mar. 28, 2023. |
| Number | Date | Country | |
|---|---|---|---|
| 20220221537 A1 | Jul 2022 | US |
| Number | Date | Country | |
|---|---|---|---|
| 62546395 | Aug 2017 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 15998525 | Aug 2018 | US |
| Child | 17581684 | US |
