刀边法测MTF

由于审稿人要求，需要在图示中加入MTF附加图注来说明成像系统的提升，本来没有计算过MTF，一直以为MTF只有用专业的鉴别率板测量不同周期的对比度来表征，没想到只要有点物或者阶跃物品，利用刀边法就可以得到成像系统的MTF，实际上我个人觉得有些像测量CTF与OTF。

由于本人使用共焦显微镜测量的图像，所以这类扫描显微程序装置还涉及到采样规律的问题，当时没有考虑到，实际上在计算MTF刀边法的时候，应当考虑到采样问题，满足尼奎斯特采样定律，既大于截止频率二倍，简单算了一下，发现自己没有达到，不过没关系，可以利用双三次样条插值的办法，先计算插值倍率，满足采样定理后再进行计算。

步骤如下：首先当然是从图片中选取合适的‘刀边’，选择之后将这一条像素点存成一个向量;对这个向量求导数，然后做傅里叶变换，截止到截止频率的点数，就是得到的MTF曲线。

Tips:Maltab选取像素只需要在imshow后面加入ginput，在图像出来时选取需要的点回车，command窗口就会出现相应坐标值，也可以imshow图像出来之后在command窗口里敲入ginput，个人推荐第二种。（ginput默认左上角为原点且行列相反）