26.06.2013 Views

Mouvement et Vidéo - Centre de Morphologie Mathématique

Mouvement et Vidéo - Centre de Morphologie Mathématique

Mouvement et Vidéo - Centre de Morphologie Mathématique

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Htotal nspl <br />

Htotal (λ/λ0, nreg) <br />

Htotal nreg <br />

λ/λ0 = 0.25 <br />

Htotal nreg <br />

λ/λ0 = 0.125 <br />

Herr Hv <br />

Htotal Hv <br />

Htotal ∆v <br />

φ00 <br />

xij φ00


R(ϑ, ρ) ρ = 0.96 ρ = 0.97 ρ = 0.98 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

P3,0 <br />

P3,1


R(ϑ, ρ) ρ = 0.96 ρ = 0.97 ρ = 0.98 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

P3,0 <br />

P3,1


R(ϑ, ρ) ρ = 0.96 ρ = 0.97 ρ = 0.98 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

P3,0 <br />

P3,1


∗ 0<br />

r∗ 0<br />

r∗ 0<br />

r∗ 0<br />

r∗ 0<br />

r∗ 0<br />

= 6.7 <br />

= 6.7 <br />

= 6.7 <br />

= 8.0 <br />

= 8.0 <br />

= 8.0 <br />

(αl)


∗ 0<br />

r∗ 0<br />

r∗ 0<br />

r∗ 0<br />

r∗ 0<br />

r∗ 0<br />

= 6.7 <br />

= 6.7 <br />

= 6.7 <br />

= 8.0 <br />

= 8.0 <br />

= 8.0 <br />

(αl)


x x<br />

v v<br />

M M<br />

s s<br />

<br />

Id <br />

diag(λ1, λ2, . . . , λn) λ1, λ2, . . . , λn<br />

t v v<br />

t M M<br />

vp M M<br />

tr M M<br />

<br />

F <br />

W <br />

Hf f<br />

∇f f<br />

˜z z<br />

1E<br />

E<br />

card A A<br />

δ <br />

f f<br />

f ⋆ g f g<br />

w w<br />

<br />

med(X) X<br />

E X X<br />

arrondi(t) t <br />

arrondi(x) x <br />

H(X) X<br />

Hd(X) X<br />

I(X; I) X Y


R <br />

Z <br />

N <br />

ℓ2(Z) <br />

L2(R)


R <br />

Z <br />

N <br />

ℓ2(Z) <br />

L2(R)


R <br />

Z <br />

N <br />

ℓ2(Z) <br />

L2(R)


R <br />

Z <br />

N <br />

ℓ2(Z) <br />

L2(R)


I(t, x)<br />

t p <br />

<br />

d<br />

∂I<br />

I(t, x(t)) = + v · ∇I,<br />

dt ∂t<br />

v(x) x(t) <br />

<br />

<br />

∂I<br />

∂t + v · ∇I = 0


I(t, x)<br />

t p <br />

<br />

d<br />

∂I<br />

I(t, x(t)) = + v · ∇I,<br />

dt ∂t<br />

v(x) x(t) <br />

<br />

<br />

∂I<br />

∂t + v · ∇I = 0


I(t, x)<br />

t p <br />

<br />

d<br />

∂I<br />

I(t, x(t)) = + v · ∇I,<br />

dt ∂t<br />

v(x) x(t) <br />

<br />

<br />

∂I<br />

∂t + v · ∇I = 0


I(t, x)<br />

t p <br />

<br />

d<br />

∂I<br />

I(t, x(t)) = + v · ∇I,<br />

dt ∂t<br />

v(x) x(t) <br />

<br />

<br />

∂I<br />

∂t + v · ∇I = 0


I(t, x)<br />

t p <br />

<br />

d<br />

∂I<br />

I(t, x(t)) = + v · ∇I,<br />

dt ∂t<br />

v(x) x(t) <br />

<br />

<br />

∂I<br />

∂t + v · ∇I = 0


f <br />

<br />

<br />

f(t) =<br />

R<br />

f(u)δ(t − u)du. <br />

<br />

<br />

<br />

f(t) =<br />

R<br />

f(ω)e iωt du. <br />

<br />

f(t) f t <br />

<br />

f(ω) f ω f(t) <br />

f f(ω) <br />

<br />

<br />

<br />

t ↦→ δ(t − u) <br />

<br />

t ↦→ e iωt


f <br />

L2 1 <br />

|f(t)| 2 dt = 1.<br />

c(f) ∆(f) <br />

<br />

c(f) = t|f(t)| 2 dt,<br />

<br />

<br />

∆(f) =<br />

(t − c(f)) |f(t)| 2 dt.<br />

∆(()f) <br />

<br />

f L2 1 <br />

∆(f)∆( f) ≥ 1<br />

. <br />

2<br />

∆( f) f <br />

<br />

<br />

<br />

<br />

Gt0,ω0,∆t(t) = Ae − (t−t0 )2<br />

2δt2 e i∆0t<br />

,<br />

A Gt0,ω0,∆t <br />

<br />

<br />

<br />

(t, ω) (c(f), c( f)) (∆(f), ∆( f)) <br />

f L2(R)


ω<br />

<br />

t<br />

ω<br />

<br />

<br />

<br />

t0 ω0 ∆t <br />

<br />

<br />

δt t0<br />

ω0 <br />

gt0,ω0 (t) = g0(t − t0)e iω0t<br />

t2<br />

− g0(t) = A0e 2∆t2 <br />

<br />

<br />

ω0∆t c0 <br />

<br />

gt0,∆t(t) = 1<br />

√ ∆t g0<br />

t − t0<br />

∆t<br />

t2<br />

− g0 = A0e 2∆t2 eic0t <br />

∆ω<br />

ω0 <br />

ω<br />

<br />

t<br />

ω<br />

<br />

<br />

<br />

<br />

<br />

<br />

t<br />

t


ψ <br />

<br />

<br />

|<br />

Cψ =<br />

R<br />

ψ(ω)| 2<br />

dω < +∞. <br />

ω<br />

Wf f <br />

<br />

f(t) = 1<br />

<br />

Wf(t, s) =<br />

Cψ<br />

<br />

R 2<br />

R<br />

f(τ) 1 <br />

τ − t<br />

√ ψ dτ. <br />

s s<br />

Wf(τ, σ) 1<br />

√ σ ψ<br />

<br />

t − τ<br />

σ<br />

dτ dσ<br />

. <br />

σ2 <br />

<br />

f L2(R) = 1<br />

Cψ<br />

<br />

R 2<br />

|Wf(t, s)| 2 dt ds<br />

. <br />

s2 Wf <br />

f<br />

<br />

<br />

<br />

<br />

<br />

<br />

gt0,∆t(t) = 1<br />

√ ∆t g<br />

t − t0<br />

∆t<br />

∆t <br />

∆t <br />

∆t = b j<br />

t0 = k∆t<br />

2 −j <br />

g0(2 j (t−2 −j k)) = g0(2 j t−k) j <br />

k L2 <br />

(ψjk)j,k∈Z ψjk(t) = 2 j/2 ψ(2 j t−k)


ψ <br />

m0 m1 <br />

k ↦→ m0[k] k ∈ Z<br />

k ↦→ m1[k] k ∈ Z<br />

ω ↦→ m0(ω) ω ↦→ m1(ω) 2π <br />

φ ψ L2(R) <br />

+∞ <br />

φ(ω) = m0(<br />

k=1<br />

ω<br />

) <br />

2k ψ(ω) = m1( ω<br />

2 ) φ( ω<br />

2<br />

) <br />

(ψjk) <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

L2(R) (Vj)j∈Z <br />

<br />

j∈Z<br />

Vj =<br />

<br />

<br />

akφjk : ak ∈ R<br />

k∈Z<br />

<br />

Vj ⊂ Vj+1<br />

<br />

Vj = {0} <br />

<br />

Vj = L2(R). <br />

j∈Z


Vj <br />

(φjk)k∈Z <br />

<br />

φjk <br />

φ φjk <br />

L2 <br />

ℓ2(Z) → L2(R)<br />

(ak)k∈Z ↦→ <br />

k∈Z<br />

akφ0k<br />

φ <br />

<br />

(φjk)k∈Z Vj+1<br />

Vj <br />

<br />

φ ∈ V1 <br />

<br />

φ(t) = 2 <br />

m0[k]φ(2t − k). <br />

k∈Z<br />

m0<br />

<br />

L2(R) <br />

φ φ 0<br />

Vj <br />

Vj+1 W0 <br />

V0 ⊕ W0 = V1. <br />

W0 ψ<br />

<br />

W0 =<br />

<br />

t ↦→ <br />

<br />

dkψ(t − k) : dk ∈ Z<br />

k∈Z<br />

ψ V1 <br />

. <br />

ψ(t) = <br />

m1[k]φ(2t − k), <br />

k∈Z<br />

m1


j <br />

<br />

Vj ⊕ Wj = Vj+1, <br />

Vj ⊕ Wj ⊕ · · · ⊕ Wj ′ −1 = Wj ′ j < j′ . <br />

j ′ +∞ j −∞ <br />

<br />

L2(R) = Vj ⊕<br />

L2(R) =<br />

+∞<br />

j=−∞<br />

+∞<br />

j ′ =j<br />

Wj<br />

Wj ′ ∀j ∈ Z <br />

<br />

<br />

<br />

Bj = {φjk : k ∈ Z} ∪ {ψj ′ k : j ′ ≥ j, k ∈ Z} <br />

B = {ψjk : j ∈ Z, k ∈ Z} <br />

<br />

<br />

<br />

f = 2<br />

f[k/2 j ] <br />

<br />

j/2<br />

f[k/2 j ]φjk<br />

k∈Z<br />

f[k/2 j ] f(k/2 j ).<br />

f Vj <br />

<br />

L < j.<br />

VL ⊕ WL ⊕ · · · ⊕ Wj−1<br />

<br />

Vj ′ Vj ′ −1 ⊕ Wj ′ −1 <br />

<br />

Vj−1 ⊕ Wj−1<br />

Vj−2 ⊕ Wj−2 ⊕ Wj−1<br />

<br />

VL ⊕ WL ⊕ WL+1 ⊕ · · · ⊕ Wj−1.


Vj+1 → Vj ⊕ Wj.<br />

<br />

ℓ2(Z) → ℓ2(Z) × ℓ2(Z)<br />

(aj+1,k)k∈Z ↦→ [(ajk)k∈Z, (djk)k∈Z].<br />

Aj Dj 2π <br />

k ↦→ ajk k ↦→ djk <br />

Aj(ω) = <br />

k∈Z<br />

Dj(ω) = <br />

k∈Z<br />

ajke −ikω<br />

djke −ikω<br />

<br />

Aj(2ω)<br />

Dj(2ω)<br />

<br />

=<br />

m0(ω) m0(ω + π)<br />

m1(ω) m1(ω + π)<br />

Aj+1(ω)<br />

Aj+1(ω + π)<br />

<br />

<br />

<br />

<br />

T (ω) =<br />

m0(ω) m0(ω + π)<br />

m1(ω) m1(ω + π)<br />

[0, 2π] [0, 2π] <br />

˜ T (ω) = T (ω) −1 2π<br />

˜m0 ˜ m1 ˜ T (ω) <br />

˜T (ω) =<br />

˜m0(ω) ˜m0(ω + π)<br />

˜m1(ω) ˜m1(ω + π)<br />

˜ φ ˜ ψ <br />

<br />

˜φ(ω) =<br />

+∞ <br />

k=1<br />

˜ψ(ω) = ˜m1<br />

˜m0<br />

ω<br />

2<br />

<br />

ω<br />

<br />

2 k<br />

˜ φ<br />

ω<br />

2<br />

<br />

<br />

<br />

<br />

<br />

˜ φ ˜ ψ j <br />

L2(R) <br />

f = <br />

k∈Z<br />

<br />

f, ˜ <br />

φjk φjk + <br />

j ′ ≥j,k∈Z<br />

<br />

f, ˜ ψj ′ <br />

k ψj ′ k


j ∈ Z f ∈ ℓ2(Z) j −∞ <br />

<br />

f = <br />

j,k∈Z<br />

<br />

f, ˜ <br />

ψjk ψjk<br />

<br />

<br />

m0 m1 ˜m0 ˜m1 <br />

<br />

{φjk : k ∈ Z} ∪ {ψjk : k ∈ Z} ↔ {φj+1,k : k ∈ Z}<br />

<br />

<br />

ajk = 2 <br />

l∈Z<br />

djk = 2 <br />

l∈Z<br />

˜m0[k]aj+1,2l−k<br />

˜m1[k]aj+1,2l−k<br />

<br />

aj+1,k = 1<br />

2<br />

<br />

m0[2l − k]ajl + m1[2l − k]djl<br />

l∈Z<br />

<br />

j = 0 j = −3<br />

N <br />

A × N<br />

A <br />

N log N<br />

<br />

<br />

ψ t ↦→ 2j/2ψ(2jt− k) <br />

j,k∈Z L2(R) φ = ˜ φ ψ = ˜ ψ <br />

<br />

ω <br />

m0 m1 <br />

|m0(ω)| 2 + |m0(ω + π)| 2 = 1 ∀ω<br />

m0(ω)m1(ω) + m0(ω + π)m1(ω + π) = 0 ∀ω<br />

|m1(ω)| 2 + |m1(ω + π)| 2 = 1 ∀ω


a−3,k<br />

d−3,k<br />

a0,k<br />

1 ↑ 2<br />

1 ↑ 2<br />

˜m0<br />

˜m1<br />

m0<br />

m1<br />

2 ↓ 1<br />

2 ↓ 1<br />

+<br />

d−2,k<br />

˜m0<br />

˜m1<br />

2 ↓ 1<br />

2 ↓ 1<br />

˜m0<br />

˜m1<br />

2 ↓ 1<br />

2 ↓ 1<br />

a−3,k<br />

d−1,k d−2,k d−3,k<br />

<br />

1 ↑ 2<br />

1 ↑ 2<br />

m0<br />

m1<br />

+<br />

d−1,k<br />

<br />

<br />

<br />

<br />

1 ↑ 2<br />

1 ↑ 2<br />

m0<br />

m1<br />

+<br />

a0,k


m0 m1 <br />

<br />

<br />

m0 m1 <br />

m1(ω) = e iω m0(ω + π)<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

L2 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

j → +∞ <br />

<br />

˜ ψ p <br />

<br />

˜ψ(t)t k dt = 0 ∀k ∈ {0, . . . , p − 1}.<br />

R<br />

˜ ψ <br />

p ω = 0 ψ <br />

p<br />

f p <br />

p I 2 −j(p+1/2)


I M <br />

|〈 ˜ ψjk, f〉| ≤ M2 −j(p+1/2) supp ˜ ψjk ⊂ I.<br />

<br />

<br />

|〈 ˜ ψjk, f〉| ≤ M2 −j(p+1/2) ψ p <br />

f r r < p<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

x = (xk)k∈Z <br />

(x2k)k∈Z (x2k+1)k∈Z xe<br />

xo <br />

h <br />

ℓ <br />

xe <br />

xo xo <br />

P <br />

h = xo − P (xe) <br />

xe h xo <br />

xo = d + P (xe) <br />

P xo,k hk h <br />

xo P <br />

<br />

(xe, d) (xe, xo)


xe <br />

h <br />

ℓ = xe + U(h) <br />

<br />

xe = ℓ + U(h) <br />

<br />

<br />

<br />

<br />

(x)<br />

(x2k)<br />

(x2k+1)<br />

P<br />

−<br />

+ (ℓk)<br />

U<br />

<br />

(hk)<br />

(ℓk)<br />

(hk)<br />

<br />

−<br />

U<br />

P<br />

+<br />

(x2k)<br />

(x2k+1)<br />

<br />

<br />

<br />

<br />

<br />

<br />

T = GL1U1 · · · LnUn, <br />

(Li)1≤i≤n (Ui)1≤i≤n <br />

G <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

(xk)


xo,k <br />

xe,k <br />

<br />

hk = x2k+1 − x2k<br />

<br />

<br />

<br />

<br />

<br />

<br />

ℓk = x2k + 1<br />

2 hk<br />

hk = x2k+1 − x2k<br />

<br />

<br />

ℓk = 1<br />

2 (x2k + x2k+1) <br />

d <br />

<br />

<br />

<br />

hk = x2k+1 − 1<br />

2 (x2k + x 2(k+1)) <br />

<br />

<br />

ℓk = x2k + 1<br />

4 (hk + hk+1) <br />

hk = x2k+1 − 1<br />

2 (x2k + x 2(k+1)) <br />

ℓk = 3<br />

4 x2k + 1<br />

4 (x2k−1 + x2k+1) − 1<br />

8 (x2k−2 + x2k+2)


yk = x2k+1 − α(x2k + x 2(k+1)) α −1.58613<br />

zk = x2k + β(x 2(k−1)+1 + x2k+1) β −0.05298<br />

h ′ k = yk − γ(zk + zk+1) γ 0.88291<br />

ℓ ′ k = zk + δ(yk−1 + yk) δ 0.44351<br />

hk = ηh ′ k<br />

ℓk = 1<br />

η ℓ′ k<br />

η 1.14960<br />

<br />

<br />

(xk)<br />

(x2k)<br />

(x2k+1)<br />

P U P U<br />

−<br />

+<br />

−<br />

+ (ℓk)<br />

<br />

<br />

<br />

<br />

<br />

y z <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

(hk)


hk = x2k+1 − ⌊ 1<br />

2 (x2k + x 2(k+1))⌋ <br />

ℓk = x2k + ⌊ 1<br />

4 (hk + hk+1)⌋


x2k+1 <br />

<br />

<br />

x2k+1 <br />

P (x2k+1) = x2k + x2k<br />

,<br />

2<br />

<br />

<br />

x2k+1 <br />

P (x2k+1) = x2k + x2k − x2(k−1) .<br />

2<br />

<br />

<br />

<br />

<br />

<br />

Gi <br />

i ≤ j ⇒ Gi ⊂ Gj<br />

Gi Gi+1 \ Gi <br />

<br />

<br />

<br />

<br />

<br />

φ ψ <br />

φℓℓ = φ ⊗ φ <br />

ψℓh = φ ⊗ ψ <br />

ψhℓ = ψ ⊗ φ <br />

ψhh = ψ ⊗ ψ. <br />

φℓℓ ψℓh ψhℓ ψhh


ma mb <br />

(f ∗ ma) ∗ mb = (f ∗ mb) ∗ ma = f ∗ (ma ∗ mb).<br />

<br />

<br />

d ψhℓ d ψℓh d ψhh<br />

<br />

P U <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

xe = (xi,j)i,j∈Z,i+j≡0 mod 2 <br />

xo = (xi,j)i,j∈Z,i+j≡1 mod 2, <br />

<br />

j −2 <br />

j − 1


j = 0 j = −1 j = −2 j = −3 j = −4<br />

j <br />

<br />

s<br />

j = 0<br />

j = −1<br />

s<br />

√ 2<br />

s<br />

2<br />

j = −2<br />

j = −3<br />

s<br />

2 √ 2<br />

s<br />

4<br />

j = −4<br />

j<br />

2 <br />

<br />

<br />

<br />

hi,j = xi,j − med(xi−1,j, xi+1,j, xi,j−1, xi,j+1) i + j ≡ 1 mod 2 <br />

ℓi,j = xi,j + 1<br />

2 med(hi−1,j, hi+1,j, hi,j−1, hi,j+1) i + j ≡ 0 mod 2 <br />

<br />

<br />

<br />

hi,j = xi,j − max(xi−1,j, xi+1,j, xi,j−1, xi,j+1) i + j ≡ 1 mod 2 <br />

ℓi,j = xi,j + max(hi−1,j, hi+1,j, hi,j−1, hi,j+1) i + j ≡ 0 mod 2


s t ⇐⇒ 0 ≤ s ≤ t t ≤ s ≤ 0 ∀s, t ∈ R, <br />

<br />

K <br />

K = med(0, min K, max K).<br />

<br />

hi,j = xi,j − {xi−1,j, xi+1,j, xi,j−1, xi,j+1} i + j ≡ 1 mod 2 <br />

ℓi,j = xi,j + 1<br />

2 {hi−1,j, hi+1,j, hi,j−1, hi,j+1} i + j ≡ 0 mod 2 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

h = xo − Padap(xe) <br />

ℓ = xe + U(h) <br />

Padap(xe)(n) = xe(k)Fn(n − k) Fn (xe)<br />

n


(x)<br />

(x2k)<br />

(x2k+1)<br />

Padap<br />

−<br />

+ (ℓk)<br />

<br />

<br />

<br />

<br />

h xe <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

˜xe <br />

xe <br />

<br />

|Padap( ˜xe) − Padap(xe)| 0.<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

U<br />

(hk)


(xk)<br />

(x2k)<br />

(x2k+1)<br />

+ (ℓk)<br />

U<br />

Padap<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

−<br />

(hk)


I(t, x) <br />

t x = (x1, x2) <br />

t x I <br />

<br />

<br />

X1(t), X2(t), X3(t) (x1(t), x2(t))<br />

<br />

v = dx(t)<br />

<br />

dt<br />

x(t) <br />

I(t, x(t))


d<br />

∂I<br />

I(t, x(t)) = + v · ∇I <br />

dt ∂t<br />

<br />

<br />

<br />

d<br />

I(t, x(t)) = 0 <br />

dt<br />

<br />

∂I<br />

+ v · ∇I = 0 <br />

∂t<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

∇I<br />

?<br />

?<br />

?<br />

?<br />

?<br />

?<br />

?<br />

t = 0 t = 1<br />

<br />

<br />

∇I t = 0 t = 1


Image<br />

échelle s<br />

vδt<br />

Image<br />

Image<br />

échelle s<br />

? Image<br />

?<br />

vδt<br />

|v| ≪ s |v| > s<br />

<br />

s <br />

v |v| s <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

εx0,t = I(t, x0 + v(x0)) − I(t + 1, x0)


G <br />

<br />

W0 <br />

v0 = argmin<br />

v∈G<br />

<br />

W0<br />

(I(t, x + v) − I(t + 1, x)) 2 d 2 x <br />

<br />

<br />

<br />

<br />

I Ĩ <br />

<br />

<br />

( Ĩ(t, x + v) − I(t + 1, x))2 d 2 x <br />

v0 = argmin<br />

v∈G W0<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

W1 <br />

W0 <br />

<br />

v0 = argmin<br />

v∈G<br />

<br />

W1<br />

( Ĩ(t, x + v) − I(t + 1, x))2 d 2 x <br />

W0 ⊂ W1


2 k <br />

2 k−1 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

M(v) =<br />

∂I<br />

∂t + v · ∇I 2 d 2 x. <br />

<br />

<br />

<br />

R(v) =<br />

v 2 Hd 2 x <br />

·H <br />

<br />

<br />

<br />

v <br />

v = argmin M(v) + λR(v).


λ<br />

<br />

<br />

<br />

<br />

<br />

F <br />

<br />

<br />

I(t, x) = I(0, x) ∗ δ(x − vt). <br />

<br />

Î(ω, ξ) ∝ Î(0, ξ)δ(tvξ − ω). <br />

<br />

<br />

<br />

t<br />

x<br />

(x, t) (ωx, τ)<br />

<br />

<br />

x t <br />

<br />

<br />

(vx, 1) <br />

vx<br />

<br />

<br />

<br />

<br />

1<br />

τ<br />

vx<br />

ωx


M(v) = <br />

x∈W (x0)<br />

ρ(x) 2 ∇I · v + ∂I 2 <br />

∂t<br />

W (x0) <br />

ρ <br />

<br />

<br />

N W (x0) t<br />

t AR 2 Av = t AR 2 b <br />

A = t ∇I(x1), . . . , ∇I(xN) <br />

<br />

R = diag(ρ(x1), . . . , ρ(xN)) <br />

b = − t ( ∂I<br />

∂t (x1), . . . , ∂I<br />

∂t (xN)) <br />

t AR 2 A 2 × 2 <br />

<br />

t AR 2 A <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

v j Br r <br />

2 j Z 2


vj−1 <br />

r<br />

B r 2 2 <br />

v {j,j−1} = v j + v j−1 <br />

{v1 } + B r<br />

<br />

2<br />

2j−1Z2 <br />

<br />

<br />

<br />

<br />

<br />

(ψ n )1≤n≤N <br />

x0 <br />

<br />

∇I · v + ∂I<br />

∂t<br />

<br />

<br />

∂I<br />

∂x1<br />

v1, ψ n x0<br />

<br />

ψ n k (x − x0)d 2 x = 0 ∀n = 1 . . . N <br />

<br />

〈f, g〉 =<br />

f(x)g(x)d 2 x<br />

ψ n x0 (x) = ψn (x − x0),<br />

<br />

∂I<br />

+ v2, ψ<br />

∂x2<br />

n <br />

∂I<br />

x0 +<br />

∂t , ψn <br />

x0 = 0 ∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

I, ∂ψn x0<br />

∂x1<br />

<br />

v1(x0) +<br />

I, ∂ψn x0<br />

∂x2<br />

<br />

v2(x0) = ∂<br />

<br />

I, ψ<br />

∂t<br />

n <br />

x0<br />

∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

|Av − b| A


vj−1 <br />

r<br />

B r 2 2 <br />

v {j,j−1} = v j + v j−1 <br />

{v1 } + B r<br />

<br />

2<br />

2j−1Z2 <br />

<br />

<br />

<br />

<br />

<br />

(ψ n )1≤n≤N <br />

x0 <br />

<br />

∇I · v + ∂I<br />

∂t<br />

<br />

<br />

∂I<br />

∂x1<br />

v1, ψ n x0<br />

<br />

ψ n k (x − x0)d 2 x = 0 ∀n = 1 . . . N <br />

<br />

〈f, g〉 =<br />

f(x)g(x)d 2 x<br />

ψ n x0 (x) = ψn (x − x0),<br />

<br />

∂I<br />

+ v2, ψ<br />

∂x2<br />

n <br />

∂I<br />

x0 +<br />

∂t , ψn <br />

x0 = 0 ∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

I, ∂ψn x0<br />

∂x1<br />

<br />

v1(x0) +<br />

I, ∂ψn x0<br />

∂x2<br />

<br />

v2(x0) = ∂<br />

<br />

I, ψ<br />

∂t<br />

n <br />

x0<br />

∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

|Av − b| A


vj−1 <br />

r<br />

B r 2 2 <br />

v {j,j−1} = v j + v j−1 <br />

{v1 } + B r<br />

<br />

2<br />

2j−1Z2 <br />

<br />

<br />

<br />

<br />

<br />

(ψ n )1≤n≤N <br />

x0 <br />

<br />

∇I · v + ∂I<br />

∂t<br />

<br />

<br />

∂I<br />

∂x1<br />

v1, ψ n x0<br />

<br />

ψ n k (x − x0)d 2 x = 0 ∀n = 1 . . . N <br />

<br />

〈f, g〉 =<br />

f(x)g(x)d 2 x<br />

ψ n x0 (x) = ψn (x − x0),<br />

<br />

∂I<br />

+ v2, ψ<br />

∂x2<br />

n <br />

∂I<br />

x0 +<br />

∂t , ψn <br />

x0 = 0 ∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

I, ∂ψn x0<br />

∂x1<br />

<br />

v1(x0) +<br />

I, ∂ψn x0<br />

∂x2<br />

<br />

v2(x0) = ∂<br />

<br />

I, ψ<br />

∂t<br />

n <br />

x0<br />

∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

|Av − b| A


vj−1 <br />

r<br />

B r 2 2 <br />

v {j,j−1} = v j + v j−1 <br />

{v1 } + B r<br />

<br />

2<br />

2j−1Z2 <br />

<br />

<br />

<br />

<br />

<br />

(ψ n )1≤n≤N <br />

x0 <br />

<br />

∇I · v + ∂I<br />

∂t<br />

<br />

<br />

∂I<br />

∂x1<br />

v1, ψ n x0<br />

<br />

ψ n k (x − x0)d 2 x = 0 ∀n = 1 . . . N <br />

<br />

〈f, g〉 =<br />

f(x)g(x)d 2 x<br />

ψ n x0 (x) = ψn (x − x0),<br />

<br />

∂I<br />

+ v2, ψ<br />

∂x2<br />

n <br />

∂I<br />

x0 +<br />

∂t , ψn <br />

x0 = 0 ∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

I, ∂ψn x0<br />

∂x1<br />

<br />

v1(x0) +<br />

I, ∂ψn x0<br />

∂x2<br />

<br />

v2(x0) = ∂<br />

<br />

I, ψ<br />

∂t<br />

n <br />

x0<br />

∀n = 1 . . . N <br />

<br />

<br />

<br />

<br />

<br />

|Av − b| A


N (F [k])0≤k


¯ J <br />

[yn, yn+1[ <br />

n <br />

<br />

[yn, yn+1[ <br />

¯ F [k] F [k] xn ∈ [yn, yn+1[ <br />

<br />

F [k] ∈ [yn, yn+1[ ⇒ ¯ F [k] = xn.<br />

<br />

∆ <br />

T θ = T/∆ <br />

θ = 2<br />

n <br />

F [k] <br />

<br />

F [k] ¯ F [k] <br />

<br />

<br />

<br />

<br />

<br />

I(x, y) <br />

Ī(x, y) mathcalI <br />

<br />

(I, Ī) = 10 log 10<br />

<br />

(x,y)∈I<br />

<br />

(x,y)∈I I(x, y)2<br />

(Ī(x, y) − I(x, y))2 .<br />

<br />

<br />

<br />

<br />

<br />

<br />

Vmax <br />

<br />

(I, Ī) = 10 log 10<br />

<br />

(x,y)∈I<br />

V 2 max card I<br />

.<br />

(Ī(x, y) − I(x, y))2


X <br />

p X <br />

H(X) <br />

H(X) = − <br />

p(x)log2p(x)<br />

x∈X<br />

<br />

<br />

<br />

p


ε Eε <br />

Eε = {v |Ax − b| < ε}


ε Eε <br />

Eε = {v |Ax − b| < ε}


ε Eε <br />

Eε = {v |Ax − b| < ε}


E N <br />

(ei)1≤i≤N <br />

I = I1 × I2 × · · · × IN,<br />

Ii R <br />

(Gr )rmin≤r≤rmax <br />

<br />

G r = 2 r Z ∩ I.<br />

p Gr Ar p m × n <br />

br p n vr p m Ar br vr Gr <br />

<br />

<br />

(A r v r )p = A r pv r p.<br />

ar <br />

<br />

r<br />

a = a r p.<br />

p∈G r<br />

<br />

A r pv r p = b r p.<br />

<br />

A r = F(v r+1 , D)<br />

b r = G(v r+1 , D)<br />

D It It+1 <br />

<br />

v rmax+1 = 0.<br />

<br />

<br />

<br />

v rmin <br />

v r rmax rmin <br />

p r<br />

argmin v r p |A r pv r p − b r p| 2 <br />

t A r pA r pv r p − t A r pb r p = 0,


t A r pA r p <br />

<br />

A r <br />

A r (v r ) = |A r v r − b r | 2 .<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

argmin<br />

vr A r (v r ) + λB r (v r ) <br />

Br <br />

<br />

<br />

<br />

B r (v r <br />

) = |∇v r i | 2 (v r i )p = v r p,i<br />

1≤i≤m<br />

<br />

λ <br />

A <br />

B <br />

<br />

<br />

B r M r H r c r <br />

(B r v r )p,i,j = ∂vr q,j<br />

(p) B<br />

∂qi<br />

r ∇<br />

M r = t A r A r<br />

H r = t B r B r H r <br />

c r = t A r b r<br />

A r B r M r H r <br />

(A r ) i′ 1 ...i′ N ,j′<br />

i1...iN ,j = 0 ∃ k ik = i ′ k<br />

(B r ) i′ 1 ...i′ N ,j′<br />

i1...iN ,j,l = 0 j = j′<br />

(M r ) i′ 1 ...i′ N ,j′<br />

i1...iN ,j = 0 ∃ k ik = i ′ k<br />

(H r ) i′ 1 ...i′ N ,j′<br />

i1...iN ,j<br />

= 0 j = j′


argmin<br />

v r<br />

|A r v r − b r | 2 + λ tv r H r v r ,<br />

argmin<br />

vr tvr r r r r t r r<br />

((M + λH )v − 2c ) + b b .<br />

λ <br />

λ <br />

<br />

λ <br />

<br />

tr M r<br />

λ0 = . <br />

tr Hr λ0 <br />

<br />

<br />

<br />

<br />

(H r ) i′ 1 ...i′ N ,j<br />

i1...iN ,j = (Hr ) (i′ 1 −i1)...(i ′ N −iN ),1<br />

0...0,1 , <br />

λ0 =<br />

=<br />

<br />

tr M r<br />

i1,...iN ,j (Hr ) i1...iN ,j<br />

i1...iN ,j<br />

tr M r<br />

mN card Gr(H r ) 0...0,1<br />

0...0,1<br />

<br />

<br />

λ0 <br />

λ/λ0 <br />

<br />

<br />

<br />

<br />

<br />

∇(A + λB) = 2[(M r + λH r )v r − c r ]<br />

<br />

(v r,0<br />

p ) <br />

A


nreg <br />

v r,i+1 = (1 − ɛ(M r + λH r ))v r,i + ɛc r <br />

ɛ <br />

<br />

0 < ɛ <<br />

2<br />

max vp(M r + λHr . <br />

)<br />

<br />

<br />

H r κ = 4N <br />

0 < ɛ <<br />

2<br />

max vp(M r ) + λκ<br />

<br />

ɛ <br />

<br />

λ/λ0 nreg <br />

<br />

<br />

<br />

<br />

M r p <br />

<br />

∂vl,k<br />

∂x = vl+1,k − vl,k.<br />

H r <br />

∂ 2 vl,k<br />

∂x 2 = vl+1,k − 2vl,k + vl−1,k<br />

<br />

<br />

l,k (M r l,k,xx + M r l,k,yy )<br />

8N r xN r y<br />

2<br />

ɛ <<br />

max vp M r + 8λ<br />

λ0 =


Nx Ny <br />

<br />

v r,i+1<br />

l,k,x<br />

v r,i+1<br />

l,k,y<br />

= vr,i<br />

l,k,x − ɛ(M r l,k,xx vr,i<br />

l,k,x + M r l,k,xy vr,i<br />

l,k,y )<br />

− ɛλ(v r,i<br />

l−1,k,x − 2vr,i<br />

l,k,x + vr,i<br />

l+1,k,x )<br />

+ ɛc r l,k,x<br />

= vr,i<br />

l,k,y − ɛ(M r l,k,yx vr,i<br />

l,k,x + M r l,k,yy vr,i<br />

l,k,y )<br />

− ɛλ(v r,i<br />

l,k−1,y − 2vr,i<br />

l,k,y + vr,i<br />

l,k+1,y )<br />

+ ɛc r l,k,y<br />

<br />

<br />

<br />

<br />

<br />

Hv <br />

<br />

Herr <br />

<br />

<br />

S <br />

Hv Herr <br />

Htotal <br />

<br />

Htotal = 0.25Hv + Herr.<br />

<br />

λ/λ0 nreg<br />

v It It−1<br />

v ˜v<br />

∆v ˜v Q1(˜v)<br />

Hv<br />

Qv(˜v) ¯˜v<br />

¯˜v ¯v<br />

I pred<br />

t<br />

It−1 ¯v<br />

Et I pred<br />

t It<br />

Et ˜ Et<br />

∆err ˜ Et Qerr( ˜ Et)<br />

Herr<br />

Qerr( ˜ Et) ¯˜ Et


¯˜ Et Ēt<br />

Īt I pred<br />

t<br />

Īt It<br />

Entrées<br />

It−1 ME v CD1 ¯v<br />

It<br />

−<br />

Et<br />

Īt<br />

I pred<br />

t<br />

+<br />

MC<br />

CD2<br />

Ēt<br />

Ēt<br />

Ev<br />

Eerr<br />

PSNR<br />

<br />

<br />

Qi ∆i <br />

<br />

<br />

× × <br />

<br />

ɛ <br />

<br />

<br />

<br />

<br />

nspl = 3 <br />

<br />

nspl = 1 <br />

<br />

Sorties


46<br />

45<br />

44<br />

43<br />

42<br />

A<br />

B<br />

A<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

B<br />

0.10 0.15 0.20 0.25<br />

A<br />

<br />

B<br />

A<br />

40<br />

39<br />

38<br />

37<br />

36<br />

A<br />

B<br />

A<br />

0.2 0.3 0.4 0.5 0.6<br />

B<br />

A<br />

<br />

= f(Htotal) λ/λ0 = 0.25 ∆v = 1 nreg = 1000<br />

∆err nspl = 1 <br />

nspl = 3<br />

A


λ/λ0<br />

λ/λ0 <br />

λ/λ0 <br />

<br />

<br />

<br />

<br />

47<br />

46<br />

45<br />

44<br />

43<br />

42<br />

A<br />

D<br />

C<br />

B<br />

A<br />

D<br />

C<br />

0.05 0.10 0.15 0.20<br />

B<br />

A<br />

<br />

D<br />

C<br />

B<br />

A<br />

D<br />

38<br />

36<br />

34<br />

32<br />

30<br />

A<br />

D<br />

C<br />

B<br />

A<br />

D<br />

C<br />

B<br />

0.6 0.8 1.0 1.2 1.4 1.6 1.8<br />

<br />

= f(Htotal) ∆v = 4 nspl = 1 ∆err <br />

nreg = 100 λ/λ0 = 0.25 <br />

nreg = 100 λ/λ0 = 0.125 nreg = 1000 λ/λ0 = 0.25<br />

nreg = 1000 λ/λ0 = 0.125<br />

<br />

<br />

<br />

<br />

<br />

λ/λ0 = 0.25 λ/λ0 = 0.125<br />

<br />

<br />

<br />

<br />

A<br />

D<br />

C<br />

B<br />

A


47<br />

46<br />

45<br />

44<br />

43<br />

42<br />

A<br />

C<br />

B<br />

A<br />

C<br />

0.05 0.10 0.15 0.20<br />

B<br />

<br />

A<br />

C<br />

B<br />

A<br />

38<br />

36<br />

34<br />

32<br />

30<br />

A<br />

C<br />

B<br />

A<br />

C<br />

0.4 0.6 0.8 1.0 1.2 1.4 1.6<br />

B<br />

A<br />

<br />

= f(Htotal) λ/λ0 = 0.25 nspl = 1 ∆v = 4<br />

∆err nreg = 0 <br />

nreg = 100 nreg = 1000<br />

38<br />

36<br />

34<br />

32<br />

A<br />

C<br />

B<br />

A<br />

C<br />

B<br />

0.4 0.6 0.8 1.0 1.2<br />

A<br />

<br />

C<br />

B<br />

A<br />

38<br />

36<br />

34<br />

32<br />

30<br />

A<br />

C<br />

B<br />

A<br />

C<br />

B<br />

0.6 0.8 1.0 1.2 1.4 1.6 1.8<br />

<br />

= f(Htotal) λ/λ0 = 0.125 nspl = 1 ∆v = 4<br />

∆err nreg = 0 <br />

nreg = 100 nreg = 1000<br />

A<br />

C<br />

C<br />

B<br />

B<br />

A<br />

A


Herr Hv <br />

Htotal <br />

<br />

<br />

∆v = 2 <br />

<br />

<br />

<br />

<br />

<br />

<br />

1.0<br />

0.8<br />

0.6<br />

0.4<br />

A<br />

C<br />

B<br />

C<br />

A<br />

B<br />

A<br />

0.05 0.10 0.15 0.20<br />

C<br />

<br />

B<br />

A<br />

0.58<br />

0.56<br />

0.54<br />

0.52<br />

0.50<br />

0.48<br />

0.46<br />

A<br />

C<br />

B<br />

C<br />

A<br />

0.05 0.10 0.15<br />

B<br />

<br />

Herr = f(Hv) λ/λ0 = 0.25 nspl = 3 ∆err = 8<br />

∆v nreg = 0 <br />

nreg = 100 nreg = 1000<br />

<br />

<br />

nreg λ/λ0 <br />

<br />

<br />

<br />

<br />

<br />

C<br />

A


1.2<br />

1.1<br />

1.0<br />

0.9<br />

0.8<br />

0.7<br />

0.6<br />

A<br />

C<br />

B<br />

C<br />

A<br />

B<br />

A<br />

C<br />

0.05 0.10 0.15 0.20<br />

<br />

B<br />

A<br />

0.68<br />

0.66<br />

0.64<br />

0.62<br />

0.60<br />

0.58<br />

0.56<br />

0.54<br />

A<br />

C<br />

B<br />

A<br />

C<br />

0.05 0.10 0.15<br />

<br />

Htotal = f(Hv) λ/λ0 = 0.25 nspl = 3 ∆err = 8<br />

∆v nreg = 0 <br />

nreg = 100 nreg = 1000<br />

47<br />

46<br />

45<br />

44<br />

43<br />

42<br />

A<br />

A<br />

B<br />

C<br />

A<br />

B<br />

0.05 0.10 0.15 0.20 0.25<br />

<br />

C<br />

A<br />

B<br />

C<br />

40<br />

38<br />

36<br />

34<br />

32<br />

A<br />

C<br />

B<br />

A<br />

0.2 0.4 0.6 0.8 1.0<br />

C<br />

B<br />

<br />

= f(Htotal) λ/λ0 = 0.25 nreg = 1000 nspl = 1<br />

∆err ∆v = 4 <br />

∆v = 2 ∆v = 1<br />

B<br />

A<br />

A<br />

A


λ/λ0 = 0.125


λ/λ0 = 0.125


λ/λ0 = 0.125


φij xij xij <br />

<br />

<br />

v(x) = <br />

i,j<br />

vij1<br />

vij2<br />

<br />

φij(x) <br />

<br />

<br />

<br />

φij


φ00 <br />

<br />

xij φ00<br />

<br />

<br />

<br />

<br />

<br />

<br />

[v] opt = argmin<br />

[v]<br />

<br />

It+1(x + v(x)) − It(x) A<br />

x<br />

<br />

·A ·1 ·2 <br />

<br />

<br />

[v] opt = argmin<br />

[v]<br />

⎛<br />

<br />

x<br />

⎝It+1<br />

<br />

x + <br />

<br />

vij1<br />

φij(x)<br />

vij2<br />

i,j<br />

⎞<br />

<br />

− It(x) ⎠<br />

2


v<br />

x <br />

<br />

<br />

<br />

<br />

<br />

<br />

χ 2 ([vijl]) <br />

<br />

χ 2<br />

<br />

[v]<br />

= <br />

⎛<br />

x<br />

⎝It+1<br />

<br />

x + <br />

<br />

vij1<br />

φij(x)<br />

vij2<br />

ij<br />

⎞<br />

<br />

− It(x) ⎠<br />

2<br />

. <br />

<br />

<br />

<br />

<br />

<br />

χ 2 ([vijl]) <br />

<br />

χ 2 ([v]) = γ − ∇χ 2 · [v] + 1<br />

2 [v] · Hχ2 · [v] <br />

<br />

[v] <br />

[v] = [v] + D −1 (−∇χ 2 ([v] )). <br />

<br />

<br />

[v] = [v] − · ∇χ 2 ([v] ), <br />

<br />

<br />

<br />

<br />

x ′ = x + v(x). <br />

x ′ <br />

<br />

χ 2 (v) = <br />

It+1(x ′ ) − It(x) .<br />

x


vijl <br />

∂χ 2<br />

∂vijl<br />

= 2 <br />

<br />

∂ 2 χ 2<br />

∂vijl∂vi ′ j ′ l ′<br />

= 2 <br />

x<br />

x<br />

(φijφi ′ j ′)(x)<br />

<br />

∂2It+1 φij(x) ∂It+1<br />

(x<br />

∂xl<br />

′ <br />

) It+1(x ′ <br />

) − It(x)<br />

∂xl∂xl ′<br />

(x ′ )<br />

<br />

It+1(x ′ <br />

) − It(x)<br />

<br />

∂It+1<br />

+<br />

∂xl<br />

∂It+1<br />

∂xl ′<br />

<br />

<br />

(x ′ <br />

)<br />

<br />

<br />

<br />

∂ 2 χ 2<br />

∂vijl∂vi ′ j ′ l ′<br />

= 2 <br />

<br />

x<br />

αkk<br />

(φijφi ′ j ′)(x)<br />

<br />

∂It+1 ∂It+1<br />

βk = ∂χ2<br />

∂vijl<br />

1<br />

′ =<br />

2<br />

∂ 2 χ 2<br />

∂vijl∂vi ′ j ′ l ′<br />

∂xl<br />

,<br />

∂xl ′<br />

<br />

(x ′ ) <br />

k k ′ ijl i ′ j ′ l ′<br />

<br />

<br />

αkk ′δvk ′ = βk, <br />

δvk = × βk<br />

<br />

<br />

<br />

1/αkk βk <br />

λ <br />

k <br />

δvk = 1<br />

λαkk<br />

βk λαkkδvk = βk. <br />

αkk <br />

αkk


[α ′ ] [α] <br />

<br />

α ′ kk = (1 + λ)αkk<br />

α ′ kk ′ = αkk ′ k = k′ .<br />

<br />

′<br />

α kk ′δvk ′ = βk<br />

<br />

λ <br />

λ [α ′ ]<br />

<br />

λ [α] [α ′ ] <br />

<br />

<br />

<br />

<br />

[v]<br />

λ<br />

<br />

δ [v] ← (Hχ 2 ([v]) + λId) −1 ∇χ 2 ([v])<br />

χ 2 ([v] + δ [v]) < χ 2 ([v]) <br />

λ ← λ/10<br />

[v] ← [v] + δ [v]<br />

<br />

λ ← 10λ<br />

<br />

χ 2 < ε1 δχ 2 < ε2<br />

<br />

t + 1 <br />

<br />

<br />

<br />

<br />

<br />

<br />

mij


q <br />

<br />

It+1(q) = <br />

mijψi(qx)ψj(qy) = <br />

i,j<br />

∇It+1(q) = <br />

i,j<br />

i,j<br />

mijΨi,j(q) <br />

mij∇Ψi,j(q) <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

s <br />

E(s) = 0.


s Sp s <br />

Sp s = E(|s|).<br />

<br />

<br />

<br />

<br />

<br />

AutoCorr(s)(x) = E(<br />

s(t)s(t + x)dt)<br />

<br />

<br />

<br />

<br />

AutoCorr(s)(x) = E(<br />

Sp s =<br />

s(ω)e ixω s(ω)dω)<br />

= F −1 E(ω ↦→ |ω| 2 ) (x)<br />

<br />

AutoCorr s<br />

f <br />

Sp s ⋆ f = E(|s f| 2 )<br />

= | f| 2 E(|s| 2 )<br />

= Sp f Sp s<br />

f s f<br />

<br />

<br />

T <br />

δT (x) = <br />

δ(x − kT ).<br />

<br />

k∈Z<br />

FδT = 1<br />

T<br />

δ 2π .<br />

T<br />

<br />

T <br />

<br />

Sp sδT = E(|s ⋆ 1<br />

T<br />

= E(ω ↦→ | 1<br />

T<br />

≤ <br />

k∈Z<br />

δ 2π |<br />

T<br />

2 )<br />

<br />

k∈Z<br />

s(ω + 2kπ<br />

T )|)<br />

E(ω ↦→ | 1 2kπ<br />

s(ω +<br />

T T )|)


[− π<br />

T<br />

π , T ] <br />

<br />

2π<br />

T <br />

<br />

T <br />

2T <br />

<br />

R+ <br />

<br />

hT <br />

[− π π<br />

T , T ] s h <br />

Sp(s ⋆ hT ) = Sp s Sp hT .<br />

h2T hT <br />

<br />

h2T (ω) = hT (2ω)<br />

h2T <br />

<br />

i <br />

Sp(s ⋆ hT ⋆ h2T ) = Sp f Sp hT Sp h2T .<br />

ki(ω) =<br />

=<br />

i<br />

hT (2 j ω)<br />

j=0<br />

i<br />

j=0<br />

h2 j T (ω).<br />

| ki(ω)| 2 ≪ 1 ω ∈ [π/2, π], <br />

i <br />

<br />

<br />

<br />

<br />

(1/4, 1/2, 1/4) <br />

<br />

Sp(h)(ω) = cos 4 ω<br />

2 .


cos4 ω<br />

2 cos4 ω <br />

[π/4, π/2] <br />

<br />

<br />

<br />

− π<br />

2<br />

− π<br />

2<br />

0<br />

b k0(ω)<br />

π<br />

2<br />

b k0(ω) <br />

0<br />

ω<br />

π<br />

2 ω<br />

− π<br />

2<br />

− π<br />

2<br />

0<br />

b k1(ω)<br />

b k1(ω) <br />

<br />

<br />

G = (1/4, 1/2, 1/4) <br />

d <br />

u <br />

<br />

<br />

It <br />

J 0 t = G ⋆ It<br />

K (J i t )1≤i≤K <br />

<br />

J i+1<br />

t<br />

= d(G ⋆ J i t ).<br />

<br />

i i − 1<br />

<br />

0<br />

π<br />

2<br />

ω<br />

π<br />

2 ω


φ(x) = 1 [−1,1](|1 − x1|)1 [−1,1](|1 − x2|),<br />

φij <br />

n × n <br />

<br />

φ n ij(x) = φ(i − x1 x2<br />

, j −<br />

n n )<br />

<br />

<br />

<br />

<br />

u <br />

φ2 00<br />

<br />

<br />

J 0 t ← G ⋆ It)<br />

i 1 K <br />

J i t ← d(u(G) ⋆ J i−1<br />

t )<br />

<br />

vK+1 ← 0<br />

i K <br />

w i ← φ 2 00 ⋆ u(vi+1 )<br />

v i ← J i t J i t+1 w i<br />

<br />

w ← v 0<br />

vopt ← It It+1 w


4 × 4 8 × 8 <br />

<br />

<br />

<br />

4 × 4 <br />

8 × 8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

4 × 4 <br />

<br />

<br />

<br />

<br />

<br />

<br />

8×8


40<br />

38<br />

36<br />

34<br />

42<br />

41<br />

40<br />

39<br />

38<br />

37<br />

36<br />

A<br />

A<br />

D<br />

C<br />

A<br />

D<br />

D<br />

C<br />

B<br />

A<br />

C<br />

A<br />

B<br />

D<br />

A<br />

C C<br />

A<br />

B<br />

5 10 15<br />

<br />

A<br />

B<br />

D<br />

D<br />

C C<br />

A<br />

A<br />

D<br />

B<br />

A<br />

A<br />

A<br />

B D<br />

5 10 15<br />

<br />

B<br />

D<br />

B<br />

A<br />

C<br />

C<br />

D<br />

A<br />

B<br />

D<br />

B<br />

41<br />

40<br />

39<br />

38<br />

37<br />

36<br />

35<br />

42<br />

40<br />

38<br />

36<br />

A<br />

A<br />

D<br />

D<br />

C<br />

C<br />

A<br />

A<br />

AB<br />

D<br />

C<br />

A<br />

A<br />

D<br />

B<br />

C<br />

A<br />

A<br />

B<br />

D<br />

5 10 15<br />

<br />

A<br />

B<br />

D<br />

C C<br />

A<br />

A<br />

B<br />

D<br />

A<br />

A<br />

D<br />

B<br />

5 10 15<br />

<br />

<br />

4 × 4 8 × 8 4 × 4<br />

8 × 8<br />

C<br />

BA<br />

C<br />

D<br />

B<br />

A<br />

D<br />

C


26.5<br />

26.0<br />

25.5<br />

25.0<br />

33.0<br />

32.5<br />

32.0<br />

31.5<br />

31.0<br />

30.5<br />

30.0<br />

A<br />

A<br />

D<br />

C C<br />

D<br />

B<br />

A A<br />

D<br />

B<br />

5 10 15<br />

<br />

A<br />

C<br />

B<br />

D<br />

5 10 15<br />

<br />

B<br />

A<br />

B<br />

D<br />

C<br />

C<br />

A<br />

A<br />

B<br />

D<br />

B<br />

A<br />

D<br />

A<br />

C<br />

B<br />

C<br />

31.0<br />

30.5<br />

30.0<br />

29.5<br />

29.0<br />

28.5<br />

33.5<br />

33.0<br />

32.5<br />

32.0<br />

31.5<br />

31.0<br />

30.5<br />

A<br />

A<br />

D<br />

B<br />

C<br />

A A A<br />

D<br />

B<br />

5 10 15<br />

<br />

AC<br />

B<br />

D D<br />

5 10 15<br />

<br />

<br />

4 × 4 8 × 8 4 × 4<br />

8 × 8<br />

A<br />

B<br />

C<br />

C<br />

A<br />

D<br />

D<br />

B<br />

B<br />

A<br />

A<br />

C<br />

C<br />

D


30<br />

29<br />

28<br />

27<br />

26<br />

34<br />

33<br />

32<br />

31<br />

30<br />

29<br />

28<br />

A<br />

A<br />

D<br />

D<br />

C<br />

A<br />

C<br />

A<br />

C<br />

B<br />

D<br />

A<br />

C<br />

A<br />

C<br />

B<br />

A<br />

D D<br />

5 10 15<br />

C C<br />

D<br />

A<br />

<br />

B<br />

A<br />

D<br />

C<br />

A<br />

D<br />

A<br />

C<br />

A<br />

BD<br />

C C<br />

5 10 15<br />

<br />

B<br />

A<br />

D<br />

B<br />

A<br />

33<br />

32<br />

31<br />

30<br />

29<br />

28<br />

34<br />

32<br />

30<br />

28<br />

A<br />

A<br />

D<br />

D<br />

C<br />

C<br />

CB<br />

A<br />

A<br />

C<br />

A<br />

D<br />

D<br />

A<br />

C<br />

A<br />

D<br />

B<br />

5 10 15<br />

C C<br />

D<br />

A<br />

<br />

BA<br />

D<br />

C<br />

A<br />

D<br />

C<br />

B<br />

A<br />

D<br />

A<br />

B<br />

C C<br />

5 10 15<br />

<br />

<br />

4×4 8×8 4×4 <br />

8 × 8<br />

D<br />

A<br />

B<br />

C<br />

A


37<br />

36<br />

35<br />

34<br />

40<br />

39<br />

38<br />

37<br />

36<br />

35<br />

34<br />

A<br />

A<br />

D<br />

D<br />

C<br />

B<br />

C<br />

A<br />

B<br />

D<br />

A<br />

C<br />

5 10 15<br />

B<br />

<br />

A<br />

B<br />

A<br />

D<br />

A<br />

D<br />

B<br />

D<br />

5 10 15<br />

<br />

C<br />

A<br />

B<br />

C<br />

D<br />

A<br />

C<br />

B<br />

A<br />

D<br />

B<br />

C<br />

A<br />

D<br />

39<br />

38<br />

37<br />

36<br />

35<br />

34<br />

40<br />

39<br />

38<br />

37<br />

36<br />

35<br />

A<br />

A<br />

D<br />

D<br />

C<br />

B<br />

C<br />

B<br />

A<br />

A<br />

D<br />

B<br />

5 10 15<br />

<br />

A<br />

B<br />

A<br />

D<br />

C<br />

C<br />

A<br />

B<br />

D<br />

A<br />

B<br />

D<br />

5 10 15<br />

<br />

<br />

4 × 4 8 × 8 <br />

4 × 4 8 × 8<br />

C<br />

D<br />

A<br />

B<br />

A<br />

C<br />

D<br />

B


29<br />

28<br />

27<br />

26<br />

25<br />

24<br />

30<br />

28<br />

26<br />

24<br />

A<br />

D<br />

A<br />

C<br />

D<br />

C<br />

B<br />

C C C<br />

D<br />

A<br />

B<br />

A<br />

A<br />

D<br />

5 10 15<br />

<br />

B<br />

A<br />

B<br />

C C<br />

B<br />

D<br />

A<br />

5 10 15<br />

<br />

A<br />

B<br />

D<br />

A<br />

A<br />

D<br />

B<br />

30<br />

28<br />

26<br />

24<br />

30<br />

28<br />

26<br />

24<br />

A<br />

A<br />

D<br />

D<br />

C<br />

C<br />

B<br />

A<br />

C<br />

A<br />

D A D A<br />

B<br />

5 10 15<br />

<br />

B<br />

A<br />

B<br />

A<br />

A<br />

B<br />

C<br />

C C<br />

D<br />

5 10 15<br />

<br />

<br />

4 × 4 8 × 8 4 × 4<br />

8 × 8<br />

B<br />

D<br />

A<br />

B


40<br />

39<br />

38<br />

37<br />

36<br />

35<br />

34<br />

33<br />

42<br />

41<br />

40<br />

39<br />

38<br />

37<br />

36<br />

C<br />

A<br />

D<br />

A<br />

D<br />

AD A<br />

B<br />

C<br />

C<br />

B<br />

A<br />

D<br />

B<br />

C<br />

A<br />

D<br />

B<br />

5 10 15<br />

B<br />

A<br />

<br />

D<br />

A<br />

C<br />

B<br />

A<br />

D<br />

5 10 15<br />

<br />

A<br />

C<br />

A<br />

B<br />

C<br />

A<br />

D<br />

A<br />

DB<br />

C<br />

B<br />

C<br />

A C<br />

A<br />

B<br />

D<br />

D<br />

B<br />

41<br />

40<br />

39<br />

38<br />

37<br />

36<br />

35<br />

42<br />

41<br />

40<br />

39<br />

38<br />

37<br />

36<br />

A<br />

D<br />

A<br />

D<br />

C<br />

C<br />

B A<br />

D<br />

C<br />

A<br />

B<br />

D<br />

A<br />

A<br />

C<br />

B<br />

D<br />

5 10 15<br />

B<br />

A<br />

<br />

D<br />

A<br />

C<br />

B<br />

A<br />

D<br />

5 10 15<br />

<br />

C<br />

A<br />

B<br />

A<br />

A<br />

D<br />

C<br />

A C<br />

B<br />

<br />

4 × 4 8 × 8 4 × 4<br />

8 × 8<br />

C<br />

B<br />

D<br />

A<br />

D<br />

B<br />

C


I2t<br />

I2t+1<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

Dnr I2t<br />

arrondi(v) <br />

x ∈ Dnr ⇐⇒ ∀y ∈ I2t+1, y − arrondi(v(y)) = x


I2t <br />

arrondi(v) <br />

Dnc I2t+1<br />

arrondi(v) <br />

<br />

<br />

x ∈ Dnc ⇐⇒ ∃y ∈ I2t, ∃x ′ ∈ I2t+1, y = x−arrondi(v(x)) = x ′ −arrondi(v(x ′ )),<br />

|I2t(y) − I2t+1(x ′ )| < |I2t(y) − I2t+1(x ′ )|. <br />

Dnr Dnc <br />

<br />

<br />

˜ J J <br />

<br />

H(x) = 1<br />

<br />

√ I2t+1(x) −<br />

2<br />

Ĩ2t(x<br />

<br />

− v(x))<br />

<br />

<br />

L(y) = √ 2I2t(y).<br />

y I2t x I2t+1 y = x−v(x)<br />

y <br />

L(y) = √ 2I2t(y) + ˜ H(y + v(x))<br />

<br />

I2t I2t+1 <br />

L H


I2t <br />

arrondi(v) <br />

Dnc I2t+1<br />

arrondi(v) <br />

<br />

<br />

x ∈ Dnc ⇐⇒ ∃y ∈ I2t, ∃x ′ ∈ I2t+1, y = x−arrondi(v(x)) = x ′ −arrondi(v(x ′ )),<br />

|I2t(y) − I2t+1(x ′ )| < |I2t(y) − I2t+1(x ′ )|. <br />

Dnr Dnc <br />

<br />

<br />

˜ J J <br />

<br />

H(x) = 1<br />

<br />

√ I2t+1(x) −<br />

2<br />

Ĩ2t(x<br />

<br />

− v(x))<br />

<br />

<br />

L(y) = √ 2I2t(y).<br />

y I2t x I2t+1 y = x−v(x)<br />

y <br />

L(y) = √ 2I2t(y) + ˜ H(y + v(x))<br />

<br />

I2t I2t+1 <br />

L H


' () * $ &<br />

!"# $%" &<br />

<br />

<br />

&' ( # %<br />

!" #$" %


) *+ , &" (<br />

"#$% &'$ (<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

!


) *+ , &" (<br />

"#$% &'$ (<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

!


) *+ , &" (<br />

"#$% &'$ (<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

!


F [k] N <br />

B = {gk}0≤k


B ∗ = {g ∗ k }0≤k


]yk, yk+1]<br />

p(x) pk ∆k <br />

<br />

∆k = yk+1 − yk,<br />

<br />

∀x ∈]yk, yk+1] p(x) = pk<br />

∆k<br />

<br />

p(x) <br />

xk<br />

[yk, yk+1] <br />

d = 1 <br />

pk∆k<br />

12<br />

<br />

<br />

<br />

∀k ∈ Z yk+1 − yk = ∆.<br />

d = ∆<br />

12 .<br />

X <br />

H(X) = −<br />

+∞<br />

−∞<br />

p(x) log 2 p(x)dx.<br />

<br />

X ¯ X <br />

H( ¯ X) ≥ Hd(X) − 1<br />

2 log 2(12d)<br />

Q <br />

¯ X RX <br />

¯ X <br />

RX = H( ¯ H) = Hd(X) − log 2(12d).<br />

<br />

d(RX) = 1<br />

12 22Hd(X) 2 −2RX .


FB[k] <br />

<br />

∆k = 12d<br />

N<br />

0 ≤ k < N<br />

d( ¯ R) = N Hd<br />

¯ −2<br />

22 2<br />

12 ¯ R<br />

,<br />

¯ Hd


D <br />

H <br />

R <br />

G <br />

F <br />

<br />

M <br />

<br />

A <br />

AM AM <br />

<br />

W <br />

∆d R/F <br />

B<br />

<br />

R = D/F<br />

∆d = 0<br />

H = Rprec<br />

W = AM<br />

<br />

<br />

W > M <br />

D/F


M<br />

AM<br />

W<br />

0<br />

∆b<br />

<br />

<br />

AM<br />

M<br />

W<br />

0<br />

∆b<br />

<br />

<br />

AM<br />

W<br />

M<br />

0<br />

Saut <strong>de</strong> trame<br />

∆b<br />

<br />

<br />

AM W < AM <br />

W > AM


W < AM <br />

AM − W <br />

W > AM <br />

W/F <br />

<br />

<br />

<br />

<br />

<br />

<br />

W ← max(W + H − D/F, 0)<br />

skip ← 1<br />

W > M <br />

W ← max(W − D/F, 0)<br />

skip ← skip + 1<br />

<br />

skip.G/F − 1 <br />

W > AM <br />

∆d ← W/F<br />

<br />

∆d ← W − AM<br />

<br />

B ← D/F − ∆d<br />

R<br />

<br />

<br />

<br />

<br />

QP <br />

<br />

QPprec <br />

<br />

<br />

<br />

<br />

W/F


Npix <br />

<br />

Ri i <br />

σi i <br />

∆i i <br />

Nmb <br />

N = NpixNmb<br />

C <br />

<br />

K <br />

<br />

<br />

Ri = Npix(K σ2 i<br />

∆ 2 i<br />

+ C).<br />

¯ Xi <br />

∆i X <br />

σ2<br />

∆ 2 i<br />

< 1<br />

2e <br />

H( ¯ Xi) = e<br />

ln2<br />

σ2 i<br />

∆2 i<br />

. <br />

K e<br />

ln2 <br />

<br />

<br />

K C <br />

<br />

i i − 1<br />

<br />

<br />

<br />

αi <br />

<br />

αi =<br />

<br />

2 R<br />

N (1 − σi) + σi<br />

1 <br />

R < 0.5N


i i Nmb<br />

ni i <br />

ni = Nmb + 1 − i.<br />

Li <br />

Li = ri − NpixniC.<br />

Si <br />

Nmb<br />

Si =<br />

<br />

αiσi.<br />

<br />

<br />

σi <br />

Li > 0 <br />

<br />

∆ ∗ i =<br />

j=i<br />

NpixKiσi<br />

Liαi<br />

<br />

∆ ∗ i <br />

<br />

C K <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

ρ R Rz<br />

Rnz <br />

<br />

<br />

<br />

Si,


(a[k], r[k]) <br />

( ¯ FB[k])1≤k≤Npix ↦→ ((a[k], r[k]))1≤k≤m<br />

a[k] k <br />

r[k] <br />

m Npix <br />

<br />

<br />

ρ = 1<br />

N<br />

Nmb mi <br />

r[k],<br />

i=0 k=0<br />

<br />

ν = 1<br />

Nmb <br />

N<br />

Qnz <br />

<br />

Qnz = 1<br />

N<br />

i=0<br />

mi<br />

Nmb mi <br />

⌈log2 |a[k]|⌉ + 1.<br />

i=0 k=0<br />

Qz <br />

<br />

Qz = 1<br />

nN<br />

Nmb mi <br />

⌈log2 |r[k]|⌉ + 1.<br />

i=0 k=0<br />

Qnz Qz ∆ <br />

<br />

<br />

ρ <br />

Qnz ρ <br />

Qnz = ϑ(1 − ρ). <br />

ϑ <br />

ϑ Qz ϑ ρ <br />

<br />

Qz = Pρ(ϑ) <br />

Pρ <br />

ρ


R(ρ) = ξ1(ρ)Qnz(ρ) + ξ2(ρ)Qz(ρ) + ξ3(ρ). <br />

ξ1 ξ2 ξ3 <br />

<br />

Pρ ξ1 ξ2 ξ3 <br />

ρ <br />

R ρ R <br />

<br />

R ρ <br />

ρ <br />

T <br />

ρ ∆ ∗ θ = T/∆ <br />

∆ ∗ ↦→ ρ <br />

<br />

<br />

<br />

<br />

ρ <br />

0.9 < ρ < 0.95.<br />

Qnz <br />

ϑ<br />

ϑ Qz<br />

ρ R<br />

ρ ↦→ R<br />

<br />

R ↦→ ρ<br />

ρ QP0 <br />

<br />

QP <br />

ρ <br />

<br />

ρ + <br />

ρ −


κ = ρ−<br />

,<br />

ρ +<br />

<br />

QP0 − 3 QP0 + 3 κ <br />

<br />

κ <br />

<br />

<br />

<br />

<br />

<br />

Rz<br />

Rnz <br />

<br />

ρ Rz <br />

<br />

FB[k] s[k] <br />

s[k] =<br />

<br />

0 ¯ FB[k] = 0<br />

1 ¯ FB[k] = 0<br />

<br />

S ρ ν = 1 − ρ<br />

<br />

H(S) = −ρ log 2 ρ + ν log 2 ν.<br />

rz <br />

rz = Rz<br />

νN<br />

<br />

S rz <br />

Rz<br />

N ≤ −ρ log 2 ρ + ν log 2 ν.<br />

x ∈ [0, 1[ −x log 2 x ≤ (1 − x) log 2 e <br />

rz = Rz<br />

νN ≤ log 2 e − log 2 ν


ρ <br />

<br />

rz <br />

rz <br />

rz <br />

<br />

<br />

<br />

<br />

<br />

<br />

T <br />

∆ θ = T/∆ <br />

<br />

<br />

<br />

νN <br />

<br />

<br />

rnz <br />

rnz = Rnz<br />

νN<br />

νN ≫ 1 νN <br />

<br />

pT (x) = 1<br />

ν p(x)1 {|x|>T }(x).<br />

XT pT <br />

<br />

<br />

rnz = HdXT − log 2 ∆<br />

∼ <br />

<br />

f ∼ g ⇐⇒ g = O(f) f = O(g)<br />

νN > 1<br />

ɛ ν < ɛ ɛ<br />

F r B [k] FB[k] <br />

<br />

k < k ′ ⇒ F r B[k] ≥ F r B[k ′ ]


|F r B[k]| ∼ k −γ(k) , <br />

γ k<br />

N <br />

rnz γ(ν) θ γ(ν) <br />

rnz <br />

γ γ > 0.5 <br />

<br />

γ = 1 <br />

1/2 < γ ≤ 1 γ <br />

<br />

T <br />

<br />

νN ∼ T −γ .<br />

T ∆ <br />

νN ∼ ∆ −γ . <br />

Rnz = Ra + Rs Ra <br />

Rs <br />

<br />

<br />

0 ≤ Rs ≤ νN.<br />

pj <br />

(j + 1/2)∆ + T <br />

lj = − log 2 pj <br />

<br />

Ha = −νN<br />

+∞<br />

j=0<br />

pj log 2 pj. <br />

nj = νNpj j |F r B [k]| <br />

|F r B<br />

<br />

[k]| ∈ [(j + θ)∆, (j + 1 + θ)∆[ <br />

1<br />

1<br />

− −<br />

nj ∼ ((j + θ)∆) γ − ((j + θ + 1)∆) γ .<br />

<br />

pj = nj<br />

νN<br />

1<br />

1<br />

−<br />

∼ (j + θ)− γ − (j + θ + 1) γ .<br />

+∞<br />

j=0 pj log 2 pj νN <br />

<br />

Ha ∼ νN


s <br />

s <br />

s = 1/2 Ra νN<br />

<br />

Rnz = Ra + Rs ∼ νN,<br />

rnz <br />

<br />

Rz <br />

Rnz <br />

rz rnz <br />

<br />

<br />

r ∗ = rz + rnz <br />

R<br />

N = Rz + Rnz<br />

N<br />

<br />

= rzν + rnzν <br />

= r ∗ (1 − ρ) <br />

ρ <br />

γ = 1/2 <br />

R<br />

N<br />

1<br />

∼ K ,<br />

∆2 <br />

<br />

<br />

r ∗ <br />

<br />

r ∗ <br />

<br />

r ∗


Qnz Qz


∆ = W − AD/F W<br />

AM <br />

<br />

<br />

<br />

<br />

<br />

ρ ↦→ R


∗<br />

r ∗ <br />

<br />

r ∗ r ∗ <br />

r ∗ t − 1 t − 2<br />

t − 3 r ∗ 1 r ∗ 2 r ∗ 3 r ∗ t <br />

˜r ∗ = 0.6r ∗ 1 + 0.3r ∗ 2 + 0.1r ∗ 3.<br />

<br />

<br />

r ∗ = r ∗ 0<br />

QP <br />

r ∗ QP r ∗ QP r ∗ QP <br />

r ∗ <br />

<br />

L <br />

<br />

N = 25344 Nmb = 99<br />

<br />

<br />

<br />

<br />

3<br />

2 Nc1 + Nmbc2<br />

Nmbc3<br />

Nmbc4<br />

NCtmn<br />

<br />

<br />

<br />

QP ↦→ ρ ϑ <br />

<br />

<br />

<br />

ϑ <br />

Qnz


Nmbh1<br />

3<br />

2 Nh2<br />

Lh3<br />

ρ − R 32h4<br />

ϑ Lh5<br />

R − QP QP0 h6<br />

Nmbh7<br />

Nmbh8<br />

NCh1 + LCh2<br />

<br />

<br />

θ <br />

ρ ↦→ R <br />

<br />

<br />

<br />

<br />

<br />

θ <br />

<br />

Nmbh1<br />

3<br />

2 Nh2<br />

Lh3<br />

ρ − QP 32h4<br />

ρ − R QP0 a6<br />

Nmba7<br />

Nmba8<br />

a6 < h6 a7 < h7 a8 < h8 NCa1 +LCa2 Ca1 < Ch1<br />

Ca2 < Ch2


ρ


QP <br />

<br />

r ∗ <br />

<br />

<br />

r ∗ 0


∗ <br />

r ∗


∗ <br />

r ∗


0.7<br />

0.6<br />

0.5<br />

0.4<br />

0.3<br />

0.2<br />

A<br />

C<br />

B<br />

C<br />

B<br />

A<br />

CB<br />

3.0 3.5 4.0 4.5 5.0 5.5 6.0<br />

R(ϑ, ρ) ρ = 0.96 ρ = 0.97 ρ = 0.98<br />

R <br />

ϑ <br />

A


∗ 0 = 6.7 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

r∗ 0 = 6.7 <br />

3.5 ϑ 5.4 ρ ↦→ R <br />

<br />

ρ <br />

R ρ <br />

<br />

<br />

<br />

ϑ


∗ 0 = 6.7 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

r∗ 0 = 8.0 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

r∗ 0 = 8.0


∗ 0 = 8.0 <br />

<br />

r ∗<br />

r ∗ 0 <br />

r ∗ <br />

<br />

<br />

r ∗ <br />

<br />

<br />

r ∗ <br />

<br />

<br />

<br />

<br />

r ∗ <br />

r ∗ <br />

r ∗ <br />

<br />

r ∗ 0<br />

<br />

r ∗ <br />

<br />

<br />

<br />

<br />

<br />

r ∗


∗ <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

r∗ <br />

<br />

r ∗ 0<br />

<br />

= 8.0


4000<br />

3000<br />

2000<br />

1000<br />

0<br />

0 50 100 150 200 250<br />

<br />

4000<br />

3000<br />

2000<br />

1000<br />

0<br />

4000<br />

3000<br />

2000<br />

1000<br />

0 50 100 150 200 250<br />

r ∗ = 6.7<br />

0<br />

0 50 100 150 200 250<br />

r ∗ = 8<br />

<br />

<br />

<br />

<br />

<br />

r ∗ QP r ∗ QP r ∗ <br />

QP r ∗ QP


3000<br />

2500<br />

2000<br />

1500<br />

1000<br />

500<br />

0<br />

0 50 100 150 200 250<br />

<br />

2500<br />

2000<br />

1500<br />

1000<br />

500<br />

0<br />

2000<br />

1500<br />

1000<br />

500<br />

0 50 100 150 200 250<br />

r ∗ = 6.7<br />

0<br />

0 50 100 150 200 250<br />

r ∗ = 8<br />

<br />

<br />

<br />

<br />

<br />

r ∗ QP r ∗ QP r ∗ <br />

QP r ∗ QP


S(v) = |err| + λ(QP )C(v),<br />

λ C(v)


C(v) <br />

v <br />

¯v


C(v) <br />

v <br />

¯v


H(X)<br />

<br />

H(X) = − <br />

p(x) log p(x) <br />

x∈X<br />

X X p <br />

<br />

<br />

p


H(X, Y ) (X, Y ) <br />

p(x, y) <br />

H(X, Y ) = − <br />

p(x, y) log p(x, y) <br />

x∈X y∈Y<br />

<br />

<br />

<br />

H(Y |X) <br />

H(X, Y ) = − <br />

p(x, y) log p(y|x) <br />

x∈X y∈Y<br />

H(X|Y ) = H(Y |X) X = ⌊Y/5⌋ <br />

Y X <br />

<br />

<br />

<br />

p q <br />

D(p q) = <br />

x∈X<br />

p(x) log p(x)<br />

q(x)<br />

<br />

<br />

q <br />

X p <br />

<br />

<br />

<br />

<br />

I(X; Y ) <br />

<br />

I(X; Y ) = I(Y ; X) = <br />

p(x, y)<br />

p(x, y) log<br />

p(x)p(y)<br />

x∈X y∈Y<br />

<br />

I(X; Y ) = H(X) − H(X|Y ) = H(Y ) − H(Y |X) = H(X) + H(Y ) −<br />

H(X, Y )


n (V i<br />

1<br />

, V i<br />

2 )0


20<br />

15<br />

10<br />

5<br />

0<br />

−5<br />

−10<br />

−10 −5 0 5 10 15 20<br />

<br />

<br />

<br />

10<br />

0<br />

−10<br />

−20<br />

−15 −10 −5 0 5 10


vx vy


140<br />

120<br />

100<br />

80<br />

60<br />

40<br />

20<br />

0<br />

0 50 100<br />

<br />

150<br />

∂vx<br />

∂x<br />

140<br />

120<br />

100<br />

80<br />

60<br />

40<br />

20<br />

0<br />

0 50 100<br />

<br />

150<br />

∂vx<br />

∂y<br />

140<br />

120<br />

100<br />

80<br />

60<br />

40<br />

20<br />

0<br />

0 50 100<br />

<br />

150<br />

∂vy<br />

∂x<br />

140<br />

120<br />

100<br />

80<br />

60<br />

40<br />

20<br />

0<br />

0 50 100<br />

<br />

150<br />

∂vy<br />

∂y<br />

∇v


(x k 1,n )0≤n≤N (x k 2,n )0≤n≤N <br />

(i, j) ∈ {(1, 2), (2, 1)} <br />

˜x k i,n,α = |x k i,2n+1 − Pα((x k i,2n))|,<br />

Pα <br />

αl = argmin ˜xi,l,α.<br />

α<br />

2l+1 (x k j,n )0≤n≤N Pαl <br />

Uαl <br />

Pαl<br />

<br />

<br />

(x k i,n )0≤n≤N (x k j,n )0≤n≤N <br />

<br />

<br />

<br />

<br />

<br />

<br />

˜xj,l,β = |x k+1<br />

j,l<br />

− Pβ((x k+1<br />

j,l ))|,<br />

Pβ β <br />

<br />

γl = 1 + 2 argmin ˜xj,l,β.<br />

β<br />

2l + 1 Pγl <br />

Uγl


P lin<br />

−1 P lin<br />

0 P lin<br />

1 <br />

P lin<br />

−1 ((x2k)) = ⌊− 1<br />

2 (x 2(k−1)) + 3<br />

2 (x2k)⌋<br />

P lin<br />

0 ((x2k)) = ⌊ 1<br />

2 (x2k) + 1<br />

2 (x 2(k+1))⌋<br />

P lin<br />

1 ((x2k)) = ⌊ 3<br />

2 (x 2(k+1) − 1<br />

2 (x 2(k+1))⌋<br />

P lin<br />

0 <br />

P lin<br />

−1 P lin<br />

1 <br />

<br />

(w t s) s∈{−1,0,1},t∈{d,b}<br />

<br />

<br />

w d −1 = w d 1 = 2 w d 0 = 1<br />

w b −1 = w b 1 = 1 w b 0 = 1<br />

2<br />

αl = argmin w<br />

α∈{−1,0,1}<br />

d α|x k i,2l+1 − Pα((x k i,2n))| + w b α<br />

αl l <br />

(xj,n) w0 <br />

αl <br />

<br />

αl <br />

U−1 U0 <br />

U1 <br />

U lin<br />

−1 ((x2k+1)) = ⌊− 1<br />

4 (x2k−3) + 3<br />

4 (x2k−1)⌋<br />

U lin<br />

0 ((x2k+1)) = ⌊ 1<br />

4 (x2k−1) + 1<br />

4 (x2k+1)⌋<br />

U lin<br />

1 ((x2k+1)) = ⌊ 3<br />

4 (x 2k+1) − 1<br />

4 (x2k+3)⌋<br />

αl ≤ 0 x k i,2l<br />

x k i,2l+1 2l + 1 <br />

<br />

2l − 3 αl−1 ≥ 0 <br />

x k i,2l x k i,2l−1 2l − 1 <br />

2l + 3


αl ≤ 0 αl > 0<br />

αl−1 ≥ 0 U lin<br />

0<br />

αl−1 < 0 U lin<br />

1<br />

U lin<br />

−1<br />

U lin<br />

0<br />

(αl)<br />

αl > 0 αl−1 < 0 <br />

<br />

2l + 1 2l − 1 <br />

<br />

<br />

(x l i,n )<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

(x1) <br />

<br />

<br />

<br />

<br />

Pα <br />

n ˜xi,α,n <br />

<br />

ρi <br />

ρi =<br />

<br />

α ′ =α (˜xi,α ′ ,n + εα ′)<br />

(A − 1) <br />

α ′(˜xi,α ′ ,n + εα ′),<br />

A <br />

εα


εα <br />

<br />

εα


s <br />

Hs ns <br />

R <br />

R = <br />

nsHs.<br />

s<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

vx <br />

120<br />

100<br />

80<br />

60<br />

40<br />

20<br />

0<br />

0 20 40 60 80 100 120<br />

vx<br />

120<br />

100<br />

80<br />

60<br />

40<br />

20<br />

0<br />

0 20 40 60 80 100 120<br />

vy


vy <br />

vy <br />

vx <br />

vx <br />

vy <br />

vy <br />

vx <br />

vx


vy <br />

vy <br />

vx <br />

vx <br />

vy <br />

vy <br />

vx <br />

vx


vy <br />

vy <br />

vx <br />

vx <br />

vy <br />

vy <br />

vx <br />

vx


1.0<br />

0.5<br />

0.0<br />

-0.5<br />

-1.0<br />

0 2 4 6 8 10<br />

<br />

1.0<br />

0.5<br />

0.0<br />

-0.5<br />

-1.0<br />

0 2 4 6 8 10


x2k+1 x2k x2k+1 <br />

P 2p − 1 2p <br />

<br />

2k + 1 <br />

x2k−2p+2, . . . , x2k, x2k+2, . . . , x2k+2p<br />

x2k+1 = P (2k + 1).<br />

<br />

Mn (n, xn) <br />

Mn (n, xn)<br />

<br />

hk = x2k+1 − x2k+1<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

2k 2(k + 1) [xk ∧ x 2(k+1), x2k ∨ x 2(k+1)]


x−4 x−2 x0 x2 x4 x6<br />

ˆx3<br />

<br />

n<br />

Pn <br />

Pn,0 <br />

⎧<br />

⎪⎨<br />

Pn,0(x) =<br />

⎪⎩<br />

x2k ∧ x 2(k+1) Pn(x) ≤ x2k ∧ x 2(k+1)<br />

x2k ∨ x 2(k+1) Pn(x) ≥ x2k ∨ x 2(k+1)<br />

Pn(x) <br />

[2k, 2(k + 1)]<br />

Pn,0 <br />

Pn,0(x) = ((x2k ∧ x 2(k+1)) ∨ Pn(x)) ∧ (x2k ∨ x 2(k+1))<br />

= med(x2k, Pn(x), x 2(k+1))<br />

Pn,0 <br />

x−4 x−2 x0 x2 x4 x6<br />

ˆx1<br />

P3,0<br />

∧ ∨ Pn,0 <br />

<br />

Pn,0


x2k, x2k+2, . . . , x2l<br />

<br />

<br />

x2k, x2k+1, x2k+2, . . . , x2l<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

x1 (1, x1) ∆1 ∆2 ∆1<br />

(M−2M2) M0 ∆2 (M0, M2) <br />

M4 M2 <br />

M0 ∆ ′ 1 ∆ ′ 2 <br />

Pk,lm <br />

Pk;l,m(2k) = x2k<br />

Pk;l,m(2m) − Pk;l,m(2l) = x2m − x2l<br />

2k + 1 <br />

Pk;k−1,k+1(x) Pn,1(x) =<br />

∨ Pk;k,k+1(x) ∧ Pk+1;k,k+2(x) ∨ Pk+1;k,k+1(x) <br />

∧ Pn(x)∨<br />

Pk;k−1,k+1(x) ∧ Pk;k,k+1(x) ∨ Pk+1;k,k+2(x) ∧ Pk+1;k,k+1(x) <br />

<br />

<br />

= med(Pk;k−1,k+1(x), med(Pk+1;k,k+2(x), Pn(x), Pk+1;k,k+1(x)), Pk+1;k,k+1(x))


x−2<br />

x−2<br />

M−2<br />

(∆1)<br />

M0<br />

x0<br />

(∆2)<br />

ˆx1<br />

M2<br />

x2<br />

M4<br />

<br />

M−2<br />

M0<br />

x0<br />

ˆx1<br />

(∆ ′ 2)<br />

M2<br />

x2<br />

(∆ ′ 1)<br />

M4<br />

<br />

P3,1<br />

x4<br />

x4


k<br />

f<br />

I ⊂ R R x ∈ I <br />

f <br />

∀ε > 0, ∃x ′ , x ′′ ∈ I, |x − x ′ | < ε, |x − x ′′ | < ε f(x ′ ) < 0, f(x ′′ ) > 0 <br />

f I ⊂ R R <br />

<br />

f I ⊂ R R k <br />

k − 1 i ∈ {0, . . . , k − 1} f (i) <br />

<br />

<br />

f <br />

f(x) =<br />

<br />

0 x = 0<br />

x 9 sin 1<br />

x<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

k k <br />

<br />

<br />

<br />

<br />

P3,1 <br />

C 4 <br />

O(2 −4j ) j <br />

<br />

k k <br />

k O(2 −jk )


k <br />

<br />

<br />

pI I = [a, b] <br />

⎧<br />

⎪⎨ a x < a<br />

pI(x) = x a ≤ x ≤ b<br />

⎪⎩<br />

b <br />

<br />

<br />

I1 = [a1, b1] <br />

I2 = [a2, b2] I1 ∩ I2 = ∅ I1 ∩ I2 <br />

<br />

I1 ∩ I2 = ∅ =⇒ pI1 ◦ pI2 = pI2 ◦ pI1 = pI1∩I2 .<br />

I1 ∩ I2 = ∅ I1<br />

I2 <br />

<br />

I1 ∩ I2 = ∅ =⇒ pI1<br />

◦ pI2 = pI2 ◦ pI1<br />

<br />

<br />

<br />

<br />

I1 I2 <br />

Rk k <br />

k <br />

intmed(I1, I2) <br />

intmed(I1, I2) = [R1(a1, b1, a2, b2), R2(a1, b1, a2, b2)].<br />

I1 I2 <br />

intmed(I1, I2) <br />

<br />

<br />

intmed(I1, I2) ∩ I1 = ∅<br />

intmed(I1, I2) ∩ I2 = ∅


k <br />

<br />

intmed(([aj, bj])1≤j≤k) = [Rk(a1, . . . , ak, b1, . . . , bk), Rk+1(a1, . . . , ak, b1, . . . , bk)].<br />

<br />

]−∞, b] [a, +∞[ <br />

<br />

<br />

<br />

Ik;d(x) Pk;k−d,k+d(x) Pk;k,k+d(x) <br />

Pn,1(x) 2k + 1 <br />

Pn,1(x) = p intmed(Ik,1(x),Ik+1,−1(x))(Pn(x))<br />

<br />

<br />

<br />

N <br />

<br />

<br />

N <br />

<br />

<br />

<br />

<br />

<br />

N <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

N


4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N <br />

<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N


6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

100<br />

50<br />

0<br />

-50<br />

-100<br />

-150<br />

-200<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N <br />

<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

100<br />

50<br />

0<br />

-50<br />

-100<br />

-150<br />

-200<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N


4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N <br />

<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

250<br />

200<br />

150<br />

100<br />

50<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N


6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

100<br />

50<br />

0<br />

-50<br />

-100<br />

-150<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N <br />

<br />

6<br />

5<br />

4<br />

3<br />

2<br />

1<br />

0<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

<br />

<br />

100<br />

50<br />

0<br />

-50<br />

-100<br />

-150<br />

0.0 0.2 0.4 0.6 0.8 1.0<br />

N


PSNR<br />

50<br />

45<br />

40<br />

35<br />

30<br />

25<br />

20<br />

15<br />

1.0 2. 5. 10.0 20. 50. 100.0 200.<br />

nombre <strong>de</strong> coefficients<br />

<br />

N


EE EO<br />

OE OO <br />

EE<br />

OE<br />

EO<br />

OO<br />

P<br />

P<br />

U<br />

U<br />

P<br />

<br />

x <br />

y <br />

<br />

⊗ I2 I2 ⊗<br />

<br />

I U<br />

0 I<br />

P<br />

U<br />

U<br />

I 0<br />

P I<br />

<br />

EE<br />

OE<br />

EO<br />

OO<br />

P<br />

P<br />

P<br />

P<br />

U<br />

<br />

<br />

<br />

<br />

<br />

<br />

P/2 <br />

<br />

U<br />

U<br />

h2k+1 = x2k+1 − ˆx2k+1/2<br />

U<br />

LL<br />

HL<br />

LH<br />

HH<br />

LL<br />

HL<br />

LH<br />

HH


EE<br />

OE<br />

EO<br />

OO<br />

P/2<br />

P/2<br />

P<br />

<br />

h2k+1 = x2k+1 − ˆx2k+1<br />

<br />

<br />

<br />

<br />

EE<br />

OE<br />

EO<br />

OO<br />

PHL<br />

PLH<br />

P<br />

PHH UHH<br />

<br />

<br />

<br />

<br />

<br />

P HL<br />

n,0 (x) = med(x2k,2l, x, x2k+2,2l)<br />

P LH<br />

n,0 (x) = med(x2k,2l, x, x2k,2l+2)<br />

P HH<br />

n,0 (x) = med(x2k,2l, x2k,2l+2, x, x2k+2,2l, x2k+2,2l+2)<br />

<br />

<br />

<br />

P/2<br />

P/2<br />

ULH<br />

UHL<br />

LL<br />

HL<br />

LH<br />

HH<br />

LL<br />

HL<br />

LH<br />

HH


Pk;l,m(2k) = x2k<br />

Pk;l,m(2m) − Pk;l,m(2l) = x2m − x2l<br />

<br />

<br />

Ik,∆(x) Pk;k−∆,k+∆(x) Pk;k,k+∆(x)<br />

<br />

∆ HH<br />

1 = (1, 1) ∆ HH<br />

2 = (−1, 1)<br />

∆ LH = (0, 1) ∆ HL = (1, 0)<br />

2k + ∆ HH<br />

1<br />

P HL<br />

n,1 (x) = p intmed(Ik,∆ HL(x),I k+∆ HL ,−∆ HL(x))(x)<br />

P LH<br />

n,1 (x) = p intmed(Ik,∆ LH (x),I k+∆ LH ,−∆ LH (x))(x)<br />

<br />

P HH<br />

n,1 (x) = pintmed(Ik,∆HH (x),I<br />

1<br />

k+∆HH 1<br />

,−∆HH (x),I<br />

1<br />

k,∆HH (x),I<br />

2<br />

k+∆HH 2<br />

,−∆HH (x))(x)<br />

2


Pk;l,m(2k) = x2k<br />

Pk;l,m(2m) − Pk;l,m(2l) = x2m − x2l<br />

<br />

<br />

Ik,∆(x) Pk;k−∆,k+∆(x) Pk;k,k+∆(x)<br />

<br />

∆ HH<br />

1 = (1, 1) ∆ HH<br />

2 = (−1, 1)<br />

∆ LH = (0, 1) ∆ HL = (1, 0)<br />

2k + ∆ HH<br />

1<br />

P HL<br />

n,1 (x) = p intmed(Ik,∆ HL(x),I k+∆ HL ,−∆ HL(x))(x)<br />

P LH<br />

n,1 (x) = p intmed(Ik,∆ LH (x),I k+∆ LH ,−∆ LH (x))(x)<br />

<br />

P HH<br />

n,1 (x) = pintmed(Ik,∆HH (x),I<br />

1<br />

k+∆HH 1<br />

,−∆HH (x),I<br />

1<br />

k,∆HH (x),I<br />

2<br />

k+∆HH 2<br />

,−∆HH (x))(x)<br />

2


Pk;l,m(2k) = x2k<br />

Pk;l,m(2m) − Pk;l,m(2l) = x2m − x2l<br />

<br />

<br />

Ik,∆(x) Pk;k−∆,k+∆(x) Pk;k,k+∆(x)<br />

<br />

∆ HH<br />

1 = (1, 1) ∆ HH<br />

2 = (−1, 1)<br />

∆ LH = (0, 1) ∆ HL = (1, 0)<br />

2k + ∆ HH<br />

1<br />

P HL<br />

n,1 (x) = p intmed(Ik,∆ HL(x),I k+∆ HL ,−∆ HL(x))(x)<br />

P LH<br />

n,1 (x) = p intmed(Ik,∆ LH (x),I k+∆ LH ,−∆ LH (x))(x)<br />

<br />

P HH<br />

n,1 (x) = pintmed(Ik,∆HH (x),I<br />

1<br />

k+∆HH 1<br />

,−∆HH (x),I<br />

1<br />

k,∆HH (x),I<br />

2<br />

k+∆HH 2<br />

,−∆HH (x))(x)<br />

2


Pk;l,m(2k) = x2k<br />

Pk;l,m(2m) − Pk;l,m(2l) = x2m − x2l<br />

<br />

<br />

Ik,∆(x) Pk;k−∆,k+∆(x) Pk;k,k+∆(x)<br />

<br />

∆ HH<br />

1 = (1, 1) ∆ HH<br />

2 = (−1, 1)<br />

∆ LH = (0, 1) ∆ HL = (1, 0)<br />

2k + ∆ HH<br />

1<br />

P HL<br />

n,1 (x) = p intmed(Ik,∆ HL(x),I k+∆ HL ,−∆ HL(x))(x)<br />

P LH<br />

n,1 (x) = p intmed(Ik,∆ LH (x),I k+∆ LH ,−∆ LH (x))(x)<br />

<br />

P HH<br />

n,1 (x) = pintmed(Ik,∆HH (x),I<br />

1<br />

k+∆HH 1<br />

,−∆HH (x),I<br />

1<br />

k,∆HH (x),I<br />

2<br />

k+∆HH 2<br />

,−∆HH (x))(x)<br />

2


Pk;l,m(2k) = x2k<br />

Pk;l,m(2m) − Pk;l,m(2l) = x2m − x2l<br />

<br />

<br />

Ik,∆(x) Pk;k−∆,k+∆(x) Pk;k,k+∆(x)<br />

<br />

∆ HH<br />

1 = (1, 1) ∆ HH<br />

2 = (−1, 1)<br />

∆ LH = (0, 1) ∆ HL = (1, 0)<br />

2k + ∆ HH<br />

1<br />

P HL<br />

n,1 (x) = p intmed(Ik,∆ HL(x),I k+∆ HL ,−∆ HL(x))(x)<br />

P LH<br />

n,1 (x) = p intmed(Ik,∆ LH (x),I k+∆ LH ,−∆ LH (x))(x)<br />

<br />

P HH<br />

n,1 (x) = pintmed(Ik,∆HH (x),I<br />

1<br />

k+∆HH 1<br />

,−∆HH (x),I<br />

1<br />

k,∆HH (x),I<br />

2<br />

k+∆HH 2<br />

,−∆HH (x))(x)<br />

2


S<br />

N : x ↦→ x + S<br />

<br />

<br />

N : x ↦→ x + SD(x)<br />

D


v + <br />

v − <br />

<br />

<br />

M + M − <br />

<br />

M + (t, x, y) = (t + 1, x + v + x (t, x, y), y + v + y (t, x, y)) <br />

M − (t, x, y) = (t − 1, x + v − x (t, x, y), y + v − y (t, x, y)) <br />

<br />

M p ⎧<br />

⎪⎨ Id p = 0<br />

= (M<br />

⎪⎩<br />

+ ) p p > 0<br />

(M − ) −p p < 0<br />

T <br />

T (a,b)(t, x, y) = (t, x + a, y + b).<br />

<br />

M0 <br />

<br />

<br />

<br />

S x <br />

<br />

V (x) = {x} + S(x) <br />

V (x) = {(t + δt, x + δx, y + δy)|(δt, δx, δy) ∈ S(x)} <br />

<br />

VM(x) = {M δt (T (δx,δy)(t, x, y))|(δt, δx, δy) ∈ S(x)}


WM(x) = {T (δx,δy)(M δt (t, x, y))|(δt, δx, δy) ∈ S(x)}, <br />

<br />

<br />

<br />

<br />

<br />

t<br />

S − T<br />

x<br />

t<br />

T − S<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

x


t<br />

ST<br />

x<br />

t<br />

T S<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

v(t, x, y) = (+∞, +∞).<br />

<br />

<br />

<br />

VM,I(t, x, y) = VM(t, x, t) ∩ I.<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

x


B <br />

<br />

<br />

<br />

<br />

A <br />

<br />

<br />

t<br />

A<br />

B<br />

<br />

A <br />

B <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

k Rk<br />

Rk(f)(x) = f(yx,k)<br />

x


{yx,i} x f(yx,i) ≤ f(yx,j) i < j <br />

k <br />

<br />

med(f)(x) = med<br />

y∈V (x) f(y)<br />

<br />

α ∈ [0, 1] <br />

k k <br />

k = ⌊α card V (x)⌋.


9 8<br />

7 4 3<br />

7 3 2<br />

<br />

<br />

α = 0.5 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

S1 S2 <br />

<br />

HOM(X) = {x|S1(x) ∈ X S2(x) ∩ X = ∅} <br />

<br />

S1(x) ∩ S2(x) = ∅ ∀x ∈ I,


S2,M S ′ 2,M <br />

S ′ 2,M(x) = S2,M(x) \ S1,M(x).<br />

<br />

S1,M S2,M S ′ 1,M S ′ 2,M <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

S1,M


t<br />

t<br />

(S1, S2)<br />

(S1,M , S ′ 2,M )<br />

t<br />

x<br />

x<br />

t<br />

t<br />

(S ′ 1,M , S ′ 2,M )<br />

?<br />

(S1,M , S2,M )<br />

(S ′ 1,M , S2,M )<br />

<br />

<br />

<br />

x<br />

x<br />

x


ε(X) = {x ∈ X|V (x) ⊂ X} <br />

ε(f)(x) = min f(y) <br />

y∈V (x)<br />

<br />

ˇ V (x) V (x) x <br />

δ(Y ) = <br />

y∈Y<br />

ˇV (y)<br />

δ(f)(y) = max<br />

x∈ ˇ f(x)<br />

V (y)<br />

<br />

<br />

<br />

Y ⊆ ε(X) ⇐⇒ δ(Y ) ⊆ X<br />

<br />

δ(Y ) = {y|V (y) ∩ Y = ∅} <br />

δ(f)(y) = max<br />

x|y∈V (x) f(x)


x ∈ I <br />

I ← +∞<br />

x ′ ∈ V (x) <br />

I(x) ← min(I(x), I(x ′ ))<br />

<br />

<br />

<br />

x ∈ I <br />

I ← −∞<br />

<br />

x ∈ I <br />

x ′ ∈ V (x) <br />

I(x ′ ) ← max(I(x ′ ), I(x))


x ∈ I <br />

I ← +∞<br />

x ′ ∈ V (x) <br />

I(x) ← min(I(x), I(x ′ ))<br />

<br />

<br />

<br />

x ∈ I <br />

I ← −∞<br />

<br />

x ∈ I <br />

x ′ ∈ V (x) <br />

I(x ′ ) ← max(I(x ′ ), I(x))


352 × 288 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

4 <br />

3


1/8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ <br />

Ĩ Ĩ <br />

<br />

<br />

f : I → J ˜ f :<br />

Ĩ → J


1/8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ <br />

Ĩ Ĩ <br />

<br />

<br />

f : I → J ˜ f :<br />

Ĩ → J


1/8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ <br />

Ĩ Ĩ <br />

<br />

<br />

f : I → J ˜ f :<br />

Ĩ → J


1/8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ <br />

Ĩ Ĩ <br />

<br />

<br />

f : I → J ˜ f :<br />

Ĩ → J


1/8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ <br />

Ĩ Ĩ <br />

<br />

<br />

f : I → J ˜ f :<br />

Ĩ → J


1/8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ <br />

Ĩ Ĩ <br />

<br />

<br />

f : I → J ˜ f :<br />

Ĩ → J


1/8 <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ <br />

Ĩ Ĩ <br />

<br />

<br />

f : I → J ˜ f :<br />

Ĩ → J


ϕ <br />

<br />

ϕ <br />

ϕ ψ <br />

<br />

ϕ ◦ ψ = ψ ◦ ϕ<br />

<br />

<br />

min max


h(i-1)<br />

h(i-2)<br />

h(i)<br />

h(i+1)<br />

h(i+2)


y I x Ĩ y ∈ V (x) <br />

<br />

Ĩ <br />

<br />

Ĩ <br />

mathcalI ˜ <br />

<br />

<br />

<br />

<br />

x <br />

<br />

<br />

<br />

<br />

<br />

g ≤ ε(f) ⇐⇒ δ(g) ≤ f <br />

<br />

S = {(0, 0), (−1, 0), (1, 0), ( 1<br />

2<br />

, 1)} <br />

, −1)}<br />

y x y ∈ V (x) ˇ S = {(0, 0), (1, 0), (−1, 0), (− 1<br />

2<br />

<br />

S <br />

<br />

<br />

g = ε(f) <br />

δ(ε(f)) ≤ f.<br />

<br />

I Ĩ <br />

<br />

ε(f)(x) = min ˜f(y)<br />

y∈V (x)<br />

δ(g)(x) = max<br />

x∈V (y) ˜g(y).<br />

<br />

Ĩ <br />

g ≤ ε( ˜ f) δ(˜g) ≤ f.


(0,0)<br />

S<br />

(0,0)<br />

ˇ S<br />

· · · · · · · ·<br />

· 3 3 6 3 3 3 ·<br />

· 4.5 6.5 8 5 2 2 ·<br />

· 3 3 4 4 3.5 4.5 ·<br />

ε(g)<br />

1 3 5 8 11 3 5 8<br />

3 3 6 7 9 3 4 7<br />

8 7 9 11 13 5 2 7<br />

9 6 3 4 9 8 10 5<br />

g<br />

· · · · · · · ·<br />

· · 6 7.25 6.5 3.5 · ·<br />

· · 8 8 8 8 · ·<br />

· · · · · · · ·<br />

δ1(ε(g))<br />

δ ˇ S ε S <br />

δ1(ε(g)) g <br />

· <br />

<br />

<br />

IV <br />

<br />

<br />

IV = <br />

V (x).<br />

x∈I<br />

g ≤ ε(f) ⇔ ∀x ∈ I, g(x) ≤ min ˜f(y)<br />

y∈V (x)<br />

⇔ δ(g) ≤ f<br />

g h<br />

<br />

∀x ∈ I, ∀y ∈ V (x), g(x) ≤ ˜ h(y), <br />

h <br />

h ≥ δ(g).<br />

˜ h<br />

h <br />

h <br />

δ(g)


h1 <br />

h2 min(h1, h2) <br />

h ≥ δ(g) <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

∧ε(hi) = ε(∧hi),<br />

<br />

<br />

<br />

<br />

<br />

<br />

I <br />

Ĩ Ĩ <br />

n <br />

<br />

1<br />

n <br />

n <br />

1<br />

n <br />

<br />

n <br />

<br />

Ĩ<br />

<br />

<br />

<br />

n <br />

<br />

<br />

<br />

<br />

<br />

<br />

<br />

n n 2


x2<br />

C1 : x1 ≥ 10 δ(g)<br />

h1(x1)<br />

h2(x1)<br />

C<br />

h2(x2)<br />

h1(x2)<br />

h1 h2 min(h1, h2)<br />

C2 : x2 ≥ 20<br />

<br />

x1<br />

x2<br />

C1 : x1 ≥ 10<br />

A<br />

C2 : x2 ≥ 20<br />

C3 : 1<br />

3 x1 + 2<br />

3 x2 ≥ 20<br />

<br />

<br />

<br />

C <br />

C <br />

h(xi) ≥ <br />

C1 C2 <br />

<br />

λixi ≥ A <br />

C1 C2 C3 <br />

h A <br />

h ≥ A ⇒ ε(h) ≥ g <br />

<br />

<br />

x1


n (xi)1≤i≤n n > 1 <br />

m m < n M +<br />

Z <br />

∀i ∈ [1, n] yi,Z = M +<br />

Z (xi) <br />

<br />

M +<br />

R yi,R xi <br />

<br />

yi,R <br />

(y j<br />

i,Z )j∈N<br />

(y0 i,Z ) <br />

n (ji,opt)1≤i≤n<br />

<br />

i1 = i2 ⇒ y ji 1 ,opt<br />

i1,Z = yji 2 ,opt<br />

i2,Z ∀(i1, i2)


y ji,opt<br />

i,Z<br />

= argmin <br />

1≤i≤n<br />

|y ji<br />

i,Z − yi,R| 2 . <br />

<br />

<br />

<br />

n (ji,opt)1≤i≤n <br />

<br />

n


y ji,opt<br />

i,Z<br />

= argmin <br />

1≤i≤n<br />

|y ji<br />

i,Z − yi,R| 2 . <br />

<br />

<br />

<br />

n (ji,opt)1≤i≤n <br />

<br />

n


k


k


N


N


ρ <br />

<br />

<br />

<br />

ρ


ρ <br />

<br />

<br />

<br />

ρ


ρ <br />

<br />

<br />

<br />

ρ

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!