% 读取视频帧序列 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形态学处理填充，有边缘求质心坐标

可以在 Canny 边缘检测后，对边缘进行形态学处理，例如使用 imfill 函数进行填充，然后使用 regionprops 函数求取连通区域的质心坐标。

代码示例：

% 读取视频帧序列
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');

    % 进行形态学处理填充
    filled_edge = imfill(canny_edge, 'holes');

    % 求取连通区域的质心坐标
    stats = regionprops(filled_edge, 'Centroid');
    centroids = cat(1, stats.Centroid);
    disp(centroids);

    % 显示处理后的帧
    imshow(filled_edge);
end