MATLAB code to extract radius of curvature from FEA images