Open access peer-reviewed chapter

Probabilistic Analysis of Projected Features in Binocular Stereo

By Lorenzo J. Tardón, Isabel Barbancho and Carlos Alberola-López

Submitted: December 5th 2011Reviewed: April 18th 2012Published: July 11th 2012

DOI: 10.5772/46027

Downloaded: 1141

1. Introduction

Some geometrical relationships between projected primitives in binocular stereo systems will be analysed in the next sections with the aim of providing a characterization from a probabilistic point of view. To this end, we will consider the parallel stereo system model and the well known pinhole camera model [1].

The characterizations that will be derived will be readily usable as valuable sources of information to solve the correspondence problem in stereo systems [2] and their nature will be that of a priori information sources in Bayesian models.

To begin with, we will introduce the stereo system model that will be used for the analysis together with the notation that will be employed and the parameters that will be necessary for the calculations. Afterwards, we will use this model to derive the joint probability density function (pdf) of the orientation of the projections on the image planes of arbitrary small edges. In this case, we will find a cumbersome expression so, then, we will focus on the derivation of a tractable pdf of a convenient function of the orientation of the projections.

Later, we will turn our attention to the so called disparity gradient, which defines important relationships between projections in stereo systems. We will find three different usable pdfs of the disparity gradient that can be used to solve the correspondence problem in parallel stereo systems. Finally, a brief summary will be drawn.

2. Geometric relationships in the parallel stereo system model

In order to perform our analysis, we consider a common model for stereo image acquisition systems. The two cameras of the stereo system are considered to be identical. These cameras are modelled using the well known pinhole camera model with focal length f, parallel optical axes and image planes defined on the same geometric plane [3], [1]. This description defines the so called parallel stereo system model. An illustration of the geometry and the projection process with this model is represented in Fig. .

For simplicity, the centre of the real world coordinate system is considered to be equidistant to the optical centres of the two cameras of the system (Cland Cr). The optical centers of the cameras are separated a distance b: the baseline. As shown in Fig. , the Xaxis is parallel to the linebase band the Zaxis is perpendicular to the image planes.

Figure 1.

Parallel stereo system model.

In Fig. , Aand Brepresent the edges of a straight segment AB¯of length δ. Ais located at (X,Y,Z)in the world coordinate system. The segment has an arbitrary orientation described by the angles αand βdefined with respect to the XYand XZplanes, respectively.

The edge points and the segment are projected onto the left and right image planes of our parallel stereo system. Thus, we find the projected points Aland Blon the left image and the projected segment δlon the same image. Also, the angle between δland the horizontal on the left image is denoted θl. Similarly, on the right image plane we find Ar, Br, δrand the angle θr.

Recall that the optical axes of the two cameras are parallel in our stereo model. Also, we consider that equally numbered horizontal lines on the two image planes comply with the epipolar constraint [4].

The segments on the image planes that correspond to the projection of the same segment in the real world are partially characterized and related by their respective orientations on the left and right images. This orientation can be analysed to be used to solve the correspondence problem in stereo systems.

Using the model selected, we will focus in the next sections on the orientation of the projection of small straight edges (δland δr). Then, we will also consider a well known feature: the disparity gradient [5], and we will show how to develop probability characterizations of this feature under different conditions [6].

3. Joint probability density function of the orientation of projected edges

Making use of the geometrical relationships established in the previous section and in Fig. , we will derive a relationship between the location and orientation of the edgel δ[7] in the real world, and the orientations of its projections described by the angles θland θr(Fig. ) in the corresponding image planes. Then, under appropriate hypotheses, we will find the description of the joint probabilistic behaviour of the projected angles.

Consider the definitions and the geometry shown in Fig. where the length of the segment δis arbitrarily small. We can write the location of the projected points in the left and right images using their coordinates on the corresponding image planes [1]. Let Bl=(Blx,Bly), then, using the geometry involved and using Band its projections as starting reference, we can write Al=(Alx,Aly)=(Blx+δlcosθl,Bly+δlsinθl).

Now, let's look at the right (r) image. Under the hypotheses described previously, and using the length of the projected edgel on the right image, δr, making use of the fact that the ycoordinates must be the same in the two images, it is simple to observe that δlsinθl=δrsinθrand, so, δr=δlsinθlsinθr.

After these observations, the coordinates of the projections of Aand Bcan be written as follows:

Ar=(Arx,Ary)=(Brx+δlsinθlsinθrcosθr,Bry+δlsinθl)Br=(Brx,Bry)uid2

But our objective must be to find the relation between the projections and the orientation of the edgel in the real world, such orientation is described by the angles αand βin Fig. . Working in this direction, the following relations can be observed:

α=arctanBz-AzBx-Axβ=arctanBy-Ay(Bx-Ax)2+(Bz-Az)2=arcsinBy-Ay(Bx-Ax)2+(By-Ay)2+(Bz-Az)2uid3

On the other hand, using the projection equations of the pinhole camera model [1], the following relations can be found:

X=-xlbxr-xl-b2Y=-ylbxr-xlZ=fbxr-xluid4

where (X,Y,Z)correspond to the coordinates of a generic point in the real world and (xl,yl), (xr,yr)correspond to its projections on the left and right images, respectively.

Now, using eqs. () to () together with eqs. () and (), it is possible to find the expressions of the following terms involved in the calculation of the projected angles:

Bx-AxBy-AyBz-Azuid5

Then, using these expressions in eq. () and writing all the terms as functions of the real world coordinates of A, the coordinates of Br, the camera parameters fand band the orientation of the projections of the edgel (θland θr), we find the equations that lead us from (α,β)to (θl,θr):

α=arctanZsin(θr-θl)Xsin(θr-θl)-b2sin(θr+θl)β=arctanbsinθlsinθr-Ysin(θr-θl)[Xsin(θr-θl)-bsin(θr+θl)]2+[Zsin(θr-θl)]2uid6

After these operations, we are ready to derive the joint pdf of the orientation of the projections of the segment: fθl,θr(θl,θr). To this end, only the pdf of (α,β)is required at this stage.

Since there is no reason to think differently, we will assume that these two parameters are independent uniform random variables (rv's) ranging from 0 to π[8]. Under these hypotheses, it is evident that the joint pdf of (α,β)is fαβ(α,β)=1π2. So, in order to derive the desired expression, we only need to calculate the modulus of the Jacobian of the transformation [9]:

|Jd|=αθlαθrβθlβθruid7

Thus, we must find the partial derivatives of αand βwith respect to θly θr. These are not simple expressions because of the functions involved. As an example, observe the result obtained for the last element of Jd:

βθr=[bsinθlcosθr-Ycos(θr-θl)]{[Xsin(θr-θl)-bsin(θr+θl)]2+Z2sin2(θr-θl)}-...{[Xsin(θr-θl)-bsin(θr+θl)]2+......[bsinθlsinθr-Ysin(θr-θl)]{[Xsin(θr-θl)-bsin(θr+θl)]......Z2sin2(θr-θl)+[bsinθlsinθr-Ysin(θr-θl)]2}......[Xcos(θr-θl)-bcos(θr+θl)]+Z2sin(θr-θl)cos(θr-θl)}[Xsin(θr-θl)-bsin(θr+θl)]2+Z2sin2(θr-θl)uid8

Since analytical expressions for all the required terms can be found by direct calculations, it is possible to obtain the desired pdf operating in the usual way [9]:

fθl,θr(θl,θr)=1π2|Jd|uid9

Unfortunately, this expression far from being simple because of the complexity of the terms involved. This fact should encourage us to search for a more usable expression capable of statistically describing a certain relation between the orientation of the projected segments. In the next section, we find such expression by using a function of cot θland cot θr.

4. Probability density function of the difference of the cot of the orientation of projected segments

A tractable expression to relate the orientation of projected segments can be found by defining a suitable function of the projected angles shown in Fig. . Let fK(k), with ka function of {θl,θr}denote such function.

More specifically, the pdf of the modulus of the difference of the cot of the projected angles in the selected binocular stereo system will be derived.

Taking into account the scene depicted in Fig. , let AB¯define, again, a straight segment with arbitrary length δ. The orientation of this segment is described by the angles αy βas shown in the figure.

Now, the location of the edges of the segment in the real world coordinate system will be written as follows:

A:(Ax,Ay,Az)=(X,Y,Z)uid10
B:(Bx,By,Bz)=(X+δcosβcosα,Y-δsinβ,Z-δcosβsinα)uid11

And taking into account the geometry selected, the coordinates of the projections of the edges of the segment can be written as:

Arx=-fAz(Ax-b2)Alx=-fAz(Ax+b2)Ary=-fAzAyAly=-fAzAyuid12
Brx=-fBz(Bx-b2)Blx=-fBz(Bx+b2)Bry=-fBzByBly=-fBzByuid13

Now, let

k=|cot(θl)-cot(θr)|uid14

Substituting the cot functions by the corresponding expressions in terms of the projections of the edges of the segment, using the projection equations () to (), multiplying by AzBz, substituting Bias a function of the coordinates of Aand dividing by cosβ, the following expression is found:

k=-bsinαZtanβ-Ysinαuid15

This expression will be used to derive the pdf of k.

To begin with, the joint pdf of kand αwill be derived. To this end, the following transformation equations will be used:

k=-bsinαZtanβ-Ysinαα=αuid16

The modulus of the Jacobian of the transformation can be easily determined:

|J|=kαkβαααβ=bsinαZsec2β(Ztanβ-Ysinα)2uid17

With all this, the joint pdf of kand αcan be readily obtained [10], [9]:

fk,α(k,α)=rf(α(kr,αr),β(kr,αr))1|Jr|uid18

where rrepresents the set of roots of the transformation of (α,β)as a function of (k,α). Two different solutions can be found for this transformation because of the modulus operation in equation ():

β=arctansinαkY+bkZ,withk=bsinαZtanβ-Ysinalphaβ=arctansinαkY-bkZ,withk=-bsinαZtanβ-Ysinalphaα=αuid19

Assuming, that the orientation angles αand βbehave as uniform random variables [8] with range (0,π)and assuming independence, it is clear that f(α,β)=1π2[9]. Then, equation () can be written, after substitution of the terms involved as:

fk,α(k,α)=1π2(Ztanβ-Ysinα)2bsinαZsec2ββ=arctansinαkY+bkZ+1π2(Ztanβ-Ysinα)2bsinαZsec2ββ=arctansinαkY-bkZuid20

Now, αand βcan be expressed in terms of αand k, making use of the following identity: sec[arctana]=1+a2. Thus, the following expression is found after some simplifications:

fk,α(k,α)=1π2bsinαk2Z1+sin2αkY+bkZ2+1π2bsinαk2Z1+sin2αkY-bkZ2uid21

Now, the last step to reach our objective is to integrate with respect to α. The two terms of the previous fdp can be integrated similarly. It will be shown how the first one is handled:

I1=α=0πbπ2k2Zsinα1+sin2αkY+bkZ2dα=cosα=x-sinαdα=dxbπ2k2Zx(α=0)x(α=π)-dx1+(1-x2)kY+bkZ2=bπ2k2Z1+kY+bkZ2x(α=0)x(α=π)-dx1-x2kY+bkZ21+kY+bkZ2=xkY+bkZ1+kY+bkZ2=ydx=1+kY+bkZ2kY+bkZdybπ2k2Z1+kY+bkZ2kY+bkZy(x(α=0))y(x(α=π))-dy1-y2=2bπ2k1+kY+bkZ2(kY+b)arctanhkY+bk2Z2+(kY+b)2uid22

The second term can be integrated likewise.

Finally, the target pdf, fk(k), can be written:

fk(k)=2bπ2k1+kY+bkZ2(kY+b)arctanhkY+bk2Z2+(kY+b)2+2bπ2k1+kY-bkZ2(kY-b)arctanhkY-bk2Z2+(kY-b)2,k>0uid23

This is the expression we were looking for. The behaviour of this function is represented in Fig. .

Figure 2.

Probability density function of the modulus of the difference of the cot of the orientation of projected segments (Y=0).

5. The disparity gradient

The disparity gradient has been successfully used in the process of establishment of the correspondence relationships in stereo vision systems. Although the probabilistic behaviour of this feature has been used previously [11], [12], the process to derive some of the pdfs related to the disparity gradient has not been detailed. In this section, we will focus on the specific procedure to find different approximations of the probabilistic characterization of the disparity gradient. Thus, we will derive several expressions of the pdf of the disparity gradientDGrepresents the random variable whereas dgrepresents a realization of DG.:

fDG(dg)uid26

We will pay attention to the assumptions required to derive the pdfs and to the approximations used in the different cases considered.

5.1. Comments on the disparity gradient

The disparity gradient has been successfully used as a source of information to solve the correspondence problem in stereo systems [13], [14], [5], [15], [16], [17], [12].

Generally speaking, the disparity gradient provides a priori information regarding how the real world scene is projected onto the image planes of a stereo system and, consequently, how different matching points in the projected images must be related in terms of geometrical (disparity related) relationships

The disparity refers to the difference between the coordinates of the projections of a certain point of the 3D world onto the image planes of a stereo system. Obviously, the disparity gradient refers to the rate of change of the disparity between nearby or related points [5].

Furthermore, it has been confirmed that the human visual system shows certain limitations related to the disparity gradient when matching stereo images [18]. More specifically, it was proved that 1 represents the limit of the disparity gradient for most of the subjects evaluated. On the other hand, other experiments were performed by other authors that showed that, under certain conditions, the disparity gradient can be over that threshold but with low probability. In fact, Pollard [19] derived a probability function for the disparity gradient in a stereo system with fixation point.

Additionally, the disparity gradient is able to consider other important constraints often employed for the analysis of three dimensional scenes such as figural continuity, ordering of projected features or continuity of the disparity gradient itself [17], [5].

5.2. Stereo system for the probabilistic analysis of the disparity gradient

In the following sections devoted to the probabilistic analysis of the disparity gradient in a parallel binocular stereo system, the specific geometry that will be considered is shown in Fig. . According to this figure, the locations in the real world of the points Aand B, that define a straight segment with its mid-point at (X0,Y0,Z0)and length 2δ, are given by the following expressions:

A=X0+δcosβcosα,Y0+δcosβsinα,Z0-δsinβB=X0-δcosβcosα,Y0-δcosβsinα,Z0+δsinβuid29

Figure 3.

Parallel binocular stereo system for the analysis of the disparity gradient.

Then, the projections of the edge points of the segment onto the right and left image planes are given by:

Ar=-fAz(Ax-b2),-fAzAyBr=-fBz(Bx-b2),-fBzByAl=-fAz(Ax+b2),-fAzAyBl=-fBz(Bx+b2),-fBzByuid31

In this scenario, the disparity gradient is defined as the quotient between the difference of disparity between the two points observed and their Cyclopean separation [19]:

dg=DifferenceofdisparityCyclopeanseparationuid32

Taking into account that the Cyclopean projections of Aand Bare given by the following equation:

Ar+Al2andBr+Bl2uid33

and using the disparity vectors associated to the points Aand Bgiven by

(Al-Ar)and(Bl-Br)uid34

respectively. Then the disparity gradient can be written as follows:

dg=2||(Ar-Br)-(Al-Bl)||||(Ar-Br)+(Al-Bl)||uid35

Now, by substitution of the expressions of Al, Bl, Arand Br, multiplying by AzBz, substituting by their expressions in terms of δ, βand Z0, after some simplifications and reordering all the terms, the following expression is found:

dg=|bsinβ|||-X0sinβ-Z0cosβcosα,-Y0sinβ-Z0cosβsinα||uid36

This is the main equation that will be used to derive different expressions of the disparity gradient in different scenarios.

The following sections describe the scenarios and the procedures issued to derive the different probability density functions.

5.3. Primitives centred in the world reference system

In our first scenario, we will be able to derive an exact analytical expression of the pdf of the disparity gradient This expression can be considered to be illustrative of the behaviour of dg. Moreover, in the next subsection, we will show how the same expression is found under different conditions and assumptions.

In this first scenario, we will assume that X0=0, Y0=0and α=0(see Fig. ). Then, the expression of the disparity gradient (eq. ()) is readily simplified to give:

dg=bZ0|tanβ|uid38

We will assume that the angle of orientation βbehaves as a uniform random variable in the range (0,π).

Paying attention to the symmetry of dg, it is possible to pose the problem in a more convenient way. Without loss of generality, the modulus of tanβin eq. () can be removed by simply allowing the random variable βto be defined as a uniform random variable in (0,π2). The application of this and other symmetry conditions that will be considered later will allow us to avoid some expressions that involve the calculation of the modulus of certain functions and thus the analysis and some of the expressions involved will remain conveniently more simple.

According to equation (), it is quite simple to obtain the derivative of the disparity gradient with respect to β. Let g(β)=dg, then g'(β)=bZ0cos2β. On the other hand, it is possible to obtain βas g-1(dg)=arctanZ0bdg. Thus, finally, the pdf of DGis directly obtained:

fDG(dg)=π2bZ01cos2ββ=g-1(dg)=2Z0πb1+tan2ββ=g-1(dg)=2Z0πbtanarctanZ0bdg2+1=
fDG(dg)=2πbZ0dg2+bZ02,dg(0,)uid39

Figure 4.

Probability density function of the disparity gradient when the primitives projected are centred in the world reference system.

Figure 5.

Distribution function of the disparity gradient when the primitives projected are centred in the world reference system.

In this expression (eq. () and Fig. ), a unilateral Cauchy probability density function should be identified. In our scenario, this Cauchy function is tuned by the parameters 0 and bZ0[20]. The distribution function can be easily found (See Fig. ):

FDG(dg)=2πarctanZ0bdg,dg(0,)uid42

5.4. Narrow field of view cameras

In this section, another step in the analysis of the behaviour of the disparity gradient will be done. We will consider a binocular stereo system with cameras of narrow field of view satisfying the epipolar constraint. This is a scenario that can be applied in numerous cases. Moreover, we can consider this scenario as a basic model for the analysis of stereo systems and suitable for practical applications.

In this scenario, the disparity gradient is given by:

dg=2||fbAz,0-fbBz,0||||-2fAz(X0+δcosβcosα),-2fAz(Y0+δcosβsinα)......--2fBz(X0-δcosβcosα),-2fBz(Y0-δcosβsinα)||uid44

After the substitution of Azand Bzby their respective expressions in terms of X0, Y0, Z0, α, βand δand reordering all the terms the following expression can be found:

dg=b2sin2β(X02+Y02)sin2β+Z02cos2β+2Z0sinβcosβ(X0cosα+Y0sinα)uid45

We will derive the desired pdf making use of this equation.

The fact that the cameras of the stereo system have a narrow field of view implies that the coordinates in the real world of the projected objects should satisfy the following condition: Z0X0,Y0. On the other hand, the angle βshould not be equal to π2(as a matter of fact, being βa continuous random variable, this conditions represents and event with zero probability).

Under the hypotheses described, removing X0and Y0from the expression of the disparity gradient, because of the narrow field approximation, and assuming that Z0Z02, the following simplified expression is found:

dgb2sin2βZ02cos2β=bsinβZ0|cosβ|uid46

In this scenario, the symmetry of the geometry and the behaviour of the random variables αand βallows us to consider the following range for the uniform random variables αand β: (-π2,π2)and (0,π2), respectively. And then, the expression of the disparity gradient can be written as:

dg=bsinβZ0cosβuid47

Now, in order to derive the behaviour of the disparity gradient, we will observe the region in which the random variable DGis smaller than a certain value dg. Then, Prob{DG<dg}is given by the probability that the random variables αand βare such that DG<dg. Let Cdgdenote the region in the α-βplane that complies with this condition:

Prob{DG<dg}=Prob{(α,β)Cdg}uid48

This probability can be easily found by integrating the joint pdf of αand βin the region Cdg:

FDG(dg)=Cdgfα,β(α,β)dαdβuid49

where, according to the selected hypotheses, the joint pdf required is given by fα,β(α,β)=2π2.

In order to define the region Cdg, eq. () must be used in order to obtain the solutions of β:

β=arctandgZ0buid50

So, the region in the α-βplane that defines Cdgis given by the following relations:

α-π2,π2β0,arctandgZ0buid51

Thus, it is possible to derive the probability distribution function of the disparity gradient solving the following integral:

FDG(dg)=α=-π2α=π2β=0β=arctandgZ0b2π2dβdαuid52

which is given by:

FDG(dg)=2πarctanZ0bdguid53

Then, the probability density function can be readily obtained:

fDG(dg)=2πbZ0dg2+bZ02uid54

Observe that, under different conditions and hypotheses, the same expressions for the behaviour of the disparity gradient as in the case of primitives centred in the world coordinate system (Sec. ) have been obtained. Of course, this fact comes from the assumption that Z0X0,Y0which asymptotically leads to the more specific case in which X0=0and Y0=0.

5.5. General case. Approximate expression

Under general conditions, a close analytic solution for the probability density function or the probability distribution function of the disparity gradient has not been found. So, we will face the derivation of an approximate solution.

To this end, consider the following approximate expression of the disparity gradient in our stereo system (Fig. ):

dg=bX02+Y02+Z02cot2β+2Z0cotβK(X0+Y0)uid56

In this expression, obtained after eq. (), the terms (X0cosα+Y0sinα)have been substituted by K(X0+Y0). Note that Kshould not modify the region in which the disparity gradient is properly defined: DG[0,). Using this idea, it is possible to arrive at the desired goal. Now the procedure is described.

We know that if β0, then dg0. So, we can find a condition to impose on Kso that maxDG. To this end, the minimum of the denominator in eq. () can be found in the usual way, deriving the expression in the square root with respect to βand finding the roots:

βX02+Y02+Z02cot2β+2Z0cotβK(X0+Y0)=0uid57
-2Z02cotβcsc2β-2Z0csc2βK(X0+Y0)=0uid58

Now, since cscβ0β, the following must be fulfilled:

Z0cotβ+K(X0+Y0)=0uid59

Thus, the following relation is found:

cotβ=-K(X0+Y0)Z0uid60

Recall that in the minimum the denominator in eq. () must be zero. Substituting cotβaccording to the previous expression in the denominator of eq. (), the following must be fulfilled:

X02+Y02+Z02-K(X0+Y0)Z02+2Z0-K(X0+Y0)Z0K(X0+Y0)=0uid61

which leads to the following expression:

K=X02+Y02(X0+Y0)2uid62

Thus, the approximation of the disparity gradient that will be used is given by:

dgbX02+Y02+Z02cot2β+2Z0X02+Y02cotβuid63

Now, the probability distribution function will be found. Consider Cdgas the region in which DG<dgand let Cdg(α,β)denote the region in the α-βplane such that DG<dg. Then, again:

FDG(dg)=Cdg(α,β)fα,β(α,β)dαdβuid64

(a)

(b)

(c)

(d)

Figure 6.

Probability distribution function of the disparity gradient {1}. General case: simulation results (solid line) and analytic approximation (dashed line).

Since DGdoes not depend on α(eq. ()), the region Cdg(α,β)can be defined as a function of β, exclusively:

FDG(dg)=Cdg(β)αfα,β(α,β)dαdβ=Cdg(β)1πdβuid66

In order to define Cdg(β), dgmust also be written as a function of β; the following result if easily obtained:

cotβ=-X02+Y02Z0±bdgZ0uid67

Let β1and β2represent the two solutions of this equation, then the region Cdg(β)is defined by the following intervals:

Cdg(β)=(-π2,min(β1,β2))(max(β1,β2),π2)uid68

With all this, the desired solution, the probability distribution function of the disparity gradient, is given by (Figs. and ):

FDG(dg)=1-1πarccot-X02+Y02Z0-bdgZ0-arccot-X02+Y02Z0+bdgZ0uid69

(e)

(f)

(g)

(h)

Figure 7.

Probability distribution function of the disparity gradient {2}. General case: simulation results (solid line) and analytic approximation (dashed line).

Note that this solution is mathematically correct, however some considerations must be taken into account so that FDG(dg)behaves as a proper probability distribution function sec. 2.2[10]. Specifically, the function arccotreturns an angular value which, ultimately, can be seen as a periodic function with period π. This means that there is an infinite number of solutions of arccot, although the main solution is often considered to be in the interval (-π2,π2). In our specific development, the function derived behaves properly if the solutions of the function arccotare selected in the range (-π,0).

After the probability distribution function (eq. ()), the probability density function (pdf) of the disparity gradient is readily found [9]:

fDG(dg)=1π2bZdg2Z02+b2+dg2(X02+Y02)dg4Z04+b4+dg4(X02+Y02)2+2dg2Z02b2......+2dg4Z02(X02+Y02)-2b2dg2(X02+Y02)uid71

which is a usable expression of the pdf of the disparity gradient that completes the analysis of the probabilistic behaviour of this parameter under the conditions and hypotheses selected.

6. Concluding summary

In this chapter, we have dealt with the probabilistic behaviour of certain relations established between the projection of features onto the image planes of a parallel stereo system. Specifically, we have considered relations between the orientation of projected edgels and the disparity gradient.

The projected edgels are simple features that can be considered in a matching stage [7]. The relation between their orientations constitutes an a priori source of information that, using the models proposed, can be used in the matching processes [21] of stereo systems. The formulae of the relation between the orientation of the projections derived are perfectly suited for application in Bayesian models for stereo matching [22].

The disparity gradient is an important parameter for stereo matching systems [14]. In this chapter, it has been analysed under different conditions to find proper probability density functions usable in a probabilistic context.

The functions derived can be used alone to match random dot stereo pairs [23], [24], [17], [25]. Also, these functions can contribute and collaborate with other matching models in the solution of the correspondence problem in stereo systems. Specifically, Bayesian approaches can be employed to solve the correspondence problem [26] using the proposed models of the disparity gradient [12].

Acknowledgements

This work was supported by the Ministerio de Economía y Competitividad of the Spanish Government under Project No. TIN2010-21089-C03-02 and Project No. IPT-2011-0885-430000.

© 2012 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Lorenzo J. Tardón, Isabel Barbancho and Carlos Alberola-López (July 11th 2012). Probabilistic Analysis of Projected Features in Binocular Stereo, Current Advancements in Stereo Vision, Asim Bhatti, IntechOpen, DOI: 10.5772/46027. Available from:

chapter statistics

1141total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Stereo Algorithm with Anisotropic Reaction-Diffusion Systems

By Atsushi Nomura, Koichi Okada, Hidetoshi Miike, Yoshiki Mizukami, Makoto Ichikawa and Tatsunari Sakurai

Related Book

First chapter

Calibration and Sensitivity Analysis of a Stereo Vision-Based Driver Assistance System

By Andras Bodis-Szomoru, Tamas Daboczi and Zoltan Fazekas

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us