evm 代码分析（一）

处理时：

• 为了放大运动，我们推荐我们新的基于相位的管道（SIGGRAPH 2013）。
• 要放大颜色，请使用线性管道（SIGGRAPH 2012）。
• 选择您想要放大的正确时间尺度。 例如，成年人的心跳往往每秒发生一次，对应于 1Hz，您可以将内容放大 0.5Hz 到 3Hz 以确保安全。 间隔越窄，放大越集中，噪声被放大的越少，但有遗漏物理现象的风险。
• 在指定时间通带时不要忘记考虑视频帧速率！ 有关示例，请参见我们的代码。

先介绍2012相位欧拉代码

reproduceResultsSiggraph12

% Generates all the results for the SIGGRAPH paper at:
% http://people.csail.mit.edu/mrub/vidmag
%
% Copyright (c) 2011-2012 Massachusetts Institute of Technology, 
% Quanta Research Cambridge, Inc.
%
% Authors: Hao-yu Wu, Michael Rubinstein, Eugene Shih
% License: Please refer to the LICENCE file
% Date: June 2012
%

clear;
dataDir = './data';
resultsDir = 'ResultsSIGGRAPH2012';

mkdir(resultsDir);

%% baby
% inFile = fullfile(dataDir,'baby.mp4'); % 路径拼接
% fprintf('Processing %s\n', inFile);
% amplify_spatial_lpyr_temporal_iir(inFile, resultsDir, 10, 16, 0.4, 0.05, 0.1);
% amplify_spatial_lpyr_temporal_iir 函数是一种放大算法，即空间滤波采用 拉普拉斯金字塔，时域滤波采用 IIR滤波器。

% Alternative processing using butterworth filter  使用巴特沃斯滤波器的替代处理
% amplify_spatial_lpyr_temporal_butter(inFile, resultsDir, 30, 16, 0.4, 3, 30, 0.1);
% 空间滤波采用 lpyr 拉普拉斯金字塔，时域滤波采用 butter 巴特沃斯滤波器

%% baby2
% inFile = fullfile(dataDir,'baby2.mp4');
% fprintf('Processing %s\n', inFile);
% amplify_spatial_Gdown_temporal_ideal(inFile,resultsDir,150,6, 140/60,160/60,30, 1);
% 空间滤波采用 Gaussian金字塔  时域滤波采用 Ideal bandpass

%% camera
% inFile = fullfile(dataDir,'camera.mp4');
% fprintf('Processing %s\n', inFile);
% amplify_spatial_lpyr_temporal_butter(inFile, resultsDir, 150, 20, 45, 100, 300, 0);
% 空间滤波采用 lpyr 拉普拉斯金字塔，时域滤波采用 butter 巴特沃斯滤波器


%% subway
% inFile = fullfile(dataDir,'subway.mp4');
% fprintf('Processing %s\n', inFile);
% amplify_spatial_lpyr_temporal_butter(inFile, resultsDir, 60, 90, 3.6, 6.2, 30, 0.3);
% 空间滤波采用 lpyr 拉普拉斯金字塔，时域滤波采用 butter 巴特沃斯滤波器

%% wrist
%% No mask is used here to generate the output video.
% inFile = fullfile(dataDir,'wrist.mp4');
% fprintf('Processing %s\n', inFile);
% amplify_spatial_lpyr_temporal_iir(inFile, resultsDir, 10, 16, 0.4, 0.05, 0.1);
% 空间滤波采用 lpyr 拉普拉斯金字塔，时域滤波采用 IIR滤波器

% Alternative processing using butterworth filter
% amplify_spatial_lpyr_temporal_butter(inFile, resultsDir, 30, 16, 0.4, 3, 30, 0.1);
% 空间滤波采用 lpyr 拉普拉斯金字塔，时域滤波采用 butter 巴特沃斯滤波器


%% shadow
% inFile = fullfile(dataDir,'shadow.mp4');
% fprintf('Processing %s\n', inFile);
% amplify_spatial_lpyr_temporal_butter(inFile, resultsDir, 5, 48, 0.5, 10, 30, 0);
% 空间滤波采用 lpyr 拉普拉斯金字塔，时域滤波采用 butter 巴特沃斯滤波器

%% guitar
% inFile = fullfile