寻找彩色图片目标物体的质心

clc,close all; %converts an rgb to hsi rgb=imread('1.bmp'); rgb = im2double(rgb); r = rgb(:, :, 1); g = rgb(:, :, 2); b = rgb(:, :, 3); %Implement the conversion equations. num = 0.5*((r - g) + (r - b)); den = sqrt((r - g).^2 + (r - b).*(g - b)); theta = acos(num./(den + eps)); H = theta; H(b > g) = 2*pi - H(b > g); H = H/(2*pi); num = min(min(r, g), b); den = r + g + b; den(den == 0) = eps; S = 1 - 3.* num./den; H(S == 0) = 0; I = (r + g + b)/3; % Combine all three results into an hsi image. hsi = cat(3, H, S, I); figure,imshow(hsi); title('HSI图像'); %figure,imshow(I); %title('HSI图像中的强度'); %figure,imshow(H); %title('HSI图像中色调'); figure,imshow(S); title('HSI图像中的饱和度'); %filter K=filter2(fspecial('average',3),S); figure,imshow(K); title('滤波后的图'); figure,imhist(K); title('滤波后的直方图') %T=0.3 K2=im2bw(K,0.3); figure,imshow(K2); title('阈值后的图像'); %marker the center L=bwlabel(K2); hold on %marker centre in K2 % for k= 1:n % [r,c]=find(L==k); % rbar = mean(r); % cbar = mean(c); % plot(cbar,rbar,'Marker','o','MarkerEdgeColor','k','MarkerFaceColor','k','MarkerSize',10); % plot(cbar,rbar,'Marker','*','MarkerEdgeColor','w'); % end %search the centre , Coordinate data in centroids s = regionprops(L, 'centroid'); centroids = cat(3, s.Centroid); centroids; % for Original picture imshow(rgb); title('标记后的结果') for k=1:8 % marker Center in 'o' text( centroids(1,1,k),centroids(1,2,k),'o','FontSize',8) % marker number eval(['text( centroids(1,1,k),centroids(1,2,k),','''', num2str(k) ,'''', ',' ,'''', 'FontSize' ,'''', ',16)']) end