matlab-对视频目标进行光流提取,对比Lucas-Kanade和Horn-Schunck-源码

clc; clear; close all; warning off; addpath(genpath(pwd)); %% Lucas Kanade directorio = dir('Combo_4_(C4)/*png'); t_global = tic; for k = 1:length(directorio)-1 k im1 = directorio(k).name; img_1 = imresize(imread(strcat('Combo_4_(C4)/',im1)),0.5); img_1 = im2double(rgb2gray(img_1)); im2 = directorio(k+1).name; img_2 = imresize(imread(strcat('Combo_4_(C4)/',im2)),0.5); img_2 = im2double(rgb2gray(img_2)); % funcion de flujo óptico de lucas Kanade (tamaño de ventana, frame_1, % frame_2, distancia señalizaciones tic; [X, Y, u, v] = lucas_kanade(11,img_1,img_2,20); toc figure(1); imshow(img_2); hold on; quiver(X, Y, u, v, 'y'); hold off; savefig(strcat('results/',im1 ,'.fig')); pause(0.01); end toc(t_global) %% Lucas Kanade mediante ecuación directorio = dir('Combo_4_(C4)/*png'); t_global = tic; for k = 1:length(directorio)-1 k im1 = directorio(k).name; img_1 = imresize(imread(strcat('Combo_4_(C4)/',im1)),0.5); img_1 = im2double(rgb2gray(img_1)); im2 = directorio(k+1).name; img_2 = imresize(imread(strcat('Combo_4_(C4)/',im2)),0.5); img_2 = im2double(rgb2gray(img_2)); % funcion de flujo óptico de lucas Kanade (tamaño de ventana, frame_1, % frame_2, distancia señalizaciones tic; [X, Y, u, v] = lucas_kanade2(11,img_1,img_2,20); toc figure(1); imshow(img_2); hold on; quiver(X, Y, u, v, 'y'); hold off; savefig(strcat('results/',im1 ,'.fig')); pause(0.001); end toc(t_global) %% Metodo H&S directorio = dir('Combo_4_(C4)/*png'); alpha = 60; it = 50; %t_global = tic; for k = 1:length(directorio)-1 k im1 = directorio(k).name; img_1 = imresize(imread(strcat('Combo_4_(C4)/',im1)),0.5); img_1 = im2double(rgb2gray(img_1)); im2 = directorio(k+1).name; img_2 = imresize(imread(strcat('Combo_4_(C4)/',im2)),0.5); img_2 = im2double(rgb2gray(img_2)); % funcion de flujo óptico de lucas Kanade (tamaño de ventana, frame_1, % frame_2, distancia señalizaciones %tic; [X, Y, u, v] = Horn_Schunck(alpha,img_1,img_2,it,20); %toc figure(1); imshow(img_2); hold on; quiver(X, Y, u, v, 'y'); hold off; savefig(strcat('results/',im1 ,'.fig')); pause(0.001); end %toc(t_global)