TP-traitement-audio-visuel/TP9/filtrage_sinogramme.m
2023-06-22 20:47:16 +02:00

17 lines
290 B
Matlab

function s = filtrage_sinogramme(s)
fft_s = fft(s, [], 1);
fft_s = fftshift(fft_s, 1);
filtre = linspace(-1, 1, size(fft_s, 1))';
filtre = abs(filtre);
filtred = filtre .* fft_s;
filtred = ifftshift(filtred, 1);
s = ifft(filtred);
s = real(s);
end