compute_point_match_stack_for_sample.m

script_folder_path = fileparts(mfilename('fullpath')) ;
memo_folder_path = fullfile(script_folder_path, 'memos') ;

raw_tiles_path = sprintf('/groups/mousebrainmicro/mousebrainmicro/data/%s/Tiling', sample_date) ;

% Build an index of the paths to raw tiles
raw_tile_index = compute_or_read_from_memo(memo_folder_path, ...
                                           sprintf('raw-tile-index-%s', sample_date), ...
                                           @()(build_raw_tile_index(raw_tiles_path)), ...
                                           do_force_computation) ;
tile_index_from_ijk1 = raw_tile_index.tile_index_from_ijk1 ;
raw_ijk1_from_tile_index = raw_tile_index.raw_ijk1_from_tile_index ;
ijk1_from_tile_index = raw_tile_index.ijk1_from_tile_index ;
xyz_from_tile_index = raw_tile_index.xyz_from_tile_index ;
relative_path_from_tile_index = raw_tile_index.relative_path_from_tile_index ;

% Display some features of the raw tile index
tile_lattice_shape = size(tile_index_from_ijk1)
tile_count = length(relative_path_from_tile_index) 

% Make a stack showing the number of point matches for each tile
match_count_from_ijk1_and_channel_index = ...
    compute_or_read_from_memo(memo_folder_path, ...
                              sprintf('point-match-stack-%s', sample_date), ...
                              @()(compute_point_match_stack(sample_date, tile_index_from_ijk1, ijk1_from_tile_index, relative_path_from_tile_index)), ...
                              do_force_computation) ;