Bonjour,
j'ai implémenté l'algorithme ci dessous, mais son exécution ne donne pas toujours le même résultat .
Qui peut voir avec moi, s'il vous plaît, où se trouve le problème
voici l'algo:
Les étapes suivantes sont exécutées
successivement :
1. Nous appliquons l’analyse en composantes principales sur l’image de test en cours.
2. Nous appliquons l’algorithme de segmentation des « nuées dynamiques » sur les données générées par l’ACP. La segmentation est en deux classes (texte original du recto ou autre).
3. Nous projetons les deux classes générées dans l’espace des couleurs d’origine.
4. Nous comparons l’histogramme logarithmique de chaque classe image. L’histogramme ayant le plus grand nombre de pixels sombres est associé à la classe image qui doit être ensuite traitée. Cette image va subir tous les traitements précédents en commençant par l’étape (1).
Voici mon code
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
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
clc; 
close all; 
clear all
[x, pathname] = uigetfile( ...
    {
    '*.*','All Files (*.*)'}, ...
    'Import first image');
x = strcat(char(pathname),char(x));
img_originale=imread(x);
 
figure, imshow(img_originale),title('originale');
 
tic
img=img_originale;
 
for cmp=1:2
    n =size(img,1);
    m =size(img,2);
 
    R=[];V=[];B=[];
    for i=1: size(img,1)
        for j=1:size(img,2)
            r=img(i,j,1);
            v=img(i,j,2);
            b=img(i,j,3);
            R=[R;r];
            V=[V;v];
            B=[B;b];
        end
    end
 
    x=[R V B];
    x=double(x);
    covariance= cov(x);
 
    [U E] = eig(covariance); %E: Valeurs et U: vecteurs propres
    E = ([E(3,3),E(2,2),E(1,1)]);  % arrange the Eignvalue in a decending order
    U = [U(:,3),U(:,2),U(:,1)];  % arrange the Eignvector to make it compatable with Eignvalue.
 
    if cmp==1
        imgR= img(:,:,1);
        imgV= img(:,:,2);
        imgB= img(:,:,3);
    else
        imgR=[];%zeros(size(img,1),size(img,2),'uint8');
        imgV=[];%zeros(size(img,1),size(img,2),'uint8');
        imgB=[];%zeros(size(img,1),size(img,2),'uint8');
        for k=1:size(avant_plan,1)
            i= avant_plan(k,1);
            j= avant_plan(k,2);
            imgR = [imgR;img(i,j,1)] ;
            imgV = [imgV;img(i,j,2)] ;
            imgB = [imgB;img(i,j,3)] ;
        end
    end
 
    maximum=max(E);
    idx=[];
    if cmp == 1
        switch ( maximum )
            case E(1)
                idx=kmeans(double(imgR(:)),2);
            case E(2)
                idx=kmeans(double(imgV(:)),2);
            case E(3)
                idx=kmeans(double(imgB(:)),2);
        end
    else
        switch ( maximum )
            case E(1)
                idx=kmeans(double(imgR),2);
            case E(2)
                idx=kmeans(double(imgV),2);
            case E(3)
                idx=kmeans(double(imgB),2);
        end
    end
 
    if cmp ==1
        idxx=[];
        X=idx;
        %if cmp ==1
        for j=1:m
            var=(n-1)*(j-1)+j;
            colonne= X(var:var+(n-1),1);
            idxx=[idxx colonne];
        end
    else
        idxx=zeros(size(img,1),size(img,2));
 
        for i=1:size(avant_plan,1)
            idxx(avant_plan(i,1),avant_plan(i,2))=idx(i);
        end
    end
    %%projetons les deux classes générées dans l'espace des couleurs d'origine.
    c1=[];
    c2=[];
    [I1 J1] = find(idxx==1);
    c1 = [I1 J1];
    [I2 J2] = find(idxx==2);
    c2 = [I2 J2];
 
    classe1=zeros(n,m,'uint8');
    classe2=zeros(n,m,'uint8').*255;
 
    for k=1:size(c1,1)
        i= c1(k,1);
        j= c1(k,2);
        classe1(i,j,1) = img_originale(i,j,1) ;
        classe1(i,j,2) = img_originale(i,j,2) ;
        classe1(i,j,3) = img_originale(i,j,3) ;
    end
      % figure,imshow(classe1),title('classe 1');
 
    for k=1:size(c2,1)
        i= c2(k,1);
        j= c2(k,2);
        classe2(i,j,1) = img_originale(i,j,1) ;
        classe2(i,j,2) = img_originale(i,j,2) ;
        classe2(i,j,3) = img_originale(i,j,3) ;
    end
    %figure,imshow(classe2),title('classe 2');
 
    switch ( maximum )
        case E(1)
            [counts1,x1] = imhist(classe1(:,:,1));
            [counts2,x2] = imhist(classe2(:,:,1));
        case E(2)
            [counts1,x1] = imhist(classe1(:,:,2));
            [counts2,x2] = imhist(classe2(:,:,2));
        case E(3)
            [counts1,x1] = imhist(classe1(:,:,3));
            [counts2,x2] = imhist(classe2(:,:,3));
    end
 
    %figure
    y1=semilogy(x1,counts1);
 
    %figure
    y2=semilogy(x2,counts2);
 
    avant_plan=[];
    if max(y1,y2)==y1
        img= classe1;
        avant_plan=c1;
        cmp
        'Classe1'
 
    else
        img= classe2;
        avant_plan=c2;
        cmp
        'Classe2'
    end
 
end
 
img_finale = zeros(size(img_originale,1),size(img_originale,2),'uint8');
 
for i=1:size(avant_plan,1)
    x=avant_plan(i,1);
    y=avant_plan(i,2);
    img_finale(x,y)= 255;
end
figure,imshow(img_finale);
 
ind1=0;sum1=0;Moy1=0;
ind1=sum(img_finale(:)==0);
 
img1=img_originale(:,:,1);
img2=img_originale(:,:,2);
img3=img_originale(:,:,3);
 
sum1=sum(img1(img_finale== 0));Moy1= sum1/ind1;
sum2=sum(img2(img_finale== 0));Moy2= sum2/ind1;
sum3=sum(img3(img_finale== 0));Moy3= sum3/ind1;
 
finale=img_originale; %zeros(h,l,'uint8');
finale (:,:,3)=round(Moy3);
finale (:,:,2)=round(Moy2);
finale (:,:,1)=round(Moy1);
 
    for i=1:size(avant_plan,1)
        x=avant_plan(i,1);
        y=avant_plan(i,2);
        finale(x,y,1)=img_originale(x,y,1);
        finale(x,y,2)=img_originale(x,y,2);
        finale(x,y,3)=img_originale(x,y,3);
    end
 
figure,imshow(finale,[]),title('image finale');