CN104819941B - 一种多波段光谱成像方法 - Google Patents
一种多波段光谱成像方法 Download PDFInfo
- Publication number
- CN104819941B CN104819941B CN201510230474.XA CN201510230474A CN104819941B CN 104819941 B CN104819941 B CN 104819941B CN 201510230474 A CN201510230474 A CN 201510230474A CN 104819941 B CN104819941 B CN 104819941B
- Authority
- CN
- China
- Prior art keywords
- data
- multiband
- image
- image sensor
- multispectral
- 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
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Color Television Image Signal Generators (AREA)
Abstract
一种多波段光谱成像方法,在由光学成像器件和彩色摄像传感器构成的光路中,加入多波段带通滤光片,所述多波段带通滤光片置于所述光学成像器件和彩色摄像传感器的光路上,并且所述彩色摄像传感器进行光谱成像形成混叠的多波段多光谱图像,通过FPGA处理模块将混叠的多波段多光谱图像,按预先建立的模型,实时还原成独立的多波段图像送入计算机。采用本发明可以得到2‑3个波段的实时多光谱成像。实现15帧/秒以上的三波段多光谱成像,满足大部分医学多光谱成像应用的要求。
Description
技术领域
本发明涉及一种多波段光谱成像方法,属于医学光谱成像技术领域,具体用于细胞成分的定位和定量分析,特别是医学领域的细胞内组分的定量分析。
背景技术
多光谱成像能够在较复杂的环境下定量分析细胞内的成分信息,如复染环境中细胞核内DNA的含量,对疾病诊断具有重要意义。
多光谱成像主要有两大类技术,扫描多光谱成像技术和凝视式多光谱成像技术。
扫描多光谱成像技术涉及逐点或逐行扫描,成像速度较慢,主要用在空间卫星上。凝视式多光谱成像目前有电调谐滤光器和旋转滤光轮两种,电调谐滤光器能够快速灵活选择需要的波段,成像较快但器件成本很高。旋转滤光轮每成像一个波段,需旋转更换一次滤光片,因此速度较慢,用于医学细胞图像扫描时,每例需要很长的时间,难以投入实际应用。
本发明提出了一种低成本,快速的多光谱成像方法,采用本方法和适当的电子技术,可以实现15帧/秒以上的三波段多光谱成像。满足大部分医学多光谱成像应用的要求。
发明内容
本发明的目的是提供一种多波段光谱成像方法,解决医学领域快速多光谱成像的问题,实现15帧/秒以上的三波段多光谱成像,满足大部分医学多光谱成像应用的要求。
本发明的技术方案是:
一种多波段光谱成像方法,在由光学成像器件和彩色摄像传感器构成的光路中,加入多波段带通滤光片,所述多波段带通滤光片置于所述光学成像器件和彩色摄像传感器的光路上,并且所述彩色摄像传感器进行光谱成像形成混叠的多波段多光谱图像,通过FPGA芯片将混叠的多波段多光谱图像,按预先建立的模型,实时还原成独立的多波段图像送入计算机。
所述光路还包括照明光源,所述照明光源置于所述光学成像器件的前端,并且在所述照明光源和所述光学成像器件之间装有样品。
所述FPGA芯片可采用DSP数字信号处理芯片代替。
FPGA芯片选用ALTERA公司的Cyclone IV E系列的EP4CE40F23C6。
所述FPGA芯片的内部包括图像数据获取与存储模块、同步解码模块、同步矩阵剥离运算模块、以及数据传输模块。
所述FPGA芯片对多路并行处理的流程包括:
a、根据彩色图像传感器的输出时序和输出信号拼接为一个128位的数据,并根据行序列识别模块,判断该行图像是G/R行还是B/G行,进而判定将此128位数据是存入存储G/R行图像的FIFO中还是存入存储B/G行图像的FIFO中;
b、在获取了一组G/R和B/G图像(即一行G/R图像和一行B/G图像)后,同步解码模块对图像数据进行插值运算,最终输出数据为256位的数据,该数据包含4组像素的R、G、B通道的信息,在同步矩阵剥离运算中,每一组的运算包含9个硬件乘法器M和3个硬件除法器,乘法器的一个乘数为计算出的图像采样值,另一个乘数为逆矩阵计算出的逆矩阵T-1的一个光强参数,由于逆矩阵中的参数都是小数,若直接使用这些小数计算会加大运算难度,因此通过还原成整数的乘除法解决该问题,最终计算出4个像素各自的三个波段的光强参数;
c、FPGA芯片中的FIFO存储器可以配置为256位数据输入,16位数据输出的模式,低位数据首先输出,高位数据后输出,同步传输模块将16位数据从FIFO中读出,然后将该数据以符合USB传输时序的方式写入到数据传输芯片CY7C68013中,然后传输到计算机。
所述模型的建立,包括以下步骤:
步骤a、将定制的只有一个波段的窄带滤光片放入光路中,窄带曲线对应多波段带通滤光片的一个通道;
步骤b、对空白样品进行多光谱成像;测出彩色图像传感器的响应x1(x1r,x1g,x1b);
重复步骤a、b,依次将定制的另一个波段的窄带滤光片放入光路中,测出响应xi(xir,xig,xib);
测量出所有波段响应xi(r,g,b),然后综合成敏感度矩阵X,X=[x1,x2,……xn];
计算得到模型M=(XtX)-1Xt。
以上所述每个窄带滤光片的滤光曲线只能覆盖多波段滤光片的一个通道。
采用本发明可以得到2-3个波段的实时多光谱成像。实现15帧/秒以上的三波段多光谱成像,满足大部分医学多光谱成像应用的要求。
附图说明
图1是本发明提出的低成本快速的多光谱成像方法框图;
图2是本发明采用的多波段带通滤光片滤光曲线;
图3是彩色摄像传感器的像素结构;
图4是摄像传感器bayes滤光膜阵列的滤光曲线;
图5是本发明用于复染环境下DNA定量测量的一个显微多光谱成像实例框图;
图6是本发明FPGA实时数据处理流程图;
图7是为充分利用FPGA并行处理能力进行的一次拼接处理8个像素示意图,由于彩色图像传感器像素结构(见图4),必须两行一起处理;
图8是彩色图像传感器每行输出的彩色通道时序图;
图9是两行、16像素、256位一起处理的数据处理示意框图;
图10是采用伊红、亮绿、福尔根染料复染时,宫颈脱落细胞在显微镜下放大成像的图像;
图11采用本发明的三波段实时光谱成像。
具体实施方式
下面结合附图对本发明的技术方案进行具体阐述,需要指出的是,本发明的技术方案不限于实施例所述的实施方式,本领域的技术人员参考和借鉴本发明技术方案的内容,在本发明的基础上进行的改进和设计,应属于本发明的保护范围。
本发明实施例所述的本发明提出的低成本快速的多波段光谱成像方法,基于彩色摄像传感器和定制的多波段(2-3个波段)带通滤光片。
提出的多光谱成像方法如图1所示:被光源照明的样品通过成像光学成像在彩色摄像传感器上,光线在到达传感器前,先通过多波段带通滤光片。多波段带通滤光片可以在图1所示的位置,也可以在光路的任何其他位置。
多波段带通滤光片的透过率按照多光谱成像的要求,设计成2-3个特定波段带通的滤光曲线,可以如图2所示。
彩色摄像传感器检测的图像数据,在DSP数字信号处理芯片或FPGA芯片中按照预先建立的模型被实时还原成多光谱图像送入计算机。
光学检测中,物质的含量是根据Lambert-Beer定律来检测的。
Lambert-Beer定律表达式为:
A=-lg(I/I0)=εbC
式中,I0为波长为λ的平行、均匀入射单色光束强度;I是透过目标后的光束强度;A为光度度;ε为待测组分的摩尔吸光系数;b为光程;C为待测组分的物质的量浓度。
因此单色光的吸光度正比于物质的含量,物质的含量越多,吸收光越多,光的透射越低。Lambert-Beer成立的前提是在检测光的波段中,待测组分的摩尔吸光系数ε为常数。理想的是单波长的单色光,实际用的是和待测组分的摩尔吸光系数变化相比波段足够窄的带通光。
彩色摄像传感器是在传感器前加一个Bayer滤光膜阵列构成。模式如图3所示,四个像素一组,一个红色R、两个绿色G和一个蓝色B。Bayer滤光膜的透过率如图4所示。
彩色摄像机为我们带来了比单色摄像机更丰富的图像信息。但从图4我们可以看出,RGB三色传感器的透过曲线波段非常宽,而且彼此重叠,因此它检测的数据并不符合Lambert-Beer定律成立的条件。而且滤光曲线的峰值往往不在我们需要的波长。
本发明设计了一个多波段带通滤光片来克服这一缺陷。多波段带通滤光片的波长中心在测量所需要的波长上,带宽按满足Lambert-Beer定律成立条件设计。
通过多波段带通滤光片后光线只剩下所需要的窄带光。剩下的问题是RGB三色传感器对所有波段的光都有响应。因此本发明在彩色摄像传感器后设计了一个由DSP或FPGA电子器件构成的部件,按照建立的模型剥离彩色摄像传感器的混合响应,得到多波段单独响应的多光谱图像。
模型的建立包括以下步骤:
步骤a、将定制的只有一个波段的窄带滤光片放入图1所示的光路,窄带曲线对应多波段带通滤光片的一个通道;
步骤b、对空白样品采用本发明的多光谱成像系统进行多光谱成像;测出RGB传感器的响应x1(x1r,x1g,x1b);
重复步骤a、b,依次将定制的另一个波段的窄带滤光片放入光路,测出响应xi(xir,xig,xib);
测量出所有波段响应xi(r,g,b),然后综合成敏感度矩阵X,X=[x1,x2,……xn];
计算得到模型M=(XtX)-1Xt;
采用本发明可以得到2-3个波段的实时多光谱成像。
如图5所示是采用本发明搭建的多光谱显微成像仪的原理框图。光源1发出的白光照亮载玻片2上染色后的脱落细胞,脱落细胞位于显微镜物镜3的焦点上,因此经过脱落细胞上的点的光线,在经过显微镜物镜3后变成平行光,再经显微镜成像镜4将细胞图像成像在彩色图像传感器6的靶面上,在彩色图像传感器6的靶面前插入了一个多波段带通滤光片5,因此此时彩色摄像机上采集到的是混叠的多波段多光谱图像。混叠的多波段多光谱图像通过FPGA芯片7或DSP数字信号处理芯片,按预先建立的模型,实时还原成独立的多波段图像送入计算机8。
混叠的多波段多光谱图像通过DSP数字信号处理芯片或FPGA芯片7和彩色图像传感器6(摄像机)集成在一起。
FPGA芯片选用ALTERA公司的Cyclone IV E系列的EP4CE40F23C6型号,该芯片含有39600个逻辑单元,1134Kbits的嵌入式存储器。数据传输芯片采用CYPRESS公司的CY7C68013,此芯片集成了USB 3.0收发器、增强型8051微处理器以及智能的SIE(SerialInterface Engine),符合USB3.0规范,理论最大数据传输速率为480Mb/s。FPGA通过模拟IIC总线与MT9T031通信,用于配置MT9T031的工作模式、输出图像大小和速率、增益以及曝光时间等。
主要的数据获取和处理过程都是通过FPGA内部的各个功能模块完成的,主要包括:图像数据获取与存储模块,同步解码模块,同步矩阵剥离运算,以及数据传输模块,如图6所示,整个硬件算法模块采用流水线操作以达到图像输出实时化。FPGA内部FIFO(FirstInput First Output)存储器的IP核模块能处理的最大位宽为256位,长位宽为实时多路并行处理多流水线的数据处理提供了可能。
如图7-9所显示了整个多路并行流水线处理流程:
a、根据彩色图像传感器的输出时序和输出信号拼接为一个128位的数据,如图7所示,并根据图8中的行序列识别模块,判断该行图像是G/R行还是B/G行,以此来判定将此128位数据是存入存储G/R行图像的FIFO中还是存入存储B/G行图像的FIFO中。
b、在获取了一组G/R和B/G图像(即一行G/R图像和一行B/G图像)后,同步解码模块对图像数据进行插值运算,最终输出数据为256位的数据,该数据包含4组像素的R、G、B通道的信息,在同步矩阵剥离运算中,每一组的运算包含9个硬件乘法器M和3个硬件除法器,乘法器的一个乘数为计算出的图像采样值,另一个乘数为逆矩阵计算出的逆矩阵T-1的一个光强参数。由于逆矩阵中的参数都是小数,若直接使用这些小数计算会加大运算难度,因此通过还原成整数的乘除法解决该问题,最终计算出4个像素各自的三个波段的光强参数。
c、FPGA芯片中的FIFO存储器可以配置为256位数据输入,16位数据输出的模式,如图9所示,低位数据首先输出,高位数据后输出,同步传输模块将16位数据从FIFO中读出,然后将该数据以符合USB传输时序的方式写入到数据传输芯片CY7C68013中,然后传输到计算机。
采用以上所述FPGA处理结构,对于500百万像素的彩色相机芯片,可以实现15帧/秒的处理和传输速度。
在样品位置分别放置中心波长470nm,560nm和650nm带宽10nm的干涉滤光片进行成像试验。未采用FPGA剥离处理时,三种波长的光在摄相机芯片RGB传感器上的响应为:
FPGA剥离后,三种波长的光在RGB三个通道上响应为:
B | G | R | |
470nm | 205 | 0 | 0 |
560nm | 0 | 231 | 0 |
650nm | 0 | 0 | 245 |
可以看出RGB三个通道上已分别对应470nm,560nm和650nm三个光谱波段的检测,对其他通道没有响应。
图10是采用伊红、亮绿、福尔根染料复染时,宫颈脱落细胞在显微镜下放大成像的图像,可以看出细胞核、细胞浆的福尔根、伊红、亮绿染色叠在一起无法进行细胞核内DNA(福尔根染料表示)的定量。
图11采用本发明的三波段实时光谱成像,从复染图像中剥离的细胞核福尔根染色图像,可以看出伊红、亮绿染料的影响已经去除,只剩下细胞核福尔根染料的图像。根据此图像,我们就可以测量细胞核内DNA的含量变化,判断细胞有无发生癌变。
Claims (7)
1.一种多波段光谱成像方法,在由光学成像器件和彩色摄像传感器构成的光路中,加入多波段带通滤光片,所述多波段带通滤光片置于所述光学成像器件和彩色摄像传感器的光路上,并且所述彩色摄像传感器进行光谱成像形成混叠的多波段多光谱图像,通过FPGA芯片将混叠的多波段多光谱图像,按预先建立的模型,实时还原成独立的多波段图像送入计算机;
所述FPGA芯片对多路并行处理的流程包括:a、根据彩色图像传感器的输出时序和输
出信号拼接为一个128位的数据,并根据行序列识别模块,判断行图像是G/R行还是
B/G行,进而判定将此128位数据是存入存储G/R行图像的FIFO中还是存入存储B/G
行图像的FIFO中;
b、在获取了一组G/R和B/G图像后,同步解码模块对图像数据进行插值运算,最终输出数据为256位的数据,该数据包含4组像素的R、G、B通道的信息,在同步矩阵剥离运算中,每一组的运算包含9个硬件乘法器M和3个硬件除法器,乘法器的一个乘数为计算出的图像采样值,另一个乘数为模型Model计算出的一个参数;
c、FPGA芯片中的FIFO存储器可以配置为256位数据输入,16位数据输出的模式,低位数据首先输出,高位数据后输出,同步传输模块将16位数据从FIFO中读出,然后将该数据以符合USB传输时序的方式写入到数据传输芯片CY7C68013中,然后传输到计算机。
2.如权利要求1所述的方法,其特征在于,所述光路还包括照明光源,所述照明光源置于所述光学成像器件的前端,并且在所述照明光源和所述光学成像器件之间装有样品。
3.如权利要求1所述的方法,其特征在于,所述FPGA芯片可采用DSP数字信号处理芯片代替。
4.如权利要求1所述的方法,其特征在于,FPGA芯片选用ALTERA公司的Cyclone IV E系列的EP4CE40F23C6。
5.如权利要求1所述的方法,其特征在于,所述FPGA芯片的内部包括图像数据获取与存储模块、同步解码模块、同步矩阵剥离运算模块、以及数据传输模块。
6.如权利要求1所述的方法,其特征在于,所述模型的建立,包括以下步骤:
步骤a、将定制的只有一个波段的窄带滤光片放入光路中,窄带曲线对应多波段带通滤光片的一个通道;
步骤b、对空白样品进行多光谱成像;测出彩色图像传感器的响应x1(x1r,x1g,x1b);重复步骤a、b,依次将定制的另一个波段的窄带滤光片放入光路中,测出另一个波段的窄带滤光片的响应xi(xir,xig,xib);
测量出所有波段响应xn(r,g,b),然后综合成敏感度矩阵X,X=[x1,x2,……xn];n≥3;
计算得到模型Model=(XtX)-1Xt。
7.如权利要求6所述的方法,其特征在于,以上所述每个窄带滤光片的滤光曲线只能覆盖多波段滤光片的一个通道。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510230474.XA CN104819941B (zh) | 2015-05-07 | 2015-05-07 | 一种多波段光谱成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510230474.XA CN104819941B (zh) | 2015-05-07 | 2015-05-07 | 一种多波段光谱成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104819941A CN104819941A (zh) | 2015-08-05 |
CN104819941B true CN104819941B (zh) | 2017-10-13 |
Family
ID=53730297
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510230474.XA Active CN104819941B (zh) | 2015-05-07 | 2015-05-07 | 一种多波段光谱成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104819941B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110530781A (zh) * | 2019-08-30 | 2019-12-03 | 武汉呵尔医疗科技发展有限公司 | 一种单帧三波段数字图像采集剥离系统 |
CN113125008A (zh) * | 2019-12-31 | 2021-07-16 | 北京微视新纪元科技有限公司 | 一种高光谱相机、光谱仪及其制作方法 |
CN111565305A (zh) * | 2020-05-12 | 2020-08-21 | Oppo(重庆)智能科技有限公司 | 一种图像传感器、信号处理方法及存储介质 |
CN114323278A (zh) * | 2021-12-23 | 2022-04-12 | 长光驰宇科技(长春)有限公司 | 一种灵巧型多光谱成像设备 |
CN118050328A (zh) * | 2024-03-26 | 2024-05-17 | 华南农业大学 | 一种单目高分多光谱的成像系统及控制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1914492A (zh) * | 2003-12-18 | 2007-02-14 | 阿瓦戈科技传感器Ip股份有限公司 | 利用混合滤光片的波长相关成像和探测方法及系统 |
CN103308466A (zh) * | 2013-06-04 | 2013-09-18 | 沈阳仪表科学研究院有限公司 | 便携式滤光片色轮型多光谱成像系统及其光谱图像处理方法 |
CN103954570A (zh) * | 2014-04-21 | 2014-07-30 | 江苏大学 | 一种基于光谱成像技术的粮食虫害程度判别方法 |
JP2014202673A (ja) * | 2013-04-08 | 2014-10-27 | キヤノン株式会社 | マルチバンド撮像装置およびその制御方法 |
JP2015072713A (ja) * | 2009-12-18 | 2015-04-16 | フラームセ・インステリング・フォール・テヒノロヒス・オンデルズーク・ナムローゼ・フェンノートシャップVlaamse Instelling Voor Technologisch Onderzoek N.V. | マルチスペクトルデータの幾何学的リファレンシング |
US9285309B2 (en) * | 2013-02-04 | 2016-03-15 | Hanwha Techwin Co., Ltd. | Image fusion method and apparatus using multi-spectral filter array sensor |
-
2015
- 2015-05-07 CN CN201510230474.XA patent/CN104819941B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1914492A (zh) * | 2003-12-18 | 2007-02-14 | 阿瓦戈科技传感器Ip股份有限公司 | 利用混合滤光片的波长相关成像和探测方法及系统 |
JP2015072713A (ja) * | 2009-12-18 | 2015-04-16 | フラームセ・インステリング・フォール・テヒノロヒス・オンデルズーク・ナムローゼ・フェンノートシャップVlaamse Instelling Voor Technologisch Onderzoek N.V. | マルチスペクトルデータの幾何学的リファレンシング |
US9285309B2 (en) * | 2013-02-04 | 2016-03-15 | Hanwha Techwin Co., Ltd. | Image fusion method and apparatus using multi-spectral filter array sensor |
JP2014202673A (ja) * | 2013-04-08 | 2014-10-27 | キヤノン株式会社 | マルチバンド撮像装置およびその制御方法 |
CN103308466A (zh) * | 2013-06-04 | 2013-09-18 | 沈阳仪表科学研究院有限公司 | 便携式滤光片色轮型多光谱成像系统及其光谱图像处理方法 |
CN103954570A (zh) * | 2014-04-21 | 2014-07-30 | 江苏大学 | 一种基于光谱成像技术的粮食虫害程度判别方法 |
Non-Patent Citations (3)
Title |
---|
Multispectral imaging using;Rui Shogenji et al;《opticas express》;20040419;第12卷(第8期);第1646-1650页第3.1-3.4、4.1节 * |
基于人造复眼的多谱成像系统;金建 等;《光学精密工程》;20141031;第22卷(第10期);第2691-2697页 * |
多波段滤光片辐亮度计的系统级定标方法研究;张艳娜 等;《光学学报》;20130831;第33卷(第8期);都0812004-1-0812004-7页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104819941A (zh) | 2015-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104819941B (zh) | 一种多波段光谱成像方法 | |
CN102279050B (zh) | 一种多光谱计算重构方法及系统 | |
CN101839979B (zh) | 作物冠层植被指数测量方法及装置 | |
CN108020519B (zh) | 一种基于颜色恒常性的虚拟多光源光谱重建方法 | |
CN106226247A (zh) | 一种基于高光谱显微成像技术的细胞检测方法 | |
CN111256826B (zh) | 显示屏色度测量方法、装置及终端设备 | |
EP3446083B1 (en) | Digital pathology color calibration | |
CN103389162B (zh) | 一种利用彩色数码相机的反射物体颜色测量仪 | |
IL141902A (en) | Spectral imaging to classify cells using internal reference | |
CN106679808B (zh) | 基于压缩光谱的关联成像系统及其方法 | |
CN111024229B (zh) | 单芯片集成型光谱成像微系统光谱数据校正方法 | |
AU700085B2 (en) | Multi-spectral segmentation for image analysis | |
Gao et al. | High-throughput fast full-color digital pathology based on Fourier ptychographic microscopy via color transfer | |
CN115100082A (zh) | 一种基于高光谱相机的高精度颜色显示系统 | |
CN106895916B (zh) | 一种单次曝光拍摄获取多光谱图像的方法 | |
US20220138986A1 (en) | Matching two color measurement devices using artificial neural network | |
CN112098415A (zh) | 一种杨梅品质无损检测方法 | |
JP4174707B2 (ja) | 分光測定システム、色再現システム | |
CN106124054A (zh) | 一种大幅面光谱分光成像测色装置 | |
Rice et al. | Hyperspectral image projector applications | |
WO2022057506A1 (zh) | 图像传感器、终端、数据处理方法、装置及存储介质 | |
CN211504402U (zh) | 一种基于渐变薄膜滤光片的宽波段高光谱相机 | |
Kim et al. | Radiometric characterization of spectral imaging for textual pigment identification | |
Sun et al. | Simulation of proposed eight-band camera for capturing multispectral images | |
CN111007020A (zh) | 双帧四光谱成像方法及应用 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |