CN106768349A - Adapt to the TDI scanning imaging system relative radiometric calibration methods of multichannel bidirectional output - Google Patents
Adapt to the TDI scanning imaging system relative radiometric calibration methods of multichannel bidirectional output Download PDFInfo
- Publication number
- CN106768349A CN106768349A CN201611035496.1A CN201611035496A CN106768349A CN 106768349 A CN106768349 A CN 106768349A CN 201611035496 A CN201611035496 A CN 201611035496A CN 106768349 A CN106768349 A CN 106768349A
- Authority
- CN
- China
- Prior art keywords
- section
- relative radiometric
- radiometric calibration
- imaging system
- multichannel
- 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 31
- 238000000034 method Methods 0.000 title claims abstract description 23
- 230000002457 bidirectional effect Effects 0.000 title claims abstract description 11
- 230000005855 radiation Effects 0.000 claims abstract description 22
- 238000012937 correction Methods 0.000 claims abstract description 8
- 238000012935 Averaging Methods 0.000 claims abstract description 5
- 241001269238 Data Species 0.000 claims description 8
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000003595 spectral effect Effects 0.000 claims description 3
- 238000003359 percent control normalization Methods 0.000 claims description 2
- 230000006978 adaptation Effects 0.000 abstract description 2
- 230000003287 optical effect Effects 0.000 description 4
- 230000000903 blocking effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000003331 infrared imaging Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000003333 near-infrared imaging Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000011896 sensitive detection 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
- 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)
- Facsimile Scanning Arrangements (AREA)
Abstract
The invention belongs to relative radiometric calibration technical field, and in particular to a kind of TDI scanning imaging system relative radiometric calibration methods of adaptation multichannel bidirectional output.The method including record greatest irradiation brightness, by interval record brightness data, calculate multichannel average radiation response data, linear segmented is carried out to collection point, 2 relative radiometric calibrations corrections are carried out to each section, record each section of head and the tail image averaging gray scale representative value, from steps such as corresponding relative radiometric calibration coefficients.Compared with prior art, the present invention solves the problems, such as multichannel TDI scanning imaging system relative radiometric calibrations, not only can control the relative radiometric calibration precision in single channel, and each interchannel relative radiometric calibration precision of recoverable.
Description
Technical field
The invention belongs to relative radiometric calibration technical field, and in particular to a kind of TDI of adaptation multichannel bidirectional output sweeps
Retouch imaging system relative radiometric calibration method.
Background technology
At present, optical remote sensing imaging system all increasingly requires possess the big breadth of acquisition, noise high at two aspects of the army and the people
Than target image ability, in order to obtain big breadth, high s/n ratio target image, TDI scanning imagery modes are generally used.TDI sweeps
Imaging system is retouched mainly comprising optical system and photoelectric detector, and one side optical system has vignetting, i.e. homogeneous target warp
Optical system conversion is crossed, energy is reallocated, cause focal plane center energy to be concentrated, gradually reduced to edge, cause Jiao
Face Energy distribution is uneven；Another aspect sensitive detection parts cause that each detection pixel is input into identical energy due to material intrinsic property
Response characteristic have differences, produce Non-uniform responsivity；Secondly TDI device data transfer rates it is high, device inside frequently with
Multichannel export structure form meets high-frequency data transfer, because different channel signals cause by different paths and different amplifiers
Image produces passage blocking artifact.The inconsistency that the comprehensive function of these three factors is produced can cause that the image for obtaining is produced very
Strong pattern noise, drastically influence image quality, be that of restriction multichannel TDI scanning imaging systems is important and crucial
Factor.Therefore, it is necessary to carry out relative radiometric calibration to multichannel TDI scanning imaging systems.
The relative radiometric calibration method of the TDI scanning imagers of the near-infrared multichannel bidirectional output of research at present, it is different
In the relative radiometric calibration method of infrared imaging system.Infrared imaging system relative radiometric calibration is not suitable for TDI type multichannels
The near infrared imaging system of output, is mainly manifested in each interchannel relative radiometric calibration poor effect, therefore, it is defeated for multichannel
The TDI scanning imaging systems for going out, interchannel relative radiometric calibration is influence picture quality most critical technology.
The content of the invention
(1) technical problem to be solved
The technical problem to be solved in the present invention is：How to overcome prior art not enough, there is provided it is defeated that one kind is applied to multichannel
The relative radiometric calibration method of the TDI scanning imaging systems for going out.
(2) technical scheme
In order to solve the above technical problems, the present invention provides a kind of TDI scanning imaging systems for adapting to multichannel bidirectional output
Relative radiometric calibration method, the method comprises the following steps：
Step 1：Imaging system alignment criteria radiation source, adjusts radiation source brightness, makes image saturation, and record greatest irradiation is bright
Degree Lmax；
Step 2：In 0-LmaxIn dynamic range, by interval delta L=Lmax/ N, N are interval quantity, gather N+1 width images, together
Shi Jilu brightness datas L1、L2……LN+1；
Step 3：Multichannel average radiation response data DN is calculated to N+1 width image1、DN2……DNN+1；Wherein, DN=L
× c, c are the spectral response of detector, obtain the brightness-average response data pair of detector；
Step 4：Linear segmented is carried out to N+1 collection point, i.e.,Its
In, m represents the section of the division of selection；When representing any input brightness L in the section, correspondence image it is flat
Equal gray scale；amRepresent that this section meets the slope of linear response；bmRepresent that this section radiates the intercept of linear response；
Step 5：According to the segmentation result in step 4,2 relative radiometric calibration corrections are carried out to each section, solved
Every section of relative radiometric calibration coefficient of all pixels of passage；
Step 6：Record each section of head and the tail image averaging gray scale representative value [Nimin, Nimax]；
Wherein, i=1,2,3 ..., which section represented；Min represents this section of image minimum gradation value；Max represents that this section is schemed
As maximum gradation value；
Step 7：During working flight, which section is imaging system belong to according to image intensity value, from corresponding relative spoke
Calibration coefficient is penetrated, realizes that multichannel scan image is corrected.
Wherein, the calibrated radiation source in the step 1 is solar simulator.
Wherein, it is segmented according to foundation in the step 4：
In the past 3 gathered datas start by least square fitting linearity curve, make inclined between matched curve and sampled point
Poor quadratic sum selects a for minimum constraintsm、bmCoefficient；
If preceding 3 gathered datas meet precision controlling constraints, it is fitted whether preceding 4 data meet precision controlling
Constraints, then pointwise is increased up i+1 data and meets precision controlling constraints, then take preceding i data as the 1st
Section；
I-th to the i-th+2 this 3 data is taken again according to the matched curve ... ... of precision controlling constraints until discontented
The 2nd section is determined untill sufficient precision controlling constraints；
Relative radiometric calibration correction section is divided in radiation full dynamic range according to same method.
Wherein, precision is segmented in the step 4 by 1% control.
(3) beneficial effect
Compared with prior art, the present invention solves the problems, such as multichannel TDI scanning imaging system relative radiometric calibrations, no
Only can control the relative radiometric calibration precision in single channel, and each interchannel relative radiometric calibration precision of recoverable.
Brief description of the drawings
Fig. 1 is multichannel TDI scanning imaging systems according to the sectional linear fitting that drafting is responded in radiation full dynamic range
Curve synoptic diagram.
Fig. 2 is the original image schematic diagram of multichannel TDI scanning imaging systems output.
After Fig. 3 is segmented precision controlling for multichannel TDI scanning imaging systems scan image according to 1%, then carry out two point phases
To the image schematic diagram after radiation calibration.
Specific embodiment
To make the purpose of the present invention, content and advantage clearer, with reference to the accompanying drawings and examples, to of the invention
Specific embodiment is described in further detail.
In order to solve the above technical problems, the present invention provides a kind of TDI scanning imaging systems for adapting to multichannel bidirectional output
Relative radiometric calibration method, the method comprises the following steps：
Step 1：Imaging system alignment criteria radiation source (solar simulator), adjusts radiation source brightness, makes image saturation, remembers
Record greatest irradiation brightness Lmax；
Step 2：In 0-LmaxIn dynamic range, by interval delta L=Lmax/ N, N are interval quantity, gather N+1 width images, together
Shi Jilu brightness datas L1、L2……LN+1；
Step 3：Multichannel average radiation response data DN is calculated to N+1 width image1、DN2……DNN+1；Wherein, DN=L
× c, c are the spectral response of detector, obtain the brightness-average response data pair of detector；
Step 4：Linear segmented is carried out to N+1 collection point, i.e.,Wherein, m represents drawing for selection
The section divided；When representing any input brightness L in the section, the average gray of correspondence image；amRepresent that the section meets line
Property response slope；bmRepresent that this section radiates the intercept of linear response；
Segmentation is according to foundation：In the past 3 gathered datas start by least square fitting linearity curve, make matched curve with
The quadratic sum of the deviation between sampled point selects a for minimum constraintsm、bmCoefficient；
If preceding 3 gathered datas meet precision controlling constraints, it is fitted whether preceding 4 data meet precision controlling
Constraints, then pointwise is increased up i+1 data and meets precision controlling constraints, then take preceding i data as the 1st
Section；
I-th to the i-th+2 this 3 data is taken again according to the matched curve ... ... of precision controlling constraints until discontented
The 2nd section is determined untill sufficient precision controlling constraints；
Relative radiometric calibration correction section is divided in radiation full dynamic range according to same method；
Accuracy constraint condition：At 4%, human eye does not see obvious difference in response to general relative radiation precision controlling substantially
It is different, but for multichannel TDI scanning imaging systems, passage blocking artifact calibration result is bad, experiment proves that, precision presses 1%
During control, passage blocking artifact calibration result is preferable；
Step 5：According to the segmentation result in step 4,2 relative radiometric calibration corrections are carried out to each section, solved
Every section of relative radiometric calibration coefficient of all pixels of passage；
Step 6：Record each section of head and the tail image averaging gray scale representative value [Nimin, Nimax]；
Wherein, i=1,2,3 ..., which section represented；Min represents this section of image minimum gradation value；Max represents that this section is schemed
As maximum gradation value；
Step 7：During working flight, which section is imaging system belong to according to image intensity value, from corresponding relative spoke
Calibration coefficient is penetrated, realizes that multichannel scan image is corrected.
With reference to specific embodiment, the present invention will be described in detail.
Embodiment
In the present embodiment：
(1) the detailed operation method of step 4
If finite number meets according to radiation linear function
f(Li)m=am·Lmi+bm (1)
In formula：
The section of m-division；
LmiI-th radiance in-m sections；
f(Li)mAny input brightness L in the-sectioniWhen, according to the image averaging gray scale that formula (1) is solved；
am- this section meets the slope of linear response；
bm- this section radiates the intercept of linear response；
Input radiation brightness L in-m sectionsiWhen, gather the average gray of image.
am, bmSolve as follows：
If each sample point data and fitting data deviation areSum of square of deviations isSeek sum of square of deviations partial derivative and obtained after arranging：
Relative deviation accordance judges as follows：
If deviation is UN, precision controlling constraints
(2) the specific calculating process of the coefficient of step 5
In m sections, the linear response of multichannel image average gray responder, using 2 relative radiometric calibration algorithms.
WithDuring for multichannel TDI scanning imaging systems in m sections, radiance is L1And L2When institute
There are the average of pixel, DNi(L1) and DNi(L2) for i-th pixel in L1、L2When response, by being calculated i-th pixel
KiAnd biValue：
Obtain the k of each pixeliAnd biAfterwards, whenWhen, then the relative radiation from m sections
Correction coefficient, each pixel is output as：
Using multichannel TDI scanning imaging systems relative radiometric calibration method of the invention, rdaiation response is carried out first
Segmentation is optimized, then different multichannel TDI scanning imaging system relative radiometric calibrations is realized under the conditions of segmentation is optimized.
The above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, on the premise of the technology of the present invention principle is not departed from, some improvement and deformation can also be made, these improve and deform
Also should be regarded as protection scope of the present invention.
Claims (4)
1. it is a kind of to adapt to the TDI scanning imaging system relative radiometric calibration methods that multichannel bidirectional is exported, it is characterised in that the party
Method comprises the following steps：
Step 1：Imaging system alignment criteria radiation source, adjusts radiation source brightness, makes image saturation, records greatest irradiation brightness
Lmax；
Step 2：In 0-LmaxIn dynamic range, by interval delta L=Lmax/ N, N are interval quantity, gather N+1 width images, while note
Record brightness data L1、L2……LN+1；
Step 3：Multichannel average radiation response data DN is calculated to N+1 width image1、DN2……DNN+1；Wherein, DN=L × c, c
It is the spectral response of detector, obtains the brightness-average response data pair of detector；
Step 4：Linear segmented is carried out to N+1 collection point, i.e.,Wherein, m represents the division of selection
Section；When representing any input brightness L in the section, the average gray of correspondence image；amRepresent that the section meets linear sound
The slope answered；bmRepresent that this section radiates the intercept of linear response；
Step 5：According to the segmentation result in step 4,2 relative radiometric calibration corrections are carried out to each section, solve every section
The relative radiometric calibration coefficient of all pixels of passage；
Step 6：Record each section of head and the tail image averaging gray scale representative value [Nimin, Nimax]；
Wherein, i=1,2,3 ..., which section represented；Min represents this section of image minimum gradation value；Max represents this section of image most
High-gray level value；
Step 7：During working flight, which section is imaging system belong to according to image intensity value, fixed from corresponding relative radiation
Mark coefficient, realizes that multichannel scan image is corrected.
2. it is as claimed in claim 1 to adapt to the TDI scanning imaging system relative radiometric calibration methods that multichannel bidirectional is exported, its
It is characterised by, the calibrated radiation source in the step 1 is solar simulator.
3. it is as claimed in claim 1 to adapt to the TDI scanning imaging system relative radiometric calibration methods that multichannel bidirectional is exported, its
It is characterised by, is segmented according to foundation in the step 4：
In the past 3 gathered datas start by least square fitting linearity curve, make deviation between matched curve and sampled point
Quadratic sum selects a for minimum constraintsm、bmCoefficient；
If preceding 3 gathered datas meet precision controlling constraints, it is fitted whether preceding 4 data meet precision controlling constraint
Condition, then pointwise is increased up i+1 data and meets precision controlling constraints, then take preceding i data as the 1st section；
I-th to the i-th+2 this 3 data is taken again according to precision controlling constraints matched curve ... ... until being unsatisfactory for essence
Determine the 2nd section untill degree control constraints condition；
Relative radiometric calibration correction section is divided in radiation full dynamic range according to same method.
4. it is as claimed in claim 1 to adapt to the TDI scanning imaging system relative radiometric calibration methods that multichannel bidirectional is exported, its
It is characterised by, precision is segmented in the step 4 by 1% control.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611035496.1A CN106768349B (en) | 2016-11-18 | 2016-11-18 | Adapt to the TDI scanning imaging system relative radiometric calibration method of multichannel bidirectional output |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611035496.1A CN106768349B (en) | 2016-11-18 | 2016-11-18 | Adapt to the TDI scanning imaging system relative radiometric calibration method of multichannel bidirectional output |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106768349A true CN106768349A (en) | 2017-05-31 |
| CN106768349B CN106768349B (en) | 2019-04-19 |
Family
ID=58971095
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611035496.1A Active CN106768349B (en) | 2016-11-18 | 2016-11-18 | Adapt to the TDI scanning imaging system relative radiometric calibration method of multichannel bidirectional output |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106768349B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107907222A (en) * | 2017-11-13 | 2018-04-13 | 国网山东省电力公司荣成市供电公司 | A kind of thermal infrared imaging electric power facility fault locator and detection method |
| CN108181005A (en) * | 2017-11-17 | 2018-06-19 | 天津津航技术物理研究所 | A kind of method and system for the debugging of TDI ccd detectors focal plane |
| US11604146B2 (en) | 2017-09-19 | 2023-03-14 | Beckman Coulter, Inc. | Analog light measuring and photon counting with a luminometer system for assay reactions in chemiluminescence measurements |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101226639A (en) * | 2008-01-29 | 2008-07-23 | 航天东方红卫星有限公司 | Relative radiometric correction method for star-load TDICCD camera |
| CN103076097A (en) * | 2013-01-06 | 2013-05-01 | 河北汉光重工有限责任公司 | Reference-source-based segmented linear non-uniform correction method |
| CN103226819A (en) * | 2013-03-29 | 2013-07-31 | 武汉大学 | Segmental counting-based relative radiation correction method |
| CN104050356A (en) * | 2014-04-08 | 2014-09-17 | 北京空间机电研究所 | Scan rate evaluation method for TDI (Transport Driver Interface) infrared detector scanning imaging system |
| US9106857B1 (en) * | 2014-05-09 | 2015-08-11 | Teledyne Dalsa, Inc. | Dynamic fixed-pattern noise reduction in a CMOS TDI image sensor |
-
2016
- 2016-11-18 CN CN201611035496.1A patent/CN106768349B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101226639A (en) * | 2008-01-29 | 2008-07-23 | 航天东方红卫星有限公司 | Relative radiometric correction method for star-load TDICCD camera |
| CN103076097A (en) * | 2013-01-06 | 2013-05-01 | 河北汉光重工有限责任公司 | Reference-source-based segmented linear non-uniform correction method |
| CN103226819A (en) * | 2013-03-29 | 2013-07-31 | 武汉大学 | Segmental counting-based relative radiation correction method |
| CN104050356A (en) * | 2014-04-08 | 2014-09-17 | 北京空间机电研究所 | Scan rate evaluation method for TDI (Transport Driver Interface) infrared detector scanning imaging system |
| US9106857B1 (en) * | 2014-05-09 | 2015-08-11 | Teledyne Dalsa, Inc. | Dynamic fixed-pattern noise reduction in a CMOS TDI image sensor |
Non-Patent Citations (3)
| Title |
|---|
| 刘亚侠: "TDI CCD遥感相机标定技术的研究", 《中国优秀博硕士学位论文全文数据库（博士）信息科技辑》 * |
| 朱建等: "一种基于两点线性校正的改进算法", 《科学技术与工程》 * |
| 王密等: "基于分段辐射校正的星载TDI-CCD成像数据辐射处理方法", 《中国科学：信息科学》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US11604146B2 (en) | 2017-09-19 | 2023-03-14 | Beckman Coulter, Inc. | Analog light measuring and photon counting with a luminometer system for assay reactions in chemiluminescence measurements |
| US11754504B2 (en) | 2017-09-19 | 2023-09-12 | Beckman Coulter, Inc. | System for analog light measuring and photon counting in chemiluminescence measurements |
| CN107907222A (en) * | 2017-11-13 | 2018-04-13 | 国网山东省电力公司荣成市供电公司 | A kind of thermal infrared imaging electric power facility fault locator and detection method |
| CN107907222B (en) * | 2017-11-13 | 2019-10-18 | 国网山东省电力公司荣成市供电公司 | A kind of thermal infrared imaging electric power facility fault locator and detection method |
| CN108181005A (en) * | 2017-11-17 | 2018-06-19 | 天津津航技术物理研究所 | A kind of method and system for the debugging of TDI ccd detectors focal plane |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106768349B (en) | 2019-04-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9787911B2 (en) | Systems and methods for photometric normalization in array cameras | |
| US11448586B2 (en) | Inspection apparatus, sensing apparatus, sensitivity control apparatus, inspection method, and program with pixel sensitivity control | |
| CN110211056B (en) | Self-adaptive infrared image de-striping algorithm based on local median histogram | |
| US8569684B2 (en) | Infrared sensor control architecture | |
| CN100399798C (en) | Apparatus and method for improving image quality in image sensor | |
| CN110686783B (en) | Blind pixel detection method and image processing method for InGaAs short-wave infrared imager | |
| CN103335716B (en) | Area array infrared camera calibration and non-uniformity correction method based on variable integration time | |
| CN106768349A (en) | Adapt to the TDI scanning imaging system relative radiometric calibration methods of multichannel bidirectional output | |
| CN101753797B (en) | Image signal processing apparatus and method | |
| CN103438900A (en) | Three-line-array camera image collaborative absolute radiometric calibration and compensation method | |
| CN102622739A (en) | Image non-uniformity correction method for Bayer filter array color camera | |
| CN110033414B (en) | Infrared image non-uniformity correction method and system based on averaging processing | |
| CN109974854A (en) | A kind of radiation correction method of frame width formula FPI high spectrum image | |
| CN110782403A (en) | Infrared image non-uniformity correction method | |
| US20220187131A1 (en) | Method and system of two-point correction based on temperature substitution | |
| US9143709B1 (en) | Non-uniformity correction (NUC) gain damping | |
| CN107633487A (en) | A kind of system-level relative radiometric correction method of aviation sweep type multispectral scanner image | |
| CN112561829B (en) | Multi-region non-uniform brightness distortion correction algorithm based on L-channel Gamma transformation | |
| CN111609937B (en) | Thermal infrared imager external field calibration method and device | |
| Chakhvashvili et al. | Comparison of reflectance calibration workflows for a UAV-mounted multi-camera array system | |
| US8792612B2 (en) | Radiation imaging system, control method for the same, and program | |
| CN115222619A (en) | Wall surface emissivity solving method based on thermal infrared image correction brightness of thermal imager | |
| CN112882081B (en) | PET detector full-channel gain calibration method | |
| CN111598815B (en) | Thermal infrared imager image non-uniformity correction method and device | |
| Feng et al. | Calculating the leaf-area based on non-loss correction algorithm |
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 |