资源描述：
《ct图像三维重建(附源码)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、结束计算重建图像和原图的性噪比，并进行输出用函数ifanbeam根据扇束投影数据重建图像，并展示用函数fanbeam进行映射，得到扇束的数据，并展示对图片信息进行预处理，并进行展示输入图片数字选择生成128的图片信息程序流图：MATLAB源码：clc;clearall;closeall;%loadmri%载入mri数据，是matlab自带库%ph=squeeze(D);%压缩载入的数据D，并赋值给phph=phantom3d(128);prompt={'EnterthePiecenum(1to128):'};%提示信息“输入1到27的片的数字”name='Inpu
2、tnumber';%弹出框名称defaultanswer={'1'};%默认数字numInput=inputdlg(prompt,name,1,defaultanswer)%弹出框，并得到用户的输入信息P=squeeze(ph(:,:,str2num(cell2mat(numInput))));%将用户的输入信息转换成数字，并从ph中得到相应的片信息Pimshow(P)%展示图片PD=250;%将D赋值为250,是从扇束顶点到旋转中心的像素距离。dsensor1=2;%正实数指定扇束传感器的间距2F1=fanbeam(P,D,'FanSensorSpacing',
3、dsensor1);%通过P，D等计算扇束的数据值dsensor2=1;%正实数指定扇束传感器的间距1F2=fanbeam(P,D,'FanSensorSpacing',dsensor2);%通过P，D等计算扇束的数据值dsensor3=0.25%正实数指定扇束传感器的间距0.25[F3,sensor_pos3,fan_rot_angles3]=fanbeam(P,D,...'FanSensorSpacing',dsensor3);%通过P，D等计算扇束的数据值，并得到扇束传感器的位置sensor_pos3和旋转角度fan_rot_angles3figure,%创
4、建窗口imagesc(fan_rot_angles3,sensor_pos3,F3)%根据计算出的位置和角度展示F3的图片colormap(hot);%设置色图为hotcolorbar;%显示色栏xlabel('FanRotationAngle(degrees)')%定义x坐标轴ylabel('FanSensorPosition(degrees)')%定义y坐标轴output_size=max(size(P));%得到P维数的最大值，并赋值给output_sizeIfan1=ifanbeam(F1,D,...'FanSensorSpacing',dsensor1,
5、'OutputSize',output_size);%根据扇束投影数据F1及D等数据重建图像figure,imshow(Ifan1)%创建窗口，并展示图片Ifan1title('图一');disp('图一和原图的性噪比为：');result=psnr1(Ifan1,P);Ifan2=ifanbeam(F2,D,...'FanSensorSpacing',dsensor2,'OutputSize',output_size);%根据扇束投影数据F2及D等数据重建图像figure,imshow(Ifan2)%创建窗口，并展示图片Ifan2disp('图二和原图的性噪比为
6、：');result=psnr1(Ifan2,P);title('图二');Ifan3=ifanbeam(F3,D,...'FanSensorSpacing',dsensor3,'OutputSize',output_size);%根据扇束投影数据F3及D等数据重建图像figure,imshow(Ifan3)%创建窗口，并展示图片Ifan3title('图三');disp('图三和原图的性噪比为：');result=psnr1(Ifan3,P);function[p,ellipse]=phantom3d(varargin)%PHANTOM3DThree-dimen
7、sionalanalogueofMATLABShepp-Loganphantom%P=PHANTOM3D(DEF,N)generatesa3Dheadphantomthatcan%beusedtotest3-Dreconstructionalgorithms.%%DEFisastringthatspecifiesthetypeofheadphantomtogenerate.%Validvaluesare:%%'Shepp-Logan'Atestimageusedwidelybyresearchersin%tomography%'ModifiedShepp-Loga
8、n'(de
