用2023a版matlab写一段将四副256*256的灰度图像I1、I2、I3、I4进行灰度加权平均图像融合后的图像进行几何质心提取固定大小的有效区域

好的，我理解了你的问题。以下是代码实现：

% 读取四幅256*256的灰度图像
I1 = imread('image1.png');
I2 = imread('image2.png');
I3 = imread('image3.png');
I4 = imread('image4.png');

% 灰度加权平均融合
alpha = 0.25; % 权重
I = alpha*I1 + alpha*I2 + alpha*I3 + alpha*I4;

% 提取有效区域
threshold = graythresh(I); % 自适应阈值
BW = imbinarize(I, threshold); % 二值化
stats = regionprops(BW, 'BoundingBox'); % 获取边界框信息
box = stats.BoundingBox; % 边界框
x = round(box(1));
y = round(box(2));
w = round(box(3));
h = round(box(4));

% 获取有效区域
ROI = I(y:y+h-1, x:x+w-1);

% 计算几何质心
props = regionprops(BW, 'Centroid');
centroid = props.Centroid; % 几何质心

% 显示结果
subplot(1,2,1);
imshow(I); title('灰度加权平均融合图像');
subplot(1,2,2);
imshow(ROI); title('有效区域');
hold on;
plot(centroid(1), centroid(2), 'r+', 'MarkerSize', 10); % 显示几何质心

请将代码保存为.m文件，并将四幅灰度图像命名为image1.png、image2.png、image3.png、image4.png放在同一目录下，然后运行即可得到灰度加权平均融合图像和有效区域，并在有效区域中显示几何质心。