GVF snake C++源码 两个版本

       今天在写算法的时候，发现在迭代的过程中，如果力场很小的话，活动轮廓会往左上角，即图像原点方向移动，想到可能是因为opencv里CvPoint是的坐标是整数类型，而在迭代的过程中，算出来的下一代坐标是double类型，在double向int转化的时候，小数点后的数都被cut了。比如轮廓上的点(10, 10)在力场的作用下移动到(9.8, 9.7)，这时候从double到int转化，结果为(9, 9)，这就会出现轮廓向左上角移动的情况。

       解决的方法有两个，一是将算法中的CvPoint都用double的类型替换，二是简单的在迭代过程中进行四舍五入。第一种方法修改代价比较高，但是结果要好于第二种方法。

=============================================================================================

开学这段时间一直在导师的指导下看GVF snake的论文，这两天参照xu的GVF snake matlab代码和crond123的源码，做了一些修改，增加了一些功能，用两种方法对图像求梯度。最终结果和xu给出的matlab结果一样。crond123的源码不能很好的收敛到凹区域，迭代次数较少，且代码有些小问题，不能成功编译。

GVF snake参数的分析这几天整理一下，再发上来。导师又给找了一些论文，接下来又要投入到endless papers当中了~~~~

http://download.csdn.net/source/2776453