Estudiante Ingenieria

f=imread("binuclear.tif");

figure, imshow(f)

[pic]

figure, hist(f(:),0:255)

[pic]

pmin=107; pmax=205;f_adj=255*(double(f)-pmin)/(pmax-pmin);

f_adj=uint8(f_adj);

figure, imshow(f_adj)

[pic]

figure, hist(f_adj(:),0:255)

[pic]

T=isodata(f_adj)

T = 140

g=(f_adj0)

w =0 0 0 0 1 0 0 0 0

0 0 1 1 1 1 1 0 0

0 1 1 1 1 1 1 1 0

0 1 1 1 1 1 1 1 01 1 1 1 1 1 1 1 1

0 1 1 1 1 1 1 1 0

0 1 1 1 1 1 1 1 0

0 0 1 1 1 1 1 0 00 0 0 0 1 0 0 0 0



sum(w(:))

ans = 49

bw_min=ordfilt2(g,1,w,"symmetric");

figure, imshow(bw_min)

[pic]bw_max=ordfilt2(bw_min,1,w,"symmetric");

figure, imshow(bw_max)

[pic]

L=bwlabel(bw_max);

figure, imshow(L,[]),colormap jet,colorbar

[pic]

stats=regionprops(L,"area");

area=[stats.Area]

area =



3 16 84 256 174 1092 244 764 886 596 41



[M,N]=size(f);

AT=M*N

AT =1228800

Area=100*area/AT

Area =

Columns 1 through 6:

2.4414e-004 1.3021e-003 6.8359e-003 2.0833e-002 1.4160e-002 8.8867e-002

Columns 7through 11:

1.9857e-002 6.2174e-002 7.2103e-002 4.8503e-002 3.3366e-003

Area'

ans =



2.4414e-004

1.3021e-003

6.8359e-0032.0833e-002

1.4160e-002

8.8867e-002

1.9857e-002

6.2174e-002

7.2103e-002

4.8503e-002

3.3366e-003



id=find(Area> 4.8501e-002 & Area