因为要模拟multi-pinhole collimator,涉及到计算projection的重叠问题。
类似小孔成像,建模是下一步工作。
假设single light source,
clear all
clc
N=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+a1).^2+(n+b1).^2).^(1/2); %圆函数关系式
i=find(D1<=r1); %单下标寻址实现语句,find()函数足条件D<-r的像素点的单下标值
I64(i)=0.4; %按要求为像素点赋值
r2=15; %圆半径大小控制变量(单位pixel)
a2=10; %圆心位置控制变量
b2=0; %圆心位置控制变量
D2=((m+a2).^2+(n+b2).^2).^(1/2); %圆函数关系式
i=find(D2<=r2);%单下标寻址实现语句,find()函数足条件D<-r的像素点的单下标值
I64(i)=0.4+I64(i); %按要求为像素点赋值
J=mat2gray(I64,[0,1]);
imshow(J);
len=length(find(I64(i)==0.8))
overlap=len/(N*N)
>>
len =
157
overlap =
0.0245