TP-traitement-audio-visuel/TP9/filtrage_sinogramme.m

17 lines
290 B
Mathematica
Raw Permalink Normal View History

2023-06-22 18:47:16 +00:00
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