% 读取视频帧序列 v = VideoReader('shoujilux7.mp4'); % 遍历视频帧序列 while hasFrame(v) % 读取当前帧并将其转换为灰度图像 frame = readFrame(v); gray_frame = rgb2gray(frame); % 进行gamma校正 gamma = 1.5; gamma_corrected = imadjust(gray_frame,[],[],gamma); % 进行高斯滤波 sigma = 1; hsize = ceil(6sigma); h = fspecial('gaussian', hsize, sigma); filtered_frame = imfilter(gamma_corrected, h); % 对滤波后的图像进行模拟退火算法改进下的Otsu阈值分割 T = graythresh(filtered_frame); [m, n] = size(filtered_frame); E = bwareaopen(im2bw(filtered_frame, T), round(mn/1000), 8); % 进行Canny边缘检测 canny_edge = edge(E, 'canny'); % 显示处理后的帧 imshow(canny_edge); end形态学处理填充，建立像素坐标系，由边缘求质心坐标，显示在图上

可以在代码的最后加上以下部分来进行形态学处理填充、建立像素坐标系、求质心坐标并显示在图上：

% 进行形态学处理填充
se = strel('disk', 5);
filled_frame = imclose(canny_edge, se);

% 建立像素坐标系
[x, y] = meshgrid(1:n, 1:m);

% 求质心坐标
centroid_x = sum(x(filled_frame)) / sum(filled_frame(:));
centroid_y = sum(y(filled_frame)) / sum(filled_frame(:));

% 在图上显示质心
hold on
plot(centroid_x, centroid_y, 'r*', 'MarkerSize', 15);
hold off

最后整个代码应该是这样的：

v = VideoReader('shoujilux7.mp4');

while hasFrame(v)
    frame = readFrame(v);
    gray_frame = rgb2gray(frame);
    
    % 进行gamma校正
    gamma = 1.5;
    gamma_corrected = imadjust(gray_frame,[],[],gamma);

    % 进行高斯滤波
    sigma = 1;
    hsize = ceil(6*sigma);
    h = fspecial('gaussian', hsize, sigma);
    filtered_frame = imfilter(gamma_corrected, h);

    % 对滤波后的图像进行模拟退火算法改进下的Otsu阈值分割
    T = graythresh(filtered_frame);
    [m, n] = size(filtered_frame);
    E = bwareaopen(im2bw(filtered_frame, T), round(m*n/1000), 8);

    % 进行Canny边缘检测
    canny_edge = edge(E, 'canny');

    % 进行形态学处理填充
    se = strel('disk', 5);
    filled_frame = imclose(canny_edge, se);

    % 建立像素坐标系
    [x, y] = meshgrid(1:n, 1:m);

    % 求质心坐标
    centroid_x = sum(x(filled_frame)) / sum(filled_frame(:));
    centroid_y = sum(y(filled_frame)) / sum(filled_frame(:));

    % 在图上显示质心
    hold on
    plot(centroid_x, centroid_y, 'r*', 'MarkerSize', 15);
    hold off
    
    % 显示处理后的帧
    imshow(canny_edge);
end

注意，由于在每个循环中都要绘制质心，因此需要在`imshow`之前加上`hold on`，在`plot`之后加上`hold off`。