图像处理二阶微分有关问题

图像处理二阶微分问题
Two-Phase Kernel Estimation for Robust 文中提到的！


在网上找的处理二阶微分的都是laplacian operator，形式是这样的
和上面提到的这个不一样啊，请问这样的代码该怎么实现，或者用什么样的算子？？？？
在网上找的处理二阶微分的都是laplacian operator，形式是这样的
和上面提到的这个不一样啊，请问这样的代码该怎么实现，或者用什么样的算子？？？？

计算微分实际就是离散，delta足够小就可以了

比如f(x)一阶微分求解

double f1(double x)
{
     return (f(x+0.0000005)-f(x))/0.0000005;
}