#+!

因为要模拟multi-pinhole collimator,涉及到计算projection的重叠问题。类似小孔成像，建模是下一步工作。假设single light source,clear allclcN=80;I64=zeros(N,N); %预定义平面164的灰度值为O[m,n]=meshgrid(linspace(-N/2,N/2-1,N));%确定坐标系及坐标原点的位置r1=15; %圆半径大小控制变量(单位pixel)a1=-10; %圆心位置控制变量b1=0; %圆心位置控制变量D1=((m+a

模拟的r射线，不用考虑散射了。function []=spheresource()h=input('source to collimator:');H=input('source to detector:');r=input('hole diameter:');d1=input('pinhole1 parameter(x y axis):');d2=input('pinhole2 parameter:');d3=input('pinhole3 parameter:');d4=input('pinhole4