CN114235153A - 一种基于数码相机的织物颜色数字化测量方法 - Google Patents

一种基于数码相机的织物颜色数字化测量方法 Download PDF

Info

Publication number
CN114235153A
CN114235153A CN202210169547.9A CN202210169547A CN114235153A CN 114235153 A CN114235153 A CN 114235153A CN 202210169547 A CN202210169547 A CN 202210169547A CN 114235153 A CN114235153 A CN 114235153A
Authority
CN
China
Prior art keywords
color
fabric
pure
measured
data
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.)
Granted
Application number
CN202210169547.9A
Other languages
English (en)
Other versions
CN114235153B (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.)
Wuhan Textile University
Original Assignee
Wuhan Textile University
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 Wuhan Textile University filed Critical Wuhan Textile University
Priority to CN202210169547.9A priority Critical patent/CN114235153B/zh
Publication of CN114235153A publication Critical patent/CN114235153A/zh
Application granted granted Critical
Publication of CN114235153B publication Critical patent/CN114235153B/zh
Priority to US17/869,771 priority patent/US11614362B1/en
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
    • 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/465Measurement of colour; Colour measuring devices, e.g. colorimeters taking into account the colour perception of the eye; using tristimulus detection
    • 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
    • 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0004Industrial image inspection
    • G06T7/001Industrial image inspection using an image reference approach
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/90Determination of colour characteristics
    • 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
    • G01J2003/466Coded colour; Recognition of predetermined colour; Determining proximity to predetermined colour
    • 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
    • G01J2003/467Colour computing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30108Industrial image inspection
    • G06T2207/30124Fabrics; Textile; Paper

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Quality & Reliability (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Image Processing (AREA)

Abstract

本发明为一种基于数码相机的织物颜色数字化测量方法,包括:制作纯色织物样本;利用分光光度计测量获得纯色织物真实颜色数据;利用数码相机拍摄获得纯色织物raw格式数字图像,提取纯色织物raw格式响应值;利用数码相机拍摄获得任意目标织物raw格式数字图像,提取目标织物待测量区域raw格式响应值;计算待测量区域与纯色织物raw响应值之间的欧氏距离和相似系数;对欧氏距离和相似系数进行归一化;计算每个纯色织物颜色数据的加权系数;利用加权系数对纯色织物颜色数据进行加权;对所有纯色织物的加权颜色数据求和,得到目标织物待测量区域的颜色值。本发明克服了分光光度计仅能进行固定孔径单点测量的局限性,适用于任意复杂图案织物。

Description

一种基于数码相机的织物颜色数字化测量方法
技术领域
本发明属于计算机数字图像处理技术领域,具体涉及一种基于数码相机的织物颜色数字化测量方法。
背景技术
颜色是织物(如色棉、纱线、面料等)的重要外观属性和功能属性。一方面,织物的颜色效果,影响消费者的喜好程度,进而影响服装的销售;另一方面,在专业应用领域,如迷彩服装的设计应用中,通过设计实现其对环境颜色的仿制,可发挥重要的隐蔽功能属性。颜色控制是纺织生产制造工艺中的关键技术,颜色测量是颜色控制的重要前提。但织物具有种类丰富性、结构形态多样性,给织物颜色的准确测量造成阻碍。
在织物颜色测量方面,传统方法主要通过人工对色——目视比色法,利用色卡进行织物颜色的定性测量和评估,但由于工作人员的年龄、性别和工作经验等主观差异性,不可避免导致颜色目视评估结果的不一致性,而且工作效率较低。随着分光光度计等专业测量设备的出现,分光光度法成为目前织物颜色测量的主流方法,分光光度法基于光学原理和色度学理论,能够准确测量织物颜色数据,但分光光度计仅能针对颜色单一的平面物体进行接触式颜色测量,无法直接面向形状不规则的物体、具有复杂图案或特殊结构的纺织面料、纤细的纱线、柔软的色棉等进行快速颜色测量,在实际纺织生产应用中仍然存在局限性,难以面向种类丰富性、结构形态多样性的织物应用。针对当前分光光度计在复杂织物颜色测量应用中的局限性,目前国内外的学术界和工业界均未提出很好的应对解决方法。
发明内容
本发明的目的是为了解决背景技术中所述问题,提出一种基于数码相机的织物颜色数字化测量方法。织物的颜色数据和其数码相机响应值具有高度关联性,因此基于数码相机成像技术对织物颜色进行测量理论可行,该技术以数字成像原理为基础,首先拍摄获得目标织物数字图像,然后依据成像模型反向计算织物的颜色数据,能够非接触地对任意形态结构和复杂图案织物对象进行客观的颜色测量,能够有效避免传统人工对色所存在的主观差异,避免目前价格昂贵的分光光度计在实际应用中的局限性,实现各类型织物颜色的快速、灵活测量与分析。
针对当前分光光度计在复杂织物颜色测量方面存在的实际应用问题,本发明基于数码相机成像技术提出了一种织物颜色测量新方法,首先收集制作纯色织物样本并测量获得其颜色数据,其次利用数码相机分别拍摄纯色织物和目标织物待测量区域数字图像,提取对应的响应值,然后计算目标织物待测量区域与纯色织物之间响应值欧氏距离和相似测度,得到纯色织物颜色数据的加权系数,最后利用加权系数对纯色织物颜色数据加权求和,得到目标织物待测量区域颜色数据。本发明的技术方案为一种基于数码相机的织物颜色数字化测量方法,具体包括以下步骤:
步骤1,制作纯色织物样本;
步骤2,利用分光光度计测量获得纯色织物的真实颜色数据;
步骤3,利用数码相机拍摄获取纯色织物raw格式数字图像,并提取纯色织物raw响应值;
步骤4,利用数码相机拍摄获取任意目标织物raw格式数字图像,并提取目标织物待测量区域raw响应值;
步骤5,计算待测量区域与所有纯色织物raw响应值之间的欧氏距离和相似系数;
步骤6,对欧氏距离和相似系数进行归一化处理;
步骤7,计算每个纯色织物颜色数据的加权系数;
步骤8,利用每个纯色织物的加权系数对其颜色数据进行加权;
步骤9,对所有纯色织物的加权颜色数据求和,得到目标织物待测量区域的颜色值。
进一步的,步骤1中,由于纯色织物样本的颜色分布直接影响目标织物颜色测量的准确性,通常制作的纯色织物样本的色域体积越大、且样本在颜色空间分布越均匀,其应用性能越好,具体制作方法如下。针对特定类型织物生产工艺:
首先,收集该类型纯色织物的历史颜色数据,并利用凸包(convex hull)算法提取织物颜色数据的外壳数据点Chull;其中,凸包算法为现有算法,不在此赘述,详见参考文献1:Graham R L. An efficient algorithm for determining the convex hull of afinite planar set[J]. Info. Pro. Lett., 1972, 1: 132-133;
其次,调研获得该类型织物生产工艺中的颜色波动误差ΔE
然后,利用发明人前期提出的色彩体系构建方法(详见参考文献2:公开号为CN113096072A的中国发明专利“一种色织纱线色彩体系构建方法”),将颜色数据的外壳数据Chull和颜色波动误差ΔE作为输入数据,生成纯色织物的理论颜色数据集;
最后,由纺织企业完成纯色织物实体样本的生产,并将纯色织物样固定于白卡纸板,最终获得纯色织物实体样本。
进一步的,步骤2中,利用分光光度计测量获得纯色织物真实颜色数据的方法具体如下:
首先,测量获得织物的CIEXYZ三刺激值数据。依据色度学理论,三刺激值的测量原理如式(1)至式(2)所示,
Figure 32894DEST_PATH_IMAGE001
, (1)
其中,
Figure 340248DEST_PATH_IMAGE002
, (2)
式中,x(λ) 、y(λ)和z(λ)为标准观察者颜色匹配函数,E(λ)为物体的光谱反射率,S(λ)为光源的相对光谱功率分布函数,λ为波长,k为调节因子,XYZ分别为织物的三刺激值数据。
然后,计算得到其对应的CIELab颜色数据。依据色度学理论,由三刺激值CIEXYZ数据计算得到其对应的CIELab颜色数据的方法如式(3)至式(4)所示,
Figure 214925DEST_PATH_IMAGE003
, (3)
其中,
Figure 797085DEST_PATH_IMAGE004
, (4)
式中,Lab分别为织物在CIELab颜色空间的亮度、红绿和黄蓝颜色值,XYZ分别为织物的三刺激值颜色数据,X n Y n Z n 分别为参考光源的三刺激值颜色数据,式(4)中HH n 分别表示织物和参考光源的CIEXYZ三刺激值。
完成纯色织物的真实颜色数据测量,得到每个纯色织物样本的颜色数据c i ,c i 为1×3的行向量,其形式如式(5)所示,
Figure 549140DEST_PATH_IMAGE005
, (5)
其中,i指示第i个纯色织物样本,i的取值为1到NN为纯色织物样本的数量,L i a i b i 分别为第i个纯色织物样本在CIELab颜色空间的亮度、红绿和黄蓝颜色数据。
进一步的,步骤3中利用数码相机拍摄获取纯色织物raw格式数字图像,并提取纯色织物raw响应值的方法如下 :
首先,在均匀照明的光照环境下,利用数码相机拍摄纯色织物的raw格式数字图像;
然后,对于每个纯色织物样本数字图像,提取其中心m×m像素区域内的所有像素的raw响应值数据,并对m×m像素点raw响应值数据求平均,得到纯色织物的raw响应值数据,如式(6)所示:
Figure 57089DEST_PATH_IMAGE006
, (6)
其中,i指示第i个纯色织物样本,j指示提取区域内第j个像素,r i,j g i,j b i,j 分别为第i个纯色样本的第j个像素红、绿、蓝三通道raw响应值,d i 为第i个纯色织物样本的raw响应值,为1×3的行向量。
进一步的,步骤4中利用数码相机拍摄获取任意目标织物raw格式数字图像,并提取目标织物待测量区域raw响应值的方法如下:
首先,在与拍摄纯色织物的光照和成像参数完全一致的条件下,利用数码相机获取任意目标织物raw格式数字图像;
然后,选取待测量区域,待测量区域如果是单个像素,则直接存储为d t ,d t 为待测量区域的raw响应值行向量,维度为1×3,待测量区域如果包含L个像素,则按照式(6)求取L个像素raw响应值的平均值,并存储为d t ,d t 的具体形式如式(7)所示。
Figure 527384DEST_PATH_IMAGE007
, (7)
其中,r t g t b t 分别为待测量区域红、绿、蓝三通道raw响应值。
进一步的,步骤5中计算待测量区域与所有纯色织物raw响应值之间的欧氏距离和相似系数的方法分别如下:
首先,计算待测量区域与所有纯色织物raw格式响应值之间的欧氏距离,并存储为行向量,如式(8)和式(9)所示,
Figure 647656DEST_PATH_IMAGE008
, (8)
Figure 254218DEST_PATH_IMAGE009
, (9)
其中,r t g t b t 分别为待测量区域红、绿、蓝三通道raw格式响应值,r i g i b i 分别为第i个纯色织物样本红、绿、蓝三通道raw格式响应值,e i 为待测量区域与第i个纯色织物raw格式响应值之间的欧氏距离,e为包含N个欧氏距离的行向量;
然后,计算待测量区域与所有纯色织物raw响应值之间的相似系数,并存储为行向量,如式(10)和式(11)所示,
Figure 935997DEST_PATH_IMAGE010
, (10)
Figure 893589DEST_PATH_IMAGE011
, (11)
其中,d t 为待测量区域raw响应值向量,d i 为第i个纯色织物样本raw响应值向量,‘·’表示向量内积,s i 为待测量区域与第i个纯色织物样本raw格式响应值之间的相似系数,s为包含N个相似系数的行向量。
进一步的,步骤6中对欧氏距离和相似系数进行归一化的方法分别如式(12)和式(13)所示:
Figure 817551DEST_PATH_IMAGE012
, (12)
Figure 278620DEST_PATH_IMAGE013
, (13)
其中,e n,i s n,i 分别为归一化的欧氏距离和相似系数,max(·)为求最大值函数,‘·/’表示相除。
进一步的,步骤7中计算每个纯色织物颜色数据加权系数的方法具体如下:
Figure 974043DEST_PATH_IMAGE014
, (14)
其中,w i 表示第i个纯色织物颜色数据的加权系数,exp(·)为指数函数。
进一步的,步骤8中利用每个纯色织物加权系数对其颜色数据进行加权的具体方法如下:
Figure 901155DEST_PATH_IMAGE015
, (15)
其中,c w,i 为第i个纯色织物的加权颜色数据,w i 表示第i个纯色织物颜色数据的加权系数,c i 为第i个纯色织物的颜色数据,L i a i b i 分别为第i个纯色织物在CIELab颜色空间的亮度、红绿和黄蓝颜色数据。
进一步的,步骤9中对所有纯色织物的加权颜色数据求和,得到目标织物待测量区域的颜色值的方法如式(16)所示:
Figure 379540DEST_PATH_IMAGE016
, (16)
其中,c t 为目标织物待测量区域的颜色值,c w,i 为第i个纯色织物的加权颜色数据,N表示纯色织物的个数,至此完成目标织物待测量区域的颜色测量。
本发明针对当前分光光度计在复杂织物颜色测量方面存在的局限性,依据织物颜色数据和其数码相机响应值之间的高度关联性,提出了基于数码相机的织物颜色数字化测量方法。能够面向形状不规则的物体、具有复杂图案或特殊结构的纺织面料、纤细的纱线、柔软的色棉等进行快速非接触式颜色测量,面向种类丰富性、结构形态多样性的织物进行应用。本发明克服了传统人工对色所存在的主观差异,解决了分光光度计在实际应用中的局限性。
附图说明
图1为本发明实施例的流程图。
图2为本发明实施例20个测量区域的颜色数据的分布。
具体实施方式
本发明技术方案具体实施时可由本领域技术人员采用计算机软件技术运行。
结合附图1,本发明实施例通过一种基于数码相机的织物颜色数字化测量方法,具体包括以下步骤:
步骤1,制作纯色织物样本;
步骤2,利用分光光度计测量获得纯色织物的真实颜色数据;
步骤3,利用数码相机拍摄获取纯色织物raw格式数字图像,并提取纯色织物raw响应值;
步骤4,利用数码相机拍摄获取任意目标织物raw格式数字图像,并提取目标织物待测量区域raw响应值;
步骤5,计算待测量区域与所有纯色织物raw响应值之间的欧氏距离和相似系数;
步骤6,对欧氏距离和相似系数进行归一化处理;
步骤7,计算每个纯色织物颜色数据的加权系数;
步骤8,利用每个纯色织物的加权系数对其颜色数据进行加权;
步骤9,对所有纯色织物的加权颜色数据求和,得到目标织物待测量区域的颜色值。
下面通过实施例具体说明每个步骤的处理过程:实施例以全棉平纹色织工艺、X-rite Color i7分光光度计、Canon 600D数码相机、以及VeriVide DigiEye灯箱为基础,对本发明方法进行测试。
步骤1中,实施例以A公司的纯色色织织物的1872个历史颜色数据为基础,利用凸壳算法提取外壳数据点共237个,该公司纯色色织生产工艺的颜色波动误差ΔE的值为3.5。按照参考文献2所示方法,共计算得到1236个纯色织物理论数据,通过打样工艺制作得到纯色织物样本1236个,将纯色织物样固定于白卡纸板,获得纯色织物实体样本。并按照相同生产工艺制作提花织物样本20个,作为待测试的目标织物。
步骤2中,纯色织物样本的真实颜色数据,采用X-rite Color i7分光光度计测量获得。实际测量时,首先测量获得纯色织物的三刺激值,具体测量原理如式(1)和式(2)所示:
Figure 616487DEST_PATH_IMAGE017
, (1)
其中,
Figure 171227DEST_PATH_IMAGE002
, (2)
式中,x(λ) 、y(λ)和z(λ)为标准观察者颜色匹配函数,E(λ)为物体的光谱反射率,S(λ)为光源的相对光谱功率分布函数,λ为波长,k为调节因子,XYZ分别为织物的三刺激值数据。获得织物的三刺激值之后,进一步利用式(3)和式(4)所示方法计算得到其对于的CIELab颜色数据,完成纯色织物的颜色测量。实施例中,计算织物的CIELab颜色数据时,采用纺织企业常用的CIED65日光光源作为参考光源。
Figure 41094DEST_PATH_IMAGE003
, (3)
其中,
Figure 572438DEST_PATH_IMAGE004
, (4)
式中,Lab分别为织物在CIELab颜色空间的亮度、红绿和黄蓝颜色值,XYZ分别为织物的三刺激值颜色数据,X n Y n Z n 分别为参考光源的三刺激值颜色数据,式(4)中HH n 分别表示织物和参考光源的CIEXYZ三刺激值。至此,完成纯色织物的真实颜色数据测量,得到每个纯色织物样本的颜色数据c i ,c i 为1×3的行向量,其形式如式(5)所示,
Figure 8099DEST_PATH_IMAGE005
, (5)
其中,i指示第i个纯色织物样本,i的取值为1到NN为纯色织物样本的数量,L i a i b i 分别为第i个纯色织物样本在CIELab颜色空间的亮度、红绿和黄蓝颜色数据。实施例中,N的值为1236.
步骤3中,完成纯色织物的真实颜色数据测量后,在VeriVide DigiEye均匀照明灯箱内,依次拍摄纯色织物样本,并存储每个织物样本的raw格式数字图像(*.cr2格式文件)。对于每个纯色织物样本数字图像,提取其中心m×m像素区域内的所有像素的raw响应值数据,并对m×m像素点raw响应值数据求平均,得到纯色织物的raw响应值数据,如式(6)所示:
Figure 465232DEST_PATH_IMAGE006
, (6)
其中,i指示第i个纯色织物样本,j指示提取区域内第j个像素,r i,j g i,j b i,j 分别为第i个纯色样本的第j个像素红、绿、蓝三通道raw响应值,d i 为第i个纯色织物样本的raw响应值,为1×3的行向量。实施例中,m的值取100。
步骤4中,完成纯色织物raw格式图像拍摄和raw响应值提取之后,在与拍摄纯色织物的光照和成像参数完全一致的条件下,利用数码相机获取目标提花织物raw格式数字图像,并选取待测量区域,提取待测量区域raw响应值。其中,待测量区域如果是单个像素,则直接存储为d t ,d t 为待测量区域的raw响应值行向量,维度为1×3,待测量区域如果包含L个像素,则按照式(6)求取L个像素raw响应值的平均值,并存储为d t ,d t 的具体形式如式(7)所示。
Figure 353553DEST_PATH_IMAGE007
, (7)
其中,r t g t b t 分别为待测量区域红、绿、蓝三通道raw响应值。实施例中,针对每个提花织物分别选择一个待测量区域,共20个待测量区域。其中,10个待测量区域为单个像素,10个待测量区域为不同大小的矩形区域,每个矩形完全位于提花织物的纯色区域内。
步骤5中,完成待测区域raw响应值提取之后,计算待测量区域与所有纯色织物raw响应值之间的欧氏距离和相似系数,方法分别如下:
首先,计算待测量区域与所有纯色织物raw格式响应值之间的欧氏距离,并存储为行向量,如式(8)和式(9)所示,
Figure 501638DEST_PATH_IMAGE008
, (8)
Figure 306652DEST_PATH_IMAGE009
, (9)
其中,r t g t b t 分别为待测量区域红、绿、蓝三通道raw格式响应值,r i g i b i 分别为第i个纯色织物样本红、绿、蓝三通道raw格式响应值,e i 为待测量区域与第i个纯色织物raw格式响应值之间的欧氏距离,e为包含N个欧氏距离的行向量;
然后,计算待测量区域与所有纯色织物raw响应值之间的相似系数,并存储为行向量,如式(10)和式(11)所示,
Figure 186883DEST_PATH_IMAGE010
, (10)
Figure 313233DEST_PATH_IMAGE011
, (11)
其中,d t 为待测量区域raw响应值向量,d i 为第i个纯色织物样本raw响应值向量,‘·’表示向量内积,s i 为待测量区域与第i个纯色织物样本raw格式响应值之间的相似系数,s为包含N个相似系数的行向量。
步骤6中,完成待测量区域与所有纯色织物raw响应值之间的欧氏距离和相似系数计算之后,对欧氏距离和相似系数进行归一化处理,方法分别如式(12)和式(13)所示:
Figure 671533DEST_PATH_IMAGE012
, (12)
Figure 331054DEST_PATH_IMAGE013
, (13)
其中,e n,i s n,i 分别为归一化的欧氏距离和相似系数,max(·)为求最大值函数,‘·/’表示相除。
步骤7中,完成欧氏距离和相似系数进行归一化之后,计算每个纯色织物颜色数据加权系数,方法如式(14)所示:
Figure 382186DEST_PATH_IMAGE014
, (14)
其中,w i 表示第i个纯色织物颜色数据的加权系数,exp(·)为指数函数。
步骤8中,完成每个纯色织物颜色数据加权系数计算之后,利用每个纯色织物加权系数对其颜色数据进行加权处理,方法如式(15)所示:
Figure 992903DEST_PATH_IMAGE015
, (15)
其中,c w,i 为第i个纯色织物的加权颜色数据,w i 表示第i个纯色织物颜色数据的加权系数,c i 为第i个纯色织物的颜色数据,L i a i b i 分别为第i个纯色织物在CIELab颜色空间的亮度、红绿和黄蓝颜色数据。
步骤9中,获得每个纯色织物的加权颜色数据之后,最后按照式(16)所示,对所有纯色织物的加权颜色数据求和,得到目标织物待测量区域的颜色值,至此完成目标织物的测量区域颜色测量。
Figure 154894DEST_PATH_IMAGE016
, (16)
其中,c t 为目标织物待测量区域的颜色值,c w,i 为第i个纯色织物的加权颜色数据,N表示纯色织物的个数。实施例中,按照上述步骤完成20个待测量区域的颜色测量,对应的颜色分布如附图2所示。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (10)

1.一种基于数码相机的织物颜色数字化测量方法,其特征在于,包括以下步骤:
步骤1,制作纯色织物样本;
步骤2,利用分光光度计测量获得纯色织物的真实颜色数据;
步骤3,利用数码相机拍摄获取纯色织物raw格式数字图像,并提取纯色织物raw响应值;
步骤4,利用数码相机拍摄获取任意目标织物raw格式数字图像,并提取目标织物待测量区域raw响应值;
步骤5,计算待测量区域与所有纯色织物raw响应值之间的欧氏距离和相似系数;
步骤6,对欧氏距离和相似系数进行归一化处理;
步骤7,计算每个纯色织物颜色数据的加权系数;
步骤8,利用每个纯色织物的加权系数对其颜色数据进行加权;
步骤9,对所有纯色织物的加权颜色数据求和,得到目标织物待测量区域的颜色值。
2.如权利要求1所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤1中,具体制作方法如下;
首先,收集某一类型纯色织物的历史颜色数据,并利用凸包算法提取织物颜色数据的外壳数据点Chull
其次,调研获得该类型织物生产工艺中的颜色波动误差ΔE
然后,利用色彩体系构建方法,将颜色数据的外壳数据Chull和颜色波动误差ΔE作为输入数据,生成纯色织物的理论颜色数据集;
最后,由纺织企业完成纯色织物实体样本的生产,并将纯色织物样固定于白卡纸板,最终获得纯色织物实体样本。
3.如权利要求1所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤2中,利用分光光度计测量获得纯色织物真实颜色数据的方法具体如下:
首先,测量获得纯色织物的CIEXYZ三刺激值数据,依据色度学理论,三刺激值的测量原理如式(1)至式(2)所示,
Figure 433432DEST_PATH_IMAGE001
, (1)
其中,
Figure 107996DEST_PATH_IMAGE002
, (2)
式中,x(λ) 、y(λ)和z(λ)为标准观察者颜色匹配函数,E(λ)为物体的光谱反射率,S(λ)为光源的相对光谱功率分布函数,λ为波长,k为调节因子,XYZ分别为纯色织物的三刺激值数据;
然后,计算得到其对应的CIELab颜色数据;依据色度学理论,由三刺激值CIEXYZ数据计算得到其对应的CIELab颜色数据的方法如式(3)至式(4)所示,
Figure 723785DEST_PATH_IMAGE003
, (3)
其中,
Figure 437269DEST_PATH_IMAGE004
, (4)
式中,Lab分别为纯色织物在CIELab颜色空间的亮度、红绿和黄蓝颜色值,XYZ分别为纯色织物的三刺激值颜色数据,X n Y n Z n 分别为参考光源的三刺激值颜色数据,式(4)中,HH n 分别表示纯色织物和参考光源的CIEXYZ三刺激值;
完成纯色织物的真实颜色数据测量,得到每个纯色织物样本的颜色数据c i ,c i 为1×3的行向量,其形式如式(5)所示,
Figure 556535DEST_PATH_IMAGE005
, (5)
其中,i指示第i个纯色织物样本,i的取值为1到NN为纯色织物样本的数量,L i a i b i 分别为第i个纯色织物样本在CIELab颜色空间的亮度、红绿和黄蓝颜色数据。
4.如权利要求1所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤3中利用数码相机拍摄获取纯色织物raw格式数字图像,并提取纯色织物raw响应值的方法如下:
首先,在均匀照明的光照环境下,利用数码相机拍摄纯色织物的raw格式数字图像;
然后,对于每个纯色织物样本数字图像,提取其中心m×m像素区域内的所有像素的raw响应值数据,并对m×m像素点raw响应值数据求平均,得到纯色织物的raw响应值数据,如式(6)所示:
Figure 933159DEST_PATH_IMAGE006
, (6)
其中,i指示第i个纯色织物样本,j指示提取区域内第j个像素,r i,j g i,j b i,j 分别为第i个纯色样本的第j个像素红、绿、蓝三通道raw响应值,d i 为第i个纯色织物样本的raw响应值,为1×3的行向量。
5.如权利要求1所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤4中利用数码相机拍摄获取任意目标织物raw格式数字图像,并提取目标织物待测量区域raw响应值的方法如下:
首先,在与拍摄纯色织物的光照和成像参数完全一致的条件下,利用数码相机获取任意目标织物raw格式数字图像;
然后,选取待测量区域,待测量区域如果是单个像素,则直接存储为d t ,d t 为待测量区域的raw响应值行向量,维度为1×3,待测量区域如果包含L个像素,则按照式(6)求取L个像素raw响应值的平均值,并存储为d t ,d t 的具体形式如式(7)所示;
Figure 770665DEST_PATH_IMAGE007
, (7)
其中,r t g t b t 分别为待测量区域红、绿、蓝三通道raw响应值。
6.如权利要求1所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤5中计算待测量区域与所有纯色织物raw响应值之间的欧氏距离和相似系数的方法分别如下:
首先,计算待测量区域与所有纯色织物raw格式响应值之间的欧氏距离,并存储为行向量,如式(8)和式(9)所示,
Figure 25191DEST_PATH_IMAGE008
, (8)
Figure 998963DEST_PATH_IMAGE009
, (9)
其中,r t g t b t 分别为待测量区域红、绿、蓝三通道raw格式响应值,r i g i b i 分别为第i个纯色织物样本红、绿、蓝三通道raw格式响应值,e i 为待测量区域与第i个纯色织物raw格式响应值之间的欧氏距离,e为包含N个欧氏距离的行向量;
然后,计算待测量区域与所有纯色织物raw响应值之间的相似系数,并存储为行向量,如式(10)和式(11)所示,
Figure 546488DEST_PATH_IMAGE010
, (10)
Figure 871290DEST_PATH_IMAGE011
, (11)
其中,d t 为待测量区域raw响应值向量,d i 为第i个纯色织物样本raw响应值向量,‘·’表示向量内积,s i 为待测量区域与第i个纯色织物样本raw格式响应值之间的相似系数,s为包含N个相似系数的行向量。
7.如权利要求6所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤6中对欧氏距离和相似系数进行归一化的方法分别如式(12)和式(13)所示:
Figure 926577DEST_PATH_IMAGE012
, (12)
Figure 551594DEST_PATH_IMAGE013
, (13)
其中,e n,i s n,i 分别为归一化的欧氏距离和相似系数,max(·)为求最大值函数,‘·/’表示相除。
8.如权利要求7所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤7中计算每个纯色织物颜色数据加权系数的方法具体如下:
Figure 535599DEST_PATH_IMAGE014
, (14)
其中,w i 表示第i个纯色织物颜色数据的加权系数,exp(·)为指数函数。
9.如权利要求1所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤8中利用每个纯色织物加权系数对其颜色数据进行加权的具体方法如下:
Figure 816539DEST_PATH_IMAGE015
, (15)
其中,c w,i 为第i个纯色织物的加权颜色数据,w i 表示第i个纯色织物颜色数据的加权系数,c i 为第i个纯色织物的颜色数据,L i a i b i 分别为第i个纯色织物在CIELab颜色空间的亮度、红绿和黄蓝颜色数据。
10.如权利要求1所述的一种基于数码相机的织物颜色数字化测量方法,其特征在于:步骤9中对所有纯色织物的加权颜色数据求和,得到目标织物待测量区域的颜色值的方法如式(16)所示:
Figure 84971DEST_PATH_IMAGE016
, (16)
其中,c t 为目标织物待测量区域的颜色值,c w,i 为第i个纯色织物的加权颜色数据,N表示纯色织物的个数,至此完成目标织物待测量区域的颜色测量。
CN202210169547.9A 2022-02-24 2022-02-24 一种基于数码相机的织物颜色数字化测量方法 Active CN114235153B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202210169547.9A CN114235153B (zh) 2022-02-24 2022-02-24 一种基于数码相机的织物颜色数字化测量方法
US17/869,771 US11614362B1 (en) 2022-02-24 2022-07-20 Method of digital measuring color of fabrics based on digital camera

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210169547.9A CN114235153B (zh) 2022-02-24 2022-02-24 一种基于数码相机的织物颜色数字化测量方法

Publications (2)

Publication Number Publication Date
CN114235153A true CN114235153A (zh) 2022-03-25
CN114235153B CN114235153B (zh) 2022-05-13

Family

ID=80747957

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210169547.9A Active CN114235153B (zh) 2022-02-24 2022-02-24 一种基于数码相机的织物颜色数字化测量方法

Country Status (2)

Country Link
US (1) US11614362B1 (zh)
CN (1) CN114235153B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115855809A (zh) * 2023-02-21 2023-03-28 武汉纺织大学 一种基于光谱重建的纺织品色牢度自动评级方法
CN115937207A (zh) * 2023-01-09 2023-04-07 江苏恒力化纤股份有限公司 一种图像式染色涤纶织物颜色测量方法
CN117073842A (zh) * 2023-07-21 2023-11-17 武汉纺织大学 一种基于纹理特征加权校正的纺织面料照相测色方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102714687A (zh) * 2010-01-19 2012-10-03 阿克佐诺贝尔国际涂料股份有限公司 从图像确定颜色的方法和系统
CN103389162A (zh) * 2012-11-26 2013-11-13 中国计量学院 一种利用彩色数码相机的反射物体颜色测量仪
CN106896069A (zh) * 2017-04-06 2017-06-27 武汉大学 一种基于彩色数码相机单幅rgb图像的光谱重建方法
CN110044485A (zh) * 2019-04-17 2019-07-23 江阴芗菲服饰有限公司 一种图像式织物颜色测量方法
CN111750994A (zh) * 2020-06-04 2020-10-09 武汉纺织大学 一种基于数码相机成像模型的光谱测量方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113096072B (zh) 2021-03-18 2022-08-30 武汉纺织大学 一种色织纱线色彩体系构建方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102714687A (zh) * 2010-01-19 2012-10-03 阿克佐诺贝尔国际涂料股份有限公司 从图像确定颜色的方法和系统
CN103389162A (zh) * 2012-11-26 2013-11-13 中国计量学院 一种利用彩色数码相机的反射物体颜色测量仪
CN106896069A (zh) * 2017-04-06 2017-06-27 武汉大学 一种基于彩色数码相机单幅rgb图像的光谱重建方法
CN110044485A (zh) * 2019-04-17 2019-07-23 江阴芗菲服饰有限公司 一种图像式织物颜色测量方法
CN111750994A (zh) * 2020-06-04 2020-10-09 武汉纺织大学 一种基于数码相机成像模型的光谱测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
万晓霞,等: "非接触式颜色测量参数的最优化研究", 《光谱学与光谱分析》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115937207A (zh) * 2023-01-09 2023-04-07 江苏恒力化纤股份有限公司 一种图像式染色涤纶织物颜色测量方法
CN115855809A (zh) * 2023-02-21 2023-03-28 武汉纺织大学 一种基于光谱重建的纺织品色牢度自动评级方法
CN117073842A (zh) * 2023-07-21 2023-11-17 武汉纺织大学 一种基于纹理特征加权校正的纺织面料照相测色方法和系统
CN117073842B (zh) * 2023-07-21 2024-03-29 武汉纺织大学 基于纹理特征加权校正的纺织面料照相测色方法和系统

Also Published As

Publication number Publication date
US11614362B1 (en) 2023-03-28
CN114235153B (zh) 2022-05-13

Similar Documents

Publication Publication Date Title
CN114235153B (zh) 一种基于数码相机的织物颜色数字化测量方法
US8743137B2 (en) Method for electronic color matching
CN105787508B (zh) 纺织品颜色识别方法及系统
CN109671113B (zh) 一种基于图像技术的色纺面料的测配色方法
CN109141640A (zh) 基于机器视觉的色差检测方法、系统、设备及存储介质
WO2003029811A1 (en) Assessing colour fastiness
JP3809838B2 (ja) 照明源固有hsv色座標を用いた画像ハイライト補正方法、画像ハイライト補正プログラム、および画像取得システム
CN115855809B (zh) 一种基于光谱重建的纺织品色牢度自动评级方法
CN101100807A (zh) 纺织品色牢度色差评级方法
CN110644216B (zh) 一种纺织品自动化色牢度评级装置
CN109582994B (zh) 一种色纺纱智能选色配色方法
CN106504241A (zh) 一种自动对色装置和方法
Gökmen et al. A non-contact computer vision based analysis of color in foods
CN110108362A (zh) 基于slic超像素分割的色差自适应在线检测方法及装置
CN111476153A (zh) 一种水果成熟度计算方法
CN114674438A (zh) 用于表征效果颜料斑点的装置和方法
CN104537639A (zh) 基于平移评价优选的光谱图像数据融合方法
Shrestha et al. Evaluation and comparison of multispectral imaging systems
JP2014187558A (ja) 画像色分布検査装置および画像色分布検査方法
CN112488997B (zh) 基于特征插值的古代绘画印刷品颜色复现检测和评价方法
KR101881693B1 (ko) 패션상품 디자인에 있어서 소재 직물 선택과 조명 조건하에서 패션상품의 비주얼 검토 장치 및 방법
MATUSİAK et al. Comparison of spectrophotometric and DigiEye colour measurements of woven fabrics
CN115880297A (zh) 基于机器视觉的被套染色质量评估方法
CN115753644B (zh) 一种基于非线性约束插值的织物颜色照相测量方法
CN108873955B (zh) 一种原液着色纺织材料的配色方法

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