医学成像技术滤波反投影
function reconstruction
clear;
close all;
N=256;%N*N大小图像
delta=3; %相当于Φ
a = 100;
b = 50;
d = 1;
%生成椭圆
for m =1:180/delta;
L=a^2*cos(m*delta)^2+b^2*sin(m*delta)^2;
for t=-N/2+1:N/2;
if (t^2 <= L)
proj(t+N/2,m)=2*a*b*sqrt(L-t^2)/L;
else
proj(t+N/2,m)=0;
end
end
end
%%%%各种滤波函数
for t = -N+1:N-1;
h1(t+N)=-2/(pi^2*d^2*(4*t^2-1)); % S-L
end;
for t=-N+1:N-1;
if t==0
h2(t+N)=1/(4*d^2); % R-L
elseif mod(t,2) == 1
h2(t+N)=-1/(t^2*pi^2*d^2);
else
h2(t+N)=0;
end;
end;
% Lewitt
esp=0.5;
for t=-N+1:N-1;
if t==0
h3(t+N)=(1-2*esp/3)/(4*d^2);
elseif mod(t,2)==1
h3(t+N)=-(1-esp)/(t^2*pi^2*d^2);
else
h3(t+N)=-esp/(t^2*pi^2*d^2);
end;
end;
%存储各种反投影结果
rProj1=zeros(N,N);
rProj2=zeros(N,N);
rProj3=zeros(N,N);
rProj4=zeros(N,N);
nProj=[];
.....................................................
相关文章
- 并发数据结构-1.1.2 阻塞技术
- 【华为云技术分享】基于Atlas 200 DK的原版YOLOv3(基于Darknet-53)实现(Python版本)
- 【华为云技术分享】云小课 | SAP容灾一点通
- 【华为云技术分享】华为云多元计算+AI 打造企业级智能数据湖
- 【6.19日云享专家技术公开课】从零开始快速搞定扫码付款
- 对一门技术的掌握程度
- Java后端整体架构技术栈及工具说明
- Atitit 2016 技术趋势与没落技术 目录 1.1. 离线优先Web应用程序(Offline first web applications)1 1.2. 依照产品而不是项目(product
- Unity技术-GameFramework文档系列(五)- 创建实体
- 【华为云技术分享】华为云多元计算+AI 打造企业级智能数据湖
- [SSD核心技术:FTL 5] 固态硬盘DRAM缓存技术详解
- 如何学习性能测试?盘点【性能测试】必须弄懂的5大技术点
- 数据技术初学者不用再东搜西找了,一篇文章解读20个最新数据技术概念
- 擎创技术流 | ClickHouse实用工具—ckman教程(13)
- 擎创技术流 | ClickHouse管理工具—ckman教程(2)ClickHouse集群