CN108896176A - A kind of Space Consistency bearing calibration of multi-optical spectrum imaging system - Google Patents
A kind of Space Consistency bearing calibration of multi-optical spectrum imaging system Download PDFInfo
- Publication number
- CN108896176A CN108896176A CN201810458025.4A CN201810458025A CN108896176A CN 108896176 A CN108896176 A CN 108896176A CN 201810458025 A CN201810458025 A CN 201810458025A CN 108896176 A CN108896176 A CN 108896176A
- Authority
- CN
- China
- Prior art keywords
- imaging system
- optical spectrum
- spectrum imaging
- response
- reflectivity
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 33
- 238000001228 spectrum Methods 0.000 title claims abstract description 33
- 238000002310 reflectometry Methods 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 18
- 230000003595 spectral effect Effects 0.000 claims description 7
- 239000003086 colorant Substances 0.000 claims description 5
- 230000001373 regressive effect Effects 0.000 claims description 5
- 241001270131 Agaricus moelleri Species 0.000 claims description 4
- 238000005286 illumination Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 210000001525 retina Anatomy 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/2823—Imaging spectrometer
- G01J2003/2826—Multispectral imaging, e.g. filter imaging
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The invention discloses a kind of Space Consistency bearing calibrations of multi-optical spectrum imaging system.Method includes：Serial grayscale paper is measured in the response of multi-optical spectrum imaging system objective table different location, the color lump response rebuild in multi-optical spectrum imaging system for reflectivity is measured simultaneously, with polynomial regression model, different grayscale papers are mapped in the response of different location using location information as independent variable to the response of center, while the response for being used for the color lump of reflectivity reconstruction is mapped to its target value using location information as independent variable.The present invention has combined Space Consistency accuracy requirement and reflectivity reconstruction precision demand in multi-optical spectrum imaging system, so that after Space Consistency corrects, the Space Consistency color difference of system greatly reduces, while the color difference typical value that reflectivity is rebuild also is reduced.
Description
Technical field
The present invention relates to the Space Consistency correction sides of multi-optical spectrum imaging technology more particularly to a kind of multi-optical spectrum imaging system
Method can make reflectivity rebuild color difference typical value and also decrease in the case where system space consistency color difference is greatly reduced.
Background technique
When on a branch of radiation of visible light to different colours object, degree of reflection of the object to different-waveband visible light
Difference, so that the stimulation that human eye retina is subject to has differences, to form human eye to the subjective feeling of different colours.Therefore,
The essence of object color measurement is the spectral reflectivity for measuring object.
Multi-optical spectrum imaging system is a kind of equipment that can measure object color and reflectivity, its working principle is that by uniform
Illumination be irradiated on the object of objective table, the visible light of different-waveband is filtered by light-splitting device, and object is acquired by camera
The response of different-waveband in body reflected light, and then collected response is rebuild to obtain object spectra reflection by reflectivity
Rate.The response referred in the present invention, the number that the not reflected rate as directly measured in multi-optical spectrum imaging system is rebuild
According to.
However the illumination due to being irradiated on object may be uneven, and in the optical filter of colour wheel and the camera lens of camera
The heart and the light transmittance of surrounding are not exactly the same, so measuring when the sample of same color is placed on objective table different location
Color may be different.
Therefore multi-optical spectrum imaging system needs to carry out Space Consistency correction, so that same color exists after corrected
The color difference for the color that different location measurement obtains is as small as possible.
A kind of method of common Space Consistency correction is the Space Consistency bearing calibration based on interpolation virtual whiteboard,
Method process is as follows：Take the grayscale paper of N number of different gray scales that it is limited to place it in system stage for each grayscale paper
Limited response for shooting of M position, one piece of virtual whiteboard is obtained using interpolation method according to limited response,
N block virtual whiteboard can be generated according to the difference of gray value, and the virtual mean value of N block is drawn high to 1.For multi-optical spectrum imaging system
Obtained response is shot, finds one piece of nearest virtual whiteboard of its response distance, and draw high to after 1 according to divided by the blank
The value of corresponding position obtains corrected response, thus completes the system space consistency school based on interpolation virtual whiteboard
Just.Although this method can improve the Space Consistency of system to a certain extent, by application method to multispectral
After imaging system is corrected, system space consistency promoted it is still limited, and this method will affect in degree together it is anti-
It penetrates rate and rebuilds color difference, do not ensure that reflectivity is rebuild color difference and reduced.
Summary of the invention
In view of the above-mentioned deficiencies in the prior art, it is an object of the present invention to provide a kind of Space Consistency of multi-optical spectrum imaging system
Bearing calibration, this method can be such that the Space Consistency color difference of system greatly reduces, while reflectivity being made to rebuild color difference typical value
Also reduced.
The purpose of the present invention is achieved through the following technical solutions：A kind of Space Consistency of multi-optical spectrum imaging system
Bearing calibration, the method includes the steps of：
(1) the response R for the color lump rebuild with multi-optical spectrum imaging system shooting for reflectivity1, while according to reflectivity
The inverse process of reconstruction is instead released from the standard reflectivity of these color lumps can make the smallest response S of reflectivity reconstruction error1, S1
As R1Regressive object value；
(2) grayscale paper for selecting N number of different gray scales, is sequentially placed multi-optical spectrum imaging system objective table for each grayscale paper
M different location, measure the grayscale paper respectively in the response R of M different location2, the response work of an optional position
For R2Regressive object value, be denoted as S2；
(3) regression model is established：Polynomial regression, regression model are individually carried out to each channel of multi-optical spectrum imaging system
For：
S=f (r, x, y)=a0+a1r+a2x+a3y+a4x2+a5y2+a6xy+…+akxmyn (1)
Wherein m and n is natural number, akFor xmynParameter before, r are the response in certain channel, and x, y are grayscale paper or color
The position of block, s is the regressive object value of grayscale paper or color lump, according to step (1) and step (2) ready R1,S1And R2,S2
Data are trained regression model, calculate the parameters a in regression modeli；
(4) arbitrary sample is returned according to its own response r and its position coordinates x, y according to what step (3) were established
Model is returned to calculate its regressand value s, to realize that Space Consistency corrects.
Further, the grayscale paper of N number of difference gray scale used in the step (2) can use the color lump of N number of different colours
To substitute.
Further, in the step (3), regression model is：
S=f (r, x, y)=a0+a1r+a2x+a3y+a4x2+a5y2+a6xy (2)
Further, the color lump shape is rectangle or circle, by the spectral reflectivity mean value generation for measuring color block areas
The spectral reflectivity of table color lump.
Further, the regressor r in the regression model, x, y are required to normalize between 0-1.
The beneficial effects of the invention are as follows：Although common Space Consistency bearing calibration can reduce more to a certain extent
The Space Consistency color difference of spectrum imaging system, but the reflectivity that will affect system rebuilds color difference, does not ensure that reflectivity
Color difference is rebuild to reduce.The color lump for being used for reflectivity reconstruction is innovatively introduced free-air correction by the method for the present invention, uses enough complexity
Polynomial regression model, grayscale paper is not only corrected to same response in the different responses of different location, and will use
In reflectivity rebuild color lump according to position information correction to can make reflectivity reconstruction the smaller response of color difference.Reach same
When reduce the effect that multi-optical spectrum imaging system Space Consistency color difference and reflectivity rebuild color difference.
Detailed description of the invention
Fig. 1 is that reflectivity rebuilds color used in multi-optical spectrum imaging system of the present invention；
Fig. 2 is the present invention for placing the plate schematic diagram with M groove of color lump, M=9 in the figure；
Fig. 3 is to change about position x, the corrected value of y with x, y in calibration model after being corrected with the method for the present invention
Surface chart, each point around mean camber are corrected value of the every training data in different location.
Specific embodiment
The specific embodiment of the invention is described further with reference to the accompanying drawing.
A kind of Space Consistency bearing calibration of multi-optical spectrum imaging system provided by the invention, comprises the steps of：
(1) the response R for the color lump rebuild as shown in Figure 1 for reflectivity with multi-optical spectrum imaging system shooting1, while root
The smallest sound of reflectivity reconstruction error can be made by instead releasing according to the inverse process that reflectivity is rebuild from the standard reflectivity of these color lumps
It should value S1, S1As R1Regressive object value；
(2) grayscale paper for selecting N number of different gray scales, by each grayscale paper be sequentially placed it is as shown in Figure 2 it is multispectral at
As M different location of system stage, M position distribution is uniformly and most wide as far as possible, measures the grayscale paper respectively in M difference
The response R of position2, the response of an optional position is as R2Regressive object value, be denoted as S2；
(3) regression model is established：Polynomial regression, regression model are individually carried out to each channel of multi-optical spectrum imaging system
For：
S=f (r, x, y)=a0+a1r+a2x+a3y+a4x2+a5y2+a6xy+…+akxmyn (1)
Wherein m and n is natural number, akFor xmynParameter before, r are the response in certain channel, and x, y are grayscale paper or color
The position of block, s is the regressive object value of grayscale paper or color lump, according to step (1) and step (2) ready R1,S1And R2,S2
Data are trained regression model, calculate the parameters a in regression modeli；
(4) arbitrary sample is returned according to its own response r and its position coordinates x, y according to what step (3) were established
Model is returned to calculate its regressand value s, to realize that Space Consistency corrects.
Further, the grayscale paper of N number of difference gray scale used in the step (2) can use the color lump of N number of different colours
To substitute.
Further, in the step (3), regression model be can simplify as following form：
S=f (r, x, y)=a0+a1r+a2x+a3y+a4x2+a5y2+a6xy (2)
Further, the color lump shape is rectangle or circle, by the spectral reflectivity mean value generation for measuring color block areas
The spectral reflectivity of table color lump.
Further, the regressor r in the regression model, x, y are required to normalize between 0-1.
Embodiment 1
Fig. 3 be using the method for the present invention to multi-optical spectrum imaging system carry out Space Consistency after, in calibration model (1) about
The corrected value of position x, y are with x, the surface chart of y variation, and each point around mean camber is every training data in different location
When corrected value.
Table 1 is to utilize Space Consistency bearing calibration proposed by the present invention and the Space Consistency based on interpolation virtual whiteboard
The correction Comparative result table of bearing calibration.It can be found that no matter Space Consistency bearing calibration proposed by the present invention is in space one
Cause property color difference result still rebuilds in color difference result in reflectivity and is all substantially better than the Space Consistency based on interpolation virtual whiteboard
Bearing calibration.
Table 1
The above is only the specific embodiment of the invention, cannot be limited the scope of the invention with this, in the art
Those skilled in the art change according to known to equivalent change made by this creation and those skilled in that art, all should still belong to
The range that the present invention covers.
Claims (5)
1. a kind of Space Consistency bearing calibration of multi-optical spectrum imaging system, which is characterized in that comprise the steps of：
(1) the response R for the color lump rebuild with multi-optical spectrum imaging system shooting for reflectivity1, while being rebuild according to reflectivity
Inverse process is instead released from the standard reflectivity of these color lumps can make the smallest response S of reflectivity reconstruction error1, S1As R1
Regressive object value；
(2) each grayscale paper, is sequentially placed the M of multi-optical spectrum imaging system objective table by the grayscale paper for selecting N number of different gray scales
A different location measures the grayscale paper in the response R of M different location respectively2, the response of an optional position is as R2
Regressive object value, be denoted as S2；
(3) regression model is established：Polynomial regression is individually carried out to each channel of multi-optical spectrum imaging system, regression model is：
S=f (r, x, y)=a0+a1r+a2x+a3y+a4x2+a5y2+a6xy+…+akxmyn (1)
Wherein m and n is natural number, akFor xmynParameter before, r are the response in certain channel, and x, y are grayscale paper or color lump
Position, s is the regressive object value of grayscale paper or color lump, according to step (1) and step (2) ready R1,S1And R2,S2Data
Regression model is trained, the parameters a in regression model is calculatedi；
(4) for arbitrary sample, according to its own response r and its position coordinates x, y, the recurrence mould established according to step (3)
Type calculates its regressand value s, to realize that Space Consistency corrects.
2. the Space Consistency bearing calibration of multi-optical spectrum imaging system according to claim 1, it is characterised in that：The step
Suddenly the grayscale paper of N number of difference gray scale used in (2) can be substituted with the color lump of N number of different colours.
3. the Space Consistency bearing calibration of multi-optical spectrum imaging system according to claim 1, it is characterised in that：The step
Suddenly in (3), regression model is：
S=f (r, x, y)=a0+a1r+a2x+a3y+a4x2+a5y2+a6xy (2)
4. the Space Consistency bearing calibration of multi-optical spectrum imaging system according to claim 1, it is characterised in that：The color
Block-shaped is rectangle or circle, and the spectral reflectivity mean value by measuring color block areas represents the spectral reflectivity of color lump.
5. the Space Consistency bearing calibration of multi-optical spectrum imaging system according to claim 1, it is characterised in that：Described time
Return the regressor r in model, x, y are required to normalize between 0-1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810458025.4A CN108896176B (en) | 2018-05-14 | 2018-05-14 | A kind of Space Consistency bearing calibration of multi-optical spectrum imaging system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810458025.4A CN108896176B (en) | 2018-05-14 | 2018-05-14 | A kind of Space Consistency bearing calibration of multi-optical spectrum imaging system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108896176A true CN108896176A (en) | 2018-11-27 |
| CN108896176B CN108896176B (en) | 2019-10-11 |
Family
ID=64343882
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810458025.4A Active CN108896176B (en) | 2018-05-14 | 2018-05-14 | A kind of Space Consistency bearing calibration of multi-optical spectrum imaging system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108896176B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109738069A (en) * | 2018-12-27 | 2019-05-10 | 浙江农林大学暨阳学院 | The method that multispectral imaging illuminates spatial heterogeneity correction |
| CN109738068A (en) * | 2018-12-25 | 2019-05-10 | 浙江农林大学暨阳学院 | A kind of correction multispectral camera response non-linear method |
| CN110017816A (en) * | 2019-04-28 | 2019-07-16 | 新疆师范大学 | It is a kind of to take photo by plane the image flame detection device of measurement suitable for low-to-medium altitude aircraft |
| CN110660112A (en) * | 2019-09-29 | 2020-01-07 | 浙江大学 | Drawing spectrum reconstruction method based on special color card and multispectral imaging |
| CN111047539A (en) * | 2019-12-27 | 2020-04-21 | 上海工程技术大学 | Fabric image color calibration algorithm based on spectral reflectivity reconstruction |
| CN113947553A (en) * | 2021-12-20 | 2022-01-18 | 山东信通电子股份有限公司 | Image brightness enhancement method and device |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1902471A (en) * | 2003-11-07 | 2007-01-24 | 帝国学院创新有限公司 | Method of spectroscopy |
| CN105067116A (en) * | 2015-07-15 | 2015-11-18 | 北京农业信息技术研究中心 | Method and system for splicing frame imaging spectral data |
| CN106458586A (en) * | 2014-06-17 | 2017-02-22 | Ay实验室有限公司 | Method for quantifying the amount of ammonium bicarbonate in a solid sample of ammonium carbamate |
| CN107271372A (en) * | 2017-06-02 | 2017-10-20 | 西北农林科技大学 | A kind of Apple Leaves chlorophyll remote sensing estimation method |
| CN107505046A (en) * | 2017-07-31 | 2017-12-22 | 征图新视（江苏）科技有限公司 | Multispectral color measurement method |
-
2018
- 2018-05-14 CN CN201810458025.4A patent/CN108896176B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1902471A (en) * | 2003-11-07 | 2007-01-24 | 帝国学院创新有限公司 | Method of spectroscopy |
| CN106458586A (en) * | 2014-06-17 | 2017-02-22 | Ay实验室有限公司 | Method for quantifying the amount of ammonium bicarbonate in a solid sample of ammonium carbamate |
| CN105067116A (en) * | 2015-07-15 | 2015-11-18 | 北京农业信息技术研究中心 | Method and system for splicing frame imaging spectral data |
| CN107271372A (en) * | 2017-06-02 | 2017-10-20 | 西北农林科技大学 | A kind of Apple Leaves chlorophyll remote sensing estimation method |
| CN107505046A (en) * | 2017-07-31 | 2017-12-22 | 征图新视（江苏）科技有限公司 | Multispectral color measurement method |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109738068A (en) * | 2018-12-25 | 2019-05-10 | 浙江农林大学暨阳学院 | A kind of correction multispectral camera response non-linear method |
| CN109738068B (en) * | 2018-12-25 | 2021-06-04 | 浙江农林大学暨阳学院 | Method for correcting nonlinearity of response value of multispectral camera |
| CN109738069A (en) * | 2018-12-27 | 2019-05-10 | 浙江农林大学暨阳学院 | The method that multispectral imaging illuminates spatial heterogeneity correction |
| CN110017816A (en) * | 2019-04-28 | 2019-07-16 | 新疆师范大学 | It is a kind of to take photo by plane the image flame detection device of measurement suitable for low-to-medium altitude aircraft |
| CN110660112A (en) * | 2019-09-29 | 2020-01-07 | 浙江大学 | Drawing spectrum reconstruction method based on special color card and multispectral imaging |
| CN111047539A (en) * | 2019-12-27 | 2020-04-21 | 上海工程技术大学 | Fabric image color calibration algorithm based on spectral reflectivity reconstruction |
| CN113947553A (en) * | 2021-12-20 | 2022-01-18 | 山东信通电子股份有限公司 | Image brightness enhancement method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108896176B (en) | 2019-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108896176B (en) | A kind of Space Consistency bearing calibration of multi-optical spectrum imaging system | |
| CN110392252B (en) | Method for generating correction model of camera to correct aberration | |
| CN111351446B (en) | Light field camera calibration method for three-dimensional topography measurement | |
| CN112669394B (en) | Automatic calibration method for vision detection system | |
| CN109827502B (en) | High-precision calibration method for line-structured light vision sensor for calibration point image compensation | |
| CN105913401B (en) | Industrial camera photography measurement image luminance compensation method | |
| CN109727290B (en) | Zoom camera dynamic calibration method based on monocular vision triangulation distance measurement method | |
| CN106705855A (en) | High-dynamic performance three-dimensional measurement method based on adaptive grating projection | |
| CN111009014A (en) | Calibration method of orthogonal spectral imaging pose sensor of general imaging model | |
| CN105466889B (en) | The acquisition method of complex organization's body surface face illuminance in a kind of spatial frequency domain imaging | |
| CN110660112B (en) | Drawing spectrum reconstruction method based on special color card and multispectral imaging | |
| CN104296968A (en) | Modulation transfer function test method of multichannel CCD | |
| CN108917640A (en) | A kind of laser blind hole depth detection method and its system | |
| CN112848281A (en) | Light compensation method for photocuring 3D printer | |
| CN117053852B (en) | Optical measurement correction method and system for underwater structure | |
| CN109253862A (en) | A kind of colour measurement method neural network based | |
| CN111047586A (en) | Pixel equivalent measuring method based on machine vision | |
| CN112734824A (en) | Three-dimensional reconstruction method based on generalized luminosity stereo model | |
| CN109323851B (en) | Terahertz focal plane response rate and response rate nonuniformity test system and method | |
| CN108226059A (en) | A kind of satellite EO-1 hyperion CO2The in-orbit Calibration Method of survey meter | |
| CN108765301B (en) | Optical system and distortion correction method and system thereof | |
| CN102419861A (en) | Color image correcting method based on topology subdivision of uniform color space | |
| CN109238461A (en) | Room lighting dazzle measurement method based on digital camera images | |
| CN113936580B (en) | Display screen brightness measurement model generation method and related device | |
| CN110095865A (en) | A kind of uniform irradiation design method of LED structure light source |
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 |