利用投影演算法來計算系統矩陣左乘和右乘

採用distance driven投影模型

mode 1: 左乘得影象向量

% 引數設定

para.ht_orig = 0; %切片初始高度，即第一（下層）切片與探測器的間距，單位是切片個數。
para.pv_ind = 14; %投影角序號，如1、3、14等。
para.nvx = 1080; %x軸方向體素個數1080
para.nvy = 2304; %y軸方向體素個數2304
para.nvz = 1;  %z軸方向體素個數
% mode = 1; % 右乘
% img_mat = ones(para.nvx*para.nvy*para.nvz,1);
mode = 2; % 左乘
img_mat = ones(1,2580480);
[sps_sys_mat,myproj_vectori] = voxel_based_ddp2_fun2(img_mat,para,mode);

>> vector2 = ones(1,2580480);

>> rslt2 = vector2*sps_sys_mat;
>> isequal(rslt2,myproj_vectori)
ans =
1

mode2： 右乘得投影向量

>> vector1 = ones(2488320,1);

>> rslt1 = sps_sys_mat*vector1;
>> isequal(rslt1,myproj_vectori)
ans =
0

0這種差異，我實在不知道是怎麼造成的。

利用投影演算法來計算系統矩陣左乘和右乘

採用distance driven投影模型 mode 1 左乘得影象向量 引數設定 para.ht orig 0 切片初始高度，即第一 下層 切片與探測器的間距，單位是切片個數。para.pv ind 14 投影角序號，如1 3 14等。para.nvx 1080 x軸方向體素個數1080 para...

計算矩陣演算法

to change this template,choose tools templates and open the template in the editor.package testmutlchat author administrator public class computematri...

基於Siddon演算法的CT正投影計算

ct技術在醫學診斷 工業無損檢測 安檢及地質勘測等方面發揮著重要的作用。ct技術的數理基礎是radon變換，在實際的ct系統中，radon正變換稱為投影資料的採集，radon逆變換稱為切片影象重建。ct投影 在優化radon變換的投影成像過程，測試radon逆變換的ct重建演算法等方面具有十分重要的...