reconstruction d'une image

[mawkli]

Bonjour,
j'ai écrit un code Matlab au but de reconstruire une image à partir de sa transformée de fourier,on fixe l'amplitude et on initialise la phase puis on la varie en posant des contraintes ,je vais mettre mon code mais il manque quelque chose au niveau des contraintes parce que,comme vous le verrez, on obtient pas un bon résultat

merci d'avance

Code :

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
A = imread('barbara.png');
A=A(1:50,1:50); 
F = fft2(A);
M = abs(F);
P = angle(F);

%M donné 
[Ma Na]=size(A);
L=10;
R=zeros(size(A)); R(Ma/2+1-L:Ma/2+1+L,Na/2+1-L:Na/2+1+L)=ones(2*L+1,2*L+1);  
R=fftshift(R);

iM=real(ifft2(exp(j*P)));
Ik1=iM;


for k=0:150

    
V=find(Ik1>255); Ik1(V)=255;
V=find(Ik1<0); Ik1(V)=128;
Ik1=Ik1/max(max(Ik1))*255; 
If_k=fft2(Ik1);
M_k1=abs(If_k); 
If_k1=M_k1.*exp(j*P); 
Ik1old=Ik1; 
Ik1=(real(ifft2(If_k1)));

%Ik1_filtered=imfilter(real(ifft2(If_k1)),1/4*ones(2,2));
%M_k1=abs(fft2(Ik1_filtered)); 
%If_k1=M_k1.*exp(j*P); 
%Ik1=(real(ifft2(If_k1)));
end

Ik1=imfilter(real(ifft2(If_k1)),1/25*ones(5,5));
figure(1);
subplot(211); imagesc(A);colormap('gray');
subplot(212); imagesc(Ik1);colormap('gray');
figure(2); imagesc(Ik1old-Ik1); colormap('gray');

[mawkli]

aucune réponse??