function loadGroovePoints(obj)
% LOAD Load segmented surface and landmnarks
  %   Load the segmented scapula surface (if exist) and the
  %   scapula landmarks (6 scapula, 5 groove, 5 pillar) from
  %   amira directory.

  SCase = obj.shoulder.SCase;

  filename = fullfile(...
    SCase.dataMatlabPath,...
    "scapulaLandmarksAuto"+obj.shoulder.side+".mat");

  load(filename, "ScapulaLandmarks");
  obj.groove = ScapulaLandmarks.groove;
  obj.groove = obj.getSortedGrooveLateralToMedial();
end