Discussion :

[carpette]

Bonjour a tous,

voila, je debute sur matlab en traitement des images, et je souhaite tout d'abord me familiariser avec en faisant des petites manipulations.
Je souhaite avant tout balayer, avec une image de taille 31 par exemple, une autre image de taille beaucoup plus grande afin de faire une correlation.
Sur ma grande image, j'ai des points blanc (je me suis arrangé pour faire des images binaire), et avec mon autre image crée, je souhaite savoir le nombre de points. Voici mon programme:

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
52
clear all
close all
clc
 
%Déclaration des variables
taille=31;
val=(taille-1)/2;
nbadventis=0; %nb adventices initialisé a zero
 
%Matrice circulaire
SE=strel('disk',16,0);
A=zeros(taille,taille);
A(val+1,val+1)=1;
imc=imdilate(A,SE);
taille_imc=size(imc)
Y=taille_imc(2);
X=taille_imc(1);
figure
subplot(1,2,1)
imshow(A)
subplot(1,2,2)
imshow(imc)
 
%%Image a traiter, variable
img=imread('mahala_filtre.bmp');
taille_img=size(img)
y=taille_img(2);
x=taille_img(1);
img1=img(:,:,1);
 
 
%%Balayage et corrélation
res1=zeros(size(x),size(y));% matrice de correlation 
init=zeros(size(x),size(y));% matrice d'initialisation 
for i=1:size(x)-taille
    for j=1:size(Y):size(y)-taille
        res1(i,j)=corr2(imc,img1); 
        if res1>=0,60
            nbadventis=nbadventis+1
            init(i,j)=init(i,j)+res(i,j);
            i=i+size(X)
        else i=i+1;
 
        end
    end
end
 
figure
subplot(1,2,1)
imshow(img1)
subplot(1,2,2)
imshow(init)

Voila, le probléme viens, je pense de mon balayage, et je me retrouve coincé.
Merci d'avance

PS: je joins la photo si vous voulez compiler et regarder mes resultats

[Invité]

Bonjour,

Déjà je repère un problème majeur: i est la variable sur laquelle tu fais ta boucle for, or tu la changes à l'intérieur: ceci n'est pas possible, l'incrémentation se fait automatiquement.
Petit exemple:

Code :

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

1
2
3
4
5
for i = 1:10
   disp(['avant: ' num2str(i)])
   i= 1;
   disp(['après: ' num2str(i)])
end

Si tu veux changer i à l'intérieur, pense à utiliser plutôt une boucle while

PS: Pense à l'indentation de ton code qui le rend plus visible: Ctrl+A , Ctrl+I

[Invité]

Ensuite,

taille_imc=size(imc)
Y=taille_imc(2);
X=taille_imc(1);

et

taille_img=size(img)
y=taille_img(2);
x=taille_img(1);

Peuvent être remplacé par:

Code :

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

1
2
[X Y] = size(imc);
[x y pasbesoin]=size(img);

Et par la suite, x, y, X, Y étant des valeurs, pourquoi mettre size(x), size(y)?

res1=zeros(size(x),size(y));% matrice de correlation
init=zeros(size(x),size(y));% matrice d'initialisation

Remplacé par:

Code :

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

1
2
res1=zeros(x,y);% matrice de correlation 
init=zeros(x,y);% matrice d'initialisation

[carpette]

Merci pour vos réponses. Alors, pour ma boucle, je souhaite avant tout de copier le point dans ma nouvelle matrice init, et le compter dans la variable nbadventis et de l'afficher. Cependant je ne vois pas comment faire, car c'est la façon de manipuler la matrice circulaire qui me pose un probléme, au niveau du balayage.