CN103630072B - The layout optimization method of video camera in two CCD camera measure system - Google Patents
The layout optimization method of video camera in two CCD camera measure system Download PDFInfo
- Publication number
- CN103630072B CN103630072B CN201310508264.3A CN201310508264A CN103630072B CN 103630072 B CN103630072 B CN 103630072B CN 201310508264 A CN201310508264 A CN 201310508264A CN 103630072 B CN103630072 B CN 103630072B
- Authority
- CN
- China
- Prior art keywords
- camera
- video camera
- point
- straight line
- centerdot
- 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
- 238000005457 optimization Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 21
- 230000002068 genetic effect Effects 0.000 claims abstract description 21
- 238000005242 forging Methods 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims abstract description 12
- 238000000605 extraction Methods 0.000 claims abstract description 9
- 239000000203 mixture Substances 0.000 claims abstract description 9
- 239000000284 extract Substances 0.000 claims abstract description 8
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 238000005259 measurement Methods 0.000 claims description 39
- 238000012360 testing method Methods 0.000 claims description 32
- 230000000007 visual effect Effects 0.000 claims description 15
- 101000857682 Homo sapiens Runt-related transcription factor 2 Proteins 0.000 claims description 12
- 102100025368 Runt-related transcription factor 2 Human genes 0.000 claims description 12
- 230000003287 optical effect Effects 0.000 claims description 8
- 230000035772 mutation Effects 0.000 claims description 7
- 230000000452 restraining effect Effects 0.000 claims description 7
- 238000001514 detection method Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 3
- 230000000694 effects Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
In two CCD camera measure system of the present invention, the layout optimization method of video camera belongs to Computer Vision Detection and field of image detection, particularly for obtain large forgings dimensional parameters two CCD camera measure system in the layout optimization method of video camera.In two CCD camera measure system, consider that the picture point that image sampling causes extracts deviation, the present invention sets up the numerical relationship model being extracted the measuring error that causes of deviation and focal length of camera, parallax range, video camera deflection angle three structural parameters by picture point, i.e. Camera composition optimized mathematical model, and utilize genetic algorithm to obtain optimum structural parameters combination.The present invention utilizes genetic algorithm to achieve the layout optimization design of video camera in two CCD camera measure system, the measuring error caused by picture point extraction error is made to reach minimum, when not demarcating video camera, namely the rational deployment that can be video camera provides effective theoretical direction and reference, there is effect of optimization good, the feature that application is strong.
Description
Technical field
The invention belongs to Computer Vision Detection and field of image detection, particularly for obtain large forgings dimensional parameters two CCD camera measure system in the layout optimization method of video camera.
Background technology
Binocular vision photogrammetry is as a kind of real-time, non-contact measurement method that measuring accuracy is high, be widely used in the numerous areas such as industrial detection, target identification, especially in the forging and stamping of measurement large forgings in real time process, in hot physical dimension, there is incomparable advantage.Many scholars have carried out large quantifier elimination around how obtaining high-precision measurement result, but these research work mainly concentrate on the raising stated accuracy of camera system and the matching precision of unique point, often have ignored the impact of structural parameters on measuring accuracy of measuring system, and structural parameters not only determine the size of apparent field, and decide the measuring accuracy of diverse location in apparent field.When utilizing two CCD camera measure system to measure, target to be measured must in apparent field.Therefore, when analytical structure parameter is on the affecting of measuring error, the research carried out under not considering apparent field's constraint all has one-sidedness.In actual measurement process, once after system calibrating, system just must keep relatively fixing, and the structural parameters of system all can not change, so before starting to demarcate, be necessary to be optimized system structure parameter.The measuring accuracy of two CCD camera measure system is improved by the installation position of reasonable Arrangement video camera.
Existing Camera composition optimization method draws general conclusion by investigating single structure parameter one by one to the impact of measuring accuracy mostly, fail each structural parameters to consider and obtain optimum structural parameters combination, or utilize first order optimization method etc. to carry out video camera to optimize distribution, but be easily absorbed in local minimum point.In fact, current most of optimization method all belongs to local optimal searching category, and its effect of optimization depends on choosing of initial value to a great extent, and genetic algorithm is a kind of heuristic random searching algorithm, has efficient global optimizing ability.In addition, the picture point seldom having focus of attention to produce by sampling extracts deviation to the impact of final measuring error.
Summary of the invention
Technical matters to be solved by this invention overcomes the deficiencies in the prior art, effective Camera composition optimization method is lacked at forging scene, and not by problems such as each structural parameters of system consider, invent a kind of layout optimization method based on video camera in the two CCD camera measure system of genetic algorithm.In two CCD camera measure system, consider that the picture point that image sampling causes extracts deviation, the present invention establishes the numerical relationship model between measuring error and focal length of camera, parallax range and video camera deflection angle three structural parameters caused by picture point extraction deviation.Distance between the Liang Tai video camera centre of perspectivity, is parallax range; The angle of two camera optical axis and Z axis, is video camera deflection angle.Under the constraints such as consideration apparent field, by the single objective programming problem that the layout optimization question variation of two CCD camera measure system is belt restraining, and use genetic algorithm to carry out global optimizing, obtain the measuring system structural parameters of one group of optimum, make the measuring error caused by picture point extraction deviation reach minimum.
The technical scheme that the present invention takes is the layout optimization method of video camera in a kind of two CCD camera measure system, in two CCD camera measure system, image sampling can cause picture point to extract deviation, sets up the measuring error and focal length of camera f, parallax range D, video camera deflection angle that are caused by picture point extraction deviation
the numerical relationship model of three structural parameters, i.e. Camera composition optimized mathematical model, and utilize genetic algorithm to obtain optimum structural parameters combination.The concrete steps that layout optimization method adopts are as follows:
Step 1: in two CCD camera measure system, sets up two Camera composition optimized mathematical models.
As shown in Figure 2, left side camera CCD1 and right camera CCD2 converges arranged in form with optical axis, the initial point O of world coordinate system OXYZ and the camera coordinate system O of left side camera
lx
ly
lz
linitial point O
loverlap.If the coordinate of a certain object point P in world coordinate system OXYZ is (X
w, Y
w, Z
w), at the camera coordinate system O of left side camera
lx
ly
lz
lwith the camera coordinate system O of right camera
rx
ry
rz
runder coordinate be respectively (X
l, Y
l, Z
l), (X
r, Y
r, Z
r), in left and right side video camera CCD1, CCD2 image planes, the image physical coordinates of picture point is respectively (x
l, y
l) and (x
r, y
r).
Object point P is at the camera coordinate system O of left and right side video camera
lx
ly
lz
l, O
rx
ry
rz
runder coordinate and the image physical coordinates of picture point there is following relation:
Object point P is at the camera coordinate system O of left and right side video camera
lx
ly
lz
l, O
rx
ry
rz
runder coordinate and its coordinate under world coordinate system OXYZ there is following relation:
Formula (3) and formula (2) are substituted into formula (1), there is following relation in the image physical coordinates obtaining the coordinate of object point P under world coordinate system OXYZ and picture point:
Simultaneous formula (1) ~ formula (5), obtains the coordinate (X of object point P under world coordinate system OXYZ
w, Y
w, Z
w) as follows:
Wherein, if the expression formula of Ψ is as follows:
Consider that sampling causes the situation of maximum picture point deviation, suppose that Pixel Dimensions is δ, then the picture point on the actual left and right camera image plane obtained and the extraction deviation existed between picture point are ideally:
ε
l=ε
r=±0.5δ(10)
Picture point physical coordinates on the left and right camera image plane that actual extracting obtains is respectively
itself and picture point physical coordinates (x ideally
l, y
l), (x
r, y
r) pass be:
The picture point physical coordinates obtained by actual extracting obtains the coordinate of object point P under world coordinate system OXYZ
as follows:
Wherein, if
expression formula as follows:
Setting up the mathematical relation that picture point extracts between the measuring error Q that causes of deviation and structural parameters is:
For obtaining the measurement space of forging ' s block dimension as shown in the hatched region in Fig. 2, it is of a size of L × W × H, L represents the size of length direction (X to), W represents the size of Width (Z-direction), H represents the size of short transverse (Y-direction), Y-axis forward be vertical paper inwards.
In XOZ plane, the angle theta of camera optical axis and the left and right visual field border line of video camera is:
Wherein, γ is effective image planes size.
The front depth of field Δ L1 of video camera is:
The rear depth of field Δ L2 of video camera is:
The depth of field Δ L of video camera is:
Wherein, F is f-number, d=|O
la|=|O
rb| be focusing from, CoC, for allowing blur circle diameter, is calculated by following formula:
CoC=a/1730(21)
Wherein, a is the catercorner length of effective image planes.
In XOZ coordinate system, build the equation of each place, border straight line, the form of straight-line equation is z=k
0x+b, wherein k
0for the slope of straight line, b is the intercept of straight line.
The equation of right camera CCD2 depth of field inner boundary l1 place straight line is:
The equation of right camera CCD2 depth of field outer boundary l2 place straight line is:
The equation of left side camera CCD1 depth of field inner boundary l3 place straight line is:
The equation of left side camera CCD1 depth of field outer boundary l4 place straight line is:
The equation of left margin l5 place, right camera CCD2 visual field straight line is:
The equation of left margin l6 place, left side camera CCD1 visual field straight line is:
The equation of right margin l7 place, right camera CCD2 visual field straight line is:
The equation of right margin l8 place, left side camera CCD1 visual field straight line is:
On measurement space length direction, the equation of left margin l9 (l9') place straight line is:
x=D/2-L/2(30)
On measurement space length direction, the equation of right margin l10 place straight line is:
x=D/2+L/2(31)
The statement of measurement space Width size W is divided into two kinds of situations: | CU|=W, | C'V'|=W.Then there is W=min{|CU|, | C'V'|}.Wherein,
According to the range of size of actual measurement demand determination measurement space, the size L>=L of length direction (X to)
0, the size W>=W of Width (Z-direction)
0, the size H>=H of short transverse (Y-direction)
0, should L=L be met
0time corresponding measurement space width
choose and be positioned at apparent field center L
0× W
0× H
0in measurement space, an equally distributed p position is as test point, and the layout scenarios of test point is as follows: according to being divided into q the test plane be equally spaced along Z-direction direction from the close-by examples to those far off, test interplanar spacing is W
0/ (q-1), each test plane is L
0× H
0region, each region is chosen r is capable, s arranges uniform test point.
Video camera adopts zoom lens, and the range of adjustment of focal distance f is [f
1, f
2]; The scope of parallax range D is [D
1, D
2]; Video camera depth of field Δ L>=W
0; Distance between reference measure system and forging arranges suitable focusing from d; Video camera deflection angle
scope be [0, θ];
The layout optimization question variation of this two CCD camera measure system is the single objective programming problem of following belt restraining:
Wherein, Σ Q is the summation of p test point measuring error.Γ is the mean value of each test point measuring error, and it is objective function to be optimized.
for the reference position along Z-direction two camera field of view laps, its value is less than focusing from d.
Step 2: use genetic algorithm to find optimum solution, namely optimum structural parameters combination.
For the single objective programming problem of the belt restraining obtained in step 1, genetic algorithm is adopted to carry out global optimizing.Detailed process is as follows:
(1) algorithm routine is write, and comprises and writes the objective function of optimization problem described in step 1 and the program file of non-linear constrain.
(2) optimized variable is focal distance f, parallax range D and video camera deflection angle
the span of each variable is respectively [f
1, f
2], [D
1, D
2] and
(3) Population Size is set, crossover probability, mutation probability and algorithm end condition.
Population Size elects PopSize as, and crossover probability elects p as
c, mutation probability elects p as
m.The stop criterion of algorithm is: iterations reaches maximum iteration time MaxGen; Or when algorithm is in the algebraically ConGen that stagnation algebraically specifies, the weighted mean change of fitness function is less than function franchise FunTol.
(4) run the program of genetic algorithm, draw the optimum solution of optimization problem, namely optimum structural parameters combination, makes the measuring error in whole measurement space reach minimum.
The invention has the beneficial effects as follows the layout optimization design utilizing genetic algorithm to achieve video camera in two CCD camera measure system, the measuring error caused by picture point extraction error is made to reach minimum, when not demarcating video camera, namely the rational deployment that can be video camera provides effective theoretical direction and reference, there is effect of optimization good, the features such as application is strong.
Accompanying drawing explanation
Fig. 1 is the schematic diagram that analog image is sampled as digital picture.A () figure is the analog image of unique point, (b) figure is corresponding digital picture after unique point sampling, wherein: 1-unique point, and 2-boost line, 3-pixel, 4-unique point place pixel.
Fig. 2 is the structural representation of two CCD camera measure system, and Fig. 3 is the distribution schematic diagram of test point.Wherein: CCD1-left side camera, CCD2-right camera, O
lx
ly
lz
lthe camera coordinate system of-left side camera, O
l-left side camera the centre of perspectivity, O
lz
l-left side camera optical axis, | O
lthe focusing of A|-left side camera is from, O
rx
ry
rz
rthe camera coordinate system of-right camera, O
r-right camera the centre of perspectivity, O
rz
r-right camera optical axis, | O
rthe focusing of B|-right camera from, OXYZ-world coordinate system, O-world coordinate system initial point, D-parallax range, f-focal length of camera,
-video camera deflection angle, the angle of θ-camera optical axis and camera field of view boundary line, L-measurement space length, W-measurement space width, Z
0-along the reference position of Z-direction two camera field of view laps, the Δ L-video camera depth of field, l1-right camera depth of field inner boundary, l2-right camera depth of field outer boundary, l3-left side camera depth of field inner boundary, l4-left side camera depth of field outer boundary, l5-right camera visual field left margin, l6-left side camera visual field left margin, l7-right camera visual field right margin, l8-left side camera visual field right margin, left margin on l9-measurement space length direction, the intersection point of C-l9 and l2, the intersection point of U-l9 and l3, the intersection point of V-l9 and l5, right margin on l10-measurement space length direction, left margin on measurement space length direction in the another kind of situation of l9'-, the intersection point of C'-l9' and l2, the intersection point of U'-l9' and l3, the intersection point of V'-l9' and l5, P
k, (i, j)test point in-kth test a plane on the i-th row jth row.
Embodiment
The specific embodiment of the present invention is further described below in conjunction with accompanying drawing and technical scheme.
As shown in Figure 1, after only having the sampling of the analog image (a) of a unique point 1, obtain the corresponding digital picture (b) of unique point 1, continuous print image planes by discrete be many square region, i.e. pixel 3.For the characteristic point positioning method of Pixel-level precision, the positional information of unique point 1 place pixel 4 can only being obtained, if think that the position of this pixel is the position of unique point 1, so extracting deviation by there is the picture point caused by image sampling.The system structure parameter that the present invention considers is the deflection angle of focal length of camera f, parallax range D and video camera
the present invention establishes the numerical relationship model of measuring error and the said structure parameter caused by picture point extraction deviation, i.e. Camera composition optimized mathematical model, under the constraint of consideration apparent field, use genetic algorithm to carry out global optimizing, obtain the measuring system structural parameters of one group of optimum.Concrete steps are as follows:
Step 1: in two CCD camera measure system, sets up two Camera composition optimized mathematical models.
The imaging device that the present invention selects is Princeton MegaPlusIIES4020 type black-white CCD video camera, and its resolution is 2048 × 2048, and Pixel Dimensions is δ=7.4 μm; Effective image planes are of a size of 15.2mm × 15.2mm, i.e. γ=15.2mm, catercorner length
for adapting to different visual field demands, two video cameras are all furnished with TamronDi-IILD zoom lens, and the range of adjustment of its focal distance f is [f
1, f
2]=[18,250], unit is mm.The range of adjustment of the f-number F of this camera lens is [3.5,6.3], and for making logical light quantity large as far as possible, to obtain image clearly, lens aperture value F selects 3.5.Suppose that the parameters of two video cameras is identical.
The derivation of aforesaid formula (1) ~ formula (9) is all be based upon on the basis of the perspective projection relation met ideally, and namely object point distributes at Spatial continual.But, the object point of Spatial continual is sampled as picture point discrete in digital camera image planes, cause and there is picture point extraction deviation, consider the situation causing maximum deviation in sampling process, then tried to achieve to exist between picture point in the left and right picture plane of actual acquisition and picture point ideally by formula (10) and extract deviation:
ε
l=ε
r=±3.7μm
Calculating permission blur circle diameter CoC by formula (21) is:
When utilizing binocular vision system to measure, target to be measured must in apparent field, and apparent field refers to the overlapped fov be in two video camera field depths, as shown in the shadow region in Fig. 2.For obtaining the measurement space of forging ' s block dimension as shown in the hatched region in Fig. 2.According to practical measurement requirement, require that measurement space need meet: L
0=3000mm, W
0=3000mm, H
0=3000mm, i.e. the size L>=3000mm of length direction (X to), the size W>=1000mm of Width (Z-direction), the size H>=3000mm of short transverse (Y-direction).
The measuring error of two CCD camera measure system is uneven in the distribution of measurement space, for without loss of generality, choose and be positioned at apparent field center 3000mm × 1000mm × 3000mm measurement space (the shade rectangular parallelepiped region shown in Fig. 3) equally distributed 180 location points as test point, be i.e. p=180.Test point is numbered according to following rule: according to being divided into 5 tests plane, i.e. q=5 being equally spaced along Z-direction direction from the close-by examples to those far off.Interplanar spacing is 1000/4=250mm, and each test plane is the region of 3000mm × 3000mm, each region is chosen 6 row × 6 and arranges uniform test point, i.e. r=6, s=6.Test point coding rule in each test plane is encoded consistent with the pixel coordinate of image, and the upper left corner is coordinate origin, and (i, j) represents the test point on the i-th row jth row, and the sequence number of row increases progressively from top to bottom, and the sequence number of row increases progressively from left to right, uses P
k, (i, j)represent the test point on the i-th row jth row in a kth test plane.Therefore, the value of all test points each coordinate components in the ideal coordinates under world coordinate system can be obtained, X
wvalue have: { D/2-1500, D/2-900, D/2-300, D/2+300, D/2+900, D/2+1500}, Y
wvalue have :-1500 ,-900 ,-300,300,900,1500}, Z
wvalue according to intersection point C (C'), U (U') and the focusing of V (V') the relative camera position from d place,
evenly choose in scope 5 be worth and be spaced apart 250mm or
evenly choose 5 in scope be worth and be spaced apart 250mm.Wherein, if the expression formula of d' is:
The Z-direction coordinate of intersection point C (C'), U (U') and V (V') is solved under XOZ coordinate system:
The Z-direction coordinate that C (C') puts is:
The Z-direction coordinate that U (U') puts is:
The Z-direction coordinate that V (V') puts is:
According to forging on-site actual situations, consideration equipment install and avoid equipment because of some constraint conditions obtaining other apart from factors such as hot large forgings are nearer and overheated as follows: the statement of the size W of measurement space Width is divided into two kinds of situations, a kind of situation is that the distance between C and U equals W, namely | CU|=W, another kind of situation be on measurement space length direction left margin l9' when apparent field's right side edge, distance between C' and V' equals W, namely | C'V'|=W, then there is W=min{|CU|, | C'V'|}.During measurement space lengthwise dimension L=3000mm, then corresponding measurement space Width size W
l=3000=min{|CU|
l=3000, | C'V'|
l=3000and W
l=3000>=W
0, therefore have | CU|
l=3000>=1000 and | C'V'|
l=3000>=1000.Consider the compactedness requirement of two CCD camera measure system, in conjunction with the restriction of video camera own dimensions and erecting bed size, parallax range D span is [D
1, D
2]=[100,2000], unit is mm.Because measurement space is arranged in the overlapped fov of two video camera field depths, therefore, video camera depth of field Δ L must be greater than the width W of measurement space
0, namely have Δ L>=1000mm.For avoiding high temperature to cause damage to measuring equipment, under measuring system distance forge press, anvil center is at least d
1=8000mm is far away, with reference to d
1choose video camera focusing from d=8000mm.For making full use of apparent field, video camera deflection angle
scope is [0, θ].By formula (17), in conjunction with focal range [f
1, f
2] and effectively image planes size γ calculate the span of θ for [1.74 °, 22.89 °].
In sum, the layout optimization question variation of this two CCD camera measure system is the single objective programming problem of following belt restraining:
Wherein, Q
k, (i, j)for the measuring error of the upper test point of the i-th row jth row in a kth test plane, Z
0for along the overlapping reference position of Z-direction two camera field of view, Z should be had
0< 8000mm.
Step 2: use genetic algorithm to find optimum solution, namely optimum structural parameters combination.
For the single objective programming problem of the belt restraining obtained in step 1, genetic algorithm is adopted to carry out global optimizing.From initial population, use genetic operator, namely selection opertor, crossover operator and mutation operator produce population of future generation, so make population evolve towards the direction of optimization solution, until meet the end condition of setting.Detailed process is as follows:
(1) algorithm routine is write, and comprises and writes the objective function of optimization problem described in step 1 and the program file of non-linear constrain.
(2) optimized variable is focal distance f, parallax range D and video camera deflection angle
each variable-value scope is respectively [18,250], [100,2000] and [0,22.89], and unit is respectively millimeter, millimeter and degree.
(3) Population Size is set, crossover probability, mutation probability and algorithm end condition
Algorithm parameter is set as PopSize=500, p
c=0.78, p
m=0.008, MaxGen=2000, ConGen=500, FunTol=10
-12, that is: Population Size elects 500 as, and crossover probability elects 0.78 as, and mutation probability elects 0.008 as.The stop criterion of algorithm is: iterations reaches maximum iteration time 2000; Or when algorithm is in the algebraically 500 that stagnation algebraically specifies, the weighted mean change of fitness function is less than function franchise 10
-12.
(4) call genetic algorithm master routine, draw the optimum solution of optimization problem, namely optimum structural parameters combination, makes the measuring error in whole measurement space reach minimum.
In step 1, the optimum solution of optimization problem is focal length of camera f=42.552mm, parallax range D=2000mm, video camera deflection angle
target function value converges to minimum value 5.834mm.Above-mentioned experiment is carried out on the computing machine of 1.99GB internal memory, 2.09GHz processor, and on average consuming time is 2.8h.
In two CCD camera measure system, image sampling can cause picture point to extract deviation, the present invention extracts mathematical relation between the measuring error that causes of deviation and structural parameters by setting up picture point, consider that measurement space, equipment are installed and avoids equipment because of apart from factors such as large forgings radiation are nearer and overheated, for video camera and the fixed two CCD camera measure system of camera lens, adopt genetic algorithm to carry out the layout optimization design of video camera, rely on the efficient global optimizing ability of genetic algorithm can obtain the structural parameters combination of global optimum.
Claims (1)
1. the layout optimization method of video camera in a two CCD camera measure system, in two CCD camera measure system, consider that the picture point that image sampling causes extracts deviation, set up the numerical relationship model being extracted the measuring error that causes of deviation and focal length of camera, parallax range, video camera deflection angle three structural parameters by picture point, i.e. Camera composition optimized mathematical model, and utilize genetic algorithm to obtain optimum structural parameters combination, it is characterized in that, the concrete steps that layout optimization method adopts are as follows:
Step 1: in two CCD camera measure system, sets up two Camera composition optimized mathematical models;
Left side camera CCD1 and right camera CCD2 converges arranged in form with optical axis, the initial point O of world coordinate system OXYZ and the camera coordinate system O of left side camera
lx
ly
lz
linitial point O
loverlap; If the coordinate of a certain object point P in world coordinate system OXYZ is (X
w, Y
w, Z
w), at the camera coordinate system O of left side camera
lx
ly
lz
lwith the camera coordinate system O of right camera
rx
ry
rz
runder coordinate be respectively (X
l, Y
l, Z
l), (X
r, Y
r, Z
r), in left and right side video camera CCD1, CCD2 image planes, the image physical coordinates of picture point is respectively (x
l, y
l) and (x
r, y
r);
Object point P is at the camera coordinate system O of left and right side video camera
lx
ly
lz
l, O
rx
ry
rz
runder coordinate and the image physical coordinates of its picture point there is following relation:
Object point P is at the camera coordinate system O of left and right side video camera
lx
ly
lz
l, O
rx
ry
rz
runder coordinate and its coordinate under world coordinate system OXYZ there is following relation:
Formula (3) and formula (2) are substituted into formula (1), there is following relation in the image physical coordinates obtaining the coordinate of object point P under world coordinate system OXYZ and picture point:
Simultaneous formula (1) ~ formula (5), obtains the coordinate (X of object point P under world coordinate system OXYZ
w, Y
w, Z
w) as follows:
Wherein, if the expression formula of Ψ is as follows:
Consider that sampling causes the situation of maximum picture point deviation, suppose that Pixel Dimensions is δ, then the picture point on the actual left and right camera image plane obtained and the extraction deviation existed between picture point are ideally:
ε
l=ε
r=±0.5δ(10)
Picture point physical coordinates on the left and right camera image plane that actual extracting obtains is respectively
itself and picture point physical coordinates (x ideally
l, y
l), (x
r, y
r) pass be:
The picture point physical coordinates obtained by actual extracting obtains the coordinate of object point P under world coordinate system OXYZ
as follows:
Wherein, if
expression formula as follows:
Setting up the mathematical relation that picture point extracts between the measuring error Q that causes of deviation and structural parameters is:
Measurement space for obtaining forging ' s block dimension is of a size of L × W × H, L represents the size of length direction (X to), W represents the size of Width (Z-direction), and H represents the size of short transverse (Y-direction), Y-axis forward be vertical paper inwards;
In XOZ plane, the angle theta of camera optical axis and the left and right visual field border line of video camera is:
Wherein, γ is effective image planes size;
The front depth of field Δ L1 of video camera is:
The rear depth of field Δ L2 of video camera is:
The depth of field Δ L of video camera is:
Wherein, F is f-number, d=|O
la|=|O
rb| be focusing from, CoC, for allowing blur circle diameter, is calculated by following formula:
CoC=a/1730(21)
Wherein, a is the catercorner length of effective image planes;
In XOZ coordinate system, build the equation of each place, border straight line, the form of straight-line equation is z=k
0x+b, wherein k
0for the slope of straight line, b is the intercept of straight line;
The equation of right camera CCD2 depth of field inner boundary l1 place straight line is:
The equation of right camera CCD2 depth of field outer boundary l2 place straight line is:
The equation of left side camera CCD1 depth of field inner boundary l3 place straight line is:
The equation of left side camera CCD1 depth of field outer boundary l4 place straight line is:
The equation of left margin l5 place, right camera CCD2 visual field straight line is:
The equation of left margin l6 place, left side camera CCD1 visual field straight line is:
The equation of right margin l7 place, right camera CCD2 visual field straight line is:
The equation of right margin l8 place, left side camera CCD1 visual field straight line is:
On measurement space length direction, the equation of left margin l9 (l9') place straight line is:
x=D/2-L/2(30)
On measurement space length direction, the equation of right margin l10 place straight line is:
x=D/2+L/2(31)
The statement of measurement space Width size W is divided into two kinds of situations: | CU|=W or | C'V'|=W, then have W=min{|CU|, | C'V'|}; Wherein,
According to the range of size of actual measurement demand determination measurement space, the size L>=L of length direction (X to)
0, the size W>=W of Width (Z-direction)
0, the size H>=H of short transverse (Y-direction)
0, should L=L be met
0time corresponding measurement space width
; Choose and be positioned at apparent field center L
0× W
0× H
0in measurement space, an equally distributed p position is as test point, and the layout scenarios of test point is as follows: according to being divided into q the test plane be equally spaced along Z-direction direction from the close-by examples to those far off, test interplanar spacing is W
0/ (q-1), each test plane is L
0× H
0region, each region is chosen r is capable, s arranges uniform test point;
Video camera adopts zoom lens, and the range of adjustment of focal distance f is [f
1, f
2]; The scope of parallax range D is [D
1, D
2]; Video camera depth of field Δ L>=W
0; Distance between reference measure system and forging arranges suitable focusing from d; Video camera deflection angle
scope be [0, θ];
The layout optimization question variation of this two CCD camera measure system is the single objective programming problem of following belt restraining:
Wherein, ∑ Q is the summation of p test point measuring error; Γ is the mean value of each test point measuring error, and it is objective function to be optimized;
for the reference position along Z-direction two camera field of view laps, its value is less than focusing from d;
Step 2: use genetic algorithm to find optimum solution, namely optimum structural parameters combination;
For the single objective programming problem of the belt restraining obtained in step 1, genetic algorithm is adopted to carry out global optimizing; Detailed process is as follows:
(1) algorithm routine is write, and comprises and writes the objective function of optimization problem described in step 1 and the program file of non-linear constrain;
(2) optimized variable is focal distance f, parallax range D and video camera deflection angle
the span of each variable is respectively [f
1, f
2], [D
1, D
2] and
(3) Population Size is set, crossover probability, mutation probability and algorithm end condition;
Population Size elects PopSize as, and crossover probability elects p as
c, mutation probability elects p as
m; The stop criterion of algorithm is: iterations reaches maximum iteration time MaxGen; Or when algorithm is in the algebraically ConGen that stagnation algebraically specifies, the weighted mean change of fitness function is less than function franchise FunTol;
(4) run the program of genetic algorithm, draw the optimum solution of optimization problem, namely optimum structural parameters combination, makes the measuring error in whole measurement space reach minimum.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310508264.3A CN103630072B (en) | 2013-10-25 | 2013-10-25 | The layout optimization method of video camera in two CCD camera measure system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310508264.3A CN103630072B (en) | 2013-10-25 | 2013-10-25 | The layout optimization method of video camera in two CCD camera measure system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103630072A CN103630072A (en) | 2014-03-12 |
CN103630072B true CN103630072B (en) | 2016-01-13 |
Family
ID=50211363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310508264.3A Active CN103630072B (en) | 2013-10-25 | 2013-10-25 | The layout optimization method of video camera in two CCD camera measure system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103630072B (en) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105430266A (en) * | 2015-11-30 | 2016-03-23 | 努比亚技术有限公司 | Image processing method based on multi-scale transform and terminal |
CN107066747A (en) * | 2017-04-24 | 2017-08-18 | 哈尔滨理工大学 | A kind of vision measurement network organizing planing method |
CN107443369A (en) * | 2017-06-25 | 2017-12-08 | 重庆市计量质量检测研究院 | A kind of robotic arm of the inverse identification of view-based access control model measurement model is without demarcation method of servo-controlling |
CN107483820A (en) * | 2017-08-23 | 2017-12-15 | 福建星瑞格软件有限公司 | A kind of quick focus method and device of the focusing of contrast formula |
CN108562265B (en) * | 2018-03-12 | 2021-01-12 | 内蒙古大学 | Measuring range estimation method of binocular stereo vision distance measuring device |
CN109003307B (en) * | 2018-06-11 | 2021-10-22 | 西北工业大学 | Underwater binocular vision measurement-based fishing mesh size design method |
CN109102548B (en) * | 2018-08-23 | 2020-12-11 | 武汉中观自动化科技有限公司 | Method and system for identifying tracking range |
CN109084688B (en) * | 2018-09-20 | 2020-09-29 | 杭州电子科技大学 | Binocular vision distance measurement method based on variable-focus camera |
CN110342252B (en) * | 2019-07-01 | 2024-06-04 | 河南启迪睿视智能科技有限公司 | Automatic article grabbing method and automatic grabbing device |
CN110851978B (en) * | 2019-11-08 | 2024-03-19 | 江苏科技大学 | Camera position optimization method based on visibility |
CN111460735B (en) * | 2020-04-08 | 2023-09-01 | 深圳市瑞立视多媒体科技有限公司 | Camera layout function optimization method based on genetic inheritance and related equipment |
CN112254672B (en) * | 2020-10-15 | 2022-02-15 | 天目爱视(北京)科技有限公司 | Height-adjustable's intelligent 3D information acquisition equipment |
CN114329854B (en) * | 2020-11-06 | 2023-05-12 | 北京航空航天大学 | Two-dimensional space vision sensor layout optimization method based on multi-target constraint |
CN112556639B (en) * | 2020-11-06 | 2023-05-09 | 广州艾目易科技有限公司 | Device and method for testing actual effective visual field range of binocular vision system |
CN114329855B (en) * | 2020-11-06 | 2023-05-12 | 北京航空航天大学 | Sensor layout optimization and rapid deployment method of wireless vision sensing network |
CN113240829B (en) * | 2021-02-24 | 2022-09-23 | 南京工程学院 | Intelligent gate passing detection method based on machine vision |
CN113532473B (en) * | 2021-06-17 | 2024-04-19 | 浙江工业大学 | Image pickup measurement error suppression method by arranging near-field stationary points |
CN115529437A (en) * | 2021-06-25 | 2022-12-27 | 青岛海信智慧生活科技股份有限公司 | Method, device, equipment and medium for determining monitoring equipment arrangement information |
CN114724323B (en) * | 2022-06-09 | 2022-09-02 | 北京科技大学 | Point distribution method of portable intelligent electronic fence for fire scene protection |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251925A (en) * | 2007-12-29 | 2008-08-27 | 中国科学院计算技术研究所 | Virtual video camera planning and distributing method and system |
CN103247053A (en) * | 2013-05-16 | 2013-08-14 | 大连理工大学 | Accurate part positioning method based on binocular microscopy stereo vision |
CN104349134A (en) * | 2013-08-07 | 2015-02-11 | 安讯士有限公司 | Method and system for selecting position and orientation for a monitoring camera |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU3249600A (en) * | 1999-03-18 | 2000-10-04 | Showbites, Inc. | Method for optimization of video coverage |
-
2013
- 2013-10-25 CN CN201310508264.3A patent/CN103630072B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101251925A (en) * | 2007-12-29 | 2008-08-27 | 中国科学院计算技术研究所 | Virtual video camera planning and distributing method and system |
CN103247053A (en) * | 2013-05-16 | 2013-08-14 | 大连理工大学 | Accurate part positioning method based on binocular microscopy stereo vision |
CN104349134A (en) * | 2013-08-07 | 2015-02-11 | 安讯士有限公司 | Method and system for selecting position and orientation for a monitoring camera |
Non-Patent Citations (1)
Title |
---|
立体视觉测量中的摄像机优化布局;陈杰春 等;《机床与液压》;20120831;第40卷(第15期);34-37 * |
Also Published As
Publication number | Publication date |
---|---|
CN103630072A (en) | 2014-03-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103630072B (en) | The layout optimization method of video camera in two CCD camera measure system | |
Chen et al. | High-accuracy multi-camera reconstruction enhanced by adaptive point cloud correction algorithm | |
CN103714535B (en) | Binocular vision measurement system camera parameter online adjustment method | |
CN101814185B (en) | Line structured light vision sensor calibration method for micro-size measurement | |
CN102063718B (en) | Field calibration and precision measurement method for spot laser measuring system | |
CN103292695A (en) | Monocular stereoscopic vision measuring method | |
CN104331896A (en) | System calibration method based on depth information | |
CN104484887B (en) | External parameters calibration method when video camera is used in combination with scanning laser range finder | |
CN104408689A (en) | Holographic-image-based streetscape image fragment optimization method | |
Zhu et al. | Camera calibration method based on optimal polarization angle | |
CN103994723A (en) | Macro and micro composite grating ruler measurement system based on vertical and horizontal transformation amplification segmentation | |
CN103632364A (en) | Camera spatial position relation calibration device in multi-camera photographing measurement system | |
CN107816942A (en) | A kind of planar dimension measurement method based on cross structure light vision system | |
CN103308000B (en) | Based on the curve object measuring method of binocular vision | |
CN110879080A (en) | High-precision intelligent measuring instrument and measuring method for high-temperature forge piece | |
CN105066962A (en) | Multiresolution large visual field angle high precision photogrammetry apparatus | |
CN107092905B (en) | Method for positioning instrument to be identified of power inspection robot | |
WO2015098222A1 (en) | Information processing device, information processing method, and program | |
CN103824277A (en) | Substation three-dimensional live-action monitoring stationing method based on nonlinear parameter optimization calibration | |
CN105469412A (en) | Calibration method of assembly error of PTZ camera | |
CN108088381A (en) | A kind of contactless minim gap method for measuring width based on image procossing | |
CN103913149A (en) | Binocular range finding system based on STM 32 single chip microcomputer and range finding method thereof | |
CN111681283A (en) | Monocular stereoscopic vision-based relative pose calculation method applied to wireless charging alignment | |
CN104048603A (en) | System and method for measuring array laser reticle and area array CCD high-resolution large visual field | |
CN102663727A (en) | Method for calibrating parameters by dividing regions in a camera based on CMM moving target |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |