CN108225570A - A kind of adaptive non-uniformity correction algorithm in short-wave infrared focal plane - Google Patents
A kind of adaptive non-uniformity correction algorithm in short-wave infrared focal plane Download PDFInfo
- Publication number
- CN108225570A CN108225570A CN201711497395.0A CN201711497395A CN108225570A CN 108225570 A CN108225570 A CN 108225570A CN 201711497395 A CN201711497395 A CN 201711497395A CN 108225570 A CN108225570 A CN 108225570A
- Authority
- CN
- China
- Prior art keywords
- focal plane
- heterogeneity
- coefficient
- correction
- short
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 30
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 11
- 239000000523 sample Substances 0.000 claims description 13
- 230000005855 radiation Effects 0.000 claims description 6
- 239000000654 additive Substances 0.000 claims description 3
- 230000000996 additive effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims 1
- 238000003491 array Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 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
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
-
- 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
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/80—Calibration
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
- Radiation Pyrometers (AREA)
Abstract
The invention discloses a kind of adaptive non-uniformity correction algorithms in short-wave infrared focal plane, the adaptive non-uniformity correction algorithm in short-wave infrared focal plane carries out regression estimates using the heterogeneity of logistics function combination short-wave infrared focal plane asymmetric model focal planes and Nonuniformity Correction coefficient is obtained, and achievees the purpose that correction.
Description
Technical field
The present invention relates to the technical fields of short-wave infrared focal plane, adaptive non-equal more particularly to short-wave infrared focal plane
Even correcting algorithm.
Background technology
The heterogeneity of short-wave infrared focal plane arrays (FPA) can be with use environment temperature, imageable target temperature, focal plane itself
Temperature, focal plane working time length variation and change, correction coefficient needs real-time update.
It is above-mentioned existing therefore, it is desirable to there is a kind of adaptive non-uniformity correction algorithm in short-wave infrared focal plane that can effectively solve
There is the defects of technology.
Invention content
The purpose of the present invention is to provide a kind of short-wave infrared focal planes based on logistics and multiple-factor calibration model
Adaptive non-uniformity correction algorithm uses logistics function combination short-wave infrared focal plane asymmetric model focal planes
Heterogeneity carry out regression estimates Nonuniformity Correction coefficient is obtained, achieve the purpose that correction.
To achieve the above object, the present invention provides a kind of short-wave infrared based on logistics and multiple-factor calibration model
The adaptive non-uniformity correction algorithm in focal plane includes the following steps：
1. two hypotheses are carried out to linear model；
2. ignoring probe unit, there are under error condition, approx think burnt flat in the smaller saturation with removal of radiation scope
Face response is linear, and under the conditions of uniform radiancy, the heterogeneity linear model of focal plane can be represented by the following formula：
Yk(i, j)=ak(i,j)*Xk(i,j)-bk(i,j)+nk(i,j)
In formula, XkWhen (i, j) represents kth frame, the true infra-red radiation degree that detects member received by of the coordinate for (i, j), ak
(i, j) and bk(i, j) is respectively at this point, the multiplicative noise variable of the probe unit and additive noise variable, nk(i, j) is reads
The electronic noise of circuit, ak(i, j) and bkThe value of (i, j) is also changed over time and is changed simultaneously；
3. the heterogeneity of focal plane 2. can be obtained mainly by the gain of multiplying property and the biasing group of additivity by step
Into influence focal plane asymmetric is mainly to bias, and with environment temperature, detector temperature and the change of working time
Change, be biased in and ceaselessly change, therefore many scene correcting algorithms occur, only to biasing ceaselessly predictive compensation, in scene
Constantly good effect is can be obtained by the case of movement；
4. the combination for the heterogeneity of bias being resolved into a few class impact factors of environment temperature, detector temperature adds one
A fixed mean value biasing forms model：
Yk(i, j)=Ak(i,j)*Xk(i,j)+Bk(i,j)*TFPA+Ck(i,j)*TENV+m
In formula, Ak(i, j) and Bk(i,j)、Ck(i, j) is respectively that focal plane asymmetric response coefficient, detector temperature are non-
Uniformity coefficient, working time heterogeneity coefficient, m are fixed mean value biasings, can be ignored in correction；
5. gray scale is continuous between assuming the pixel of infrared focal plane imaging, then in infrared image between adjacent picture elements
Difference is very little in the statistical significance of certain frame number, it means that two adjacent picture elements are almost on the histogram of time
Equal；
6. the true response approximation of single probe unit is equal to according to the hypothesis of step 5. by neighborhood detection member response
Mean value Y, so as to the estimated value returned as logistics：
Then the logistics function models of nonuniformity correction estimated value are
7. according to above-mentioned nonuniformity correction model, the logistics function models of nonuniformity correction are：
8. using maximal possibility estimation, the parameter of Nonuniformity Correction, such as following formula is obtained：
Solve Ak(i, j) and Bk(i,j)、Ck(i, j) focal plane asymmetric response coefficient, detector temperature heterogeneity
Coefficient, working time heterogeneity coefficient, specific solve can use newton-La Feisen alternative manners to solve.
Preferably, two hypotheses of the linear model include：(1) response of each probe unit exists in focal plane
It is stable on time, will not changes because of working time length；(2) working region of each probe unit is linear.
For the deficiency that this non-homogeneous feature of short-wave infrared focal plane and shutter correct, the present invention provides one kind to be based on
The adaptive non-uniformity correction algorithm in short-wave infrared focal plane of logistics and multiple-factor calibration model, the algorithm Comprehensive
The considerations of arrived each heteropical impact factor of environment temperature, explorer response, detector temperature, so as to solving well
The uncertainty of shortwave focal plane asymmetric correction, quantitative test the result shows that, this method can reach 0.001% it is non-
Even property correction is remaining, significantly improves correction accuracy, improves the picture quality of infrared imaging, expands being applicable in for alignment technique
Range.
Specific embodiment
The adaptive non-uniformity correction algorithm in short-wave infrared focal plane based on logistics and multiple-factor calibration model includes
Following steps：
1. two hypotheses of traditional linear model：First, the response of each probe unit is in time in focal plane
It is stable, will not changes because of working time length；Second is that the working region of each probe unit is linear；
2. ignoring probe unit, there are under error condition, approx think burnt flat in the smaller saturation with removal of radiation scope
Face response is linear, and under the conditions of uniform radiancy, the heterogeneity linear model of focal plane can be represented by the following formula：
Yk(i, j)=ak(i,j)*Xk(i,j)-bk(i,j)+nk(i,j)
In formula, XkWhen (i, j) represents kth frame, the true infra-red radiation degree that detects member received by of the coordinate for (i, j), ak
(i, j) and bk(i, j) is respectively at this point, the multiplicative noise variable of the probe unit and additive noise variable, nk(i, j) is reads
The electronic noise of circuit, ak(i, j) and bkThe value of (i, j) is also changed over time and is changed simultaneously；
3. the heterogeneity of focal plane 2. can be obtained mainly by the gain of multiplying property and the biasing group of additivity by step
Into influence focal plane asymmetric is mainly to bias, and with environment temperature, detector temperature and the change of working time
Change, be biased in and ceaselessly change, therefore many scene correcting algorithms occur, only to biasing ceaselessly predictive compensation, in scene
Constantly good effect is can be obtained by the case of movement；
4. the combination for the heterogeneity of bias being resolved into a few class impact factors of environment temperature, detector temperature adds one
A fixed mean value biasing forms model：
Yk(i, j)=Ak(i,j)*Xk(i,j)+Bk(i,j)*TFPA+Ck(i,j)*TENV+m
In formula, Ak(i, j) and Bk(i,j)、Ck(i, j) is respectively that focal plane asymmetric response coefficient, detector temperature are non-
Uniformity coefficient, working time heterogeneity coefficient, m are fixed mean value biasings, can be ignored in correction；
5. gray scale is continuous between assuming the pixel of infrared focal plane imaging, then in infrared image between adjacent picture elements
Difference is very little in the statistical significance of certain frame number, it means that two adjacent picture elements are almost on the histogram of time
Equal；
6. the true response approximation of single probe unit is equal to according to the hypothesis of step 5. by neighborhood detection member response
Mean value Y, so as to the estimated value returned as logistics：
Then the logistics function models of nonuniformity correction estimated value are
7. according to above-mentioned nonuniformity correction model, the logistics function models of nonuniformity correction are：
8. using maximal possibility estimation, the parameter of Nonuniformity Correction, such as following formula is obtained：
Solve Ak(i, j) and Bk(i,j)、Ck(i, j) focal plane asymmetric response coefficient, detector temperature heterogeneity
Coefficient, working time heterogeneity coefficient, specific solve can use newton-La Feisen alternative manners to solve.
In practical application, with gradient ascent algorithm, iteration continues always.Gradient ascent algorithm is compared with Newton iteration,
Convergence rate is slow, because gradient ascent algorithm is single order convergence, and Newton iteration belongs to second order convergence.
Above-mentioned algorithm practical application, the logistics estimation windows of correction are 3 × 3 statistical window, are risen by gradient
Method, generally converges to optimal in 100 frame image post-equalization coefficients, and image is more uniform after correction, and calibration result is more apparent, and should
Algorithm is not in exaggerated correction when scene stillness or slow movement, leads to " ghost due to using multivariate calibration model
Shadow " without scene motion, calculates simply convenient for real-time parallel processing, breaches the practicability limit of adaptive nonuniformity correction
System
It is last it is to be noted that：The above embodiments are merely illustrative of the technical solutions of the present invention, rather than its limitations.To the greatest extent
Pipe is with reference to the foregoing embodiments described in detail the present invention, it will be understood by those of ordinary skill in the art that：It is still
It can modify to the technical solution recorded in foregoing embodiments or which part technical characteristic is equally replaced
It changes；And these modifications or replacement, the essence for various embodiments of the present invention technical solution that it does not separate the essence of the corresponding technical solution
God and range.
Claims (2)
1. a kind of adaptive non-uniformity correction algorithm in short-wave infrared focal plane, which is characterized in that include the following steps：
1. two hypotheses are carried out to linear model；
2. ignoring probe unit, there are under error condition, approx think focal plane sound in the smaller saturation with removal of radiation scope
Should be linear, under the conditions of uniform radiancy, the heterogeneity linear model of focal plane can be represented by the following formula：
Yk(i, j)=ak(i,j)*Xk(i,j)-bk(i,j)+nk(i,j)
In formula, XkWhen (i, j) represents kth frame, the true infra-red radiation degree that detects member received by of the coordinate for (i, j), ak(i,
And b j)k(i, j) is respectively at this point, the multiplicative noise variable of the probe unit and additive noise variable, nk(i, j) is reading circuit
Electronic noise, ak(i, j) and bkThe value of (i, j) is also changed over time and is changed simultaneously；
3. it 2. can show that the heterogeneity of focal plane is mainly made of the gain of multiplying property and the biasing of additivity by step, shadow
Mainly biasing for focal plane asymmetric is rung, and with environment temperature, detector temperature and the variation of working time, biasing
Ceaselessly changing, therefore many scene correcting algorithms occurring, only to biasing ceaselessly predictive compensation, constantly moved in scene
In the case of can be obtained by good effect；
4. the combination for the heterogeneity of bias being resolved into a few class impact factors of environment temperature, detector temperature is solid plus one
Fixed mean value biasing forms model：
Yk(i, j)=Ak(i,j)*Xk(i,j)+Bk(i,j)*TFPA+Ck(i,j)*TENV+m
In formula, Ak(i, j) and Bk(i,j)、Ck(i, j) is respectively that focal plane asymmetric response coefficient, detector temperature are non-homogeneous
Property coefficient, working time heterogeneity coefficient, m are fixed mean value biasings, can be ignored in correction；
5. gray scale is continuous between assuming the pixel of infrared focal plane imaging, then the difference in infrared image between adjacent picture elements
It is very little in the statistical significance of certain frame number, it means that two adjacent picture elements are almost equal on the histogram of time
's；
6. the true response approximation of single probe unit is equal to according to the hypothesis of step 5. by neighborhood detection member response mean value
Y, so as to the estimated value returned as logistics：
Then the logistics function models of nonuniformity correction estimated value are
7. according to above-mentioned nonuniformity correction model, the logistics function models of nonuniformity correction are：
8. using maximal possibility estimation, the parameter of Nonuniformity Correction, such as following formula is obtained：
Solve Ak(i, j) and Bk(i,j)、Ck(i, j) focal plane asymmetric response coefficient, detector temperature heterogeneity coefficient,
Working time heterogeneity coefficient, specific solve can use newton-La Feisen alternative manners to solve.
2. the adaptive non-uniformity correction algorithm in short-wave infrared focal plane as described in claim 1, it is characterised in that：It is described linear
Two hypotheses of model include：(1) in focal plane the response of each probe unit be in time it is stable, will not be because of
Working time length and change；(2) working region of each probe unit is linear.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711497395.0A CN108225570B (en) | 2017-12-31 | 2017-12-31 | Short wave infrared focal plane self-adaptive non-uniform correction algorithm |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711497395.0A CN108225570B (en) | 2017-12-31 | 2017-12-31 | Short wave infrared focal plane self-adaptive non-uniform correction algorithm |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108225570A true CN108225570A (en) | 2018-06-29 |
| CN108225570B CN108225570B (en) | 2020-01-14 |
Family
ID=62644960
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711497395.0A Active CN108225570B (en) | 2017-12-31 | 2017-12-31 | Short wave infrared focal plane self-adaptive non-uniform correction algorithm |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108225570B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111678605A (en) * | 2020-06-17 | 2020-09-18 | 北京中云微迅信息技术有限公司 | Infrared temperature measurement compensation method based on thermopile electricity |
| CN111998961A (en) * | 2020-09-23 | 2020-11-27 | 国科天成（北京）科技有限公司 | Infrared focal plane non-uniformity correction method based on pixel point temperature drift estimation |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030122077A1 (en) * | 2001-05-11 | 2003-07-03 | Kaufman Charles S. | Method and apparatus for temperature compensation of an uncooled focal plane array |
| US20060231760A1 (en) * | 2005-04-15 | 2006-10-19 | Korea Advanced Institute Of Science And Technology | Compensation circuit for compensating non-uniformity according to change of operating temperature of bolometer |
| CN101776486A (en) * | 2009-12-31 | 2010-07-14 | 华中科技大学 | Method for correcting non-uniformity fingerprint pattern on basis of infrared focal plane |
| CN102855610A (en) * | 2012-08-03 | 2013-01-02 | 南京理工大学 | Method for correcting infrared image heterogeneity by using parameter correctness factor |
| KR101282567B1 (en) * | 2013-04-25 | 2013-07-04 | 주식회사 콕스 | Apparatus and method for amending non-uniformity of handheld thermal imaging camera |
| CN103604503A (en) * | 2013-06-17 | 2014-02-26 | 中国航天科工集团第三研究院第八三五八研究所 | High-precision non-uniformity correction method for dynamic adjustment of integration time |
| CN104268870A (en) * | 2014-09-24 | 2015-01-07 | 北京津同利华科技有限公司 | Short-wave infrared focal plane non-uniformity correction algorithm based on wavelet transformation histogram |
-
2017
- 2017-12-31 CN CN201711497395.0A patent/CN108225570B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030122077A1 (en) * | 2001-05-11 | 2003-07-03 | Kaufman Charles S. | Method and apparatus for temperature compensation of an uncooled focal plane array |
| US20060231760A1 (en) * | 2005-04-15 | 2006-10-19 | Korea Advanced Institute Of Science And Technology | Compensation circuit for compensating non-uniformity according to change of operating temperature of bolometer |
| CN101776486A (en) * | 2009-12-31 | 2010-07-14 | 华中科技大学 | Method for correcting non-uniformity fingerprint pattern on basis of infrared focal plane |
| CN102855610A (en) * | 2012-08-03 | 2013-01-02 | 南京理工大学 | Method for correcting infrared image heterogeneity by using parameter correctness factor |
| KR101282567B1 (en) * | 2013-04-25 | 2013-07-04 | 주식회사 콕스 | Apparatus and method for amending non-uniformity of handheld thermal imaging camera |
| CN103604503A (en) * | 2013-06-17 | 2014-02-26 | 中国航天科工集团第三研究院第八三五八研究所 | High-precision non-uniformity correction method for dynamic adjustment of integration time |
| CN104268870A (en) * | 2014-09-24 | 2015-01-07 | 北京津同利华科技有限公司 | Short-wave infrared focal plane non-uniformity correction algorithm based on wavelet transformation histogram |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111678605A (en) * | 2020-06-17 | 2020-09-18 | 北京中云微迅信息技术有限公司 | Infrared temperature measurement compensation method based on thermopile electricity |
| CN111998961A (en) * | 2020-09-23 | 2020-11-27 | 国科天成（北京）科技有限公司 | Infrared focal plane non-uniformity correction method based on pixel point temperature drift estimation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108225570B (en) | 2020-01-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108846805B (en) | Infrared thermal image two-point non-uniform correction method based on scene self-adaption | |
| CN109272520B (en) | Self-adaptive infrared focal plane non-uniform correction method combining motion guidance and edge detection | |
| CN107255521B (en) | A kind of Infrared Image Non-uniformity Correction method and system | |
| CN106017695B (en) | Adaptive infrared asymmetric correction method based on state estimation | |
| CN109813442B (en) | Multi-frame processing-based internal stray radiation non-uniformity correction method | |
| US10803557B2 (en) | Non-uniformity correction method for infrared image based on guided filtering and high-pass filtering | |
| CN106197673B (en) | A kind of adaptive wide temperature range non-uniform correction method and system | |
| WO2013084782A1 (en) | Image processing device | |
| CN109741267A (en) | Infrared Image Non-uniformity Correction method based on three sides filtering and neural network | |
| CN102968776B (en) | Linear filter and non-linear filter combined heterogeneity correction method | |
| CN108225570A (en) | A kind of adaptive non-uniformity correction algorithm in short-wave infrared focal plane | |
| CN107271046B (en) | A kind of asymmetric correction method and device based on partial statistics correlation | |
| CN108665425A (en) | Infrared Image Non-uniformity Correction method based on interframe registration and adaptive step | |
| WO2014173464A1 (en) | Noise correction in an infra-red imaging system | |
| CN110033414A (en) | A kind of Infrared Image Non-uniformity Correction method and system based on equalization processing | |
| US11195254B2 (en) | Interframe registration and adaptive step size-based non-uniformity correction method for infrared image | |
| EP3785220A1 (en) | A method for confident registration-based non-uniformity correction using spatio-temporal update mask | |
| CN113432723B (en) | Image processing method, system and computer system for weakening stray radiation | |
| CN108225571A (en) | A kind of adaptive non-uniformity correction algorithm in short-wave infrared focal plane | |
| CN103997611B (en) | Method for suppressing image fixed-pattern noise based on noise template | |
| CN112393804B (en) | Image correction method and device | |
| CN103868601A (en) | Bilateral total variation regularization correction method for non-uniform response of IRFPA detector | |
| CN111998961B (en) | Infrared focal plane non-uniformity correction method based on pixel point temperature drift estimation | |
| Geng et al. | An adjacent differential statistics method for IRFPA nonuniformity correction | |
| CN110850500B (en) | Infrared image multi-section single-point correction parameter correction method |
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 |