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'); % 形态学膨胀 se = strel('disk', 2); dilated_edge = imdilate(canny_edge, se); % 连通域分析 stats = regionprops('table', dilated_edge, 'Area', 'Centroid'); % 筛选面积最大的连通区域 [~, idx] = max(stats.Area); centroid = stats.Centroid(idx, :); % 显示帧和质心 imshow(dilated_edge); hold on; plot(centroid(1), centroid(2), 'r+', 'MarkerSize', 10); hold off; end建立坐标系，实时显示质心坐标

好的，可以使用 MATLAB 的图像处理工具箱来实现这个功能。代码如下：

v =Reader('shoujilux7.mp4'); 
figure; % 新建一个图像窗口
axis([0 v.Width 0 v.Height]); % 设置坐标系范围
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', 2);
    dilated_edge = imdilate(canny_edge, se);

    % 连通域分析
    stats = regionprops('table', dilated_edge, 'Area', 'Centroid');

    % 筛选面积最大的连通区域
    [~, idx] = max(stats.Area);
    centroid = stats.Centroid(idx, :);

    % 显示帧和质心
    imshow(frame);
    hold on;
    plot(centroid(1), centroid(2), 'r+', 'MarkerSize', 10);
    hold off;
    drawnow; % 实时刷新图像窗口
end

需要注意的地方是，由于每一帧的图像大小可能不同，因此需要在代码中动态设置坐标系范围。另外，在循环中加入 `drawnow` 命令可以实现实时刷新图像窗口。