Измерение пульса человека - портирование готового кода с Matlab

dataDir = './data';
vidFile = fullfile(dataDir,'face.mp4');
% compute Gaussian blur stack
disp('Spatial filtering...')
vid = VideoReader(vidFile);
% vid = videoinput('winvideo', 1);
alpha = 50;
samplingRate = 30;
chromAttenuation = 1;
fl = 50/60;
fh = 60/60;
level = 4;
startIndex = 1;
len = vid.NumberOfFrames;
endIndex = len-10;
 
% Extract video info
vidHeight = vid.Height;
vidWidth = vid.Width;
nChannels = 3;
temp = struct('cdata', zeros(vidHeight, vidWidth, nChannels, 'uint8'), 'colormap', []);
 
% firstFrame
temp.cdata = read(vid, startIndex);
[rgbframe, ~] = frame2im(temp);
rgbframe = im2double(rgbframe);
frame = rgb2ntsc(rgbframe);
 
blurred = blurDnClr(frame,level);
% create pyr stack
GDown_stack = zeros(endIndex - startIndex +1, size(blurred,1),size(blurred,2),size(blurred,3));
k = 1;
for i=startIndex+1:endIndex
    disp(k);
    k = k+1;
    temp.cdata = read(vid, i);
    [rgbframe,~] = frame2im(temp);
    rgbframe = im2double(rgbframe);
    frame = rgb2ntsc(rgbframe);
    blurred = blurDnClr(frame,level);
    GDown_stack(k,:,:,:) = blurred;
    filtered_stack = ideal_bandpassing(GDown_stack, 1, fl, fh, samplingRate);
    %% amplify
    filtered_stack(:,:,:,1) = filtered_stack(:,:,:,1) .* alpha;
    filtered_stack(:,:,:,2) = filtered_stack(:,:,:,2) .* alpha .* chromAttenuation;
    filtered_stack(:,:,:,3) = filtered_stack(:,:,:,3) .* alpha .* chromAttenuation;
    
    filtered = squeeze(filtered_stack(k,:,:,:));
    filtered = imresize(filtered,[vidHeight vidWidth]);
    filtered = filtered+frame;
    filtered = ntsc2rgb(filtered);
%     imshow(blurred, 'InitialMagnification',800);
    imshow(blurred);
    drawnow;
end

function filtered = ideal_bandpassing(input, dim, wl, wh, samplingRate)
    if (dim > size(size(input),2))
        error('Exceed maximum dimension');
    end
    input_shifted = shiftdim(input,dim-1);
    Dimensions = size(input_shifted);
    n = Dimensions(1);
    dn = size(Dimensions,2);
    Freq = 1:n;
    Freq = (Freq-1)/n*samplingRate;
    mask = Freq > wl & Freq < wh;
    Dimensions(1) = 1;
    mask = mask(:);
    mask = repmat(mask, Dimensions);
    F = fft(input_shifted,[],1);
    F(~mask) = 0;
    filtered = real(ifft(F,[],1));
    filtered = shiftdim(filtered,dn-(dim-1));
end