实验二扇形束投影数据采集与滤波反投影重建实验 (2).doc_第1页
实验二扇形束投影数据采集与滤波反投影重建实验 (2).doc_第2页
实验二扇形束投影数据采集与滤波反投影重建实验 (2).doc_第3页
实验二扇形束投影数据采集与滤波反投影重建实验 (2).doc_第4页
全文预览已结束

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

南昌大学实验报告学生姓名: 胡文松 学 号: 6103413007 专业班级: 生医131班 实验类型:验证综合设计创新 实验日期: 20160516 实验成绩: 实验二、扇形束投影数据采集与滤波反投影重建实验一、实验目的及要求:1、 用MATLAB中的fanbeam函数获取Shepp-Logan头模型的扇形束投影数据;2、 显示扇形束投影数据并和平行束投影数据进行比较;3、 用MATLAB中的ifanbeam函数实现图像重建;4、 改变投影参数,观察对重建图像的影响。二、实验基本原理:fanbeam 函数,计算扇形束投影。句法F = fanbeam(I,D)F = fanbeam(.,param1,val1,param1,val2,.)F,sensor_positions,fan_rotation_angles = fanbeam(.)F = fanbeam(I,D) 计算图像I的扇形束投影 (正弦图) , D 是扇形束顶点到旋转中心的距离。F的每一列表示一个角度下的扇形采样数据。探测器圆弧状排布,探测器间相隔一度,旋转360度 0:359 。F = fanbeam(.,param1,val1,param1,val2,.) 指定参数控制扇形投影FanRotationIncrement 指定旋转角度增量,单位为度,默认值为1度。 FanSensorGeometry 指定探测器的排布方式 arc探测器圆弧状排布(默认值) line -探测器等距离线性排布FanSensorSpacing 指定探测器间的间隔,探测器圆弧状排布的情况,单位为度,默认值为1。探测器等距离线性排布,单位为像素。 F,sensor_positions,fan_rotation_angles = fanbeam(.) 返回探测器位置和旋转角度信息 ifanbeam 函数,由扇形束投影数据重建图像句法I = ifanbeam(F,D)I = ifambeam(.,param1,val1,param2,val2,.)I,H = ifanbeam(.)I = ifanbeam(F,D) 由投影数据F重建图像I。D是扇形束顶点到旋转中心的距离 I = ifanbeam(.,param1,val1,param2,val2,.) 指定参数控制图像重建。 FanCoverage 指定旋转角度范围 cycle旋转360度 0,360). minimal 重建图像所需的最小旋转角度。FanRotationIncrement 同fanbeamFanSensorGeometry 同fanbeamFanSensorSpacing 同fanbeamOutputSize 指定重建图像矩阵的大小,假如不指定,ifanbeam自动决定矩阵大小(扇形投影重建图像和原始图像矩阵大小可能不同)。三、主要仪器设备及实验耗材:具有XP或2000系统,并装有MATLAB系统的PC机。四、实验内容或步骤以及显示结果:1用MATLAB图像处理工具箱的phantom生成Shep-Logan头模型; P=phantom(256); imshow(P)2用MATLAB中的fanbeam函数获得头模型的扇形束和平行束投影数据;D = 200; F = fanbeam(P,D,FanRotationIncrement,0.5,FanSensorSpacing,0.5);theta=0:1:179;R4,xp=radon(P,theta);%平行束投影数据采集显示投影数据:figure.imshow(F,); title(扇形束投影数据) figure,imagesc(theta,xp,R4);xlabel(theta);ylabel(xprime);title(平行束180个角度投影数据);3用MATLAB中的ifanbeam函数对获得的投影数据进行滤波反投影重建,获得Shepp-Logan模型的重建图像。I = ifanbeam(F,D, FanSensorSpacing,0.5,FanRotationIncrement,0.5);I4=iradon(R4,1); %平行束投影数据重建显示重建图像:figure .imshow(I); title(扇形束投影数据重建) figure .imshow(I4); title (平行束180个角度滤波重建图像);五、实验小结通过实验结果可知,平行束通过180度扫描产生一个峰值,扇形束扫描360度出现两个峰值改变探测器的旋转角度或者排布方

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论