Phase segmentation method, device and computer equipment for magnetic resonance electrocardiogram gating

Abstract

The present application relates to a magnetic resonance electrocardiogram (MRI) gating phase segmentation method, apparatus, and computer equipment, wherein the method comprises: acquiring first MRI waveform data, performing phase analysis on the first MRI waveform data, and determining phase analysis parameters; synchronously and continuously acquiring second MRI waveform data and corresponding MRI data; performing phase segmentation on at least one second cycle in the second MRI waveform data according to the phase analysis parameters, and determining the phase encoding corresponding to each second cycle; encoding the synchronously acquired MRI data according to the phase encoding to obtain target images of different phases. Through the present application, the problem of artifacts generated when the heart rate of the acquired subject is irregular is solved, the generation of artifacts is reduced, and the quality of the target image is improved.

Description

Technical Field

[0001] The present application relates to the field of medical imaging technology, and in particular to a method, device and computer equipment for phase segmentation of magnetic resonance electrocardiogram (MRI) gating. Background Art

[0002] With the development of medical imaging technology, various medical imaging devices are widely used in clinical practice to perform medical imaging on subjects for diagnosis and treatment. The heart is an important component of the human body and a rapidly moving organ. Physiologically, different moments in the cardiac cycle or different states of motion of the heart are referred to as phases. When imaging at different phases of the cardiac cycle, images of the heart are obtained under different working conditions. For example, imaging at the end of systole will obtain an image of the heart in full contraction; imaging at the end of diastole will obtain an image of the heart in full relaxation. Based on the end-diastolic and end-systolic volumes of the heart, important cardiac hemodynamic indicators such as ejection fraction can be calculated, providing a basis for evaluating cardiac function. This is a basic requirement for clinical cardiac imaging.

[0003] MRI ECG gating uses the ECG R wave as a trigger, starting data acquisition after a certain period of time. This "certain period" is called the time of delay (TD), and its value is affected by the subject's heart rate.

[0004] Currently, in retrospectively acquired MRI ECG-gated imaging schemes, a large amount of data is collected non-selectively and continuously over multiple cardiac cycles, and then the same level of each cardiac cycle is combined to form an image. For example, within each cardiac cycle, images are continuously collected for eight time periods, and then the data for the first time period of each cardiac cycle is filled into K-space a, and the data for the second time period of each cardiac cycle is filled into K-space b, and so on. Ultimately, eight images of different time periods of the cardiac cycle are obtained, and continuous playback can form a movie image, showing the heart movement during a cardiac cycle. However, when the subject's heart rate is irregular, the R wave triggering time is different. After the same delay time (TD), the position and morphology of the heart change, resulting in artifacts.

[0005] Regarding the problem in related technologies that artifacts may be generated when the heart rate of the subject being collected is irregular, no effective solution has been proposed so far. Summary of the Invention

[0006] In this embodiment, a magnetic resonance electrocardiogram (MRI) gating phase segmentation method, apparatus, and computer device are provided to solve the problem in related technologies that artifacts may be generated when the heart rate of the subject being acquired is irregular.

[0007] In a first aspect, this embodiment provides a method for phase segmentation of magnetic resonance electrocardiography gating, comprising:

[0008] collecting first electrocardiogram waveform data, performing time phase analysis on the first electrocardiogram waveform data, and determining time phase analysis parameters;

[0009] synchronously and continuously acquiring second electrocardiogram waveform data and corresponding magnetic resonance data;

[0010] Performing time phase segmentation on at least one second cycle in the second ECG waveform data according to the time phase analysis parameter, and determining a time phase code corresponding to each second cycle; wherein the second cycle is a period between two adjacent first target waves in the second ECG waveform data;

[0011] The synchronously acquired magnetic resonance data are coded according to the phase coding, and the magnetic resonance data with the same phase coding are reconstructed to obtain target images of different phases.

[0012] In some embodiments, the collecting of first ECG waveform data and performing phase analysis on the first ECG waveform data includes:

[0013] collecting first electrocardiogram waveform data, performing first target wave detection on the first electrocardiogram waveform data, and determining a first period between two adjacent first target waves in the first electrocardiogram waveform data;

[0014] performing a second target wave detection on the electrocardiogram data corresponding to the first cycle;

[0015] The first cycle is divided based on the detection result of the second target wave to obtain a first interval matrix.

[0016] In some embodiments, determining the phase analysis parameters includes:

[0017] Each waveform parameter in the first interval matrix is analyzed based on the first cycle to determine a phase analysis parameter corresponding to each waveform parameter in the first interval matrix.

[0018] In some embodiments, the first target wave is an R wave; the second target wave is a T wave and a P wave; and the first cycle is a first RR interval.

[0019] In some embodiments, the method further comprises:

[0020] After synchronously and continuously acquiring the second ECG waveform data and the corresponding magnetic resonance data, first target wave detection is performed on the second ECG waveform data to determine a second period between every two adjacent first target waves in the second ECG waveform data.

[0021] In some embodiments, performing phase segmentation on at least one second cycle in the second ECG waveform data according to the phase analysis parameter and determining a phase code corresponding to each second cycle includes:

[0022] According to the phase analysis parameters, at least one second cycle in the second electrocardiogram waveform data is segmented to obtain a phase interval matrix, and the phase interval matrix is calculated to obtain a phase code corresponding to the second cycle.

[0023] In some embodiments, encoding the synchronously acquired magnetic resonance data according to the phase encoding, and reconstructing the magnetic resonance data with the same phase encoding to obtain target images of different phases, includes:

[0024] Determine the coding range based on preset requirements;

[0025] According to the phase coding, the synchronously acquired magnetic resonance data within the coding range are coded, and the magnetic resonance data with the same phase coding are reconstructed to obtain target images of different phases.

[0026] In some embodiments, encoding the synchronously acquired magnetic resonance data within the encoding range according to the phase encoding, and reconstructing the magnetic resonance data with the same phase encoding to obtain target images of different phases, includes:

[0027] Based on the coding range, corresponding phase codes are assigned to the magnetic resonance data of different time periods, and the magnetic resonance data with the same phase code are filled into the corresponding position of the K space for image reconstruction to obtain target images of different phases.

[0028] In a second aspect, a magnetic resonance electrocardiogram (MRI)-gated phase segmentation device is provided in this embodiment, comprising: a first acquisition module, a second acquisition module, a first processing module, and a second processing module;

[0029] The first acquisition module is configured to acquire first ECG waveform data, perform time phase analysis on the first ECG waveform data to obtain a first interval matrix, and determine a time phase analysis parameter based on the first interval matrix;

[0030] The second acquisition module is used to synchronously and continuously acquire second ECG waveform data and corresponding magnetic resonance data;

[0031] The first processing module is configured to perform time phase segmentation on at least one second cycle in the second ECG waveform data according to the time phase analysis parameter, and determine a time phase code corresponding to each second cycle; wherein the second cycle is a period between two adjacent first target waves in the second ECG waveform data;

[0032] The second processing module is configured to perform encoding processing on the synchronously acquired magnetic resonance data according to the phase encoding, and to perform image reconstruction on the magnetic resonance data with the same phase encoding to obtain target images of different phases.

[0033] In a third aspect, a computer device is provided in this embodiment, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the magnetic resonance electrocardiogram (MRI) gating phase segmentation method described in the first aspect is implemented.

[0034] In a fourth aspect, a storage medium is provided in this embodiment, on which a computer program is stored. When the program is executed by a processor, the phase segmentation method of magnetic resonance electrocardiography gating described in the first aspect is implemented.

[0035] Compared with the related art, the magnetic resonance electrocardiogram gating phase segmentation method, device and computer equipment provided in this embodiment collect first electrocardiogram waveform data, perform phase analysis on the first electrocardiogram waveform data, and determine the phase analysis parameters; synchronously and continuously collect second electrocardiogram waveform data and corresponding magnetic resonance data; perform phase segmentation on at least one second cycle in the second electrocardiogram waveform data according to the phase analysis parameters, and determine the phase encoding corresponding to each second cycle; encode the synchronously collected magnetic resonance data according to the phase encoding, and reconstruct the magnetic resonance data with the same phase encoding to obtain target images of different phases, which solves the problem of artifacts generated when the heart rate of the collected object is irregular, reduces the generation of artifacts, and can improve the quality of the target image.

[0036] The details of one or more embodiments of the present application are set forth in the following drawings and description to make other features, objects, and advantages of the present application more readily apparent. BRIEF DESCRIPTION OF THE DRAWINGS

[0037] The drawings described herein are used to provide a further understanding of the present application and constitute a part of the present application. The illustrative embodiments of the present application and their descriptions are used to explain the present application and do not constitute an improper limitation on the present application. In the drawings:

[0038] Figure 1 This is a hardware structure block diagram of a terminal device for a magnetic resonance electrocardiogram (MRI) gating phase segmentation method provided in one embodiment of the present application;

[0039] Figure 2 This is a flow chart of a method for phase segmentation of magnetic resonance electrocardiography gating provided in one embodiment of the present application;

[0040] Figure 3 is a schematic diagram of electrocardiogram waveform data provided by an embodiment of the present application;

[0041] Figure 4 1 is a flow chart of a method for phase segmentation of magnetic resonance electrocardiography gating provided in one embodiment of the present application;

[0042] Figure 5 This is a structural block diagram of a magnetic resonance electrocardiogram (MRI) gated phase segmentation device provided in one embodiment of the present application.

[0043] In the figure: 210, a first acquisition module; 220, a second acquisition module; 230, a first processing module; 240, a second processing module. DETAILED DESCRIPTION

[0044] In order to more clearly understand the purpose, technical solutions and advantages of the present application, the present application is described and illustrated below in conjunction with the accompanying drawings and embodiments.

[0045] Unless otherwise defined, the technical terms or scientific terms involved in this application should have the general meaning understood by people with ordinary skills in the technical field to which this application belongs. The words "one", "an", "a", "the", "these" and the like in this application do not indicate quantitative restrictions, and they can be singular or plural. The terms "include", "comprise", "have" and any variants thereof involved in this application are intended to cover non-exclusive inclusions; for example, a process, method and system, product or device comprising a series of steps or modules (units) is not limited to the listed steps or modules (units), but may include unlisted steps or modules (units), or may include other steps or modules (units) inherent to these processes, methods, products or devices. The words "connect", "connected", "coupled" and the like involved in this application are not limited to physical or mechanical connections, but may include electrical connections, whether direct or indirect. The "plurality" involved in this application refers to two or more. "And / or" describes the relationship between related objects, indicating that three possible relationships exist. For example, "A and / or B" can represent: A exists alone, A and B exist simultaneously, and B exists alone. Generally, the character " / " indicates that the related objects are in an "or" relationship. The terms "first," "second," "third," etc. used in this application are only used to distinguish similar objects and do not represent a specific ordering of the objects.

[0046] The method embodiment provided in this embodiment can be executed in a terminal, a computer or a similar computing device. For example, running on a terminal, Figure 1 FIG. 1 is a block diagram of the hardware structure of the terminal of the magnetic resonance electrocardiogram gating phase segmentation method of this embodiment. Figure 1 As shown, the terminal may include one or more ( Figure 1 Only one is shown) a processor 102 and a memory 104 for storing data, wherein the processor 102 may include but is not limited to a processing device such as a microprocessor MCU or a programmable logic device FPGA. The above terminal may also include a transmission device 106 and an input and output device 108 for communication functions. It will be understood by those skilled in the art that Figure 1 The structure shown is only for illustration and does not limit the structure of the above terminal. Figure 1 More or fewer components than shown, or with Figure 1 Different configurations shown.

[0047] The memory 104 can be used to store computer programs, for example, software programs and modules of application software, such as the computer program corresponding to the phase segmentation method for MRI ECG gating in this embodiment. The processor 102 executes the computer programs stored in the memory 104 to perform various functional applications and data processing, thereby implementing the above-mentioned method. The memory 104 may include high-speed random access memory and may also include non-volatile memory, such as one or more magnetic storage devices, flash memory, or other non-volatile solid-state memory. In some examples, the memory 104 may further include memory remotely located relative to the processor 102, and these remote memories may be connected to the terminal via a network. Examples of such networks include, but are not limited to, the Internet, corporate intranets, local area networks, mobile communication networks, and combinations thereof.

[0048] The transmission device 106 is used to receive or send data via a network. The network may include a wireless network provided by the terminal's telecommunications provider. In one embodiment, the transmission device 106 includes a network interface controller (NIC), which can be connected to other network devices via a base station to enable communication with the Internet. In another embodiment, the transmission device 106 may be a radio frequency (RF) module, which is used to communicate with the Internet wirelessly.

[0049] In this embodiment, a magnetic resonance electrocardiogram (MRI) gating phase segmentation method is provided. Figure 2 FIG. 1 is a flow chart of the magnetic resonance electrocardiogram (MRI) gating phase segmentation method of this embodiment, as shown in FIG. Figure 2 As shown, the process includes the following steps:

[0050] Step S210, collecting first ECG waveform data, performing time phase analysis on the first ECG waveform data, and determining time phase analysis parameters;

[0051] Step S220, synchronously and continuously acquiring second ECG waveform data and corresponding magnetic resonance data;

[0052] Step S230, performing time phase segmentation on at least one second cycle in the second ECG waveform data according to the time phase analysis parameter, and determining a time phase code corresponding to each second cycle; wherein the second cycle is the period between two adjacent first target waves in the second ECG waveform data;

[0053] Step S240 , encoding the synchronously acquired magnetic resonance data according to the phase encoding, and reconstructing the magnetic resonance data with the same phase encoding to obtain target images of different phases.

[0054] It should be noted that if Figure 3 As shown, a normal ECG waveform data for a cardiac cycle consists of the P wave, QRS complex, and T wave. In addition to heart rate, cardiac activity is more importantly reflected in the P, QRS, and T waves within the ECG waveform data, as well as their various parameters. The P wave, QRS complex, and T wave are the most important characteristic waves. These, along with the PR interval, QT interval, and ST interval formed from them, are the most important ECG characteristic signals, reflecting multiple aspects of the cardiac conduction system and the presence of pathological changes in the heart itself. The QRS complex reflects the potential changes during ventricular myocardial depolarization. A typical QRS complex consists of three closely connected waves: the first downward wave is called the Q wave (low voltage, called the q wave); the first upward wave is called the R wave (low voltage, called the r wave); and the downward wave following the R (or r) wave is called the S wave (low voltage, called the s wave). These three waves are closely connected, with a total duration of no more than 0.1 second. Since they all reflect ventricular activation, they are collectively referred to as the QRS complex. The common QRS complex may present various waveforms such as qRs, qR, Qr, RS, rS, QS, Rs, etc. in different leads.

[0055] In this embodiment, it can be divided into two processes.

[0056] One process is to determine the phase analysis parameters through phase analysis during the self-learning process.

[0057] Specifically, the first ECG waveform data can be obtained by collecting the subject by an ECG acquisition device, and the first ECG waveform data can include multiple cardiac cycles; however, in order to reduce the impact of data differences between multiple cardiac cycles on the results, the preferred first ECG waveform data only includes one cardiac cycle. Using algorithms such as wavelet algorithms, pattern matching algorithms, and artificial neural network algorithms, the first ECG waveform data is subjected to phase analysis to obtain a first interval matrix. Different first interval matrices can be divided based on the detection of different first target waves. Preferably, taking the R wave as the first target wave as the detection object, it can be considered that the first interval matrix is composed of the RS interval, ST interval, T wave, TP interval, P wave, PQ interval, and QR interval; the first interval matrix can be expressed as: [RS ST T TP P PQ QR]. For another example, if the T wave is the first target wave to be detected, the first interval matrix can be considered to be composed of the T wave, TP interval, P wave, PQ interval, QR interval, RS interval, and ST interval; the first interval matrix can be expressed as: [RS ST T TP P PQ QR RS ST]. There is no limitation on the division of the first interval matrix and the second interval matrix described below, and examples are not given here one by one.

[0058] The phase analysis parameters can be determined based on the first interval matrix. For example, the first interval matrix can be transformed using the preset correction parameters and the first cycle in the first ECG waveform data to obtain the phase analysis parameters. The preset correction parameters and the first cycle have corresponding correction formulas for each waveform interval in the first interval matrix. For example, the correction formula can be Xindex = X / RR 1 / 2 ; The correction formula can be Xindex=X / RR 1 / 3 ; The correction formula may be Xindex=X+0.154(1-RR); The correction formula may be Xindex=X+0.00175([60 / RR]-60); and other correction formulas. Wherein, X is the interval parameter of each waveform interval in the first interval matrix before the transformation; Xindex represents the parameter of each waveform interval in the first interval matrix after the transformation, that is, the phase analysis parameter. Phase analysis can obtain appropriate phase analysis parameters according to the different adjustments of each person's heart rate (that is, standardize different heart rates), eliminate the influence of one's own heart rate changes, and thereby reduce the artifacts of the target image. In other embodiments, there is no limitation on the correction formula.

[0059] The other process is the retrospective acquisition process, in which the phase analysis parameters are applied to the retrospective acquisition process to ultimately obtain target images at different phases, thereby reducing the generation of artifacts and improving the quality of the target image.

[0060] Specifically, during the retrospective acquisition process, it is necessary to synchronously and continuously acquire the second ECG waveform data and the corresponding magnetic resonance data. The magnetic resonance data is a magnetic resonance image, and the optional magnetic resonance data may also be K-space data, etc. The second ECG waveform data may span multiple cardiac cycles, and the second ECG waveform data of each cardiac cycle may also be composed of the RS interval, ST interval, T wave, TP interval, P wave, PQ interval, and QR interval. The acquired magnetic resonance data may be collected at different time periods, such as: eight consecutive time periods between the ST interval and the TP interval, eight consecutive time periods between the RR interval, etc.

[0061] Based on the phase analysis parameters determined during the self-learning process and the second cycle between two adjacent first target waves in the second ECG waveform data, at least one second cycle in the second ECG waveform data is phase-segmented according to the phase analysis parameters, and a phase code corresponding to each second cycle is determined. Phase segmentation is required for each second cycle acquired continuously. For example, taking the first target wave as the R wave and the second cycle as the second RR interval, the phase code for the T wave in the second RR interval can be: Tindex - 0.2 * (1000 - RR), where 0.2 and 1000 are empirical values that can be adjusted based on the characteristics of the subject. The phase code for each waveform in the second RR interval can be derived using a formula similar to the above. The phase code can be obtained using the phase analysis parameters and the second RR interval. The synchronously acquired magnetic resonance data is then encoded based on the obtained phase code to obtain target images at different phases.

[0062] According to the above steps, by acquiring first ECG waveform data, performing phase analysis on the first ECG waveform data, and determining phase analysis parameters; synchronously and continuously acquiring second ECG waveform data and corresponding magnetic resonance data; performing phase segmentation on at least one second cycle in the second ECG waveform data according to the phase analysis parameters, and determining the phase encoding corresponding to each second cycle; encoding the synchronously acquired magnetic resonance data according to the phase encoding to obtain target images of different phases; utilizing the phase analysis parameters to match the time when the subject's R wave is triggered when the heart rate is irregular, thereby solving the problem of artifacts generated when the subject's heart rate is irregular because the position and morphology of the heart change after the same delay time (TD), thereby reducing the generation of artifacts and improving the quality of the target image.

[0063] Combine Figure 4 , the above steps are described in detail below:

[0064] In the application, when performing phase analysis on the first ECG waveform data, the first target wave detection included and the process of performing first target wave detection on the second ECG waveform data can be the same. Here, the phase analysis on the first ECG waveform data is described in detail.

[0065] In some embodiments, collecting the first ECG waveform data and performing phase analysis on the first ECG waveform data in step S210 includes the following steps:

[0066] Step S211, collecting first ECG waveform data, performing first target wave detection on the first ECG waveform data, and determining a first period between two adjacent first target waves in the first ECG waveform data;

[0067] Step S212, performing a second target wave detection on the ECG data corresponding to the first cycle;

[0068] Step S213 , dividing the first period based on the detection result of the second target wave to obtain a first interval matrix.

[0069] Specifically, the first target wave can be any specified wave among the P wave, Q wave, R wave, S wave, or T wave. Depending on the specified wave, the first cycle will also be adaptively adjusted. The second target wave can be the T wave or the P wave, and the second target wave can also be adjusted to other waves in the first cycle based on the first target wave.

[0070] Taking the first target wave as the R wave; the second target wave as the T wave and the P wave; and the first cycle as the first RR interval as an example, the above process can be: collecting the first ECG waveform data, performing R wave detection on the first ECG waveform data, and determining the first RR interval between two adjacent R waves in the first ECG waveform data; performing T wave and P wave detection on the first RR interval; segmenting the first RR interval based on the detection results of the T wave and the P wave to obtain a first interval matrix.

[0071] In which, an R-wave detection algorithm is used to perform R-wave detection on the first ECG waveform data to detect the position of the R wave in each cardiac cycle; then, based on the positions of two adjacent R waves in the first ECG waveform data, the first RR interval can be quickly and accurately detected. In which, the R-wave detection algorithm can be a Pan-Tompkins algorithm, which can not only detect the position of the R wave, but also detect parameters related to the QRS complex. In other embodiments, other R-wave detection algorithms can also be used, such as: based on a binary spline 4-layer wavelet transform, using the maximum and minimum method in the 3-layer detail coefficients to detect the position of the R wave in each cardiac cycle; there is no limitation on other R-wave detection algorithms.

[0072] Using a detection algorithm similar to the R-wave detection algorithm, T-wave and P-wave detection can be performed on the first RR interval, determining the positions of the T and P waves. The first RR interval is then segmented based on the positions of the T and P waves, resulting in a first interval matrix, which can be expressed as: [RS ST T TP P PQ QR]. This allows accurate detection of R-wave and other waveform parameters even when an irregular heartbeat results in varying R-wave triggering times.

[0073] For example, R wave detection is performed on the collected first ECG waveform data, and the first RR interval of the first ECG wave is obtained to be 1000ms. Then T wave and P wave detection are performed, and the first interval matrix after segmentation is obtained as [RS ST T TP P PQQR] = [50ms 100ms 300ms 300ms 100ms 100ms 50ms].

[0074] In some embodiments, determining the phase analysis parameters in step S210 includes the following steps:

[0075] Each waveform parameter in the first interval matrix is analyzed based on the first RR interval, and a phase analysis parameter corresponding to each waveform parameter in the first interval matrix is determined.

[0076] Specifically, each waveform parameter in the first interval matrix corresponds to a phase analysis parameter, which can be obtained by analyzing each waveform parameter in the first interval matrix based on the first RR interval. Each waveform parameter has a corresponding phase analysis parameter (index); for example, the phase analysis parameter RSindex of the RS interval is RS itself. The phase analysis parameter STindex of the ST interval is ST+0.1*(1000-RR1), where RR1 is the first RR interval; the phase analysis parameter Tindex of the T wave is T+0.2*(1000-RR1); the phase analysis parameter TP index of the TP interval is TP+0.5*(1000-RR1); the phase analysis parameter Pindex of the P wave is P+0.1*(1000-RR1); the phase analysis parameter PQindex of the PQ interval is PQ+0.1*(1000-RR1); the phase analysis parameter QRindex of the QR interval is QR itself. Similar formulas, where 0.1, 0.2, 0.5 and 1000 are empirical values, can be adjusted according to the nature of the subject. Then it can be considered that the phase analysis parameter matrix index corresponding to the first interval matrix can be expressed as: [RSindex STindex Tindex TPindex Pindex PQindeQRindex]=[RS ST+0.1*(1000-RR1) T+0.2*(1000-RR1) TP+0.5*(1000-RR1) P+0.1*(1000-RR1) PQ+0.1*(1000-RR1) QR]. For example, the phase analysis parameter matrix index corresponding to the first interval matrix obtained by the above segmentation is: [RSindex STindex Tindex TPindex Pindex PQinde QRindex] = [50ms 100ms+0.1*(1000-1000ms) 300ms+0.2*(1000-1000ms) 300ms+0.5*(1000-1000ms)100ms+0.1*(1000-1000ms) 100ms+0.1*(1000-1000ms) 50ms] = [50ms 100ms 300ms300ms 100ms 100ms 50ms]. This embodiment can ensure that each waveform parameter of the first interval matrix in each cardiac cycle matches the phase analysis parameter without involving a complex algorithm process.

[0077] In some embodiments, the magnetic resonance electrocardiography gating phase segmentation method further includes the following steps:

[0078] After synchronously and continuously acquiring the second ECG waveform data and the corresponding magnetic resonance data, first target wave detection is performed on the second ECG waveform data to determine a second period between every two adjacent first target waves in the second ECG waveform data.

[0079] The process of performing the first target wave detection on the second ECG waveform data can be the same as the above process, which will not be repeated here.

[0080] In some embodiments, step S230 of performing phase segmentation on at least one second cycle in the second ECG waveform data according to the phase analysis parameter and determining the phase code corresponding to each second cycle includes the following steps:

[0081] Step S231 , dividing at least one second cycle in the second ECG waveform data according to the phase analysis parameters to obtain a phase interval matrix, and calculating the phase interval matrix to obtain a phase code corresponding to the second cycle.

[0082] Specifically, since the length of the collected cycles may vary in each cardiac cycle, this embodiment uses the phase analysis parameters to segment at least one second cycle in the second ECG waveform data, thereby accurately segmenting cycles of different lengths and obtaining a phase interval matrix. The phase interval matrix is calculated to obtain a phase code corresponding to the second cycle for rapid encoding processing.

[0083] Segmentation refers to segmenting the second period into corresponding time phase interval matrices using time phase analysis parameters after obtaining the periodic data of the second period, thereby improving segmentation efficiency and ensuring the accuracy of the segmentation results.

[0084] The phase code can be obtained by segmenting the second period, that is, the elements in the phase interval matrix. For example, the phase code corresponding to STnew is STindex-0.1*(1000-RR).

[0085] Take the second cycle as the second RR interval as an example to illustrate:

[0086] R wave detection is performed on the second RR interval of the second cycle, and the period of the second RR interval is obtained as RR2. The phase code is the elements in the phase interval matrix obtained by segmenting the second RR interval. The phase interval matrix can be expressed as: [RSnew STnew Tnew TPnew Pnew PQnew QRnew] = [RSindex STindex-0.1*(1000-RR2)Tindex-0.2*(1000-RR2) TPindex-0.5*(1000-RR2) Pindex-0.1*(1000-RR2) PQinde-0.1*(1000-RR2) QRindex]. For example, according to the above-mentioned phase analysis parameter matrix index: [50ms 100ms300ms 300ms 100ms 100ms 50ms], the second RR interval of RR2=800ms is divided, and the obtained phase interval matrix is: [RSnew STnew Tnew TPnew Pnew PQnew QRnew]=[50ms 100ms-0.1*(1000-800ms) 300ms-0.2*(1000-800ms) 300ms-0.5*(1000-800ms) 100ms-0.1*(1000-800ms)100ms-0.1*(1000-800ms) 50ms]=[50ms 80ms 260ms 200ms 80ms 80ms 50ms].

[0087] In some embodiments, step S240 of encoding the synchronously acquired magnetic resonance data of different time periods according to the phase coding, and reconstructing the magnetic resonance data with the same phase coding to obtain target images of different phases includes the following steps:

[0088] Determine the coding range based on preset requirements;

[0089] According to the phase coding, the synchronously acquired magnetic resonance data within the coding range are coded and processed, and the magnetic resonance data with the same phase coding are image reconstructed to obtain target images of different phases.

[0090] The encoding range is determined based on the preset requirements. For example, if the preset requirement is the ST to TP interval, the encoding range is: multiple consecutive periods from the ST to TP interval. If the preset requirement is the RR interval, the encoding range is: multiple consecutive periods of the RR interval. The multiple consecutive periods here can be 8 consecutive periods, or any other number of periods, without limitation.

[0091] According to the phase coding, the synchronously acquired magnetic resonance data within the above-mentioned coding range is coded and processed, thereby assigning corresponding phase coding to the magnetic resonance data of different time periods; the magnetic resonance data with the same phase coding are image reconstructed to obtain target images of different phases.

[0092] Since the magnetic resonance data can be magnetic resonance images, K-space data, etc.; Figure 4 As shown, taking the magnetic resonance data as magnetic resonance images as an example, the magnetic resonance images are acquired synchronously and continuously; according to the phase encoding and the encoding range, the magnetic resonance images acquired synchronously and located within the encoding range are subjected to image phase encoding, and the magnetic resonance images with the same phase encoding are reconstructed to obtain target images of different phases.

[0093] Specifically, based on the above coding range, corresponding phase codes can be assigned to MRI images of different time periods, and the data corresponding to the MRI images with the same phase codes can be filled into the corresponding positions of the K space for image reconstruction to obtain target images of different phases.

[0094] For example: for 8 consecutive time periods with a coding range of RR intervals, a corresponding time phase code (RSindex) is assigned to the first time period, a corresponding time phase code (STindex-0.1*(1000-RR2)) is assigned to the second time period, a corresponding time phase code (Tindex-0.2*(1000-RR2)) is assigned to the third time period, a corresponding time phase code (TPindex-0.5*(1000-RR2)) is assigned to the fourth and fifth time periods, a corresponding time phase code (Pindex-0.1*(1000-RR2)) is assigned to the sixth time period, a corresponding time phase code (PQinde-0.1*(1000-RR2)) is assigned to the seventh time period, and a corresponding time phase code (QRindex) is assigned to the eighth time period. The data corresponding to the magnetic resonance image with the same phase encoding is then filled into the corresponding position of the K space for image reconstruction to obtain target images of different phases; phase encoding is used to reduce the impact of arrhythmia on imaging, and the target image can be adjusted based on the encoding range, which is convenient for use.

[0095] It should be noted that the steps shown in the above process or the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0096] This embodiment also provides a magnetic resonance electrocardiogram (MRI) gated phase segmentation device, which is used to implement the above-mentioned embodiments and preferred embodiments. Details already described will not be repeated here. The terms "module," "unit," "subunit," etc. used below may refer to a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation using hardware, or a combination of software and hardware, is also possible and contemplated.

[0097] Figure 5 This is a structural block diagram of the magnetic resonance electrocardiogram gated phase segmentation device of this embodiment, as shown in FIG. Figure 5 As shown, the device includes: a first acquisition module 210, a second acquisition module 220, a first processing module 230 and a second processing module 240;

[0098] A first acquisition module 210 is configured to acquire first ECG waveform data, perform time phase analysis on the first ECG waveform data, and determine time phase analysis parameters;

[0099] A second acquisition module 220 is used to synchronously and continuously acquire second ECG waveform data and corresponding magnetic resonance data;

[0100] a first processing module 230 configured to perform time phase segmentation on at least one second cycle in the second ECG waveform data according to the time phase analysis parameter, and determine a time phase code corresponding to each second cycle; wherein the second cycle is a period between two adjacent first target waves in the second ECG waveform data;

[0101] The second processing module 240 is configured to perform encoding processing on the synchronously acquired magnetic resonance data according to the phase encoding, and to perform image reconstruction on the magnetic resonance data with the same phase encoding to obtain target images of different phases.

[0102] The above device solves the problem of artifacts generated when the heart rate of the collected object is irregular, reduces the generation of artifacts, and improves the quality of the target image.

[0103] In some embodiments, the first acquisition module 210 is further configured to acquire first ECG waveform data, perform first target wave detection on the first ECG waveform data, and determine a first period between two adjacent first target waves in the first ECG waveform data;

[0104] Performing a second target wave detection on the electrocardiogram data corresponding to the first cycle;

[0105] The first period is segmented based on the detection result of the second target wave to obtain a first interval matrix.

[0106] In some embodiments, the first acquisition module 210 is further configured to analyze each waveform parameter in the first interval matrix based on the first cycle, and determine a phase analysis parameter corresponding to each waveform parameter in the first interval matrix.

[0107] In some embodiments, the first target wave is the R wave; the second target wave is the T wave and the P wave; and the first cycle is the first RR interval.

[0108] In some of these embodiments, Figure 5 On the basis of , it also includes a phase analysis module;

[0109] The phase analysis module is used to perform first target wave detection on the second ECG waveform data after synchronously and continuously acquiring the second ECG waveform data and the corresponding magnetic resonance data, and determine the second period between every two adjacent first target waves in the second ECG waveform data.

[0110] In some embodiments, the first processing module 230 is further configured to segment at least one second cycle in the second ECG waveform data according to the phase analysis parameters to obtain a phase interval matrix, and calculate the phase interval matrix to obtain a phase code corresponding to the second cycle.

[0111] In some of the embodiments, the second processing module 240 is further configured to determine a coding range according to a preset requirement;

[0112] According to the phase coding, the synchronously acquired magnetic resonance data within the coding range are coded and processed, and the magnetic resonance data with the same phase coding are image reconstructed to obtain target images of different phases.

[0113] In some embodiments, the second processing module 240 is further configured to assign corresponding phase codes to magnetic resonance data of different time periods based on a coding range, and to fill magnetic resonance data with the same phase code into corresponding positions in the K space for image reconstruction to obtain target images of different phases.

[0114] It should be noted that the above modules can be functional modules or program modules, and can be implemented through software or hardware. For modules implemented through hardware, the above modules can be located in the same processor; or the above modules can be located in different processors in any combination.

[0115] This embodiment further provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor is configured to run the computer program to execute the steps in any one of the above method embodiments.

[0116] Optionally, the computer device may further include a transmission device and an input / output device, wherein the transmission device is connected to the processor, and the input / output device is connected to the processor.

[0117] Optionally, in this embodiment, the processor may be configured to execute the following steps through a computer program:

[0118] S1, collecting first ECG waveform data, performing time phase analysis on the first ECG waveform data, and determining time phase analysis parameters;

[0119] S2, synchronously and continuously acquiring second ECG waveform data and corresponding magnetic resonance data;

[0120] S3, performing time phase segmentation on at least one second cycle in the second ECG waveform data according to the time phase analysis parameter, and determining a time phase code corresponding to each second cycle; wherein the second cycle is a period between two adjacent first target waves in the second ECG waveform data;

[0121] S4, encoding the synchronously acquired magnetic resonance data according to the phase encoding, and reconstructing the magnetic resonance data with the same phase encoding to obtain target images of different phases.

[0122] It should be noted that, for specific examples in this embodiment, reference may be made to the examples described in the above embodiments and optional implementation modes, and will not be repeated in this embodiment.

[0123] In addition, in conjunction with the MRI ECG-gated phase segmentation method provided in the above embodiments, a storage medium may also be provided in this embodiment to implement the method. The storage medium stores a computer program; when the computer program is executed by a processor, the computer program implements any of the MRI ECG-gated phase segmentation methods in the above embodiments.

[0124] It should be understood that the specific embodiments described herein are only used to explain this application and are not used to limit it. Based on the embodiments provided in this application, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of this application.

[0125] Obviously, the accompanying drawings are merely examples or embodiments of the present application. A person skilled in the art can also apply the present application to other similar situations based on these drawings without inventive effort. Furthermore, it is understandable that, although the work involved in this development process may be complex and lengthy, certain design, manufacturing, or production changes based on the technical content disclosed in this application are merely routine technical means for a person skilled in the art and should not be considered to constitute a deficiency in the disclosure of the present application.

[0126] The term "embodiment" as used in this application refers to specific features, structures, or characteristics described in conjunction with the embodiment that can be included in at least one embodiment of this application. The appearance of this phrase in various places in the specification does not necessarily mean that the embodiment is the same, nor does it mean that it is mutually exclusive with other embodiments and is independent or optional. It is understood, either explicitly or implicitly, by those skilled in the art that the embodiments described in this application can be combined with other embodiments when there is no conflict.

[0127] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of patent protection. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present application, and these modifications and improvements fall within the scope of protection of the present application. Therefore, the scope of protection of the present application shall be determined by the appended claims.

Claims

1. A magnetic resonance electrocardiogram (MRI) gating phase segmentation method, characterized in that: include: collecting first electrocardiogram waveform data, performing time phase analysis on the first electrocardiogram waveform data, and determining time phase analysis parameters; The phase analysis is used to adjust and obtain appropriate phase analysis parameters according to the different heart rates of each person; the phase analysis parameters are the parameters after the transformation of each waveform interval in the first interval matrix obtained by the phase analysis; synchronously and continuously acquiring second electrocardiogram waveform data and corresponding magnetic resonance data; Performing phase segmentation on at least one second cycle in the second ECG waveform data according to the phase analysis parameters, and determining a phase code corresponding to each second cycle; wherein the second cycle is a cycle between two adjacent first target waves in the second ECG waveform data; wherein performing phase segmentation on at least one second cycle in the second ECG waveform data according to the phase analysis parameters comprises: using the phase analysis parameters to segment at least one second cycle in the second ECG waveform data, thereby achieving accurate segmentation of cycles of different lengths and obtaining a phase interval matrix; The synchronously acquired magnetic resonance data are coded according to the phase coding, and the magnetic resonance data with the same phase coding are reconstructed to obtain target images of different phases.

2. The magnetic resonance electrocardiography gating phase segmentation method according to claim 1, characterized in that: The collecting of the first electrocardiogram waveform data and performing phase analysis on the first electrocardiogram waveform data includes: collecting first electrocardiogram waveform data, performing first target wave detection on the first electrocardiogram waveform data, and determining a first period between two adjacent first target waves in the first electrocardiogram waveform data; performing a second target wave detection on the electrocardiogram waveform data corresponding to the first cycle; The first cycle is divided based on the detection result of the second target wave to obtain a first interval matrix.

3. The magnetic resonance electrocardiography gating phase segmentation method according to claim 2, characterized in that: Determining the phase analysis parameters includes: Each waveform parameter in the first interval matrix is analyzed based on the first cycle to determine a phase analysis parameter corresponding to each waveform parameter in the first interval matrix.

4. The magnetic resonance electrocardiography gating phase segmentation method according to claim 1, characterized in that: The method further comprises: After synchronously and continuously acquiring the second ECG waveform data and the corresponding magnetic resonance data, first target wave detection is performed on the second ECG waveform data to determine a second period between every two adjacent first target waves in the second ECG waveform data.

5. The magnetic resonance electrocardiography gating phase segmentation method according to claim 4, characterized in that: The step of performing time phase segmentation on at least one second cycle in the second electrocardiogram waveform data according to the time phase analysis parameter and determining a time phase code corresponding to each second cycle further includes: The phase interval matrix is calculated to obtain a phase code corresponding to the second period.

6. The magnetic resonance electrocardiography gating phase segmentation method according to claim 1, characterized in that: The encoding process is performed on the synchronously acquired magnetic resonance data according to the phase encoding, and image reconstruction is performed on the magnetic resonance data having the same phase encoding to obtain target images of different phases, including: Determine the coding range based on preset requirements; According to the phase coding, the synchronously acquired magnetic resonance data within the coding range are coded, and the magnetic resonance data with the same phase coding are reconstructed to obtain target images of different phases.

7. The magnetic resonance electrocardiography gating phase segmentation method according to claim 6, characterized in that: The step of encoding the synchronously acquired magnetic resonance data within the encoding range according to the phase encoding, and reconstructing the magnetic resonance data with the same phase encoding to obtain target images of different phases includes: Based on the coding range, corresponding phase codes are assigned to the magnetic resonance data of different time periods, and the magnetic resonance data with the same phase code are filled into the corresponding position of the K space for image reconstruction to obtain target images of different phases.

8. A magnetic resonance electrocardiogram (MRI) gating phase segmentation device, characterized in that: include: a first acquisition module, a second acquisition module, a first processing module, and a second processing module; The first acquisition module is configured to acquire first ECG waveform data, perform time phase analysis on the first ECG waveform data to obtain a first interval matrix, and determine time phase analysis parameters based on the first interval matrix; the time phase analysis is configured to obtain appropriate time phase analysis parameters based on the different heart rates of each person; the time phase analysis parameters are parameters obtained by transforming each waveform interval in the first interval matrix obtained by the time phase analysis; The second acquisition module is used to synchronously and continuously acquire second ECG waveform data and corresponding magnetic resonance data; The first processing module is configured to perform phase segmentation on at least one second cycle in the second ECG waveform data according to the phase analysis parameters, and determine a phase code corresponding to each second cycle; wherein the second cycle is a cycle between two adjacent first target waves in the second ECG waveform data; wherein performing phase segmentation on at least one second cycle in the second ECG waveform data according to the phase analysis parameters comprises: using the phase analysis parameters to segment at least one second cycle in the second ECG waveform data, thereby achieving accurate segmentation of cycles of different lengths and obtaining a phase interval matrix; The second processing module is configured to perform encoding processing on the synchronously acquired magnetic resonance data according to the phase encoding, and to perform image reconstruction on the magnetic resonance data with the same phase encoding to obtain target images of different phases.

9. A computer device comprising a memory and a processor, characterized in that: A computer program is stored in the memory, and the processor is configured to run the computer program to perform the steps of the magnetic resonance electrocardiography gating phase segmentation method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the magnetic resonance electrocardiography gating phase segmentation method according to any one of claims 1 to 7 are implemented.

Images