用2023a版matlab写一段将四副256*256的灰度图像I1、I2、I3、I4进行灰度加权平均图像融合后的图像进行几何质心提取固定大小的有效区域，其中有效区域为将标准图像内切圆中的内接正方形区域确定为标准图像的有效区域

好的，以下是一段用 Matlab 实现的代码，可以实现将四幅图像进行灰度加权平均后，提取有效区域的几何质心。

% 读取四个灰度图像
I1 = imread('image1.jpg');
I2 = imread('image2.jpg');
I3 = imread('image3.jpg');
I4 = imread('image4.jpg');

% 灰度加权平均图像融合
w1 = 0.25;
w2 = 0.25;
w3 = 0.25;
w4 = 0.25;
I_merge = w1*I1 + w2*I2 + w3*I3 + w4*I4;

% 确定有效区域
center_x = 128; % 圆心x坐标
center_y = 128; % 圆心y坐标
radius = 100; % 圆的半径
side = 2*radius/sqrt(2); % 正方形的边长
x1 = center_x - radius; % 正方形左上角x坐标
y1 = center_y - radius; % 正方形左上角y坐标

% 提取有效区域
I_effective = I_merge(y1:y1+side-1, x1:x1+side-1);

% 计算有效区域的几何质心
[x, y] = meshgrid(x1:x1+side-1, y1:y1+side-1);
x_weighted = sum(x(:).*double(I_effective(:)))/sum(double(I_effective(:)));
y_weighted = sum(y(:).*double(I_effective(:)))/sum(double(I_effective(:)));

% 在原图上显示有效区域和几何质心
figure;
imshow(I_merge);
hold on;
rectangle('Position', [x1, y1, side, side], 'EdgeColor', 'r');
plot(x_weighted, y_weighted, 'gx', 'LineWidth', 2, 'MarkerSize', 10);

需要注意的是，该代码中的四个图像需要与该 Matlab 脚本文件在同一目录下，并且需要将`image1.jpg`、`image2.jpg`、`image3.jpg`、`image4.jpg`替换为实际的图像文件名。此外，该代码中的有效区域为正方形，而不是圆形，因为只有正方形可以使用矩阵索引方便地提取有效区域。如果需要圆形有效区域，可以先提取正方形有效区域，然后使用`imellipse`函数创建一个内接圆形，再使用`createMask`方法将圆形区域提取出来。