Add a new main file ("main_imregistration.m") to demo MATLAB embedded image registration method: good for sparse features and low quality volumetric images, but only works for very simple deformations.
 Welcome to compare this code with our proposed ALDVC code, which combines both advantages of Local Subvolume DVC (fast, compute in parallel) and Global DVC (guarantee kinematic compatibility, see https://www.mathworks.com/matlabcentral/fileexchange/77019-augmented-lagrangian-digital-volume-correlation-aldvc.
 % ===================================
% ****** ATTENTION ******
% The "x,y,z" or "1-,2-,3-" coordinates in the ALDVC code always correspond to the 1st, 2nd and 3rd indices of Matlab workspace variable. For example, p_meas(:,1) and p_meas(:,2) are the x- & y-coordinates of scattered points.
% This is a little different from some MATLAB image processing functions. % For example, if a 3D image has size MxNxL, in this code, we always have the image size_x=M, size_y=N, size_z=L. If you use some Matlab computer vision/image post-processing function, for example, 'imagesc3D', or 'imshow3D', or 'surf', it will reads size_x=N, size_y=M, size_z=L.
% Please pay attention to this difference.
 I also attach my EASF webinar to introduce AL-DIC/DVC algorithm and review local and global DIC/DVC methods: