matlab中心型程序,MatLab计算图形中心(2)

clear;

f0=imread('1.jpg');

f1=double(f0);

f2=(f1(:,:,1)+f1(:,:,2)+f1(:,:,3))/3; %转化为灰度图像.

subplot(2,2,1);imshow(f0);title('彩色图像');

subplot(2,2,2);imshow(uint8(f2));title('灰度图像');

[m,n]=size(f2);

for i=1:m %二值化处理

for j=1:n

if f2(i,j)<180

f(i,j)=1;

else

f(i,j)=0;

end

end

end

subplot(2,2,3);imshow(f,[0,1]);title('二值图像(黑白图像)');

x=0;

s=0;

for i=1:m

for j=1:n

x=x+j*f(i,j);

s=s+f(i,j);

end

end

x=x/s;

y=0;

for i=1:m

for j=1:n

y=y+i*f(i,j);

end

end

y=y/s;

hold on

subplot(2,2,3);plot(x,y,'+r');%质心显示

hold on

subplot(2,2,3);plot(n/2,m/2,'+b');