Discussion :

[matriax]

Bonjour,

J'ai écrit une partie d'un algorithme dans le but de tester si chaque pixel d'une image (fillimage20) fait partie d'un cercle ou pas.
1)Si le pixel est dans le cercle je souhaiterais le remplacer par la valeur d'un pixel d'une autre image (img02) (en utilisant les mêmes coordonnées i,j entre chaque image)
2)Si ce n'est pas le cas je voudrais garder la valeur du pixel de l'image originale (fillimage20).

Malheureusement, cet algorithme semble remplacer tous les pixels par la valeur des pixels de limage img02 (pour résumer fillimage20 devient img02). Je suppute que mon erreur est probablement liée à une mauvaise indexation de la boucle. En dépit de mes efforts je n'arrive point à discerner mon erreur.
Si quelqu'un à la gentillesse de m'aider, cela serait vraiment sympa.
Merci par avance.

Voici donc ma partie de l'algorithme ou j'ai un peu du mal. (Pour information, j'utilise le théorème de notre bon vieux Pythagore afin de comparer les distances entre le pixel i,j, le centre du cercle et le rayon du cercle):

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
fillimage20=logical(imfill(img800,'holes'));
[centersBrighti, radiiBrighti] = imfindcircles(fillimage,[Rcellmini Rcellmaxi],'ObjectPolarity','bright');
 
 
[L,C]=size(fillimage20);
Numberofcircle=size(radiiBrighti,1);
 
if Numberofcircle~=0;
    centersBright2=[centersBrighti, radiiBrighti];
    sortcentersBrighti=centersBright2(:,[1 2]); 
    sortradiiBrighti=centersBright2(:,3);
 
    for i=1:L
        for j=1:C
            for nan=1:Numberofcell; 
                if img800(i,j)==0 && (((i-sortcentersBrighti(nan,2)))^2)+((j-sortcentersBrighti(nan,1))^2)<=((sortradiiBrighti(nan,1))^2;
                    fillimage20(i,j)=img02(i,j);
                end
            end
        end 
    end 
end

[Jerome Briot]

Les parenthèses sont mal placées dans la condition suivante :

Code :

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

(((i-sortcentersBrighti(nan,2)))^2)+((j-sortcentersBrighti(nan,1))^2)<=((sortradiiBrighti(nan,1))^2

Après correction (et simplification) :

Code :

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

((i-sortcentersBrighti(nan,2))^2+(j-sortcentersBrighti(nan,1))^2)<=sortradiiBrighti(nan,1)^2

Évite d'utiliser nan comme nom de variable car c'est déjà la nom d'une fonction MATLAB.

Tu peux également utiliser une variable intermédiaire pour améliorer la lisibilité de ton code.

Voici une version corrigée :

Code :

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

1
2
3
4
5
6
7
8
9
10
    for i = 1:L
        for j = 1:C
            for n = 1:Numberofcell; 
                dij = (i-sortcentersBrighti(n,2))^2+(j-sortcentersBrighti(n,1))^2;
                if img800(i,j)==0 && dij<=sortradiiBrighti(n,1)^2
                    fillimage20(i,j) = img02(i,j);
                end
            end
        end 
    end

[matriax]

Merci énormément pour ton aide et les petites astuces.
Grace a tes conseils je vais pouvoir continuer mon projet.
Merci encore.
Un nouvel utilisateur conquis