投影交集overlapping仿真

因为要模拟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