clc,clear all,close all
%% Import parameters
selpath_LA = uigetdir(['C:\Users\sachal.hussain3\OneDrive - Alma Mater Studiorum Università di Bologna\' ...
    'FATA project\Regional segmentation\Contraction parameters\Left atrium_contraction\Cardiac cycle implementation']);
LA_files = dir(fullfile(selpath_LA,'*.mat'));
%% Evaluate volumetric change
for i = 1:length(LA_files)                                          % To load all stl LA_files
    filename = fullfile(selpath_LA,LA_files(i).name);
    load(filename)
    clearvars -except mv wb selpath_LA LA_files i LA_len
    LA_len(i) = norm(minus(mv,wb))
    clear mv wb
end
plot(LA_len)

   
