How to divide a binary image circularly into "n" non overlapping equal size sectors around its center of gravity to count foreground pixels in each sector?

How to divide a binary image circularly into "n" non overlapping equal size sectors around its center of gravity to count foreground pixels in each sector?

Answer 1

[Matt J]

Use NDGRID and FIND,

[M,N]=size(binaryImage);

[X,Y]=ndgrid((1:M)-centerX,(1:N)-centerY)
[i,j]=find(binaryImage);

Then use cart2pol to convert to polar coordinates,

   theta = cart2pol(X(i),Y(j));

and finally apply histcounts() to theta.