test.m

这是我们组参加2008年全国大学生数学建模时A题(数码相机定位)所用到得程序

clear
M=getM;
I=imread('shumo.bmp');
level = graythresh(I);
BW = im2bw(I,level);%转化为二值图象：0黑1白
[rowNum,columnNum]=size(BW);
Black=[];
for i=1:1024
    for j=1:768
        if(BW(j,i)==0)
            Black=[Black;i,j];%找出所有的像点，即所有的黑点
        end
    end
end
r=size(Black,1);              
disp('所有特征点在像上的像素点个数为：');
disp(r);
num=0;
for i=1:r
    aim=getObjectPoint(M,Black(i,1),Black(i,2));
    x=aim(1);
    y=aim(2);
    if(((x-12)^2+(y-12)^2<=144)|((x-42)^2+(y-12)^2<=144)|((x-112)^2+(y-12)^2<=144)|((x-12)^2+(y-112)^2<=144)|((x-112)^2+(y-112)^2<=144))
        num=num+1;%统计像素点落在A~E五个圆内的数量
    end
end
disp('经统计后像素点落在A~E五个圆内的数量为：');
disp(num);
disp('误码率为：');
disp((r-num)/r);