%% Identify Vertebrae and Show Deviations
[x, y] = ginput;
plot(x, y, 'ro', 'LineWidth', 2, 'MarkerSize', 15);
points = [x, y];
points = round(points);
xCenter = points(find(y==max(y)), 1);
xNew = x - xCenter;
subplot(1, 2, 2);
stem(y, xNew);
camroll(-90);
yticks([]);
xticks([]);
xlim([min(y)-50 max(y)+50]);
if max(abs(xNew)) >= min(abs(xNew))
   ylim([-max(abs(xNew))-50 max(abs(xNew))+50]);
else
   ylim([-min(abs(xNew))-50 min(abs(xNew))+50]);
end
title('Deviations From Normal Spine');
pause(2.0);