说明~

在phi的迭代计算中，如果只更新曲线附近的phi会带来比较好的效果；否则会使得整张图片的phi值都为0。另外在更新公式中，图像力需要归一化。这是两点和原文中不同的地方。