Matlab 水平集实现
我有一个关于级别集实现的问题 在文章“距离正则化水平集演化及其在图像分割中的应用”中,李春明,徐晨阳,你可以找到一个扩散方程: 抱歉,我不允许发布图像: 14 pdf第4页 对于dps=1-1/s,实现为[MATLAB]:Matlab 水平集实现,matlab,image-processing,image-segmentation,Matlab,Image Processing,Image Segmentation,我有一个关于级别集实现的问题 在文章“距离正则化水平集演化及其在图像分割中的应用”中,李春明,徐晨阳,你可以找到一个扩散方程: 抱歉,我不允许发布图像: 14 pdf第4页 对于dps=1-1/s,实现为[MATLAB]: distRegTerm = 4*del2(phi)-curvature; 其中: [phi_x,phi_y]=gradient(phi); s=sqrt(phi_x.^2 + phi_y.^2); smallNumber=1e-10; Nx=phi_x./(s+sma
distRegTerm = 4*del2(phi)-curvature;
其中:
[phi_x,phi_y]=gradient(phi);
s=sqrt(phi_x.^2 + phi_y.^2);
smallNumber=1e-10;
Nx=phi_x./(s+smallNumber);
Ny=phi_y./(s+smallNumber);
curvature=div(Nx,Ny);
function f = distReg_p2(phi)
[phi_x,phi_y]=gradient(phi);
s=sqrt(phi_x.^2 + phi_y.^2);
a=(s>=0) & (s<=1);
b=(s>1);
ps=a.*sin(2*pi*s)/(2*pi)+b.*(s-1);
dps=((ps~=0).*ps+(ps==0))./((s~=0).*s+(s==0));
f = div(dps.*phi_x - phi_x, dps.*phi_y - phi_y) + 4*del2(phi);
这是可以的,因为对于dp,方程是15
我不明白为什么10年了
其中:
[phi_x,phi_y]=gradient(phi);
s=sqrt(phi_x.^2 + phi_y.^2);
smallNumber=1e-10;
Nx=phi_x./(s+smallNumber);
Ny=phi_y./(s+smallNumber);
curvature=div(Nx,Ny);
function f = distReg_p2(phi)
[phi_x,phi_y]=gradient(phi);
s=sqrt(phi_x.^2 + phi_y.^2);
a=(s>=0) & (s<=1);
b=(s>1);
ps=a.*sin(2*pi*s)/(2*pi)+b.*(s-1);
dps=((ps~=0).*ps+(ps==0))./((s~=0).*s+(s==0));
f = div(dps.*phi_x - phi_x, dps.*phi_y - phi_y) + 4*del2(phi);
ps是16
代码是:
distRegTerm=distReg_p2(phi);
其中:
[phi_x,phi_y]=gradient(phi);
s=sqrt(phi_x.^2 + phi_y.^2);
smallNumber=1e-10;
Nx=phi_x./(s+smallNumber);
Ny=phi_y./(s+smallNumber);
curvature=div(Nx,Ny);
function f = distReg_p2(phi)
[phi_x,phi_y]=gradient(phi);
s=sqrt(phi_x.^2 + phi_y.^2);
a=(s>=0) & (s<=1);
b=(s>1);
ps=a.*sin(2*pi*s)/(2*pi)+b.*(s-1);
dps=((ps~=0).*ps+(ps==0))./((s~=0).*s+(s==0));
f = div(dps.*phi_x - phi_x, dps.*phi_y - phi_y) + 4*del2(phi);
我不明白这个函数的最后一行
谢谢请尝试在计算科学交换或信号处理上提问。