MATLAB中不用循环生成圆盘（圆形）/圆环掩膜矩阵

           由于在计算中，经常需要生成圆盘（or叫圆形？）掩膜，所以特意写了一个生成圆盘的matlab函数，当然经过修改，也可以生成圆环函数。

       譬如，我们想在一个矩阵里提取以某一点为圆心，半径r内的元素，让这个圆范围外的元素都为0时，这就需要一个圆盘掩膜。但是按普通人的第一直觉都是利用对行和列的两重循环，把矩阵内的每一个元素都遍历一遍，看是否在圆内，如果在的话，赋值为1，不在的话，赋值为0 。对于这种传统做法，明显采用了两重循环，会大大降低matlab程序的效率，那有没有不用循环的方法呢？

       下面就讲如何不用循环生成圆盘掩膜：

        假设矩阵大小为512*512 ，取矩阵最中间一点为圆心，半径r=50。先用函数meshgrid生成矩阵对应的每一个坐标，然后用circle=x.^2+y.^2生成代表半径的平方的矩阵，最后用find函数找到符合要求的元素，赋值为1，不符合要求的赋值为0 。具体程序如下：

 

m=512; %矩阵的函数
n=512; %矩阵的列数
r=50;   %生成圆的半径
m1=-m/2:m/2-1;   %把圆心变到矩阵的中间
n1=-n/2:n/2-1;
[x,y]=meshgrid(m1,n1);
circle=x.^2+y.^2;   %计算出每一点到圆心的距离的平方

circ_mask=zeros(m,n);  
circ_mask(find(circle<=r*r))=1;  %找到圆内的元素，并复制为1
circ_mask(find(circle>r*r))=0;   %找到圆外的元素，并复制为0

figure;
imshow(255*circ_mask);

 

    程序运行结果如下：

     另外，如果要生成圆环的话，只需要把find函数那里做相应的改变即可。