- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
计算机图像处理技术结课(报告)
共9页第 PAGE 11页
卷积反投影图像重建
1 反投影重建基本介绍 [1]
设待重建图像为,它的二维傅氏变换为。根据中心切片定理,可通过在不同视角下的投影的一维傅氏变换求得。即:
待建图像:
(1.1)
因为,所以有:
同时:
(1.2)
先来看该式的第二个积分:
(1.3)
式中:
(1.4)
式(3.10)的物理意义是投影经过传递函数为的滤波器后得到的修正后的投影在满足时的值。将(3.11)代入(3.8),得到:
(1.5)
称为滤波反投影方程,其物理意义是经过给定点的所有滤波后的投影在~范围内的累加—反投影重建,得出点的像素值。
可见,滤波(卷积)反投影算法的具体包含三大步:
(1) 把在固定视角下测得的投影经过滤波,得到滤波后的投影;
(2) 对每一个,把反投影于满足的射线上的所有各点;
(3) 将步骤(2)中的反投影值对所有进行累加(积分),得到重建后的图像。
2 重建流程
2.1 首先我们利用phantom()函数产生一个头部幻影图像,用以检测二维重建算法,代码如下:
I=phantom(256);
subplot(2,2,1)
imshow(I,[]);
title(256*256原始图像);
效果图如图1所示, 图1
为一个大椭圆和几个小椭圆。
2.2 初始参数设置
重建采用的是平移加旋转的扫描方式,射线源在某一角度下水平移动,将物体全部照射后旋转一角度,如此重复,在这个过程中探测器相应地运动以接收X射线。根据此原理,将重建程序的初始参数设置如下:
[N,N]=size(I);
z=2*ceil(norm(size(I)-floor((size(I)-1)/2)-1))+3;% radon变换默认平移点数/角度
Nt=360; % 角度采样点数
Nd=N; % 平移数
x=pi/180; % 角度增量
d=N/Nd; % 平移步长
theta = 1:Nt;
a=zeros(N);
2.3 产生无噪声投影数据
[R,xp] = radon(I,theta);
e=floor((z-Nd)/2)+2;
R=R(e:(Nd+e-1),:);
R1=reshape(R,256,360);
radon(I,theta)产生I投影,默认z点/角度,即使指定N点也是z点.所以为避免重建图像放大或缩小,下面计算取投影时需补偿,补偿量e,如对256的图像,补偿为55,即pm的第55个点作为计算用的第一个投影。
2.4 添加噪声并将所有噪声平行投影进行显示
[mm,nn]=size(R1);
di=lognrnd(0,0.15,mm,nn);
R1= 10*(R1-min(R1(:)))/( max(R1(:))-min(R1(:)));
I0 = 1.5e5; % incident photons; decrease this for simulating low dose scans
rand(state, 0), randn(state, 0);
yi= poissrnd(I0 * di.*exp(-R1))+3*randn(size(R1));
if any(yi(:) == 0)
warn(%d of %d values are 0 in sinogram!, ...
sum(yi(:)==0), length(yi(:)));
end
R1 = log(I0 ./ max(yi,0.01)); % noisy sinogram
R1=max(R1,0);
% 显示
ff=2;
uu=22000;
v=ff*exp(R1/uu);
subplot(2,2,2)
imagesc(R1);
title(256*360有噪声平行投影);
colormap(gray)
colorbar
Q=reshape(R1,256,360);
效果图如图-2: 图-2
2.5 滤波器的选择与设计
最基本的从投影重建图像的滤波器:1971年提出的R-L重建滤波器(下图中实线表示)。
图-3 R-L 滤波器示意图
空域表达式为:
(2.1)
其中B为截至频率。若,d为空间采样间隔,可解出离散的滤波器空间脉冲响应:
(2.2)
MATLAB 代码如下:
g=-(Nd/2-1):(Nd/2);
for i=1:256
if g(i)==0
hl(i)=1/(4*d^2);
您可能关注的文档
- 东师《英语阅读(一)16秋在线作业2.doc
- 辽宁鞍山玉佛苑英文导游词.doc
- QFD(质量功能配置).doc
- 酒店加盟品牌优势对比一览表(图表).docx
- 大学体验英语综合教程第三版第三单元1.ppt
- 寒假班中级口语.doc
- 柴移英文说明书.doc
- 外贸英语:海关检查.docx
- 双语:从二次元到三次元-你懂弹幕吗.docx
- 【斯芬克艺术职业规划】艺术留学职业规划-新时代下会有哪些新鲜的设计角色出现.doc
- GB/T 34877.4-2024工业风机 标准实验室条件下风机声功率级的测定 第4部分:声强法.pdf
- 中国国家标准 GB/T 34877.4-2024工业风机 标准实验室条件下风机声功率级的测定 第4部分:声强法.pdf
- 中国国家标准 GB/T 5578-2024固定式发电用汽轮机规范.pdf
- GB/T 5578-2024固定式发电用汽轮机规范.pdf
- 《GB/T 5578-2024固定式发电用汽轮机规范》.pdf
- 《GB/T 4340.1-2024金属材料 维氏硬度试验 第1部分:试验方法》.pdf
- GB/T 4340.1-2024金属材料 维氏硬度试验 第1部分:试验方法.pdf
- 中国国家标准 GB/T 4340.1-2024金属材料 维氏硬度试验 第1部分:试验方法.pdf
- GB/T 43995-2024数字航天摄影测量 空中三角测量规范.pdf
- 中国国家标准 GB/T 43995-2024数字航天摄影测量 空中三角测量规范.pdf
文档评论(0)