Bezdepbonusu7Программирование: статьи и ответы
Меню

Вопрос

Изображение истории движения в Matlab

Автор вопроса: Hind FawziОтветов: 1

Я работаю над проектом о распознавании действий с использованием изображений истории движения в Matlab. Я новичок в этой области. Я сделал вычитание фона, используя метод разности кадров, чтобы получить изображения, на которых есть только движущийся человек. Теперь я хочу вычислить MHI. Я нашел следующий код для MHI. Я не понял, что такое fg{1} и как его использовать. Любая помощь будет оценена. Спасибо.

vid= VideoReader('PS7A1P1T1.avi');
n = vid.NumberOfFrames;

fg = cell(1, n);

 for i = 1:n
      frame = read(vid,i);
      frame = rgb2gray(frame);
      fg{i} = frame;
 end
%---------------------------------------------------------------
%background subtraction using frame differencing method
thresh = 25;
bg = fg{1};           % read in 1st frame as background frame
% ----------------------- set frame size variables -----------------------
fr_size = size(bg);
width = fr_size(2);
height = fr_size(1);
% --------------------- process frames -----------------------------------
for i = 2:n

    fr = fg{i};       % read in frame

    fr_diff = abs(double(fr) - double(bg));  % cast operands as double to avoid negative overflow

    for j=1:width                 % if fr_diff > thresh pixel in foreground
        for k=1:height
            if ((fr_diff(k,j) > thresh))
                fg {i}(k,j) = fr(k,j);
            else
                fg {i}(k,j) = 0;
            end
        end
    end

    bg = fr;

    imshow(fg{i})

end

out = MHI(fg);

//------------------------------------------------------- функция MHI = MHI(fg)

% Initialize the output, MHI a.k.a. H(x,y,t,T)
MHI = fg;

% Define MHI parameter T
T = 15; % # of frames being considered; maximal value of MHI.

% Load the first frame
frame1 = fg{1};

% Get dimensions of the frames
[y_max x_max] = size(frame1);

% Compute H(x,y,1,T) (the first MHI)
MHI{1} = fg{1} .* T;

% Start global loop for each frame
for frameIndex = 2:length(fg)

    %Load current frame from image cell
    frame = fg{frameIndex};

    % Begin looping through each point
    for y = 1:y_max
        for x = 1:x_max
            if (frame(y,x) == 255)
                MHI{frameIndex}(y,x) = T;
            else
                if (MHI{frameIndex-1}(y,x) > 1)
                    MHI{frameIndex}(y,x) = MHI{frameIndex-1}(y,x) - 1;
                else
                    MHI{frameIndex}(y,x) = 0;
                end
            end
        end
    end
end
  • Если я помог, подумайте о том, чтобы принять мой ответ, чтобы сообщить сообществу, что вам больше не нужна помощь. Удачи rayryeng 28.12.2015

computer-visionmatlabvideo-processingimage-processingmotion

Ответы

Ответ 1 · рейтинг 1 · rayryeng

fg{1}, скорее всего, является первым кадром видео в градациях серого. Учитывая ваши комментарии, вы используете класс VideoReader для чтения во фреймах. Таким образом, считывайте каждый кадр отдельно, преобразуйте в оттенки серого, а затем размещайте в ячейке массива ячеек. Когда вы закончите, вызовите скрипт.

Вот код, измененный из ваших комментариев, чтобы удовлетворить эту задачу:

vid = VideoReader('PS7A1P2T1.avi');
n = vid.NumberOfFrames;

fg = cell(1, n);

for i = 1:n
   frame = read(vid,i);
   frame = rgb2gray(frame);
   fg{i} = frame;
end

Затем вы можете вызвать сценарий MHI:

out = MHI(fg);

Перевод вопроса и ответов сообщества Stack Overflow, лицензия CC BY-SA. Оригинал на Stack Overflow.

Похожие материалы