CN108389238A - A kind of analysis method of colored drawing class historical relic hybrid pigment - Google Patents
A kind of analysis method of colored drawing class historical relic hybrid pigment Download PDFInfo
- Publication number
- CN108389238A CN108389238A CN201810259228.0A CN201810259228A CN108389238A CN 108389238 A CN108389238 A CN 108389238A CN 201810259228 A CN201810259228 A CN 201810259228A CN 108389238 A CN108389238 A CN 108389238A
- Authority
- CN
- China
- Prior art keywords
- pigment
- spectrum
- colored drawing
- historical relic
- end member
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Spectrometry And Color Measurement (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
The invention discloses a kind of analysis methods of colored drawing class historical relic hybrid pigment, and the broadband continuous reflection spectrum of pending colored drawing class cultural artifact surface area-of-interest is acquired first with high light spectrum image-forming and field spectroradiometer;Then the broadband continuous reflection spectrum is considered as migration fractionation information, realizes that the pigment in same color region decomposes using Algorithms of Non-Negative Matrix Factorization, obtains the reflection spectrum curve of different pigment;The reflection spectrum curve of the different pigment is matched with the spectrum in pigment library of spectra again, obtains the type of hybrid pigment;The abundance of various species pigment is obtained using the Ratio spectrum derivative based on mode, and then obtains the type and ratio of pending colored drawing class cultural artifact surface color composition pigment.The above method can permanently retain the pigment digital information of colored drawing class historical relic, contribute to the digital virtual reparation to historical relic and entity reparation.
Description
Technical field
The present invention relates to historical relic's protections and identification technology field more particularly to a kind of analysis of colored drawing class historical relic hybrid pigment
Method.
Background technology
Currently, in historical relic's protection and identification technology field, hybrid pigment analysis is to excavate ancient book to draw inherent value, just for it
True analysis and assessment provide the technology on scientific basis basis.However thousand paper longevity, thin,tough silk longevity 800, even if without environment, artificial etc.
The influence of factor, ancient book draw it is now mouldy, damage by worms, stain, the various situations such as breakage it is also very universal, in order to preferably retain and
The value of the ancient painting and calligraphy pieces of displaying, generally requires to carry out carrying out color reparation and duplication to it, this is just it needs to be determined that the kind of pigment
Class.
There are many science and technology to be applied to the pigment category authentication of colored drawing class historical relic in the prior art, such as X-ray diffraction point
The technologies such as analysis, x-ray fluorescence analysis, Raman spectrum analysis, electronic micro-analysis, near-infrared spectrum analysis, by measuring pigment
The features such as element, structure and external pattern determine the classification of substance, but the method for the prior art in most cases needs
Want sample preparation that could accurately be measured, and process of the sample preparation inherently to a secondary damage of historical relic, since historical relic can not
Regenerated particularity should be lossless for technology used by Analysis of Cultural Relic, and also lack in the prior art to mixing face
Expect the solution of analysis.
Invention content
The object of the present invention is to provide a kind of analysis method of colored drawing class historical relic hybrid pigment, this method can be retained permanently
The pigment digital information of colored drawing class historical relic contributes to the digital virtual reparation to historical relic, and effectively shortens cultural relic entity
Repair time has larger economic benefit and social benefit to the digital protection and reparation of colored drawing class historical relic.
The purpose of the present invention is what is be achieved through the following technical solutions:
A kind of analysis method of colored drawing class historical relic hybrid pigment, the method includes:
Step 1 acquires pending colored drawing class cultural artifact surface area-of-interest using high light spectrum image-forming and field spectroradiometer
Broadband continuous reflection spectrum;
The broadband continuous reflection spectrum is considered as migration fractionation information by step 2, is realized using Algorithms of Non-Negative Matrix Factorization
The pigment in same color region decomposes, and obtains the reflection spectrum curve of different pigment;
Step 3 matches the reflection spectrum curve of the different pigment with the spectrum in pigment library of spectra, is mixed
Close the type of pigment;
Step 4 recycles the Ratio spectrum derivative based on mode to obtain the abundance of various species pigment, and then obtains and wait locating
Manage the type and ratio of colored drawing class cultural artifact surface color composition pigment.
As seen from the above technical solution provided by the invention, the above method can permanently retain the face of colored drawing class historical relic
Material digitalized information contributes to the digital virtual reparation to historical relic, and effectively shortens repair time of cultural relic entity, to colored drawing
The digital protection of class historical relic has larger economic benefit and social benefit with reparation.
Description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, required use in being described below to embodiment
Attached drawing be briefly described, it should be apparent that, drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill in field, without creative efforts, other are can also be obtained according to these attached drawings
Attached drawing.
Fig. 1 is the analysis method flow diagram of colored drawing class historical relic hybrid pigment provided in an embodiment of the present invention.
Specific implementation mode
With reference to the attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based on this
The embodiment of invention, every other implementation obtained by those of ordinary skill in the art without making creative efforts
Example, belongs to protection scope of the present invention.
The embodiment of the present invention is described in further detail below in conjunction with attached drawing, is implemented as shown in Figure 1 for the present invention
The analysis method flow diagram for the colored drawing class historical relic hybrid pigment that example provides, the method includes:
Step 1 acquires pending colored drawing class cultural artifact surface area-of-interest using high light spectrum image-forming and field spectroradiometer
Broadband continuous reflection spectrum;
In this step, the coverage area of broadband continuous reflection spectrum is 350nm-2500nm, and acquisition environment is darkroom,
The light source carried using instrument, measurement are repeatedly averaged, to deep or light region of differing more than at same color acquisition at least two
Reflectance spectrum;The broadband continuous reflection spectrum is that the spectral reflectance of a variety of pure pigment by synthesizing color mixes.
The broadband continuous reflection spectrum is considered as migration fractionation information by step 2, is realized using Algorithms of Non-Negative Matrix Factorization
The pigment in same color region decomposes, and obtains the reflection spectrum curve of different pigment;
Here, it is non-negative that Algorithms of Non-Negative Matrix Factorization, which requires all elements in matrix, in the case where not considering error,
Matrix is expressed as it is that two nonnegative matrixes are multiplied completely as a result, just because of decomposition the result is that non-negative feature, thus
There is specific physical significance in practical application, since high-spectral data is a kind of record to atural object reflectivity data, meet
The requirement of nonnegativity, so as to the Endmember extraction by Algorithms of Non-Negative Matrix Factorization for high-spectral data.
In the specific implementation, Algorithms of Non-Negative Matrix Factorization NMF is indicated using following formula:
Xn×m=An×rSr×m+En×m
Wherein,For the mixed spectra matrix of broadband continuous reflection spectrum to be decomposed;WithFor two end member matrixes after decomposition;En×mIt is expressed as error matrix;R values are end member number;And | | E | | to the greatest extent may be used
Energy is small, here x ∈ R+Indicate x >=0;
In the influence for not considering error, the problem of for solving two unknown quantity A and S by a known quantity X, often turn
Turn to following formula to solve, i.e.,:
min||X-AS||;
Further, two end member matrix As and S are solved by objective function and selection Optimality Criteria, it is defined
Object function includes minimum Euclidean distance and K-L divergences, wherein:
The minimum Euclidean distance is expressed as:
The K-L divergences are expressed as:
Wherein i (i=1,2 ... n) be mixed spectra matrix X line number (i.e. spectrum channel), maximum value is exactly n, and n is
Total points of a measured curve of spectrum, it is related with instrument performance, it is total wave band number of instrument;J (j=1,2 ...
M) be mixed spectra matrix X row number, the maximum value of j is m, and m, which is represented, selects several mixed reflection spectrum to participate in operations.
In conjunction with the object function, continuous alternating iteration calculates the value of end member matrix A and S, until being optimal solution, changes
It is expressed as formula:
From the point of view of the formula of Non-negative Matrix Factorization, which is to search out one group of base vector so that object function is minimum
Change, for the separation of hybrid pigment spectrum, each pure pigment component spectrum is exactly such one group of base vector, thus is solved
Go out this group of base vector, is equivalent to obtain the endmember spectra for constituting mixed spectra.
Step 3 matches the reflection spectrum curve of the different pigment with the spectrum in pigment library of spectra, is mixed
Close the type of pigment;
Step 4 recycles the Ratio spectrum derivative based on mode to obtain the abundance of various species pigment, and then obtains this and wait for
Handle the type and ratio of colored drawing class cultural artifact surface color composition pigment.
Here, the ratio spectro-derivative spectrophotometry based on mode is built upon on the basis of original Ratio spectrum derivative, it is also contemplated that
There are the factor of non-linear wave band, by improving the choosing method of characteristic wave bands in ratio spectro-derivative spectrophotometry, introduce mode this
Statistic concept for characterizing its general characteristic, and then determines abundance.
In the specific implementation, above-mentioned obtain the process of the abundance of various species pigment using the Ratio spectrum derivative based on mode
For:
Statistics section is determined according to the relationship between value of chromatism and pigment abundance first;
Then by the Abundances of the various species pigment solved based on Ratio spectrum derivative according to constraint item non-negative and more than 1
After part rejecting abnormalities value, the number fallen in statistics section is calculated, and preliminary Abundances model is regarded as into section where it
It encloses;
It is again proportional spacing with step-length 0.01, builds a series of simulation mixed spectras, calculating simulation spectrum and measured spectra
Related coefficient, the highest one group of end member ratio of related coefficient is the final Abundances of various species pigment.
It is built upon on the basis of mixed spectra meets linear mixed model in the specific implementation, derivative ratio solution mixes algorithm,
And meet the constraints that end member abundance is non-negative and the sum of end member abundance is 1.When mixture only includes two kinds of mineral constituent,
And when not considering limits of error item, the model of the Ratio spectrum derivative is expressed as:
r(λi)=F1×r1(λi)+F2×r2(λi)
Wherein, i is spectrum channel, r (λi) be mixed spectra in λiThe reflectivity of wavelength location, F1、F2For the rich of end member
Degree, r1(λi)、r2(λi) it is respectively end member F1、F2It is corresponding in λiThe reflectivity of wavelength location;
By end member F2It is considered as interfering component, above-mentioned expression formula both sides while divided by F2The reflectivity of endmember spectra, obtains
Again to r (λi) first derivation, it obtains
It can be seen that at this time that derivative spectrum has been and F by the formula2The content of end member is unrelated, and then both sides are removed simultaneously again
WithObtain F1The abundance of end member;
Again by F1End member passes through Initial abundance value of the abundance result as various species pigment after second order derivation, specifically
Expression formula is:
Ratio spectro-derivative spectrophotometry solution procedure is can be seen that, it is simple clear, complicated interative computation is not needed, it can be fast
Speed, the inverting for accurately carrying out endmember spectra abundance.
It is worth noting that, the content not being described in detail in the embodiment of the present invention belongs to professional and technical personnel in the field's public affairs
The prior art known.
In conclusion the method that the embodiment of the present invention is provided can permanently retain the pigment digital information of historical relic, have
Help the digital virtual reparation to historical relic, improve the artistic expression effect of historical relic, is conducive to the on-line off-line digitlization of historical relic
Displaying;In addition, entity reparation of the abundance inverting of pigment to historical relic, especially to its panchromatic process be capable of providing science according to
According to can effectively shorten repair time of cultural relic entity, there are a large amount of collection colored drawing class cultural relic digitalizations to protect and repair China
There is larger economic benefit and social benefit again.
The foregoing is only a preferred embodiment of the present invention, but scope of protection of the present invention is not limited thereto,
Any one skilled in the art is in the technical scope of present disclosure, the change or replacement that can be readily occurred in,
It should be covered by the protection scope of the present invention.Therefore, protection scope of the present invention should be with the protection model of claims
Subject to enclosing.
Claims (5)
1. a kind of analysis method of colored drawing class historical relic hybrid pigment, which is characterized in that the method includes:
Step 1, the wide wave that pending colored drawing class cultural artifact surface area-of-interest is acquired using high light spectrum image-forming and field spectroradiometer
Section continuous reflection spectrum;
The broadband continuous reflection spectrum is considered as migration fractionation information by step 2, is realized using Algorithms of Non-Negative Matrix Factorization identical
The pigment of color region decomposes, and obtains the reflection spectrum curve of different pigment;
Step 3 matches the reflection spectrum curve of the different pigment with the spectrum in pigment library of spectra, obtains mixing face
The type of material;
Step 4 recycles the Ratio spectrum derivative based on mode to obtain the abundance of various species pigment, and then it is pending to obtain this
Colored drawing class cultural artifact surface color forms the type and ratio of pigment.
2. the analysis method of colored drawing class historical relic hybrid pigment according to claim 1, which is characterized in that in step 1,
The coverage area of the broadband continuous reflection spectrum is 350nm-2500nm;And the broadband continuous reflection spectrum be by
The spectral reflectance for synthesizing a variety of pure pigment of color mixes.
3. the analysis method of colored drawing class historical relic hybrid pigment according to claim 1, which is characterized in that in step 2,
The Algorithms of Non-Negative Matrix Factorization is indicated using following formula:
Xn×m=An×rSr×m+En×m
Wherein,For the mixed spectra matrix of broadband continuous reflection spectrum to be decomposed;With
For two end member matrixes after decomposition;En×mIt is expressed as error matrix;R values are end member number;
Further, two end member matrix As and S, defined target are solved by objective function and selection Optimality Criteria
Function includes minimum Euclidean distance and K-L divergences, wherein:
The minimum Euclidean distance is expressed as:
The K-L divergences are expressed as:
Wherein i (i=1,2 ... be n) mixed spectra matrix X line number, maximum value is n, and n is a measured spectrum song
Total points of line;J (j=1,2 ... m) be mixed spectra matrix X row number, maximum value is m, and m represents several mixing of selection
Reflectance spectrum participates in operation;
In conjunction with the object function, continuous alternating iteration calculates the value of end member matrix A and S, and until being optimal solution, iteration is public
Formula is expressed as:
4. the analysis method of colored drawing class historical relic hybrid pigment according to claim 1, which is characterized in that in step 4, described
The abundance that various species pigment is obtained using the Ratio spectrum derivative based on mode, is specifically included:
Statistics section is determined according to the relationship between value of chromatism and pigment abundance first;
Then the Abundances of the various species pigment solved based on Ratio spectrum derivative are picked according to non-negative and more than 1 constraints
After exceptional value, the number fallen in statistics section is calculated, and preliminary Abundances range is regarded as into section where it;
It is again proportional spacing with step-length 0.01, builds a series of simulation mixed spectras, the phase of calculating simulation spectrum and measured spectra
Relationship number, the highest one group of end member ratio of related coefficient is the final Abundances of various species pigment.
5. the analysis method of colored drawing class historical relic hybrid pigment according to claim 4, which is characterized in that the Ratio spectrum derivative
Model be expressed as:
r(λi)=F1×r1(λi)+F2×r2(λi)
Wherein, i is spectrum channel, r (λi) be mixed spectra in λiThe reflectivity of wavelength location, F1、F2For the abundance of end member, r1
(λi)、r2(λi) it is respectively end member F1、F2It is corresponding in λiThe reflectivity of wavelength location;
By end member F2It is considered as interfering component, above-mentioned expression formula both sides while divided by F2The reflectivity of endmember spectra, obtains
Again to r (λi) first derivation, it obtains
Then again both sides simultaneously divided byObtain F1The abundance of end member;
By F1End member passes through Initial abundance value of the abundance result as various species pigment after second order derivation, expression
For:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810259228.0A CN108389238A (en) | 2018-03-27 | 2018-03-27 | A kind of analysis method of colored drawing class historical relic hybrid pigment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810259228.0A CN108389238A (en) | 2018-03-27 | 2018-03-27 | A kind of analysis method of colored drawing class historical relic hybrid pigment |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108389238A true CN108389238A (en) | 2018-08-10 |
Family
ID=63072671
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810259228.0A Pending CN108389238A (en) | 2018-03-27 | 2018-03-27 | A kind of analysis method of colored drawing class historical relic hybrid pigment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108389238A (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297918A (en) * | 2018-08-21 | 2019-02-01 | 广西科技大学 | A method of detection Acid red 26 |
CN109523600A (en) * | 2018-11-09 | 2019-03-26 | 广东省博物馆 | A kind of historical relic's protection reparation recruits color method |
CN110441244A (en) * | 2019-07-26 | 2019-11-12 | 北京建筑大学 | A kind of spectrum segmentation coloring earth recognition methods for taking Absorption Characteristics into account |
CN113160077A (en) * | 2021-04-08 | 2021-07-23 | 武汉纺织大学 | High-fidelity digital restoration method for color of faded colored mural |
CN115128017A (en) * | 2022-08-11 | 2022-09-30 | 天津大学 | Colored drawing cultural relic illumination damage assessment method and monitoring method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102193090A (en) * | 2010-03-19 | 2011-09-21 | 复旦大学 | Mixed pixel decomposition method for remote sensing images |
CN105224915A (en) * | 2015-09-07 | 2016-01-06 | 河海大学 | A kind of mixed pixel of hyper-spectral image decomposition method |
-
2018
- 2018-03-27 CN CN201810259228.0A patent/CN108389238A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102193090A (en) * | 2010-03-19 | 2011-09-21 | 复旦大学 | Mixed pixel decomposition method for remote sensing images |
CN105224915A (en) * | 2015-09-07 | 2016-01-06 | 河海大学 | A kind of mixed pixel of hyper-spectral image decomposition method |
Non-Patent Citations (2)
Title |
---|
侯妙乐 等: "《高光谱成像技术在彩绘文物分析中的研究综述》", 《光谱学与光谱分析》 * |
吕书强 等: "《基于众数的比值导数法在混合颜料解混中的研究》", 《光散射学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109297918A (en) * | 2018-08-21 | 2019-02-01 | 广西科技大学 | A method of detection Acid red 26 |
CN109523600A (en) * | 2018-11-09 | 2019-03-26 | 广东省博物馆 | A kind of historical relic's protection reparation recruits color method |
CN109523600B (en) * | 2018-11-09 | 2021-08-31 | 广东省博物馆 | Cultural relic protection repair complementary color matching method |
CN110441244A (en) * | 2019-07-26 | 2019-11-12 | 北京建筑大学 | A kind of spectrum segmentation coloring earth recognition methods for taking Absorption Characteristics into account |
CN113160077A (en) * | 2021-04-08 | 2021-07-23 | 武汉纺织大学 | High-fidelity digital restoration method for color of faded colored mural |
CN113160077B (en) * | 2021-04-08 | 2022-05-13 | 武汉纺织大学 | High-fidelity digital restoration method for color of faded colored mural |
CN115128017A (en) * | 2022-08-11 | 2022-09-30 | 天津大学 | Colored drawing cultural relic illumination damage assessment method and monitoring method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108389238A (en) | A kind of analysis method of colored drawing class historical relic hybrid pigment | |
CN107643267B (en) | A kind of lossless comprehensive recognition methods of ancient wall pigment based on visible spectrum imaging | |
WO2016000088A1 (en) | Hyperspectral waveband extraction method based on optimal index factor-correlation coefficient method | |
KR101725700B1 (en) | Texture analysis of a painted surface using specular angle data | |
CN105051762B (en) | System and method for determining formulation for coating material | |
AU2014346520B2 (en) | Texture analysis of a coated surface using pivot-normalization | |
CN103612483B (en) | A kind of printing ink color matching method based on Spectral matching | |
Nidamanuri et al. | Use of field reflectance data for crop mapping using airborne hyperspectral image | |
US10969952B2 (en) | Color and texture match ratings for optimal match selection | |
CN108427934A (en) | A kind of Hyperspectral imaging mixed pixel decomposition method | |
CN113189021A (en) | Method for identifying rock color based on spectrum | |
CN109946241A (en) | A kind of classification of soils method calculating imaging system based on EO-1 hyperion | |
CN104268896B (en) | Hyper spectrum dimensionality reduction matching method and system based on spectrum sampling histogram | |
CN110070004A (en) | A kind of field hyperspectrum Data expansion method applied to deep learning | |
CN110926608A (en) | Spectrum reconstruction method based on light source screening | |
CN115597726A (en) | Color temperature detection method and device and electronic equipment | |
Marchant et al. | Spectral invariance under daylight illumination changes | |
He et al. | Method for grade identification of tobacco based on machine vision | |
CN108519157B (en) | A kind of generation of metamerism spectrum and evaluation method and system for light source detection | |
CN111750993B (en) | Open measurement environment spectrum measurement method based on imaging condition correction | |
Kuntz et al. | Karlsruhe optimized and precise radiative transfer algorithm: III. ADDLIN and TRANSF algorithms for modeling spectral transmittance and radiance | |
Contreras et al. | Geochemical and Hyperspectral Data Fusion for Drill-Core Mineral Mapping | |
CN113191259B (en) | Dynamic data expansion method for hyperspectral image classification and image classification method | |
Djiongo Kenfack et al. | Quaternion-based texture analysis of multiband satellite images: application to the estimation of aboveground biomass in the east region of Cameroon | |
Tominaga | Estimation Method of Fluorescent Donaldson Matrices Based on Multispectral Imaging Data |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180810 |
|
RJ01 | Rejection of invention patent application after publication |