CN103134789A - Spectrum recover method based on Laplacian-Markov field - Google Patents
Spectrum recover method based on Laplacian-Markov field Download PDFInfo
- Publication number
- CN103134789A CN103134789A CN2012104756108A CN201210475610A CN103134789A CN 103134789 A CN103134789 A CN 103134789A CN 2012104756108 A CN2012104756108 A CN 2012104756108A CN 201210475610 A CN201210475610 A CN 201210475610A CN 103134789 A CN103134789 A CN 103134789A
- Authority
- CN
- China
- Prior art keywords
- spectrum
- standard deviation
- laplacian
- overbar
- iteration
- 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
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000001237 Raman spectrum Methods 0.000 claims abstract description 22
- 239000011159 matrix material Substances 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 6
- 230000003595 spectral effect Effects 0.000 claims description 40
- 238000011084 recovery Methods 0.000 claims description 9
- 238000005316 response function Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008859 change Effects 0.000 abstract description 11
- 238000010606 normalization Methods 0.000 abstract description 4
- 238000004458 analytical method Methods 0.000 abstract description 3
- 238000001514 detection method Methods 0.000 abstract description 3
- 239000000463 material Substances 0.000 abstract description 3
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 23
- 239000008103 glucose Substances 0.000 description 23
- 238000004422 calculation algorithm Methods 0.000 description 13
- 238000001069 Raman spectroscopy Methods 0.000 description 6
- 230000001629 suppression Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 125000002496 methyl group Chemical group [H]C([H])([H])* 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
The invention discloses a spectrum recover method based on the Laplacian-Markov field. The a spectrum recover method based on the Laplacian-Markov field comprises the following steps: conducting normalization processing on discretized spectrum to obtain normalized spectrum intensity; calculating a first derivative of the normalized spectrum intensity; calculating a neighborhood standard deviation of each normalized sectrum fi, finding out a maximum standard deviation and a minimum standard deviation, constructing a weighting matrix according to maximum standard deviation and the minimum standard deviation; and using the Laplacian-Markov field, and solving a raman spectrum deconvolution and split Bregman iteration method. The spectrum recover method based on the Laplacian-Markov field suppresses noise in places gentle in intensity change, adopts weak reatrain to store detail of the spectrum in places drastic in intensity change, achieves balance of recover resolution ratio and noise compress capability, introduces the split Bregman iteration method to overcome indifferentiable problem of conventional gradient descenting method, and provides technical support for application of the raman spectrum in quality detection, material analysis and the like.
Description
Technical Field
The invention belongs to the technical field of Raman spectrum data processing, and particularly relates to a method for recovering a degenerated Raman spectrum by utilizing a Laplace-Markov Laplacian-Markov) field, a deconvolution and splitting Bregman iteration method, which provides technical support for the Raman spectrum in the aspects of quality detection, material analysis and the like.
Background
In recent years, the research on raman spectroscopy technology is active, so that raman spectroscopy is widely applied to the fields of food quality monitoring, chemical rapid detection, material analysis, biomedical science and the like. However, during the measurement process of the spectrometer, the spectral data is often affected by factors such as natural broadening of spectral lines, collisional broadening, low-pass characteristics of detectors and circuits, and the like, which causes strong noise and sub-band overlapping phenomena of the spectral data. These noise interferences cause difficulties in the spectral application and can even lead to serious identification and estimation errors, which greatly limits the application of the spectrum. Therefore, it is essential to preprocess data and improve the resolution of the spectrum.
At present, there are many methods for recovering the degraded spectrum, and these methods can be mainly classified into two main categories: curve fitting and spectral deconvolution. The curve fitting method considers that the spectrum curve is formed by combining a Gaussian function and a Lorentzian function, namely, the measured spectrum can be obtained by fitting the functions. The method can well solve the problem of subband overlapping, namely, overlapping peaks can be fitted, but the method is extremely sensitive to noise. Spectral deconvolution methods, mainly to recover the degraded spectrum by different constraints, are High Order Statistics (HOS), Tikhonov regularization, homomorphic filtering, etc. In these methods, strong constraints can well suppress raman spectral noise (Tikhonov regularization), but too strong constraints can cause that severely degraded spectral details cannot be recovered. The weaker constraint can better recover the degraded details (higher order statistics), but the noise suppression capability is reduced. Therefore, such methods are difficult to balance in recovering spectral resolution and noise suppression capability.
Disclosure of Invention
The invention provides a spectrum recovery method based on a Laplacian-Markov field, which solves the problem that Raman spectrum data measured by a spectrometer is often influenced by random noise, instrument errors and the like to cause low resolution.
A spectrum recovery method based on a Laplacian-Markov field is carried out according to the following steps:
(1) for discretized spectrumNormalized processing is carried out to obtain normalized spectral intensity fiI is 1,2,3 … n, n is the total number of spectra;
(2) calculating normalized spectral intensity fiFirst derivative f ofi′;
(3) Calculating each spectral point fiNeighborhood of [ f ]i-3,fi-2,fi-1,fi,fi+1,fi+2,fi+3]Standard deviation of (S)iFrom which the maximum standard deviation S is foundmaxAnd minimum standard deviation SminConstructing a weighting matrix Q of n x n dimensions, the diagonal elements of the weighting matrix QSetting the other elements to be 0, wherein ln is a natural logarithm, and e is a natural constant;
(4) and (3) solving the Raman spectrum by adopting a splitting iterative method:
(41) solving parameters by an initialization splitting method: initial value f of the spectrum0F, the initial value b of the over-noise-removal compensation 00, substitute variable initial value d 00, and 0 is the iteration number k;
(42) updating calculation fk+1And dk+1:
Where max is the maximum value, H is the instrument response function, HTIs the transpose of H, Δ t is the time step, n × n dimensional difference matrix of the spectrumNormalized spectral intensity seti ═ 1,2, …, n, regularization parameters λ ∈ (0,1) and β ∈ (0, 1);
(43) update calculation bk+1=bk+((Df)k+1-dk+1);
(44) Judgment of fk+1Whether or not an iteration stop condition (| f) is satisfiedk+1-fk‖/‖fk|) > epsilon, epsilon is an iteration stop threshold, if not, then return to step (42), otherwise output raman spectrum fk+1。
The instrument response function H is a Gaussian function or a Lorentzian function.
The technical effects of the invention are as follows: first, a method of spectral deconvolution with data weighting is proposed. For the raman spectrum data curve, the influence of interference on the data itself is not the same for all data, and the influence of noise is larger for the wave number with gentle intensity change and smaller for the wave number with severe intensity change. The data weighting adopts different constraints according to different intensity changes, and adopts stronger constraint at the place with gentle intensity change to restrain good noise; the details of the good spectrum are preserved by taking weaker constraints where the intensity changes are severe. This achieves a balance between the ability to restore resolution and noise suppression. Secondly, a Laplacian-Markov prior is applied in the solution of the spectrum deconvolution, and is a weak smoothness prior which can well keep the details of the spectrum, so that the Laplacian-Markov prior is first applied in the deconvolution of the Raman spectrum data. Thirdly, a splitting Bregman iteration method is introduced to solve the deconvolution model, the splitting Bregman iteration method overcomes the infinitesimal problem of the conventional gradient descent method, and compared with the conventional method, the method has extremely high algorithm solving speed and high practical value.
Drawings
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a flow chart of a Laplacian-Markov spectrum deconvolution algorithm;
FIG. 3 is a diagram of a mathematical model of spectral degradation;
FIG. 4 shows glucose (C)15H29NO7) Actually measuring a Raman spectrum curve;
FIG. 5 shows glucose (C)15H29NO7) Actually measuring a spectrum curve after normalization of the Raman spectrum curve;
FIG. 6 is a neighborhood standard deviation S output spectral curve for each spectral point;
FIG. 7 is a spectral curve of each spectral point normalized by the neighborhood standard deviation S;
FIG. 8 is a diagram of kernel functions of different shapes having the same full width at half maximum;
FIG. 9 shows glucose (C)15H29NO7) Result graphs of different iteration times of the spectrum; wherein,
FIG. 9(a) shows glucose (C)15H29NO7) Performing spectrum curve original drawing;
FIG. 9(b) shows glucose (C)15H29NO7) 3 times of iteration result graphs of the spectral curve;
FIG. 9(C) shows glucose (C)15H29NO7) 5 times of iteration result graphs of the spectrum curve;
FIG. 9(d) shows glucose (C)15H29NO7) 7 times of iteration result graphs of the spectral curve;
FIG. 9(e) shows glucose (C)15H29NO7) A spectrum curve 10 times iteration result graph;
FIG. 9(f) shows glucose (C)15H29NO7) A spectrum curve 20 times iteration result graph;
FIG. 9(g) shows glucose (C)15H29NO7) A spectrum curve 30 times iteration result graph;
FIG. 9(h) shows glucose (C)15H29NO7) A spectrum curve 40 times of iteration result graphs;
FIG. 9(i) is glucose (C)15H29NO7) A spectrum curve 50 times iteration result graph;
FIG. 9(j) is glucose (C)15H29NO7) A spectrum curve 60 times iteration result graph;
FIG. 9(k) shows glucose (C)15H29NO7) A spectrum curve 70 times of iteration result graphs;
FIG. 9(l) is glucose (C)15H29NO7) A spectrum curve 80 times iteration result graph;
FIG. 9(m) shows glucose (C)15H29NO7) A spectrum curve 90 times iteration result graph;
FIG. 9(n) is glucose (C)15H29NO7) A spectrum curve 100 times iteration result graph;
FIG. 10 is a comparison of the spectrum recovery results of degraded glucose by different methods, wherein 10(a) is a graph of the actually measured spectrum of a glucose instrument, 10(b) is a graph of the experimental result of the TRSD method, and 10(c) is a graph of the experimental result of the HOS method, and 10(d) is a graph of the experimental result of the LM-QSD method;
Detailed Description
The present invention is described in further detail below with reference to the attached drawings and examples.
For the raman spectrum data curve, the influence of interference on the data itself is not the same for all data, and the influence of noise is larger for the wave number with gentle intensity change and smaller for the wave number with severe intensity change. If the adaptive processing is performed according to the difference of the intensity change, the balance of the resolution restoring capability and the noise suppression capability can be achieved. Therefore, the invention provides a spectrum deconvolution method with data weighting, wherein the data weighting adopts different constraints according to different intensity changes, and adopts stronger constraint at a place with mild intensity change to inhibit noise; the details of the good spectrum are preserved by taking weaker constraints where the intensity changes are severe. On the other hand, a Laplacian-Markov prior is a weak smoothness prior. The smoothness of the prior constrained Raman spectrum is adopted, the Laplacian-Markov prior is applied to deconvolution of Raman spectrum data for the first time, and an experimental result also proves that the method is very effective for recovering the degraded spectrum.
Most raman spectral data measured by raman spectroscopy can be described as a convolution of the true signal with the corresponding function of the measurement instrument, as shown in fig. 2.
g(v)=Hf(v)+n(v),
Wherein, f (v) is a real signal, g (v) is an observed signal measured by a measuring instrument, and H is a response function of the measuring instrument, also called a fuzzy core, and mainly comprises an internal linear function and an instrument broadening function.
As shown in FIG. 3, FIG. 3 is glucose (C)6H12O6) The measured Raman spectrum curve has the abscissa Wavenumber of the curve as the wave number, the ordinate Raman Intensity as the Raman spectrum Intensity, and the resolution ratio of the spectrum curve is 1cm-1The curve is 2850 to 2900cm-1The wavenumber range has severe sub-band overlapping phenomena and is affected by noise, and the spectrum can be used for testing the resolution recovery capability of the algorithm on the spectrum curve and the noise suppression capability.
Referring to fig. 1, the method of the present invention specifically comprises the following steps:
(1) degraded spectrum to discretizationPreprocessing is performed to find the values of the points where the spectral intensity is the greatest and least, denoted as fmaxAnd fminNormalizing the spectra in order to adapt the algorithm to different spectral curves, i.e.
As shown in FIG. 4, the glucose concentration is 2850 to 2900cm-1The maximum value of the wavelength intensity in the wavenumber range is 15, but the maximum intensity variation range is greatly different for different wavenumber ranges and different substances, so that the normalization process is required, and the result of the normalization process performed on fig. 4 is shown in fig. 5.
(2) Normalized spectral intensity fiThe first derivative and the second derivative of (a),
fi′=(fi+1-fi)/2,
(3) and solving the weighting matrix Q. Q has dimension n, for each spectral point fiFind its neighborhood (f)i-3,fi-2,fi-1,fi,fi+1,fi+2,fi+3) Standard deviation of S. When i is 1, fiIs taken as (f)1,f1,f1,f1,f2,f3,f4) (ii) a When i is 2, fiIs taken as (f)1,f1,f1,f2,f3,f4,f5) (ii) a When i is 3, fiIs taken as (f)1,f1,f2,f3,f4,f5,f6). Taking the maximum value and the minimum value of S as S respectivelymax、SminThen the diagonal element Q of the weighting matrix QiiCan be expressed as:
wherein ln is a natural logarithm, and e is a natural constant, and the value is 2.7.
Let sigmai=qiiThen Q can be expressed as:
the standard deviation S is solved for the spectrum data in fig. 5, the solving result is shown in fig. 6, and it can be seen from comparison between fig. 5 and fig. 6 that the intensity change of the spectrum data is more drastic, the larger the value of the corresponding point S is, the noise is in an oscillation change rule, so the standard deviation S value is smaller, and different weights can be given to the spectrum recovery and the denoising according to S according to the characteristic of the spectrum. Normalizing S to obtain qiiThereby obtaining a matrix Q, the result of which is shown in fig. 7.
(4) Setting an instrument response function H, the instrument response function usually has two selectable functions, a Gaussian function and a Lorentz function, as shown in fig. 8, which are a Gaussian function curve (Gaussian shape) and a Lorentz function curve (Lorentz shape), respectively, and it can be seen that although the shapes of the two functions are different, the half-widths and the heights are the same, and glucose (C) is considered6H12O6) The symmetry of the molecule, a gaussian function was chosen as the kernel function for this example, with the width of the narrowest peak as the width of this gaussian kernel.
(5) And setting a splitting method to solve parameters. Lambda is a regular parameter, the function of lambda lies in balancing a data item and smoothness prior, the value range is (0,1), beta is the regular parameter in the splitting method, the value range is (0,1), k is the iteration number, and d is used for replacing fiThe' is a first derivative of the spectrum, b is an over-denoising compensation value, the value range is (0,1), and the actual value is an empirical value. Initialization of the spectral initial value f0F, the initial value b of the over-noise-removal compensation 00, substitute variable initial value d0=0。
(6) Solving for f using a split Bregman iteration methodk+1And dk+1To f fork+1Linear solution by iterative Gauss-Seidel, HTIs the transpose of matrix H. The delta t is a step of time,energy functional, Df, is the derivation of the matrix.
Wherein the n × n dimensional difference matrix of the spectrumNormalized spectral intensity seti=1,2,…,n。
If the difference is (Df)k+1+bkIf the | is smaller than the threshold lambda/beta, the noise point is considered to be the noise point, and the intensity of the noise point is replaced by 0; if it is larger than the threshold lambda/beta, it is considered as a minutia point of the spectrum, it is smoothness-constrained with a Laplacian-Markov prior. Then-update b, bk+1=bk+((Df)k+1-dk+1) Judgment fk+1Whether or not an iteration stop condition [ (iif) is satisfiedk+1-fk‖/‖fkII > epsilon, epsilon is a threshold value for stopping iteration, and the value range is (10)-9~10-7) If the condition is not met, the loop step (6) is carried out to continue solving, and if the condition is met, the Raman spectrum f is outputk+1。
Setting glucose (C)6H12O6) Process parameters, initialize f first0=g,b0=0,d0When the value is 0, two regular parameters in the splitting method take values respectively as follows: α =0.05 and β = 0.1. And (3) carrying out iterative solution, wherein the intermediate output result of the solution is shown in fig. 9, and as the iteration times are increased, the curve changes less and less, and finally converges when 100 times are carried out.
In order to verify the effectiveness of the algorithm of the present invention, the present invention has been compared with a Tikhonov Regularized Spectrum Deconvolution (TRSD) (see J.Ottaway, J.H.Kalivas, and E.Andries, "Spectral Multivariate Calibration with wavelet reconstruction Using Variants of" applied. Spectroscopy.64, 1388-1395 (2010)), a method of Higher Order Statistics (HOS) (see J.Yuan, Z.Hu, G.Wang, and Z.Xu, "structured high-order statistical breakdown of Spectral data," Chin.Opt.Lett.3,552-555(2005), in short, the name of the method of the present invention is-ordered Spectral data), and the length of the spectrum can be calculated as a ratio of the length of the spectrum to the length of the spectrum, as shown in the longitudinal length chart, and the length of the LM can be calculated as a ratio of the length of the longitudinal spectrum, as shown in the chart of the longitudinal length of the curve, and the length of the LM can be calculated as a ratio of the length of the longitudinal spectrum, as shown in the chart 1; compared with the TRSD algorithm and the HOS algorithm, the LM-QSD algorithm is faster than the TRSD algorithm and the HOS algorithm, and is faster than the TRSD algorithm and the HOS algorithm by more than 4 times, which is mainly benefited by the splitting method introduced by the invention.
TABLE 1 CPU time consumption comparison units(s) for three algorithms
Spectral curve | TRSD | HOS | LM-QSD |
D(+)-Glucopyranose | 15.365s | 17.524s | 4.365s |
methyl | 25.625s | 20.152s | 5.574s |
It will be understood by those skilled in the art that the foregoing is only a preferred embodiment of the present invention, and is not intended to limit the invention, and that any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (3)
1. A spectrum recovery method based on a Laplacian-Markov field is carried out according to the following steps:
(1) for discretized spectrumNormalized processing is carried out to obtain normalized spectral intensity fiI is 1,2,3 … n, n is the total number of spectra;
(2) calculating normalized spectral intensity fiFirst derivative f ofi′;
(3) Calculating each spectral point fiNeighborhood of fi-3,fi-2,fi-1,fi,fi+1,fi+2,fi+3]Standard deviation of (S)iFrom which the maximum standard deviation S is foundmaxAnd minimum standard deviation SminConstructing a weighting matrix Q of n x n dimensions, the diagonal elements of the weighting matrix QSetting the other elements to be 0, wherein ln is a natural logarithm, and e is a natural constant;
(4) and (3) solving the Raman spectrum by adopting a splitting iterative method:
(41) solving parameters by an initialization splitting method: initial value f of the spectrum0F, the initial value b of the over-noise-removal compensation00, substitute variable initial value d00, and 0 is the iteration number k;
(42) updating calculation fk+1And dk+1:
Where max is the maximum value, H is the instrument response function, HTIs the transpose of H, Δ t is the time step, n × n dimensional difference matrix of the spectrumNormalized spectral intensity seti ═ 1,2, …, n, regularization parameters λ ∈ (0,1) and β ∈ (0, 1);
(43) update calculation bk+1=bk+((Df)k+1-dk+1);
(44) Judgment of fk+1Whether or not an iteration stop condition (| f) is satisfiedk+1-fk‖/‖fk|) > epsilon, epsilon is an iteration stop threshold, if not, then return to step (42), otherwise output raman spectrum fk+1。
2. The method for spectrum recovery according to claim 1, wherein the instrument response function H is a gaussian function or a lorentzian function.
3. A spectrum recovery system based on Laplacian-Markov field comprises
A first module for discretizing the spectrumNormalized processing is carried out to obtain normalized spectral intensity fiI is 1,2,3 … n, n is the total number of spectra;
a second module for calculating the normalized spectral intensity fiFirst derivative f ofi′;
A third module for calculating each spectral point fiNeighborhood of [ f ]i-3,fi-2,fi-1,fi,fi+1,fi+2,fi+3]Standard deviation of (S)iFrom which the maximum standard deviation S is foundmaxAnd minimum standard deviation SminConstructing a weighting matrix Q of n x n dimensions, the diagonal elements of the weighting matrix QSetting the other elements to be 0, wherein ln is a natural logarithm, and e is a natural constant;
a fourth module for solving the raman spectrum using a split iteration method, comprising:
the first submodule is used for initializing the splitting method to solve parameters: initial value f of the spectrum0F, the initial value b of the over-noise-removal compensation00, substitute variable initial value d00, and 0 is the iteration number k;
a second submodule for updating the calculation fk+1And dk+1:
Where max is the maximum value, H is the instrument response function, HTIs the transpose of H, Δ t is the time step, n × n dimensional difference matrix of the spectrumNormalized spectral intensity seti ═ 1,2, …, n, regularization parameters λ ∈ (0,1) and β ∈ (0, 1);
a third submodule for updating the calculation bk+1=bk+((Df)k+1-dk+1);
A fourth sub-module for judging fk+1Whether or not iteration is satisfiedStop condition (| f)k+1-fk‖/‖fk|) > epsilon, epsilon is an iteration stop threshold, if not, then return to step (42), otherwise output raman spectrum fk+1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210475610.8A CN103134789B (en) | 2012-11-21 | 2012-11-21 | Spectrum recover method based on Laplacian-Markov field |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210475610.8A CN103134789B (en) | 2012-11-21 | 2012-11-21 | Spectrum recover method based on Laplacian-Markov field |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103134789A true CN103134789A (en) | 2013-06-05 |
CN103134789B CN103134789B (en) | 2015-05-20 |
Family
ID=48494899
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210475610.8A Active CN103134789B (en) | 2012-11-21 | 2012-11-21 | Spectrum recover method based on Laplacian-Markov field |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103134789B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104267018A (en) * | 2014-10-27 | 2015-01-07 | 武汉四方光电科技有限公司 | Method for processing gas concentration signal in Raman gas analyzer |
CN104457986A (en) * | 2014-10-23 | 2015-03-25 | 南京邮电大学 | Spectral resolution enhancing method based on self-adaptation regularization |
CN105354799A (en) * | 2015-09-22 | 2016-02-24 | 华中师范大学 | Blind deconvolution infrared spectrogram super-resolution restoration method |
CN105654447A (en) * | 2016-03-04 | 2016-06-08 | 华中师范大学 | Telemetering spectrum noise suppression method suitable for infrared small target recognition |
CN106991651A (en) * | 2016-12-12 | 2017-07-28 | 南昌大学 | Fast imaging method and system based on synthesis analysis deconvolution network |
CN107067377A (en) * | 2017-03-02 | 2017-08-18 | 中国科学院遥感与数字地球研究所 | The method and device that a kind of shadow Detection of high spectrum image is recovered with spectrum |
CN107966428A (en) * | 2016-10-19 | 2018-04-27 | 西派特(北京)科技有限公司 | A kind of method for improving Miniature Raman spectrometer resolution ratio |
CN109187488A (en) * | 2018-09-30 | 2019-01-11 | 姚志湘 | Raman spectra qualitative for different resolution compares processing method |
CN111104876A (en) * | 2019-12-04 | 2020-05-05 | 南京邮电大学 | Infrared spectrum deconvolution method based on neural network |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002012862A (en) * | 2000-06-29 | 2002-01-15 | Kasei Optonix Co Ltd | Zinc sulfide phosphor and method for evaluating luminance of the phosphor. |
JP2009207697A (en) * | 2008-03-04 | 2009-09-17 | Kao Corp | Method for measuring fiber orientation of hair |
-
2012
- 2012-11-21 CN CN201210475610.8A patent/CN103134789B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002012862A (en) * | 2000-06-29 | 2002-01-15 | Kasei Optonix Co Ltd | Zinc sulfide phosphor and method for evaluating luminance of the phosphor. |
JP2009207697A (en) * | 2008-03-04 | 2009-09-17 | Kao Corp | Method for measuring fiber orientation of hair |
Non-Patent Citations (2)
Title |
---|
杨怀栋等: "盲目反卷积光谱图超分辨复原算法", 《光谱学与光谱分析》 * |
袁玉峰等: "化学计量学结合拉曼光谱在生物材料检测中的应用", 《光谱实验室》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104457986A (en) * | 2014-10-23 | 2015-03-25 | 南京邮电大学 | Spectral resolution enhancing method based on self-adaptation regularization |
CN104267018A (en) * | 2014-10-27 | 2015-01-07 | 武汉四方光电科技有限公司 | Method for processing gas concentration signal in Raman gas analyzer |
CN105354799A (en) * | 2015-09-22 | 2016-02-24 | 华中师范大学 | Blind deconvolution infrared spectrogram super-resolution restoration method |
CN105654447A (en) * | 2016-03-04 | 2016-06-08 | 华中师范大学 | Telemetering spectrum noise suppression method suitable for infrared small target recognition |
CN107966428A (en) * | 2016-10-19 | 2018-04-27 | 西派特(北京)科技有限公司 | A kind of method for improving Miniature Raman spectrometer resolution ratio |
CN106991651A (en) * | 2016-12-12 | 2017-07-28 | 南昌大学 | Fast imaging method and system based on synthesis analysis deconvolution network |
CN106991651B (en) * | 2016-12-12 | 2019-10-29 | 南昌大学 | Fast imaging method and system based on synthesis analysis deconvolution network |
CN107067377A (en) * | 2017-03-02 | 2017-08-18 | 中国科学院遥感与数字地球研究所 | The method and device that a kind of shadow Detection of high spectrum image is recovered with spectrum |
CN107067377B (en) * | 2017-03-02 | 2019-10-29 | 中国科学院遥感与数字地球研究所 | A kind of method and device of the shadow Detection of high spectrum image and spectrum recovery |
CN109187488A (en) * | 2018-09-30 | 2019-01-11 | 姚志湘 | Raman spectra qualitative for different resolution compares processing method |
CN111104876A (en) * | 2019-12-04 | 2020-05-05 | 南京邮电大学 | Infrared spectrum deconvolution method based on neural network |
Also Published As
Publication number | Publication date |
---|---|
CN103134789B (en) | 2015-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103134789B (en) | Spectrum recover method based on Laplacian-Markov field | |
Liu et al. | Efficient blind signal reconstruction with wavelet transforms regularization for educational robot infrared vision sensing | |
Liu et al. | Fast blind instrument function estimation method for industrial infrared spectrometers | |
CN105069758B (en) | A kind of high spectrum image denoising method based on robust low-rank tensor | |
Gholami et al. | A balanced combination of Tikhonov and total variation regularizations for reconstruction of piecewise-smooth signals | |
Liu et al. | A MAP-based algorithm for spectroscopic semi-blind deconvolution | |
CN102184533B (en) | Non-local-restriction-based total variation image deblurring method | |
Liu et al. | NIRS feature extraction based on deep auto-encoder neural network | |
Mohimani et al. | Complex-valued sparse representation based on smoothed ℓ 0 norm | |
Ebadi et al. | A review of applying second-generation wavelets for noise removal from remote sensing data | |
Liu et al. | Multi-order blind deconvolution algorithm with adaptive Tikhonov regularization for infrared spectroscopic data | |
Liu et al. | Blind spectrum reconstruction algorithm with L0-sparse representation | |
Shukla et al. | Generalized fractional filter-based algorithm for image denoising | |
Chen et al. | Hyperspectral imagery denoising using a spatial-spectral domain mixing prior | |
CN105354799A (en) | Blind deconvolution infrared spectrogram super-resolution restoration method | |
Urbas et al. | Automated spectral smoothing with spatially adaptive penalized least squares | |
Cai et al. | Hyperspectral image denoising using adaptive weight graph total variation regularization and low-rank matrix recovery | |
CN109839362B (en) | Infrared spectrum quantitative analysis method based on progressive denoising technology | |
CN116883720A (en) | Fruit and vegetable pesticide residue detection method and system based on spatial spectrum attention network | |
Cohen et al. | Modes of homogeneous gradient flows | |
CN104457986B (en) | A kind of spectrum simulation method based on adaptive regularization | |
Qin et al. | Improved deep residual shrinkage network on near infrared spectroscopy for tobacco qualitative analysis | |
Liu et al. | Estimation of chlorophyll content in maize canopy using wavelet denoising and SVR method | |
CN109143341A (en) | Reduced-rank filtering method based on Hampel norm | |
CN111598795A (en) | Hyperspectral image denoising method and device based on tensor ring decomposition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |