灰度图像生成程序matlab

clear; close all; clc; rgb = imread('测试图像.jpg'); % Read the color image if size(rgb, 3) == 3 I = rgb2gray(rgb); else I = rgb; end A = sum(I(:)); B = sum(I(:)); I1 = [1:numel(A)]; AS = A .* I1; X = sum(AS) / sum(A); J1 = [1:numel(B)]; BS = B .* J1; Y = sum(BS) / sum(B); X1 = max(1, round(X - 1000)); X2 = min(size(I, 2), round(X + 1000)); Y1 = max(1, round(Y - 1000)); Y2 = min(size(I, 1), round(Y + 1000)); S1 = (X1:1:X2); C1 = (Y1:1:Y2); figure(1); imshow(I, []); title('Center of gravity') hold on; plot(X, Y, 'r+'); hold on; plot(X1, C1, 'r'); hold on; plot(X2, C1, 'r'); hold on; plot(S1, Y1, 'r'); hold on; plot(S1, Y2, 'r'); FIN = I(Y1:Y2, X1:X2); figure(2); imshow(FIN, []); title('The extraction results'); hold on; plot(1000, 1000, 'r+');