2 .2.2. Алгоритм вычисления характеристик гистограммы расстояний и ее построение.

I. Расчет расстояний между клетками.

n=1;

u=1:(length(centroids1(:,1))-1);

ras=zeros(sum(u),1);

for i=1:(length(centroids1(:,1))-1)

for j=i+1:length(centroids1(:,1))

ras(n)=sqrt((centroids1(i,1)-centroids1(j,1))^2+((centroids1(i,2)-centroids1(j,2))^2));

n=n+1;

end

end

II. Построение гистограммы расстояний

[M,N]=size(I(:,:,1));

k=1:sqrt((M-1)^2+(N-1)^2);

h=hist(ras,k);

figure;plot(h);

title('Гистограмма распределения расстояний между парами олигодендроцитов');

xlabel('Расстояние в пикселях');

ylabel('Частота расположения олигодендроцитов на данном расстоянии');

n=1;

ras2=zeros(length(centroids1),1);

ho=[];

k=1:sqrt((M-1)^2+(N-1)^2)/10;

for i=1:length(centroids2)

for j=1:length(centroids1)

ras2(n)=sqrt((centroids2(i,1)-centroids1(j,1))^2+((centroids2(i,2)-centroids1(j,2))^2));

n=n+1;

end

h2=hist(ras2,k*10);

ho=[ho;h2]; %#ok<AGROW>

n=1;

end

figure;imshow(ho,[]);

hold on

Y(1:length(ho(1,:)))=round(1+length(ho(:,1))/4);

X=1:length(ho(1,:));

plot(X,Y,'r-');

Y(1:length(ho(1,:)))=round(length(ho(:,1))-length(ho(:,1))/4);

plot(X,Y,'r-');

hold off

III. Усреднение полученной гитограммы.

i=round(1+length(ho(:,1))/4):round(length(ho(:,1))-length(ho(:,1))/4);

sum(ho(i,:));

hosr=sum(ho(i,:))/length(i);

hx=1:10:sqrt((M-1)^2+(N-1)^2)-10;

figure;plot(hx,hosr);

title('Усредненная гистограмма расположения олигодендроцитов отностельно нейронов');

xlabel('Расстояние в пикселях');ylabel('Средняя частота расположения');

3 Экспериментальные исследования

17