[go: up one dir, main page]

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 PDF

Info

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
Application number
CN201810212643.0A
Other languages
Chinese (zh)
Other versions
CN110276725B (en
Inventor
吕骏
郑亦嘉
张珏
方竞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking University
Original Assignee
Peking University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Peking University filed Critical Peking University
Priority to CN201810212643.0A priority Critical patent/CN110276725B/en
Publication of CN110276725A publication Critical patent/CN110276725A/en
Application granted granted Critical
Publication of CN110276725B publication Critical patent/CN110276725B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/02Neural networks
    • G06N3/04Architecture, e.g. interconnection topology
    • G06N3/045Combinations of networks
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; 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

一种自由呼吸下腹部磁共振成像伪影的快速去除方法A fast method for removing artifacts in free-breathing lower abdomen MRI

技术领域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)

1. a kind of quick minimizing technology for freely breathing lower abdomen magnetic resonance imaging artifact, which is characterized in that by following steps group At:
1) acquisition subject has the image data I0 of jail-bar artifacts under the free respiratory state of abdomen;
2) collected I0 is pre-processed to obtain L subgraph;
3) L subgraph is subjected to artifact removal by the good calibration model M of precondition respectively;
4) the L subgraph splicing after removing artifact is known as image I1, as the output for going pseudo- movie queen.
2. a kind of quick minimizing technology for freely breathing lower abdomen magnetic resonance imaging artifact as described in claim 1, feature Be, when image data I0 of the acquisition subject under the free respiratory state of abdomen used magnetic resonance imaging sequence For 3D Golden Angle radial direction stack-of-stars (SOS) scanning sequence;The image data I0 is adopted by 5-20 times of drop Sample obtains.
3. a kind of quick minimizing technology for freely breathing lower abdomen magnetic resonance imaging artifact as described in claim 1, feature It is, the preprocess method includes: that use side length for the square of 5-10 pixel be sliding window, and traversal image I0 is obtained L sub-regions, the degree of overlapping between adjacent subarea domain are 70-85%.
4. a kind of quick minimizing technology for freely breathing lower abdomen magnetic resonance imaging artifact as described in claim 1, feature It is, the calibration model M is the Nonlinear Mapping obtained by supervised learning, by stacking convolution autoencoder network structure It builds;Wherein the correct labeled data of the sample of calibration model M is the artifact free image that full K spatial sampling obtains, the training sample of model M What this part 5%-20% K spatial sampling from corresponding correct labeled data obtained has artifacts.
5. calibration model M as claimed in claim 4, which is characterized in that the magnetic resonance imaging sequence of the correct labeled data of sample For it is fully sampled to complete the space K with 3D Golden Angle radial direction stack-of-stars (SOS) scanning sequence.
CN201810212643.0A 2018-03-15 2018-03-15 Method for rapidly removing lower abdominal magnetic resonance imaging artifacts in free breathing Active CN110276725B (en)

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)

* Cited by examiner, † Cited by third party
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

Patent Citations (9)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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