如何在Matlab中实现梯度矢量流分割图像？

Question

我想应用“梯度矢量流”来检测图像中的轮廓。在我的问题中，我有一个由Canny边缘检测器提取的边缘图和一个二进制图像作为初始轮廓。初始轮廓的作用是定位活动轮廓模型。 如何实现渐变矢量流算法，以使其在Matlab中接收边缘图和二进制图像作为初始轮廓？ 谢谢，如果您能指导我解决这个问题。