matlab识别中心点

clear all;
clc
dir_name = 'E:\HUST\group_medicine_big_data\source_images_and_corresponding_ROI';
dir_root = dir('E:\HUST\group_medicine_big_data\source_images_and_corresponding_ROI');
x = zeros(length(dir_root),1);
y = zeros(length(dir_root),1);
z(1)={'aa'};
z(2)={'bb'};
for ii = 3:length(dir_root)
%     pic_num = '37';
%     suffix = '.png';
%     pic_name = [pic_num suffix];
    pic_name = dir_root(ii).name;
    z(ii) = {pic_name};
    full_pic_name = [dir_name '\' pic_name];
    pic = imread(full_pic_name);
    [m,n,~] = size(pic);
    R = pic(:,:,1);
    G = pic(:,:,2);
    B = pic(:,:,3);
    for i =1:m
        for j = 1:n
            if R(i,j) == 255 && G(i,j) == 128 && B(i,j) == 0
                x1 = i;
                y1 = j;
                break
            end
        end
    end

    for i = m:-1:1
        for j = n:-1:1
            if R(i,j) == 255 && G(i,j) == 128 && B(i,j) == 0
                x2 = i;
                y2 = j;
                break
            end
        end
    end

    x(ii) = (x1+x2)/2;
    y(ii) = (y1+y2)/2;
    


end
z = z';
xy = cat(2,x,y); 
xlswrite('E:\HUST\group_medicine_big_data\1.csv',z,'data','A');
xlswrite('E:\HUST\group_medicine_big_data\1.csv',x,'data','B');
xlswrite('E:\HUST\group_medicine_big_data\1.csv',y,'data','C');