Discussion :

[assoumagh]

bonjour
je suis entrain d'implémenter un algorithme de flot optique donc j'ai besoin de calculer le gradient d'une image:seulement le gradient(derivée partielle par rapport de la luminance par rapport à x et à y et à t)mon problème c'est comment calculé la dérivée de la luminance par rapport à t.


merci pour votre aide

[Jerome Briot]

Plusieurs discussions ont déjà été ouvertes à ce sujet par l'utilisateur hanane78 ces derniers jours sur les forums MATLAB et Algorithmes...

[hanane78]

La dérivée de I par rapport au temps n'est autre que la différence entre l'image 1 et l'image 2

It= I1 - I2
et pour calculer la dérivée par rapport à x et y
vous pouvez utiliser cette code

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
 
[Ix, Iy]=gradient(double(I));

mais fait attention pour faire le suivi par le flot optique vous devez résoudre l'équation de flot optique. Et cette équation ne se résout pas de manière directe, il faut prendre en compte certaines contraintes

Je pense que j'ai été un peu claire

Bon Courage

[assoumagh]

bonjour,
merci oui j'ai déjà utilisé I1-I2
mais dans matlab est il possible de créer une matrice dont chaque élément est une matrice
c'est à dire par exemple en C on peut créer une liste chainée remplit par les image donc chaque élément de la liste est une image(matrice de luminance)
ce que je veux faire c'est de prendre une séquence d'image et non pas deux images seulement.
une autre question:tu parles des contraintes oui j'ai utilisé des contraintes je ne sais pas s'il y a d'autres enfaite j'ai trouvé que je dois calculer Ix et Iy et It puis la valeur moyenne de Vx et Vy puis une formule qui détermine la vitesse voir la pièce jointe.
si t'as déjà implémenter en matlab tu veux que je t'envoie mon code pour me le corriger j'ai pas trouvé de résultat

une dernière question:pouvez vous me donner deux images consécutives pour vérifier mon programme
merci d'avance

[assoumagh]

ci dessous ce j'ai pu faire avec matlab mais j'ai pas eu de résultat
pouvez vous m'aidez pour corriger ce programme meme s'il faut changer les entées et donner de nouvelles images faites le seulement envoyez moi ces images:
j'attends vos aides c'est trés urgent .....
merci d'avance

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
clear
lamda=2;n0=50;
Nfr=240; % numero du premier frame à convertir
 
nbfr=2; % nombre de frames à convertir
 
entree=fopen('foreman.qcif','r');
 
sortie=fopen('foreman230.dat','w');
 
fseek(entree,Nfr*99*384,0);
xprec=fread(entree,[176 144],'uchar');
for k_fr=1: nbfr
 
%k_fr
x=fread(entree,[176 144],'uchar');
yr=fread(entree,[88 72],'uchar');
yb=fread(entree,[88 72],'uchar');
[Ix,Iy]=gradient(x);
 
for i=1: 144
    for j=1: 176
        [Ix,Iy]=gradient(x);
        It = x - xprec ; 
        vx(i,j)= 0;
        vy(i,j)= 0;
    end
end
xprec=x;
end
n=1;
while(n<n0)
    for k_fr=1: nbfr
        for i=1: 144
            for j=1: 176
                vxb= 1/4*(vx(j-1,i)+vx(j+1,i)+vx(j,i-1)+vx(j,i+1));
                vyb= 1/4*(vy(j-1,i)+vy(j+1,i)+vy(j,i-1)+vy(j,i+1));
                alpha=lamda*(x*vxb+y*vyb+It)/(1+lamda*((Ix*Ix)+(Iy*Iy)));
                vx(i,j)=vxb-alpha*Ix;
                vy(i,j)=vyb-alpha*Iy;
                hold on;
                h=quiver(i,j,vx,vy);
                hold off;
            end
        end
    k_fr=k_fr+1;
    end
n=n+1;
end
 
fclose('all');

[hanane78]

la formule qui existe dans le fichier joint je ne comprends pas

ce que je connais c l'equation de flot optique voir l'image dans le fichier ci joint
cette equation n'est pas resoudre d'ne maniere simple parceque il contient de inconnu u et v .

pour mon application j'ai cherche un motif de reference dans la 1er image apres je fixé un point ds le 1 ere image et je cherche son deplacement ds le second image cette procedure fait pour tout les image de la sequence
chaque fois je trouve mon motif dans l'image suivi
donc le deplacement de ce pixel n'es autre que les inconnu u et v qu'on existe dans l'equation

voir des documentation sur Lukas-Kanade flow
quant a votre code j'arrive pas a le comprendre maintenant j'ai pas matlab sur ce Pc pour le tester