xingxin.m

求取一幅灰度图像中的光斑形心的位置的Matlab的源代码。算法使用的是最大类间方差法编写的程序

tic
clc
clear   

I=imread('1.bmp');
subplot(2,2,1);
imshow(I);
title('原始图像'); 
I=rgb2gray(I);
subplot(2,2,2);
imshow(I);
title('灰度图像'); 

[Row,Col]=size(I);
Row
Col
I1=medfilt2(I);
I=I1;
subplot(2,2,3);
imshow(I); 
title('中值滤波');
[ni,x] = imhist(I,255);

pi=ni/(Row*Col);

ave_gray=0;
for i=1:255
    ave_gray=ave_gray+i*pi(i);
end

t=0;
for k=1:254
    w0=0;
    w1=0;
    u0=0;
    u1=0;
   
	for i=1:k
        w0=w0+pi(i);
        if w0==0
            w0=0.001;
        end
        u0=u0+i*pi(i);        
	end
    u0=u0/w0; 
    
	for i=k+1:255
        w1=w1+pi(i);
        if w1==0
            w1=0.001;
        end
        u1=u1+i*pi(i);
	end
    u1=u1/w1;
   
    d_jian(k)=w0*(u0-ave_gray)^2+w1*(u1-ave_gray)^2;
    if d_jian(k)>t
        t=d_jian(k);
        m=k;
    end
end
m

hgray=0;
x(Row)=0;

for i=30:Row-30
    for j=30:Col-30
        if I(i,j) >m
            %I(i,j) = 0;
            hgray=hgray+i*double(I(i,j));
            x(i)=x(i)+double(I(i,j));
            
            %else
            %I(i,j) =255;
        end
        xc=fix(hgray/x(i));
        xc
    end
end

subplot(2,2,4);
imshow(I);
title(' 最大类间方差法分割 ')