MATLAB，3x3mm的光感接收器上，40X40的房间，每个房间（细胞）75X75，有效接收区域（活跃区域）65X65。已知光子落点坐标（x，y）有9000个。

Question

三个要求
1号:不管活性区域,进入每个细胞的光子数。（落入75x75的光子数）

2号:进入每个细胞活性区域的光子数。（落入65X65的光子数）

3号:PD适用值(0.563)来获取并选择值

最后,为了做3号,要经过1号和2号。
下面的代码已经完成三号，我的问题是一号和二号要求代码怎么加入下面的代码里。我自己也写了一下不过我刚开始学，希望可以给我一个参考。
这里是示意图：

%%光子消减
clear all；clc；
％ 重置
数据 = cell2mat(readcell("COUNTS1"));
%读取数据
all_length = 3；gapd_length = 0.025；active_length = 0.015；inactive_length = 0.005；
% 单元格大小和活动区域设置
gapd_active_count = zeros([all_length/gapd_length all_length/gapd_length]);
%根据单元格大小创建数组
对于我 = 1：大小（数据，1）
    % 重复进行到输PHOTON 的位置
     temp_x = floor(data(i,1)/gapd_length);temp_y = floor(data(i,2)/gapd_length);
     %用PHOTON 检测细胞位置
     如果 ((data(i,1)-temp_x*gapd_length)>=inactive_length)&&((data(i,1)-temp_x*gapd_length)<=(inactive_length+active_length)) && ((data(i,2)- temp_y*gapd_length)>=inactive_length)&&((data(i,2)-temp_y*gapd_length)<=(inactive_length+active_length))
         % CELL 判断是否为活动区域，进行条件语句
         gapd_active_count(temp_y+1, temp_x+1) = gapd_active_count(temp_y+1, temp_x+1) + 0.563;
         %如果为真，则添加 0.563 的 PD 值
     end
gapd_active_count(gapd_active_count>1) = 1;
%由于值不能超过1，所以超过1就设置为1
总和（总和（gapd_active_count））
% 添加所有落入活动区域的值

Answer 1

这样应该可以吧

data = rand(10000,2)*3
[cell_count,area_count]=get_map(data)

h =heatmap(cell_count);
title(h,'cell count')
h = heatmap(area_count);
title(h,'area count')

function [cell_count,area_count]=get_map(data)
w = 0.075;
area_b = [0.005, w-0.005];
cell_count = zeros(40,40);
area_count = zeros(40,40);
re = mod(data,w);
location = int32((data-re)/w);
for i=1:size(location,1)
    x = location(i,1)+1;
    y = location(i,2)+1;
    cell_count(x, y) = cell_count(x, y) + 1;
    if all(re(i)>area_b(1)) && all(re(i)<area_b(2))
        area_count(x,y) = area_count(x,y) + 1;
    end
end
end