基于泊松模型的EM-ML算法:
详细推导过程可参照网络资源:
MLEM算法全过程推导_百度文库:https://wenku.baidu.com/view/be21f4520a4c2e3f5727a5e9856a561252d321f3.html
这里只做简单介绍,如有问题请指正。谢谢!
1、迭代公式
图片来源:https://wenku.baidu.com/view/74774f53dd36a32d7275813e.html
2、计算流程
(1)给定图像初值;
(2)对当前图像进行投影运算;
(3)将(2)中的投影数据与真实的图像投影数据进行比较(每个向量元素相除);
(4)将比值反投影到图像空间,获得修正因子;
(5)对图像进行修正;
(6)重复(2)-(5)的过程,直到满足停机准则。
3、MatLab代码
clc
clear
H = [1 0 1 0;0 1 0 1;1 1 0 0;0 1 1 0;0 0 1 1];
p = [2.1;3.2;2.5;1.5;2.2];
x = [1;1;1;1];%初始化
% EM-ML
for i = 1:1:100
p1 = H*x;%投影运算
ratio = p./p1;%投影数据与真实投影比较
den = diag(H'*ratio);%分子
num = sum(H);%分母
x = den*x./num';%图像修正
end