CV模型论文公式转换代码

cv模型能量函数：

又可写成：

C1、C2公式：

代码：

mu=1；nu=0.003 * 255 * 255；epison=1；step=0.1；pi=3.14159265；

 

//Img：原图 LSF：生成图

 

Mat Drc = (epison / pi) / (epison*epison+ LSF.mul(LSF)); //Dirac 函数

//mul函数：mul会计算两个Mat矩阵对应位的乘积，所以要求参与运算的矩阵A的行列和B的行列数一致。计算结果是跟A或B行列数一致的一个Mat矩阵。

 

Mat Hea = 0.5*(1 + (2 / pi)*atan(LSF/epison)); //Heaviside 函数

//atan：自定义函数计算矩阵的反三角函数

 

Mat Length = nu*Drc.mul(cur);

 

Mat La