博主github:https://github.com/MichaelBeechan
博主CSDN:https://blog.csdn.net/u011344545
function [u, v] = optic_flow_brox(img1, img2)
img1 = imread('frame07.png');
img2 = imread('frame08.png');
alpha = 30.0 ; % Global smoothness variable.
gamma = 80.0 ; % Global weight for derivatives.
[ht, wt, dt] = size( img1 ) ;
num_levels = 40 ; % for face; Number of pyramid levels. Can also be calculated automatically
%power(0.95,40) = 0.1285
im1_hr = gaussianRescaling( img1, power( 0.95, num_levels ) ) ; % Rescaling images
im2_hr = ga