201223-记录numpy和Matlab下制作圆形掩膜

一. 引言

图像处理中change需要使用numpy或Matlab制作圆形掩膜，记录如下。
参考网址：Mask a circular sector in a numpy array

1.numpy版代码如下：

mask=np.zeros((1,256,256,1))
H, W=256,256
x, y = np.ogrid[:H, :W]
cx, cy = H/2, W/2
radius = int(np.random.uniform(0.9, 0.98) * H/2)
r2 = (x - cx) * (x - cx) + (y - cy) * (y - cy)
circmask = r2 >= radius * radius
circmask[:int(H/4), :] = 0
circmask[int(H*3/4):, :] = 0
mask[:, circmask, :] = 1

2.matlab版代码如下：

mask= zeros(512,512);
x=1:512;
y=x;
[X,Y]=meshgrid(x,y);
A=9*(X-12).^2+9*(X-12).*(Y-10)+3*(Y-10).^2;
mask(A<49)=1;