⛄一、获取代码方式
获取代码方式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 备注
简介此部分摘自互联网,仅供参考,若侵权,联系删除