-
大小: 4KB文件類型: .m金幣: 1下載: 0 次發(fā)布日期: 2021-01-03
- 語言: Matlab
- 標(biāo)簽:
資源簡介
關(guān)于使用FBP算法進(jìn)行圖像的重建算法的Matlab程序,對于大家很有幫助,是自己研究出來的!!
代碼片段和文件信息
%%
clc;
clear;
close?all;
tic;?%計時開始
imageN=128;???????%圖像大小imageN*imageN,即64*64
I=?phantom(imageN);%讀取Shepp-Logan頭像
[imageNimageN]=size(I);
subplot(221)%最后顯示的圖有2行2列,現(xiàn)在顯示的是第一個
imshow(I);
title(‘64*64原始圖像‘);%顯示原始圖像
%%
%4096*4096的矩陣
nx=imageN;???????%圖像的寬64
ny=imageN;???????%圖像的高64
na=128;??????????%投影角度64
nb=128;??????????%每個角度下的射線數(shù)64
mask=true(nxny);??%預(yù)先生成64*64的單位矩陣,加快速度
ray_pix?=?1;???????%?像素大小
x?=?[0:nx-1]?-?(nx-1)/2;
y?=?[0:ny-1]?-?(ny-1)/2;
[xy]?=?ndgrid(x?y);?%相對于中心的坐標(biāo)轉(zhuǎn)換
x?=?x(mask(:));
y?=?y(mask(:));
np?=?length(x); ?
angle?=?[0:na-1]‘/na?*?pi;
tau?=?cos(angle)?*?x‘?+?sin(angle)?*?y‘;
tau?=?tau?/?ray_pix;
tau?=?tau?+?(nb+1)/2; %?射線編號從1開始
ibl?=?floor(tau); ????%?左半部
%?floor?為向后取整1.6-1;??1.2-1;
val?=?1?-?(tau-ibl); %?左半部的加權(quán)
評論
共有 條評論