summation of a matrix inside circle

Question

JK 2021-11-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1583739-summation-of-a-matrix-inside-circle

编辑： DGM 2021-11-11

采纳的回答： DGM

Hello,

I have 2D NxN matrix array and I want to get summation of data inside circle from center.

The size of circle is varied from smallest i.e. one pixel at the center to largest i.e. N/2

I want to plot summation w.r.t. circle.

Thanks,

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

DGM 2021-11-10

Is N always odd? Is it always even? Could it be either?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

DGM 2021-11-10

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1583739-summation-of-a-matrix-inside-circle#answer_828739

编辑：DGM 2021-11-11

在 MATLAB Online 中打开

This should work regardless of whether N is odd or even.

s = 100;

A = ones(s); % test array

% build distance map

rmax = s/2 - 0.5;

x = linspace(-rmax,rmax,s);

D = sqrt(x.^2 + x.'.^2);

r = 0:ceil(D(1,floor(s/2)));

circint = zeros(1,numel(r));

for rind = 1:numel(r)

circint(rind) = sum(A(D<=r(rind)));

end

plot(r,circint)

2 个评论
显示无隐藏无

JK 2021-11-11

any thought on even case?

DGM 2021-11-11

编辑：DGM 2021-11-11

I modified it. Bear in mind that there is no center pixel in an array with even geometry. The center four pixels are 1/sqrt(2) away from the center. Consequently, the first element of the output is zero.

Similarly, the pixels straddling the midpoints of the edges of a 10x10 array are not at a radius of 5px. They're at a radius of 4.5227px. As the radius vector in this example is integer-valued, you can either choose to stop short at 4px or go a bit too far (5px). I chose to go to 5px. If you want to stop short in the even case, you can remove the ceil() from the definition of r.

请先登录，再进行评论。