【图像重建】OSEM算法图像重建【含Matlab源码 3573期】

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【图像重建】基于matlab OSEM算法图像重建【含Matlab源码 3573期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab图像处理（初级版）

备注：
点击上面蓝色字体付费专栏Matlab图像处理（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab图像处理（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分代码

clc;
clear;
N=200; % 图像大小
N2=N^2;
disp(‘Generating phantom image’);
I=phantom(N);
Image=reshape(I’,N2,1); % 产生实时图像数据

theta=linspace(0,180,61);
theta=theta(1:60); % 投影角度

%% ========= 产生投影数据 ========= %%
disp(‘产生投影数据’);
P_num=185; % 探测器通道个数
P=Projection_Function(theta,N,P_num); % 产生投影数据
M=P_num*length(theta); % 投影射线的总条数
P=reshape(P,M,1);

%% ========= 获取系统矩阵 ========= %%
disp(‘获取系统矩阵’)
delta=1; % 网格大小
[W_ind,W_dat]=SystemMatrix(theta,N,P_num,delta);

%% ========= 划分子集 ========= %%
disp(‘划分子集’)
L=10; % 子集个数
T=length(theta)/L; % 每个子集包含的角度数
theta_seq=reshape(1:length(theta),L,T); % 产生角度编号矩阵5D
W_index=zeros(L,TP_num); % 包含L个子集的射线的行号
for i=1:L
temp=P_numtheta_seq(i,:);
for j=1:T
W_index(i,(j-1)P_num+1 :jP_num)=temp(j)-P_num+1:temp(j);
end
end

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1]姜虹.主成分分析的图像压缩与重构[J].电子设计工程. 2012,20(05)

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除