主动轮廓模型——Snake分割算法 matlab源码

学习图像分割算法，在网上找到的关于主动轮廓模型的实现代码，自己简化总结了一下，在这里和大家分享，欢迎提问

进入正题：

snake是一种能量最小的曲线，表示为v(s) = (x(s), y(s)), s为归一化的曲线长度，s∈[0, 1]。

能量函数由曲线内部能量和外部约束力（图像力）组成， 表示为 Esnake =  ∫Esnake(v(s))ds

                                                                                = ∫( Eint(v(s))  +  Eimage(v(s)) )ds

                                                                                      曲线内部能量        图像力         

 内部能量分为弹性势能和弯曲势能两部分：

                                                                Eint = (α(s)|vs(s)|^2 + β(s)|Vss(s)|^2）/2

                                                                              弹