我有一个能量方程,例如
该算法由 Lankton 实现,您可以在 code 下载代码和图像。我想基于该代码来绘制能量函数。请注意,F 是在该代码中计算的。我的目标能量数字如
我试图通过该代码实现它。但这不是正确答案
Energy=[];
%--main loop
for its = 1:max_its % Note: no automatic convergence test
%-- get the curve's narrow band
idx = find(phi <= 1.2 & phi >= -1.2)';
[y x] = ind2sub(size(phi),idx);
%-- get windows for localized statistics
xneg = x-rad; xpos = x+rad; %get subscripts for local regions
yneg = y-rad; ypos = y+rad;
xneg(xneg<1)=1; yneg(yneg<1)=1; %check bounds
xpos(xpos>dimx)=dimx; ypos(ypos>dimy)=dimy;
%-- re-initialize u,v,Ain,Aout
u=zeros(size(idx)); v=zeros(size(idx));
Ain=zeros(size(idx)); Aout=zeros(size(idx));
F_energy=zeros(size(idx));
%-- compute local stats
for i = 1:numel(idx) % for every point in the narrow band
img = I(yneg(i):ypos(i),xneg(i):xpos(i)); %sub image
P = phi(yneg(i):ypos(i),xneg(i):xpos(i)); %sub phi
upts = find(P<=0); %local interior
Ain(i) = length(upts)+eps;
u(i) = sum(img(upts))/Ain(i);
vpts = find(P>0); %local exterior
Aout(i) = length(vpts)+eps;
v(i) = sum(img(vpts))/Aout(i);
F_energy(i)=sum((img(upts)-u(i)).^2)+sum((img(vpts)-v(i)).^2); %% Compute the first term in (5) without integrate
end
%-- get image-based forces
F = -(u-v).*(2.*I(idx)-u-v);
% Compute the second term in (5)
u=phi<=0;
bw2=bwperim(u);
Length_phi=sum(sum(bw2));
Energy=[Energy (sum(F_energy(:))+alpha.*Length_phi)];
end
也许这是一项艰巨的任务,因为能量函数是如此复杂。但是,除能源项外,所有事情都由上述代码实现。希望您能理解并帮助我绘制能量函数。预先感谢
这是我的图结果。然而,它与论文结果并不相似。我的结果是接近零点的最小能量。但论文结果不是。我的代码中发生了什么。
最佳答案
您确定您的参数与原始论文中使用的参数相似吗?我观察到每次迭代中的能量至少取决于两件事:
该论文确实证实了这种关系:
和
下图显示了我使用您的代码设法实现的目标:
请注意,在原始纸张中,X 轴单位为秒。代码在每次迭代中为我们提供能量。如果不知道原始计算(在论文中描述)中一次迭代的持续时间,我们无法真正比较这些图。然而,我的更类似于原始的。
这是初始化掩码的代码(对应于绘图):
I = imread('Mushroom.png'); %-- load the image
m = false(size(I,1),size(I,2)); %-- create initial mask
m(60:100,15:80) = true; %-- initial mask coordinates
最大迭代次数:400
半径:20
希望我有所帮助。