上面程序有问题,重新上传
if abs(dp) >= abs(vdp)
hi = table(:,:,1) + table(:,:,2)./dp + table(:,:,3)./(dp.^2);
else
hi = table(:,:,4) + table(:,:,5)./vdp + table(:,:,6)./(vdp.^2);
end
hi = 1./(hi+eps);
hi(Ostu_mask(:,:,j)==0) = inf;
if abs(dp) >= abs(vdp)
Heig(:,:,j) = hi;
else
Heig(:,:,j) = max_h - hi;
end |