function f=ispectremat(Ihat,Omega);

%          f=ispectremat(Ihat,Omega);
%
%Calcule, étant donnée une image  
%             Ihat(omega)=[:, omega=-Omega:Omega/K:Omega-Omega/K]  
%les valeurs de l'antécédent de Ihat (ligne par ligne) par la prise de
%spectre aux points t=k pi/Omega, k=-K:1:K-1. 

[M,MM]=size(Ihat);
MMM=fix(MM/2);
w=-ones(M,MM);
for k=1:2:MM
  w(:,k)=-w(:,k);
end
IIhat=Ihat';
UU=MM*ifft(IIhat);
U=UU';
U=w.*U;
U1=U(:,1:MMM);
U2=U(:,MMM+1:MM);
U=[U2,U1];
f=(Omega/(MM*pi))*U; 