% Copyright (c) 2019 Lara Raad <lara.raad@upf.edu>,
% Enric Meinhardt-Llopis <enric.meinhardt@cmla.ens-cachan.fr>
%
% This program is free software: you can redistribute it and/or modify
% it under the terms of the GNU Affero General Public License as
% published by the Free Software Foundation, either version 3 of the
% License, or (at your option) any later version.
%
% You should have received a copy of the GNU Affero General Public License
% along with this program. If not, see <http://www.gnu.org/licenses/>.

% laplace beltrami flow interpolation
% function []=laplace_beltrami_flow_interpolation(input, output, lambda, g)
% input_flow: texture sample path
% mask: inpainting mask path
% guide: guiding image path
% output_flow: synthesized texture path
% lambda: anisotropic weight
% g: metric

function inpaint_OF_with_laplace_beltrami_flow_interpolation(flow_in,...
                                          guide, mask, lambda, w, flow_out)

% values of lambda and w are chars --> convert to double
lambda = str2double(lambda);
w = str2double(w);
                                      
kappa = im2double(imread(mask)); % read inpainting mask
kappa = kappa(:,:,1); % 
kappa(kappa~=0) = 1; % values ~= 0 are pixels to inpaint
g = im2double(imread(guide)); % read guiding image
v = readFlowFile(flow_in); % read ground truth flow
u = laplace_beltrami_interpolation(v,g,kappa,lambda,w); % inpaint missing regions
writeFlowFile(u,flow_out)