CN113516725B - 基于fpga飞焦点模式下的暗电流智能处理方法 - Google Patents

基于fpga飞焦点模式下的暗电流智能处理方法 Download PDF

Info

Publication number
CN113516725B
CN113516725B CN202110331749.4A CN202110331749A CN113516725B CN 113516725 B CN113516725 B CN 113516725B CN 202110331749 A CN202110331749 A CN 202110331749A CN 113516725 B CN113516725 B CN 113516725B
Authority
CN
China
Prior art keywords
dark current
fpga
ram
data
detector
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.)
Active
Application number
CN202110331749.4A
Other languages
English (en)
Other versions
CN113516725A (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.)
Minfound Medical Systems Co Ltd
Original Assignee
Minfound Medical Systems 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 Minfound Medical Systems Co Ltd filed Critical Minfound Medical Systems Co Ltd
Priority to CN202110331749.4A priority Critical patent/CN113516725B/zh
Publication of CN113516725A publication Critical patent/CN113516725A/zh
Application granted granted Critical
Publication of CN113516725B publication Critical patent/CN113516725B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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/008Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/94Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR 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/10081Computed x-ray tomography [CT]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D10/00Energy efficient computing, e.g. low power processors, power management or thermal management

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Transforming Light Signals Into Electric Signals (AREA)

Abstract

本发明属于图像处理技术领域,具体涉及基于FPGA飞焦点模式下的暗电流智能处理方法。包括FPGA、上位机、滑环和探测器,所述FPGA内含有RAM A和RAM B;包括如下步骤:在探测器开始采集曝光数据前,先对探测器开启飞焦点采数模式,采集探测器飞焦点模式下,两个焦点情况下各个像素点的暗电流大小,将这些值存放在FPGA内的两个RAM中,当开始采集曝光数据后,根据控制板切换球管焦点位置信息,将探测器采集到的数据减去相应像素点的暗电流,最终将处理后的数据打包输出给上位机。本发明具有能够使采集到的飞焦点数据均一性提高,进而提升图像质量的特点。

Description

基于FPGA飞焦点模式下的暗电流智能处理方法
技术领域
本发明属于图像处理技术领域,具体涉及基于FPGA飞焦点模式下的暗电流智能处理方法。
背景技术
目前,CT系统为了获得更高的采样率以及更高的分辨率,通常采用的一种方法是飞焦点模式。飞焦点模式即控制焦点来回切换从而获得2倍或多倍的采样效率,从而使得图像的分辨率更高。
但是,在飞焦点模式下探测器在焦点切换后产生的暗电流会随着焦点的变化而变化,并且每次扫描暗电流会随着探测器温湿度等环境大小的不同而产生变化,如果对暗电流不做相对应的处理,则在飞焦点模式下采集到的数据进行图像重建会产生伪影,影响图像质量。
因此,设计一种能够使采集到的飞焦点数据均一性提高,进而提升图像质量的暗电流处理方法,就显得十分必要。
例如,申请号为CN201910627524.6的中国发明专利所述的一种z向飞焦点扫描方式和图像重建方法,z向飞焦点扫描方式包括Z向飞焦点采样、调节相邻两个焦点的位置、轴扫模式下旋转机架轴180度、Z向飞焦点控制、设计线性变化的斜率;其中图像重建方法包括对两个焦点的数据分别进行重新排列,使其排列成两个平行束;在角度方向进行插值,使两个飞焦点的扫描数据角度相同;在径向方向进行插值,使两个焦点的非等距离平行束变为等距离的平行束;对两个平行束分别进行滤波操作;进行反投影计算。虽然上述的Z向飞焦点轴扫模式和图像重建方法可以在Z向任意位置重建而得到层厚完全相同的图像;同时此特殊轴扫模式和重建方法增加了Z向的采样点数因而减少了Z向混叠伪影,但是其缺点在于,上述的Z向飞焦点轴扫模式和图像重建方法并不能用于处理在飞焦点模式下探测器在焦点切换后产生的暗电流。
发明内容
本发明是为了克服现有技术中,在飞焦点模式下,两个焦点对应暗电流大小不一致,导致飞焦点模式下采集到的数据进行图像重建时会产生伪影,影响图像质量的问题,提供了一种能够使采集到的飞焦点数据均一性提高,进而提升图像质量的基于FPGA飞焦点模式下的暗电流智能处理方法。
为了达到上述发明目的,本发明采用以下技术方案:
基于FPGA飞焦点模式下的暗电流智能处理方法,包括FPGA、上位机、滑环和探测器,所述FPGA内含有RAM A和RAM B;包括如下步骤:
S1,所述上位机发送开启飞焦点模式至FPGA,并开启数据采集模式;
S2,所述上位机对滑环发出旋转指令,所述滑环开始旋转并发送滑环脉冲信号至FPGA;
S3,所述上位机对探测器发送数据采集指令,所述探测器开始采集暗电流数据;
S4,所述FPGA获得探测器采集到的暗电流数据,并根据焦点位置将暗电流数据存放在RAM A和RAM B中;
S5,所述上位机开启曝光模式,并发送曝光指令至FPGA,所述FPGA接收到曝光指令后,锁存住当前时刻下的探测器飞焦点模式下各个像素点的暗电流值;
S6,所述探测器将曝光的数据发送至FPGA,并对所述曝光的数据进行暗电流处理,最终将处理后的数据打包输出给上位机。
作为优选,步骤S2还包括如下步骤:
S21,所述FPGA通过滑环接收到滑环上的孔位信号,转换成采样脉冲;
S22,所述FPGA根据采样脉冲在孔位间插入一个虚拟脉冲,用以切换焦点位置。
作为优选,步骤S3包括如下步骤:
S31,所述FPGA开始根据滑环上的采样点切换飞焦点位置,所述探测器采集暗电流数据,采集的每个像素点的暗电流数据按照view方向依次排列并发送到FPGA中。
作为优选,步骤S4包括如下步骤:
S41,根据焦点的位置采用PING-PONG RAM模式将飞焦点模式下两个焦点的暗电流数据分别存在RAM A和RAM B中;
其中,所述滑环脉冲采集的暗电流数据放置在RAM A中,所述虚拟脉冲采集的暗电流数据放置在RAM B中。
作为优选,所述RAM A和RAM B中的暗电流数据在每次对应的暗电流数据采集后进行刷新。
作为优选,步骤S6中所述暗电流处理,包括如下步骤:
所述FPGA根据焦点位置读取RAM A或RAM B内锁存的探测器飞焦点模式下各个像素点的暗电流值,并将所述曝光的数据各个像素点减去对应像素点的暗电流值后获得的处理数据。
本发明与现有技术相比,有益效果是:(1)本发明通过飞焦点模式下的暗电流处理,使采集到的飞焦点数据均一性提高,进而提升图像质量;(2)本发明采用每次扫描前采集本次扫描条件下每个像素点飞焦点模式下的两个暗电流值,并存储在FPGA中的方式,来提高采集数据的准确性。
附图说明
图1为本发明基于FPGA飞焦点模式下的暗电流智能处理方法的一种流程图;
图2为本发明基于FPGA飞焦点模式下的暗电流智能处理方法的一种原理框图;
图3为本发明中采集飞焦点暗电流阶段的一种原理示意图;
图4为本发明中曝光模式下采数阶段的一种原理示意图;
图5为无飞焦点暗电流处理下的数据的一种效果示意图;
图6为经过本发明基于FPGA飞焦点模式下的暗电流智能处理方法处理后的数据的一种效果示意图。
具体实施方式
为了更清楚地说明本发明实施例,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
实施例1:
如图1和图2所示,基于FPGA飞焦点模式下的暗电流智能处理方法,包括FPGA、上位机、滑环和探测器,所述FPGA内含有RAM A和RAM B;包括如下步骤:
S1,所述上位机发送开启飞焦点模式至FPGA,并开启数据采集模式;
S2,所述上位机对滑环发出旋转指令,所述滑环开始旋转并发送滑环脉冲信号至FPGA;
S3,所述上位机对探测器发送数据采集指令,所述探测器开始采集暗电流数据;
S4,所述FPGA获得探测器采集到的暗电流数据,并根据焦点位置将暗电流数据存放在RAM A和RAM B中;
S5,所述上位机开启曝光模式,并发送曝光指令至FPGA,所述FPGA接收到曝光指令后,锁存住当前时刻下的探测器飞焦点模式下各个像素点的暗电流值;
S6,所述探测器将曝光的数据发送至FPGA,并对所述曝光的数据进行暗电流处理,最终将处理后的数据打包输出给上位机。
进一步的,步骤S2还包括如下步骤:
S21,所述FPGA通过滑环接收到滑环上的孔位信号,转换成采样脉冲;
S22,所述FPGA根据采样脉冲在孔位间插入一个虚拟脉冲,用以切换焦点位置。
所述滑环设有机架,滑环位于机架内上,通过上位机对滑环发出的旋转指令,使得滑环旋转的同时带动了机架的旋转。
进一步的,步骤S3包括如下步骤:
S31,所述FPGA开始根据滑环上的采样点切换飞焦点位置,所述探测器采集暗电流数据,采集的每个像素点的暗电流数据按照view方向依次排列并发送到FPGA中。
所述步骤S3探测器采集的暗电流数据,是指在飞焦点模式下,两个焦点情况下各个像素点的暗电流大小。所述view方向指数据查看的方向。
进一步的,步骤S4包括如下步骤:
S41,根据焦点的位置采用PING-PONG RAM模式将飞焦点模式下两个焦点的暗电流数据分别存在RAM A和RAM B中;
其中,所述滑环脉冲采集的暗电流数据放置在RAM A中,所述虚拟脉冲采集的暗电流数据放置在RAM B中。
具体如图3所示,图中PKG代表数据库文件,图中包含多个RAM A和RAMB,FPGA芯片发送采集积分时间和飞焦点位置信息,并根据焦点位置选择对应的暗电流数据,利用采数脉冲(滑环脉冲和虚拟脉冲)将采集到的暗电流数据(PKG文件)存放在RAM A或RAM B中。其中,一个PKG文件对应一个RAM。
进一步的,所述RAM A和RAM B中的暗电流数据在每次对应的暗电流数据采集后进行刷新。所述过程用于确保数据的实时准确性,保证最后暗电流处理的结果更加准确。
另外,步骤S5的具体过程如图4所示,在曝光采数模式下,当所述FPGA接收到曝光指令后,锁存住当前时刻下的探测器飞焦点模式下各个像素点的暗电流值(即图4中的各个PKG文件),此时的暗电流值最接近本次扫描的两个飞焦点位置下暗电流大小。每个RAM对应一个探测器曝光的数据和一个PKG文件,方便进行后续的数据处理操作。
进一步的,步骤S6中所述暗电流处理,包括如下步骤:
所述FPGA根据焦点位置读取RAM A或RAM B内锁存的探测器飞焦点模式下各个像素点的暗电流值,并将所述曝光的数据各个像素点减去对应像素点的暗电流值后获得的处理数据。
本发明为了解决飞焦点模式下,两个焦点对应暗电流大小不一致的情况,在探测器开始采集曝光数据前,先对探测器开启飞焦点采数模式,采集探测器飞焦点模式下,两个焦点情况下各个像素点的暗电流大小,将这些值存放在FPGA内的两个RAM中,当开始采集曝光数据后,根据控制板切换球管焦点位置信息,将探测器采集到的数据减去相应像素点的暗电流,使得采集上的飞焦点数据均一性提高,图像质量提升。本发明中的FPGA位于对应的数据采集板上。
具体对比结果如图5和图6所示,图5为无飞焦点暗电流处理下的图像中各个像素点的灰度值大小;图6为经过本发明基于FPGA飞焦点模式下的暗电流智能处理方法处理后的图像中各个像素点的灰度值大小。从图5和图6中虚线框部分数据的比较可以看出,经过本发明方法处理后的各个像素点灰度值的均一性,有了明显的提高。经过FPGA处理后的数据明显优于无暗电流处理的数据。
本发明通过飞焦点模式下的暗电流处理,使采集到的飞焦点数据均一性提高,进而提升图像质量;本发明采用每次扫描前采集本次扫描条件下每个像素点飞焦点模式下的两个暗电流值,并存储在FPGA中的方式,来提高采集数据的准确性。
以上所述仅是对本发明的优选实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。

Claims (6)

1.基于FPGA飞焦点模式下的暗电流智能处理方法,包括FPGA、上位机、滑环和探测器,所述FPGA内含有RAM A和RAM B;其特征在于,包括如下步骤:
S1,所述上位机发送开启飞焦点模式至FPGA,并开启数据采集模式;
S2,所述上位机对滑环发出旋转指令,所述滑环开始旋转并发送滑环脉冲信号至FPGA;
S3,所述上位机对探测器发送数据采集指令,所述探测器开始采集暗电流数据;
S4,所述FPGA获得探测器采集到的暗电流数据,并根据焦点位置将暗电流数据存放在RAM A和RAM B中;
S5,所述上位机开启曝光模式,并发送曝光指令至FPGA,所述FPGA接收到曝光指令后,锁存住当前时刻下的探测器飞焦点模式下各个像素点的暗电流值;
S6,所述探测器将曝光的数据发送至FPGA,并对所述曝光的数据进行暗电流处理,最终将处理后的数据打包输出给上位机。
2.根据权利要求1所述的基于FPGA飞焦点模式下的暗电流智能处理方法,其特征在于,步骤S2还包括如下步骤:
S21,所述FPGA通过滑环接收到滑环上的孔位信号,转换成采样脉冲;
S22,所述FPGA根据采样脉冲在孔位间插入一个虚拟脉冲,用以切换焦点位置。
3.根据权利要求1所述的基于FPGA飞焦点模式下的暗电流智能处理方法,其特征在于,步骤S3包括如下步骤:
S31,所述FPGA开始根据滑环上的采样点切换飞焦点位置,所述探测器采集暗电流数据,采集的每个像素点的暗电流数据按照view方向依次排列并发送到FPGA中;
探测器采集的暗电流数据,是指在飞焦点模式下,两个焦点情况下各个像素点的暗电流大小;所述view方向指暗电流数据查看的方向。
4.根据权利要求2所述的基于FPGA飞焦点模式下的暗电流智能处理方法,其特征在于,步骤S4包括如下步骤:
S41,根据焦点的位置采用PING-PONG RAM模式将飞焦点模式下两个焦点的暗电流数据分别存在RAM A和RAM B中;
其中,所述滑环脉冲采集的暗电流数据放置在RAM A中,所述虚拟脉冲采集的暗电流数据放置在RAM B中。
5.根据权利要求4所述的基于FPGA飞焦点模式下的暗电流智能处理方法,其特征在于,所述RAM A和RAM B中的暗电流数据在每次对应的暗电流数据采集后进行刷新。
6.根据权利要求1所述的基于FPGA飞焦点模式下的暗电流智能处理方法,其特征在于,步骤S6中所述暗电流处理,包括如下步骤:
所述FPGA根据焦点位置读取RAM A或RAM B内锁存的探测器飞焦点模式下各个像素点的暗电流值,并将所述曝光的数据各个像素点减去对应像素点的暗电流值后获得的处理数据。
CN202110331749.4A 2021-03-29 2021-03-29 基于fpga飞焦点模式下的暗电流智能处理方法 Active CN113516725B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110331749.4A CN113516725B (zh) 2021-03-29 2021-03-29 基于fpga飞焦点模式下的暗电流智能处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110331749.4A CN113516725B (zh) 2021-03-29 2021-03-29 基于fpga飞焦点模式下的暗电流智能处理方法

Publications (2)

Publication Number Publication Date
CN113516725A CN113516725A (zh) 2021-10-19
CN113516725B true CN113516725B (zh) 2023-11-28

Family

ID=78062146

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110331749.4A Active CN113516725B (zh) 2021-03-29 2021-03-29 基于fpga飞焦点模式下的暗电流智能处理方法

Country Status (1)

Country Link
CN (1) CN113516725B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114334588B (zh) * 2021-12-31 2024-03-22 明峰医疗系统股份有限公司 一种x射线管灯丝预热方法、系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008187255A (ja) * 2007-01-26 2008-08-14 Fujifilm Corp 撮像装置および画像補正方法
CN103366389A (zh) * 2013-04-27 2013-10-23 中国人民解放军北京军区总医院 Ct图像重建方法
JP2016100879A (ja) * 2014-11-26 2016-05-30 オリンパス株式会社 撮像装置、画像処理方法
CN107465874A (zh) * 2017-09-07 2017-12-12 维沃移动通信有限公司 一种暗电流处理方法及移动终端
CN110264541A (zh) * 2019-07-12 2019-09-20 四川明峰医疗科技有限公司 一种z向飞焦点扫描方式和图像重建方法
CN110363825A (zh) * 2019-07-12 2019-10-22 四川明峰医疗科技有限公司 一种z向高分辨率的CT扫描方式和图像重建方法
CN110477943A (zh) * 2019-08-26 2019-11-22 广东明峰医疗科技有限公司 一种利用飞焦点方式实现双能ct的设备和方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6754298B2 (en) * 2002-02-20 2004-06-22 The Regents Of The University Of Michigan Method for statistically reconstructing images from a plurality of transmission measurements having energy diversity and image reconstructor apparatus utilizing the method

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008187255A (ja) * 2007-01-26 2008-08-14 Fujifilm Corp 撮像装置および画像補正方法
CN103366389A (zh) * 2013-04-27 2013-10-23 中国人民解放军北京军区总医院 Ct图像重建方法
JP2016100879A (ja) * 2014-11-26 2016-05-30 オリンパス株式会社 撮像装置、画像処理方法
CN107465874A (zh) * 2017-09-07 2017-12-12 维沃移动通信有限公司 一种暗电流处理方法及移动终端
CN110264541A (zh) * 2019-07-12 2019-09-20 四川明峰医疗科技有限公司 一种z向飞焦点扫描方式和图像重建方法
CN110363825A (zh) * 2019-07-12 2019-10-22 四川明峰医疗科技有限公司 一种z向高分辨率的CT扫描方式和图像重建方法
CN110477943A (zh) * 2019-08-26 2019-11-22 广东明峰医疗科技有限公司 一种利用飞焦点方式实现双能ct的设备和方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Evaluation of two iterative techniques for reducing metal artifacts in computed tomography;F.Edward Boas等;《Radiology》;第259卷(第3期);第894-902 *
两通道按列输出CMOS图像的非均匀性及校正方法;杨国鹏等;《光学精密工程》;第24卷(第8期);第2087-2093页 *
医用数字化 X 射线影像探测器的质量评价探讨;李宝良等;《中国医疗器械信息》;第36-39页 *

Also Published As

Publication number Publication date
CN113516725A (zh) 2021-10-19

Similar Documents

Publication Publication Date Title
US9224193B2 (en) Focus stacking image processing apparatus, imaging system, and image processing system
US7834317B2 (en) Scanning electron microscope and system for inspecting semiconductor device
US20070019847A1 (en) Radiation image processing apparatus, image processing system, radiation image processing method, storage medium, and program
JP5320329B2 (ja) Sem式欠陥観察装置および欠陥画像取得方法
CN113516725B (zh) 基于fpga飞焦点模式下的暗电流智能处理方法
FR2705479A1 (fr) Procédé et appareil de reconnaissance de formes.
US20240046620A1 (en) Fully automated sem sampling system for e-beam image enhancement
FR2791776A1 (fr) Procede et appareil a contraste de tension pour inspecter des semi-conducteurs en utilisant un faisceau de particules a basse tension
US9784698B2 (en) Method for correcting timing skew in X-ray data read out of an X-ray detector in a rolling shutter mode
CN112882057B (zh) 一种基于插值的光子计数非视域三维成像超分辨方法
JP2003023572A (ja) パネル検出器のピクセル置換の方法及び装置
JPH11352011A (ja) フラットパネル表示器検査用画像取得方法、フラットパネル表示器検査用画像取得装置
JP2022528751A (ja) 4d stemを用いて元素組成を判定するシステム、装置、および方法
CN115248197A (zh) 一种三维成像装置和成像方法
JP2005317818A (ja) パターン検査装置およびパターン検査方法
CN116008794A (zh) 高压断路器分合闸速度特性的测试装置、系统、方法
US20220343463A1 (en) Changing the size of images by means of a neural network
TW201035559A (en) Patterned wafer inspection system using a non-vibrating contact potential difference sensor
CN107242878A (zh) 放射线成像装置和操作该装置的方法、放射线成像系统
US7831083B1 (en) Image quality monitoring for substrate inspection
JP5455694B2 (ja) 荷電粒子線装置
CN105161147A (zh) 一种压水堆乏燃料元件三维中子照相无损检测方法
JP3661633B2 (ja) 電子ビームテストシステム及び電子ビームテスト方法
CN101465941B (zh) 全角度图像螺旋失真的自动校正装置及校正方法
JP3633545B2 (ja) 電子ビームテストシステム及び電子ビームテスト方法

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