As seen in the paper entitled 2D peripheral image quality metrics with different types of multifocal contact lenses. https://www.nature.com/articles/s41598-019-54783-x
Use as FIG_EstimationSF() to obtain the same image as the one printed in the article.
You can modify the image used in the subplot b, by calling to the function FIG_EstimationSF(imageDir), input should be readable by imread function in matlab.
If you want to add values to the function to calculate the SFcutoff limit, look at the function estimateSFCutOFF.
help FIG_EstimationSF
-----