【原创】matlab每日练习：用matlab…

【原创】matlab每日练习：用matlab作圆（转载请注明出处：http://blog.sina.com.cn/hapcat）

  作者：曾维 

    在图像识别中，圆是最基本的图形之一。怎么样在一个二维矩阵中标识出一个圆，及相关算法是本文将要介绍的。这将为进一步学习它的逆过程，即圆的识别，圆心坐标的确立，作一个铺垫。

1. 产生一个单像素圆环

思路：圆的参数方程 其中(a,b)为圆心坐标，r为圆的半径，设单位为像素，则r 即为半径的像素个数。

x-a=rcos(w)

y-b=rsin(w)

估计一个圆上的点的个数，圆的周长为2*pi*r,即每个扫描的小角度增量为(2*pi)/2*pi*r=1/r（单位为弧度）。

角度如果画一个整圆,则w=0 到w=2*pi

代码如下

 

r=100;%半径

O=[100 100];%圆心

Bkground=zeros(2*r,2*r);%画布

dat=1/r;

for i=0:1:2*pi*r

%画图,加一是因为matlab的%数组就从1开始的.如a(1),没有a(0)

Bkground(int32(r*cos(i*dat))+O(1)+1,int32(r*sin(i*dat))+O(2)+1)=1;

end

imshow(Bkground);

 

结果如下图

                    图1-1

2. 建立一个数组,保存圆周的坐标数据.

思路:画一个圆,圆周上的像素个数为int32(2*pi*r). 则保存这写点的坐标需要的数组为int32(2*pi*r) 2个单元.

代码如下

r=100;%半径

C_data=zeros(int32(2*pi*r),2);% 创建用于保存圆的坐标的数组

O=[1 1];%注意圆心

dat=1/r;

for i=1:1:2*pi*r

%画图,加一是因为matlab的%数组就从1开始的.如a(1),没有a(0)

C_data(i,1)=int32(r*cos(i*dat))+O(1);

C_data(i,2)=int32(r*sin(i*dat))+O(2);

end

%根据C_data作一个圆试试

Bkground=zeros(300,300);%黑色画布

for i=1:1:2*pi*r

%圆心坐标为(120,120).在matlab中数组标为Bkground (121,121)

Bkground(120+C_data(i,1),120+C_data(i,2))=1; 

end

imshow(Bkground);

运行结果如下图,

图1-2