CN107505046B - 多光谱颜色测量方法 - Google Patents

多光谱颜色测量方法 Download PDF

Info

Publication number
CN107505046B
CN107505046B CN201710637410.0A CN201710637410A CN107505046B CN 107505046 B CN107505046 B CN 107505046B CN 201710637410 A CN201710637410 A CN 201710637410A CN 107505046 B CN107505046 B CN 107505046B
Authority
CN
China
Prior art keywords
parameters
matrix
color
spectrum
image
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
CN201710637410.0A
Other languages
English (en)
Other versions
CN107505046A (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.)
Focusight Technology Co Ltd
Original Assignee
Focusight 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 Focusight Technology Co Ltd filed Critical Focusight Technology Co Ltd
Priority to CN201710637410.0A priority Critical patent/CN107505046B/zh
Publication of CN107505046A publication Critical patent/CN107505046A/zh
Application granted granted Critical
Publication of CN107505046B publication Critical patent/CN107505046B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J3/50Measurement of colour; Colour measuring devices, e.g. colorimeters using electric radiation detectors
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/46Measurement of colour; Colour measuring devices, e.g. colorimeters
    • G01J3/463Colour matching
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Color Image Communication Systems (AREA)

Abstract

本发明涉及一种多光谱颜色测量方法,包含颜色测量过程中参数的计算方法以及利用参数实施颜色测量的方法;首先计算出颜色测量中使用到的参数,主要包括三个部分:色差仪参数的计算,图像的预处理过程,颜色测量模型的建立;其次,应用计算得到参数,实施颜色测量。本发明解决了现有技术中数字相机无法对光谱值进行测量;而光谱仪虽然能够对光谱值进行测量,但是其测量幅面有限的问题;通过数字相机检测光谱反射率,相对于传统的方法,同时兼顾了两个优点,能够对光谱放射率进行测量,同时能过实现全幅面的检测。

Description

多光谱颜色测量方法
技术领域
本发明涉及印刷品表面颜色测量领域,尤其是一种通过4镜头的数字相机采集的待检测品的RGB图像,测量相应的光谱反射率和Lab颜色值的多光谱颜色测量方法。
背景技术
传统印刷产品表面的颜色质量测量方法有目视法、光电积分法和分光光度法。
目视法是利用人眼观察标准印刷品的颜色和样品印刷品的颜色差别,通常需要在国际色彩联盟规定的标准光源下进行。但是目视法主要借助检验工的颜色观察经验进行评判,受主观因素的影响过大,测量效率低下,不能形成统一的颜色质量评价标准。
光电积分法就是一种颜色测量仪器,通过探测器接收印刷品表面的光谱能量,和其自身的光谱响应进行积分测量。虽然光电积分法测量精度高、测量速度快,在印刷行业和纺织印染行业的颜色控制中得到广泛的应用,但是光电积分法无法获得印刷品表面的光谱组成,对于同色异谱的产品无法识别。
分光光度法是对印刷品表面的光谱反射特性进行测量,根据得到的光谱数据计算在CIE标准光源下对应的三刺激值。其测量方法是根据分光色散系统对印刷品表面的光谱进行机械扫描,测量各个波段对应的光谱能量,由此测量整个光谱功率分布,虽然这种方法测量精度高,但是测量速度慢、测量效率低。为了解决测量速度的问题提出用多通道技术平行探测全波段光谱能量分布的方法,利用快速存取和分组处理的方法实现了时间分辨和光谱分辨之间的兼顾。但是分光光度法只能实现印刷品的局部区域的颜色检测,通过在印刷品表面逐行扫描对印刷品进行检测,对于全幅面的检测无法达到。
发明内容
本发明要解决的技术问题是:提出一种多光谱颜色测量方法,通过数字相机的RGB值,建立RGB图像和光谱反射率之间的关系,并利用这个关系通过RGB图像计算出光谱值。
本发明所采用的技术方案为:一种多光谱颜色测量方法,包括以下步骤,
1)计算出颜色测量中使用到的参数;包括色差仪参数的计算、图像的预处理过程以及颜色测量模型的建立;
2)应用计算得到参数,实施颜色测量;
所述的步骤1)中,色差仪参数的计算包括光谱到XYZ的转换参数计算以及XYZ值到Lab的转换参数;图像的预处理过程包括白平衡参数的计算以及图像配准参数的计算;颜色测量模型的建立包括RGB到光谱的转换参数计算。
进一步的说,本发明所述的光谱到XYZ的转换参数计算的方式为:选取色卡图,对色卡图的不同颜色区进行测量,得到光谱数据的矩阵和XYZ数据的矩阵;采用线性回归法和最小二乘法求得光谱值与XYZ值之间的关系公式为:Q=(SPT*SP)-1*SPT*C;其中,SP为光谱数据的矩阵,C为XYZ数据的矩阵。
再进一步的说,本发明所述的白平衡参数的计算方式为:从采集到的图像中提取白平衡纸区域作为白平衡图像,再分别计算白平衡图像中每一列的三个通道平均值。
再进一步的说,本发明所述的RGB到光谱的转换参数计算方式为:提取每个色块中心区域的RGB平均值,记为矩阵D;通过最小二乘法求解得到RGB到光谱的转换参数矩阵为:
Figure BDA0001365157570000021
其中,SP为光谱数据的矩阵,
Figure BDA0001365157570000031
为矩阵D进行扩展后得到的另一个矩阵。
再进一步的说,本发明所述的步骤2)中,应用计算得到参数,计算出检测区域的光谱反射率和Lab值的方法是:
A、利用相机采集待检测品的图像,对采集到的图像进行预处理;
B、在预处理后的图像上选定检测区域,提取检测区域内的平均值,将所有通道的数据记为矩阵d;
C、对提取的通道数据进行扩展,记为另一个的矩阵dm
D、通过已计算得到RGB到光谱的转换参数P以及光谱到XYZ的转换参数Q,计算检测区域的光谱Spm以及XYZ值Cm如下式:
Spm=dm*P
Cm=Spm*Q
E、根据标准公式计算从XYZ计算得到Lab值。
本发明的有益效果是:解决了现有技术中数字相机无法对光谱值进行测量;而光谱仪虽然能够对光谱值进行测量,但是其测量幅面有限的问题;通过数字相机检测光谱反射率,相对于传统的方法,同时兼顾了两个优点,能够对光谱放射率进行测量,同时能过实现全幅面的检测。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明颜色测量参数计算流程图;
图2是本发明颜色测量流程图;
图3是本发明相机成像示意图;
图4是本发明包含4个镜头的图像;
图5是本发明白平衡纸区域图像;
图6是本发明提取红色矩形框内的RGB平均值图像。
具体实施方式
现在结合附图和优选实施例对本发明作进一步详细的说明。这些附图均为简化的示意图,仅以示意方式说明本发明的基本结构,因此其仅显示与本发明有关的构成。
一种多光谱颜色测量方法,包含颜色测量过程中参数的计算方法以及利用参数实施颜色测量的方法。
首先计算出颜色测量中使用到的参数,主要包括三个部分:色差仪参数的计算(包括:光谱到XYZ的转换参数计算,XYZ值到Lab的转换参数),图像的预处理过程(包括:白平衡参数的计算,图像配准参数的计算),颜色测量模型的建立(RGB到光谱的转换参数计算)。参数的计算流程如图1所示。
其次,应用计算得到参数,实施颜色测量,颜色测量流程如图2所示。
1、色差仪参数的计算
以X-rite SP62色差仪为例,该色差仪在400nm-700nm光谱范围内以10nm等间隔采样,对物体进行光谱反射率测量。因此印刷品表面的光谱反射率由31个光谱参数构成,即[S1 S2 … S31];同时可以得到XYZ值和L*a*b*值。
1)光谱到XYZ的转换参数计算
采用线性回归法求解光谱值与XYZ值之间的关系,设Q为31*3的参数矩阵。例如采用chromasens公司的色卡图,该色卡图包含345个不同颜色的色块区域(如图4),对色卡图的345个不同颜色区进行测量,得到光谱数据的矩阵如下:
Figure BDA0001365157570000051
其中,i表示第i个测量区域,j表示第j个光谱值。
同时,得到XYZ数据的矩阵如下:
Figure BDA0001365157570000052
其中,i表示第i个测量区域
由线性回归法,有如下关系公式:
SP*Q=C
用最小二乘法求解,如式:
Q=(SPT*SP)-1*SPT*C
2)XYZ值到Lab的转换参数
XYZ值到Lab值的转换采用标准公式,如下:
L=116*(f(Y/Yn))-16
a=500[f(X/Xn)-f(Y/Yn)]
b=200[f(Y/Yn)-f(Z/Zn)]
其中,
Figure BDA0001365157570000053
分别带入XYZ值到Lab值,求出转换参数Xn,Yn,Zn。
2、图像的预处理过程
本实施例基于多光谱成像系统,由多通道数字相机,镜头,滤光片和光源构成,如图3所示。12通道的线扫描CCD相机是Chromasens的TruePixa相机,包含四个镜头,每个镜头前面安装不同的滤光片(分别为:蓝色、普通透明、浅黄、绿色),每个镜头成像后都能得到一幅三通道的RGB图像,因此每种颜色会产生12个不同的刺激响应。包含四个镜头的图像图4所示。
1)白平衡参数的计算
本实施例采用的线扫描CCD相机,对其采集的图像(如图5)的每一列的每一个通道计算一个白平衡参数。
以纯正白色(白平衡纸)为基准,以光源的色温为尺度,通过自动或手动调整白平衡,对偏离真实的景物色彩进行补偿和改正,消除偏色,以贴近人的视觉标准。
在检品的顶端放置一条白平衡纸,相机采集到的图像如图4所示。从采集到的图像中提取白平衡纸区域作为白平衡图像,再分别计算白平衡图像中每一列的三个通道(R通道、G通道或B通道)平均值。设定白平衡的目标值(取值在180至220之间),则对应每一列的一个通道(R通道、G通道或B通道)白平衡参数为目标值除以对应列的平均值。
2)几何校正参数的计算
十二通道的CCD相机因为安装精度,视场的变化等因素的影响,所以四个镜头得到的四幅图像的位置会有一些差异。在处理多光谱图像时需要对四幅图像配准,即寻找四幅图像之间的对应关系,对四幅图像进行平移几何校正。
本实施例采用互相关度量法计算两个图像之间的互相关函数,根据互相关函数的峰值位置确定平移变换关系。
3)利用计算得到的白平衡参数和几何校正参数对原采集图像进行白平衡和图像
3、颜色测量模型的建立
颜色测量模型的建立,实际上建立RGB图像(4个镜头,共十二通道)和光谱反射率之间的关系。
1)RGB数据提取以及光谱数据采集
通过TruePixa相机采集色卡图像,得到4个镜头的RGB图像,对图像进行预处理后,提取图像每一个色块的中心区域的平均值,图6显示其中两个镜头的色块中心区域(红色矩形块所示)。提取每个色块中心区域的RGB平均值,共四个镜头,记为一个1*12的矩阵:
[R1 G1 B1 R2 G2 B2 R3 G3 B3 R4 G4 B4]
提取345个色块的中心区域的12通道RGB平均值,记为:
Figure BDA0001365157570000071
其中,i代表第几个色块
例如用X-rite SP62色差仪对345个色块进行测量,得到光谱数据记为:
Figure BDA0001365157570000072
其中,i表示第i个色块,j表示第j个光谱值。
2)RGB到光谱的转换参数计算
本实施例采用了一种二次多项式回归法建立RGB与光谱之间的关系,12通道的数据可以看作一个1*12的矩阵,对这个1*12的矩阵进行扩展,增加一个常数项1,同时增加每一个通道的数据值的平方项,得到一个维度为1*25的矩阵。按照这种方式对12*345矩阵D进行扩展得到一个25*345的矩阵,如下:
Figure BDA0001365157570000081
记P为25*31的RGB到光谱的转换参数矩阵,有式:
Figure BDA0001365157570000082
用最小二乘法求解,如式:
Figure BDA0001365157570000083
4、颜色测量的实现
颜色测量的实现方法如图2所示,通过truePIXA相机采集待检测品的图像,利用本方法计算出检测区域的光谱反射率和Lab值。
1)利用truePIXA相机采集待检测品的图像,对采集到的图像进行预处理
2)在预处理后的图像上选定检测区域,提取检测区域内的平均值,将12个通道的数据记为一个1*12的矩阵。
3)采用上文提到的方式对提取12通道数据进行扩展,记为一个1*25的矩阵dm
4)已计算得到RGB到光谱的转换参数P以及光谱到XYZ的转换参数Q,可计算检测区域的光谱Spm以及XYZ值Cm如下式:
Spm=dm*P
Cm=Spm*Q
根据标准公式计算从XYZ计算得到Lab值。
以上说明书中描述的只是本发明的具体实施方式,各种举例说明不对本发明的实质内容构成限制,所属技术领域的普通技术人员在阅读了说明书后可以对以前所述的具体实施方式做修改或变形,而不背离本发明的实质和范围。

Claims (2)

1.一种多光谱颜色测量方法,其特征在于:包括以下步骤,
1)计算出颜色测量中使用到的参数;包括色差仪参数的计算、图像的预处理过程以及颜色测量模型的建立;
2)应用计算得到参数,实施颜色测量;
所述的步骤1)中,色差仪参数的计算包括光谱到XYZ的转换参数计算以及XYZ值到Lab的转换参数;图像的预处理过程包括白平衡参数的计算以及图像配准参数的计算;颜色测量模型的建立包括RGB到光谱的转换参数计算;
所述的RGB到光谱的转换参数计算方式为:提取每个色块中心区域的RGB平均值,记为矩阵D;通过最小二乘法求解得到RGB到光谱的转换参数矩阵为:
Figure FDA0002455309720000011
其中,SP为光谱数据的矩阵,
Figure FDA0002455309720000012
为矩阵D进行扩展后得到的另一个矩阵;
所述的步骤2)中,应用计算得到参数,计算出检测区域的光谱反射率和Lab值的方法是:
A、利用相机采集待检测品的图像,对采集到的图像进行预处理;
B、在预处理后的图像上选定检测区域,提取检测区域内的平均值,将所有通道的数据记为矩阵d;
C、对提取的通道数据进行扩展,记为另一个的矩阵dm;
D、通过已计算得到RGB到光谱的转换参数P以及光谱到XYZ的转换参数Q,计算检测区域的光谱Spm以及XYZ值Cm如下式:
Spm=dm*P
Cm=Spm*Q
根据标准公式计算从XYZ计算得到Lab值;
所述的光谱到XYZ的转换参数计算的方式为:选取色卡图,对色卡图的不同颜色区进行测量,得到光谱数据的矩阵和XYZ数据的矩阵;采用线性回归法和最小二乘法求得光谱值与XYZ值之间的关系公式为:Q=(SPT*SP)-1*SPT*C;其中,SP为光谱数据的矩阵,C为XYZ数据的矩阵。
2.如权利要求1所述的一种多光谱颜色测量方法,其特征在于:所述的白平衡参数的计算方式为:从采集到的图像中提取白平衡纸区域作为白平衡图像,再分别计算白平衡图像中每一列的三个通道平均值。
CN201710637410.0A 2017-07-31 2017-07-31 多光谱颜色测量方法 Active CN107505046B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710637410.0A CN107505046B (zh) 2017-07-31 2017-07-31 多光谱颜色测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710637410.0A CN107505046B (zh) 2017-07-31 2017-07-31 多光谱颜色测量方法

Publications (2)

Publication Number Publication Date
CN107505046A CN107505046A (zh) 2017-12-22
CN107505046B true CN107505046B (zh) 2020-06-12

Family

ID=60689077

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710637410.0A Active CN107505046B (zh) 2017-07-31 2017-07-31 多光谱颜色测量方法

Country Status (1)

Country Link
CN (1) CN107505046B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108319993A (zh) * 2017-12-26 2018-07-24 钟建哲 一种布料加工监测管理系统
CN108896176B (zh) * 2018-05-14 2019-10-11 浙江大学 一种多光谱成像系统的空间一致性校正方法
CN108896177B (zh) * 2018-05-29 2021-12-31 中国科学院光电研究院 多通道复用的光谱图像获取方法
CN108645515A (zh) * 2018-06-14 2018-10-12 征图新视(江苏)科技有限公司 基于多光谱同色同谱的颜色测量系统
CN109886351B (zh) * 2019-03-04 2021-10-15 北京麦飞科技有限公司 高光谱数据和高分辨率图像融合方法
CN110276129B (zh) * 2019-06-21 2022-11-29 贵州大学 基于燃烧控制系统控制燃烧炉供给燃料当量比的方法、系统及装置
CN110487406A (zh) * 2019-08-26 2019-11-22 北京麦飞科技有限公司 基于多项式拟合的高光谱辐射定标方法及系统
CN111798442B (zh) * 2020-07-15 2023-11-28 河南科技大学 一种白度测量方法及白度测量中的参数标定方法
CN117730655B (zh) * 2024-02-20 2024-05-14 云南省农业科学院质量标准与检测技术研究所 水稻种子活力定量分析方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103575673A (zh) * 2013-09-24 2014-02-12 岑夏凤 可测量色彩情感的分光光度计
CN104849220A (zh) * 2015-06-09 2015-08-19 武汉大学 一种平面式文物光谱图像获取方法
CN106896069A (zh) * 2017-04-06 2017-06-27 武汉大学 一种基于彩色数码相机单幅rgb图像的光谱重建方法
CN106895916A (zh) * 2017-01-09 2017-06-27 浙江大学 一种单次曝光拍摄获取多光谱图像的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7034942B2 (en) * 2002-01-07 2006-04-25 Hewlett-Packard Development Company, L.P. Color measurement device and color measurement method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103575673A (zh) * 2013-09-24 2014-02-12 岑夏凤 可测量色彩情感的分光光度计
CN104849220A (zh) * 2015-06-09 2015-08-19 武汉大学 一种平面式文物光谱图像获取方法
CN106895916A (zh) * 2017-01-09 2017-06-27 浙江大学 一种单次曝光拍摄获取多光谱图像的方法
CN106896069A (zh) * 2017-04-06 2017-06-27 武汉大学 一种基于彩色数码相机单幅rgb图像的光谱重建方法

Also Published As

Publication number Publication date
CN107505046A (zh) 2017-12-22

Similar Documents

Publication Publication Date Title
CN107505046B (zh) 多光谱颜色测量方法
CN107144353B (zh) 一种基于数码相机的纺织品色差测量方法
Leon et al. Color measurement in L∗ a∗ b∗ units from RGB digital images
RU2251084C2 (ru) Способ подбора цвета посредством устройства формирования электронного отображения
CN108020519B (zh) 一种基于颜色恒常性的虚拟多光源光谱重建方法
CN111256826B (zh) 显示屏色度测量方法、装置及终端设备
CN103954634A (zh) 一种印刷品在线质量检测系统
CN110017904B (zh) 一种基于ccd相机的多光谱辐射测温方法
JP2008528341A (ja) 画像内カラー計測を用いたウェブ印刷機のカラー制御
KR100238960B1 (ko) 색감측정장치
US9194746B1 (en) System for measuring deviation of printed color at a selected location on a moving substrate from a target color
CN101294867B (zh) Ccd成像设备相对光谱响应特性的标定方法及系统
JP2000102042A (ja) 撮像装置の分光感度特性測定方法および撮像データ構成方法
CN204831751U (zh) 一种工业彩色相机颜色响应测试与校正装置
CN111579505A (zh) 一种基于数字图像处理的卷烟包装材料色差检测方法
CN110044485B (zh) 一种图像式织物颜色测量方法
CN112488997B (zh) 基于特征插值的古代绘画印刷品颜色复现检测和评价方法
Picollo et al. Application of hyper-spectral imaging technique for colorimetric analysis of paintings
CN105628206A (zh) 一种测定不同位置茶树叶片的颜色的方法
CN114910416A (zh) 一种高光谱颜色测量定标系统及定标方法
CN211504402U (zh) 一种基于渐变薄膜滤光片的宽波段高光谱相机
CN110726536B (zh) 一种彩色数码反射显微镜颜色校正方法
Gebejes et al. Color and image characterization of a three CCD seven band spectral camera
CN111007020A (zh) 双帧四光谱成像方法及应用
JP2002323376A (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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 213161 No. 258-6 Jinhua Road, West Taihu Lake, Wujin District, Changzhou City, Jiangsu Province

Applicant after: Zhengtu Xinshi (Jiangsu) Science and Technology Co., Ltd.

Address before: 213161 No. 258-6 Jinhua Road, West Taihu Lake, Wujin District, Changzhou City, Jiangsu Province

Applicant before: Sign new map (Jiangsu) Technology Co. Ltd.

CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Wang Yansong

Inventor after: Du Weidong

Inventor after: Wang Zheng

Inventor after: Zhang Chune

Inventor before: Wang Yansong

Inventor before: Du Weidong

Inventor before: Wang Zheng

GR01 Patent grant
GR01 Patent grant