CN104408709A - Rapid correction method of image with non-uniform gray scale of linear-array CCD camera - Google Patents
Rapid correction method of image with non-uniform gray scale of linear-array CCD camera Download PDFInfo
- Publication number
- CN104408709A CN104408709A CN201410597643.9A CN201410597643A CN104408709A CN 104408709 A CN104408709 A CN 104408709A CN 201410597643 A CN201410597643 A CN 201410597643A CN 104408709 A CN104408709 A CN 104408709A
- Authority
- CN
- China
- Prior art keywords
- image
- array ccd
- ccd camera
- line
- gray
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Transforming Light Signals Into Electric Signals (AREA)
- Image Processing (AREA)
Abstract
The invention relates to a rapid correction method of an image with a non-uniform gray scale of a linear-array CCD camera. The linear-array CCD camera collects an image of a to-be-measured object, the image is corrected line by line, and then a corrected image is formed by splicing and combination; during each line scanning period, gray levels of all pixels of a line image of the to-be-measured object image are arranged at a corresponding linear-array CCD pixel one-dimensional direction to form gray level distribution function; and correction is carried out based on a population-mean-based moving average idea. According to the invention, the correction method has the high real-time performance. Defects of high cost, low precision, and short service life and the like of hardware correction as well as defects of complex algorithm, poor real-time performance, requirement of complex mathematic model establishment, and requirement of parameter calibration of the traditional software correction method can be overcome. Therefore, with comprehensive consideration of the total gray mean value and the local gray mean value, gray compensation and correction of the dark parts of two sides of the CCD imaging are realized, thereby realizing real-time gray correction of the linear-array CCD.
Description
Technical field
The invention belongs to technical field of image processing, relate to the irregular quick antidote of a kind of linear array CCD camera gradation of image.
Background technology
CCD (Charge Couple Device) camera is divided into face battle array and linear array two kinds, and wherein linear array CCD camera is very suitable for the image acquisition of high-speed moving object, is now widely used in the fields such as Machine Vision Detection.But under the demand of high precision and Large visual angle image acquisition, the imaging pixel number of linear array CCD camera is also increase thereupon, the loss that directly results in luminous energy arrival imaging plane significantly increases, make the gray scale of the output image of CCD camera there is the dark irregular phenomenon of illumination in middle bright both sides, seriously have impact on the operations such as follow-up image analysis processing.Therefore irregular to the gray scale exported image carries out corresponding corrective operations and seems very necessary.
In actual applications, as shown in Figure 1, linear array CCD camera imaging pixel arranges along straight line, and due to optical energy loss, the optical axis place energy of the imaging plane of optical imaging system is maximum, and both sides light energy losses is larger.Now, even if adopt parallel uniform source of light to be radiated on measured object, the light intensity arrived on imaging plane also there will be the phenomenon pockety because of optical energy loss, causes collected image intensity value to present the darker irregular result in middle bright both sides, as shown in Figure 2.
At present, the antidote irregular for linear array CCD camera output image gray scale is divided into two kinds: (1) hardware is corrected: carry out sectional-regulated to illumination system, the LED light source generally selecting sectional to regulate carries out illumination compensation to two ends, though this method does not increase the excessive data processing time, there is Automatic adjusument difference and the precision deficiency such as low; (2) software is corrected: mainly contain histogram equalization, based on frequency domain filtering and with reference to antidotes such as standard specimens, these class methods need to set up complicated mathematical model usually, its correlation parameter also needs to demarcate in advance, is difficult to meet the demand such as high real-time and self-adaptation in practical application.
Summary of the invention
The present invention is that the deficiency overcoming the existence of above-mentioned prior art proposes the irregular quick antidote of a kind of linear array CCD camera gradation of image, the method is simple to operate, adaptivity good, without the need to setting up complicated mathematical model and parameter calibration operation, can the intensity profile of efficient, flexible ground orthotic line array CCD camera image irregular.The present invention exports the imaging law of a line view data combining linear array CCD camera each scan period, by the real-time process to progressive image data, propose a kind of irregular antidote of real-time gray scale in conjunction with population mean and the average thought of micro-slip.
The irregular quick antidote of a kind of linear array CCD camera gradation of image of the present invention, it is the good antidote of a kind of real-time, overcoming algorithm complexity, poor real that the high cost that hardware corrects, low precision and the life-span not enough and traditional software antidote such as low exists, complicate mathematical model need set up and carry out on the shortcoming bases such as parameter calibration, by comprehensively taking into account overall gray level average and local gray average, realize correcting the grey level compensation of CCD imaging both sides compared with dark-part, thus achieve the real-time gray-level registration of line array CCD.
The quick antidote that linear array CCD camera imaging gray scale is irregular, linear array CCD camera gathers measured object image, corrects this image line by line, then pieces together the image after rectification; Each line period, the gray-scale value of each pixel of row image of measured object image presses the one-dimensional square of the pixel of corresponding line array CCD to arrangement, forms a grey value profile function, is designated as F (x); Because optical energy loss imaging obviously can present the dark irregular phenomenon of illumination in middle bright both sides, need correct F (x), obtaining the grey value profile function after correcting is F ' (x), and its concrete steps of correcting are as follows:
(1) solve the overall gray level mean value of scan line image corresponding to F (x), and be designated as G
Wherein, N is the number of pixels of current line image;
(2) F (x) is solved in interval
sliding average, and be designated as M (x)
Wherein, m is forward slip step-length, and n slides backward step-length,
total step-length step=m+n+1,0<step<N/4; I is temporary variable,
(3) the grey value profile function F after obtaining rectification is solved ' (x)
F’(x)=G+(F(x)-M(x))。
As above the irregular quick antidote of a kind of linear array CCD camera gradation of image described in 1, m+1=n.
The irregular quick antidote of a kind of linear array CCD camera gradation of image as above, pixel number >=1024 of linear array CCD camera.
Beneficial effect
1, method calculating is simple, real-time is good;
2, compared to the hardware type antidote that traditional employing additional light source compensates, described method has that self-adaptation is good, precision is high and without advantages such as restrictions in serviceable life;
3, method takes into account the inner link of population mean and micro-slip mean value, and its Output rusults has very strong adaptivity.
Accompanying drawing explanation
Fig. 1 is line array CCD optical energy loss schematic diagram
Fig. 2 is the irregular schematic diagram of line array CCD imaging intensity profile
Fig. 3 is the schematic diagram of running mean thought for correcting that the present invention is based on population mean
Fig. 4 is web original image
Web image after correcting when Fig. 5 is step-length step=192
Fig. 6 is the irregular comparison diagram of intensity profile before and after correcting
Embodiment
Below in conjunction with embodiment, set forth the present invention further.Should be understood that these embodiments are only not used in for illustration of the present invention to limit the scope of the invention.In addition should be understood that those skilled in the art can make various changes or modifications the present invention, and these equivalent form of values fall within the application's appended claims limited range equally after the content of having read the present invention's instruction.
The quick antidote that a kind of linear array CCD camera imaging gray scale is irregular, linear array CCD camera gathers measured object image, correct this image line by line, then piece together the image after rectification, Fig. 3 is the schematic diagram of running mean thought for correcting that the present invention is based on population mean; Each line period, the gray-scale value of each pixel of row image of measured object image presses the one-dimensional square of the pixel of corresponding line array CCD to arrangement, forms a grey value profile function, is designated as F (x); Because optical energy loss imaging obviously can present the dark irregular phenomenon of illumination in middle bright both sides, affect image quality, need correct F (x), obtaining the grey value profile function after correcting is F ' (x), and its concrete steps of correcting are as follows:
(1) solve the overall gray level mean value of scan line image corresponding to F (x), and be designated as G
Wherein, N is the number of pixels of current line image;
(2) F (x) is solved in interval
sliding average, and be designated as M (x)
Wherein, m is forward slip step-length, and n slides backward step-length,
total step-length step=m+n+1,0<step<N/4; I is temporary variable,
(3) the grey value profile function F after obtaining rectification is solved ' (x)
F’(x)=G+(F(x)-M(x))。
The irregular quick antidote of a kind of linear array CCD camera gradation of image as above, m+1=n.
The irregular quick antidote of a kind of linear array CCD camera gradation of image as above, pixel number >=1024 of linear array CCD camera.
Embodiment 1
Fig. 4 does not utilize the method to carry out the web original image corrected, and this image is in middle bright as can be seen from this figure, and the irregular distribution characteristics of gray scale that both sides are dark, now collected picture quality is poor.Utilize the method to be corrected below, be that the row image progressive that CCD in each line period collects is corrected in practical application, then piece together a width complete image.Here for be contrasted with former figure, first split line by line by the web original image after collection, then correct line by line, finally piece together a width complete image again, concrete steps are as follows:
The gray-scale value of each pixel of row image of the measured object image after segmentation is pressed the one-dimensional square of the pixel of corresponding line array CCD to arrangement, form a grey value profile function, be designated as F (x), then correct F (x), obtaining the grey value profile function after correcting is F ' (x).
(1) solve the overall gray level mean value of scan line image corresponding to F (x), and be designated as G
Wherein, N is the number of pixels of current line image, now N=6144;
(2) F (x) is solved in interval
sliding average, and be designated as M (x)
Wherein, m is forward slip step-length, and n slides backward step-length,
total step-length step=m+n+1,0<step<N/4; I is temporary variable,
now the size of step=192, m=95, n=96, step depends primarily on the numerical value of m or n, empirically chooses the best in conjunction with measured object Texture eigenvalue, excessive or too smallly all can cause by the change of adopting thing texture.
(3) the grey value profile function F after obtaining rectification is solved ' (x)
F’(x)=G+(F(x)-M(x))。
Irregular distribution function F (x) of gray scale before rectification and irregular distribution function F ' (x) of gray scale after correcting are as shown in Figure 6, can find out that the irregular phenomenon of original gray scale has greatly improved, finally each row Image Mosaic after correcting is become the web image after a complete rectification, as shown in Figure 5.
Claims (3)
1. the quick antidote that linear array CCD camera imaging gray scale is irregular, is characterized in that: linear array CCD camera gathers measured object image, then corrects this image line by line, finally pieces together the image after rectification;
In each line period, the gray-scale value of each pixel of row image of measured object image presses the one-dimensional square of the pixel of corresponding line array CCD to arrangement, forms a grey value profile function, is designated as F (x); Then correct F (x), obtaining the grey value profile function after correcting is F ' (x), and its concrete steps of correcting are as follows:
(1) solve the overall gray level mean value of scan line image corresponding to F (x), and be designated as G
Wherein, N is the number of pixels of current line image;
(2) F (x) is solved in interval
sliding average, and be designated as M (x)
Wherein, m is forward slip step-length, and n slides backward step-length, m,
total step-length step=m+n+1,0<step<N/4; I is temporary variable,
(3) the grey value profile function F after obtaining rectification is solved ' (x)
F’(x)=G+(F(x)-M(x))。
2. the irregular quick antidote of a kind of linear array CCD camera gradation of image according to claim 1, is characterized in that, m+1=n.
3. the irregular quick antidote of a kind of linear array CCD camera gradation of image according to claim 1, is characterized in that, pixel number >=1024 of linear array CCD camera.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410597643.9A CN104408709B (en) | 2014-10-30 | 2014-10-30 | A kind of irregular quick antidote of linear array CCD camera gradation of image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410597643.9A CN104408709B (en) | 2014-10-30 | 2014-10-30 | A kind of irregular quick antidote of linear array CCD camera gradation of image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104408709A true CN104408709A (en) | 2015-03-11 |
CN104408709B CN104408709B (en) | 2017-11-17 |
Family
ID=52646338
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410597643.9A Active CN104408709B (en) | 2014-10-30 | 2014-10-30 | A kind of irregular quick antidote of linear array CCD camera gradation of image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408709B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631808A (en) * | 2015-12-23 | 2016-06-01 | 浙江越剑机械制造有限公司 | Method for correcting image distortion caused by cloth shaking in automatic cloth inspection machine |
CN105741250A (en) * | 2016-02-04 | 2016-07-06 | 东华大学 | Quadratic interpolation based image correction method for automatic cloth inspecting machine with non-uniform cloth travel speed |
CN105761222A (en) * | 2016-02-04 | 2016-07-13 | 东华大学 | Image correction method for non-uniform cloth feeding speed of automatic cloth inspecting machine based on Newton interpolation method |
CN105761221A (en) * | 2016-02-04 | 2016-07-13 | 东华大学 | Image correction method for non-uniform cloth feeding speed of automatic cloth inspecting machine based on natural spline interpolation method |
CN105761220A (en) * | 2016-02-04 | 2016-07-13 | 东华大学 | Image correction method for non-uniform cloth feeding speed of automatic cloth inspecting machine based on linear interpolation method |
CN107610116A (en) * | 2017-09-15 | 2018-01-19 | 天津工业大学 | Conveyer belt On-line Fault Detection inhomogeneous illumination method for correcting image |
CN105741301B (en) * | 2016-02-04 | 2018-10-23 | 东华大学 | The irregular image correction method of automatic cloth inspecting machine walk cloth speed based on Lagrange's interpolation |
CN109615579A (en) * | 2018-11-27 | 2019-04-12 | 成都圭目机器人有限公司 | A kind of illumination compensation algorithm applied to the area array cameras for using LED spotlight |
CN118376157A (en) * | 2024-06-21 | 2024-07-23 | 钛玛科(北京)工业科技有限公司 | Method, device and equipment for detecting size of collimated beam-expanding laser |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070177019A1 (en) * | 2006-01-25 | 2007-08-02 | Jurgen Reichenbach | Method and device for the correction of the brightness of a raw image generated by a sensor matrix |
CN102938137A (en) * | 2012-10-25 | 2013-02-20 | 苏州有色金属研究院有限公司 | Dynamic non-uniformity correction method for linear scanned image based on image sequence analysis |
CN103955897A (en) * | 2014-04-22 | 2014-07-30 | 中国资源卫星应用中心 | CDD image aberration elimination method based on nearest statistics |
-
2014
- 2014-10-30 CN CN201410597643.9A patent/CN104408709B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070177019A1 (en) * | 2006-01-25 | 2007-08-02 | Jurgen Reichenbach | Method and device for the correction of the brightness of a raw image generated by a sensor matrix |
CN102938137A (en) * | 2012-10-25 | 2013-02-20 | 苏州有色金属研究院有限公司 | Dynamic non-uniformity correction method for linear scanned image based on image sequence analysis |
CN103955897A (en) * | 2014-04-22 | 2014-07-30 | 中国资源卫星应用中心 | CDD image aberration elimination method based on nearest statistics |
Non-Patent Citations (4)
Title |
---|
尹楠: "基于区域特征的线阵CCD成像校正方法的研究", 《中国优秀硕士学位论文全文数据库 信息科技辑(月刊)》 * |
王雪晶: "彩色线阵CCD遥感图像的校正", 《中国优秀博硕士学位论文全文数据库(博士) 信息科技辑 (半年刊)》 * |
袁金国: "《遥感图像数字处理》", 28 February 2006 * |
裴益轩 等: "滑动平均法的基本原理及应用", 《火炮发射与控制学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105631808A (en) * | 2015-12-23 | 2016-06-01 | 浙江越剑机械制造有限公司 | Method for correcting image distortion caused by cloth shaking in automatic cloth inspection machine |
CN105631808B (en) * | 2015-12-23 | 2018-10-19 | 浙江越剑智能装备股份有限公司 | A kind of method of pattern distortion caused by cloth cover shake in correction automatic cloth inspecting machine |
CN105761221A (en) * | 2016-02-04 | 2016-07-13 | 东华大学 | Image correction method for non-uniform cloth feeding speed of automatic cloth inspecting machine based on natural spline interpolation method |
CN105761222A (en) * | 2016-02-04 | 2016-07-13 | 东华大学 | Image correction method for non-uniform cloth feeding speed of automatic cloth inspecting machine based on Newton interpolation method |
CN105761220A (en) * | 2016-02-04 | 2016-07-13 | 东华大学 | Image correction method for non-uniform cloth feeding speed of automatic cloth inspecting machine based on linear interpolation method |
CN105741250A (en) * | 2016-02-04 | 2016-07-06 | 东华大学 | Quadratic interpolation based image correction method for automatic cloth inspecting machine with non-uniform cloth travel speed |
CN105761222B (en) * | 2016-02-04 | 2018-10-23 | 东华大学 | The irregular image correction method of automatic cloth inspecting machine walk cloth speed based on Newton interpolating method |
CN105761221B (en) * | 2016-02-04 | 2018-10-23 | 东华大学 | The irregular image correction method of automatic cloth inspecting machine walk cloth speed based on natural spline interpolation method |
CN105761220B (en) * | 2016-02-04 | 2018-10-23 | 东华大学 | The irregular image correction method of automatic cloth inspecting machine walk cloth speed based on linear interpolation method |
CN105741250B (en) * | 2016-02-04 | 2018-10-23 | 东华大学 | The irregular image correction method of automatic cloth inspecting machine walk cloth speed based on quadratic interpolattion |
CN105741301B (en) * | 2016-02-04 | 2018-10-23 | 东华大学 | The irregular image correction method of automatic cloth inspecting machine walk cloth speed based on Lagrange's interpolation |
CN107610116A (en) * | 2017-09-15 | 2018-01-19 | 天津工业大学 | Conveyer belt On-line Fault Detection inhomogeneous illumination method for correcting image |
CN109615579A (en) * | 2018-11-27 | 2019-04-12 | 成都圭目机器人有限公司 | A kind of illumination compensation algorithm applied to the area array cameras for using LED spotlight |
CN109615579B (en) * | 2018-11-27 | 2023-04-14 | 成都圭目机器人有限公司 | Illumination compensation method applied to area-array camera using LED spotlight |
CN118376157A (en) * | 2024-06-21 | 2024-07-23 | 钛玛科(北京)工业科技有限公司 | Method, device and equipment for detecting size of collimated beam-expanding laser |
Also Published As
Publication number | Publication date |
---|---|
CN104408709B (en) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408709A (en) | Rapid correction method of image with non-uniform gray scale of linear-array CCD camera | |
CN105960586B (en) | System and method for image acquisition using supervised high quality imaging | |
CN105046655B (en) | A kind of automatic sharpening method of video image and device | |
CN102438107B (en) | Method for processing microscope video image | |
CN101819024B (en) | Machine vision-based two-dimensional displacement detection method | |
CN102519398B (en) | Optical fiber base pin concentricity measuring method based on visual measurement and system thereof | |
CN103226819A (en) | Segmental counting-based relative radiation correction method | |
CN109872309B (en) | Detection system, method, device and computer readable storage medium | |
CN101865671B (en) | Projection three-dimensional measurement method | |
CN103456021B (en) | A kind of Fabric Defect detection method based on morphological analysis | |
CN104978748A (en) | Liquid crystal screen defect detection method based on local pixel values | |
CN104122271A (en) | Automated optical inspection (AOI)-based bullet apparent defect detection method | |
CN109191520B (en) | Plant leaf area measuring method and system based on color calibration | |
CN102938137A (en) | Dynamic non-uniformity correction method for linear scanned image based on image sequence analysis | |
CN103698341A (en) | System for detecting cotton rolling quality and method for detecting cotton rolling quality based on image | |
CN102306378B (en) | Image enhancement method | |
CN104504722A (en) | Method for correcting image colors through gray points | |
CN113487563B (en) | EL image-based self-adaptive detection method for hidden cracks of photovoltaic module | |
CN104021531A (en) | Improved method for enhancing dark environment images on basis of single-scale Retinex | |
CN103475821A (en) | Adjustment method based on automatic integration time of near infrared camera | |
CN104361566B (en) | A kind of image processing method in optimization dark portion region | |
CN107133937B (en) | A kind of self-adapting enhancement method of infrared image | |
CN103514588B (en) | Image enchancing method and system | |
CN103500457A (en) | Method of color cast detection of video image | |
CN104820970A (en) | Infrared image relative radiation correction method based on on-orbit classified statistic |
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 |