鱼眼镜头校正matlab程序

包含测试图像、测试结果 采用matlab实现鱼眼镜头标定 方法简答，可靠，标定效果可以

function[RAB]=circ（xyN）
x1=0;
x2=0;
x3=0;
y1=0;
y2=0;
y3=0;
x1y1=0;
x1y2=0;
x2y1=0;
for i =1:N
x1=x1+x（i）;
x2=x2+x（i）*x（i）;
x3=x3+x（i）*x（i）*x（i）;
y1=y1+y（i）;
y2=y2+y（i）*y（i）;
y3=y3+y（i）*y（i）*y（i）;
x1y1=x1y1+x（i）*y（i）;
x1y2=x1y2+x（i）*y（i）*y（i）;
x2y1=x2y1+x（i）*x（i）*y（i）;
end
C=N*x2-x1*x1;
D=N*x1y1-x1*y1;
E=N*x3+N*x1y2-（x2+y2）*x1;
G=N*y2-y1*y1;
H=N*x2y1+N*y3-（x2+y2）*y1;
a=（H*D-E*G）/（C*G-D*D）;
b=（H*C-E*D）/（D*D-G*C）;
c=-（a*x1+b*y1+x2+y2）/N;

A=a/（-2）; %x坐标
B=b/（-2）; %y坐标
R=sqrt（a*a+b*b-4*c）/2;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     276879  2012-08-01 10:26  校正\060-1.fig

     文件     415798  2012-07-18 09:51  校正\060.bmp

     文件     268030  2012-08-01 09:39  校正\060.fig

     文件     415798  2012-07-18 09:51  校正\080.bmp

     文件     415798  2012-07-18 09:51  校正\090l.bmp

     文件     415798  2012-07-18 09:51  校正\090r.bmp

     文件     283827  2012-08-01 09:34  校正\090r.fig

     文件     415798  2012-07-20 16:14  校正\90.bmp

     文件        533  2012-07-30 15:32  校正\circ.m

     文件    4001078  2012-08-01 19:52  校正\correct_image_061g.bmp

     文件       3203  2012-07-28 22:38  校正\fitellipse.m

     文件        621  2012-08-01 09:57  校正\markpoint.m

     文件       3368  2012-08-04 10:40  校正\rec.asv

     文件       3399  2012-08-04 11:16  校正\rec.m

     文件       2067  2012-08-01 15:54  校正\rec2.asv

     文件       2067  2012-08-01 15:55  校正\rec2.m

     文件       3777  2012-07-27 15:30  校正\shiyan5.m

     文件       1626  2012-08-01 11:04  校正\test.asv

     文件       1668  2012-08-01 16:43  校正\test.m

     文件       4018  2012-08-01 10:09  校正\test2.asv

    ..A.SH.     36864  2012-07-18 09:51  校正\Thumbs.db

     文件       1538  2012-08-02 16:59  校正\trapezial_correction-wrong.m

     文件       3308  2012-08-02 17:23  校正\trapezial_correction.asv

     文件       1695  2012-08-04 10:09  校正\trapezial_correction.m

     文件       1480  2012-07-26 20:43  校正\Untitled.m

     文件       4075  2012-07-25 15:44  校正\Untitled3.m

     文件       3845  2012-07-23 08:58  校正\Untitled5.m

     文件       3914  2012-07-23 15:35  校正\Untitled6.m

     文件     128822  2012-07-29 20:25  校正\截图.jpg

     目录          0  2013-04-09 09:07  校正

............此处省略3个文件信息