基于泊松模型的EM-ML算法简单描述（MatLab代码）

基于泊松模型的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