METHOD AND APPARATUS FOR DETERMINING RESPIRATION PHASE, MAGNETIC RESONANCE IMAGING METHOD AND SYSTEM
20210052187 ยท 2021-02-25
Assignee
Inventors
- Jin Qiang He (Shenzhen, CN)
- De He Weng (Shenzhen, CN)
- Shu Qun Xie (Shenzhen, CN)
- Fang Dong (Shenzhen, CN)
Cpc classification
A61B5/7285
HUMAN NECESSITIES
A61B5/055
HUMAN NECESSITIES
A61B5/7264
HUMAN NECESSITIES
A61B5/7275
HUMAN NECESSITIES
A61B5/7292
HUMAN NECESSITIES
A61B5/7271
HUMAN NECESSITIES
A61B5/0816
HUMAN NECESSITIES
International classification
A61B5/055
HUMAN NECESSITIES
A61B5/00
HUMAN NECESSITIES
Abstract
The present disclosure provides techniques for determining a respiration phase by extracting a distance characteristic value, a score characteristic value, and an area characteristic value from the respiration signal, wherein the distance characteristic value, the score characteristic value and the area characteristic value are used to indicate waveform variation between two adjacent waveforms in the respiration signal. The techniques include training a respiration signal model according to the distance characteristic value, the score characteristic value, and the area characteristic value to determine the respiration phase of the respiration signal using the respiration signal model.
Claims
1. A method for determining a respiration phase of a respiration signal, comprising: extracting, via one or more processors, a distance characteristic value, a score characteristic value, and an area characteristic value from the respiration signal, the distance characteristic value, the score characteristic value, and the area characteristic value indicating a waveform variation between two adjacent waveforms in the respiration signal; training, via one or more processors, a respiration signal model according to the distance characteristic value, the score characteristic value, and the area characteristic value; determining, via one or more processors, the respiration phase of the respiration signal using the respiration signal model; and performing, via a magnetic resonance imaging (MRI) apparatus, magnetic resonance imaging of an examination region of an examination subject using the determined respiration phase.
2. The method as claimed in claim 1, wherein the distance characteristic value represents a ratio of (i) a width of a second-half-waveform of a first waveform of the two adjacent waveforms in the respiration signal to (ii) a width of a first-half-waveform of a second waveform of the two adjacent waveforms in the respiration signal.
3. The method as claimed in claim 1, wherein the score characteristic value represents a ratio of (i) a width of a rising waveform portion leading to a crest point of a first waveform of the two adjacent waveforms in the respiration signal to (ii) a width of a falling waveform portion leading to a trough point of the first waveform.
4. The method as claimed in claim 1, wherein the area characteristic value represents a ratio of (i) an area formed by two crest points of the two adjacent waveforms in the respiration signal and a trough point between the two adjacent waveforms in the respiration signal to (ii) an area formed by a second waveform of the two adjacent waveforms in the respiration signal.
5. The method as claimed in claim 1, further comprising: prior to extracting the distance characteristic value, the score characteristic value, and the area characteristic value from the respiration signal: acquiring a crest point of a first waveform and a crest point of a second waveform of the two adjacent waveforms in the respiration signal as a first maximum point and a second maximum point, respectively; acquiring a trough point between the two adjacent waveforms in the respiration signal as a first minimum point, and acquiring a trough point of the second waveform after the second maximum point as a second minimum point; and determining time points corresponding to the first maximum point, the first minimum point, the second maximum point, and the second minimum point, respectively, as a first maximum time point, a first minimum time point, a second maximum time point, and a second minimum time point.
6. The method as claimed in claim 5, wherein extracting the distance characteristic value from the respiration signal comprises: calculating a ratio of (i) a difference value between the first minimum time point and the first maximum time point to (ii) the difference value between the second maximum time point and the first minimum time point; and extracting the distance characteristic value as the calculated ratio.
7. The method as claimed in claim 5, wherein extracting the score characteristic value from the respiration signal comprises: determining, based on a predetermined amplitude ratio parameter, the first maximum point, and the first minimum point, (i) a start time point of a rising waveform portion leading to the first maximum point of the first waveform, and (ii) a start time point of a falling waveform portion leading to the first minimum point of the first waveform, as (i) a first time point, and (ii) a second time point, respectively, for determining an end-of-expiration start time point of the respiration phase; calculating a ratio of (i) a difference value between the first maximum time point and the first time point to (ii) a difference value between the first minimum time point and the second time point; and extracting the score characteristic value as the calculated ratio.
8. The method as claimed in claim 7, wherein the first time point is determined by evaluating:
mag(t_1)=mag(t_max1)*p+mag(t_min1)*(1p), and wherein the second time point is determined by evaluating:
mag(t_2)=mag(t_max1)*(1p)+mag(t_min1)*p, wherein t_1 represents the first time point, t_max1 represents the first maximum time point, p represents an amplitude ratio parameter, t_min1 represents the first minimum time point, and mag represents an amplitude function.
9. The method as claimed in claim 8, wherein the amplitude ratio parameter p is equal to 80%.
10. The method as claimed in claim 5, wherein extracting the area characteristic value from the respiration signal comprises: calculating a ratio of (i) a first polygon area determined by the first maximum point, the second maximum point, and the first minimum point to (ii) a second polygon area determined by the first minimum point, the second maximum point, and the second minimum point; and extracting the area characteristic value as the calculated ratio.
11. An apparatus for determining a respiration phase of a respiration signal, comprising: extraction circuitry configured to extract a distance characteristic value, a score characteristic value, and an area characteristic value from the respiration signal, the distance characteristic value, the score characteristic value, and the area characteristic value indicating waveform variation between two adjacent waveforms in the respiration signal; and determining circuitry configured to (i) train a respiration signal model according to the distance characteristic value, the score characteristic value, and the area characteristic value, and (ii) determine the respiration phase of the respiration signal using the respiration signal model, wherein a magnetic resonance imaging (MRI) apparatus uses the determined respiration phase of the respiration signal to perform magnetic resonance imaging of an examination region of an examination subject.
12. The apparatus as claimed in claim 11, wherein the extraction circuitry is further configured to: acquire a crest point of a first waveform and a crest point of a second waveform of the two adjacent waveforms of the respiration signal as a first maximum point and a second maximum point, respectively; acquire a trough point between the two adjacent waveforms of the respiration signal as a first minimum point, and to acquire a trough point of the second waveform after the second maximum point as a second minimum point; and determine time points corresponding to the first maximum point, the first minimum point, the second maximum point, and the second minimum point, respectively, as a first maximum time point, a first minimum time point, a second maximum time point, and a second minimum time point.
13. The apparatus as claimed in claim 12, wherein the extraction circuitry is further configured to: calculate a ratio of (i) a difference value between the first minimum time point and the first maximum time point, to (ii) a difference value between the second maximum time point and the first minimum time point; and extract the distance characteristic value as the calculated ratio.
14. The apparatus as claimed in claim 12, wherein the extraction circuitry is further configured to: determine, based on a predetermined amplitude ratio parameter, the first maximum point, and the first minimum point, (i) a start time point of a rising waveform portion leading to the first maximum point of the first waveform, and (ii) a start time point of a falling waveform portion leading to the first minimum point of the first waveform, as (i) a first time point, and (ii) a second time point, respectively, to determine an end-of-expiration start time point of the respiration phase; calculate a ratio of (i) a difference value between the first maximum time point and the first time point, and (ii) a difference value between the first minimum time point and the second time point; and extract the score characteristic value as the calculated ratio.
15. The apparatus as claimed in claim 12, wherein the extraction circuitry is further configured to: calculate a ratio of (i) a first polygon area determined by the first maximum point, the second maximum point, and the first minimum point, and (ii) a second polygon area determined by the first minimum point, the second maximum point, and the second minimum point; and extract the area characteristic value as the calculated ratio.
16. A non-transitory readable storage medium having a program stored thereon that, when executed by one or more processors, cause the one or more processors to: extract a distance characteristic value, a score characteristic value, and an area characteristic value from a respiration signal, the distance characteristic value, the score characteristic value, and the area characteristic value indicating a waveform variation between two adjacent waveforms in the respiration signal; train a respiration signal model according to the distance characteristic value, the score characteristic value, and the area characteristic value; determine a respiration phase of the respiration signal using the respiration signal model; and perform a magnetic resonance imaging of an examination region of an examination subject using the determined respiration phase.
Description
BRIEF DESCRIPTION OF THE DRAWINGS/FIGURES
[0048] The drawings described here are used to provide a better understanding of the present disclosure and form part of the present disclosure; the schematic embodiments of the present disclosure and the descriptions thereof are used to explain the present disclosure, but do not constitute an improper limitation of the present disclosure. In the drawings:
[0049]
[0050]
[0051]
[0052]
[0053]
[0054]
[0055]
[0056]
[0057]
[0058]
DETAILED DESCRIPTION
[0059] To enable those skilled in the art to better understand the present disclosure, embodiments of the present disclosure are described clearly and completely below in conjunction with the drawings accompanying the present disclosure; obviously, the embodiments which are described are merely some, and not all, of the embodiments of the present disclosure. All other embodiments obtained by those skilled in the art on the basis of the embodiments in the present disclosure without expending creative effort shall fall within the scope of protection of the present disclosure.
[0060] According to a natural understanding of respiratory motion, in general, a respiratory motion cycle consists of three phases: aspiration, expiration, and end-of-expiration.
[0061] An exemplary embodiment of the present disclosure provides a method for determining a respiration phase of respiratory motion.
[0062] Step S22: extracting respiration characteristic values from the respiration signal.
[0063] Respiration characteristic values such as a distance characteristic value, a score characteristic value, and an area characteristic value are extracted from the respiration signal. These respiration characteristic values are used to indicate waveform variation between two adjacent waveforms in the respiration signal. Thus, by finding by comparison the waveform variation between two adjacent waveforms, using differences in the shape of aspiration and expiration curves, it is possible to determine an EOE phase in one waveform.
[0064] Step S24: training a respiration signal model according to the extracted respiration characteristic values, and using the respiration signal model to determine a respiration phase of the respiration signal.
[0065] After these respiration characteristic values have been extracted from the respiration signal, e.g. a Pilot Tone respiration signal, they may be used to train a Support Vector Machine (SVM) model. When the SVM model has been trained, this model can be used to determine a respiration phase of the respiration signal.
[0066]
[0067] Step S32: determining extremum points of two adjacent waveforms of the respiration signal and corresponding time points.
[0068]
[0069] First, a crest point of a first waveform and a crest point of a second waveform of the two waveforms are acquired as a first maximum point max1 and a second maximum point max2, respectively. Next, a trough point between the two waveforms is acquired as a first minimum point min1, and a trough point immediately after the second maximum point max2 is acquired as a second minimum point min2.
[0070] Next, time points corresponding to the first maximum point max1, the first minimum point min1, the second maximum point max2, and the second minimum point min2, respectively, are determined, i.e. a first maximum time point t_max1, a first minimum time point t_min1, a second maximum time point t_max2, and a second minimum time point t_min2.
[0071] Step S34: determining a distance characteristic value.
[0072] After the extremum points and the time points corresponding thereto have been determined, the corresponding time points are used to calculate the distance characteristic value.
[0073]
f1=(t_min1t_max1)/(t_max2t_min1).Eqn. 1:
[0074] The formula above may also be expressed as f1=D1/D2, wherein f1 denotes the distance characteristic value.
[0075] Step S36: determining a score characteristic value.
[0076] After the extremum points and the time points corresponding thereto have been determined, the corresponding time points are used to calculate the score characteristic value.
[0077]
[0078] Based on a predetermined amplitude ratio parameter p, the first maximum point max1, and the first minimum point min1, a start time point t_1 (not shown) of a rising waveform portion N1 leading to the first maximum point max1 on the first waveform is determined, and a start time point t_2 (not shown) of a falling waveform portion N2 leading to the first minimum point min1 on the first waveform is determined. As discussed herein, t_1 is referred to as the first time, and t_2_ is referred to as the second time.
[0079] In an exemplary embodiment of the present disclosure, the first time t_1 and the second time t_2 are calculated according to the following two Equations, respectively:
mag(t_1)=mag(t_max1)*p+mag(t_min1)*(1p); andEqn. 2:
mag(t_2)=mag(t_max1)*(1p)+mag(t_min1)*p.Eqn. 3:
[0080] wherein mag(t) is the amplitude of the respiration signal as a function of time t, and p is the predetermined amplitude ratio parameter.
[0081] For example, supposing that the predetermined amplitude ratio parameter p is 80%, then the first time t_1 and second time t_2 should satisfy:
mag(t_1)=mag(t_max1)*0.8+mag(t_min1)*0.2; andEqn. 2:
mag(t_2)=mag(t_max1)*0.2+mag(t_min1)*0.8.Eqn. 3:
[0082] After the possible EOE phase start time points, i.e. the first time point t_1 and second time point t_2, have been determined, the ratio of the difference value between the first maximum time point t_max1 and first time point t_1 to the difference value between the first minimum time point t_min1 and second time point t_2 is calculated, and the result is taken to be the score characteristic value. The following equation is used to calculate the score characteristic value:
f2=(t_max1t_1)/(t_min1t_2).Eqn. 4:
[0083] Step S38: determining an area characteristic value.
[0084] After the extremum points and the time points corresponding thereto have been determined, the corresponding time points are used to calculate the area characteristic value.
[0085]
[0086] In an embodiment of the present disclosure, the following equation is used to calculate the area characteristic value:
f3=s(t_max1,t_max2)/s(t_min1,t_min2)Eqn. 5:
[0087] f3 denotes the area characteristic value, and s denotes the area of the polygon enclosed by the extremum points; in fact, f3 is the ratio S1/S2 of the two areas S1 and S2 shown in
[0088] The embodiments of the present disclosure do not impose restrictions regarding the order of steps S34-S38 above, which may be performed in a different order or performed in parallel (e.g. simultaneously).
[0089] An embodiment of the present disclosure further provides an apparatus for determining a respiration phase of a respiration signal.
[0090] The extraction module 82 is configured to extract respiration characteristic values from the respiration signal for the purpose of training an SVM model. These respiration characteristics include but are not limited to, e.g. a distance characteristic value, a score characteristic value, and an area characteristic value.
[0091] In an exemplary embodiment of the present disclosure, the distance characteristic value is calculated by the following formula:
f1=(t_min1t_max1)/(t_max2t_min1).Eqn. 6:
[0092] f1 denotes the distance characteristic value, t_min1 denotes a trough point between two adjacent waveforms in the respiration signal, t_max1 denotes a crest point of a first waveform, and t_max2 denotes a crest point of a second waveform.
[0093] In an exemplary embodiment of the present disclosure, suppose that the time points of 80% of the amplitudes of the maximum point and minimum point are a first time t_1 and a second time t_2 respectively, i.e. the following conditions are met:
mag(t_1)=mag(t_max1)*0.8+mag(t_min1)*0.2 and mag(t_2)=mag(t_max1)*0.2+mag(t_min1)*0.8.
[0094] After the first time t_1 and second time t_2 have been calculated by the above equations, the formula f2=(t_max1t_1)/(t_min1t_2) is used to calculate the score characteristic value, and the equation f3=s(t_max1, t_max2)/s(t_min1, t_min2) is used to calculate the area characteristic value, wherein f2 denotes the score characteristic value, f3 denotes the area characteristic value, s(t1, t2) is defined as the area of a closed polygon formed between the waveform and all time points from time t1 to t2, e.g. s(t_max1, t_max2) denotes the area of the closed polygon formed by reach of the time points from time t_max1 to t_max2, and s(t_min1, t_min2) denotes the area of the closed polygon formed by each of the time points from time t_min1 to t_min2.
[0095] After these respiration characteristic values have been extracted from the respiration signal, the determining module 84 uses these respiration characteristic values to train the SVM model, i.e. a respiration signal model. When the SVM model has been trained, the model can be used to determine a respiration phase.
[0096]
[0097] A schematic embodiment of the present disclosure further provides an MRI method. The MRI method is an image outputting method. First, the method for determining a respiration phase of a respiration signal according to the present disclosure is used to determine a respiration phase. Second, magnetic resonance imaging of an examination region of an examination subject is performed according to the determined respiration phase. Specifically, an MRI system is triggered to perform a scan in the EOE phase of the respiration signal, thereby avoiding motion artifacts in the MR image, and it is thus possible to generate an image of living tissue with high precision.
[0098] A schematic embodiment of the present disclosure further provides an MRI system, comprising the apparatus for determining a respiration phase provided in the present disclosure and an imaging apparatus, which is not shown in the Figures for purposes of brevity but may operate in accordance with known MR imaging systems. For instance, after an examination subject has been placed in the MRI system, the apparatus for determining a respiration phase provided in the present disclosure is used to determine a respiration phase of a respiration signal, in particular an EOE phase. Next, the MRI system is triggered at a start time point of the EOE phase to transmit an MR signal and perform a scan, and gradient coils are used to modify an external magnetic field, i.e. a body layer of the examination subject is selected while generating MR signal region encoding. An MR signal is then reconstructed, e.g. by means of a Fourier transform, to generate an image of the selected body layer, for the purpose of medical diagnosis. MR signal generation and detection is achieved using a high-frequency system; this system comprises a transmitting antenna for sending high-frequency excitation pulses into the body of the examination subject, and a receiving antenna for receiving resonance signals modulated by the examination subject.
[0099]
[0100] It should be noted that the CPU 1010 may comprise one or more processors; the one or more processors and/or other data processing circuits can in general be referred to as apparatuses for determining a respiration phase in the present disclosure. The data processing circuit may be completely or partially embodied as software, hardware, firmware or any other combination. In addition, the data processing circuit may be a single independent processing module, or completely or partially incorporated in any one of the other components in the computing apparatus 100.
[0101] The storage unit 1040 may be used to store software programs and modules of application software, such as a program instruction/data storage apparatus corresponding to a main control instruction computing method described in the present disclosure; the CPU 1010 realizes the abovementioned main control instruction computing method by running the software programs and modules stored in the storage unit 1040. The storage unit 1040 may comprise a non-volatile memory, such as one or more magnetic storage apparatuses, internal memories, or other non-volatile solid-state memories. In some practical examples, the storage unit 1040 may further comprise a memory disposed remotely from the CPU 1010; these remote memories may be connected to the computing apparatus 100 via a network. Real examples of the abovementioned network include but are not limited to the internet, intranets, local area networks, mobile communication networks and combinations thereof.
[0102] Again, an embodiment of the present disclosure further provides a non-transitory readable storage medium having a program stored thereon and, when the program is executed, a processor is caused to execute the method provided in the present disclosure for determining a respiration phase of a respiration signal. The non-transitory readable storage medium may be identified with the storage unit 1040 or a different storage unit not shown in the Figures, and the processor may be identified with the CPU 1010 or a different processor not shown in the Figures. The computing apparatus 100 may form part of and/or be in communication with a magnetic resonance scanning system, such as a control computer or other computing components that may be used by a MRI scanner to perform MR imaging of an examination region of an examination subject according to the determined respiration phase. Thus, the non-transitory readable storage medium may store instructions that are executed by the CPU 1010 or other suitable components of a MRI scanner or apparatus to perform the methods as discussed herein, which include, once determined, using the determined respiration phase to perform an MR imaging of an examination region of an examination subject.
[0103] The communication unit 1060 is used for receiving or sending data via a network. Specific examples of the abovementioned network may include a wireless network provided by a communication supplier of the computing apparatus 100. In an example, the communication unit 1060 comprises a network adapter (Network Interface Controller, NIC), which may be connected to other network equipment via a base station and can communicate with the internet. In an example, the communication unit 1060 may be a Radio Frequency (RF) module used to communicate with the Internet wirelessly.
[0104] In an embodiment of the present disclosure, a respiration signal model is trained on the basis of a machine learning method such as an SVM method by extracting characteristic values of a respiration signal such as a Pilot Tone respiration signal, and a respiration phase of the respiration signal is further determined on the basis of the respiration signal model, resulting in the beneficial effects of robustness, high performance, and support for Pilot Tone applications.
[0105] Only example embodiments of the present disclosure are described above. It should be pointed out that those skilled in the art can make improvements and modifications without departing from the principle of the present disclosure, and that these improvements and modifications should also fall within the scope of protection of the present disclosure.
[0106] The various functional blocks, apparatuses, modules, units, components of physical or functional units, etc., as shown in the drawings and described herein may be implemented via any suitable number and type of computer processors, hardware components, the execution of software algorithms, or combinations thereof, and thus may alternatively be referred to as a unit, system, circuitry, processor(s), or a device.