一、理论基础

计算机层析成像技术(CT)是近十几年发展起来的一种新的非接触无损检测技术，它具有检测精度高、重建图像无影像重叠、空间分辨率和密度分辨率高、可以直接进行数字化处理等优点,现已被广泛应用于航空、航天、机械、公安、海关、医疗等诸多领域。滤波反投影算法是目前比较常用的CT图像重建算法，它速度快，图像质量好。但在二维平面上，Radon变换不具有局部性。于是，寻找一种既能减少辐射剂量，又能重建感兴趣区域的局部图像重建算法，引起了人们的极大兴趣，这就是所说的局部CT。最后利用MATLAB对滤波反投影算，多尺度全局重建和多尺度局部重建进行了仿真，从而验证算的正确性。

其中O(x，y)代表原始图像，R(x，y)代表重建图像。当然孤立地去看相对误差并不是一个好的衡量标准，它必须同人们的视觉评估结合起来进行考虑。

用S-L头模型进行计算机仿真研究的主要好处之一是可以获得该模型投影数据的解析表达式。一个中心位置在原点且未经旋转的椭圆，其长轴与X轴重合，短轴与Y轴重合。假设椭圆内的密度为r、椭圆外密度为零，则该椭圆图像可用以下方程表示：

在研究从投影重建图像的算法时，为了比较客观的评价各种重建算法的有效性，人们常选用公认的Sheep-Logan 头模型作为研究对象。该模型由10个位置、大小、方向、密度各异的椭圆组成，象征一个脑断层图像。图3-2是S-L头模型的灰度显示图像。

二、核心程序

clc;
clear;
close all;I=phantom(256);  %生产头部模型图
figure(1);
imshow(I);       %显示图像
IMG=double(I);   %双精度显示
%=====step1================================
THETA=0:179;
PR=touyin(IMG,THETA) %进行投影
figure(2);
imshow(PR,[]);
%=====step2================================
filtPR=lvbo(PR);%进行滤波
figure(3);
imshow(filtPR,[]);
%=====step3================================
BPI=fantouyin(filtPR,THETA);%反投影
figure(4)
imshow(BPI,[]);imwrite(BPI,'BPI.bmp');

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:CountIMG = fft(PR(:,i));fimg = IMG.*f';g(:,i) = ifft(fimg);
end
filtPR= real(g);

clc;
clear;
close all;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)
imshow(cod_a,map);subplot(222)
imshow(cod_h,map);subplot(223)
imshow(cod_v,map);subplot(224)
imshow(cod_d,map);Z1=idwt2(cod_a,cod_h,cod_v,cod_d,'db1');figure(3);
imshow(Z1,map);for i=1:256for j=1:256error(i,j)=(I(i,j)-Z1(i,j))^2/I(i,j)^2;end
endfor i=1:256for j=1:256if I(i,j)==0error(i,j)=255;    endend
endfor i=1:256for j=1:256if error(i,j)>255;error(i,j)=255;    endend
endfigure(4);
imshow(error,map);

三、仿真测试结果

基于shepp-Logan模型的多尺度局部重建的MATLAB仿真与分析

其仿真图像为：

图1 原始仿真图（左）和局部数据仿真图（右）

其小波系数仿真如下所示：

图2 小波系数

最后由局部数据重建的图像如下所示：

图3 局部数据重建效果

基于shepp-Logan模型的多尺度全局重建的MATLAB仿真与分析

我们首先用一层小波变换得到分解后的小波系数。仿真结果如下所示：

图4 小波系数

然后通过二维反小波变化，可以得到如下仿真结果：

图5 小波重建结果

A09-02

基于Radon滤波反投影算法的CT图像重建matlab仿真相关推荐

【ISAR成像定标方法（3）—基于SGP4模型的空间目标定标方法MATLAB仿真】
目录前提介绍基于SGP4模型的转速估计基于SGP4模型的空间目标定标仿真实验结语前提介绍本章内容简介:本文研究了使用双行轨道报和SGP4模型估计空间LEO目标位置信息,并根据几何关系推测目 ...
matlab三相电源含谐波,谐波电能计量算法的研究及MATLAB仿真-毕业论文（设计）.doc...
毕业设计(论文) 题目: 谐波电能计量算法的研究及MATLAB 仿真学院: 机电工程学院专业: 电子信息工程班级: 电信学号: 姓名: 指导教师: 教务处制 2011 年 5 ...
【图像重建】基于布雷格曼迭代（bregman alteration）算法集合ART算法实现CT图像重建附matlab代码
1 简介 Fluorescence diffuse optical tomography (fDOT) is a noninvasive imaging technique that makes it ...
【回归预测】基于粒子滤波实现锂离子电池寿命预测附matlab代码
1 内容介绍随着现代生产生活对系统设备可靠性.安全性要求的提高,从成本.可靠性的角度考虑,电子系统正逐步由原来的定期维修变成视情维修(CBM, Condition Based Maintenanc ...
图像相减的matlab仿真及光栅滤波法,图像相减的MATLAB 仿真及光栅滤波法实验实现...
图像相减的MATLAB 仿真及光栅滤波法实验实现毕业设计(论文)中期报告题目图像相减的MATLAB仿真及光栅滤波法实验实现院(系)光电学院专业光信息科学与技术班级090106姓名陈凤学号090106 ...
基于直方图均衡化的水下激光图像处理（MATLAB仿真）
原文:http://blog.csdn.net/sunlinyi66/article/details/39215447 基于直方图均衡化的水下图像增强(MATLAB仿真) 可以对RGB三个通道进行直方 ...
【散射点模型以及雅克42飞机实测数据RD算法ISAR成像的MATLAB仿真实验】
基于散射点模型进行仿真实验分析,探讨距离压缩效果和RD成像算法的有效性. 如图2-6 a)所示为模拟空间站目标的散射点模型,该模型由模拟的空间站四块大太阳能板以及实验核心舱组成.仿真无噪声环境下的实验 ...
图像相减的matlab仿真及光栅滤波法,图像相减的MATLAB 仿真及光栅滤波法实验实现.doc...
您所在位置:网站首页 > 海量文档 &nbsp>&nbsp计算机&nbsp>&nbsp图形图像图像相减的MATLAB 仿真及光栅滤波法实验实现.do ...
基于LMS自适应滤波器的QPSK信号均衡器matlab仿真
目录 1.算法概述 2.仿真效果 3.MATLAB仿真源码 1.算法概述在自适应滤波器设计中,最小均方(Least Mean Square,LMS)算法使用随机梯度下降的方法实现代价函数的最小化,具 ...
matlab电流检测,基于d-q变换的谐波电流检测方法与MATLAB仿真验证
描述引言目前,谐波抑制的一个重要趋势是采用有源电力滤波器(Active Power Filter-APF).而该滤波器性能的好坏与它所采用的谐波电流检测方法有很大关系.因此,如何实时准确地检测出非 ...

基于Radon滤波反投影算法的CT图像重建matlab仿真

一、理论基础

二、核心程序

三、仿真测试结果

基于Radon滤波反投影算法的CT图像重建matlab仿真相关推荐

最新文章

热门文章

基于Radon滤波反投影算法的CT图像重建matlab仿真

一、理论基础

二、​​​​​​​核心程序

三、仿真测试结果

基于Radon滤波反投影算法的CT图像重建matlab仿真相关推荐

最新文章

热门文章

二、核心程序