Отфильтрованный BackProjection Matlab

Я уже реализовал код для BackProjection, и он работает. Теперь я должен реализовать код для отфильтрованной версии, используя сначала фильтр ramlak, а затем окно Хэмминга. Моя синограмма имеет размер 144х180. Моя идея

filter=zeros(1,144)';

filter(1:73)=[0:1/72:1];

filter(74:end)=[1-1/72:-1/72:1/72];

чтобы иметь такой фильтр введите описание изображения здесь

а затем для каждой проекции с

filtered_proj=fft(s).*filter;
proj=real(ifft(filtered_proj));

а затем выполнить алгоритм обратного проецирования

Не могли бы вы сказать мне, если это правильно?

Спасибо

0 ответов

Другие вопросы по тегам