CN112304987B - 基于光子计数能谱ct的含能材料等效原子序数测量方法 - Google Patents
基于光子计数能谱ct的含能材料等效原子序数测量方法 Download PDFInfo
- Publication number
- CN112304987B CN112304987B CN202011120012.XA CN202011120012A CN112304987B CN 112304987 B CN112304987 B CN 112304987B CN 202011120012 A CN202011120012 A CN 202011120012A CN 112304987 B CN112304987 B CN 112304987B
- Authority
- CN
- China
- Prior art keywords
- atomic number
- equivalent atomic
- photon counting
- energy spectrum
- sample
- 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
- 239000000463 material Substances 0.000 title claims abstract description 26
- 238000001228 spectrum Methods 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000000203 mixture Substances 0.000 claims description 3
- WHBHBVVOGNECLV-OBQKJFGGSA-N 11-deoxycortisol Chemical compound O=C1CC[C@]2(C)[C@H]3CC[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 WHBHBVVOGNECLV-OBQKJFGGSA-N 0.000 claims 1
- 238000002591 computed tomography Methods 0.000 abstract description 6
- 239000012535 impurity Substances 0.000 description 6
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000003550 marker Substances 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 229920002292 Nylon 6 Polymers 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000001419 dependent effect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 238000009659 non-destructive testing Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- -1 polytetrafluoroethylene Polymers 0.000 description 2
- 229920001343 polytetrafluoroethylene Polymers 0.000 description 2
- 239000004810 polytetrafluoroethylene Substances 0.000 description 2
- 239000012925 reference material Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 241000243251 Hydra Species 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000011521 glass Substances 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- QRXWMOHMRWLFEY-UHFFFAOYSA-N isoniazide Chemical compound NNC(=O)C1=CC=NC=C1 QRXWMOHMRWLFEY-UHFFFAOYSA-N 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
- G01N23/046—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material using tomography, e.g. computed tomography [CT]
-
- 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/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pulmonology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Molecular Biology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
本发明提供一种基于光子计数能谱CT的含能材料等效原子序数测量方法。该方法只需利用光子计数能谱CT对待测样品和3种已知材料的标定物进行CT扫描重建,进而由低、高能量区间CT图像的相对比值图像拟合出等效原子序数与相对比值的关系曲线,利用该关系曲线即可进一步计算得到待测样品和标定物的等效原子序数图像。该方法不依赖于光子计数能谱CT的专业知识,实施简便,具有很好的鲁棒性和通用性,可以大大降低现有方法的设备要求和算法复杂度。
Description
技术领域
本发明涉及一种基于光子计数能谱CT的含能材料等效原子序数测量方法,属于工业CT无损检测技术领域。
背景技术
计算机断层成像(Computed Tomography,CT)在医学诊断及工业无损检测中有着广泛的应用。CT图像可看作是材料的线性衰减系数分布,该系数取决于X射线能量和材料特性(质量密度和等效原子序数)。在实际应用中,X射线管产生的光子能谱很宽(这与传统CT理论中X射线源只能发出特定能量光子的假设相违背),积分型探测器对整个能谱的光子进行积分,导致质量密度和等效原子序数均不相同的材料在给定的能谱下可能表现为相同的线性衰减系数,这最终导致传统CT无法有效完成材质的识别。以含能材料检测为例,在实际生产中可以经常发现含能材料样品中含有高线性衰减系数的杂质,其衰减值远高于正常含能材料,甚至高于金属铝。从CT理论分析,高线性衰减系数杂质可能是由于含能材料高聚集造成的(等效原子序数与其他接近,密度远高于正常含能材料),或是由于实际生产过程中混入了高原子序数的杂质。由于传统CT无法有效获得的质量密度和等效原子序数,因此很难判断该类杂质的组成并针对杂质成分改进生产工艺。
作为传统CT的补充,双能CT和能谱CT利用不同能谱的X射线对物体进行CT扫描。结合双能CT和多能CT分解算法,可以获得被扫描物体的有效原子序数和电子密度的分布。但实际应用中上述方法存在能谱测量困难,参考材料标定过程复杂,且可能放大噪声等缺点。据此基于同步辐射以氧气为标准元素的迭代方法被提出,但是该技术设备要求高,难以满足实际需要。后续基于光子计数探测器的标定方法相继被提出,但上述方法必须设置很窄的能量区间(0.5keV或1keV),严重降低了光子计数探测器有效能量区间内接收的光子数,导致CT扫描十分耗时。此外,待测材料的原子序数精度仍可能受重建误差、探测器响应准确性、散射效应等因素影响。
现有的技术缺点主要包括:
(1)双能CT和能谱CT存在能谱测量困难、算法复杂度高、参考材料标定过程复杂,且可能放大噪声等问题。
(2)基于迭代的同步辐射测量方法设备要求高,难以广泛应用。
(3)基于光子计数探测器的标定的方法扫描效率低,难以满足工业需要,且精度易受重建误差、探测器响应准确性、散射效应等因素影响。
综上所述,现有的技术方法存在算法复杂、设备要求高、扫描效率低等问题,难以满足实际含能材料生成中的高效检测需求。
发明内容
为提高实际含能材料生成中等效原子序数的识别精度,并降低设备要求和算法复杂度,本发明提供一种基于光子计数能谱CT的等效原子序数测量方法。
本发明解决其技术问题所采用的技术方案包括以下步骤:
步骤1:获取3个材质组成已知的标定物的等效原子序数Zk,eff,k为标定物序号;
步骤2:利用光子计数能谱CT,将待测样品和3个标定物同时进行CT扫描,并重建低能区CT图像SL(x,y)和高能区CT图像SH(x,y);
步骤3:将SL(x,y)除以SH(x,y),得到待测样品和标定物的相对比值图像SR(x,y);
步骤4:利用标定物的等效原子序数及其相对比值,拟合等效原子序数与相对比值的关系曲线;
步骤5:利用SR(x,y)和步骤4的关系曲线,计算得到待测样品和标定物的等效原子序数图像Zeff(x,y)。
在上述步骤1中,设第k个标定物由原子序数为Z1,Z2,…,Zn的原子组成,根据Glasser公式,其等效原子序数可表示为
其中ai为标定物中第i种原子的质量百分比,Ai为第i种原子的原子质量,ni为第i种原子的原子数目,Zi为第i种原子的原子序数。
在上述步骤3中,当实际图像中的噪声较大时,需要先对SL(x,y)和SH(x,y)进行降噪。为避免异常点,对于SH(x,y)中灰度值小于T的点不进行计算,直接将SR(x,y)中的对应点置为0,其中T为灰度阈值。
在上述步骤4中,利用标定物的等效原子序数及其相对比值,拟合等效原子序数与相对比值的关系曲线为其中Zk,eff是步骤1中获取的标定物等效原子序数,R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数。
在上述步骤5中,根据步骤4的关系曲线R,可以得到等效原子序数图像计算公式再根据SR(x,y)即可计算得到待测样品和标定物的等效原子序数图像Zeff(x,y),其中a、b、c的取值与步骤4中相同。
本发明提供的上述方法,所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用,即此时只需对待测样品(无标定物)执行步骤2、步骤3和步骤5,即可得到待测样品的等效原子序数图像。
本发明的有益效果是:该方法在实际应用中不依赖于光子计数能谱CT的专业知识,实施简便,具有很好的鲁棒性和通用性,可以大大降低现有方法的设备要求和算法复杂度,并且所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用。
下面结合附图和实施例对本发明做进一步说明。
附图说明
图1为本发明算法示意图。
具体实施方式
利用微焦点光子计数能谱CT(光子计数探测器为XCounter Hydra FX50)应用本发明,执行以下步骤:
步骤1:选择3个材质组成已知的标定物,依次为尼龙6、聚四氟乙烯、铝,获取其等效原子序数Zk,eff,k为标定物序号,根据公式:
可计算得到3种标定物的等效原子序数分别为6.6008、8.6163和13.0000。
步骤2:利用微焦点光子计数能谱CT,选择射线源管电压100kV和管电流70μA,同时扫描标定物和待测样品,扫描几何参数为:射线源到探测器距离420mm,射线源到旋转中心距离175mm,探测器像素尺寸为0.10mm,探测器采集速度为20帧/s,20帧图像求和得到一幅投影;低能量阈值为40keV,高能量阈值为60keV;单幅投影大小为1024×64,投影幅数为900幅,重建矩阵大小为1024×1024×64。利用FDK算法分别重建低能区CT图像SL(x,y)和高能区CT图像SH(x,y)。
步骤3:将SL(x,y)除以SH(x,y),得到待测样品和标定物的相对比值图像SR(x,y)。考虑到实际图像中噪声较大,先对SL(x,y)和SH(x,y)利用BM3D算法进行降噪。为避免异常点,对于SH(x,y)中灰度值小于T的点不进行计算,直接将SR(x,y)中的对应点置为0,其中T为灰度阈值,本实例中T取0.01。
步骤4:根据SR(x,y)获得尼龙6、聚四氟乙烯、铝3种标定物的相对比值分别为1.0912、1.1684和1.4854。获取3种标定物的相对比值时,考虑到SR(x,y)中存在噪声,分别取各标定物的多个像素的灰度平均值。利用标定物的等效原子序数和相对比值拟合关系曲线其中R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数,拟合可得关系曲线参数为a=6.375,b=23390,c=125100。
步骤5:根据步骤4的关系曲线R,可以得到等效原子序数图像计算公式再根据SR(x,y)即可计算得到待测样品和标定物的等效原子序数图像Zeff(x,y),其中a、b、c的取值与步骤4中相同。
本发明提供的上述方法,所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用,即此时只需对待测样品(无标定物)执行步骤2、步骤3和步骤5,即可得到待测样品的等效原子序数图像。
Claims (2)
1.一种基于光子计数能谱CT的含能材料等效原子序数测量方法,其特征在于包括下述步骤:
步骤1:获取3个材质组成已知的标定物的等效原子序数Zk,eff,k为标定物序号,设第k个标定物由原子序数为Z1,Z2,…,Zn的原子组成,其等效原子序数表示为
其中ai为标定物中第i种原子的质量百分比,Ai为第i种原子的原子质量,ni为第i种原子的原子数目,Zi为第i种原子的原子序数;
步骤2:利用光子计数能谱CT,将待测样品和3个标定物同时进行CT扫描,并重建低能区CT图像SL(x,y)和高能区CT图像SH(x,y);
步骤3:将SL(x,y)除以SH(x,y),得到待测样品和标定物的相对比值图像SR(x,y);
步骤4:利用标定物的等效原子序数Zk,eff及其相对比值,拟合等效原子序数与相对比值的关系曲线R为相应标定物在SR(x,y)中的相对比值,a、b、c为关系曲线参数;
步骤5:利用SR(x,y)和步骤4的关系曲线R,计算得到待测样品和标定物的等效原子序数图像其中a、b、c的取值与步骤4中相同。
2.根据权利要求1所述的一种基于光子计数能谱CT的含能材料等效原子序数测量方法,其特征在于:该方法所得到的等效原子序数图像计算公式,在同一光子计数能谱CT和相同扫描参数的情况下可以重复利用,即此时只需对待测样品执行步骤2、步骤3和步骤5,即可得到待测样品的等效原子序数图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011120012.XA CN112304987B (zh) | 2020-10-19 | 2020-10-19 | 基于光子计数能谱ct的含能材料等效原子序数测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011120012.XA CN112304987B (zh) | 2020-10-19 | 2020-10-19 | 基于光子计数能谱ct的含能材料等效原子序数测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112304987A CN112304987A (zh) | 2021-02-02 |
CN112304987B true CN112304987B (zh) | 2024-04-12 |
Family
ID=74328283
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011120012.XA Active CN112304987B (zh) | 2020-10-19 | 2020-10-19 | 基于光子计数能谱ct的含能材料等效原子序数测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112304987B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114403912A (zh) * | 2022-01-17 | 2022-04-29 | 明峰医疗系统股份有限公司 | 一种基于深度学习的能谱ct有效原子序数估计方法 |
CN114813798B (zh) * | 2022-05-18 | 2023-07-07 | 中国工程物理研究院化工材料研究所 | 用于表征材料内部结构及成分的ct检测装置和成像方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2014268157A1 (en) * | 2011-07-29 | 2014-12-11 | The Johns Hopkins University | Differential phase contrast X-ray imaging and X-ray Illumination Methods |
CN104346820A (zh) * | 2013-07-26 | 2015-02-11 | 清华大学 | 一种x光双能ct重建方法 |
CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
CN105559813A (zh) * | 2014-11-04 | 2016-05-11 | 株式会社东芝 | 医用图像诊断装置以及医用图像处理装置 |
CN106353828A (zh) * | 2015-07-22 | 2017-01-25 | 清华大学 | 在安检系统中估算被检查物体重量的方法和装置 |
CN107356615A (zh) * | 2016-05-10 | 2017-11-17 | 清华大学 | 一种用于双能x射线ct的方法和系统 |
CN111707688A (zh) * | 2020-06-18 | 2020-09-25 | 清华大学 | 光子计数能谱ct成像中自适应能谱优化方法及其应用 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7190757B2 (en) * | 2004-05-21 | 2007-03-13 | Analogic Corporation | Method of and system for computing effective atomic number images in multi-energy computed tomography |
US10451568B2 (en) * | 2014-08-22 | 2019-10-22 | Canon Medical Systems Corporation | Photon counting X-ray CT apparatus |
-
2020
- 2020-10-19 CN CN202011120012.XA patent/CN112304987B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2014268157A1 (en) * | 2011-07-29 | 2014-12-11 | The Johns Hopkins University | Differential phase contrast X-ray imaging and X-ray Illumination Methods |
CN104346820A (zh) * | 2013-07-26 | 2015-02-11 | 清华大学 | 一种x光双能ct重建方法 |
CN104970815A (zh) * | 2014-04-04 | 2015-10-14 | 曹红光 | 基于光栅相位衬度和光子计数的x射线成像系统及方法 |
CN105559813A (zh) * | 2014-11-04 | 2016-05-11 | 株式会社东芝 | 医用图像诊断装置以及医用图像处理装置 |
CN106353828A (zh) * | 2015-07-22 | 2017-01-25 | 清华大学 | 在安检系统中估算被检查物体重量的方法和装置 |
CN107356615A (zh) * | 2016-05-10 | 2017-11-17 | 清华大学 | 一种用于双能x射线ct的方法和系统 |
CN111707688A (zh) * | 2020-06-18 | 2020-09-25 | 清华大学 | 光子计数能谱ct成像中自适应能谱优化方法及其应用 |
Non-Patent Citations (3)
Title |
---|
"Precise material identification method based on a photon counting technique with correction of the beam hardening effect in X-ray spectra";N. Kimoto et al.;《Applied Radiation and Isotopes》(第124期);16-26 * |
"基于光子计数能谱CT的含能材料等效原子序数测量方法";杨亚飞等;《光谱学与光谱分析》;第42卷(第5期);1400-1406 * |
"多能谱X 射线成像技术及其在CT 中的应用";郝佳等;《CT 理论与应用研究》;第20卷(第1期);141-150 * |
Also Published As
Publication number | Publication date |
---|---|
CN112304987A (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9687207B2 (en) | Pre-reconstruction calibration, data correction, and material decomposition method and apparatus for photon-counting spectrally-resolving X-ray detectors and X-ray imaging | |
US9025815B2 (en) | System and method for multi-material correction of image data | |
Beyer et al. | The use of X-ray CT for attenuation correction of PET data | |
EP2377467A1 (en) | System and method for determining the composition of an object | |
CN112304987B (zh) | 基于光子计数能谱ct的含能材料等效原子序数测量方法 | |
JP2016147055A (ja) | X線ct装置及び補正方法 | |
JP2016127907A (ja) | 医用情報処理装置、x線ct装置、医用情報処理方法、及び医用情報処理プログラム | |
WO2009102990A2 (en) | System and method for quantitative imaging of chemical composition to decompose multiple materials | |
CN110084864A (zh) | 一种基于能谱ct的电子密度图像重建方法 | |
CN108603849B (zh) | 用于确定x射线管的有效能谱的装置 | |
US11879856B2 (en) | Two-step material decomposition calibration method for a full size photon counting computed tomography system | |
Clements et al. | Multi-contrast CT imaging using a high energy resolution CdTe detector and a CZT photon-counting detector | |
US7050530B2 (en) | Method and apparatus for the spatially-resolved determination of the element concentrations in objects to be examined | |
US10799192B2 (en) | Method and apparatus for partial volume identification from photon-counting macro-pixel measurements | |
Chi et al. | Experimental feasibility of dual-energy computed tomography based on the Thomson scattering X-ray source | |
JP2022113115A (ja) | ビームハードニング補正方法、x線ct装置及びビームハードニング補正プログラム | |
Jumanazarov et al. | Significance of the spectral correction of photon counting detector response in material classification from spectral x-ray CT | |
Baek et al. | Diagnostic X-ray spectra detection by Monte Carlo simulation | |
Bisogni et al. | Experimental study of beam hardening artifacts in photon counting breast computed tomography | |
CN106373100B (zh) | 一种ct图像的硬化伪影校正方法 | |
Rosso et al. | Preliminary study of the advantages of X-ray energy selection in CT imaging | |
US11300696B2 (en) | Radiation imaging device and photon counting type detector calibration method | |
EP4400872A1 (en) | Method for calibrating detector and apparatus for calibrating detector | |
US20240230934A1 (en) | Pcct energy calibration from x-ray tube spectra using a neural network | |
US12102471B2 (en) | Iterative method for material decomposition calibration in a full size photon counting computed tomography system |
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 |