CN110276725A - A fast method for removing artifacts in free-breathing lower abdomen MRI - Google Patents
A fast method for removing artifacts in free-breathing lower abdomen MRI Download PDFInfo
- Publication number
- CN110276725A CN110276725A CN201810212643.0A CN201810212643A CN110276725A CN 110276725 A CN110276725 A CN 110276725A CN 201810212643 A CN201810212643 A CN 201810212643A CN 110276725 A CN110276725 A CN 110276725A
- Authority
- CN
- China
- Prior art keywords
- artifact
- magnetic resonance
- image
- resonance imaging
- lower abdomen
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及一种自由呼吸下腹部磁共振成像伪影的快速去除方法,包含采集带伪影的图像数据部分、预处理得到若干子图部分、通过校正模型消除伪影部分以及L个子图拼接得到输出图像部分。本发明能够实时去除由于高倍亚采样图像中带有的条纹伪影,具有临床潜在应用价值。
The present invention relates to a method for rapidly removing artifacts in MRI of the lower abdomen under free breathing, which includes collecting image data with artifacts, preprocessing to obtain several subimages, eliminating artifacts through a correction model, and splicing L subimages to obtain Output image parts. The invention can remove the streak artifact in the high-magnification sub-sampling image in real time, and has clinical potential application value.
Description
技术领域technical field
本发明涉及一种自由呼吸下腹部磁共振成像伪影的快速去除方法,属于图像处理技术领域。The invention relates to a method for quickly removing artifacts of magnetic resonance imaging of the lower abdomen during free breathing, and belongs to the technical field of image processing.
背景技术Background technique
目前,磁共振成像技术在腹部病变的临床诊断和治疗中起着至关重要的作用。但是磁共振成像时间长,腹部会受呼吸运动的影响。呼吸运动会导致器官发生运动和形变,从而导致图像中产生呼吸运动伪影,使得图像的分辨率和信噪比降低,在图像指导的介入治疗过程中则会出现静态指导信息和运动结构的位置不一致的现象。近年来,有研究者提出一种可以在受试者自由呼吸下成腹部图像的序列,即3D黄金角径向stack-of-stars(SOS)序列。然而在高倍亚采样的情况下,重建图像有明显的条纹伪影。通常采用基于压缩感知的方法来消除该类伪影,然而由于该方法需要迭代计算,重建速度慢,不具有临床应用的可行性。因此,亟需一种自由呼吸下腹部磁共振成像伪影的快速去除方法。Currently, magnetic resonance imaging plays a vital role in the clinical diagnosis and treatment of abdominal lesions. However, MRI takes a long time, and the abdomen will be affected by breathing movements. Respiratory movement will cause movement and deformation of organs, which will lead to respiratory movement artifacts in the image, which will reduce the resolution and signal-to-noise ratio of the image. In the process of image-guided interventional therapy, there will be inconsistencies between the static guidance information and the position of the moving structure. The phenomenon. In recent years, some researchers have proposed a sequence that can form abdominal images under the subject's free breathing, that is, the 3D golden angle radial stack-of-stars (SOS) sequence. However, in the case of high subsampling, the reconstructed image has obvious streak artifacts. A method based on compressed sensing is usually used to eliminate such artifacts. However, since this method requires iterative calculations, the reconstruction speed is slow, and it is not feasible for clinical application. Therefore, there is an urgent need for a rapid removal method for MRI artifacts in the freely breathing lower abdomen.
发明内容Contents of the invention
为了能够解决受试者在自由呼吸下成腹部磁共振成像过程中,消除条纹伪影速度慢,难以临床应用的不足,本发明的目的在于提供一种可以在自由呼吸下降采样的磁共振图像条纹伪影的快速去除方法。In order to be able to solve the deficiencies of slow elimination of streak artifacts and difficulty in clinical application during the subject’s abdominal magnetic resonance imaging under free breathing, the purpose of the present invention is to provide a magnetic resonance image streak that can be down-sampled during free breathing A quick way to remove artifacts.
为了实现上诉目的,本发明的技术方案如下:In order to realize the object of appeal, the technical scheme of the present invention is as follows:
一种自由呼吸下腹部磁共振成像伪影的快速去除方法,由以下步骤组成:A method for rapidly removing artifacts in MRI of the freely breathing lower abdomen, consisting of the following steps:
1)采集受试者在腹部自由呼吸状态下带有条纹伪影的图像数据I0;优选的,所述的采集的磁共振扫描序列为3D黄金角径向stack-of-stars(SOS)扫描序列;优选的,所述的图像数据I0是通过5-20倍的降采样获得。1) Collect the image data I0 of the subject with streak artifacts in the free breathing state of the abdomen; preferably, the magnetic resonance scan sequence collected is a 3D golden angle radial stack-of-stars (SOS) scan sequence ; Preferably, the image data I0 is obtained by downsampling by 5-20 times.
2)将采集到的I0进行预处理得到L个子图像,以便将降采样的图像数据I0,输入到矫正模型M中;优选的,所述的预处理方法采用边长为5-10个像素的方块为滑动窗口,遍历图像I0得到L个子区域,相邻子区域之间的重叠度为70-85%。2) Preprocessing the collected I0 to obtain L sub-images, so that the down-sampled image data I0 is input into the correction model M; The square is a sliding window, L sub-regions are obtained by traversing the image I0, and the overlapping degree between adjacent sub-regions is 70-85%.
3)将L个子图像分别通过事先训练好的校正模型M进行伪影去除,优选的,所述的校正模型M是经过有监督学习得到的非线性映射,由堆叠卷积自编码网络构建;其中校正模型M的样本正确标注数据(ground true)为全K空间采样得到的无伪影图像,模型M的训练样本来自相应正确标注数据的5%-20%部分K空间采样得到的有伪影图像;优选的,样本正确标注数据的磁共振扫描序列为用3D黄金角径向stack-of-stars(SOS)扫描序列,完成K空间全采样。3) performing artifact removal on the L sub-images respectively through a pre-trained correction model M, preferably, the correction model M is a non-linear mapping obtained through supervised learning, and is constructed by a stacked convolutional autoencoder network; wherein The correctly labeled data (ground true) of the corrected model M is an artifact-free image obtained by sampling the full K-space, and the training sample of the model M comes from an artifact-free image obtained by sampling 5%-20% of the corresponding correctly labeled data in the K-space ; Preferably, the magnetic resonance scanning sequence in which the sample is correctly marked with data is a 3D golden angle radial stack-of-stars (SOS) scanning sequence to complete the full sampling of K space.
4)将伪影去除后的L个子图拼接称为图像I1,作为去伪影后的输出。4) The splicing of L subimages after artifact removal is called image I1, which is used as the output after artifact removal.
综上所述,一种自由呼吸下腹部磁共振成像伪影的快速去除方法,主要创新之处在于:通过把全采样数据做为正确标准数据(ground true)以及把降采样数据做为样本输入,构建堆叠卷积自编码网络(Stacked Convolutional Auto-encoders,SCAE),得到了输入输出图像数据之间的非线性映射,进而通过该网络,可以去除由于降采样而造成的条纹伪影,实现在自由呼吸下降采样的磁共振图像条纹伪影的快速去除。To sum up, a fast removal method for MRI artifacts in free breathing lower abdomen, the main innovation is: by using the full sampling data as the correct standard data (ground true) and the down sampling data as the sample input , build a stacked convolutional auto-encoder network (Stacked Convolutional Auto-encoders, SCAE), and get the nonlinear mapping between the input and output image data, and then through this network, the streak artifacts caused by downsampling can be removed, and the Fast Removal of Streak Artifacts in Free-breathing Downsampled Magnetic Resonance Images.
附图说明Description of drawings
图1为本发明的整体流程图;Fig. 1 is the overall flowchart of the present invention;
图2为本发明中一个重建结果图Fig. 2 is a reconstruction result figure among the present invention
具体实施方式Detailed ways
下面通过具体实施例对本发明进行说明,但本发明并不局限于此。The present invention will be described below through specific examples, but the present invention is not limited thereto.
如图1所示,本发明的一种自由呼吸下腹部磁共振成像伪影的快速去除方法,首先采集受试者在腹部自由呼吸状态下带有条纹伪影的图像数据I0;进而将采集到的I0进行预处理得到L个子图像;然后将L个子图像分别通过事先训练好的校正模型M进行伪影去除;最终将伪影去除后的L个子图拼接称为图像I1,作为去伪影后的输出。As shown in Fig. 1, a kind of fast removal method of the MRI artifact of the free breathing lower abdomen of the present invention, at first collect the image data I0 of the subject with the streak artifact under the free breathing state of the abdomen; and then collect the collected I0 is preprocessed to obtain L sub-images; then the L sub-images are respectively removed from the artifacts through the pre-trained correction model M; finally, the L sub-images after the artifact removal are spliced and called image I1, which is used as the post-artifact removal Output.
本实施例中采集受试者在腹部自由呼吸状态下采用的磁共振扫描序列为3D黄金角径向stack-of-stars(SOS)扫描序列,进而通过5倍的降采样获得带伪影的图像数据I0。In this embodiment, the magnetic resonance scanning sequence used by the subject in the free breathing state of the abdomen is a 3D golden angle radial stack-of-stars (SOS) scanning sequence, and then the image with artifacts is obtained by downsampling by 5 times Data I0.
本实施例中预处理方法,采用边长为7个像素的方块为滑动窗口,遍历图像I0后得到L个子领域,且相邻区域之间的重叠度为80%。In the preprocessing method of this embodiment, a square with a side length of 7 pixels is used as a sliding window, and L sub-areas are obtained after traversing the image I0, and the overlapping degree between adjacent areas is 80%.
本实施例中的校正模型M,采用堆叠卷积自编码网络,该网络的样本正确标注数据为全K空间采样得到的无伪影图像,该网络的训练样本来自相应正确标注数据的20%部分K空间采样得到的有伪影图像,样本正确标注数据的磁共振扫描序列为用3D黄金角径向stack-of-stars(SOS)扫描序列,完成K空间全采样。此外该网络由9个隐藏层组成,包括:The correction model M in this embodiment adopts a stacked convolutional autoencoder network. The correctly labeled data of the sample of the network is an artifact-free image obtained by sampling in the full K space, and the training samples of the network come from 20% of the corresponding correctly labeled data. The image with artifacts obtained from K-space sampling, and the MRI scan sequence in which the sample is correctly labeled with data is a 3D golden angle radial stack-of-stars (SOS) scan sequence to complete the full sampling of K-space. In addition, the network consists of 9 hidden layers, including:
1)1个滤波器尺寸为128×5×5的卷积层;1) 1 convolutional layer with a filter size of 128×5×5;
2)1个2×2的最大池化层;2) One 2×2 maximum pooling layer;
3)64个滤波器尺寸5×5的卷积层;3) 64 convolutional layers with a filter size of 5×5;
4)1个2×2的最大池化层;4) One 2×2 maximum pooling layer;
5)64个滤波器尺寸5×5的反卷积层;5) 64 deconvolution layers with a filter size of 5×5;
6)1个2×2的上采样层;6) A 2×2 upsampling layer;
7)128个滤波器尺寸为5×5的反卷积层;7) 128 deconvolution layers with a filter size of 5×5;
8)1个2×2的上采样层;8) A 2×2 upsampling layer;
9)1个滤波器尺寸为5×5的反卷积层。9) 1 deconvolution layer with filter size 5×5.
图2显示为采用本发明得到自由呼吸下腹部磁共振成像伪影的去除图像;图像给出了三个典型受试者的结果图,左列为输入的带有条纹伪影的训练图像,中间列为堆叠卷积自编码网络(SCAE)输出的结果图,右列参考金标准图像(ground true)。结果表明本方法可以较好的删除条纹伪影,与金标准的误差仅出现在腹部边界,集中在脂肪高亮信号区域。Fig. 2 shows to adopt the present invention to obtain the removal image of MRI artifact of free breathing lower abdomen; The image provides the result figure of three typical subjects, and the left column is the training image with streak artifact of input, and the middle The column is the result image output by the stacked convolutional autoencoder network (SCAE), and the right column refers to the gold standard image (ground true). The results show that this method can remove streak artifacts well, and the error with the gold standard only appears in the abdominal border, concentrated in the fat highlight signal area.
Claims (5)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810212643.0A CN110276725B (en) | 2018-03-15 | 2018-03-15 | Method for rapidly removing lower abdominal magnetic resonance imaging artifacts in free breathing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810212643.0A CN110276725B (en) | 2018-03-15 | 2018-03-15 | Method for rapidly removing lower abdominal magnetic resonance imaging artifacts in free breathing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110276725A true CN110276725A (en) | 2019-09-24 |
| CN110276725B CN110276725B (en) | 2022-03-25 |
Family
ID=67958442
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810212643.0A Active CN110276725B (en) | 2018-03-15 | 2018-03-15 | Method for rapidly removing lower abdominal magnetic resonance imaging artifacts in free breathing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110276725B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101995561A (en) * | 2010-10-22 | 2011-03-30 | 南方医科大学 | Reconstruction method of periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) magnetic resonance data based on image domain overlapping |
| CN102005031A (en) * | 2010-11-03 | 2011-04-06 | 宁波鑫高益磁材有限公司 | Method and device for eliminating motion artifact of K spacial sampled data in MRI system |
| CN102232830A (en) * | 2010-04-30 | 2011-11-09 | 西门子(深圳)磁共振有限公司 | Method for water-fat separation through magnetic resonance imaging |
| US20140100124A1 (en) * | 2012-10-04 | 2014-04-10 | Asuragen, Inc. | Diagnostic mirnas for differential diagnosis of incidental pancreatic cystic lesions |
| US20150276909A1 (en) * | 2014-03-25 | 2015-10-01 | Beth Israel Deaconess Medical Center, Inc. (Bidmc, Inc.) | Mri methods and apparatus for flexible visualization of any subset of an enlarged temporal window |
| CN107072595A (en) * | 2013-12-31 | 2017-08-18 | 威斯康星州医药大学股份有限公司 | Adaptive Replanning Based on Multimodal Imaging |
| CN107154025A (en) * | 2016-03-04 | 2017-09-12 | 北京大学 | A kind of blood flow artifact minimizing technology being imaged for arteria carotis magnetic resonance vascular wall |
| US20170322274A1 (en) * | 2016-05-05 | 2017-11-09 | Siemens Healthcare Gmbh | Multi-Echo Pseudo-Golden Angle Stack of Stars Thermometry with High Spatial and Temporal Resolution Using k-Space Weighted Image Contrast |
| CN107507148A (en) * | 2017-08-30 | 2017-12-22 | 南方医科大学 | The method that the down-sampled artifact of MRI is removed based on convolutional neural networks |
-
2018
- 2018-03-15 CN CN201810212643.0A patent/CN110276725B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102232830A (en) * | 2010-04-30 | 2011-11-09 | 西门子(深圳)磁共振有限公司 | Method for water-fat separation through magnetic resonance imaging |
| CN101995561A (en) * | 2010-10-22 | 2011-03-30 | 南方医科大学 | Reconstruction method of periodically rotated overlapping parallel lines with enhanced reconstruction (PROPELLER) magnetic resonance data based on image domain overlapping |
| CN102005031A (en) * | 2010-11-03 | 2011-04-06 | 宁波鑫高益磁材有限公司 | Method and device for eliminating motion artifact of K spacial sampled data in MRI system |
| US20140100124A1 (en) * | 2012-10-04 | 2014-04-10 | Asuragen, Inc. | Diagnostic mirnas for differential diagnosis of incidental pancreatic cystic lesions |
| CN107072595A (en) * | 2013-12-31 | 2017-08-18 | 威斯康星州医药大学股份有限公司 | Adaptive Replanning Based on Multimodal Imaging |
| US20150276909A1 (en) * | 2014-03-25 | 2015-10-01 | Beth Israel Deaconess Medical Center, Inc. (Bidmc, Inc.) | Mri methods and apparatus for flexible visualization of any subset of an enlarged temporal window |
| CN107154025A (en) * | 2016-03-04 | 2017-09-12 | 北京大学 | A kind of blood flow artifact minimizing technology being imaged for arteria carotis magnetic resonance vascular wall |
| US20170322274A1 (en) * | 2016-05-05 | 2017-11-09 | Siemens Healthcare Gmbh | Multi-Echo Pseudo-Golden Angle Stack of Stars Thermometry with High Spatial and Temporal Resolution Using k-Space Weighted Image Contrast |
| CN107507148A (en) * | 2017-08-30 | 2017-12-22 | 南方医科大学 | The method that the down-sampled artifact of MRI is removed based on convolutional neural networks |
Non-Patent Citations (2)
| Title |
|---|
| CHRISTOPHER M. RANK等: ""Artifact model-based respiratory motion compensation (MoCo) for simultaneous PET/MR based on strongly undersampled radial MR data"", 《2014 IEEE NUCLEAR SCIENCE SYMPOSIUM AND MEDICAL IMAGING CONFERENCE (NSS/MIC)》 * |
| 肖韬辉等: ""深度学习的快速磁共振成像及欠采样轨迹设计"", 《中国图象图形学报》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110276725B (en) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111081354B (en) | System and method for denoising medical images through deep learning network | |
| CN110461228B (en) | Improving the quality of medical images using multi-contrast and deep learning | |
| Yu et al. | Computed tomography super-resolution using convolutional neural networks | |
| CN106204489B (en) | The single image super resolution ratio reconstruction method converted in conjunction with deep learning and gradient | |
| US11741580B2 (en) | Machine learning processing of contiguous slice image data | |
| CN112634197A (en) | Method for establishing intracranial angiography enhanced three-dimensional model | |
| CN104931903B (en) | A kind of magnetic resonance removes the method and device of motion artifacts | |
| KR102331104B1 (en) | Novel image reconstruction system for nuclear medicine through training the neural network for improving the spatial resolution and image quality simultaneously based on structural image of phantoms | |
| CN114419015A (en) | Brain function fusion analysis method based on multi-modal registration | |
| CN114758028B (en) | A high-performance reconstruction method for undersampled magnetic resonance images based on optimized regularization and CNN | |
| CN118261798B (en) | A DWI zero-sample super-resolution reconstruction method based on dynamic early stopping | |
| CN113592745B (en) | Unsupervised MRI image restoration method based on antagonism domain self-adaption | |
| CN112037304A (en) | Two-stage edge-enhanced QSM reconstruction method based on SWI phase images | |
| CN112508868A (en) | Intracranial blood vessel comprehensive image generation method | |
| CN114170152A (en) | Method for establishing simulated three-dimensional intracranial vascular stenosis analysis model | |
| Li et al. | Wtfusion: Wavelet-assisted transformer network for multisensor image fusion | |
| Malczewski | PET image reconstruction using compressed sensing | |
| CN110276725A (en) | A fast method for removing artifacts in free-breathing lower abdomen MRI | |
| CN112508881A (en) | Intracranial blood vessel image registration method | |
| CN114240841A (en) | Establishment method of simulated three-dimensional vascular stenosis analysis model | |
| CN108986027A (en) | Depth image super-resolution reconstruction method based on improved joint trilateral filter | |
| Mirza et al. | Super resolution mri via upscaling diffusion bridges | |
| CN116165587A (en) | Quantitative magnetic sensitive imaging method based on deep learning | |
| CN112508877A (en) | Blood vessel comprehensive image generation method | |
| CN115581449A (en) | A multi-b-value combined high-definition diffusion magnetic resonance intelligent imaging method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |