CN107643267B - 一种基于可见光谱成像的古代壁画颜料无损全面识别方法 - Google Patents
一种基于可见光谱成像的古代壁画颜料无损全面识别方法 Download PDFInfo
- Publication number
- CN107643267B CN107643267B CN201710781294.XA CN201710781294A CN107643267B CN 107643267 B CN107643267 B CN 107643267B CN 201710781294 A CN201710781294 A CN 201710781294A CN 107643267 B CN107643267 B CN 107643267B
- Authority
- CN
- China
- Prior art keywords
- pigment
- sample
- mural
- ancient
- spectral
- 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.)
- Expired - Fee Related
Links
- 239000000049 pigment Substances 0.000 title claims abstract description 243
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000001429 visible spectrum Methods 0.000 title claims abstract description 51
- 238000003384 imaging method Methods 0.000 title claims abstract description 20
- 230000003595 spectral effect Effects 0.000 claims abstract description 116
- 238000002310 reflectometry Methods 0.000 claims abstract description 57
- 239000002245 particle Substances 0.000 claims abstract description 28
- 238000010422 painting Methods 0.000 claims abstract description 17
- 238000001228 spectrum Methods 0.000 claims abstract description 15
- 238000009826 distribution Methods 0.000 claims abstract description 13
- 230000011218 segmentation Effects 0.000 claims abstract description 12
- 239000000126 substance Substances 0.000 claims abstract description 11
- 238000010521 absorption reaction Methods 0.000 claims description 40
- 239000011159 matrix material Substances 0.000 claims description 31
- 239000003973 paint Substances 0.000 claims description 30
- 238000012549 training Methods 0.000 claims description 28
- 230000004044 response Effects 0.000 claims description 24
- 239000012463 white pigment Substances 0.000 claims description 18
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000010606 normalization Methods 0.000 claims description 14
- 239000000203 mixture Substances 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 11
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000009877 rendering Methods 0.000 claims description 7
- 238000005192 partition Methods 0.000 claims description 5
- 230000009466 transformation Effects 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 claims description 3
- 238000002156 mixing Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 abstract description 8
- 239000007787 solid Substances 0.000 abstract 1
- 239000000523 sample Substances 0.000 description 78
- 238000004458 analytical method Methods 0.000 description 15
- 239000011499 joint compound Substances 0.000 description 9
- 239000010410 layer Substances 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 239000010440 gypsum Substances 0.000 description 7
- 229910052602 gypsum Inorganic materials 0.000 description 7
- 239000013074 reference sample Substances 0.000 description 7
- 238000011160 research Methods 0.000 description 6
- 238000005070 sampling Methods 0.000 description 6
- 235000008733 Citrus aurantifolia Nutrition 0.000 description 4
- 235000011941 Tilia x europaea Nutrition 0.000 description 4
- 239000006229 carbon black Substances 0.000 description 4
- 239000004571 lime Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 239000003086 colorant Substances 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 239000005995 Aluminium silicate Substances 0.000 description 2
- 108091005950 Azurite Proteins 0.000 description 2
- 244000025254 Cannabis sativa Species 0.000 description 2
- 235000012766 Cannabis sativa ssp. sativa var. sativa Nutrition 0.000 description 2
- 235000012765 Cannabis sativa ssp. sativa var. spontanea Nutrition 0.000 description 2
- 229920000742 Cotton Polymers 0.000 description 2
- 238000001069 Raman spectroscopy Methods 0.000 description 2
- 235000012211 aluminium silicate Nutrition 0.000 description 2
- 235000009120 camo Nutrition 0.000 description 2
- 235000005607 chanvre indien Nutrition 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 239000011487 hemp Substances 0.000 description 2
- 238000003709 image segmentation Methods 0.000 description 2
- NLYAJNPCOHFWQQ-UHFFFAOYSA-N kaolin Chemical compound O.O.O=[Al]O[Si](=O)O[Si](=O)O[Al]=O NLYAJNPCOHFWQQ-UHFFFAOYSA-N 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- GWBUNZLLLLDXMD-UHFFFAOYSA-H tricopper;dicarbonate;dihydroxide Chemical compound [OH-].[OH-].[Cu+2].[Cu+2].[Cu+2].[O-]C([O-])=O.[O-]C([O-])=O GWBUNZLLLLDXMD-UHFFFAOYSA-H 0.000 description 2
- 241000283707 Capra Species 0.000 description 1
- -1 Dunhuang slub Substances 0.000 description 1
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 description 1
- 241001676573 Minium Species 0.000 description 1
- 241000209140 Triticum Species 0.000 description 1
- 235000021307 Triticum Nutrition 0.000 description 1
- 229910000004 White lead Inorganic materials 0.000 description 1
- 238000002441 X-ray diffraction Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- IRERQBUNZFJFGC-UHFFFAOYSA-L azure blue Chemical compound [Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Na+].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[S-]S[S-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] IRERQBUNZFJFGC-UHFFFAOYSA-L 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000001055 blue pigment Substances 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 229910052951 chalcopyrite Inorganic materials 0.000 description 1
- DVRDHUBQLOKMHZ-UHFFFAOYSA-N chalcopyrite Chemical compound [S-2].[S-2].[Fe+2].[Cu+2] DVRDHUBQLOKMHZ-UHFFFAOYSA-N 0.000 description 1
- 229910052956 cinnabar Inorganic materials 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000003292 glue Substances 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 239000001056 green pigment Substances 0.000 description 1
- 238000000227 grinding Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 239000011133 lead Substances 0.000 description 1
- HTUMBQDCCIXGCV-UHFFFAOYSA-N lead oxide Chemical compound [O-2].[Pb+2] HTUMBQDCCIXGCV-UHFFFAOYSA-N 0.000 description 1
- 239000001035 lead pigment Substances 0.000 description 1
- YEXPOXQUZXUXJW-UHFFFAOYSA-N lead(II) oxide Inorganic materials [Pb]=O YEXPOXQUZXUXJW-UHFFFAOYSA-N 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000001053 orange pigment Substances 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 229910052957 realgar Inorganic materials 0.000 description 1
- 239000001054 red pigment Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000010902 straw Substances 0.000 description 1
- 239000002345 surface coating layer Substances 0.000 description 1
- 238000004876 x-ray fluorescence Methods 0.000 description 1
- 239000001052 yellow pigment Substances 0.000 description 1
Landscapes
- Spectrometry And Color Measurement (AREA)
Abstract
一种基于可见光谱成像的古代壁画颜料无损全面识别方法,包括依据壁画使用的基础颜料和双常数Kubelka‑Munk理论构建全面的壁画颜料数据库,作为壁画颜料识别的参考依据;获取古代壁画的可见光谱图像;对可见光谱图像进行渲染得到对应的彩色RGB图像;利用超像素分割算法对彩色RGB图像进行精细分割;计算每个子区域的平均光谱反射率;分别提取数据库样本和子区域颜料的光谱特征;构建壁画颜料识别方法,识别每个子区域颜料的物质种类,若为单色颜料,进一步识别其平均粒径大小,若为混合颜料,则判定为最小粒径颜料,并进一步计算每种颜料的浓度比例;获得每种颜料在壁画中的二维分布图;最终完成包括颜料物质、粒径、浓度及分布的古代壁画颜料全面无损识别。
Description
技术领域
本发明属于历史文化遗产保护技术领域,具体涉及一种基于可见光谱成像的古代壁画颜料无损全面识别方法。
背景技术
我国拥有大量的古代彩绘壁画,如敦煌莫高窟、陕西榆林和新疆克孜尔壁画等,它们是我国古代文明的重要组成部分,这些古代壁画蕴含着我国古代政治、经济、文化、科学技术及生产工艺发展等重要历史信息,承载着灿烂文明,传承历史文化,是研究我国古代人文学、社会学、美学以及自然科学的珍贵资源。然而受自然环境与人为因素的多重影响,我国各地大部分壁画颜料层都产生了多种病害问题,导致其可显现的历史价值信息正在快速消减,若不及时采取科学合理的保护措施,将会对国家和民族造成不可挽回的损失。
天然矿物颜料是古代壁画彩色图案层的主要物质成分,是壁画显色的唯一物质基础,因此,对壁画颜料进行科学的识别分析是深入了解这些历史文物的重要前提,也是开展古代壁画保护的重要基础。与此同时,对彩绘文物颜料分析和表征将有助于文物保护工作者了解古代艺术家创作壁画的颜料使用情况、分析学习其绘画技法、掌握彩绘文物当前的颜色保存状态,无论对于采用真实颜料对彩绘文物的实物修复、复原临摹,或是对于彩绘文物数字化的虚拟修复与再现等历史文化传承工作而言,都具有非常重要的意义。
在古代壁画颜料识别方面,传统的颜料识别多基于有损采样分析方法进行,如采用X-射线衍射和X-射线荧光分析仪、扫描电子显微镜、能量色谱仪、偏光显微镜、傅里叶变换红外光谱仪、拉曼光谱仪、拉曼显微光谱仪、激光诱导击穿光谱仪以及光学显微镜等仪器和设备,分析采样物质的物质成分信息,此类采样分析方法对宝贵的历史文物造成不可逆转的损害,并且仅能进行单点的物质成分识别,无法进行壁画全幅面颜料分布以及颜料物理属性的全面分析,具有较大局限性。近年来,原位无损分析技术越来越被文物保护工作者所重视,随之产生了诸如光纤反射光谱仪以及便携式拉曼纤维光谱仪等设备,使得颜料的分析无需破坏式采样,该类技术降低了损害文物的风险,但仪器在使用过程中仍需要与文物本体接触,且仅能进行单点分析,在壁画颜料的分析应用领域仍具有一定局限性。除此之外,当前的颜料分析主要集中于颜料物质成分的识别,而忽略了影响壁画颜色的颜料粒径大小和混合颜料中各个基础颜料浓度等因素,对于壁画颜料的分析不够全面。
对于以上古代壁画颜料识别存在的局限性问题,目前文博界和学术界均未提出合理有效的解决方法。光谱成像技术具有“谱图合一”的优势,本发明依据日趋成熟的光谱成像技术,以及矿物颜料物质、粒径属性等与光谱反射率之间关联性,提出了一种基于可见光谱成像的古代壁画颜料无损全面识别方法,该方法无需与壁画接触,即可实现包括颜料物质类别、平均粒径大小、混合浓度及空间分布的无损全面识别,有效地解决了当前分析技术对于古代壁画颜料分析不全面的局限性问题。
发明内容
本发明的目的是为了解决背景技术中所述问题,提出一种基于可见光谱成像的古代壁画颜料无损全面识别方法。
本发明的技术方案为一种基于可见光谱成像的古代壁画颜料无损全面识别方法,具体包括以下步骤:
步骤1,获得古代壁画基础颜料并构建壁画颜料数据库,所述数据库中包含基础颜料样本、基础颜料样本光谱数据以及依据双常数Kubelka-Munk理论构建的混合颜料样本的光谱数据;
步骤2,获取古代壁画可见光谱图像;
步骤3,对古代壁画可见光谱图像进行渲染,获得其对应的彩色RGB图像;
步骤4,对渲染得到的古代壁画彩色RGB图像进行超像素分割,得到分割图像,并获得彩色图像各子分区的标记数值;
步骤5,依据古代壁画彩色图像子分区标记及可见光谱图像,提取每个子区域各像素光谱反射率,并计算每个子区域的平均光谱反射率;
步骤6,提取数据库中样本及古代壁画彩色图像每个子区域的光谱反射率特征;
步骤7,依据步骤6中获得的数据库中样本和古代壁画彩色图像每个子区域光谱反射率特征,通过特征匹配方式对每个子区域颜料物质属性进行识别;
步骤8,若通过步骤7识别出的颜料为单色颜料,则进一步对该颜料平均粒径大小进行分析,实现颜料平均粒径大小的识别;
步骤9,若通过步骤7识别出的颜料为混合颜料,则判定基础颜料为最小粒径等级,并进一步对该混合颜料中各基础颜料的浓度进行分析,得到各基础颜料的最优浓度值;
步骤10,依据步骤8和步骤9中每个子区域颜料识别结果,对每种基础颜料分布进行标记,获得古代壁画中各基础颜料的二维分布图;
步骤11,完成古代壁画颜料识别,获得包含颜料物质属性、粒径大小、浓度比例及空间分布的古代壁画颜料无损全面识别结果。
而且,步骤1中基础颜料样本的光谱数据采用分光光度计测量获得,混合颜料样本的光谱数据由双常数Kubelka-Munk理论计算得到,具体实现方式如下,
A.制作用于求解基础颜料单位浓度吸收和散射系数的样本,包括基础颜料纯色样本,一个参考白色颜料纯色样本,基础颜料纯色样本按照特定质量比例与参考白色颜料混合的样本;
B.求解上述样本的吸收散射比,
利用式一求解上述所有样本的吸收散射比,
其中,表示样本的吸收散射比,Rλ表示利用分光光度计测量获得的样本光谱反射率,λ表示波长;
C.求解基础颜料的单位浓度吸收和散射系数;
设定参考白色颜料在整个光谱波长范围内的单位散射系数为1,
sλ,w=1, 式二
那么,参考白色颜料的吸收系数为,
其中,为参考白色颜料的吸收散射比;然后利用式四计算用于求解基础颜料单位浓度吸收和散射系数的中间值,
其中,sλ,0表示用于求解基础颜料单位浓度吸收和散射系数的中间值,表示混合颜料吸收散射比,表示纯色颜料吸收散射比;
那么基础颜料的单位浓度吸收和散射系数按式五和式六求解获得,
其中,ctint表示混合颜料中基础颜料所占质量比,sλ,paint表示基础颜料的单位浓度散射系数,kλ,paint表示基础颜料的单位浓度吸收系数;
D.计算混合样本光谱反射率,完成数据库构建,
利用基础颜料的单位浓度吸收和散射系数,按照式七和式八,计算任意一个按给定浓度混合得到的样本光谱反射率,
其中,表示混合颜料样本的吸收散射比,ci表示任意一个基础颜料样本浓
度,且 表示混合颜料总吸收系数,表示混合颜料总散射系数,
n表示混合颜料成分数量,Rλ,mixture表示混合颜料的光谱反射率。
而且,步骤2中获得古代壁画可见光谱的实现方式为,
相同均匀照明拍摄环境和相同拍摄参数条件下,分别拍摄获得古代壁画和光谱重建训练样本集的数字图像,提取训练样本集数字响应信号,并利用分光光度计测量获得训练样本集光谱数据,利用训练样本数据计算由数字响应值到光谱反射率转换的变化矩阵,然后利用转换矩阵对古代壁画进行光谱重建,获得敦煌壁画的可见光谱图像,具体方法如下:
步骤2.1,首先将古代壁画各像素和训练样本数字响应值按式九所示三阶多项式进行扩展
d=[1 r g b rg rb gb r2 g2 b2 rg2 r2g rb2 r2b gb2 g2b r3 g3 b3],式九
其中,d表示一个像素或者一个训练样本响应值扩展后向量,r、g和b分别表示数码相机R-通道、G-通道和B-通道的响应值;
步骤2.2,利用训练样本集和伪逆方法求解用于重建古代壁画可见光谱图像的转换矩阵,如式十所示,
Q=RtrainDtrain + 式十
其中,Q表示光谱重建转换矩阵,Rtrain表示训练样本集的光谱反射率矩阵,Dtrain表示训练样本集的扩展响应矩阵,上标‘+’表示伪逆运算符号;
步骤2.3,利用式十求得的光谱重建转换矩阵及古代壁画数字响应扩展矩阵进行壁画的光谱重建,获得古代壁画可见光谱图像,如式十一所示,
Rmural=QDmural 式十一
其中,Rmural表示重建得到的古代壁画可见光谱图像,Dmural表示古代壁画数字响应扩展矩阵。
而且,步骤3中对古代壁画可见光谱图像进行渲染,获得其对应的彩色RGB图像的实现方式如下,
步骤3.1,对于古代壁画可见光谱图像的任意一个像素点光谱反射率,首先在设定参考计算光源和标准观察者条件下计算其对应的三刺激值XYZ,如式十二所示,
其中,x(λ),y(λ),和z(λ)表示CIE标准观察者三刺激值函数,E(λ)表示CIE标准照明体相对光谱功率,r(λ)表示古代壁画可见光谱图像的任意一个像素点光谱反射率,k是调整因子,X、Y和Z表示计算得到的三刺激值;
步骤3.2,依据计算得到的三刺激值计算相应参数条件下的RGB值,如式十三所示,
其中,R、G和B为渲染后彩色图像的颜色值,M为对应三刺激值计算条件下由XYZ到RGB转换的系数矩阵。
而且,步骤4中采用基于k-means的简单线性迭代聚类方法对渲染得到的古代壁画彩色图像进行超像素分割。
而且,步骤6中提取数据库中样本及古代壁画彩色图像每个子区域光谱反射率特征的实现方式如下,
步骤6.1,首先对数据库中样本光谱反射率进行去均值归一化处理,对于任意一个样本,归一化方法如式十四所示,
其中,rnomalized表示样本归一化光谱反射率,ri表示样本某一波长处反射率值,表示样本反射率的平均值,l表示采用波段数量;
步骤6.2,求解数据库中样本光谱反射率的一阶和二阶导数,并利用式十四对数据库中样本光谱反射率的一阶导数和二阶导数进行归一化处理;
步骤6.3,将归一化处理之后的样本光谱反射率、光谱反射率一阶导数及二阶导数数据组合成为一维向量,即得到样本的光谱特征数据;
步骤6.4,针对壁画彩色图像每个子区域的平均光谱反射率,重复步骤6.1-6.3,获得古代壁画彩色图像每个子区域光谱特征数据。
本发明从壁画颜料的物质及物理属性与可见光谱反射率的关联性角度出发,提出了将壁画可见光谱图像获取技术与彩色图像分割技术以及双常数Kubelka-Munk理论等进行有机结合,构建合理有效的古代壁画颜料鉴别方法,实现了古代壁画颜料的多属性无损全面识别。本发明相对于传统及现有古代壁画颜料识别方法而言:1)克服了传统及现有技术在壁画颜料分析中对壁画造成的潜在危害,真正实现了壁画颜料的无损分析。2)弥补了传统及现有技术单点分析的缺陷性,实现了壁画颜料的全幅面识别。3)克服了传统及现有技术仅对颜料物质属性识别的局限性,实现了颜料粒径大小及浓度识别。由于本发明技术方案具有重要应用意义,受到多个项目支持:1.国家自然基金项目61275172,2.国家自然基金项目61575147,3.国家973基础研究子项目2012CB725302。对本发明技术方案进行保护,将对我国相关行业竞争国际领先地位具有重要意义。
附图说明
图1为本发明实施例的流程图。
具体实施方式
本发明技术方案具体实施时可由本领域技术人员采用计算机软件技术运行。结合附图,提供本发明实施例具体描述如下。
如图1所示,实施例提供了一种基于可见光谱成像的古代壁画颜料无损全面识别方法,可以有效地解决当前方法存在的局限性,完善古代壁画颜料识别的全面性,用于古代壁画保护领域。实施例以敦煌莫高窟249窟北壁中部经变中的飞天壁画为实验对象,创建敦煌壁画颜料样本数据库,获取249窟北壁中部经变中的飞天壁画光谱图像,对本发明方法进行说明。需要说明的是,本发明并不仅仅局限于上述实验对象的应用支持,对于任意的古代壁画对象本发明同样适用。
实施例主要包括以下步骤:
步骤1,调研总结敦煌壁画中的基础颜料种类及壁画创作技法,获得基础颜料,掌握样本制作方法;依据双常数Kubelka-Munk理论及获得的敦煌壁画基础颜料,构建敦煌壁画颜料数据库。
通过传统采样分析方法对敦煌壁画颜料研究成果的调研整理可知,敦煌壁画主要使用红、橙、黄、绿、蓝、黑和白等色系颜料,其中红色颜料主要包括土红、朱砂和朱磦等颜料,橙色主要包括铅丹和雄黄等颜料,黄色颜料主要为雌黄、石黄、土黄和密陀僧,绿色颜料主要为石绿和氯铜矿,蓝色颜料包括石青、群青和青金石等,黑色颜料为炭黑,白色颜料主要为石膏、高岭土和铅白等,加上金粉以及敦煌细泥两种金色和土色,共计约20余种颜料。在实际的敦煌壁画创作中,古代画家首先用混有麦秸及少量砂砾的粗泥将岩壁填平,待粗泥层干燥后再附一层混有棉麻纤维的细泥,然后在干燥细泥层表面用石膏、高岭土等白色颜料打底,最后以牛皮胶调和以上所述颜料进行彩绘,完成壁画创作。
本实施例中所述壁画基础颜料均为敦煌研究院提供。此外,敦煌研究院美术所专家教授如何按照敦煌壁画创作工艺制作样本方法,为最大程度接近壁画实际创作例程,本实施例选取2.5mm厚的工业装裱卡纸作为基底,裁切为2cm大小方形块,揭去其表面覆膜层之后形成毛刺表面作为粗泥层,然后在基底上刷涂取自敦煌大泉河床且掺揉了棉麻纤维的细泥,待细泥层干燥后刷涂一层石膏打底,最后针对上述基础颜料,利用羊毫白云笔在干燥白底表面进行绘制,待其自然干燥,完成样本制作。
数据库中包含基础颜料样本和依据双常数Kubelka-Munk理论构建的混合颜料样本的光谱数据。其中,基础颜料样本按照上述方法制作,其光谱数据采用GretagMacbethColor-Eye7000A分光光度计测量获得,混合颜料样本的光谱数据由双常数Kubelka-Munk理论计算得到。在计算混合样本光谱数据时,首先以基础颜料中最白的颜料为参考,采用相对值方法求解基础颜料的单位浓度吸收和散射系数,然后计算混合颜料样本光谱数据,方法详述如下:
A.制作用于求解基础颜料单位浓度吸收和散射系数的样本。
除作参考颜料的白色颜料之外,对于任意一个基础颜料,按步骤1制作方法,分别制作一个纯色样本(用下标masstone进行标示)和一个按照特定质量比例与参考白色颜料混合的样本(用下标tint进行标示),另外单独制作一个参考白色颜料纯色样本(用下标w进行标示),利用GretagMacbeth Color-Eye 7000A分光光度计测量获得所有样本的光谱反射率。
B.求解样本的吸收散射比。
利用式一求解上述所有样本的吸收散射比,如下:
其中,表示样本的吸收散射比,Rλ表示利用分光光度计测量获得的样本光谱反射率,λ表示波长。
C.求解基础颜料的单位浓度吸收和散射系数。
设定参考白色颜料在整个光谱波长范围内的单位散射系数为1,即:
sλ,w=1,式二
那么,参考白色颜料的吸收系数为:
其中,为参考白色颜料的吸收散射比。然后利用式四计算用于求解基础颜料单位浓度吸收和散射系数的中间值:
其中,sλ,0表示用于求解基础颜料单位浓度吸收和散射系数的中间值,表示混合颜料吸收散射比,表示纯色颜料吸收散射比。
那么基础颜料的单位浓度吸收和散射系数则可按式五和式六求解获得:
其中,ctint表示混合颜料中基础颜料所占质量比,sλ,paint表示基础颜料的单位浓度散射系数,kλ,paint表示基础颜料的单位浓度吸收系数。
D.计算混合样本光谱反射率,完成数据库构建。
利用基础颜料的单位浓度吸收和散射系数,按照式七和式八,计算任意一个按给定浓度混合得到的样本光谱反射率:
其中,表示混合颜料样本的吸收散射比,ci表示任意一个基础颜料样本浓度,且 表示混合颜料总吸收系数,表示混合颜料总吸收系数,n表示混合颜料成分数量,Rλ,mixture表示混合颜料的光谱反射率。
通过以上步骤求解,即可完成敦煌壁画颜料数据库的构建,作为壁画颜料识别的参考依据。实施例中采用铅白颜料作为参考白颜料,用于求解其它基础颜料的单位浓度吸收和散射系数。实施例中所有基础颜料的tint混合颜料样本配比均为0.55参考白颜料加0.45基础颜料。此外,依据对敦煌壁画颜料使用情况的研究文献调查可知,敦煌壁画中的颜料最多为三成分混合,因此混合颜料样本包含二成分颜料混合和三成分颜料混合两种形式,其中二成分颜料混合均采用0.1A:0.9B、0.3A:0.7B、0.5A:0.5B、0.7A:0.3B、0.9A:0.1B五种配比形式,三成分颜料混合均采用0.1A:0.1B:0.8C、0.1A:0.8B:0.1C、0.8A:0.1B:0.1C、0.2A:0.2B:0.6C、0.2A:0.6B:0.2C、0.6A:0.2B:0.2C、0.3A:0.3B:0.4C、0.3A:0.4B:0.3C、0.4A:0.3B:0.3C、0.1A:0.3B:0.6C、0.3A:0.1B:0.6C、0.6A:0.3B:0.1C、0.26:0.1B:0.3C、0.1A:0.6B:0.3C、0.3A:0.6B:0.1C十五种配比形式,其中A、B和C分别表示一种基础颜料。
步骤2,获取敦煌壁画可见光谱图像。
利用佳能700D数码相机在相同均匀照明拍摄环境和相同拍摄参数条件下,分别拍摄敦煌壁画和光谱重建训练样本集的数字图像,本实施例中采用敦煌彩绘文物数字化保护色卡作为训练样本集,有关色卡信息参见文献:梁金星,万晓霞,刘强,等.敦煌彩绘文物数字化保护色卡制作方法研究[J].敦煌研究,2016(2):117-124.,在MATLAB中提取训练样本集数字响应信号,利用GretagMacbeth Color-Eye 7000A分光光度计测量获得训练样本集光谱数据,利用训练样本数据计算由数字响应值到光谱反射率转换的变化矩阵,然后利用转换矩阵对敦煌壁画进行光谱重建,获得敦煌壁画的可见光谱图像。具体方法如下所述:
A.首先将敦煌壁画各像素和训练样本数字响应值按式九所示三阶多项式进行扩展
d=[1 r g b rg rb gb r2 g2 b2 rg2 r2g rb2 r2b gb2 g2b r3 g3 b3],式九
其中,d表示一个像素或者一个训练样本响应值扩展后向量,r、g和b分别表示数码相机R-通道、G-通道和B-通道的响应值。
B.利用训练样本集和伪逆方法求解用于重建敦煌壁画可见光谱图像的转换矩阵,如式十所示:
Q=RtrainDtrain + 式十
其中,Q表示光谱重建转换矩阵,Rtrain表示训练样本集的光谱反射率矩阵,Dtrain表示训练样本集的扩展响应矩阵,上标‘+’表示伪逆运算符号。
C.利用式十求得的光谱重建转换矩阵及敦煌壁画数字响应扩展矩阵进行壁画的光谱重建,获得敦煌壁画可见光谱图像,如式十一所示:
Rmural=QDmural 式十一
其中,Rmural表示重建得到的敦煌壁画可见光谱图像,Dmural表示敦煌壁画数字响应扩
展矩阵。
步骤3,对敦煌壁画可见光谱图像进行渲染,获得其对应的彩色RGB图像,以应用于后续彩色图像超像素分割处理。具体渲染方法如下:
A.对于敦煌壁画可见光谱图像的任意一个像素点光谱反射率,首先在设定参考计算光源和标准观察者条件下计算其对应的三刺激值XYZ,如式十二所示:
其中,x(λ),y(λ),和z(λ)表示CIE标准观察者三刺激值函数,E(λ)表示CIE标准照明体相对光谱功率,r(λ)表示敦煌壁画可见光谱图像的任意一个像素点光谱反射率,k是调整因子,X、Y和Z表示计算得到的三刺激值。
B.依据计算得到的三刺激值计算相应参数条件下的RGB值,如式十三所示:
其中,R、G和B为渲染后彩色图像的颜色值,M为对应三刺激值计算条件下由XYZ到RGB转换的系数矩阵。
本实施例中,以上计算均可通过MATLAB颜色科学工具箱完成,分别选用CIE D50标准光源和CIE1931标准观察者函数作为计算参数,获得敦煌壁画可见光谱图像在上述参数条件下的渲染彩色图像。
步骤4,对渲染得到的敦煌壁画彩色图像进行超像素分割,得到分割图像,获得彩色图像各子分区的标记数值。
本实施例中采用基于k-means的简单线性迭代聚类方法,进行敦煌壁画彩色渲染图像项进行超像素分割,该方法相对于其他图像分割方法而言,具有边缘轮廓保持好、计算效率高和分割精度高等优点,由于该方法为现有技术,此处不作赘述,具体原理参见文献:Achanta R,Shaji A,Smith K,et al.SLIC Superpixels Compared to State-of-the-ArtSuperpixel Methods[J].IEEE Transactions on Pattern Analysis&MachineIntelligence,2012,34(11):2274-2282.。完成壁画彩色图像分割之后,即可获得各子分区的标记数值。
步骤5,依据敦煌壁画彩色图像子分区标记及可见光谱图像,提取每个子区域各像素光谱反射率,并计算每个子区域的平均光谱反射率。
本实施例中,设置超像素分割子区域个数为1500,实际子分区个数为1455针对每个子区域,依据该区域像素在彩色图像中的位置坐标及壁画可见光谱图像,提取每个子区域各像素光谱反射率,并计算每个子区域的平均光谱反射率,用于下一步的颜料识别。
步骤6,提取数据库样本及待识别颜料光谱(即敦煌壁画彩色图像的每个子区域)反射率特征值,首先求解数据库样本和待识颜料光谱反射率的一阶和二阶导数,然后对光谱反射率、一阶导数和二阶导数进行去均值归一化处理,最后将归一化处理之后的光谱反射率、光谱反射率一阶导数、光谱反射率二阶导数组合成为一维向量,获得数据库样本及待识别颜料的光谱特征。具体方法如下
A.首先对数据库样本和待识别颜料光谱反射率进行去均值归一化处理,对于任意一个样本,归一化方法如式十四所示:
其中,rnomalized表示样本归一化光谱反射率,ri表示样本某一波长处反射率值,表示样本反射率的平均值,l表示采用波段数量。
B.同样利用式十四对数据库样本和待识别颜料光谱反射率的一阶导数和二阶导数进行归一化处理。
C.将归一化处理之后的样本光谱反射率、光谱反射率一阶导数及二阶导数数据组合成为一维向量,即得到样本及待识别颜料的光谱特征数据。
通过以上步骤,即可完成数据库样本和待识别颜料光谱特征的提取。本实施例光谱反射率波长范围为400nm-700nm,采样间隔为10nm,光特征数据为31维归一化光谱反射率加30维归一化一阶导数和29维归一化二阶导数书,光谱特征数据共90维。
步骤7,敦煌壁画颜料物质成分识别。
依据步骤6中获得的数据库样本和待识别颜料光谱特征,通过特征匹配方式进行每个子区域颜料物质属性的识别,具体方法为:以待识别颜料与数据库中参考样本光谱特征数据之间的欧式距离和夹角为基础,构建匹配度误差评价指标。其中欧式距离反映待识别颜料特征数据与参考样本特征数据之间的叠合程度,夹角反映待识别颜料特征数据与参考样本特征数据之间的形状相似度。针对每一个待识别颜料,分别计算其与每一个参考样本的特征数据欧式距离和夹角,对欧氏距离和夹角进行最大值归一化处理,然后相乘得到匹配误差,待识别颜料的特征数据与参考样本的数据匹配误差越小,证明颜料光谱特征数据的匹配度越高,二者物质属性越接近,以此确定待识别颜料的物质成分属性。
本实施例中敦煌壁画颜料的鉴别结果主要包括石青、石绿、土红、炭黑、敦煌细泥以及石膏等,包括石青与石绿的混合、石膏与敦煌细泥的混合、炭黑与石膏的混合等形式。
步骤8,敦煌壁画单色颜料平均粒径大小分析。
若通过步骤7识别出的颜料为单色颜料,则进一步对该颜料平均粒径大小进行分析,具体方法为:针对基础颜料数据库中具有粒径等级的每一种颜料,按照传统研磨和筛选工艺,制备不同平均粒径等级的颜料,利用激光粒径测量仪测量获得每种等级颜料的粒径分布及平均粒径值,同时按照壁画创作工艺涂制颜料样本,利用GretagMacbeth Color-Eye7000A分光光度计测量获得每个粒径等级的光谱反射率,通过数学拟合建模方式,建立颜料平均粒径和颜料样本平均光谱反射率之间的函数关系,利用待识别颜料的平均光谱反射率,即可实现颜料平均粒径大小的识别。
本实施例中石青、石绿和土红的颜料平均粒径识别结果均与数据库中各自的最小号11#等级颜料最为接近,对应平均颜料粒径大小计算值约为25微米。敦煌细泥、炭黑与石膏等颜料非矿物颜料,无粒径等级之分,无需粒径大小的识别。
步骤9,敦煌壁画混合颜料中个基础颜料浓度分析。
若通过步骤7识别出的颜料为混合颜料,则判定基础颜料为最小粒径等级,并进一步对该混合颜料中各基础颜料的浓度进行分析,具体方法为:依据步骤2计算得到的基础颜料单位浓度吸收和散射系数,将基础颜料浓度设置为变量,以参考样本的各基础颜料浓度为初始浓度值,结合双常数Kubelka-Munk理论,以计算得到的混合颜料光谱与实际混合颜料之间的光谱均方根误差作为目标函数,利用非负最小二乘约束优化方法,通过光谱均方根最小化求解得到各基础颜料的最优浓度值。目标函数如式十五所示:
其中,表示在变换基础颜料浓度时,由式七和式八计算得到的估计混合颜料光谱,rmixture为待识别混合颜料的光谱。本实施例中混合颜料中个基础颜料浓度在壁画不同区域并不完全相同,此处不作具体说明。
步骤10,依据步骤8和步骤9中每个子区域颜料识别结果,对每种基础颜料分布进行标记,获得敦煌壁画中各基础颜料的二维分布图。
步骤11,完成敦煌壁画颜料识别,获得包含颜料物质属性、粒径大小、浓度比例及空间分布的敦煌壁画颜料无损全面识别结果。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (6)
1.一种基于可见光谱成像的古代壁画颜料无损全面识别方法,其特征在于,包括以下步骤:
步骤1,获得古代壁画基础颜料并构建壁画颜料数据库,所述数据库中包含基础颜料样本、基础颜料样本光谱数据以及依据双常数Kubelka-Munk理论构建的混合颜料样本的光谱数据;
步骤2,获取古代壁画可见光谱图像;
步骤3,对古代壁画可见光谱图像进行渲染,获得其对应的彩色RGB图像;
步骤4,对渲染得到的古代壁画彩色RGB图像进行超像素分割,得到分割图像,并获得彩色图像各子分区的标记数值;
步骤5,依据古代壁画彩色图像子分区标记及可见光谱图像,提取每个子区域各像素光谱反射率,并计算每个子区域的平均光谱反射率;
步骤6,提取数据库中样本及古代壁画彩色图像每个子区域的光谱反射率特征;
步骤7,依据步骤6中获得的数据库中样本和古代壁画彩色图像每个子区域光谱反射率特征,通过特征匹配方式对每个子区域颜料物质属性进行识别;
步骤8,若通过步骤7识别出的颜料为单色颜料,则进一步对该颜料平均粒径大小进行分析,实现颜料平均粒径大小的识别;
步骤9,若通过步骤7识别出的颜料为混合颜料,则判定基础颜料为最小粒径等级,并进一步对该混合颜料中各基础颜料的浓度进行分析,得到各基础颜料的最优浓度值;
步骤10,依据步骤8和步骤9中每个子区域颜料识别结果,对每种基础颜料分布进行标记,获得古代壁画中各基础颜料的二维分布图;
步骤11,完成古代壁画颜料识别,获得包含颜料物质属性、粒径大小、浓度比例及空间分布的古代壁画颜料无损全面识别结果。
2.如权利要求1所述的一种基于可见光谱成像的古代壁画颜料无损全面识别方法,其特征在于:步骤1中基础颜料样本的光谱数据采用分光光度计测量获得,混合颜料样本的光谱数据由双常数Kubelka-Munk理论计算得到,具体实现方式如下,
A.制作用于求解基础颜料单位浓度吸收和散射系数的样本,包括基础颜料纯色样本,一个参考白色颜料纯色样本,基础颜料纯色样本按照特定质量比例与参考白色颜料混合的样本;
B.求解上述样本的吸收散射比,
利用式一求解上述所有样本的吸收散射比,
其中,表示样本的吸收散射比,Rλ表示利用分光光度计测量获得的样本光谱反射率,λ表示波长;
C.求解基础颜料的单位浓度吸收和散射系数;
设定参考白色颜料在整个光谱波长范围内的单位散射系数为1,
sλ,w=1, 式二
那么,参考白色颜料的吸收系数为,
其中,为参考白色颜料的吸收散射比;
然后利用式四计算用于求解基础颜料单位浓度吸收和散射系数的中间值,
其中,sλ,0表示用于求解基础颜料单位浓度吸收和散射系数的中间值,表示混合颜料吸收散射比,表示纯色颜料吸收散射比;
那么基础颜料的单位浓度吸收和散射系数按式五和式六求解获得,
其中,ctint表示混合颜料中基础颜料所占质量比,sλ,paint表示基础颜料的单位浓度散射系数,kλ,paint表示基础颜料的单位浓度吸收系数;
D.计算混合样本光谱反射率,完成数据库构建,
利用基础颜料的单位浓度吸收和散射系数,按照式七和式八,计算任意一个按给定浓度混合得到的样本光谱反射率,
其中,表示混合颜料样本的吸收散射比,ci表示任意一个基础颜料样本浓度,且表示混合颜料总吸收系数,表示混合颜料总散射系数,
n表示混合颜料成分数量,Rλ,mixture表示混合颜料的光谱反射率。
3.如权利要求1或2所述的一种基于可见光谱成像的古代壁画颜料无损全面识别方法,其特征在于:步骤2中获得古代壁画可见光谱的实现方式为,
相同均匀照明拍摄环境和相同拍摄参数条件下,分别拍摄获得古代壁画和光谱重建训练样本集的数字图像,提取训练样本集数字响应信号,并利用分光光度计测量获得训练样本集光谱数据,利用训练样本数据计算由数字响应值到光谱反射率转换的变化矩阵,然后利用转换矩阵对古代壁画进行光谱重建,获得敦煌壁画的可见光谱图像,具体方法如下:
步骤2.1,首先将古代壁画各像素和训练样本数字响应值按式九所示三阶多项式进行扩展
d=[1 r g b rg rb gb r2 g2 b2 rg2 r2g rb2 r2b gb2 g2b r3 g3 b3], 式九
其中,d表示一个像素或者一个训练样本响应值扩展后向量,r、g和b分别表示数码相机R-通道、G-通道和B-通道的响应值;
步骤2.2,利用训练样本集和伪逆方法求解用于重建古代壁画可见光谱图像的转换矩阵,如式十所示,
Q=RtrainDtrain + 式十
其中,Q表示光谱重建转换矩阵,Rtrain表示训练样本集的光谱反射率矩阵,Dtrain表示训练样本集的扩展响应矩阵,上标‘+’表示伪逆运算符号;
步骤2.3,利用式十求得的光谱重建转换矩阵及古代壁画数字响应扩展矩阵进行壁画的光谱重建,获得古代壁画可见光谱图像,如式十一所示,
Rmural=QDmural 式十一
其中,Rmural表示重建得到的古代壁画可见光谱图像,Dmural表示古代壁画数字响应扩展矩阵。
4.如权利要求3所述的一种基于可见光谱成像的古代壁画颜料无损全面识别方法,其特征在于:步骤3中对古代壁画可见光谱图像进行渲染,获得其对应的彩色RGB图像的实现方式如下,
步骤3.1,对于古代壁画可见光谱图像的任意一个像素点光谱反射率,首先在设定参考计算光源和标准观察者条件下计算其对应的三刺激值XYZ,如式十二所示,
其中,x(λ),y(λ),和z(λ)表示CIE标准观察者三刺激值函数,E(λ)表示CIE标准照明体相对光谱功率,r(λ)表示古代壁画可见光谱图像的任意一个像素点光谱反射率,k是调整因子,X、Y和Z表示计算得到的三刺激值;
步骤3.2,依据计算得到的三刺激值计算相应参数条件下的RGB值,如式十三所示,
其中,R、G和B为渲染后彩色图像的颜色值,M为对应三刺激值计算条件下由XYZ到RGB转换的系数矩阵。
5.如权利要求4所述的一种基于可见光谱成像的古代壁画颜料无损全面识别方法,其特征在于:步骤4中采用基于k-means的简单线性迭代聚类方法对渲染得到的古代壁画彩色图像进行超像素分割。
6.如权利要求5所述的一种基于可见光谱成像的古代壁画颜料无损全面识别方法,其特征在于:步骤6中提取数据库中样本及古代壁画彩色图像每个子区域光谱反射率特征的实现方式如下,
步骤6.1,首先对数据库中样本光谱反射率进行去均值归一化处理,对于任意一个样本,归一化方法如式十四所示,
其中,rnomalized表示样本归一化光谱反射率,ri表示样本某一波长处反射率值,表示样本反射率的平均值,l表示采用波段数量;
步骤6.2,求解数据库中样本光谱反射率的一阶和二阶导数,并利用式十四对数据库中样本光谱反射率的一阶导数和二阶导数进行归一化处理;
步骤6.3,将归一化处理之后的样本光谱反射率、光谱反射率一阶导数及二阶导数数据组合成为一维向量,即得到样本的光谱特征数据;
步骤6.4,针对壁画彩色图像每个子区域的平均光谱反射率,重复步骤6.1-6.3,获得古代壁画彩色图像每个子区域光谱特征数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710781294.XA CN107643267B (zh) | 2017-09-01 | 2017-09-01 | 一种基于可见光谱成像的古代壁画颜料无损全面识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710781294.XA CN107643267B (zh) | 2017-09-01 | 2017-09-01 | 一种基于可见光谱成像的古代壁画颜料无损全面识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107643267A CN107643267A (zh) | 2018-01-30 |
CN107643267B true CN107643267B (zh) | 2019-06-25 |
Family
ID=61110143
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710781294.XA Expired - Fee Related CN107643267B (zh) | 2017-09-01 | 2017-09-01 | 一种基于可见光谱成像的古代壁画颜料无损全面识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107643267B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108267426B (zh) * | 2018-02-06 | 2023-09-01 | 中国科学院西安光学精密机械研究所 | 基于多光谱成像的绘画颜料识别系统及方法 |
CN108520488B (zh) * | 2018-04-10 | 2022-06-24 | 深圳劲嘉集团股份有限公司 | 一种重构光谱并进行复制的方法以及电子设备 |
CN109974625B (zh) * | 2019-04-08 | 2021-02-09 | 四川大学 | 一种基于色相优化灰度的彩色物体结构光三维测量方法 |
CN110020998A (zh) * | 2019-04-12 | 2019-07-16 | 曲阜师范大学 | 基于光谱的褪变色中国书画颜色复原方法 |
CN110441244A (zh) * | 2019-07-26 | 2019-11-12 | 北京建筑大学 | 一种顾及吸收特征的光谱分段矿物颜料识别方法 |
CN110740537B (zh) * | 2019-09-30 | 2021-10-29 | 宁波燎原照明集团有限公司 | 一种博物馆文物的光照系统自适应调节的系统 |
CN111008653B (zh) * | 2019-11-18 | 2023-04-18 | 西安建筑科技大学 | 一种混合颜料信息解混的聚类优化方法 |
CN114519814A (zh) * | 2020-10-31 | 2022-05-20 | 华为技术有限公司 | 一种图像识别的方法及电子设备 |
CN112561812B (zh) * | 2020-12-07 | 2023-07-25 | 西安建筑科技大学 | 一种基于颜料褪色光谱色差模型的壁画颜色复原方法 |
CN112632790B (zh) * | 2020-12-29 | 2021-12-21 | 愉悦家纺有限公司 | 多维度耦合-叠加复合混色模型与渐变色谱矩阵算法构建应用 |
CN113160077B (zh) * | 2021-04-08 | 2022-05-13 | 武汉纺织大学 | 一种褪变色彩色壁画颜色高保真数字化修复方法 |
CN113792710A (zh) * | 2021-11-15 | 2021-12-14 | 滨州学院 | 光谱重建方法、装置及电子设备 |
CN118362513B (zh) * | 2024-06-19 | 2024-09-03 | 天津大学 | 对彩绘文物进行材料属性分析的方法、装置及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006053070A (ja) * | 2004-08-12 | 2006-02-23 | Kyoto Univ | 色素の推定方法 |
CN102062725A (zh) * | 2010-11-24 | 2011-05-18 | 西安交通大学 | 一种馆藏壁画表面颜料特征的无损分析方法 |
CN102305769A (zh) * | 2011-06-09 | 2012-01-04 | 天津大学 | 应用于中国古画修复的多光谱抠图方法 |
CN102854147A (zh) * | 2012-08-08 | 2013-01-02 | 北京建筑工程学院 | 一种基于高光谱数据的壁画底稿信息提取方法 |
CN103528968A (zh) * | 2013-11-01 | 2014-01-22 | 上海理工大学 | 基于迭代阈值法的反射率光谱重建方法 |
CN106469258A (zh) * | 2016-09-28 | 2017-03-01 | 武汉大学 | 一种基于双常数Kubelka‑Munk理论的有色纤维混合配色方法 |
-
2017
- 2017-09-01 CN CN201710781294.XA patent/CN107643267B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006053070A (ja) * | 2004-08-12 | 2006-02-23 | Kyoto Univ | 色素の推定方法 |
CN102062725A (zh) * | 2010-11-24 | 2011-05-18 | 西安交通大学 | 一种馆藏壁画表面颜料特征的无损分析方法 |
CN102305769A (zh) * | 2011-06-09 | 2012-01-04 | 天津大学 | 应用于中国古画修复的多光谱抠图方法 |
CN102854147A (zh) * | 2012-08-08 | 2013-01-02 | 北京建筑工程学院 | 一种基于高光谱数据的壁画底稿信息提取方法 |
CN103528968A (zh) * | 2013-11-01 | 2014-01-22 | 上海理工大学 | 基于迭代阈值法的反射率光谱重建方法 |
CN106469258A (zh) * | 2016-09-28 | 2017-03-01 | 武汉大学 | 一种基于双常数Kubelka‑Munk理论的有色纤维混合配色方法 |
Non-Patent Citations (3)
Title |
---|
基于可见光谱特征提取的敦煌壁画颜料识别方式研究;补雅晶等;<文物保护与考古科学>;20170630;第29卷(第3期);第45-51页 |
基于可见光谱的古代壁画颜料无损鉴别方法;梁金星等;<光谱法与光谱分析>;20170831;第37卷(第8期);第2519-2626页 |
敦煌壁画颜料颜色数据库构建方法;梁金星等;<敦煌研究>;20170228(第1期);第132-140页 |
Also Published As
Publication number | Publication date |
---|---|
CN107643267A (zh) | 2018-01-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107643267B (zh) | 一种基于可见光谱成像的古代壁画颜料无损全面识别方法 | |
Attas et al. | Near-infrared spectroscopic imaging in art conservation: investigation of drawing constituents | |
Cucci et al. | Remote-sensing hyperspectral imaging for applications in archaeological areas: Non-invasive investigations on wall paintings and on mural inscriptions in the Pompeii site | |
EP1814044B1 (en) | Paint pigment identification method, identification system, identification program, and recording medium thereof | |
Escadafal | Remote sensing of soil color: principles and applications | |
EP2833327B1 (en) | Method and system for digitally generating appearance data | |
Lau et al. | Multispectral and hyperspectral image analysis of elemental and micro-Raman maps of cross-sections from a 16th century painting | |
CN106052877B (zh) | 壁画绘画材料多光谱图像色彩标准色卡及其制作方法 | |
Navas et al. | Benefits of applying combined diffuse reflectance FTIR spectroscopy and principal component analysis for the study of blue tempera historical painting | |
JPH11211569A (ja) | メタリツク塗色の分類整理方法 | |
CN108389238A (zh) | 一种彩绘类文物混合颜料的分析方法 | |
Dal Fovo et al. | Insights into the stratigraphy and palette of a painting by Pietro Lorenzetti through non-invasive methods | |
Xu et al. | Can deep learning assist automatic identification of layered pigments from XRF data? | |
CN110441244A (zh) | 一种顾及吸收特征的光谱分段矿物颜料识别方法 | |
CN112488997B (zh) | 基于特征插值的古代绘画印刷品颜色复现检测和评价方法 | |
CN113160077B (zh) | 一种褪变色彩色壁画颜色高保真数字化修复方法 | |
Zhang et al. | Fast identification of mural pigments at Mogao Grottoes using a LIBS-based spectral matching algorithm | |
Creagh et al. | Recent and future developments in the use of radiation for the study of objects of cultural heritage significance | |
Lyu et al. | Quantitative analysis of mixed pigments for Chinese paintings using the improved method of ratio spectra derivative spectrophotometry based on mode | |
CN112561812B (zh) | 一种基于颜料褪色光谱色差模型的壁画颜色复原方法 | |
Ciortan et al. | Tensor decomposition for painting analysis. Part 1: pigment characterization | |
Delaney et al. | Visible and Infrared Reflectance Imaging Spectroscopy of Paintings and Works on Paper | |
CN118362513B (zh) | 对彩绘文物进行材料属性分析的方法、装置及设备 | |
Deborah et al. | Old Man in Warnemünde (1907) colouring palette: A case study on the use of hyperspectral imaging for pigment identification | |
Lau et al. | Artworks and Cultural Heritage Materials: Using Multivariate Analysis to Answer Conservation Questions |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190625 Termination date: 20200901 |