CN110363825A - 一种z向高分辨率的CT扫描方式和图像重建方法 - Google Patents

一种z向高分辨率的CT扫描方式和图像重建方法 Download PDF

Info

Publication number
CN110363825A
CN110363825A CN201910627521.2A CN201910627521A CN110363825A CN 110363825 A CN110363825 A CN 110363825A CN 201910627521 A CN201910627521 A CN 201910627521A CN 110363825 A CN110363825 A CN 110363825A
Authority
CN
China
Prior art keywords
focus
focuses
resolution
winged
rays
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
CN201910627521.2A
Other languages
English (en)
Other versions
CN110363825B (zh
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.)
Sichuan Mingfeng Medical Technology Co Ltd
Original Assignee
Sichuan Mingfeng Medical Technology Co Ltd
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 Sichuan Mingfeng Medical Technology Co Ltd filed Critical Sichuan Mingfeng Medical Technology Co Ltd
Priority to CN201910627521.2A priority Critical patent/CN110363825B/zh
Publication of CN110363825A publication Critical patent/CN110363825A/zh
Application granted granted Critical
Publication of CN110363825B publication Critical patent/CN110363825B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/005Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/002D [Two Dimensional] image generation
    • G06T11/003Reconstruction from projections, e.g. tomography
    • G06T11/006Inverse problem, transformation from projection-space into object-space, e.g. transform methods, back-projection, algebraic methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2211/00Image generation
    • G06T2211/40Computed tomography
    • G06T2211/421Filtered back projection [FBP]

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Medical Informatics (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Pulmonology (AREA)
  • Algebra (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明提供了一种z向高分辨率的CT扫描方式和图像重建方法,涉及医学影像技术领域,其中z向高分辨率的CT扫描方式包括z向飞焦点采样、调节相邻两个焦点的位置、选定最佳螺距,螺旋扫模式下旋转机架轴180度;其中图像重建方法包括以下:对两个焦点的数据分别进行重新排列,使其排列成两个平行束,在角度方向进行插值,使两个飞焦点的扫描数据角度相同,在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束,对两个平行束分别进行滤波操作,进行反投影计算。因此本发明在z向两个焦点之间切换的技术和特殊的螺旋扫描螺距相结合,以达到z向4倍采样率,减少混叠伪影使z向分辨率得到提升。

Description

一种z向高分辨率的CT扫描方式和图像重建方法
技术领域
本发明涉及医学影像技术领域,尤其涉及一种z向高分辨率的CT扫描方式和图像重建方法。
背景技术
第三代CT系统的主要组成部分包括Tube(球管),Collimator(限束器),Detector(探测器)。X射线球管发出X光,经过限束器限制形成一个锥形的光束。锥形光束照射到探测器上经过探测器转化为电信号并经过数据采集和转换单元转换为数字信息存储在图像处理系统中。图像处理系统经过一系列的校正算法和图像重建算法生成图像显示在显示器上。从成本和技术成熟程度来说,主流探测器都采用了由很多个探测器模块排列在一个弧形或者多边形的面上来形成整个探测器,每个探测器模块中规则排列一个探测器单元形成的矩阵。
现有探测器像素设计一般都为在扫描平面内(X方向)探测器像素较大,而z向探测器像素较小。在扫描平面内像素较大是为了兼顾几何探测效率和空间分辨率而设计。扫描平面内可以使用探测器1/4偏移,主射线和相隔180度机架旋转角的共轭射线结合做到2倍采样率。或者进一步和X射线焦点平面内摆动技术结合做到4倍采样率。这些技术可以做到扫描平面内的奈奎斯特采样频率大幅提升增加了平面内的空间分辨率并减少了混叠伪影。
虽然探测器像素在z方向较小,但是相比人体组织的细微结构仍然偏大,在螺旋扫描时会由于在z向采样率不够造成z向空间分辨率较低并且由此造成的螺旋扫描风车状伪影增加。可以通过进一步减少探测器像素在z方向的大小来提高采样率,但是减少探测器像素大小会造成探测器像素之间的缝隙在探测器总面积的占比提高。探测器像素之间的缝隙是不能吸收和转化X光的,因此单纯减少探测器像素大小会造成探测器几何效率过低、图像噪声增大、软组织的低密度分辨率降低。
现有技术中为了提高z向分辨率的方式包括:仅减少探测器像素在z方向的大小、单独使用特殊的螺旋扫描螺距、单独使用z方向飞焦点技术,但上述方法分别存在以下缺陷:
1、减少探测器像素在z方向的大小会导致探测器排与排之间的缝隙占探测器总面积的比例增加,因而探测器几何效率减少导致图像噪声增加而低密度分辨率降低。
2、使用特殊的螺旋扫描螺距可以增加z向采样率,但在扫描平面内此措施有效的范围较小,即射线有效交错的平面内扫描FOV减少。
3、使用z向飞焦点技术可以使采样加倍,但是z向分辨率提高的比例理论上不超过两倍。当然,为了进一步提高z向采样率,理论上可以采用z向4个焦点切换,即,可以在z向设置4个焦点位置,而使相邻的4个采样时刻采用不同的4个焦点位置。调节焦点位置可以使相邻采样时刻射线在机架旋转轴上的采样间隔等于探测器像素的1/4因而使采样率提升到原来的4倍。但是为了减少平面内角度方向采样造成的混叠伪影,每个焦点在机架旋转一周需要的采样数目有限制,使用4个焦点会造成采样数目提升到4倍,对探测器的数据采集、滑环的数据传输、重建的计算资源消耗是很大的挑战,技术难度大幅增加,因而成本会大大增加。
基于此,本案由此产生。
发明内容
为了解决现有技术中存在的上述缺陷,本发明提供了一种z向高分辨率的CT扫描方式和图像重建方法,可以达到z向4倍采样率,进而提升z向采样的奈奎斯特频率为现有技术的两倍,从而减少混叠伪影,使z向分辨率得到提升。
为了实现上述目的,本发明采取的技术方案如下:
一种z向高分辨率的CT扫描方式和图像重建方法,
其中z向高分辨率的CT扫描方式包括以下:
(1)z向飞焦点采样:相邻的两个采样采用不同的两个焦点,相应得到从焦点到探测器单元像素中心的两组主射线;
(2)调节相邻两个焦点的位置,使得主射线之间的距离为采用一个焦点时的一半;
(3)选定最佳螺距,螺旋扫模式下旋转机架轴180度:相应得到两组共轭射线;
对于一个N排探测器,所述最佳螺距为M/2N,其中N和M均为自然数,M小于2N;
其中图像重建方法包括以下:
(1)对两个焦点的数据分别进行重新排列,使其排列成两个平行束;
(2)在角度方向进行插值,使两个飞焦点的扫描数据角度相同;
(3)在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束;
(4)对两个平行束分别进行滤波操作;
(5)进行反投影计算。
进一步的,所述图像重建方法的步骤(5)中,反投影时根据需要重建的像素点找出与其最接近的两条主射线和两条共轭射线,在此四条射线中找出距离此像素点最近的两条射线进行线性插值,得到穿过被重建点的射线数值进行反投影。
进一步的,在机架旋转过程中,奇数编号的采样对应一个焦点得到采样组一,偶数编号的采样对应另外一个焦点得到采样组二。
进一步的,所述z向飞焦点扫描方式的步骤中,在机架扫描过程中患者床和球管相位位置保持不变。
本发明的原理及效果:本发明提出的方法是结合飞焦点扫描模式并且使用特殊的螺距的扫描方法,以及在重建阶段的滤波反投影阶段结合两条主射线和两条共轭射线,在最紧邻的两条射线间进行Z向插值的图像重建方法。结合飞焦点和特殊螺距的方式可以使Z向采样率提高到4倍减少混叠伪影。在最紧邻的两条射线间进行Z向插值减少了因为Z向插值所对应的平滑滤波效应造成的Z向空间分辨率的降低。因此本发明在z向两个焦点之间切换的技术和特殊的螺旋扫描螺距相结合,以达到z向4倍采样率,进而提升了z向采样的奈奎斯特频率为现有技术的两倍,从而减少混叠伪影(螺旋扫描中的风车伪影),使z向分辨率得到提升。
附图说明
图1为传统的没有Z向飞焦点的采样中焦点与探测器相对位置示意图(探测器侧视图);
图2为飞焦点模式示意图(探测器侧视图);
图3为螺距为2/16的主射线和共轭射线相对位置示意图(探测器侧视图);
图4为螺距为5/32的z向飞焦点模式扫描主射线和共轭射线相对位置示意图(探测器侧视图)。
标注说明:0度位置1,180度位置2,机架旋转轴3,探测器阵列4,主射线5,共轭射线6。
具体实施方式
为了使本发明的技术手段及其所能达到的技术效果,能够更清楚更完善的披露,兹提供了以下实施例,并结合附图作如下详细说明:
本实施例的一种z向高分辨率的CT扫描方式和图像重建方法,
其中z向高分辨率的CT扫描方式包括以下:
(1)z向飞焦点采样:相邻的两个采样采用不同的两个焦点,相应得到从焦点到探测器单元像素中心的两组主射线;在机架旋转过程中,奇数编号的采样对应一个焦点得到采样组一,偶数编号的采样对应另外一个焦点得到采样组二。z向飞焦点扫描方式的步骤中,在机架扫描过程中患者床和球管相位位置保持不变。
(2)调节相邻两个焦点的位置,使得主射线之间的距离为采用一个焦点时的一半;
(3)选定最佳螺距,螺旋扫模式下旋转机架轴180度:相应得到两组共轭射线;
对于一个N排探测器,最佳螺距为M/2N,其中N和M均为自然数,M小于2N。
如图2所示,采用Z向飞焦点时,在相邻的两个采样采用不同的两个焦点。在机架旋转过程中,奇数编号的采样对应一个焦点而偶数编号的采样对应另外一个焦点。相邻两个焦点的位置调节使得焦点与每个探测器像素中心的连线之间的距离减少到一半,因此Z向采样率提升一倍而使混叠伪影减少表现出来Z向的主观可分辨的线对数目增加。
螺旋扫描时,相隔180度旋转角的两个射线称为互相共轭。如在机架旋转中心附近最紧邻的位置,0度机架旋转角在探测器中心的射线和180度机架旋转角时在探测器中心的射线互为共轭。一条射线称为主射线,另外一条称为共轭射线。当结合Z向飞焦点扫描模式时,0度附近有两个焦点对应的两条主射线。在180度的位置有同样两个焦点对应的两条共轭射线。图3展示了两条主射线和两条共轭射线的相对位置。图3中显示了一个螺旋扫描的螺距为2/16的示例。可以看出对于任意螺距并不能保证共轭射线在Z向交错排列产生4倍的采样率。
为此,本发明提出采用特殊的螺距的Z向飞焦点模式螺旋扫描来增加Z向采样率。对于一个N排探测器,本发明提出使用的最佳螺距为
本发明提出采用特殊的螺距的Z向飞焦点模式螺旋扫描来增加Z向采样率。对于一个N排探测器,本发明提出使用的最佳螺距为M/2N,N和M都是自然数,M小于2N。采用此螺距将使两条主射线和两条共轭射线在Z向相互交错以增加Z向采样率。
本实施例对16排探测器采用螺距为5/32进行螺旋扫描(其中M为5,N为16),如图4所示。
其中图像重建方法包括以下:
(1)对两个焦点的数据分别进行重新排列,使其排列成两个平行束;
(2)在角度方向进行插值,使两个飞焦点的扫描数据角度相同;
(3)在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束;
(4)对两个平行束分别进行滤波操作;
(5)进行反投影计算。反投影时根据需要重建的像素点找出与其最接近的两条主射线和两条共轭射线,在此四条射线中找出距离此像素点最近的两条射线进行线性插值,得到穿过被重建点的射线数值进行反投影。
以上内容是结合本发明的优选实施方式对所提供技术方案所作的进一步详细说明,不能认定本发明具体实施只局限于上述这些说明,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。

Claims (4)

1.一种z向高分辨率的CT扫描方式和图像重建方法,其特征在于:
其中z向高分辨率的CT扫描方式包括以下:
(1)z向飞焦点采样:相邻的两个采样采用不同的两个焦点,相应得到从焦点到探测器单元像素中心的两组主射线;
(2)调节相邻两个焦点的位置,使得主射线之间的距离为采用一个焦点时的一半;
(3)选定最佳螺距,螺旋扫模式下旋转机架轴180度:相应得到两组共轭射线;
对于一个N排探测器,所述最佳螺距为M/2N,其中N和M均为自然数,M小于2N;
其中图像重建方法包括以下:
(1)对两个焦点的数据分别进行重新排列,使其排列成两个平行束;
(2)在角度方向进行插值,使两个飞焦点的扫描数据角度相同;
(3)在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束;
(4)对两个平行束分别进行滤波操作;
(5)进行反投影计算。
2.如权利要求1所述的一种z向飞焦点扫描方式和图像重建方法,其特征在于:所述图像重建方法的步骤(5)中,反投影时根据需要重建的像素点找出与其最接近的两条主射线和两条共轭射线,在此四条射线中找出距离此像素点最近的两条射线进行线性插值,得到穿过被重建点的射线数值进行反投影。
3.如权利要求1所述的一种z向高分辨率的CT扫描方式和图像重建方法,其特征在于:在机架旋转过程中,奇数编号的采样对应一个焦点得到采样组一,偶数编号的采样对应另外一个焦点得到采样组二。
4.如权利要求1所述的一种z向飞焦点扫描方式和图像重建方法,其特征在于:所述z向飞焦点扫描方式的步骤中,在机架扫描过程中患者床和球管相位位置保持不变。
CN201910627521.2A 2019-07-12 2019-07-12 一种z向高分辨率的CT扫描方式和图像重建方法 Active CN110363825B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910627521.2A CN110363825B (zh) 2019-07-12 2019-07-12 一种z向高分辨率的CT扫描方式和图像重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910627521.2A CN110363825B (zh) 2019-07-12 2019-07-12 一种z向高分辨率的CT扫描方式和图像重建方法

Publications (2)

Publication Number Publication Date
CN110363825A true CN110363825A (zh) 2019-10-22
CN110363825B CN110363825B (zh) 2022-12-20

Family

ID=68218939

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910627521.2A Active CN110363825B (zh) 2019-07-12 2019-07-12 一种z向高分辨率的CT扫描方式和图像重建方法

Country Status (1)

Country Link
CN (1) CN110363825B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113017663A (zh) * 2021-01-28 2021-06-25 明峰医疗系统股份有限公司 一种ct扫描数据采集方法及系统、ct扫描仪
CN113516725A (zh) * 2021-03-29 2021-10-19 明峰医疗系统股份有限公司 基于fpga飞焦点模式下的暗电流智能处理方法
CN116483025A (zh) * 2023-04-23 2023-07-25 赛诺威盛科技(北京)股份有限公司 飞焦点模式下的数据采集系统、方法、电子设备及介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5377250A (en) * 1992-08-07 1994-12-27 General Electric Company Reconstruction method for helical scanning computed tomography apparatus with multi-row detector array
US5987157A (en) * 1996-04-10 1999-11-16 Siemens Aktiengesellschaft Image reconstruction method for a multi-row detector computed tomography apparatus operable in a spiral/helical scan mode
CN102727230A (zh) * 2011-04-02 2012-10-17 沈阳东软医疗系统有限公司 Ct扫描图像重建方法及装置
CN103479379A (zh) * 2013-09-06 2014-01-01 沈阳东软医疗系统有限公司 一种倾斜螺旋扫描的图像重建方法及装置
CN106097411A (zh) * 2016-06-29 2016-11-09 沈阳开普医疗影像技术有限公司 Ct扫描机模式、图像重建方法及高分辨ct扫描机
CN106683142A (zh) * 2016-12-28 2017-05-17 上海联影医疗科技有限公司 一种动态螺距ct图像重建方法及装置
CN108720863A (zh) * 2018-02-12 2018-11-02 沈阳晟诺科技有限公司 一种焦点切换式ct机、扫描方法及图像重建方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5377250A (en) * 1992-08-07 1994-12-27 General Electric Company Reconstruction method for helical scanning computed tomography apparatus with multi-row detector array
US5987157A (en) * 1996-04-10 1999-11-16 Siemens Aktiengesellschaft Image reconstruction method for a multi-row detector computed tomography apparatus operable in a spiral/helical scan mode
CN102727230A (zh) * 2011-04-02 2012-10-17 沈阳东软医疗系统有限公司 Ct扫描图像重建方法及装置
CN103479379A (zh) * 2013-09-06 2014-01-01 沈阳东软医疗系统有限公司 一种倾斜螺旋扫描的图像重建方法及装置
CN106097411A (zh) * 2016-06-29 2016-11-09 沈阳开普医疗影像技术有限公司 Ct扫描机模式、图像重建方法及高分辨ct扫描机
CN106683142A (zh) * 2016-12-28 2017-05-17 上海联影医疗科技有限公司 一种动态螺距ct图像重建方法及装置
CN108720863A (zh) * 2018-02-12 2018-11-02 沈阳晟诺科技有限公司 一种焦点切换式ct机、扫描方法及图像重建方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
霍修坤等: "基于单螺旋CT原始数据的三维图像重建的插补算法", 《中国医疗器械杂志》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113017663A (zh) * 2021-01-28 2021-06-25 明峰医疗系统股份有限公司 一种ct扫描数据采集方法及系统、ct扫描仪
CN113017663B (zh) * 2021-01-28 2024-01-16 明峰医疗系统股份有限公司 一种ct扫描数据采集方法及系统、ct扫描仪
CN113516725A (zh) * 2021-03-29 2021-10-19 明峰医疗系统股份有限公司 基于fpga飞焦点模式下的暗电流智能处理方法
CN113516725B (zh) * 2021-03-29 2023-11-28 明峰医疗系统股份有限公司 基于fpga飞焦点模式下的暗电流智能处理方法
CN116483025A (zh) * 2023-04-23 2023-07-25 赛诺威盛科技(北京)股份有限公司 飞焦点模式下的数据采集系统、方法、电子设备及介质
CN116483025B (zh) * 2023-04-23 2024-03-22 赛诺威盛科技(北京)股份有限公司 飞焦点模式下的数据采集系统、方法、电子设备及介质

Also Published As

Publication number Publication date
CN110363825B (zh) 2022-12-20

Similar Documents

Publication Publication Date Title
CN110264541B (zh) 一种z向飞焦点扫描方式和图像重建方法
CN110363825A (zh) 一种z向高分辨率的CT扫描方式和图像重建方法
US9008402B2 (en) X-ray computed tomography apparatus
CN108720863B (zh) 一种焦点切换式ct机、扫描方法及图像重建方法
US7864916B2 (en) X-ray computer tomography apparatus
CN108283503B (zh) 一种ct机、扫描方法及图像重建方法
JP5727277B2 (ja) X線ct装置
EP1713392A1 (en) Computed tomography imaging with pixel staggering and focal spot modulation
Hagen et al. Cycloidal computed tomography
US9095259B2 (en) Method and system for high resolution nutated slice reconstruction using quarter detector offset
US6535572B2 (en) Methods and apparatus for compensating computed tomographic channel ganging artifacts
JP2002209878A (ja) 追加のコンピュータ断層撮影モードを提供する方法及び装置
JP5725885B2 (ja) X線ct装置
JP2010075553A (ja) X線コンピュータ断層撮影装置
CN113017668B (zh) 一种基于非柱面探测器的ct图像重建方法及ct扫描仪
CN110335671A (zh) 一种可调制的ct探测器数据压缩采集方式
CN100536778C (zh) 带有旋转探测模块的计算机断层摄影成像
US9784694B2 (en) X-ray computed tomography apparatus and reconstruction processing method
JP5602293B2 (ja) X線コンピュータ断層撮影装置
JP4643332B2 (ja) X線コンピュータ断層撮影装置
JP2007117717A (ja) X線コンピュータ断層撮影装置
CN115227275A (zh) 一种旋转阻挡器和基于该旋转阻挡器的散射伪影校正方法
JP2010279427A (ja) X線コンピュータ断層撮影装置
JP2006255187A (ja) X線コンピュータ断層撮影装置

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