变分法水平集matlab代码

function phi = evolution(I, phi, mu, lambda, v, delta_t, epsilon, numIter,Sigma) % evolution_withoutedge(I, phi0, mu, nu, lambda_1, lambda_2, delta_t, delta_h, epsilon, numIter); % input: % I: input image % phi0: level set function to be updated % mu: weight for length term % nu: weight for area term, default value 0 % lambda_1: weight for c1 fitting term % lambda_2: weight for c2 fitting term % delta_t: time step % epsilon: parameter for computing smooth Heaviside and dirac function % numIter: number of iterations % output: % phi: updated level set function % % created on 04/26/2004 % author: Chunming Li % email: li_chunming@hotmail.com % Copyright (c) 2004-2006 by Chunming Li I=BoundMirrorExpand(I); phi=BoundMirrorExpand(phi); for k = 1 : numIter g=edge_detector_g(Sigma,I); %高斯卷积得到g % detail_phi=sec_der( phi ); %用拉普拉斯模板效果很差 detail_phi=4*del2(phi); %拉普拉斯算子 delta_phi = Delta(phi,epsilon); %狄拉克函数 [nx,ny] = curvature(phi); Curv=div(nx,ny); %曲率 div_g_phi=Sec_energy_div(delta_phi,g,nx,ny,Curv ); %div(g*▽Φ/|▽Φ|) % updating the phi function phi=phi+delta_t*(mu.*(detail_phi-Curv)+lambda.*div_g_phi+v*g.*delta_phi); end phi=BoundMirrorShrink(phi);