Как правильно добавить пуассоновский шум на КТ-синограмму в Matlab.

Question

theta = linspace(0,360,900);
img = im2double(phantom(512));
figure; imshow(img)
c = radon(img, theta);
c = imnoise(c, 'poisson');
figure;imshow(c,[]);
c = iradon(c, theta);
figure; imshow(c);

Я хочу добавить шум Пуассона к фантомному изображению. Я делаю радоновое преобразование на изображении, а затем добавляю шум к синограмме с помощью imnoise (), а затем делаю обратное радоновое преобразование. Но я не могу получить правильную шумную синограмму с imnoise ().

Answer 1

есть некоторый коэффициент масштабирования при использовании imnoise со статистикой Пуассона. Согласно документации Matlab:

"Если I имеет двойную точность, то значения входного пикселя интерпретируются как средства распределений Пуассона, увеличенных на 1e12. Например, если входной пиксель имеет значение 5.5e-12, то соответствующий выходной пиксель будет быть сгенерированы из распределения Пуассона со средним значением 5,5, а затем уменьшены до 1e12. "

попробуйте вместо

 c = 1e12.*imnoise(c./1e12, 'poisson');