Chapter 2 CFA Interpolation for Color Reproduction of Bayer Mosaic Images
2.5 Summary
A novel hard-decision CFA interpolation procedure has been developed based on the spectral-spatial correlation of a Bayer mosaic image. The proposed HPHD CFA interpolation method effectively reconstructs fine detail features in both edge and texture regions of demosaiced images. One merit of the proposed HPHD algorithm is that it can combine with many existing image interpolation methods such as decision-based algorithm (set α =1), edge-directed interpolation, adaptive interpolation, linear interpolation, etc for improved performance. Moreover, the proposed heterogeneity-projection scheme provides an efficient method for decision-based algorithms to make accurate direction-selection before performing interpolation.
In order to reconstruct the demosaiced images with fewer color artifacts, a novel CDEA CFA interpolation method is then proposed to combine with the HPHD algorithm. The proposed CDEA algorithm effectively reduces color artifacts in both smooth and edge regions of demosaiced images. Furthermore, any existing image interpolation method can be
combined with the proposed CDEA algorithm to reconstruct the green channel. A new edge-adaptive interpolation method is also presented by adopting the hard-decision rule from HPHD algorithm to reconstruct the green channel from CFA samples. In Chapter 5, the performance of the proposed HPHD-CDEA method will be compared with three renowned CFA interpolation methods. Experimental results will show that HPHD-CDEA method not only outperforms all of them in PSNR (dB) and S-CIELAB ∆E*ab measures, but also gives superior demosaiced fidelities in visual comparison.
Chapter 3
Robust Mobile Robot Visual Tracking Control Based on a Dual-Jacobian Visual Interaction Model
3.1 Introduction
From the literature discussed in Chapter 1, we have noted that a challenge in mobile robot visual tracking control design is to develop a visual tracking control system to track a dynamic moving target based on a stability criterion and overcome the uncertainties encountered in practical systems. This problem motivates us to derive a new model for designing a robust VTC to solve the visual tracking problem of dynamic moving target and overcome the internal disturbances of practical system (such as system parametric uncertainty and velocity quantization uncertainty). To achieve this, this chapter presents a novel dual-Jacobian visual interaction model to help the design of a robust VTC for a wheeled mobile robot equipped with a tilt camera. The proposed design enhances various image-tracking applications using an on-board monocular camera, such as human-robot interaction and surveillance. Based on Lyapunov theory, the proposed control scheme not only possesses some degree of robustness against parametric uncertainty, but also overcomes the external uncertainty caused by velocity quantization noise. Moreover, the proposed controller fully works in image space; hence the computational complexity and the effects of sensor/camera modeling errors can be greatly reduced. The main differences between the proposed VTC and other existent approaches are summarized as follows:
1) The proposed dual-Jacobian visual interaction model considers not only the effect of
mobile robot motion, but also the effect of target motion. Thus, based on the proposed model, the visual tracking control problem of a unicycle-modeled mobile robot for tracking a dynamic moving target can be solved with exponential convergence using a single controller. Moreover, the proposed model also considers the kinematics of a tilt camera platform mounted on the mobile robot. Therefore, the applicability of the proposed method is greatly increased.
2) The proposed visual tracking control system not only possesses some degree of robustness against the system model uncertainties, but also overcomes the unmodelled quantization effect in the velocity commands and the occlusion effect during visual tracking process.
This advantage enhances the reliability of the proposed method in practical applications.
3) The proposed visual tracking control system works fully in image space. Therefore, compared with position-based [23], homography-based [13, 14], and epipole-based [4, 15, 16, 17] visual tracking control approaches, the computational complexity and the sensor/camera modeling errors can be reduced due to the advantages of image-based visual servo control [2].
The basic assumptions of the proposed method are listed as follows:
1) The on-board camera is supposed to be a calibrated pinhole camera. Because the proposed VTC possesses some degree of robustness against parametric uncertainty, a simple linear camera calibration method [44] can be used to estimate the intrinsic parameters of the camera.
2) The width of target is supposed to be a priori known constant in order to simplify the depth estimation problem. However, this is not a necessary assumption for the proposed method.
Any algorithm or sensor which provides the depth information can be utilized to combine with the proposed method.
The rest of this chapter is organized as follows. Section 3.2 describes the system modeling of visual tracking control problem and the proposed dual-Jacobian visual interaction model accordingly. Section 3.3 presents the results of VTC design. In Section 3.4, the robustness of the control system against the system model uncertainty is analyzed and discussed. Section 3.5 presents the design of the robust control law to overcome the velocity quantization error encountered in practical systems. Section 3.6 summarizes the contributions of this design. Experimental results will be reported in Chapter 5. Several interesting experimental observations will be presented and discussed.
3.2 Camera-Object Visual Interaction Model
This section derives the visual interaction model between a mobile robot and a dynamic moving target. We first introduce the kinematics model of wheeled mobile robot and target used in this design. The mathematic derivations of the proposed model are then presented and explained.
3.2.1 Kinematics Model of Wheeled Mobile Robot and Target
Figure 3-1 shows the model of wheeled mobile robot and target considered in the nonholonomic visual tracking control problem. The wheeled mobile robot equips with a tilt camera to track a dynamic motion target, which is supposed to be a well-recognizable object with appropriate dimensions in the image plane and can only translate with respect to the robot. The tilt camera is mounted on top of the mobile robot and its optical-axis faces the target of interest, for instance, a human face. Figure 3-1(a) illustrates a model of the wheeled mobile robot and the target in the world coordinate frame Ff (see Fig. 3-2), in which the motion of the target is supposed to be holonomic such that
t f t
f
V
X
& = , (3.1)
(a) (b)
Fig. 3-1: (a) A model of the wheeled mobile robot and the target in the world coordinate frame. (b) Side view of the wheeled mobile robot with a tilt camera mounted on top of it.
where
X
tf =[x
tfy
tfz
tf]T andV
ft =[v
xfv
yfv
zf]Tdenote, respectively, the position and velocity of target in the world coordinates.Figure 3-1(b) is the side view of the scenario under consideration, in which the tilt angle φ gives the relationship between the camera coordinate frame
F and the mobile coordinate
c frameF . The kinematics of the wheeled mobile robot is described by [45]
m⎥⎥
⎥
⎦
⎤
⎢⎢
⎢
⎣
⎡
=
m f m f
m f m f m f
v v X
θ θ
cos 0 sin
& , and
m t
m f m f
w w
=
= φ θ
&
&
, (3.2)
where
X
mf =[x
mfy
mfz
mf]T is the position of mobile robot in the world coordinates, ),
(θmf φ are the orientation angle of mobile robot and the tilt angle of onboard camera, wtm is the tilt velocity of the camera, and (
v
mf,w
mf) are the linear and angular velocities of mobile robot. In practice, (v
mf,w
mf ) can be used to calculate the velocity of each wheel of the mobile robot such that2
where (vlm,
v ) are the left- and right-wheel velocities, respectively, and D represents the
mr distance between the two drive wheels. In the rest of this chapter, the target model (3.1) and mobile robot model (3.2) will be utilized to derive the visual interaction model and to design the visual tracking control system.3.2.2 Coordinate Transformation from World Frame into Camera Frame
Figure 3-2 illustrates the relationship between coordinate systems, namely the world, camera and image coordinate frames. Let
X
f =X
tf −X
mf denote the related position between mobile robot and the target in the world coordinate frame. In order to describe a mobile robot interacting with the target in the image coordinate frame, a visual interaction model has been derived by transferring the kinematics of Xf from the world coordinate frame into the image coordinate frame. This subsection presents the transformation of the kinematics of Xf from the world coordinate frame into the camera coordinate frame.As shown in Fig. 3-2, Xc =[xc yc zc]T denotes the related position in the camera coordinate frame and can be calculated by the coordinate transformation such that
δY
is the distance between the center of robot head and the onboard camera. Because δY is a constant translation vector, the derivative of (3.4) becomes
Fig. 3-2: World, camera and image coordinate frames of robotic visual interaction.
Substituting (3.1), (3.2) and (3.4) into (3.5), the kinematics of the interaction between robot and target in camera frame can be obtained by taking the derivative of Xc such that
[ ]
⎥⎥ the following, the kinematics model (3.6) will be used to derive the interaction model in the image coordinate frame.
3.2.3 Coordinate Transformation from Camera Frame into Image Frame
In this subsection, the related position Xc is transformed into the image coordinate frame for deriving the visual interaction model based on (3.6). We first define the system state in the image coordinate frame for the controller design. Fig. 3-3 illustrates the definition of observed system state in the image plane. In Fig. 3-3,
x and
iy are, respectively, the horizontal and
i vertical position of the centroid of target in the image plane, andd is the width of target in
x the image plane. Similar to the human’s visual tracking behavior, the purpose of the visual tracking control design is to control the centroid position and width of target from an initial state into the desired state in the image plane.In the following, the visual interaction model is derived by (3.6) and the selected system state. Based on the pinhole camera model, the diffeomorphism (please see [46] for the specification) in the image plane can be defined by the standard projection equations [44]
such that:
Fig. 3-3: The definition of observed and desired system states in the image plane.
where (fx,fy) represent fixed focal length along the image x-axis and y-axis, respectively [47], and W denotes the actual width of the target. By taking the derivative of (3.7), the kinematic relationship between image and camera coordinate frames can be found such that
c
Substituting (3.6) and (3.7) into (3.8), the kinematic relationship between robot and target in the image coordinate frame can be modeled by quasi-linear parameter-varying (Quasi-LPV) description [48] such that
i
⎥⎥ elements of system matrix
A and vector
iC are time-varying function dependent on the
i robot’s pose and target’s velocity; and the elements of control matrixB are time-varying
i function dependent on the robot’s pose and current system state.3.2.4 Dual-Jacobian Visual Interaction Model
The visual interaction model (3.9) indicates that the elements of system matrix
A and
i vectorC are function of target velocity. Thus, expression (3.9) can be rewritten such that
iu
Expression (3.10) shows that the visual interaction model consists of two parts: first, the effect of target motion
[ ]
i ft& . Thus, (3.10) can be rewritten as a dual-Jacobian equation such
Fig. 3-4: Depicts the concept of dual-Jacobian visual interaction model (3.11).
that
u V
X X
X
&i = &it+ &im =J
i ft +B
i , (3.11)where matrix
J , which is named target image Jacobian, transfers the target velocity
iV
ft into target image velocity X⁢ matrix Bi, which denotes robot image Jacobian, transfers the mobile robot control velocityu into robot image velocity
X&im. In other words, the image velocity X&i is caused by the combination of target image velocity X&it and robot imagevelocity X&im. Figure 3-4 shows the concept of dual-Jacobian equation (3.11). Therefore, the visual interaction between robot and target in the image coordinate frame can be modeled as a dual-Jacobian visual interaction model (3.11), which combines the motion effect of mobile robot with moving target together.
Remark 3.1: The scalars kx =fx zc and ky =fy zc in (3.7) depend on the depth information between camera and target. The estimation of depth information is a demanding task in visual tracking control design; especially only one camera is used. Thus, an algorithm or sensor which provides the depth information is usually required during visual tracking
process. In order to simplify the depth estimation problem, an alternative is to assume that the width of target is known a priori. Therefore, the scalars
k and
x ky can be calculated usingthe state variable
d directly based on the fact that
x kx =fx zc =dx W and ky =kxfy fx, where W denotes the width of target for a specific target.3.3 Visual Tracking Controller (VTC) Design
In this section, a visual tracking control law based on the proposed dual-Jacobian visual interaction model (3.11) for tracking a target of interest in the image plane is derived by exploiting feedback linearization and pole placement approaches.
3.3.1 Error Coordinate Transformation
In order to control the system from an initial state to the desired state, an error-state model will be helpful for us to design the tracking controller. Define the error state in the image plane such that
[ ] [
i i i i x x]
TT e e e
e
x y d x y d
X
= = x − * y − * d − * , (3.12)where Xi =
[
xi yi dx]
T is the vector of fixed desired states in the image plane;[
i i x]
Ti x y d
X* = * * * is the vector of estimated states from the VSE (see Chapter 4). Based on the error state (3.12), the dynamic error state model in the image plane can be derived directly by taking the derivative of (3.12) such that
u V X
X
X
&e =− &it− &im =−J
i ft −B
i . (3.13) With the new coordinateX , the visual tracking control problem is transformed into a
e stability problem. IfX converges to zero, then the visual tracking control problem is
e solved.3.3.2 Visual Feedback Control Design
Based on the dynamic error state model (3.13), we choose the feedback linearization control law such that
) where Kg is an 3-by-3 gain matrix. Substituting (3.14) into (3.13) yields
e g
e
X
X
& =−K
. (3.15) Next, we choose the gain matrix such that) Suppose that the initial error-state
X
e(t
0) is within the image plane. Then expression (3.17) indicates that the following inequality:) value of the vector X. From (3.19), it is clear that the system error-state satisfies
) control problem is solved. Summarizing the above discussions, we obtain the following theorem.
Theorem 3.1: Suppose the the initial system state Xi is within the image plane. Let
0
( 1 2 3 > be three positive constants. Consider the closed-loop visual tracking system (3.13). If the matrix
B
i is nonsingular, then the closed-loop visual tracking system (3.13) is exponentially stable by using the control law)
Proof: Consider the closed-loop visual tracking system (3.13). We first define a
positive-definite Lyapunov function associated with the system error-state )
Taking the derivative of (3.21) yields
) that if f(u)>0 then the equilibrium point of (3.13) is asymptotically stable. Substituting the control law (3.20) into f(u), we then have
e symmetric positive definite (SPD) matrix, the following inequality holds:
0
where )λmin(
A
denotes the minimum eigenvalue of matrix A. Expression (3.24) concludes that the closed-loop visual tracking system (3.13) is exponentially stable and hence completesthe proof. ■
Remark 3.2: Although the proposed image control law (3.20) results in a smooth convergence in the image plane, it still has to proof that the robot should have followed the target. The discussion of this problem is presented in Appendix D.
3.3.3 Singularity Analysis
The feedback linearization control law (3.20) poses a singularity problem of matrix
B
i. LetB denote an element of matrix
mnB
i corresponding to the m-th row and n-th column.By directly computing, the determinant of matrix
B
i is given by32 Based on (3.25), the singularity condition of matrix
B
i can be found such thatφ following geometric relationship:
y
equals to 0 or π. The physical meaning of this is that the target is directly above or directlyFig. 3-5: Physical meaning of the singularity condition (3.26).
below the robot, and the robot will be unable to approach the target in any way due to insufficient degrees-of-freedom. Therefore, the robot will stop tracking temporarily under such circumstances.
Remark 3.3: The proposed visual interaction model (3.11) poses a question that the derivation of dx =kxW is not always hold for a cylindrical target used in the system modeling. As shown in Fig. 3-6, parameters W, fx and zc remain the same, but the cylindrical target is shifted along xc. The projection dx is not the same and includes an error δ . Because the camera is dx supposed to have a limited field of view, the error δ would be small and could be modeled dx as a system uncertainty. In the next section, the robustness analysis is presented in order to handle this kind of uncertainty.
3.4 Robustness Against System Model Uncertainty
In this section, the robustness of the VTC (3.20) is investigated against model uncertainties from camera parameters (fx,fy) , robot parameters (θmf ,φ) , and target parameters )(W,
x
&it,y
&it,d
&xt , etc. Consider the following closed-loop visual tracking system with parametric uncertainties:Fig. 3-6: Projection error in dx.
u V
u V
X
&e =−J
i ft −B
i =−(J
i+δJ
i) ft −(B
i +δB
i) , (3.30) whereδJ
i and δB
i are unknown bounded matrix-disturbances. Recall the positive-definite Lyapunov function defined in (3.21), the derivative of (3.21) with parametric uncertainties becomes) (
| )]
( ) ( [ ] [
|
|(3.30)
X X
(3.30)X V X u f u f u V
(3.13)f u
V
& = eT &e =− eTJ
i ft + eTB
i =− +δ = & −δ , (3.31)where δ
f
(u
)=X
eTδJ
iV
ft +X
eTδB
iu
is unknown. Assume that δJ
i and δB
i are bounded and there exists two positive valuesC and
MC such that
Ne N e
M X C X
C u
f( )≤ 2 +
δ . (3.32)
We now introduce the following definition.
Definition 3.1: The system (3.13) is said to be practically stable at the origin if a compact set
S in
ℜn exists containing the origin such that for eachX
e(t
0)∈S
,X
e(t
)∈S
for allt
0t
≥ .In practical applications, the practical stability problem can be more easily handled. Now, the
main result is presented as follows.
Theorem 3.2: Consider the closed-loop visual tracking system (3.13) with unknown bounded
parametric uncertainties δ
J
i and δB
i defined in (3.30). LetC and
MC be two positive
N values defined in (3.32) and assume that the target’s motion velocity Vft ≥0 is bounded and small enough. Choose the controller u as given in expression (3.20) with gain matrix0
Then the origin is asymptotically stable under the condition Vft =0, and the origin is
practically stable under the condition Vft >0. ■ Proof: Choose the controller u in (3.20) with parametric uncertainties such that
)
maximum eigenvalue of matrix A, and A denotes the Euclidean norm value of the matrix
A. From (3.31), (3.33) and (3.34), it follows that
.
Expression (3.35) implies that if [λmin(
K
g)−C
M]X
e −C
N >0 can be guaranteed, then 0|(3.13)<
V&
is satisfied and thus the system has the robust property w.r.t. the parametric uncertainties.From (3.34), the condition [λmin(
K
g)−C
M]X
e −C
N >0 can be rewritten such that[
g i i] [
g]
NN M g
e C C C
X >[λmin(K )− ]−1 = ρλ(K )− δB B−1 −1 λmax(K )−1 , (3.36) where ( )
[
( )]
min( )1
max g g
g
K K
K
λ λρλ = − . In general, we have that δBiBi−1 <1 . By
assumption (i) of the theorem, expression (3.36) follows that when ρλ(Kg)→1 satisfies,
>0
∃ε such that
[ ] [ ]
{
ζ ∈ℜ ζ < ρλ − δ λ +ε}
→ n g i i− − g − N
e C
X : (K ) B B 1 1 max(K ) 1 . (3.37) Because of CN = δJi −δBiBi−1Ji Vft , the condition Vft =0 follows that
ε
e <
X
andV
&|(3.13)≤−[
ρλ(K
g)− δB
iB
i−1] X
e 2 <0. (3.38) Thus, there exists a constant δ >0, δ <ε, ∃T such thatδ ε ⇒ <
< ( ) )
(
t
0X t
X
e e ∀t
≥t
0 +T
, (3.39) which means that the origin is asymptotically stable under the condition Vft =0.On the other hand, the condition Vft >0 follows that when ρλ(Kg)→1 satisfies,
>0
∃ε such that
0
|(3.13)<
V&
for all Xe <[
ρλ(Kg)− δBiBi−1]
−1[
λmax(Kg)]
−1CN +ε . (3.40) Since Vft is supposed to be small enough, CN ∝ Vft is sufficiently small. By assumption (ii) of the theorem, there exists a sufficiently large λmax(Kg)and a constant δ1 >0 such thatc
0X
e < , where c0 =[
ρλ(Kg)− δBiBi−1]
−1δ1+ε. Choose the setS
={
ζ ∈ℜn : ζ ≤c
1}
,where 0
c
1 > satisfiesc
1 > . Then, any trajectory that begins at anyc
0X
e(t
0)∈S
will remain in S for allt
≥ under the control law (3.20). Thus, the origin is practically stablet
0 under the condition Vft >0. This completes the proof. ■In Chapter 5, the result of Theorem 3.2 will be validated by practical experiments.
Further, in realization of the control schemes, it was noted that the quantization error in velocity commands degrade the performance of the controller and might make the system unstable. It is therefore interesting to study the robustness issues related to velocity quantization uncertainty. In the following section, a robust control law based on Lyapunov’s direct method will be derived to overcome the velocity quantization uncertainty in practical control systems.
3.5 Robustness Against Velocity Quantization Error
When tracking a target, it is desirable for the robot to have a smooth motion in human-robot interaction. But in such circumstances, one will face the problem caused by velocity quantization error in practical implementation. In this section, a robust control law is derived to eliminate the velocity quantization error encountered in practical control systems based on the dynamic error state model defined in (3.13). To do so, a stability necessary condition (SNC) is first derived for ensuring asymptotic stability and practical stability of the closed-loop visual tracking system through Lyapunov’s direct method. The robust control law is then proposed to guarantee that the visual tracking system satisfies the SNC and hence complete the controller design.
3.5.1 Stability necessary condition (SNC)
Digital control systems usually have uniform quantization errors due to the finite-length
effects on the sample values [49]. In other words, the ideal (theoretical) control command u is quantized such that
u u
u = +δ , (3.41) where u denotes the practical (actual) control command sent to the robot actuator, and δ u represents the uniform quantization error encountered in the system. Thus, in practice, (3.22) becomes
0 ) (
| )]
( ) ( ) ( [ ) (
| (3.13)
) 41 . 3 (3.13)
( =−
f u
=−f u
−f u
+f u
=V
−f u
<V
& δ δ δ & δ δ , (3.42)where δf(δu)= XeTBiδu−XeTδBiδu≅ XeTBiδu . Expression (3.42) shows that if the controller u satisfies the assumptions (i) and (ii) of Theorem 3.2, then the equilibrium point of the system (3.13) still can be unstable under the condition δ
f
(δu
)<0. Therefore, one has the following SNC in practical control implementation:SNC: Consider the closed-loop visual tracking system (3.13) with unknown bounded parametric uncertainties δ
J
i and δB
i defined in (3.30). Let δ denote the velocity output u quantization error in practical systems. Suppose that the controller u given in (3.20) satisfies the assumptions (i) and (ii) of Theorem 3.2. Then, the result of Theorem 3.2 holds under thecondition δf(δu)= XeTBiδu ≥0 ■
3.5.2 Proposed robust control law
SNC implies that a practical system may become unstable under the condition 0
) (
u
<f
δδ . Our goal is to design a robust control law which not only guarantees SNC to be always satisfied but also increases the convergence rate of the control system.
First, we expand δ
f
( uδ ) such thatFirst, we expand δ