CN102397075B - Method for eliminating chromatic aberration caused by imaging environment and inspecting stability of imaging environment - Google Patents

Method for eliminating chromatic aberration caused by imaging environment and inspecting stability of imaging environment Download PDF

Info

Publication number
CN102397075B
CN102397075B CN2010102807248A CN201010280724A CN102397075B CN 102397075 B CN102397075 B CN 102397075B CN 2010102807248 A CN2010102807248 A CN 2010102807248A CN 201010280724 A CN201010280724 A CN 201010280724A CN 102397075 B CN102397075 B CN 102397075B
Authority
CN
China
Prior art keywords
value
image
axle
color
formula
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.)
Expired - Fee Related
Application number
CN2010102807248A
Other languages
Chinese (zh)
Other versions
CN102397075A (en
Inventor
江炯宏
张恒鸿
林康平
陈瑞照
Original Assignee
江炯宏
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 江炯宏 filed Critical 江炯宏
Priority to CN2010102807248A priority Critical patent/CN102397075B/en
Publication of CN102397075A publication Critical patent/CN102397075A/en
Application granted granted Critical
Publication of CN102397075B publication Critical patent/CN102397075B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)
  • Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)

Abstract

The invention discloses a method for eliminating chromatic aberration caused by an imaging environment and inspecting stability of the imaging environment. A method for inspecting errors caused by an environment factor and a personal factor is established. The method comprises the following steps of: setting an image capturing range, capturing an image, converting a color space, calculating an average color length value, calculating a rating value X and judging the influence degree of variance, and eliminating the chromatic aberration caused by the environment factor and the personal factor by using a specified measurement distance and angle. By the method, the difference degree of the chromatic aberration caused by the environment factor and the personal factor during pick-up can be displayed under the condition of the same object to be measured and imaging environment so as to judge whether the stability of the imaging environment meets the requirement or not.

Description

Get rid of the method that imaging circumstances causes aberration and check imaging circumstances degree of stability
Technical field
The present invention is a kind of method that promotes the imaging circumstances degree of stability, and espespecially a kind of various factors of imaging circumstances of getting rid of cause the method for aberration.
Background technology
21st century diagnostic imaging medical science development, be high-resolution video diagnostic technology, the medical environment of high-speed image Transfer Technology and more convenience-for-people and environmental protection, not only there is the digitized X-ray photography to replace traditional shadowgraph, to promote the quality of image and to reduce patient's radiation dose, high-resolution color digital image acquiring technology is occupied considerable status in application medical science, especially the Evolution of Tongue Inspection of TCM instrument is the tongue color of utilizing the high-resolution color digital camera shooting patient of modernization science and technology, the numerical data that obtains can be utilized the computer image treatment technology, revise its color contrast scope or utilize different compute modes to strengthen raw video, so that being provided, the accurate image of doctor inspects pattern, the doctor only need to input patient's data such as case history number before terminating machine, just can make diagnosis according to this inspection image immediately, so can reduce patient's waiting time, to promote work efficiency and to reduce human cost.
But while utilizing digital imaging device acquisition image, may exist the error that gauge itself, anthropic factor or other environmental factors produce, for example same operator carry out the multiimage while capturing in same place in same sample, this repetitive still can produce the variation on numerical value, and this variation is survey tool itself causes when measuring variation; And the different operating personnel are when same place is measured same sample, the variation that the measurement variation of gained causes for artificial factor, or be subject to the variation of form and aspect, lightness and the chroma of other environmental factors such as light when image capture, abovely all can produce aberration and cause erroneous judgement.
Therefore, the problems referred to above, will be at the people of this art where the shoe pinches to be solved.
Summary of the invention
Because the described deficiency of above-mentioned prior art, main purpose of the present invention is to provide a kind of imaging circumstances of getting rid of to cause the method for aberration and check imaging circumstances degree of stability, can check and get rid of imaging circumstances and cause the factor of aberration.
For achieving the above object, the invention provides a kind of imaging circumstances of getting rid of and cause the method for aberration, it can get rid of the error that envirment factor or anthropic factor produce, and its method step is:
(1) the capturing images scope is set, and image capture unit is according to two above anchor points setting up on determinand, shooting distance and the angle of calibrating this device and determinand;
(2) capturing images, be to utilize image capture unit acquisition fixedly the image vision range of determinand is to obtain this pixel rgb value, and same area repeats to capture the action of image on determinand, to obtain the pixel rgb value of same position more than two groups;
(3) color space conversion, convert the pixel rgb value of same position more than two groups of this image to the Lab color space by one group of color conversion formula;
(4) calculate the average color length value, utilize a formula to ask the color length value of each position, its formula is as follows:
Figure BSA00000268391100021
Again the average color length value that represents this image is tried to achieve in the Δ E addition of each position
Figure BSA00000268391100022
Its
Figure BSA00000268391100023
Necessary≤4, it is acceptable average color length value scope;
(5) calculate comparation and assessment value X, this computing formula is:
Figure BSA00000268391100024
σ wherein 2=σ 1 2+ σ 2 2
σ wherein 2For total variation, σ 1 2Variance, σ 2 for the envirment factor generation 2For the variance that artificial error produces, A is the yield control constant, and it is worth between 1 to 7, and the impact that the larger expression envirment factor of X value and personal error produce is larger;
(6) judge the influence degree of variance, when comparation and assessment value X is less than or equal to 10%, can get rid of the error that envirment factor or anthropic factor produce with this understanding, if comparation and assessment value X, greater than 10%, must re-start the capturing images scope, set.
The present invention also provides a kind of method of checking the imaging circumstances degree of stability, and its method step is:
(1) image values acquisition, be to utilize the image of the fixing object to be measured of image capture unit acquisition to obtain its pixel rgb value, and same area repeats to capture the action of image on object to be measured, to obtain the pixel rgb value of same position more than two groups;
(2) color space conversion, convert the pixel rgb value of same position more than two groups of this image to the Lab color space by one group of color conversion formula;
(3) calculate the average color length value, utilize a formula to ask the color length value of each position, its formula is as follows:
Figure BSA00000268391100031
Again the average color length value that represents this image is tried to achieve in the color length value Δ E addition of each position
Figure BSA00000268391100032
Its average color length value For acceptable scope;
(4) calculate comparation and assessment value X, this computing formula is:
Wherein, σ 2=σ 1 2+ σ 2 2
σ wherein 2For total variation, σ 1 2Variance, σ 2 for the envirment factor generation 2For the variance that artificial error produces, A is the yield control constant, and it is worth between 1 to 7, if the impact that the larger expression envirment factor of X value and personal error produce is larger.
Wherein image capture unit is photographing unit, ccd camera or network video device.
due to image capture unit in acquisition during image, being subject to envirment factor and anthropic factor impact, can to produce color poor, this envirment factor may comprise the variation of image capture unit itself, the color phase of external light, lightness and chroma, the distance of object to be measured and image capture unit and angle, and anthropic factor may comprise operator produce when measuring variation, focusing, the phenomenons such as hands shake, can judge the error degree that environmental factors or anthropic factor produce image according to the size that comparation and assessment are worth X by method of the present invention, and then while get rid of taking because of the variation of envirment factor and anthropic factor, obtain the image close with former object color.
Description of drawings
Fig. 1 is that the present invention gets rid of the flow chart one that imaging circumstances causes aberration and checks the method for imaging circumstances degree of stability;
Fig. 2 is that the present invention gets rid of the flowchart 2 that imaging circumstances causes aberration and checks the method for imaging circumstances degree of stability.
Description of reference numerals
A1~A6: flow chart S1~S4: flow chart
The specific embodiment
Reach in order to achieve the above object effect, the accompanying drawing of now arranging in pairs or groups, the preferred embodiment that develops simultaneously, so that the people of the common knowledge of tool can be with regard to the enforcement according to this of projects in this field.
At first, see also Fig. 1 and Fig. 2, it causes the flow chart of the method for aberration and check imaging circumstances degree of stability for the present invention gets rid of imaging circumstances, know in figure and point out that a kind of imaging circumstances of getting rid of causes the method for aberration and check imaging circumstances degree of stability, the error that is produced to control imaging circumstances and anthropic factor, its method step is:
(1) the capturing images scope is set, and image capture unit can be according to two above anchor points setting up on determinand, shooting distance and the angle of calibrating this device and determinand;
(2) image values acquisition, to utilize a high-resolution color digitized video photographing unit, ccd (Charge-coupled Device, charge coupled cell) image of camera or network video device acquisition determinand, with pixel RGB (RGB) value that obtains this image, and same area repeats to capture image on determinand, to obtain the pixel rgb value of same position more than two groups;
(3) color space conversion, convert the pixel rgb value of same position more than two groups of this image to the Lab color space by one group of color conversion formula and (with three mutually perpendicular coordinate axess, represent a color space, the L axle represents lightness, the a axle is red-green axle, the b axle is Huang-blue axle), this group formula comprises the pixel rgb value is converted to the XYZ value, and the XYZ value converts f to x, f y, f zAfter be converted to again the Lab color space, its color conversion formula is as follows:
(a) first the pixel rgb value is converted to the XYZ value, the pixel rgb value calculates respectively the XYZ value according to the parameter of following different proportion:
X=0.490R+0.310G+0.200B,Y=0.177R+0.812G+0.011B,
Z=0.010G+0.990B;
(b) the XYZ value converts the Lab color space to, and its method is as follows:
(b-1) try to achieve x ry rz r, at first setting a reference white color numerical value is (X rY rZ r), then with the following formula of XYZ value substitution of step (a) calculating,
x r = X X r , y r = Y Y r , Z r = Z Z r ,
Make its (X, Y, Z) can be exchanged into (x ry rz r);
(b-2) try to achieve f x, f y, f z, with above-mentioned (x ry rz r) the following formula of according to condition substitution to be to produce (f x, f y, f z)
f x = x r 3 x r > 0.00886 kx r + 16 116 x r ≤ 0.00886 k = 903.3
f y = y r 3 y r > 0.00886 ky r + 16 116 y r ≤ 0.00886 k = 903.3
f z = z r 3 z r > 0.00886 kz r + 16 116 z r ≤ 0.00886 k = 903.3
(b-3) try to achieve (L, a, b) value, finally again with (f x, f y, f z) change into (L, a, b) value by following formula
L=116f y-16,a=500(f x-f y),b=200(f y-f z)
Can convert every group of XYZ value to numerical value (L, a, b) by above-mentioned formula;
(4) calculate the average color length value, utilize a formula to ask the color length of identical measurement point, its formula is as follows:
Its (L 0, a 0, b 0)
Be a reference value, (L, a, b) is other measured value, then the average color length value is tried to achieve in the color length Δ E addition of diverse location
Figure BSA00000268391100062
To represent this image, its average color length value
Figure BSA00000268391100063
Necessary≤4, it is acceptable scope;
(5) calculate comparation and assessment value X, this computing formula is:
Figure BSA00000268391100064
Wherein, σ 2=σ 1 2+ σ 2 2
σ wherein 2For total variation, σ 1 2Variance, σ 2 for the envirment factor generation 2For the variance that artificial error produces, A is the yield control constant, and it is worth between 1 to 7, if the impact that the larger expression envirment factor of X value and personal error produce is larger;
(6) judge the influence degree of variance, when comparation and assessment value X is less than or equal to 10%, can get rid of the error that envirment factor or anthropic factor produce with this understanding, if comparation and assessment value X, greater than 10%, must re-start the capturing images scope, set.
A kind ofly check the method practical operation of imaging circumstances degree of stability as follows, selected 10 samples, and 3 survey crews of assigned at random are measured, after everyone is measured its result after 100 times machine calculates as calculated the average color length value insert the form that designs in advance with 3 operators everyone with every 6 seconds once frequency and captured image in 10 minutes, take with cable release or computer control, ANOVA (variance analysis) analytical table is as shown in table one as calculated.
Table one
The ANOVA analytical table
Figure BSA00000268391100065
N (sample number)=10, p (operational personnel)=3, k (number of repetition)=10,
σ 1 2 = σ R 2 = MS R = 0.001992 ;
σ 2 2 = σ o 2 = ( MS O - MS PO ) / n × k = 0.009876 ;
σ 2=σ 1 22 2=0.011868;
The A coefficient is made as (being that fraction defective is below 1%) substitution formula under 5.15 specification and calculates in the hope of comparation and assessment value X;
X = 5.15 σ 4 × 100 % = 14.02576 % .
A kind of imaging circumstances of getting rid of causes the method for aberration, adopt above-mentioned steps and by formula, calculate comparation and assessment value X, when comparation and assessment value X is less than or equal to 10%, can get rid of the error that imaging circumstances or anthropic factor produce with this understanding,, if X is greater than 10% for the comparation and assessment value, must re-starts the capturing images scope and set; The present embodiment calculates comparation and assessment value X=14.02576%>10% through formula,, therefore must get back to step 1, re-starts the capturing images scope and sets, until comparation and assessment value X, less than 10%, sets with the relevant parameter of getting rid of aberration that imaging circumstances is produced.
The described determinand of the present embodiment, be the three-dimensional prosthese that constant color and labelling are coated in a surface, it is characterized in that using a kind of imaging circumstances of getting rid of to cause the method for aberration and check imaging circumstances degree of stability to carry out a series of repetition captures in order to analyze the degree of stability of imaging circumstances for this solid prosthese.
Comprehensively above-mentioned, be more than one of them preferred embodiment of the present invention, not in order to limit to protection scope of the present invention, protection scope of the present invention must be defined and is as the criterion with the scope of aforementioned claim.

Claims (3)

1. get rid of imaging circumstances and cause the method for aberration for one kind, it can get rid of the error that envirment factor or anthropic factor produce, and it is characterized in that, its method step is:
(1) the capturing images scope is set, and image capture unit is according to two above anchor points setting up on determinand, shooting distance and the angle of calibrating this device and determinand;
(2) capturing images, be to utilize image capture unit acquisition fixedly the image of determinand is to obtain its pixel rgb value, and same area repeats to capture the action of image on determinand, to obtain the pixel rgb value of same position more than two groups;
(3) color space conversion, convert the pixel rgb value of same position more than two groups of this image to the Lab color space by one group of color conversion formula;
(4) calculate the average color length value, utilize a formula to ask the color length value of each position, its formula is as follows:
ΔE = ( L - L 0 ) 2 + ( a - a 0 ) 2 + ( b - b 0 ) 2 , Again the average color length value that represents this image is tried to achieve in the Δ E addition of each position
Figure FDA00003003411200012
Its
Figure FDA00003003411200013
Necessary≤4, it is acceptable average color length value scope, and wherein Δ E is color length, and L represents lightness, L 0Be the reference value of lightness, a is red-green axle, a 0Be the reference value of red-green axle, the b axle is Huang-blue axle, b 0Axle is the reference value of Huang-blue axle;
(5) calculate comparation and assessment value X, computing formula is:
X = Aσ 4 × 100 % , σ wherein 2=σ 1 2+ σ 2 2
σ wherein 2For total variation, σ 1 2Variance, σ 2 for the envirment factor generation 2For the variance that artificial error produces, A is the yield control constant, and it is worth between 1 to 7, and the impact that the larger expression envirment factor of X value and personal error produce is larger;
(6) judge the influence degree of variance, when comparation and assessment value X is less than or equal to 10%, can get rid of the error that envirment factor or anthropic factor produce with this understanding, if comparation and assessment value X, greater than 10%, must re-start the capturing images scope, set.
2. a method of checking the imaging circumstances degree of stability, is characterized in that, its method step is:
(1) image values acquisition, be to utilize the image of the fixing object to be measured of image capture unit acquisition to obtain its pixel rgb value, and same area repeats to capture the action of image on object to be measured, to obtain the pixel rgb value of same position more than two groups;
(2) color space conversion, convert the pixel rgb value of same position more than two groups of this image to the Lab color space by one group of color conversion formula;
(3) calculate the average color length value, utilize a formula to ask the color length value of each position, its formula is as follows:
ΔE = ( L - L 0 ) 2 + ( a - a 0 ) 2 + ( b - b 0 ) 2 , Again the average color length value that represents this image is tried to achieve in the color length value Δ E addition of each position
Figure FDA00003003411200022
Its average color length value
Figure FDA00003003411200023
For acceptable scope, wherein Δ E is color length, and L represents lightness, L 0Be the reference value of lightness, a is red-green axle, a 0Be the reference value of red-green axle, the b axle is Huang-blue axle, b 0Axle is the reference value of Huang-blue axle;
(4) calculate comparation and assessment value X, computing formula is:
X = Aσ 4 × 100 % , Wherein, σ 2=σ 1 2+ σ 2 2
σ wherein 2For total variation, σ 1 2Variance, σ 2 for the envirment factor generation 2For the variance that artificial error produces, A is the yield control constant, and it is worth between 1 to 7, if the impact that the larger expression envirment factor of X value and personal error produce is larger.
3. the method for check imaging circumstances degree of stability as claimed in claim 2, is characterized in that, image capture unit is photographing unit, ccd camera or network video device.
CN2010102807248A 2010-09-14 2010-09-14 Method for eliminating chromatic aberration caused by imaging environment and inspecting stability of imaging environment Expired - Fee Related CN102397075B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010102807248A CN102397075B (en) 2010-09-14 2010-09-14 Method for eliminating chromatic aberration caused by imaging environment and inspecting stability of imaging environment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010102807248A CN102397075B (en) 2010-09-14 2010-09-14 Method for eliminating chromatic aberration caused by imaging environment and inspecting stability of imaging environment

Publications (2)

Publication Number Publication Date
CN102397075A CN102397075A (en) 2012-04-04
CN102397075B true CN102397075B (en) 2013-11-13

Family

ID=45879987

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010102807248A Expired - Fee Related CN102397075B (en) 2010-09-14 2010-09-14 Method for eliminating chromatic aberration caused by imaging environment and inspecting stability of imaging environment

Country Status (1)

Country Link
CN (1) CN102397075B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1365045A (en) * 2001-01-11 2002-08-21 金宝电子工业股份有限公司 Confirming method for automatic focusing accuracy

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7388656B2 (en) * 2005-12-30 2008-06-17 Yan Liu Method and system for color grading of gemstones
US8406514B2 (en) * 2006-07-10 2013-03-26 Nikon Corporation Image processing device and recording medium storing image processing program
US20090207624A1 (en) * 2008-02-15 2009-08-20 Acumen, Inc. Headlight assembly permitting compensation for visibility changes
US8831343B2 (en) * 2009-01-19 2014-09-09 Dolby Laboratories Licensing Corporation Image processing and displaying methods for devices that implement color appearance models

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1365045A (en) * 2001-01-11 2002-08-21 金宝电子工业股份有限公司 Confirming method for automatic focusing accuracy

Also Published As

Publication number Publication date
CN102397075A (en) 2012-04-04

Similar Documents

Publication Publication Date Title
US8319857B2 (en) Apparatus and method for correcting digital color photographs
CN103969269B (en) Method and apparatus for geometric calibration CT scanner
JP5576993B2 (en) Image measuring apparatus, image measuring method, and image measuring system
US9542758B2 (en) Image processing apparatus and light source identification method
CN104729718B (en) A kind of processing system and method for infrared imaging system NETD
CN104572538A (en) K-PLS regression model based traditional Chinese medicine tongue image color correction method
CN106296696A (en) The conforming processing method of color of image and image capture device
CN105547342A (en) Industrial lens testing device and method on the basis of liquid crystal display panel
JP6113319B2 (en) Image color distribution inspection apparatus and image color distribution inspection method
JP5611022B2 (en) Three-dimensional measuring apparatus and three-dimensional measuring method
JP5841091B2 (en) Image color distribution inspection apparatus and image color distribution inspection method
JP2022137109A (en) Image display system and image display method
JP4830068B2 (en) Two-dimensional colorimeter and spectral sensitivity correction method
US8411147B2 (en) Method for eliminating chromatic aberration caused by an imaging environment and for testing stability of the imaging environment, and chromatic aberration calibration device for use with the same
CN102397075B (en) Method for eliminating chromatic aberration caused by imaging environment and inspecting stability of imaging environment
JP2015066192A (en) Tooth colorimetric method, tooth colorimetric device, and tooth colorimetric program
JP3972647B2 (en) Diagnostic imaging apparatus, diagnostic imaging system, and diagnostic imaging method
TW201731460A (en) Hemoglobin and hematocrit analyzer and analyzing method thereof
JP2006067423A (en) Image quality quantitative evaluation method and its device
MX2021010009A (en) Film thickness gauge by near-infrared hyperspectral imaging.
CN113758683B (en) Camera system post-irradiation spectrum degradation evaluation method based on average color saturation
JP2013092372A (en) Spectral image processor, spectral image processing method and computer program
JP2004086712A (en) Image processor
CN107392955A (en) A kind of depth of field estimation device and method based on brightness
Samei et al. Method for in-field evaluation of the modulation transfer function of electronic display devices

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20131113

Termination date: 20140914

EXPY Termination of patent right or utility model