>> i=imread('04c10.tif');
>> level = graythresh(i);
>> BW = im2BW(i,level);
>> figure, imshow(BW);
>> [A,B]=find(I(i,j)<level);
??? Index exceeds matrix dimensions.
>> [i,j]=size(I);
>> for i=0:225
for j=0:225
[A,B]=find(I(i,j)<level);
[C,D]=find(I(i,j)>level);
I0=mean(men(I(C,D)));
I0=0
M=I0+log(I(A(i,j),B(i,j))/I0);
end
end
>>