mainVase.m

clear;
clc;
close all;
addpath(genpath('./'));

% parameters
filename              = 'vaseDiffuse';
folderData            = './data';
folderReconstruction  = './reconstructions';
load(fullfile(folderData, filename));

whetherConfocal   = true;       % whether confocal setting (colocated virtual source and detecotr)


%% ===== detecting discontinuities in transients =====
% ----- parameters -----
discontDetectionPara.expCoeff             = [0];          % model the exponential falloff of the SPAD signal
discontDetectionPara.sigmaBlur            = [0.75];       % Difference of Gaussian, standard deviation
discontDetectionPara.numOfDiscont         = 2;            % number of discontinuities per transient
discontDetectionPara.convolveTwoSides     = false;        % convolve transient with DoG filter in both sides (for detecting local minimum/maximum)
discontDetectionPara.whetherSortDisconts  = false;        % whether sort discontinuities

% ----- pathlength discontinuties visualization  -----
whetherVisualizePDSurface       = true;
whetherVisualizePDIndivisually  = false;
visualizaRange                  = 1:2000;

detectDiscontinuity;


%% ===== Fermat Flow (sphere-ray intersection) =====
% ----- parameters for computing x and y derivatives -----
pathSurfaceDerivativePara.detGridSize              = detGridSize;
pathSurfaceDerivativePara.planeFittingRange        = 5;           % local 5*5 patch for estimating x and y derivatives, odd number, at least 3
pathSurfaceDerivativePara.spatialSigma             = 8;           % bilateral filtering, spatial gaussian blur kernel size
pathSurfaceDerivativePara.diffSigma                = 10;          % bilateral filtering, range gaussian blur kernel size
pathSurfaceDerivativePara.fitErrorThreshold        = 5e-4;        % reconstruction threshold, larger the value, looser the constraint
pathSurfaceDerivativePara.curvatureRatioThreshold  = 0;           % reconstruction threshold, smaller the value, looser the constraint

% ----- x, y derivatives visualization -----
whetherVisualizexDyDzD = true;

fermatFlowReconstruction;