celllbp.m

基于的LBP算法的细胞分类算法

close all;
clc;
clear;
Btime=clock
cellfilename = [];
FileName='cell000.bmp';
%将图像读入一个文件中。
%**********************
number=60;
	for j=1:number
           if j < 10
              FileName(7) = int2str(j);
				  cellfilename  = [cellfilename;FileName];
           elseif j < 100
              FileName_temp = int2str(j);
              FileName(6) = FileName_temp(1);
              FileName(7) = FileName_temp(2);
				  cellfilename  = [cellfilename;FileName];
              elseif j<1000
              FileName_temp = int2str(j);
              FileName(5) = FileName_temp(1);
              FileName(6) = FileName_temp(2);
              FileName(7) = FileName_temp(3);
				  cellfilename  = [cellfilename;FileName];
              else
                  FileName_temp = int2str(j);
              FileName(4) = FileName_temp(1);
              FileName(5) = FileName_temp(2);
              FileName(6) = FileName_temp(3);
              FileName(7) = FileName_temp(4);
     		      cellfilename  = [cellfilename;FileName];              
           end
     end
CellHst = [];

for i=1:number
    Data = imread(cellfilename (i,:));
    GrayImageData = rgb2gray(Data);
    Hst = LBP_riu81(GrayImageData);
    ImageData = reshape(Hst,10,1);
    CellHst = [CellHst,ImageData];
    %clear Data;
end

CellTrainHst = [CellHst(:,1:20),CellHst(:,31:50)];
CellTestHst = [CellHst(:,21:30),CellHst(:,51:60)];
SumTemp = [];
for Loop1 = 1:20
    for Loop2 = 1:40
        Temp1 = CellTestHst(:,Loop1)-CellTrainHst(:,Loop2);
        Temp2 = CellTestHst(:,Loop1)+CellTrainHst(:,Loop2);
        Temp3 = Temp1.* Temp1;
        Temp4 = Temp3./ Temp2;
        SumTemp = [ SumTemp,sum(Temp4)];
    end 
end

for Loop3 = 1:20
    a = min(SumTemp(1,40*(Loop3 - 1) + 1 : 40*Loop3));
    for i = 1:40
        if SumTemp(1 ,40*(Loop3 - 1) + i) > a
            SumTemp(1 ,40*(Loop3 - 1) + i) = 0;
        else SumTemp(1 ,40*(Loop3 - 1) + i) = 1;
        end
    end
end

SumTempNew = reshape(SumTemp,40,20);

Result1 = SumTempNew(1:20,1:10);
Result2 = SumTempNew(21:40,11:20);
Result = sum(sum(Result1)) + sum(sum(Result2));

Class = zeros(1,20);
for j=1:20
    for i=1:40
        if SumTempNew(i,j)==1
            Class(1,j)=i;
        else
        end
    end
end

fprintf('********** The Performance of  Testing is %f **********\n',Result/20);
fprintf('********** The Performance of  Testing Class One is %f **********\n',sum(sum(Result1))/10);
fprintf('********** The Performance of  Testing Class Two is %f **********\n',sum(sum(Result2))/10);