1.完整项目描述和程序获取
>面包多安全交易平台:https://mbd.pub/o/bread/Y56XmZ5r
>如果链接失效，可以直接打开本站店铺搜索相关店铺：
>如果链接失效，程序调试报错或者项目合作也可以加微信或者QQ联系。
2.部分仿真图预览
3.算法概述
计算机层析成像技术（Computed Tomography，简称CT）是利用具有一定能量的射线源（X射线，γ射线）对物体进行断层扫描，并根据物体外部的探测器获得的物理量（指物质对射线的衰减系数）生成的一维投影数据，通过特定的重建算法，得到所扫描断层的二维图像。CT断层图像具有无影像重叠、空间和密度分辨率高、可直接进行数字化处理等优点，通过近十几年的发展已成为非接触无损检测的主流技术，是关键部件检测、机械仿型设计、安全检查等方面强有力的手段，并广泛应用于航空、航天、机械、汽车、船舶、公安等领域。此外，在医学领域，医用CT机已成为疾病诊断的重要辅助工具，对病灶的位置及病变程度的良好再现已使其成为医学影像领域里的主要医疗器械。CT机有三个组成部分。一是数据采集系统，二是图像重建（主要是图像重建算法），三是图像后处理及显示系统。CT图像重建算法主要分为两类，即变换法和级数展开法。变换法包括滤波反投影算法、傅立叶变换法以及ρ滤波法（rho-filtered layergrams），小波变换法；级数展开法主要是代数重建算法。滤波反投影是比较常用的变换算法，它具有速度快，空间和密度分辨率好的优点。
4.部分源码
function filtPR=lvbo(PR);
n = size(PR,1);
sideSize = n;
a = 1;
[Length, Count] = size(PR);
w = [-pi:(2*pi)/Length:pi-(2*pi)/Length];
rn1 = abs(2/a*sin(a.*w./2));
rn2 = sin(a.*w./2);
rd = (a*w)./2;
r = rn1*(rn2/rd)^2;
f = fftshift(r);
for i = 1:Count
IMG = fft(PR(:,i));
fimg = IMG.*f';
g(:,i) = ifft(fimg);%利用FFT,IFFT进行滤波
end
filtPR= real(g);
............................................................
I=phantom(256); %生产头部模型图
figure(1);
imshow(I); %显示图像
IMG=double(I); %双精度显示
[cod_a,cod_h,cod_v,cod_d,map]=wtest(IMG);%一层小波变换
figure(2);
subplot(221)%显示分解系数1
imshow(cod_a,map);
subplot(222)%显示分解系数2
imshow(cod_h,map);
subplot(223)%显示分解系数3
imshow(cod_v,map);
subplot(224)%显示分解系数4
imshow(cod_d,map);
Z1=idwt2(cod_a,cod_h,cod_v,cod_d,'db1');%重建
figure(3);
imshow(Z1,map);
for i=1:256
for j=1:256
error(i,j)=(I(i,j)-Z1(i,j))^2/I(i,j)^2;%误差求解
end
end
for i=1:256
for j=1:256
if I(i,j)==0
error(i,j)=255;
end
end
end
for i=1:256
for j=1:256
if error(i,j)>255;
error(i,j)=255;
end
end
end
09_002_m
