在2D阵列中集成圆形

Question

我试图整合叠加在像素阵列上的圆形孔径区域（见下图）。但是，我需要确定圆孔内的通量（面积）的分数，并在圆的边界上的每个正方形像素中省略圆孔外的任何东西。

我如何在numpy / python中对此进行编码，以便准确测量圆圈内的通量？

Answer 1

使用微积分计算圆内每个像素的比例。 （将圆的等式积分在每个的左右边界之间。）

Answer 2

1. 在您选择的图像编辑器中的半径黑色背景上绘制一个白色圆圈，然后保存输出的位图
2. 使用del在您的代码中加载图片，并将像素值除以255，以便在scipy.misc.imload
中设置面具
3. 计算该掩码与要整合的数据的乘积之和