CN104143180B - A kind of many alignment moveout scan subpixel image non-uniform correction methods - Google Patents
A kind of many alignment moveout scan subpixel image non-uniform correction methods Download PDFInfo
- Publication number
- CN104143180B CN104143180B CN201410320930.5A CN201410320930A CN104143180B CN 104143180 B CN104143180 B CN 104143180B CN 201410320930 A CN201410320930 A CN 201410320930A CN 104143180 B CN104143180 B CN 104143180B
- Authority
- CN
- China
- Prior art keywords
- image
- width
- row
- subpixel
- detector array
- 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.)
- Active
Links
Landscapes
- Image Processing (AREA)
Abstract
A kind of many alignment moveout scan subpixel image non-uniform correction methods, (1) construct many alignment moveout scan detection devices；(2) each detector array is scanned imaging according to the sampling interval in scanning direction, and imaging every time respectively obtains NtGroup view data, proceeds to step (3) immediately；(3) respectively by the N of each detector arraytGroup view data forms a frame detection image after being processed；(4) splicing of frame detection image is obtained two width subpixel images；(5) two width subpixel images are carried out intersection by row and is spliced to form image Ip1；(6) obtain arranging the image I to after nonuniformity correctionp2；(7) rebuild two width to complete to arrange the subpixel image to Nonuniformity Correction；(8) two width subpixel images are rotated by 90 ° at same direction respectively, obtain image Ip3；(9) two width are rebuild and completes subpixel image of the row to nonuniformity correction.(10) above-mentioned two width subpixel image opposite direction is rotated by 90 °.
Description
Technical field
The invention belongs to image processing field, is related to a kind of many alignment moveout scan subpixel image nonuniformity correction sides
Method.
Background technology
The application of remote sensing images is more and more extensive, but due to being subject to sensitive detection parts itself, process conditions, extraneous input etc.
The impact of composite factor, the output responsibility of each detection unit are not quite identical, will appear as regular band on image
And blind element, i.e. heterogeneity phenomenon.Heteropical occur generating picture quality very big impact, and then have impact on follow-up
Image application.Therefore, the heterogeneity of detector is always imaging system needs one of major issue of solution.General employing
Multinomial be fitting process, Nonuniformity Correction is carried out based on the method such as statistic law, eliminate or weaken heterogeneity phenomenon, be easy to follow-up
Data processing.For many alignment moveout scan images, in addition to it will eliminate the heterogeneity of each alignment image, to also
The gray scale difference for reducing Whole Response difference between alignment image as far as possible and causing, consequently facilitating follow-up frame matching etc. is processed.
Content of the invention
It is an object of the invention to provide a kind of heterogeneity school suitable for many alignment moveout scan sub-pixed mapping composographs
Correction method, solves the problems, such as the heterogeneity phenomenon of many alignment images, while the overall intensity eliminated between alignment image is poor.
The present invention technical solution be：A kind of many alignment moveout scan subpixel image non-uniform correction methods, step
Rapid as follows：
(1) many alignment moveout scan detection devices are constructed, and the device includes the detection of optical system, sweep mechanism and many alignments
Device；Described sweep mechanism includes putting mirror and its drive shaft；Described many detector array be two detector array, alignment visit
Survey device and adopt NtIndividual detection array composition, the corresponding instantaneous field of view of pixel is IFOV, and two neighboring detection array is arranged in parallel,
Vertical scanning direction staggers 1/N successivelytIndividual pixel, and arrange detection array scanning direction, in a sampling length sample St
Secondary；The sampling length is the corresponding instantaneous field of view of pixel；Described NtIt is more than or equal to 2；The StSpan is St≥2；
(2) scene in visual field is imaged in focal plane by optical system and sweep mechanism together, and drive shaft drives pendulum mirror rotation
Turn, the scene imaging in linear field with tandem two detector array on the inswept focal plane of certain speed, two
Detector array is successively imaged to same position scene in visual field, imaging time intervalIt is each alignment detection
N in devicetIndividual detection array is imaged simultaneously, obtains NtGroup view data；Step (3) is proceeded to afterwards immediately；Each detector array
Be scanned imaging according to the sampling number corresponding sampling interval is arranged in step (1) in scanning direction, imaging every time is obtained respectively
Arrive NtGroup view data, proceeds to step (3) immediately after obtaining data；
(3) respectively by the N of two detector arraytGroup view data forms a frame detection figure after carrying out alignment splicing
Picture；
(4), after completing default sampling number in a scanning direction, the frame that each detector array is correspondingly obtained is detected
Image carries out splicing according to the time and obtains two width subpixel images；
(5) two width subpixel images are carried out intersection splicing by row, forms the new stitching image I of a widthp1；
(6) to stitching image Ip1Every string image carry out following process respectively, obtain arranging after the completion of all column processing to
Image I after nonuniformity correctionp2；
(6.1) I is obtainedp1Per string image IiStatistic histogram Hi；
(6.2) for Ip1In every string image Ii, according to IiAnd its neighborhood NmThe rectangular histogram of row image, calculates Gauss and adds
Power expects rectangular histogram HGi；
(6.3) according to HiAnd HGiTo IiCarry out histogram specification conversion；
(7) row according to two width images in step (5) intersect the order of splicing, from Ip2The middle row accordingly that extract respectively are schemed
Picture, rebuilds two width and completes to arrange the subpixel image to Nonuniformity Correction；
(8) row will be completed to be rotated by 90 ° respectively to two width subpixel images after nonuniformity correction at same direction, repeats to walk
Suddenly (5)～(6), image I after being correctedp3；
(9) row according to two width images in step (5) intersect the order of splicing, from Ip3The middle row accordingly that extract respectively are schemed
Picture, rebuilds two width and completes subpixel image of the row to nonuniformity correction；
(10) row will be completed to be rotated by 90 ° to two width subpixel images of nonuniformity correction by the opposite direction that step (8) rotates,
Obtain the two width subpixel images for completing row, column both direction nonuniformity correction.
Present invention advantage compared with prior art is：
(1) this method considers the Fringe Characteristics of noise, and image is processed by column, has considered when prostatitis and week
Enclose six neighborhood row half-tone informations.The characteristics of for fringes noise, Local treatment is carried out to image by the way of row are divided, hard
In terms of part realization, in order to improve processing speed, can be with parallel processing.
(2) as the explorer response of many detector array is inconsistent, two alignment institutes are caused to there is overall ash into image
Degree is poor, treated by the present method after, the gray scale difference between alignment image can be effectively reduced so that many alignment images are in same ash
Degree scope, is conducive to the precision for improving subsequent treatment.
(3) method proposed by the present invention can be used for any two alignment or multiple alignment figures in multi-thread scanning imaging system
The Nonuniformity Correction of picture.
Description of the drawings
Fig. 1,2 are many alignment moveout scan detection device two ways schematic diagrams of the invention；
Fig. 3 is the detector array two field picture splicing schematic diagram of the present invention；
Fig. 4 is that two width sub-pixed mapping of many alignment moveout scans row of the present invention intersect splicing schematic diagram；
Specific embodiment
Below in conjunction with the accompanying drawings and example elaborates to the present invention.A kind of many alignment moveout scan subpixel images are non-
Even bearing calibration, step are as follows：
(1) many alignment moveout scan detection devices are constructed, and the device includes that optical system 1, sweep mechanism 2 and many alignments are visited
Survey device 3；Described Scan Architecture includes putting mirror and its drive shaft；Described many detector array be two detector array, line
Row detector adopts NtIndividual detection array composition, the corresponding instantaneous field of view of pixel is IFOV, two neighboring detection array parallel
Row, stagger 1/N successively in vertical scanning directiontIndividual pixel, and detection array is set in scanning direction, in a sampling length
Sampling StSecondary；The sampling length is the corresponding instantaneous field of view of pixel；Many alignment moveout scan detection device sweep mechanisms
Angular scanning speed beApart from d between two of which detector array, the target minimum movement speed of detection
Degree vmin, optical system focal length f, the ground sampled distance GSD of detector array；Described NtIt is more than or equal to 2；The StValue model
Enclose for St≥2；Below with NtIllustrate as a example by=2.
What Fig. 1 was given is front end scanning probe device；Incident illumination comprising target and the emittance information of background is through pendulum
Focal plane is converged to through optical system 1 after mirror reflection, the picture of scenery is formed, drive shaft drives pendulum mirror according to default angular speed
Rotation, makes the picture of scenery each detector array inswept successively.When the picture of scenery is with the inswept one of line of certain speed
During row detector, detector is sampled to the picture of scenery.What Fig. 2 was given is rear-end scanning detection device.Comprising target and the back of the body
The incident illumination of the emittance information of scape converges to pendulum mirror through optical system 1, reflexes to focal plane through putting mirror, forms scenery
Picture.Drive shaft drives pendulum mirror to rotate according to default angular speed, makes the picture of scenery each detector array inswept successively.Work as scape
During the picture of thing is inswept with certain speed one of detector array, detector is sampled to the picture of scenery.
In this example, optical system 1 is the optical system of typical Cassegrain form, is made up of primary mirror and secondary mirror, incident illumination
Line is incided on detector array after primary mirror and secondary mirror reflecting focal.
(2) scene in visual field is imaged in focal plane by optical system 1 and sweep mechanism 2 together, and drive shaft drives pendulum mirror
Rotation, the scene imaging in linear field with tandem two detector array on the inswept focal plane of certain speed, two
Individual detector array is successively imaged to same position scene in visual field, imaging time intervalIt is that each alignment is visited
The N surveyed in devicetIndividual detection array is imaged simultaneously, obtains NtGroup view data；Step (3) is immediately entered after obtaining view data；In
This simultaneously, each detector array still according to arranging in the sampling length sampling number in step (1), in scanning side
To imaging is scanned, imaging every time respectively obtains NtGroup view data, proceeds to step immediately after being similarly obtained view data
(3)；For example, S can be sett>=2, then can achieve detector array more than 2 times over-samplings in a scanning direction；
(3) respectively by the N of two detector arraytGroup view data forms a frame detection figure after carrying out alignment splicing
Picture, splicing is as shown in figure 3, by NtGroup pattern intersects splicing.Sub-pixed mapping frame detection image is obtained by splicing,
Realize stretching of the target in vertical scanning direction.
(4) after default sampling number, the frame detection image that each detector array is correspondingly obtained entered according to the time
Row splicing obtains two width subpixel images；Default sampling number can use 200～300 rows, and increasing default sampling number can increase
Scanogram details, improves the precision that successive image is processed；The efficiency that default sampling number can improve data processing is reduced, because
This default sampling number can be adjusted according to practical situation.
(5) two width subpixel images are carried out intersection splicing by row, forms the new stitching image I of a widthp1, two width sub-pixed mappings
The row of image intersect splicing as shown in figure 4, being two width subpixel images above in figure, lower section is spliced image；
(6) I is obtainedp1The statistic histogram H of each columni；
(7) for Ip1In every string image Ii, according to IiAnd its neighborhood NmThe rectangular histogram of (general 8-12 row) row image,
Calculate Gauss weighting and expect rectangular histogram HGi, such as formula (1),
H in formula (1)GiExpectation rectangular histogram for the i-th row image；Hi+j, j ∈ [- Nm,Nm] it is Ip1In i-th row and its surrounding Nm
The grey level histogram of neighborhood image；(σ is j) gaussian weighing function, can be calculated by formula (2) g：
In formula (2), σ is Gaussian function standard deviation, and x is distance of other image row to present image row；
(8) according to HiAnd HGiTo IiCarry out histogram specification conversion；
(9) by Ip1Every string image repeat step (6)～(8), obtain arranging the image I to after nonuniformity correctionp2；
(10) row according to two width images in step (5) intersect the order of splicing, from Ip2The middle row accordingly that extract respectively are schemed
Picture, rebuilds two width and completes to arrange the subpixel image to Nonuniformity Correction；
(11) row will be completed to be rotated by 90 ° respectively to two width subpixel images after nonuniformity correction at same direction, is repeated
Step (5)～(9), image I after being correctedp3；
(12) row according to two width images in step (5) intersect the order of splicing, from Ip3The middle row accordingly that extract respectively are schemed
Picture, rebuilds two width and completes subpixel image of the row to nonuniformity correction.
(13) row will be completed to rotate to two width subpixel images of nonuniformity correction by the opposite direction that step (11) rotates
90 °, obtain the two width subpixel images for completing row, column both direction nonuniformity correction.
Unspecified part of the present invention belongs to general knowledge as well known to those skilled in the art.
Claims (1)
1. a kind of many alignment moveout scan subpixel image non-uniform correction methods, it is characterised in that step is as follows：
(1) many alignment moveout scan detection devices are constructed, and the device includes optical system, sweep mechanism and many detector array；
Described sweep mechanism includes putting mirror and its drive shaft；Described many detector array be two detector array, alignment detect
Device adopts NtIndividual detection array composition, the corresponding instantaneous field of view of pixel is IFOV, and two neighboring detection array is arranged in parallel, is hanging down
Straight scanning direction is staggered 1/N successivelytIndividual pixel, and arrange detection array scanning direction, in a sampling length sample St
Secondary；The sampling length is the corresponding instantaneous field of view of pixel；Described NtIt is more than or equal to 2；The StSpan is St≥2；
(2) scene in visual field is imaged in focal plane by optical system and sweep mechanism together, and drive shaft drives pendulum mirror rotation, line
Scene imaging in visual field is visited with tandem two detector array on the inswept focal plane of certain speed, two alignments
Survey device to be successively imaged same position scene in visual field, imaging time intervalIt is in each detector array
NtIndividual detection array is imaged simultaneously, obtains NtGroup view data；Obtain NtThe place of step (3) is immediately performed after group view data
Reason；And obtaining NtAfter group view data, according to arranging in step (1), sampling number is corresponding to be adopted each detector array
Sample is spaced in scanning direction and is scanned imaging, and imaging every time respectively obtains NtGroup view data, is immediately performed after obtaining data
The process of step (3)；
(3) respectively by the N of two detector arraytGroup view data forms a frame detection image after carrying out alignment splicing；
(4) after completing default sampling number in a scanning direction, frame detection image that each detector array is correspondingly obtained
Splicing is carried out according to the time and obtains two width subpixel images；
(5) two width subpixel images are carried out intersection splicing by row, forms the new stitching image I of a widthp1；
(6) to stitching image Ip1Every string image carry out following process respectively, obtain arranging after the completion of all column processing to non-
Image I after even correctionp2；
(6.1) I is obtainedp1Per string image IiStatistic histogram Hi；
(6.2) for Ip1In every string image Ii, according to IiAnd its neighborhood NmThe rectangular histogram of row image, calculates the Gauss weighting phase
Hope rectangular histogram HGi；
(6.3) according to HiAnd HGiTo IiCarry out histogram specification conversion；
(7) row according to two width images in step (5) intersect the order of splicing, from Ip2Middle extract corresponding row image respectively, weight
Build two width to complete to arrange the subpixel image to Nonuniformity Correction；
(8) row will be completed to be rotated by 90 ° respectively to two width subpixel images after nonuniformity correction at same direction, repeat step
(5)～(6), image I after being correctedp3；
(9) row according to two width images in step (5) intersect the order of splicing, from Ip3Middle extract corresponding row image respectively, weight
Build two width and complete subpixel image of the row to nonuniformity correction；
(10) row will be completed to be rotated by 90 ° by the opposite direction that step (8) rotates to two width subpixel images of nonuniformity correction, is obtained
Complete two width subpixel images of row, column both direction nonuniformity correction.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410320930.5A CN104143180B (en) | 2014-07-04 | 2014-07-04 | A kind of many alignment moveout scan subpixel image non-uniform correction methods |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410320930.5A CN104143180B (en) | 2014-07-04 | 2014-07-04 | A kind of many alignment moveout scan subpixel image non-uniform correction methods |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104143180A CN104143180A (en) | 2014-11-12 |
| CN104143180B true CN104143180B (en) | 2017-03-15 |
Family
ID=51852348
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410320930.5A Active CN104143180B (en) | 2014-07-04 | 2014-07-04 | A kind of many alignment moveout scan subpixel image non-uniform correction methods |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104143180B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106408553B (en) * | 2015-07-29 | 2019-10-22 | 北京空间飞行器总体设计部 | The infrared linear array detector target response analysis method tiltedly swept |
| CN110580692B (en) * | 2019-09-11 | 2022-03-25 | 北京空间飞行器总体设计部 | Method for correcting radiation consistency of multi-line time difference scanning image |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101893552A (en) * | 2010-07-06 | 2010-11-24 | 西安电子科技大学 | Hyperspectral imager and imaging method based on compressive sensing |
| CN103607547A (en) * | 2013-12-09 | 2014-02-26 | 江苏思特威电子科技有限公司 | Pixel imaging device for mirror image and imaging method for mirror image |
-
2014
- 2014-07-04 CN CN201410320930.5A patent/CN104143180B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101893552A (en) * | 2010-07-06 | 2010-11-24 | 西安电子科技大学 | Hyperspectral imager and imaging method based on compressive sensing |
| CN103607547A (en) * | 2013-12-09 | 2014-02-26 | 江苏思特威电子科技有限公司 | Pixel imaging device for mirror image and imaging method for mirror image |
Non-Patent Citations (1)
| Title |
|---|
| 增强现实运动头部目标跟踪中的误差消除方法;王卉;《计算机仿真》;20140531;第428-431页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104143180A (en) | 2014-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102360079B (en) | Laser range finder and working method thereof | |
| CN103592030A (en) | Snapshot integral view field imaging full-polarization hyperspectral detection device | |
| TWI480578B (en) | Method for detecting optical center of wide-angle lens, and optical center detection apparatus | |
| US10462344B2 (en) | Imaging system having multiple imaging sensors and an associated method of operation | |
| CN106679676B (en) | A kind of monoscopic multifunctional optical sensor and implementation method | |
| CN107341780A (en) | A kind of Infrared images pre-processing bearing calibration | |
| CN102509304A (en) | Intelligent optimization-based camera calibration method | |
| US20240169505A1 (en) | Optical fiber image calibration method and apparatus, and computer- readable storage medium | |
| JP2017184217A5 (en) | ||
| CN104143180B (en) | A kind of many alignment moveout scan subpixel image non-uniform correction methods | |
| WO2020118896A1 (en) | Image repair method, image repair system, and flat panel detector | |
| CN103930773A (en) | Method for correcting timing skew in X-ray data read out of X-ray detector in rolling shutter mode | |
| CN105466563A (en) | Infrared thermal imager inhomogeneity two-point correction module and method | |
| CN104159007A (en) | Point target searching detection device based on multi-linear time lag extension sampling | |
| CN104156977A (en) | Point target movement velocity detection method based on multiple linear moveout scanning, extending and sampling | |
| Urban et al. | LaFiDa—A laserscanner multi-fisheye camera dataset | |
| CN103004322B (en) | Corn ear trait detecting device | |
| US11902496B2 (en) | Depth sensing camera system | |
| CN104143179B (en) | Method for enhancing moving target through multi-linear-array time difference scanning expansion sampling | |
| CN108742531B (en) | Imaging correction method based on large-range OCT scanning | |
| CN108507539B (en) | Optical camera single-line array push-broom mode equal ground resolution imaging method | |
| WO2023005942A1 (en) | Laser scanning imaging method and system, storage medium, and computer program | |
| CN104143196B (en) | A kind of point target detection method based on many alignment moveout scans extension sampling | |
| CN104143187B (en) | A kind of many alignment moveout scan extension sampling subpixel image method for registering | |
| Eglitis et al. | Baldone 1.2 m telescope plate archive-Hidden reserves of the FON project |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |