Discussion :

[Clotho]

Salut tout le monde.

J'avance (tout) doucement dans la création d'un algorithme de traitement d'images, et me retrouve confronté à un problème que j'ai du mal à définir, et donc du mal à résoudre.

Pour la culture: je veux faire un traitement d'image sur une acquisition video en 640*480. Pour diminuer la charge de travail, j'ai fait une ROI circulaire centrée. Sur cette matrice, je veux appliquer une suite de masques (des cercles concentriques de dimension croissante).

J'ai réussi à:
modéliser mon masque "tout ou rien" de ma vidéo.
modéliser ma matrice de masquage, comportant 20 cercles.
appliquer la multiplication de l'une par l'autre, afin de trouver mes portions de masque à appliquer dans mon algo ultérieurement.

Ma matrice principale fait 480*640
mes matrices de masque font 201*201


voici le code (non fonctionnel, il manque les fichiers initiaux de masquage et de max détectés)

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
load masquage.mat %matrice de masquage
load coordonnees_max.mat %max détectés via un autre script
 
nb_max=size(maxi,1);
masquelong=mask.Length;
nb_images=masquelong*nb_max;
 
masque_roi=cell(nb_max,masquelong); %prealloc
champ_cam_vide=zeros(480,640); %prealloc
 
%% création de la ROI centrale.
 
champ_cam=zeros(480,640);
 
struc_centr=strel('disk',200,0);
roi_chem=getnhood(struc_centr); %créé une matrice de 401*401 avec un cercle ne comportant que des 1, et l'extérieur ne comportant que des 0
 
%on veut insérer un cadre de 401*401 dans une matrice de 480*640, et ce, de façon (+/-) centrée.
 
x_add=(480-400)/2;
y_add=(640-400)/2;
 
champ_cam(x_add:x_add+400,y_add:y_add+400)=roi_chem; %matrice de 640*480 avec un cercle rempli de 1 au milieu
 
clear roi_chem struc_chem x_add y_add %coup de balai
 
%% application du masque en fonction des maxima détectés.
for i=1:nb_max
    coord_x_max=max(i,1); %extraction des x et des y de chaque max détecté
    coord_y_max=max(i,2);
 
    for j=1:masquelong
        imagemasque=getdatasamples(mask,j); %lecture pas à pas du masque à appliquer
        masque_roi{i,j}=champ_cam(coord_x_max+100:coord_x_max+300,coord_y_max+100:coord_y_max+300).*imagemasque; %création d'un masque sur une ROI, mise en mémoire en tableau de cellules
 
 
    end
end
 
linearisation=reshape(masque_roi,1,nb_images);

Je bute dans la boucle "for" imbriquée: je voudrai intégrer chacune de mes valeurs de masque_roi{i,j} dans une matrice nulle de 480/640 mais je n'arrive pas à conceptualiser la chose au delà. Je n'arrive pas à faire d'indexation de type

Code :

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

roi_finale(x_add:x_add+200,y_add:y_add+200)=masque_roi{i,j}

et encore moins

Code :

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

roi_finale{i,j}(x_add:x_add+200,y_add:y_add+200)=masque_roi{i,j}[/

Je ne sais même pas comment chercher une réponse dans l'aide ou le forum a ce sujet.

Si vous arrivez à me comprendre (sait-on jamais ) est-ce faisable? si oui, comment?

[Jean Dumoncel]

Bonjour,

la ligne

Code :

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

roi_finale(x_add:x_add+200,y_add:y_add+200)=masque_roi{i,j}

renvoie une erreur? si oui laquelle?

Quant à la deuxième ligne je ne l'ai pas comprise, il manque un bout en fin de ligne, non?

Code :

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

roi_finale{i,j}(x_add:x_add+200,y_add:y_add+200)=masque_roi{i,j}[/

[Clotho]

Hello,

Alors en répondant du bas vers le haut: la dernière ligne avec [/ c'est une erreur au moment d'écrire mon message, je n'y ai pas prêté attention, c'est la fin de la balise code...

pour

Code :

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

roi_finale(x_add:x_add+200,y_add:y_add+200)=masque_roi{i,j}

il fallait en fait lire:

Code :

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

roi_finale(coord_x_max+200:coord_x_max+200,coord_y_max+200:coord_y_max+200)=masque_roi{i,j}

Les message d'erreur est différent selon que je fasse une préallocation de roi_finale en tant que tableau de cellules ou en tant que matrices nulle. De base il est:

??? Cell contents assignment to a non-cell array object.

Error in ==> roi_cheminee at 35
roi_finale{i,j}(coord_x_max:coord_x_max+200,coord_y_max:coord_y_max+200)=masque_roi{i,j};

Et si je déclare roi_finale en tant que cell array j'ai un résultat que je ne m'explique pas, en lieu et place d'une matrice de 480/640 j'ai ma matrice "masque_roi".

[Clotho]

Pour faire un exemple plus simple, dans le cas ou j'ai une matrice comme ceci:

Code :

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

1
2
A=zeros(5) 
B=[0 1;1 1]

Je veux que selon un vecteur C [ 1 2] (par exemple) ma matrice B soit introduite dans les valeurs de A afin d'avoir une matrice:

Code :

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

 final=[0 0 0 0 0; 0 0 1 0 0; 0 1 1 0 0;0 0 0 0 0; 0 0 0 0]

Ca sera plus joli dans matlab.

J'arrive a faire une opération "simple" comme ca, comme je fais dans la première partie de mon code mais n'arrive pas à l'introduide dans ma boucle, a partir de tableau de cellules et du vecteur C variable.

[Invité]

Avec ton exemple simple, cela donnerait quelque chose comme:

Code :

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

1
2
3
4
5
A=zeros(5);
B=[0 1;1 1];
[m n] = size(B);
i = 2; j = 2;
A(i:i+m-1,j:j+n-1) = B

[Clotho]

En fait d'avoir simplifié mon gloubiboulga infâme et en mixant ce que tu m'as donné et mon modèle initial (oui, je tiens à mes boucles for, n'ayant jamais fait d'informatique jusqu'à y a trois mois je vois pas trop comment les éviter) j'ai réussi a sortir un "squelette" fonctionnel que je m'en vais adapter de ce pas, la bave aux lèvres et la bière à la main!

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
A=zeros(5);
B_1=[0 1;1 1];
B_2=[0 0;1 0];
C=[1 2; 2 1];
D=cat(3,B_1,B_2);
E=cell(1,2);
 
[m n o]=size(D);
 
for i=1:2
 
    x=C(i,1);
    y=C(i,2);
 
    for j=1:o
        A(x:x+m-1,y:y+n-1)=D(:,:,j);
        E{i,j}=A;
    end
 
end

Merci!