clear;
clc;
I=imread('4556.jpg');
J=rgb2gray(I);
[m,n]=size(J);
BW=J;
gmax =0;
K=60;
tic;
for i=1:m
for j=1:n
if BW(i,j)>gmax
gmax=BW(i,j);
end
if BW(i,j)<K
BW(i,j)=0;
end
end
end
BW0=im2double(BW);
y0=zeros(1,n);
for j=1:n
x=0;
x1=0;
y1=0;
for i=1:m
x=x+BW0(i,j);
y1=y1+BW0(i,j)*i;
end
y0(j)=y1/x;
end
toc;
imshow(BW);
hold on
plot(1:n,y0,'r')