Black-blood cine imaging method, apparatus, and computer device

Abstract

The application relates to a black-blood cine imaging method, device and computer equipment. The method comprises the following steps: acquiring K-space data and a physiological signal of a target part; a K-space where the K-space data is located comprises a center region and a peripheral region, a sampling rate of the center region is greater than that of the peripheral region; according to a state of the physiological signal, the K-space data is subjected to rearrangement processing to obtain rearranged K-space data; and image reconstruction is performed on the rearranged K-space data to obtain a black-blood cine image of the target part. The method can improve the speed of black-blood cine imaging.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of magnetic resonance imaging, in particular to a black-blood cine imaging method, device and computer equipment. BACKGROUND

[0002] The black-blood cine imaging method can suppress the blood flow signal of the imaging region, increase the contrast between the blood flow and the vessel wall, and has a wide range of applications in clinical practice.

[0003] In the prior art, three-dimensional gradient echo is usually used for black-blood cine imaging. However, three-dimensional gradient echo imaging is susceptible to the interference of factors such as magnetic susceptibility, and the imaging effect depends on the suppression effect on blood and fat.

[0004] Echo imaging is susceptible to factors such as magnetic susceptibility, and the imaging effect depends on the suppression effect on blood and fat. Spin echo imaging is less sensitive to magnetic field inhomogeneity, and can avoid the above problems and be suitable for imaging of blood vessels and organs in various parts of the body. However, spin echo imaging usually takes a long time to collect, and it is difficult to realize fast black-blood cine imaging.

[0005] Therefore, the current black-blood cine imaging technology has the problem of slow imaging speed. SUMMARY

[0006] Therefore, it is necessary to provide a black-blood cine imaging method, device, computer equipment and computer readable storage medium capable of fast imaging in view of the above technical problems.

[0007] In a first aspect, the present application provides a black-blood cine imaging method. The method comprises:

[0008] Obtaining K-space data and physiological signals of a target part; the K-space in which the K-space data is located comprises a center region and a peripheral region, and the sampling rate of the center region is greater than that of the peripheral region;

[0009] According to the state of the physiological signals, the K-space data is rearranged to obtain rearranged K-space data;

[0010] 5The rearranged K-space data is subjected to image reconstruction to obtain a black-blood cine image of the target part.

[0011] In one embodiment, the K-space data and physiological signals of the target part are obtained by:

[0012] According to a preset sampling trajectory, the target part is sampled in the K-space to obtain the K-space data.

[0013] In one embodiment, the K-space data is obtained by sampling the target part in the K-space according to a preset sampling trajectory, comprising:

[0014] According to the sampling trajectory, full sampling is performed on the target site in a central region of the K-space, and down sampling is performed on the target site in a peripheral region of the K-space, to obtain the K-space data.

[0015] In one of the embodiments, the down sampling of the target site in the peripheral region of the K-space comprises:

[0016] In the peripheral region of the K-space, the target site is subjected to random down sampling; the sampling rate of the random down sampling is less than the Nyquist sampling rate.

[0017] In one of the embodiments, the obtaining of the K-space data and the physiological signal of the target site further comprises:

[0018] determining a sampling time of the K-space data of the target site;

[0019] According to the sampling time, a physiological signal of the target site is collected.

[0020] In one of the embodiments, the rearrangement of the K-space data according to the state of the physiological signal to obtain the rearranged K-space data comprises:

[0021] According to the state of the physiological signal, a state stage corresponding to the physiological signal is determined;

[0022] According to the state stage, the K-space data is rearranged to obtain the rearranged K-space data.

[0023] In one of the embodiments, the physiological signal comprises an electrocardiogram signal, a respiration signal and / or a motion signal, the state stage comprises a heart stage corresponding to the electrocardiogram signal, a respiration stage corresponding to the respiration signal and / or a motion stage corresponding to the motion signal; the rearrangement of the K-space data according to the state of the physiological signal to obtain the rearranged K-space data further comprises:

[0024] determining the heart stage, the respiration stage and / or the motion stage;

[0025] According to the heart stage, the respiration stage and / or the motion stage, the K-space data is rearranged to obtain the rearranged K-space data.

[0026] In one of the embodiments, the image reconstruction of the rearranged K-space data to obtain the black-blood cine image of the target site comprises:

[0027] performing Fourier transform on the rearranged K-space data to obtain transformed K-space data;

[0028] performing image reconstruction on the transformed K-space data to obtain the black-blood cine image of the target part.

[0029] In a second aspect, the present application further provides a black-blood cine imaging device. The device comprises:

[0030] an acquisition module configured to acquire K-space data of a target part and a physiological signal; the K-space data is located in a K-space, which comprises a central region and a peripheral region, and a sampling rate of the central region is greater than that of the peripheral region;

[0031] an arrangement module configured to perform rearrangement processing on the K-space data according to a state of the physiological signal to obtain rearranged K-space data;

[0032] a reconstruction module configured to perform image reconstruction processing on the rearranged K-space data to obtain a black-blood cine image of the target part.

[0033] In a third aspect, the present application further provides a computer device. The computer device comprises a memory and a processor, the memory stores a computer program, and the processor implements the following steps when executing the computer program:

[0034] acquiring K-space data of a target part and a physiological signal; the K-space data is located in a K-space, which comprises a central region and a peripheral region, and a sampling rate of the central region is greater than that of the peripheral region;

[0035] performing rearrangement processing on the K-space data according to a state of the physiological signal to obtain rearranged K-space data;

[0036] performing image reconstruction on the rearranged K-space data to obtain a black-blood cine image of the target part.

[0037] In a fourth aspect, the present application further provides a computer readable storage medium. The computer readable storage medium stores a computer program, and the computer program is executed by a processor to implement the following steps:

[0038] acquiring K-space data of a target part and a physiological signal; the K-space data is located in a K-space, which comprises a central region and a peripheral region, and a sampling rate of the central region is greater than that of the peripheral region;

[0039] performing rearrangement processing on the K-space data according to a state of the physiological signal to obtain rearranged K-space data;

[0040] reconstructing the rearranged K-space data to obtain the black-blood cine image of the target part.

[0041] The black-blood cine imaging method, device, computer device and storage medium can perform full sampling in the center region of the K-space and down-sampling in the peripheral region, retain low-frequency information in the center region and reduce high-frequency information in the peripheral region, reduce the number of K-space sampling on the basis of ensuring the imaging quality, and improve the speed of black-blood cine imaging.

[0042] Moreover, the K-space data is rearranged according to the state of the physiological signal, and the image is reconstructed based on the rearranged K-space data, so that the image can be reconstructed according to the physiological characteristics of the target part, and the reconstructed cine image can reflect the tissue structure, motion characteristics or strain characteristics of the target part under different physiological conditions. BRIEF DESCRIPTION OF DRAWINGS

[0043] Figure 1 A flowchart of a black-blood cine imaging method in an embodiment is shown in FIG. 1.

[0044] Figure 2 A schematic diagram of K-space in which non-uniform down-sampling is performed according to a non-Cartesian trajectory in an embodiment is shown in FIG. 3.

[0045] Figure 3 A schematic diagram of K-space in which non-uniform down-sampling is performed according to a non-Cartesian trajectory in an embodiment is shown in FIG. 3.

[0046] Figure 4 A schematic diagram of monitoring a three-dimensional spin echo sequence according to a physiological signal in an embodiment is shown in FIG. 5.

[0047] Figure 5 A schematic diagram of rearranging K-space data according to a physiological signal in an embodiment is shown in FIG. 6.

[0048] Figure 6 A flowchart of a magnetic resonance black-blood cine imaging method in an embodiment is shown in FIG. 7.

[0049] Figure 7 A block diagram of a black-blood cine imaging device in an embodiment is shown in FIG. 8.

[0050] Figure 8 An internal structure diagram of a computer device in an embodiment is shown in FIG. 9. DETAILED DESCRIPTION

[0051] In order to make the purposes, technical solutions and advantages of the present application clearer, the present application will be further described in detail below with reference to the drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and not intended to limit the present application.

[0052] In one embodiment, as shown in Figure 1 A black-blood cine imaging method is provided, and the embodiment is exemplified by the method applied to a terminal. It should be understood that the method can also be applied to a server, and can also be applied to a system including a terminal and a server, and can be realized through the interaction of the terminal and the server. In the embodiment, the method includes the following steps.

[0053] In step S110, K-space data and physiological signals of a target site are acquired. The K-space in which the K-space data is located includes a central region and a peripheral region, and the sampling rate of the central region is greater than that of the peripheral region.

[0054] The target site can be a site for magnetic resonance imaging. For example, a blood vessel site.

[0055] The K-space data can be data stored in the K-space, which is used to describe the spatial frequency information of a magnetic resonance image.

[0056] The physiological signals can be, but are not limited to, electrocardiogram signals, respiratory signals or motion signals.

[0057] In a specific implementation, the magnetic resonance device can be controlled to use a spin echo sequence to scan the target site to obtain the K-space data of the target site, and the K-space data of the target site is sent to the terminal. The physiological signal acquisition device can also be controlled to acquire the physiological signals of the target site, and the physiological signals of the target site are also sent to the terminal, so that the terminal can acquire the K-space data and the physiological signals of the target site.

[0058] The K-space data and the physiological signals can be acquired synchronously.

[0059] The K-space in which the K-space data is located can be divided into a central region close to the center and a peripheral region away from the center. The sampling rate of the K-space data in the central region can be full sampling, and the sampling rate in the peripheral region can be reduced sampling.

[0060] In practical applications, the magnetic resonance device can use a three-dimensional fast spin echo sequence to scan a blood vessel part to obtain K-space data, wherein the K-space data acquisition trajectory can be but is not limited to a Cartesian trajectory, a radial trajectory, a spiral trajectory, a golden angle radial acquisition trajectory, or a spiral stack trajectory, and the K-space data is fully sampled in the K-space center region and is randomly down-sampled with an appropriate down-sampling multiple in the K-space peripheral region. For different blood vessel parts, different physiological signal acquisition devices can also be used to synchronously acquire physiological signals. Typically, for heart blood vessels, an electrocardiogram monitoring device can be used to synchronously acquire an electrocardiogram signal; for lung blood vessels, a respiration monitoring device can be used to synchronously acquire a respiration signal; and for neck blood vessels, a motion monitoring device can be used to synchronously acquire a swallowing motion signal.

[0061] Figure 2 A schematic diagram of K-space with non-uniform down-sampling according to a Cartesian trajectory is provided. According to Figure 2 , the K-space can include a center region 201 and a peripheral region 202, and using a Cartesian sampling trajectory, the center region 201 is fully sampled and the peripheral region 202 is down-sampled. Due to the different sampling rates of the center region 201 and the peripheral region 202, the sampling rate of the entire K-space is non-uniform. By reasonably setting the sampling rate of the peripheral region 202, the average sampling rate of the entire K-space can be much smaller than the Nyquist sampling rate. Since the center region of the K-space is low-frequency information and the peripheral region is high-frequency information, down-sampling in the peripheral region can reduce the sampling time while ensuring the overall contrast.

[0062] Figure 3 A schematic diagram of K-space with non-uniform down-sampling according to a non-Cartesian trajectory is provided. According to Figure 3 , the sampling trajectory can also be a radial trajectory, a spiral trajectory, a golden angle radial acquisition trajectory, or a spiral stack trajectory, etc.

[0063] Figure 4 A schematic diagram of monitoring a three-dimensional spin echo sequence according to a physiological signal is provided. According to Figure 4 , while using a three-dimensional fast spin echo sequence to scan, an electrocardiogram monitoring device can be used to synchronously acquire an electrocardiogram signal, a respiration monitoring device can be used to synchronously acquire a respiration signal, or a motion monitoring device can be used to synchronously acquire a motion signal.

[0064] In step S120, the K-space data is rearranged according to the state of the physiological signal to obtain rearranged K-space data.

[0065] The state of the physiological signal can be, but is not limited to, a state of an electrocardiogram signal, a respiratory signal, or a motion signal. For example, for an electrocardiogram signal, the state can be a state of different cardiac cycles; for a respiratory signal, the state can be a state of an inhalation phase, an exhalation phase, and a plateau phase; and for a swallowing motion signal, the state can be a normal or abnormal swallowing state.

[0066] In specific implementations, the terminal can determine the state of the physiological signal based on the acquired physiological signal, and rearrange the K-space data based on the state of the physiological signal to obtain rearranged K-space data.

[0067] In actual applications, the terminal can divide the acquired electrocardiogram signal into multiple cardiac phases (cardiac phases) according to cardiac cycles, and can also divide the acquired respiratory signal into an inhalation phase, an exhalation phase, and a plateau phase, and the like according to a respiratory wave to obtain respiratory phases. For the acquired motion signal, the terminal can also identify a normal or abnormal motion state from the motion signal, and determine a motion phase based on the identified motion state. In the time dimension, the terminal can re-allocate the K-space data, and rearrange the K-space data according to the cardiac phase, the respiratory phase, or the motion phase at the time of acquisition to obtain rearranged K-space data.

[0068] Figure 5 An example of rearranging K-space data based on a physiological signal is provided. According to the example, Figure 5 For a heart blood vessel, the K-space data can be rearranged according to the cardiac phase at the time of acquisition; for a lung blood vessel, the K-space data can be rearranged according to the respiratory phase at the time of acquisition; and for a neck blood vessel, the K-space data can be rearranged according to the motion phase and the cardiac phase at the time of acquisition.

[0069] At step S130, the rearranged K-space data is subjected to image reconstruction to obtain a black-blood cine image of the target site.

[0070] The black-blood cine image can be a dynamic series of images composed of multiple consecutive magnetic resonance images.

[0071] In specific implementations, the rearranged K-space data can be subjected to Fourier transform to obtain transformed K-space data, and the transformed K-space data can be subjected to image reconstruction to obtain a cine image of the target site.

[0072] The image reconstruction method can be, but is not limited to, a deep learning-based, machine learning-based, compressed sensing-based, parallel imaging-based, or low-rank matrix-based image reconstruction method.

[0073] For example, for K-space data rearranged according to cardiac staging, a black-blood cine image of a cardiac blood vessel can be obtained through image reconstruction; for K-space data rearranged according to respiratory staging, a black-blood cine image of a pulmonary blood vessel can be obtained through image reconstruction; for K-space data rearranged according to neck movement staging and cardiac staging, a black-blood cine image of a neck blood vessel can be obtained through image reconstruction.

[0074] The above black-blood cine imaging method can obtain K-space data of a target part and a physiological signal, wherein a K-space in which the K-space data is located includes a central region and a peripheral region, a sampling rate of the central region is greater than a sampling rate of the peripheral region, the K-space data is rearranged according to a state of the physiological signal to obtain rearranged K-space data, and a black-blood cine image of the target part is obtained through image reconstruction of the rearranged K-space data. The central region of the K-space can be fully sampled, and the peripheral region can be down-sampled. The low-frequency information of the central region is retained, and the high-frequency information of the peripheral region is reduced. The number of K-space samplings is reduced on the basis of ensuring imaging quality, and the speed of black-blood cine imaging is improved.

[0075] Moreover, the K-space data is rearranged according to the state of the physiological signal, and image reconstruction is performed based on the rearranged K-space data. The image reconstruction can be performed according to the physiological characteristics of the target part, so that the reconstructed cine image can reflect the tissue structure, movement characteristics or strain characteristics of the target part under different physiological conditions.

[0076] In one embodiment, the above step S110 can specifically include: sampling a target part in a K-space according to a preset sampling trajectory to obtain K-space data.

[0077] In a specific implementation, the terminal can control the magnetic resonance device to use a spin echo sequence, sample a target part according to a preset sampling trajectory of a K-space, and obtain K-space data. The sampling trajectory can be, but is not limited to, any one of a Cartesian trajectory, a radial trajectory, a spiral line trajectory, a golden angle radial acquisition trajectory or a spiral line stack trajectory. The sampling trajectory can also be a non-uniform down-sampling sampling trajectory that fully samples in a central region and down-samples in a peripheral region. By reasonably setting the sampling rate of the down-sampling in the peripheral region, the average sampling rate of the K-space data can be much smaller than the Nyquist sampling rate.

[0078] In this embodiment, the K-space data is obtained by sampling a target part in a K-space according to a preset sampling trajectory. The image reconstruction can be performed using less K-space data on the basis of ensuring imaging quality, and the speed of black-blood cine imaging is improved.

[0079] In one embodiment, the step of sampling the target region in the K-space according to the preset sampling trajectory to obtain the K-space data can specifically include: performing full sampling on the target region in a center region of the K-space and performing down-sampling on the target region in a peripheral region of the K-space according to the sampling trajectory to obtain the K-space data.

[0080] In a specific implementation, the terminal can control the magnetic resonance device to use a spin echo sequence, perform full sampling on the target region in a center region of the K-space and perform down-sampling on the target region in a peripheral region of the K-space according to a preset sampling trajectory to obtain K-space data. The sampling rate of the down-sampling in the peripheral region can be set to be far less than the Nyquist sampling rate.

[0081] In this embodiment, by performing full sampling on the target region in a center region of the K-space and performing down-sampling on the target region in a peripheral region of the K-space according to the sampling trajectory to obtain the K-space data, the image reconstruction can be performed using less K-space data on the basis of ensuring the imaging quality, and the speed of the black-blood cine imaging is improved.

[0082] In one embodiment, the step of performing down-sampling on the target region in the peripheral region of the K-space can specifically include: performing random down-sampling on the target region in the peripheral region of the K-space; and the sampling rate of the random down-sampling is less than the Nyquist sampling rate.

[0083] In a specific implementation, the random down-sampling can be performed on the target region in the peripheral region of the K-space, and the sampling rate of the down-sampling in the peripheral region can be set to be less than the Nyquist sampling rate. Further, the sampling rate of the down-sampling in the peripheral region can be set to be far less than the Nyquist sampling rate.

[0084] In this embodiment, by performing random down-sampling on the target region in the peripheral region of the K-space, the image reconstruction can be performed using less K-space data on the basis of ensuring the imaging quality, and the speed of the black-blood cine imaging is improved.

[0085] In one embodiment, the step S110 can specifically further include: determining a sampling time of the K-space data of the target region; and collecting a physiological signal of the target region according to the sampling time.

[0086] The physiological signal includes at least one of an electrocardiogram signal, a respiration signal, and a motion signal.

[0087] In a specific implementation, the terminal can obtain the sampling time of the K-space data of the target region, and control a physiological signal collection device to collect the physiological signal of the target region according to the obtained sampling time.

[0088] In actual application, the terminal can send an instruction of collecting the physiological signal of the target part to the physiological signal collection device at the sampling time of the K-space data, and the physiological signal collection device collects the physiological signal of the target part upon receiving the instruction. The sampling time can also be preset, and the magnetic resonance device and the physiological signal collection device are controlled to synchronously collect the K-space data and the physiological signal of the target part at the preset sampling time, and the collected K-space data and physiological signal are respectively sent to the terminal.

[0089] In this embodiment, the K-space data and the physiological signal of the target part are synchronously collected by determining the sampling time of the K-space data of the target part and collecting the physiological signal of the target part according to the sampling time, which facilitates monitoring the black-blood cine image according to the physiological signal of the target part, and further comprehensively analyzing the tissue structure, motion characteristics or strain characteristics of the target part.

[0090] In one embodiment, the above step S120 can specifically include: determining a state staging corresponding to the physiological signal according to the state of the physiological signal; and rearranging the K-space data according to the state staging to obtain the rearranged K-space data.

[0091] The state staging can be a staging corresponding to the state of the physiological signal, which can be but is not limited to a cardiac phase, a respiratory phase or a motion phase.

[0092] In specific implementation, the terminal can determine the state staging of the physiological signal according to the state of the physiological signal, and can rearrange the K-space data according to the state staging after obtaining the K-space data to obtain the rearranged K-space data.

[0093] In actual application, for the electrocardiogram signal, the electrocardiogram signal can be staged according to a cardiac cycle to obtain a cardiac phase, and the K-space data can be rearranged according to the cardiac phase; for the respiratory signal, the respiratory signal can be staged according to an inspiration phase, an expiration phase and a plateau period to obtain a respiratory phase, and the K-space data can be rearranged according to the respiratory phase; for the motion signal, the normal or abnormal motion state can be identified, the motion signal can be staged according to the normal or abnormal motion state to obtain a motion phase, and the K-space data can be rearranged according to the cardiac phase and the motion phase.

[0094] In the embodiment, the state period corresponding to the physiological signal is determined according to the state of the physiological signal, the K-space data is rearranged according to the state period, and the rearranged K-space data is obtained. The K-space data can be rearranged according to the state of the physiological signal of the target part, so that the reconstructed cine image can reflect the tissue structure, motion characteristics or strain characteristics of the target part under different physiological signal states, and the target part can be comprehensively analyzed.

[0095] In one embodiment, the physiological signal includes an electrocardiogram signal, a respiratory signal and / or a motion signal, and the state period includes a heart period corresponding to the electrocardiogram signal, a respiratory period corresponding to the respiratory signal and / or a motion period corresponding to the motion signal. The step S120 can further include determining the heart period, the respiratory period and / or the motion period, and rearranging the K-space data according to the heart period, the respiratory period and / or the motion period to obtain the rearranged K-space data.

[0096] In a specific implementation, the terminal can divide the collected electrocardiogram signal into multiple heart periods (heart phases) according to the heart cycle, divide the collected respiratory signal into respiratory periods such as the inspiration phase, the expiration phase and the plateau period according to the respiratory wave, and / or identify normal or abnormal motion states from the collected motion signal, determine the motion period according to the normal or abnormal motion state, and rearrange the K-space data according to the determined heart period, respiratory period and / or motion period to obtain the rearranged K-space data.

[0097] In the embodiment, the heart period, the respiratory period and / or the motion period are determined, and the K-space data is rearranged according to the heart period, the respiratory period and / or the motion period to obtain the rearranged K-space data. The K-space data can be rearranged according to the state of the physiological signal of the target part, so that the reconstructed cine image can reflect the tissue structure, motion characteristics or strain characteristics of the target part under different physiological signal states, and the target part can be comprehensively analyzed.

[0098] In one embodiment, the step S130 can specifically include performing Fourier transform on the rearranged K-space data to obtain transformed K-space data, and performing image reconstruction on the transformed K-space data to obtain the black-blood cine image of the target part.

[0099] The image reconstruction includes but is not limited to image reconstruction based on deep learning, machine learning, compressed sensing, parallel imaging or low-rank matrix.

[0100] In a specific implementation, for the rearranged K-space data, the terminal can perform Fourier transform on the K-space data in the time dimension to convert the K-space data to sparse space to obtain transformed K-space data, and perform image reconstruction on the transformed K-space data to obtain a black-blood cine image of a target site for different physiological signals. The image reconstruction method can be, but is not limited to, a deep learning-based, machine learning-based, compressed sensing-based, parallel imaging-based, or low-rank matrix-based image reconstruction method.

[0101] In actual applications, a dynamic image down-sampling reconstruction method based on space-time information can be used to perform Fourier transform on the rearranged K-space data in the time dimension to convert the K-space data to sparse space, and use a deep learning-based, machine learning-based, compressed sensing-based, parallel imaging-based, or low-rank matrix-based method to reconstruct the K-space data to realize black-blood cine imaging under different physiological signal dimensions.

[0102] The dynamic image down-sampling reconstruction method can use a low-rank matrix-based reconstruction method to arrange the data into a matrix and use the low-rank characteristics of the matrix as a new constraint condition for reconstruction.

[0103] The deep learning-based reconstruction algorithm can use the similarity of adjacent time frame images to use a neural network to de-alias the reconstructed image to realize high-quality reconstruction of the cine image.

[0104] In this embodiment, Fourier transform is performed on the rearranged K-space data to obtain transformed K-space data, and image reconstruction is performed on the transformed K-space data to obtain a black-blood cine image of a target site. The rearranged K-space data can be reconstructed to realize black-blood cine imaging for different physiological signals, so that the reconstructed black-blood cine image can reflect the tissue structure, motion characteristics, or strain characteristics of the target site under different physiological signal states, facilitating comprehensive analysis of the target site.

[0105] To facilitate a person skilled in the art to thoroughly understand the embodiments of the present application, the following will be described in conjunction with a specific example.

[0106] Figure 6 A flowchart of a magnetic resonance black-blood cine imaging method is provided. According to Figure 6 , the method comprises the following steps:

[0107] In step S310, a three-dimensional fast spin echo sequence is used to perform non-uniform down-sampling to quickly image a target region.

[0108] Specifically, a three-dimensional fast spin echo sequence can be used, a preparation pulse is applied, black-blood multi-contrast imaging is achieved, the fast spin echo sequence uses an excitation pulse and a refocusing pulse to realize one-time excitation and acquisition of multiple echoes, and multiple data lines of K-space are filled, which can be suitable for any acquisition trajectory, including 3D Cartesian trajectory, 3D radial trajectory, 3D spiral trajectory, koosh ball (golden angle radial acquisition trajectory), stack of spiral (spiral stack trajectory), etc.; combined with non-uniform down-sampling, the sparse characteristics or low-rank characteristics of the signal are used to implement full sampling in the central region of K-space and random down-sampling in the peripheral region of K-space using an appropriate down-sampling multiple, so as to realize fast imaging.

[0109] In step S320, physiological information such as electrocardiogram, respiration, and motion is synchronously acquired.

[0110] Specifically, electrocardiogram, respiration, and motion physiological signals can be synchronously acquired by using electrocardiogram gating, respiration gating, motion monitoring, etc.; the acquired electrocardiogram signal can be divided into multiple cardiac sub-stages (cardiac phases) according to the cardiac cycle, the acquired respiration signal can be divided into inspiration phase, expiration phase, and plateau phase according to the respiration wave, and in addition, abnormal motion signals can be extracted from the acquired motion signal, and the motion sub-stage can be determined according to the abnormal motion signals.

[0111] In step S330, K-space sampling points are re-distributed according to the cardiac sub-stage, respiration sub-stage, and motion sub-stage in which the sampling points are located when the sampling points are acquired.

[0112] Specifically, the K-space sampling points can be re-distributed in the time dimension, and the sampling points can be rearranged into a new K-space according to the cardiac sub-stage, respiration sub-stage, and motion sub-stage in which the sampling points are located when the sampling points are acquired.

[0113] In step S340, a dynamic image down-sampling reconstruction method based on space-time information is used to realize black-blood movie imaging of magnetic resonance under different motion dimensions.

[0114] Specifically, the dynamic image down-sampling reconstruction method based on space-time information can be used to perform Fourier transform in the time dimension, convert the data of the new K-space to a sparse space, and use deep learning, machine learning, compressed sensing, parallel imaging, low-rank matrix, etc. to reconstruct, so as to realize movie imaging of black-blood images under different motion dimensions.

[0115] The reconstruction method based on low-rank matrix can arrange the data of the new K-space into a matrix, and use the low-rank characteristics as a new constraint condition for reconstruction.

[0116] The deep learning-based reconstruction algorithm uses the similarity of adjacent time frame images, uses a neural network to de-alias the reconstructed images, and thus realizes high-quality reconstruction of cine imaging.

[0117] The above-mentioned magnetic resonance black-blood cine imaging method realizes high spatial resolution imaging by using non-uniform down-sampling, realizes high temporal resolution imaging by using physiological signal monitoring, and realizes three-dimensional black-blood cine imaging by combining with the fast spin echo technique, so as to comprehensively analyze the blood vessel shape, the tube wall feature, the plaque strain or the organ structure, and the blood flow strain.

[0118] In the non-uniform down-sampling, information far less than the Nyquist sampling rate is used to obtain high image resolution image reconstruction through post-processing reconstruction under high down-sampling multiple.

[0119] Moreover, multi-directional physiological feature acquisition can be realized by combining with electrocardiogram gating and motion monitoring, and information such as electrocardiogram, respiration and motion can be fully acquired, so that discrete signals under different features are fully extracted in the transform domain for reconstruction, the motion state is decomposed, and high temporal resolution imaging is realized.

[0120] Further, black-blood cine imaging can be realized by combining with the fast spin echo technique, and the advantages such as high robustness and high signal-to-noise ratio can be used to realize comprehensive analysis of the whole body blood vessels and various organs.

[0121] It should be understood that, although each step in the flowchart involved in each embodiment as described above is displayed in sequence according to the arrow, these steps are not necessarily executed in sequence according to the arrow. Unless otherwise specified herein, the execution of these steps is not strictly limited in sequence, and these steps can be executed in other sequences. Moreover, at least part of the steps in the flowchart involved in each embodiment as described above can include multiple steps or multiple stages, which are not necessarily executed at the same time, but can be executed at different times, and the execution sequence of these steps or stages is not necessarily sequential, but can be executed in rotation or alternation with at least part of other steps or steps or stages in other steps.

[0122] Based on the same inventive concept, the embodiments of the present application also provide a black-blood cine imaging device for realizing the above-mentioned black-blood cine imaging method. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme described in the above-mentioned method, and therefore the specific limitations in one or more black-blood cine imaging device embodiments provided below can refer to the limitations of the black-blood cine imaging method described above, which will not be repeated here.

[0123] In one embodiment, as Figure 7As shown, a black-blood cine imaging apparatus is provided, comprising: an acquisition module 410, an arrangement module 420 and a reconstruction module 430, wherein:

[0124] The acquisition module 410 is configured to acquire K-space data and a physiological signal of a target site; a K-space where the K-space data is located comprises a center region and a peripheral region, and a sampling rate of the center region is greater than a sampling rate of the peripheral region.

[0125] The arrangement module 420 is configured to perform rearrangement processing on the K-space data according to a state of the physiological signal, to obtain rearranged K-space data.

[0126] The reconstruction module 430 is configured to perform image reconstruction processing on the rearranged K-space data, to obtain a black-blood cine image of the target site.

[0127] In one embodiment, the acquisition module 410 is further configured to perform sampling processing on the target site in the K-space according to a preset sampling trajectory, to obtain the K-space data.

[0128] In one embodiment, the acquisition module 410 is further configured to perform full sampling processing on the target site in the center region of the K-space and down-sampling processing on the target site in the peripheral region of the K-space according to the sampling trajectory, to obtain the K-space data.

[0129] In one embodiment, the acquisition module 410 is further configured to perform random down-sampling processing on the target site in the peripheral region of the K-space; a sampling rate of the random down-sampling processing is less than a Nyquist sampling rate.

[0130] In one embodiment, the acquisition module 410 is further configured to determine a sampling time of the K-space data of the target site; and collect the physiological signal of the target site according to the sampling time.

[0131] In one embodiment, the arrangement module 420 is further configured to determine a state staging corresponding to the physiological signal according to the state of the physiological signal; and perform rearrangement processing on the K-space data according to the state staging, to obtain the rearranged K-space data.

[0132] In one embodiment, the arrangement module 420 is further configured to determine the cardiac staging, the respiratory staging and / or the motion staging; and perform rearrangement processing on the K-space data according to the cardiac staging, the respiratory staging and / or the motion staging, to obtain the rearranged K-space data.

[0133] In one embodiment, the reconstruction module 430 is further configured to perform a Fourier transform on the rearranged K-space data to obtain transformed K-space data; and to perform image reconstruction on the transformed K-space data to obtain the black-blooded movie image of the target area.

[0134] Each module in the aforementioned black-blooded film imaging device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the memory of a computer device as software, so that the processor can call and execute the operations corresponding to each module.

[0135] In one embodiment, a computer device is provided, which may be a terminal, and its internal structure diagram may be as follows: Figure 8 As shown, the computer device includes a processor, memory, input / output interface, communication interface, display unit, and input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are also connected to the system bus via the input / output interface. The processor provides computing and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system and computer programs. The internal memory provides an environment for the operation of the operating system and computer programs stored in the non-volatile storage media. The input / output interface is used for exchanging information between the processor and external devices. The communication interface is used for wired or wireless communication with external terminals; wireless communication can be achieved through Wi-Fi, mobile cellular networks, NFC (Near Field Communication), or other technologies. When the computer program is executed by the processor, it implements a black-blooded film imaging method. The display unit of the computer device is used to form a visually visible image and can be a display screen, a projection device, or a virtual reality imaging device. The display screen can be an LCD screen or an e-ink screen. The input device of the computer device can be a touch layer covering the display screen, or buttons, trackballs, or touchpads set on the casing of the computer device, or external keyboards, touchpads, or mice, etc.

[0136] Those skilled in the art will understand that Figure 8 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0137] In an embodiment, a computer device is also provided, comprising a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps in the above-mentioned method embodiments when executing the computer program.

[0138] In an embodiment, a computer readable storage medium is provided, which stores a computer program, and the computer program is executed by a processor to implement the steps in the above-mentioned method embodiments.

[0139] It should be noted that the user information (including but not limited to user equipment information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or authorized by all parties, and the collection, use and processing of related data need to comply with relevant laws, regulations and standards of relevant countries and regions.

[0140] Those skilled in the art can understand that all or part of the processes in the above-mentioned embodiment methods can be completed by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer readable storage medium, and when the computer program is executed, the processes of the above-mentioned embodiments of the methods can be included. Any reference to memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (Read-Only Memory, ROM), magnetic tape, floppy disk, flash memory, optical storage, high-density embedded non-volatile memory, resistive memory (ReRAM), magnetoresistive random access memory (Magnetoresistive Random Access Memory, MRAM), ferroelectric memory (Ferroelectric Random Access Memory, FRAM), phase change memory (Phase Change Memory, PCM), graphene memory, etc. Volatile memory can include random access memory (Random Access Memory, RAM) or external cache memory, etc. As an illustration but not limitation, RAM can be in various forms, such as static random access memory (Static Random Access Memory, SRAM) or dynamic random access memory (Dynamic Random Access Memory, DRAM), etc. The database involved in the embodiments provided in the present application can include at least one of a relational database and a non-relational database. The non-relational database can include a distributed database based on a block chain, etc., without being limited thereto. The processor involved in the embodiments provided in the present application can be a general-purpose processor, a central processing unit, a graphics processing unit, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, etc., without being limited thereto.

[0141] Any combination of the technical features of the above embodiments can be made. In order to make the description simple, all possible combinations of the technical features in the above embodiments are not described, however, as long as the combination of the technical features does not exist contradictory, it should be considered as the scope of the present application.

[0142] The above embodiments only express several implementation manners of the present application, and the description is more specific and detailed, but it should not be understood as a limitation on the scope of the patent of the present application. It should be pointed out that for ordinary skilled in the art, without departing from the concept of the present application, a number of modifications and improvements can be made, which are within the scope of protection of the present application. Therefore, the protection scope of the present application should be subject to the appended claims.

Claims

1. A black-blood cine-imaging method, characterized by, The method comprises: acquiring K-space data and physiological signals of a target site; the K-space in which the K-space data is located comprises a central region and a peripheral region, a sampling rate of the central region is greater than a sampling rate of the peripheral region, and the physiological signals comprise electrocardiogram signals, respiration signals, and / or motion signals; according to a state of the physiological signals, dividing the electrocardiogram signals into a plurality of cardiac phases according to a cardiac cycle, dividing the respiration signals into respiration phases of an inhalation phase, an exhalation phase, and a plateau period according to a respiration wave, and / or determining a motion phase according to a normal or abnormal motion state in the motion signals; performing rearrangement processing on the K-space data according to the cardiac phases, the respiration phases, and / or the motion phase, to obtain rearranged K-space data; performing image reconstruction on the rearranged K-space data, to obtain a black-blood cine image of the target site.

2. The method of claim 1, wherein, The acquiring of the K-space data and the physiological signals of the target site comprises: performing sampling processing on the target site in the K-space according to a preset sampling trajectory, to obtain the K-space data.

3. The method of claim 2, wherein, The performing of the sampling processing on the target site in the K-space according to the preset sampling trajectory, to obtain the K-space data, comprises: performing full sampling processing on the target site in the central region of the K-space according to the sampling trajectory, and performing down-sampling processing on the target site in the peripheral region of the K-space, to obtain the K-space data.

4. The method of claim 3, wherein, The performing of the down-sampling processing on the target site in the peripheral region of the K-space comprises: performing random down-sampling processing on the target site in the peripheral region of the K-space; a sampling rate of the random down-sampling processing is less than a Nyquist sampling rate.

5. The method of claim 1, wherein, The acquiring of the K-space data and the physiological signals of the target site further comprises: determining a sampling time of the K-space data of the target site; collecting the physiological signals of the target site according to the sampling time.

6. The method of claim 1, wherein, The performing of the image reconstruction on the rearranged K-space data, to obtain the black-blood cine image of the target site, comprises: performing Fourier transform on the rearranged K-space data, to obtain transformed K-space data; performing image reconstruction on the transformed K-space data, to obtain the black-blood cine image of the target site.

7. A black-blood cine-imaging apparatus, characterized by comprising: The device comprises: an acquisition module, configured to acquire K-space data and physiological signals of a target site; the K-space in which the K-space data is located comprises a central region and a peripheral region, a sampling rate of the central region is greater than a sampling rate of the peripheral region, and the physiological signals comprise electrocardiogram signals, respiration signals, and / or motion signals; a rearrangement module, configured to divide the electrocardiogram signals into a plurality of cardiac phases according to a cardiac cycle according to a state of the physiological signals, divide the respiration signals into respiration phases of an inhalation phase, an exhalation phase, and a plateau period according to a respiration wave, and / or determine a motion phase according to a normal or abnormal motion state in the motion signals, perform rearrangement processing on the K-space data according to the cardiac phases, the respiration phases, and / or the motion phase, to obtain rearranged K-space data; A reconstruction module is configured to perform image reconstruction processing on the rearranged K-space data to obtain a black-blood cine image of the target part.

8. The apparatus of claim 7, wherein, The acquisition module is further configured to sample the target part according to a preset sampling trajectory to obtain the K-space data. 9.A computer device, comprising a memory and a processor, wherein the memory stores a computer program, and the computer device is configured to perform the method according to any one of claims 1-8 when the computer program is executed by the processor. The processor, when executing the computer program, implements the steps of the method in any one of claims 1 to 6.

10. A computer-readable storage medium having stored thereon a computer program, characterized in that, The computer program, when executed by the processor, implements the steps of the method in any one of claims 1 to 6.

Images