CN105136308A - Adaptive correction method under variable integral time of infrared focal plane array - Google Patents

Adaptive correction method under variable integral time of infrared focal plane array Download PDF

Info

Publication number
CN105136308A
CN105136308A CN201510272545.2A CN201510272545A CN105136308A CN 105136308 A CN105136308 A CN 105136308A CN 201510272545 A CN201510272545 A CN 201510272545A CN 105136308 A CN105136308 A CN 105136308A
Authority
CN
China
Prior art keywords
coefficient
integral time
focal plane
plane array
gain
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
Application number
CN201510272545.2A
Other languages
Chinese (zh)
Other versions
CN105136308B (en
Inventor
张守荣
王�华
张大鹏
张旭
刘涛
潘卫军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Space Research Mechanical and Electricity
Original Assignee
Beijing Institute of Space Research Mechanical and Electricity
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Space Research Mechanical and Electricity filed Critical Beijing Institute of Space Research Mechanical and Electricity
Priority to CN201510272545.2A priority Critical patent/CN105136308B/en
Publication of CN105136308A publication Critical patent/CN105136308A/en
Application granted granted Critical
Publication of CN105136308B publication Critical patent/CN105136308B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Transforming Light Signals Into Electric Signals (AREA)
  • Radiation Pyrometers (AREA)
  • Photometry And Measurement Of Optical Pulse Characteristics (AREA)

Abstract

Provided is an adaptive correction method under variable integral time of an infrared focal plane array, comprising the following steps: first, setting different integral time for an infrared focal plane array, and calculating the average image data output collected by detection units through a two-point temperature calibration approach; working out the corresponding gain calibration coefficient and the corresponding offset calibration coefficient of each effective pixel point under different integral time conditions; then, drawing a line chart indicating the change of the gain calibration coefficient and the offset calibration coefficient of each effective pixel point with the integral time; obtaining the gain fitting coefficient and the offset fitting coefficient of each integral time interval through least square fitting; and finally, calculating the offset correction coefficient and the gain correction coefficient, and correcting an image acquired by the infrared focal plane array. As the correction coefficient and the integral time are in one-to-one correspondence relationship, the influence of change in integral time to the uniformity of the focal plane array can be isolated. Therefore, the uniformity of infrared images is improved significantly, and the uniformity of images is not affected by integral time modification.

Description

A kind of infrared focal plane array becomes the self-adapting correction method under integral time
Technical field
The present invention relates to a kind of infrared focal plane array self-adapting correction method, particularly a kind of infrared focal plane array becomes the self-adapting correction method under integral time, belongs to infrared imagery technique field.
Background technology
Refer to that infrared focal plane array probe unit produces the time of stored charge when being subject to extraneous radiation integral time.For infrared focal plane array image-forming system, the physical parameters such as output voltage, responsiveness, noise and specific detecivity are all in close relations with integral time, and suitable integral time can the overall performance of elevator system, otherwise, then suppress system performance.
In the practical application of staring infrared focal plane imaging system, often need to revise the detector integrates time according to Target Infrared Radiation intensity.And after carrying out Nonuniformity Correction to image, even if respond well, if change the integral time of detector, due to the discrete feature of its sensing circuit, device, amplifying circuit and the difference of process aspect, the homogeneity of image still can obviously be deteriorated.This essence of phenomena is because each probe unit output circuit causes the nonuniformity that integral time responds.
The application of integral time is regulated to be divided into two kinds, one is integral time is fixed gear place value, another kind is integral time is dynamic on-fixed gear value, for the first application can gear each integral time adopt existing engineering uses more and than be easier to realize method---two point correction method calculates respective correction coefficient, but this method is for rear one application, obviously no longer applicable.For the second application, existing method is according to demarcation interval section integral time, and the correction coefficient that the correction coefficient on putting the integral time in every section of interval all uses put interval certain integral time interior replaces.But existing method just weakens amendment integral time to the heteropical impact of pixel by Stepwise calibration integral time, but still exist the heterogeneity that integral time responds, so calibration result is not very desirable due to each probe unit output circuit.
Summary of the invention
The technical matters that the present invention solves is: overcome the deficiencies in the prior art, and present invention achieves a kind of in difference section integral time, correction coefficient is along with the asymmetric correction method of linear change integral time.Because correction coefficient and integral time are one_to_one corresponding, isolate the impact of the change focal plane Array Uniformity of integral time.Not only make the homogeneity of infrared image be improved by a relatively large margin, and amendment image conformity integral time is unaffected.
The technology used in the present invention solution: a kind of infrared focal plane array becomes the self-adapting correction method under integral time, and step is as follows:
(1) infrared focal plane array is aimed at black matrix by optical system, make blackbody radiation be full of the whole visual field of infrared focal plane array, and uniform irradiation is on infrared focal plane array;
(2) integral time of infrared focal plane array is set to INT 1;
(3) blackbody temperature is set to T 1, record the view data that each unit of infrared focal plane array seeker collects, gather the view data of Q frame, calculate T 1the averaged image data that at temperature, (i, j) individual probe unit collects exports and then try to achieve the averaged image data that all probe units collect and export
(4) blackbody temperature is set to T 2, utilize method in step (3), calculate T 2the averaged image data that at temperature, (i, j) individual probe unit collects exports and then try to achieve the averaged image data that all probe units collect and export
(5) choose P effectively pixel point, i.e. P effective probe unit, utilize the data in step (3) and step (4), calculating each effective pixel point is INT in integral time 1gain calibration coefficient under condition and biased calibration coefficient;
(6) integral time of the middle infrared focal plane array of modify steps (2), repeat step (3) ~ step (5), try to achieve each effective pixel point integral time after the modification corresponding gain calibration coefficient and biased calibration coefficient under condition;
(7) utilize the gain calibration coefficient of each effective pixel point of the Plotting data in step (6) and biased calibration coefficient respectively with the distribution plan of integral time, by the integral time of order from small to large, the adjacent gain calibration coefficient corresponding to each effective pixel point is connected with biased calibration coefficient straight line in distribution plan, obtains the gain calibration coefficient of each effective pixel point and biased calibration coefficient respectively with the broken line graph of change integral time;
(8) according to the turning point number n in step (7) middle polyline figure, n-1 section will be divided into integral time interval, in every section of interval gain calibration coefficient and biased calibration coefficient and integral time linear;
(9) in every section of integration time interval that step (8) is determined, choose n1 integral time point, use least square fitting to obtain gain fitting coefficient R ', the S ' of this section of integration time interval and biased fitting coefficient R ", S ";
(10) the gain fitting coefficient R ', the S ' that utilize step (9) to obtain and biased fitting coefficient R ", S ", calculate bias correction coefficient and gain correction coefficient;
(11) correct image that the bias correction coefficient that obtains in step (10) and gain correction coefficient obtain infrared focal plane array is utilized.
Blackbody temperature is set to T in (3) by described step 1, record the view data that each unit of infrared focal plane array seeker collects, gather the view data of Q frame, calculate T 1the averaged image data that at temperature, (i, j) probe unit collects exports specifically by formula:
V ij T 1 = Σ V ijm T 1 / Q
Provide, wherein for T 1the m frame image data that at temperature, (i, j) individual probe unit collects;
And then try to achieve the averaged image data that all probe units collect and export specifically by formula:
V ‾ ij T 1 = 1 M × N Σ V ij T 1
Provide, wherein, M, N are respectively line number and the columns of infrared focal plane array.
Use gain fitting coefficient R ', the S ' of least square fitting this section of integration time interval and biased fitting coefficient R ", S " in described step (9); Be specially:
Use gain fitting coefficient R ', the S ' of least square fitting this section of integration time interval specifically by formula:
R′*Int 1+S′＝K 1
R′*Int 2+S′＝K 2
......
R′*Int n1+S′＝K n1
Provide, wherein Int 1, Int 2..., Int n1and K 1, K 2... K n1n1 the integral time being respectively selected puts corresponding integral time and gain calibration coefficient;
Use the biased fitting coefficient R ", S " of least square fitting this section of integration time interval specifically by formula:
R″*Int 1+S″＝B 1
R″*Int 2+S″＝B 2
......
R″*Int n1+S″＝B n1
Provide, wherein Int 1, Int 2..., Int n1and B 1, B 2..., B n1n1 the integral time being respectively selected puts corresponding integral time and biased calibration coefficient.
The gain fitting coefficient R ', the S ' that utilize step (9) to obtain in described step (10) and biased fitting coefficient R ", S ", calculate bias correction coefficient and gain correction coefficient; Specifically by formula: use gain correction coefficient and bias correction COEFFICIENT K, B, wherein, K, B are specifically by formula:
K＝R′*Int+S′
B＝R″*Int+S″
Provide, wherein, K and B is respectively gain correction coefficient and bias correction coefficient, and Int is the current integration time, and R ' and S ' is gain fitting coefficient, and R " and S " is biased fitting coefficient.
Described in described step (3), the span of Q is: Q>=30.
In described step (5), the span of P is: P>=5.
In described step (9), n1 is more than or equal to 3.
The present invention's beneficial effect is compared with the conventional method:
(1) the present invention is by calculating the two point correction Coefficient Fitting under different integral time, find out the linear relationship of correction coefficient and integral time, compare the correction coefficient of point each integral time in existing method interval all by the method that the correction coefficient put certain integral time replaces, error is less, correction coefficient and integral time are one_to_one corresponding, match completely, isolate the impact of the change focal plane Array Uniformity of integral time like this, infrared focal plane array homogeneity obtains larger raising.
(2) when the present invention revises integral time, correction coefficient can recalculate, and the correction coefficient calculated matched with the current integration time, so image conformity is unaffected.And existing method revises integral time in interval, correction coefficient remains constant, and obvious image conformity can be deteriorated.
(3) the present invention is by the variation relation drawing different pixel point fitting coefficient and put integral time, can obtain linear segmented integral time very intuitively.Compare existing manual segmentation method, more intuitively convenient, and be not easy to make mistakes, correction accuracy is high.
(4) what the present invention mainly comprised is multiplication and additive operation, and calculated amount is little, can realize real time correction.
Accompanying drawing explanation
Fig. 1 chooses effective pixel dot gains calibration coefficient and the graph of a relation of integral time;
Fig. 2 chooses the graph of a relation that effective pixel point is biased calibration coefficient and integral time;
Fig. 3 is the inventive method realization flow figure;
Fig. 4 is the image before infrared camera high-gain, 4ms correction;
Fig. 5 is infrared camera high-gain, 4ms uses the image after the correction of existing method;
Fig. 6 is infrared camera high-gain, 4ms uses the image after method of the present invention correction.
Embodiment
Below in conjunction with accompanying drawing, the specific embodiment of the invention is described further.
One, Method And Principle
Adopt two point correction method, calculate calibration coefficient on some T integral time (integral time point choose will cover scope whole integral time), choose calibration coefficient and the graph of a relation of integral time that several effective pixel points draw the pixel point chosen, as depicted in figs. 1 and 2, wherein Fig. 1 chooses effective pixel dot gains calibration coefficient and the graph of a relation of integral time; Fig. 2 chooses the graph of a relation that effective pixel point is biased calibration coefficient and integral time; As can be seen from Fig. 1 and Fig. 2:
In different integration time interval, the biased and gain calibration coefficient of each pixel point and the relation linearly changed integral time.Again choose pixel point, curve plotting figure, the result obtained also is the same.Therefore, according to correction coefficient with piecewise linearity change integral time principle, least square fitting curve can be used in each correction coefficient linearity range, finds out every section of integration time interval, the linear relationship of calibration coefficient and integral time.When revising integral time, according to the interval at place integral time, the fitting coefficient of present segment just can be obtained by look-up table, the correction coefficient under the current integration time can be calculated according to formula (2) (3), then complete correction calculation according to formula (1).
Adopt correction coefficient K, B all to calculate along with the change of integral time, calibration model is
Y＝K*X+B(1)
K＝R′*Int+S′(2)
B＝R″*Int+S″(3)
The gain fitting coefficient that wherein R ', the S ' calculated gains that is present segment integration time interval is corresponding, R ", S " is biased corresponding biased fitting coefficient for the calculating of current integration time interval.
Two point temperature correction method be conduct a research the earliest, one of method the most ripe, generally use the infrared radiation of black matrix as input during two point correction, this input is uniform, but the heterogeneity due to infrared eye causes the uneven of output, the average generally output calibration uneven for each pixel exported to whole focal plane is gone, like this, according to the real output value of average and each point, Nonuniformity Correction coefficient k can be calculated i,j, b i,j.When inputting uneven, the output after can going out correct with nonuniformity correction formulae discovery.Peg method can represent with formula (4)
V′ i，j(Φ)＝k i,jV i,j(Φ)+b i,j(4)
Wherein: V i,j(Φ) for infrared focal plane array is when blackbody radiation degree is Φ, the output that (i, j) probe unit collects, k i,j, b i,jfor two point correction coefficient, V ' i, j(Φ) for each probe unit exports the output after two point correction.
According to detector responding range, choose Φ l, Φ htwo different irradiance carry out calibrating the response V obtaining detector cells i,j(Φ l), V i,j(Φ h). Φ respectively l, Φ hthe response mean value of each detector cells under homogeneous radiation.
Namely
V i , j ‾ ( Φ L ) = 1 M × N Σ V i , j ( Φ L ) - - - ( 5 )
V i , j ‾ ( Φ H ) = 1 M × N Σ V i , j ( Φ H ) - - - ( 6 )
By
V i , j ‾ ( Φ L ) = k i , j V i , j ( Φ L ) + b i , j - - - ( 7 )
V i , j ‾ ( Φ H ) = k i , j V i , j ( Φ H ) + b i , j - - - ( 8 )
Wherein, M, N are respectively line number and the columns of infrared focal plane array, and gain calibration coefficient and biased calibration coefficient can be calculated by following formula:
k i , j = V i , j ‾ ( Φ H ) - V i , j ‾ ( Φ L ) V i , j ( Φ H ) - V i , j ( Φ L ) b i , j = V i , j ( Φ H ) V i , j ‾ ( Φ L ) - V i , j ( Φ L ) V i , j ‾ ( Φ H ) V i , j ( Φ H ) - V i , j ( Φ L ) - - - ( 9 )
Two, detailed design
Be illustrated in figure 3 the inventive method specific implementation process flow diagram, mainly comprise two parts: design factor and correction calculation, the Detailed operating procedures of design factor is as follows:
(1) infrared focal plane array is aimed at black matrix by optical system, make blackbody radiation be full of the whole visual field of infrared focal plane array, and uniform irradiation is on infrared focal plane array;
(2) integral time of focal plane is set to INT 1;
(3) blackbody temperature is set to T 1, record the view data that each unit of infrared focal plane array seeker collects, gather the view data of Q frame, calculate T 1the averaged image data that at temperature, (i, j) individual probe unit collects exports and then try to achieve the averaged image data that all probe units collect and export specifically by formula:
V ij T 1 = Σ V ijm T 1 / Q
Provide, wherein for T 1the m frame image data that at temperature, (i, j) individual probe unit collects;
And then try to achieve the averaged image data that all probe units collect and export specifically by formula:
V ‾ ij T 1 = 1 M × N Σ V ij T 1
Provide, wherein, M, N are respectively line number and the columns of infrared focal plane array, and the span of described Q is: Q>=30.
(4) blackbody temperature is set to T 2, utilize method in step (3), calculate T 2the averaged image data that at temperature, (i, j) individual probe unit collects exports and then try to achieve the averaged image data that all probe units collect and export
(5) choose P effectively pixel point, i.e. P effective probe unit, utilize the data in step (3) and step (4), calculating each effective pixel point is INT in integral time 1gain calibration coefficient under condition and biased calibration coefficient;
(6) integral time of the middle infrared focal plane array of modify steps (2), repeat step (3) ~ step (5), try to achieve each effective pixel point integral time after the modification corresponding gain calibration coefficient and biased calibration coefficient under condition; The span of described P is: P>=5;
(7) utilize the gain calibration coefficient of each effective pixel point of the Plotting data in step (6) and biased calibration coefficient respectively with the distribution plan of integral time, by the integral time of order from small to large, the adjacent gain calibration coefficient corresponding to each effective pixel point is connected with biased calibration coefficient straight line in distribution plan, obtains the gain calibration coefficient of each effective pixel point and biased calibration coefficient respectively with the broken line graph of change integral time;
(8) according to the turning point number n in step (7) middle polyline figure, n-1 section will be divided into integral time interval, in every section of interval gain calibration coefficient and biased calibration coefficient and integral time linear;
(9) in every section of integration time interval that step (8) is determined, choose n1 integral time point, use least square fitting to obtain gain fitting coefficient R ', the S ' of this section of integration time interval and biased fitting coefficient R ", S "; Described n1 is more than or equal to 3;
Use gain fitting coefficient R ', the S ' of least square fitting this section of integration time interval specifically by formula according to formula (10):
R′*Int 1+S′＝K 1
R′*Int 2+S′＝K 2(10)
......
R′*Int n1+S′＝K n1
Provide, wherein Int 1, Int 2..., Int n1and K 1, K 2... K n1n1 the integral time being respectively selected puts corresponding integral time and gain calibration coefficient;
Use the biased fitting coefficient R ", S " of least square fitting this section of integration time interval specifically by formula:
R″*Int 1+S″＝B 1
R″*Int 2+S″＝B 2
......
R″*Int n1+S″＝B n1
Provide, wherein Int 1, Int 2..., Int n1and B 1, B 2..., B n1n1 the integral time being respectively selected puts corresponding integral time and biased calibration coefficient.
(10) the gain fitting coefficient R ', the S ' that utilize step (9) to obtain and biased fitting coefficient R ", S ", calculate bias correction coefficient and gain correction coefficient;
(11) correct image that the bias correction coefficient that obtains in step (10) and gain correction coefficient obtain infrared focal plane array is utilized.
When revising integral time, re-addressing the biased and gain fitting coefficient under the current integration time, calculating the correction coefficient under the renewal current integration time, complete new correction calculation.
Use FPGA to realize, mainly contain two kinds of computings: addition and multiplication.For the ease of hardware implementing, when meeting precision, correction coefficient is normalized to fixed-point number, to adopt fixed point totalizer and fixed-point multiplication device.Consider the fixed-point integer computing can only carrying out finite length in FPGA, therefore K, B must be integers, and simultaneously for ensureing computational accuracy, so just need first to COEFFICIENT K, B amplifies, and is reduced after having calculated again, and shift right operation can be utilized to replace division arithmetic.When completing Y=KX+B and calculating, it should be noted that B is the data with symbol, sign bit can be set in the most significant digit of B, complete the computing of addition or subtraction.Finally also need to carry out anti-spilled process to the result calculated, (Y>2 when namely there is overflow m-1, m is the valid data position of view data), result is set to 2 m-1, when overflowing under appearance, result is set to 0.
Embodiment
At different temperature, revise integral time, before measuring and calculation correction, existing bearing calibration correction is rear and use the rear image non-uniform residual volume of the inventive method correction.As shown in Table 1 to Table 3, wherein, table 1 is for correcting front original image heterogeneity residual volume statistical form, and table 2 is heterogeneity residual volume statistical form after existing method corrects, and table 3 is rear heterogeneity residual volume statistical form for the present invention corrects for test result;
Table 1
T/Int 0.1 1 2 4 6 10 15
20 0.0646 0.0920 0.1403 0.0946 0.0984 0.0978 0.0721
22 0.0653 0.0886 0.1318 0.0904 0.0938 0.0941 0.0685
27 0.0667 0.0829 0.1149 0.0847 0.0883 0.0885 0.0637
32 0.0680 0.0793 0.1019 0.0818 0.0850 0.0858 0.0391
35 0.0695 0.0770 0.0958 0.08083 0.0839 0.0846 0.0073
37 0.0698 0.0751 0.0923 0.0804 0.0832 0.0843 0.0054
42 0.0703 0.0740 0.08651 0.0802 0.0833 0.0656 0.0050
47 0.0702 0.0732 0.0835 0.0804 0.0836 0.0139 0.0048
52 0.0707 0.0740 0.0824 0.0812 0.0851 0.0038 0.0042
57 0.0695 0.0749 0.0824 0.0822 0.0579 0.0038 0.0041
62 0.0687 0.0760 0.0834 0.0827 0.0119 0.0036 0.0042
67 0.0676 0.0782 0.0853 0.0773 0.0035 0.0036 0.0042
Table 2
T/Int 0.1 1 2 4 6 10 15
20 0.0806 0.0412 0.0482 0.0498 0.0186 0.0417 0.0171
22 0.0790 0.0345 0.0211 0.0427 0.0022 0.0366 0.0154
27 0.0740 0.0300 0.0234 0.0352 0.0031 0.0303 0.0265
32 0.0672 0.0283 0.0322 0.0295 0.0034 0.0250 0.0151
35 0.0698 0.0248 0.0265 0.0264 0.0035 0.0227 0.0070
37 0.0621 0.0213 0.0247 0.0245 0.0040 0.0206 0.0050
42 0.0563 0.0189 0.0147 0.0206 0.0032 0.0348 0.0043
47 0.0538 0.0165 0.0126 0.0179 0.0028 0.0105 0.0040
52 0.0535 0.0138 0.0115 0.0147 0.0009 0.0030 0.0040
57 0.0454 0.0136 0.0125 0.0133 0.0411 0.0029 0.0041
62 0.0448 0.0122 0.0110 0.0113 0.0445 0.0029 0.0040
67 0.0394 0.0121 0.0110 0.0211 0.0445 0.0029 0.0040
Table 3
T/Int 0.1 1 2 4 6 10 15
20 0.0202 0.0187 0.0225 0.0210 0.0175 0.0121 0.0028
22 0.0200 0.0115 0.0138 0.0100 0.0021 0.0066 0.0054
27 0.0184 0.0097 0.0121 0.0099 0.0030 0.0053 0.0029
32 0.0182 0.0096 0.0110 0.0112 0.0034 0.0049 0.0030
35 0.0190 0.0107 0.0106 0.0122 0.0035 0.0023 0.0030
37 0.0195 0.0098 0.0106 0.0131 0.0040 0.0086 0.0031
42 0.0205 0.0084 0.0097 0.0098 0.0032 0.0009 0.0031
47 0.0213 0.0066 0.0096 0.0088 0.0027 0.0016 0.0031
52 0.0224 0.0066 0.0093 0.0066 0.0001 0.0016 0.0030
57 0.0218 0.0064 0.0094 0.0075 0.0005 0.0014 0.0030
62 0.0216 0.0063 0.0094 0.0070 0.0012 0.0013 0.0030
67 0.0207 0.0036 0.0088 0.0079 0.0014 0.0013 0.0030
After can finding out use bearing calibration of the present invention from table 1 and table 3 contrast, image is compared the heterogeneity remnants before correction and is reduced 1 order of magnitude, and calibration result is obvious, and after correcting, the non-homogeneous remnants of image are very little.
Contrast as can be seen from table 2 and table 3, use bearing calibration of the present invention to compare existing method, image non-uniform remnants reduce a lot, substantially in the value of ppt, existing method is compared in bearing calibration of the present invention, and image conformity obtains raising by a relatively large margin.
Accompanying drawing 4 is the image before infrared camera high-gain 4ms corrects, Fig. 5 is the image after using existing method to correct, Fig. 6 is the image after using method of the present invention to correct, by the contrast of three width pictures, can find out, before correcting, image conformity is very poor, there is heavier band interference, use the band in method of the present invention and existing method equal energy removal of images, but after using existing method, parts of images is partially bright on the lower, there is the problem that upper and lower two parts gray-scale value difference is larger, and using method of the present invention to there is not this phenomenon, integral image homogeneity is better.Comprehensive above analysis, after known use bearing calibration of the present invention, non-homogeneous remnants are minimum, reduce 1 order of magnitude before comparatively correcting, and after correcting, image display effect is optimum.
The unspecified part of the present invention belongs to general knowledge as well known to those skilled in the art.

Claims (7)

1. infrared focal plane array becomes the self-adapting correction method under integral time, it is characterized in that step is as follows:
(1) infrared focal plane array is aimed at black matrix by optical system, make blackbody radiation be full of the whole visual field of infrared focal plane array, and uniform irradiation is on infrared focal plane array;
(2) integral time of infrared focal plane array is set to INT 1;
(3) blackbody temperature is set to T 1, record the view data that each unit of infrared focal plane array seeker collects, gather the view data of Q frame, calculate T 1the averaged image data that at temperature, (i, j) individual probe unit collects exports and then try to achieve the averaged image data that all probe units collect and export
(4) blackbody temperature is set to T 2, utilize method in step (3), calculate T 2the averaged image data that at temperature, (i, j) individual probe unit collects exports and then try to achieve the averaged image data that all probe units collect and export
(5) choose P effectively pixel point, i.e. P effective probe unit, utilize the data in step (3) and step (4), calculating each effective pixel point is INT in integral time 1gain calibration coefficient under condition and biased calibration coefficient;
(6) integral time of the middle infrared focal plane array of modify steps (2), repeat step (3) ~ step (5), try to achieve each effective pixel point integral time after the modification corresponding gain calibration coefficient and biased calibration coefficient under condition;
(7) utilize the gain calibration coefficient of each effective pixel point of the Plotting data in step (6) and biased calibration coefficient respectively with the distribution plan of integral time, by the integral time of order from small to large, the adjacent gain calibration coefficient corresponding to each effective pixel point is connected with biased calibration coefficient straight line in distribution plan, obtains the gain calibration coefficient of each effective pixel point and biased calibration coefficient respectively with the broken line graph of change integral time;
(8) according to the turning point number n in step (7) middle polyline figure, n-1 section will be divided into integral time interval, in every section of interval gain calibration coefficient and biased calibration coefficient and integral time linear;
(9) in every section of integration time interval that step (8) is determined, choose n1 integral time point, use least square fitting to obtain gain fitting coefficient R ', the S ' of this section of integration time interval and biased fitting coefficient R ", S ";
(10) the gain fitting coefficient R ', the S ' that utilize step (9) to obtain and biased fitting coefficient R ", S ", calculate bias correction coefficient and gain correction coefficient;
(11) correct image that the bias correction coefficient that obtains in step (10) and gain correction coefficient obtain infrared focal plane array is utilized.
2. a kind of infrared focal plane array according to claim 1 becomes the self-adapting correction method under integral time, it is characterized in that: blackbody temperature is set to T in (3) by described step 1, record the view data that each unit of infrared focal plane array seeker collects, gather the view data of Q frame, calculate T 1the averaged image data that at temperature, (i, j) probe unit collects exports specifically by formula:
V ij T 1 = Σ V ijm T 1 / Q
Provide, wherein for T 1the m frame image data that at temperature, (i, j) individual probe unit collects;
And then try to achieve the averaged image data that all probe units collect and export specifically by formula:
V ‾ ij T 1 = 1 M × N Σ V ij T 1
Provide, wherein, M, N are respectively line number and the columns of infrared focal plane array.
3. a kind of infrared focal plane array according to claim 1 becomes the self-adapting correction method under integral time, it is characterized in that: use gain fitting coefficient R ', the S ' of least square fitting this section of integration time interval and biased fitting coefficient R ", S " in described step (9); Be specially:
Use gain fitting coefficient R ', the S ' of least square fitting this section of integration time interval specifically by formula:
R′*Int 1+S′＝K 1
R′*Int 2+S′＝K 2
......
R′*Int n1+S′＝K n1
Provide, wherein Int 1, Int 2..., Int n1and K 1, K 2... K n1n1 the integral time being respectively selected puts corresponding integral time and gain calibration coefficient;
Use the biased fitting coefficient R ", S " of least square fitting this section of integration time interval specifically by formula:
R″*Int 1+S″＝B 1
R″*Int 2+S″＝B 2
......
R″*Int n1+S″＝B n1
Provide, wherein Int 1, Int 2..., Int n1and B 1, B 2..., B n1n1 the integral time being respectively selected puts corresponding integral time and biased calibration coefficient.
4. a kind of infrared focal plane array according to claim 1 becomes the self-adapting correction method under integral time, it is characterized in that: the gain fitting coefficient R ', the S ' that utilize step (9) to obtain in described step (10) and biased fitting coefficient R ", S ", calculate bias correction coefficient and gain correction coefficient; Specifically by formula: use gain correction coefficient and bias correction COEFFICIENT K, B, wherein, K, B are specifically by formula:
K＝R′*Int+S′
B＝R″*Int+S″
Provide, wherein, K and B is respectively gain correction coefficient and bias correction coefficient, and Int is the current integration time, and R ' and S ' is gain fitting coefficient, and R " and S " is biased fitting coefficient.
5. a kind of infrared focal plane array according to claim 1 becomes the self-adapting correction method under integral time, it is characterized in that: described in described step (3), the span of Q is: Q>=30.
6. a kind of infrared focal plane array according to claim 1 becomes the self-adapting correction method under integral time, it is characterized in that: in described step (5), the span of P is: P>=5.
7. a kind of infrared focal plane array according to claim 1 becomes the self-adapting correction method under integral time, it is characterized in that: in described step (9), n1 is more than or equal to 3.
CN201510272545.2A 2015-05-25 2015-05-25 A kind of infrared focal plane array becomes the self-adapting correction method under the time of integration Active CN105136308B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510272545.2A CN105136308B (en) 2015-05-25 2015-05-25 A kind of infrared focal plane array becomes the self-adapting correction method under the time of integration

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510272545.2A CN105136308B (en) 2015-05-25 2015-05-25 A kind of infrared focal plane array becomes the self-adapting correction method under the time of integration

Publications (2)

Publication Number Publication Date
CN105136308A true CN105136308A (en) 2015-12-09
CN105136308B CN105136308B (en) 2018-05-01

Family

ID=54721741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510272545.2A Active CN105136308B (en) 2015-05-25 2015-05-25 A kind of infrared focal plane array becomes the self-adapting correction method under the time of integration

Country Status (1)

Country Link
CN (1) CN105136308B (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106644078A (en) * 2016-12-30 2017-05-10 电子科技大学 Terahertz image non-uniformity correction method
CN106768385A (en) * 2016-12-13 2017-05-31 航天海鹰光电信息技术（天津）有限公司 A kind of detector heterogeneity degeneration restorative procedure
CN107421643A (en) * 2017-07-19 2017-12-01 沈阳上博智像科技有限公司 The bearing calibration of infrared image and device
CN108120507A (en) * 2016-11-30 2018-06-05 北京航天计量测试技术研究所 A kind of thermal infrared imager high dynamic range adaptive regulation method
CN108921796A (en) * 2018-06-07 2018-11-30 西安电子科技大学 A kind of Infrared Image Non-uniformity Correction method based on deep learning
CN110686781A (en) * 2019-04-12 2020-01-14 福建鼎泰康医疗设备有限公司 Temperature calibration method and device
CN111780876A (en) * 2020-06-30 2020-10-16 平安国际智慧城市科技股份有限公司 Temperature measurement method, temperature measurement device, electronic equipment and computer readable storage medium
CN113473114A (en) * 2021-07-09 2021-10-01 中国人民解放军91550部队 Method for calibrating pixel non-uniformity of polarization camera with focal plane
CN114397024A (en) * 2021-11-17 2022-04-26 中国航空工业集团公司洛阳电光设备研究所 Temperature segment dividing and switching method based on variable integration time multi-segment correction algorithm
CN117571138A (en) * 2024-01-15 2024-02-20 中国科学院长春光学精密机械与物理研究所 Inter-chip non-uniformity correction method based on infrared detector

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013055274A2 (en) * 2011-10-14 2013-04-18 Flir Systems Ab Method for adjusting the integration time for an ir detector, and an ir camera for implementing the method
CN103604503A (en) * 2013-06-17 2014-02-26 中国航天科工集团第三研究院第八三五八研究所 High-precision non-uniformity correction method for dynamic adjustment of integration time
CN104251742A (en) * 2014-09-30 2014-12-31 成都市晶林科技有限公司 Two-point linear and target and environment-based binary non-linear infrared detector heterogeneity correction method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2013055274A2 (en) * 2011-10-14 2013-04-18 Flir Systems Ab Method for adjusting the integration time for an ir detector, and an ir camera for implementing the method
CN103604503A (en) * 2013-06-17 2014-02-26 中国航天科工集团第三研究院第八三五八研究所 High-precision non-uniformity correction method for dynamic adjustment of integration time
CN104251742A (en) * 2014-09-30 2014-12-31 成都市晶林科技有限公司 Two-point linear and target and environment-based binary non-linear infrared detector heterogeneity correction method

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108120507A (en) * 2016-11-30 2018-06-05 北京航天计量测试技术研究所 A kind of thermal infrared imager high dynamic range adaptive regulation method
CN108120507B (en) * 2016-11-30 2020-06-05 北京航天计量测试技术研究所 Thermal infrared imager high dynamic range self-adaptive adjustment method
CN106768385A (en) * 2016-12-13 2017-05-31 航天海鹰光电信息技术（天津）有限公司 A kind of detector heterogeneity degeneration restorative procedure
CN106768385B (en) * 2016-12-13 2019-04-16 天津津航技术物理研究所 A kind of detector heterogeneity degeneration restorative procedure
CN106644078A (en) * 2016-12-30 2017-05-10 电子科技大学 Terahertz image non-uniformity correction method
CN107421643A (en) * 2017-07-19 2017-12-01 沈阳上博智像科技有限公司 The bearing calibration of infrared image and device
CN108921796A (en) * 2018-06-07 2018-11-30 西安电子科技大学 A kind of Infrared Image Non-uniformity Correction method based on deep learning
CN108921796B (en) * 2018-06-07 2021-09-03 西安电子科技大学 Infrared image non-uniformity correction method based on deep learning
CN110686781B (en) * 2019-04-12 2024-04-09 福建鼎泰康医疗设备有限公司 Temperature calibration method and device
CN110686781A (en) * 2019-04-12 2020-01-14 福建鼎泰康医疗设备有限公司 Temperature calibration method and device
CN111780876A (en) * 2020-06-30 2020-10-16 平安国际智慧城市科技股份有限公司 Temperature measurement method, temperature measurement device, electronic equipment and computer readable storage medium
CN111780876B (en) * 2020-06-30 2023-12-01 平安国际智慧城市科技股份有限公司 Temperature measurement method, temperature measurement device, electronic equipment and computer readable storage medium
CN113473114B (en) * 2021-07-09 2023-05-05 中国人民解放军91550部队 Method for calibrating pixel non-uniformity of focal plane splitting polarization camera
CN113473114A (en) * 2021-07-09 2021-10-01 中国人民解放军91550部队 Method for calibrating pixel non-uniformity of polarization camera with focal plane
CN114397024A (en) * 2021-11-17 2022-04-26 中国航空工业集团公司洛阳电光设备研究所 Temperature segment dividing and switching method based on variable integration time multi-segment correction algorithm
CN114397024B (en) * 2021-11-17 2023-09-29 中国航空工业集团公司洛阳电光设备研究所 Temperature segment dividing and switching method based on variable integral time multi-segment correction algorithm
CN117571138A (en) * 2024-01-15 2024-02-20 中国科学院长春光学精密机械与物理研究所 Inter-chip non-uniformity correction method based on infrared detector
CN117571138B (en) * 2024-01-15 2024-03-19 中国科学院长春光学精密机械与物理研究所 Inter-chip non-uniformity correction method based on infrared detector

Also Published As

Publication number Publication date
CN105136308B (en) 2018-05-01

Similar Documents

Publication Publication Date Title
CN105136308A (en) Adaptive correction method under variable integral time of infrared focal plane array
CN101776486B (en) Method for correcting non-uniformity fingerprint pattern on basis of infrared focal plane
CN101776487B (en) Infrared focal plane non-uniformity fingerprint extraction and image correction method
KR102022812B1 (en) Image Contrast Reinforcement
CN103335716B (en) Area array infrared camera calibration and non-uniformity correction method based on variable integration time
CN103335724B (en) Calibration-based scene self-adaption IRFPA heterogeneity correction method
CN106782283A (en) The gamma curve correcting method and display device of a kind of display device
CN103604503A (en) High-precision non-uniformity correction method for dynamic adjustment of integration time
CN100585634C (en) A kind of relative radiation correction method of star-load TDICCD camera
CN104599632B (en) Improve the method and device of OLED display brightness uniformities
CN104269157B (en) A kind of OLED screen gamma correction method
CN108846805A (en) A kind of thermal-induced imagery two o'clock non-uniform correction method based on scene adaptive
CN106197690B (en) Image calibrating method and system under the conditions of a kind of wide temperature range
CN106855435A (en) Heterogeneity real-time correction method on long wave linear array infrared camera star
CN105719298A (en) Edge detection technology based line diffusion function extracting method
CN106017695A (en) Adaptive infrared non-uniformity calibration method based on motion state estimation
CN101674398A (en) Real-time CCD digital image data dark current correction method
CN102663714A (en) Saliency-based method for suppressing strong fixed-pattern noise in infrared image
CN105737990A (en) Detector temperature-based infrared image heterogeneity correction method and system
CN102855610A (en) Method for correcting infrared image heterogeneity by using parameter correctness factor
CN102589707B (en) Real-time compensation method of non-uniformity correction residual errors of infrared focal plane array detector
CN106644078A (en) Terahertz image non-uniformity correction method
CN103258334A (en) Method of estimating scene light source colors of color image
CN102410844B (en) Method and device for correcting non-uniformity of image of high-dynamic star sensor
CN115687850A (en) Method and device for calculating irrigation water demand of farmland

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant