www.gusucode.com > 界面看成高一维空间中某一函数ψ源码程序 > levelsetcode/binaryfit.m
function [C1,C2]= binaryfit(phi,U,epsilon) % [C1,C2]= binaryfit(phi,U,epsilon) computes c1 c2 for optimal binary fitting % input: % U: input image % phi: level set function % epsilon: parameter for computing smooth Heaviside and dirac function % output: % C1: a constant to fit the image U in the region phi>0 % C2: a constant to fit the image U in the region phi<0 % % created on 04/26/2004 % author: Chunming Li % email: li_chunming@hotmail.com % Copyright (c) 2004-2006 by Chunming Li H = Heaviside(phi,epsilon); %compute the Heaveside function values a= H.*U; numer_1=sum(a(:)); denom_1=sum(H(:)); C1 = numer_1/denom_1; b=(1-H).*U; numer_2=sum(b(:)); c=1-H; denom_2=sum(c(:)); C2 = numer_2/denom_2;