CN111965691B - 一种pet中的时间游走修正方法 - Google Patents
一种pet中的时间游走修正方法 Download PDFInfo
- Publication number
- CN111965691B CN111965691B CN202010960063.7A CN202010960063A CN111965691B CN 111965691 B CN111965691 B CN 111965691B CN 202010960063 A CN202010960063 A CN 202010960063A CN 111965691 B CN111965691 B CN 111965691B
- Authority
- CN
- China
- Prior art keywords
- time
- correction
- energy
- signal
- signals
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 16
- 238000013508 migration Methods 0.000 title claims abstract description 14
- 230000005012 migration Effects 0.000 title claims abstract description 14
- 238000013178 mathematical model Methods 0.000 claims abstract description 10
- 230000000630 rising effect Effects 0.000 claims description 9
- 239000002245 particle Substances 0.000 claims description 5
- 230000002285 radioactive effect Effects 0.000 claims description 4
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 8
- 238000012360 testing method Methods 0.000 abstract description 8
- 238000004422 calculation algorithm Methods 0.000 abstract description 3
- 238000012795 verification Methods 0.000 abstract description 2
- 238000002600 positron emission tomography Methods 0.000 description 12
- 230000000694 effects Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001934 delay Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000000700 radioactive tracer Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
- G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
- G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
- G01T1/2914—Measurement of spatial distribution of radiation
- G01T1/2985—In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Medical Informatics (AREA)
- Radiology & Medical Imaging (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Biophysics (AREA)
- Optics & Photonics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Nuclear Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明提供了一种PET中的时间游走修正方法,涉及TOF‑PET探测器技术领域,针对单阈值LED电路,本发明提出了一种通过建立简单的数学模型和算法,基于符合数据中能量信息对时间差直接进行在线Time Walk修正的方案。通过获得符合数据、建立数学模型、扫描权重因子的方式,对信号的原时间差进行修正,得到新的时间差;同时还根据得到的权重因子,提出在数据获取***中建立在线修正查找表的方案。本发明利用实际PET***的测量数据进行了测试验证,测试结果表明,通过本发明time walk修正方法,整个***的时间分辨性能得以显著提升。
Description
技术领域
本发明涉及一种PET中的时间游走修正方法,属于TOF-PET探测器技术领域。
背景技术
基于飞行时间(Time-of Flight,TOF)成像原理的正电子发射计算机断层显像(Positron Emission Tomography,PET)是利用晶体(比如BGO或者Lyso等)将放射源释放的伽马射线转换为可见光,然后利用光电转换装置(SiPM或者PMT等)将可见光转换为模拟电信号。前端电路对该模拟电信号进行高精度的能量和时间测量,然后将测量结果送入后端数据获取***(Data acquisition,DAQ)中进行符合判选。最终利用有效信号在符合线上的两个事件的时间差,重建放射源图像。其中,时间测量分辨率性能是决定最终图像质量的最关键因素之一。影响时间测量分辨率的因素很多,比如探测器本身的信噪比,时间-数字转换装置的精度等等。其中,基于前沿定时(Leading Edge Discrimination,LED)技术的时间检出电路以其结构简单和较低噪声等优点在PET***中得以广泛应用。然而,LED由于采用了固定的甄别阈值,造成不同能量信号之间过阈时间不同产生游走效应,称为Time Walk(时间游走)。Time Walk会影响***的时间分辨率性能。
目前,Time Walk的修正较多采用多阈甄别实现,这意味着更复杂的电路结构和更高的成本。
基于此,做出本申请。
发明内容
为了解决现有技术中存在的上述缺陷,本发明提供了一种PET中的时间游走修正方法。
为了实现上述目的,本发明采取的技术方案如下:
一种PET中的时间游走修正方法,包括如下步骤:
(1)获得符合数据:PET***中的探测器模块测量放射源释放的伽马光子的能量、到达时间和位置信息,得到单粒子数据,对单粒子数据中的能量、时间和位置进行判选,得到代表一对伽马光子产生物理信息的符合数据,该符合数据包括一对伽玛光子两个信号的时间差、两个信号的能量值、以及两个信号的位置;
(2)建立数学模型:根据符合数据建立时间游走修正公式,其中该公式中包含权重因子;
(3)扫描权重因子:对根据修正公式对权重因子进行扫描,然后对原时间差进行修正,得到新的时间差;通过对新的时间差数据进行统计分析,得到修正后的***时间分辨率,当修正后的时间分辨率最小时,得到时间游走修正的最佳权重因子。
进一步地,为了计算更加快速且节省资源,还包括步骤(4)建立时间差修正LUT表:在获得最佳权重因子后,将修正公式中的修正系数进行参数化,得到修正查找表LUT。
更进一步地,为了简化修正过程,所述步骤(4)具体为:以某一参考信号能量的到达时间作为参考时间,然后根据修正公式对任意能量信号进行修正时间值计算,形成能量-修正值的LUT表;符合数据中的两个事件根据各自的能量信息查阅LUT得到各自时间修正值,然后做差值,进而得到修正后的时间信息。
进一步地,所述步骤(2)具体包括如下:一对伽马光子两个信号系核信号,因此核信号数学模型的建立主要依据核信号一阶上升沿数学表达式,其信号幅度y表达式如下:
其中,A为常数,表征信号上升沿最终的最大值,τ为信号上升沿特征时间,t为时间;
能量分别是E1和E2的两个不同能量信号,在单一甄别阈值Vth处产生不同的过阈时间分别为T1和T2,根据公式(1),可以计算其时间差T2-T1即时间游走如下式:
其中,k为常数,与信号形状属性有关,A1和A2分别表征两个信号的幅度信息;
对于核信号来讲,在幅度线性良好的***可以近似能量正比与幅度,因此根据(2)式可得到以下修正公式:
其中,ΔT为修正前两个信号的时间差,δt为修正值,ΔT′为修正后两个信号的时间差,ωj为权重因子,j为扫描次序。
进一步地,所述权重因子与信号的形状、过阈点的位置相关,假定所有探测器单元信号特征相同,有相同的权重因子。
本发明的原理和有益技术效果:
(1)本发明针对单阈值LED(前沿定时)电路,采用一种通过建立简单的数学模型和算法,基于符合数据中能量信息对时间差直接进行在线Time Walk修正的方案不需要复杂的电路,成本低。
(2)本发明还根据得到的权重因子,提出在数据获取***中建立在线修正查找表的方案,使计算更加快速,同时可以节省资源。
(3)本发明还利用实际PET***的测量数据进行了初步验证。测试结果表明,通过time walk修正,整个***的时间分辨性能得以显著提升。
附图说明
图1为PET***示意图;
图2为核信号time walk效应示意图;
图3为time walk修正过程示意图;
图4为time walk修正权重因子和时间分辨率扫描测试结果;
图5为time walk修正前后***时间分辨率对比测试结果。
具体实施方式
为了使本发明的技术手段及其所能达到的技术效果,能够更清楚更完善的披露,兹提供了以下实施例,并结合附图作如下详细说明:
本实施例的一种PET中的时间游走修正方法,具体实施步骤主要分为以下几步:获得符合数据、建立数学模型、扫描权重因子、以及建立时间差修正LUT表。下面具体进行描述:
(1)获得符合数据:获得的符合数据来自PET***中的数据获取***(Data-acquisition,DAQ)。如图1,一个PET***一般由多个探测器模块按照环状排列组成。放射源(示踪剂)释放一对几乎同时的背靠背伽马光子(能量为511keV),探测器模块测量伽马光子的能量、到达时间和位置信息,得到单粒子数据(singles)。由于噪声、散射等效应,探测器模块得到的能量谱线呈现一定分布。一般来讲,对于能量分辨率性能较好的***,能量在425-650keV范围内(能量窗)的信号被认为是有效的。由于伽马光子几乎是同时产生的,两个相对的探测器测量到的伽马光子到达时间差必须在一定范围内才能被认为是同一物理事件而非噪声,这个范围称为时间窗(time window)。另外,在实际测量中,根据位置信息推算出的伽马光子响应线(Line-of-Response,LOR)必须在测量视窗(Field-of-view,FOV)范围内,否则这对数据肯定是来自不同的物理事件或者噪声。因此,根据探测器输出的singles数据对能量、时间和位置进行上述判选,可以得到代表一对伽马光子产生物理信息的符合数据(coin)。根据这些符合数据中的时间差信息,TOF-PET通过重建算法即可得到放射源图像。符合数据一般的格式如下表所示:
表格1一般符合数据格式
时间差 | 事件1能量 | 事件2能量 | 事件1位置 | 事件2位置 |
为了提高***的时间分辨,必须要进行探测器单元的时间延迟修正(time delaycorrection)和过阈时间游走修正(time walk correction)。time delay是指由于各个探测器单元本身会时间延迟不同导致的***时间差分布展宽。通过对每一个探测器单元的过阈时间延迟的平均值可以对其进行修正。time walk是指由于信号能量大小不同带来的过阈时间游走效应,即不同能量的信号过阈时间不同,这也会导致时间性能恶化。本发明主要讨论time walk的修正方案,且所使用的测试数据均已经提前进行了timedeleycorrection。
(2)建立数学模型:一对伽马光子两个信号系核信号,核信号数学模型的建立主要依据核信号一阶上升沿数学表达式,其信号幅度y表达式如下:
其中,A为常数,表征信号上升沿最终的最大值,τ为信号上升沿特征时间,t为时间;
如图2所示,当两个具有不同能量的信号(能量分别是E1和E2),其在单一甄别阈值Vth处产生不同的过阈时间T1和T2。根据公式(1),可以计算其时间差T2-T1即时间游走如下式:
其中,k为常数,与信号形状等属性有关。A1和A2分别表征两个信号的幅度信息;
对于核信号来讲,在幅度线性良好的***可以近似能量正比与幅度(E∝A),因此根据(2)式可得到以下修正公式:
其中,ΔT为修正前两个信号的时间差,δt为修正值,ΔT′为修正后两个信号的时间差,ωj为权重因子,j为扫描次序;
权重因子ωj与信号的形状,过阈点的位置等物理性质相关,如何确定最佳的ωj是修正效果的关键。在本发明中,为了简化设计,假定所有探测器单元信号特征相同,有相同的权重因子。
(3)扫描权重因子:根据公式(3),对ωj进行扫描,然后对原时间差数据进行修正,得到新的时间差数据{ΔT′}。通过对新的时间差数据进行统计分析,得到修正后的***时间分辨率。当修正后的时间分辨率最小时,也就是***时间性能最佳时所对应的ωj值,就是time walk修正的最佳权重因子ωbest。
(4)建立时间差修正LUT表是在获得最佳权重因子ωbest后,在DAQ中将公式中的修正系数δt进行参数化,得到修正查找表(Look-up-table,LUT),这样会比使用数字运算处理器(Digital Signal Processing,DSP)直接进行计算更加快速和节省资源。为了简化修正过程,本实施例以511keV能量信号(参考信号能量Eref)的到达时间作为参考时间tref,然后根据修正公式(3)对任意能量信号进行修正时间值计算,形成能量-修正值LUT。符合数据中的两个事件根据各自的能量信息查阅LUT得到各自时间修正值,然后做差值,进而可以得到修正后的时间信息。
修正过程如下式所示,具体实施步骤如图3所示。
ΔT′=ΔT+Δtbest(4)
本实施例在Minfound公司的Scintcare TOF-PET中进行了初步测试验证。其中,时间测量的最小单元(bin size)为15.625ps。图4是对权重因子进行的扫描,扫描结果显示,存在最优修正因子ωbest。图5是在最佳权重因子ωbest条件下的time walk修正前后时间谱结果。测试结果显示,经过time walk修正,***时间分辨率(半高全宽,FWHM)从约450ps优化到约370ps,时间性能得以大大提升。以上内容是结合本发明的优选实施方式对所提供技术方案所作的进一步详细说明,不能认定本发明具体实施只局限于上述这些说明,对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (2)
1.一种PET中的时间游走修正方法,其特征在于包括如下步骤:
(1)获得符合数据:PET***中的探测器模块测量放射源释放的伽马光子的能量、到达时间和位置信息,得到单粒子数据,对单粒子数据中的能量、时间和位置进行判选,得到代表一对伽马光子产生物理信息的符合数据,该符合数据包括一对伽玛光子两个信号的时间差、两个信号的能量值、以及两个信号的位置;
(2)建立数学模型:根据符合数据建立时间游走修正公式,其中该公式中包含权重因子,具体包括如下:
一对伽马光子两个信号系核信号,因此核信号数学模型的建立主要依据核信号一阶上升沿数学表达式,其信号幅度y表达式如下:
其中,A为常数,表征信号上升沿最终的最大值,τ为信号上升沿特征时间,t为时间;
能量分别是E1和E2的两个不同能量信号,在单一甄别阈值Vth处产生不同的过阈时间分别为T1和T2,根据公式(1),可以计算其时间差T2-T1即时间游走如下式:
其中,k为常数,与信号形状属性有关,A1和A2分别表征两个信号的幅度信息;对于核信号来讲,在幅度线性良好的***可以近似能量正比与幅度,因此根据(2)式可得到以下修正公式:
(3)扫描权重因子:对根据修正公式对权重因子进行扫描,然后对原时间差进行修正,得到新的时间差;通过对新的时间差数据进行统计分析,得到修正后的***时间分辨率,当修正后的时间分辨率最小时,得到时间游走修正的最佳权重因子;
(4)建立时间差修正LUT表:在获得最佳权重因子后,将修正公式中的修正系数进行参数化,得到修正查找表LUT;具体为:以某一参考信号能量的到达时间作为参考时间,然后根据修正公式对任意能量信号进行修正时间值计算,形成能量-修正值的LUT表;符合数据中的两个事件根据各自的能量信息查阅LUT得到各自时间修正值,然后做差值,进而得到修正后的时间信息。
2.如权利要求1所述的一种PET中的时间游走修正方法,其特征在于:所述权重因子与信号的形状、过阈点的位置相关,假定所有探测器单元信号特征相同,有相同的权重因子。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010960063.7A CN111965691B (zh) | 2020-09-14 | 2020-09-14 | 一种pet中的时间游走修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010960063.7A CN111965691B (zh) | 2020-09-14 | 2020-09-14 | 一种pet中的时间游走修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111965691A CN111965691A (zh) | 2020-11-20 |
CN111965691B true CN111965691B (zh) | 2022-12-23 |
Family
ID=73393138
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010960063.7A Active CN111965691B (zh) | 2020-09-14 | 2020-09-14 | 一种pet中的时间游走修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111965691B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112882084B (zh) * | 2021-01-18 | 2022-05-24 | 明峰医疗***股份有限公司 | 非线性在线修正方法、***及计算机可读存储介质 |
CN112907664B (zh) * | 2021-02-08 | 2023-04-25 | 明峰医疗***股份有限公司 | 图像像素位置查找方法、***及计算机可读存储介质 |
US11841472B2 (en) | 2021-04-16 | 2023-12-12 | Shanghai United Imaging Microelectronics Technology Co., Ltd. | Devices, systems, and methods for time correction |
CN113835112B (zh) * | 2021-09-23 | 2024-01-02 | 明峰医疗***股份有限公司 | 一种pet中基于fpga的能量修正***及修正方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101238392A (zh) * | 2005-08-04 | 2008-08-06 | 皇家飞利浦电子股份有限公司 | 用于pet的模块化信号处理主链 |
CN108061912A (zh) * | 2017-08-30 | 2018-05-22 | 深圳联影医疗科技有限公司 | 一种pet成像方法、***及探测器 |
CN110327067A (zh) * | 2019-06-10 | 2019-10-15 | 东软医疗***股份有限公司 | 图像重建方法、装置、终端设备及pet*** |
CN111221027A (zh) * | 2020-03-24 | 2020-06-02 | 明峰医疗***股份有限公司 | 用于tof-pet前沿甄别减少时间游走的电路和方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7818047B2 (en) * | 2001-11-09 | 2010-10-19 | Nova R&D, Inc. | X-ray and gamma ray detector readout system |
CN104656115B (zh) * | 2013-11-19 | 2018-04-24 | 苏州瑞派宁科技有限公司 | 一种时间标记组合的方法与*** |
US10527740B2 (en) * | 2014-04-03 | 2020-01-07 | Siemens Medical Solutions Usa, Inc. | Silicon photomultiplier based TOF-PET detector |
CN105212957B (zh) * | 2015-08-25 | 2018-01-16 | 浙江大学 | 一种基于TV Merge的晶体级PET***时间修正方法 |
CN106491153B (zh) * | 2016-12-29 | 2017-10-27 | 上海联影医疗科技有限公司 | 一种pet散射校正方法、pet成像方法及pet成像*** |
US10962662B2 (en) * | 2018-12-06 | 2021-03-30 | Canon Medical Systems Corporation | Method and apparatus to calibrate timing of a detector using multi-channel events |
CN109581461B (zh) * | 2018-12-19 | 2020-05-12 | 中国科学技术大学 | 核脉冲能量测量方法及*** |
-
2020
- 2020-09-14 CN CN202010960063.7A patent/CN111965691B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101238392A (zh) * | 2005-08-04 | 2008-08-06 | 皇家飞利浦电子股份有限公司 | 用于pet的模块化信号处理主链 |
CN108061912A (zh) * | 2017-08-30 | 2018-05-22 | 深圳联影医疗科技有限公司 | 一种pet成像方法、***及探测器 |
CN110327067A (zh) * | 2019-06-10 | 2019-10-15 | 东软医疗***股份有限公司 | 图像重建方法、装置、终端设备及pet*** |
CN111221027A (zh) * | 2020-03-24 | 2020-06-02 | 明峰医疗***股份有限公司 | 用于tof-pet前沿甄别减少时间游走的电路和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111965691A (zh) | 2020-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111965691B (zh) | 一种pet中的时间游走修正方法 | |
US8907290B2 (en) | Methods and systems for gain calibration of gamma ray detectors | |
CN106716179B (zh) | 数字正电子发射断层摄影中的飞行时间校准 | |
US10078009B2 (en) | Pulse processing circuit with correction means | |
US7414246B2 (en) | Achieving accurate time-of-flight calibrations with a stationary coincidence point source | |
Imperiale et al. | On nuclear spectrometry pulses digital shaping and processing | |
US7820975B2 (en) | Achieving accurate time-of-flight calibrations with a stationary coincidence point source | |
US10930031B2 (en) | Synchronization method and system for single event list type data in PET imaging | |
Gundacker et al. | On the comparison of analog and digital SiPM readout in terms of expected timing performance | |
WO2014209972A1 (en) | Normalization correction for multiple-detection enhanced emission tomography | |
JP6381644B2 (ja) | タイムラベル組合せ方法及びシステム | |
WO2014001926A1 (en) | Digital positron emission tomography (dpet) energy calibration method | |
US20160192896A1 (en) | A hybrid method based on simulation and experimental data to normalize pet data | |
WO2011111447A1 (ja) | 放射線計測装置、核医学診断装置 | |
CN111493915B (zh) | 一种用于pet的时间校正方法 | |
CN109581110A (zh) | 判断探测器增益状态的方法及探测器增益平衡校正方法 | |
Rinkel et al. | Equalization method for Medipix3RX | |
CN109507716B (zh) | 一种获取闪烁体探测器能量信息的方法 | |
Schug et al. | Crystal delay and time walk correction methods for coincidence resolving time improvements of a digital-silicon-photomultiplier-based PET/MRI insert | |
CN112883027B (zh) | Pet探测器能量修正方法、***及计算机可读存储介质 | |
CN113376681A (zh) | 基于晶体侧面耦合SiPM的PET探测器及其读出方法 | |
Son et al. | Development and performance evaluation of a time-of-flight positron emission tomography detector based on a high-quantum-efficiency multi-anode photomultiplier tube | |
Torres-Espallardo et al. | Effect of inter-crystal scatter on estimation methods for random coincidences and subsequent correction | |
Korpar et al. | Study of a Cherenkov TOF-PET module | |
Barrio et al. | Time and Energy Characterization of Semi-Monolithic Detectors with Different Treatments and SiPMs Suitable for Clinical Imaging |
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 |