• 大小: 138KB
    文件类型: .rar
    金币: 1
    下载: 0 次
    发布日期: 2021-06-12
  • 语言: Matlab
  • 标签: 水平集  图像分割  

资源简介

完整的基于水平集的图像分割的matlab代码和注释

资源截图

代码片段和文件信息

function [kappa] = curvature(phi dx dy)
%
% function [kappa] = curvature(phi dx dy)
%
% Calculates考虑 the curvature曲率 of the function phi(xy)
%
% Author: Baris Sumengen  sumengen@ece.ucsb.edu
% http://vision.ece.ucsb.edu/~sumengen/
%



dx2 = dx*dx;
dy2 = dy*dy;

phi_xp1 = zeros(size(phi));
phi_xm1 = zeros(size(phi));
phi_yp1 = zeros(size(phi));
phi_ym1 = zeros(size(phi));

phi_xp1(:2:end) = phi(:1:end-1);
phi_xp1(:1) = 2*phi_xp1(:2)-phi_xp1(:3);
phi_xm1(:1:end-1) = phi(:2:end);
phi_xm1(:end) = 2*phi_xm1(:end-1)-phi_xm1(:end-2);

phi_yp1(2:end:) = phi(1:end-1:);
phi_yp1(1:) = 2*phi_yp1(2:)-phi_yp1(3:);
phi_ym1(1:end-1:) = phi(2:end:);
phi_ym1(end:) = 2*phi_ym1(end-1:)-phi_ym1(end-2:);

phi_x = (phi_xm1 - phi_xp1)/(2*dx);
phi_y = (phi_ym1 - phi_yp1)/(2*dy);
phi_xx = (phi_xm1 - 2*phi + phi_xp1)/dx2;
phi_yy = (phi_ym1 - 2*phi + phi_yp1)/dy2;

dummy1 = circshift(phi [11]);
dummy1(12:end) = 2*dummy1(22:end)-dummy1(32:end);
dummy1(:1) = 2*dummy1(:2)-dummy1(:3);

dummy2 = circshift(phi [-1-1]);
dummy2(end1:end-1) = 2*dummy2(end-11:end-1)-dummy2(end-21:end-1);
dummy2(:end) = 2*dummy2(:end-1)-dummy2(:end-2);

dummy3 = circshift(phi [1-1]);
dummy3(11:end-1) = 2*dummy3(21:end-1)-dummy3(31:end-1);
dummy3(:end) = 2*dummy3(:end-1)-dummy3(:end-2);

dummy4 = circshift(phi [-11]);
dummy4(end2:end) = 2*dummy4(end-12:end)-dummy4(end-22:end);
dummy4(:1) = 2*dummy4(:2)-dummy4(:3);

phi_xy = (dummy1 + dummy2 - dummy3 - dummy4)/(4*dx*dy);

abs_grad_phi_sq = (phi_x.*phi_x + phi_y.*phi_y);
abs_grad_phi_cube = abs_grad_phi_sq.^1.5;

kappa = (phi_xx.*phi_y.*phi_y - 2.*phi_y.*phi_x.*phi_xy ...
    + phi_yy.*phi_x.*phi_x)./(abs_grad_phi_cube + (abs_grad_phi_cube == 0));




 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1794  2007-09-20 13:01  LevelSetMethods\curvature.m

     文件        624  2005-05-05 11:06  LevelSetMethods\der_ENO1_minus.m

     文件        618  2005-05-06 07:07  LevelSetMethods\der_ENO1_plus.m

     文件       1250  2005-05-05 11:06  LevelSetMethods\der_ENO2_minus.m

     文件       1246  2005-05-05 11:08  LevelSetMethods\der_ENO2_plus.m

     文件       1525  2005-05-05 11:08  LevelSetMethods\der_ENO3_minus.m

     文件       1521  2005-05-05 11:08  LevelSetMethods\der_ENO3_plus.m

     文件       1602  2005-05-05 11:05  LevelSetMethods\der_WENO_minus.m

     文件       1600  2005-05-05 11:05  LevelSetMethods\der_WENO_plus.m

     文件      11010  2005-05-12 07:10  LevelSetMethods\evolve2D.m

     文件       1749  2005-05-05 11:28  LevelSetMethods\evolve_kappa.m

     文件       1463  2005-05-05 11:26  LevelSetMethods\evolve_normal_ENO1.m

     文件       1461  2005-05-05 11:27  LevelSetMethods\evolve_normal_ENO2.m

     文件       1463  2005-05-05 11:27  LevelSetMethods\evolve_normal_ENO3.m

     文件       1302  2005-05-05 11:31  LevelSetMethods\evolve_normal_vector_ENO1.m

     文件       1730  2005-05-05 11:32  LevelSetMethods\evolve_normal_vector_ENO1_SD.m

     文件       1304  2005-05-05 11:31  LevelSetMethods\evolve_normal_vector_ENO2.m

     文件       1640  2005-05-05 11:32  LevelSetMethods\evolve_normal_vector_ENO2_SD.m

     文件       1304  2005-05-05 11:30  LevelSetMethods\evolve_normal_vector_ENO3.m

     文件       1642  2005-05-05 11:32  LevelSetMethods\evolve_normal_vector_ENO3_SD.m

     文件       1302  2005-05-05 11:30  LevelSetMethods\evolve_normal_vector_WENO.m

     文件       1640  2005-05-05 11:32  LevelSetMethods\evolve_normal_vector_WENO_SD.m

     文件       1465  2005-05-05 11:30  LevelSetMethods\evolve_normal_WENO.m

     文件        674  2005-05-05 11:24  LevelSetMethods\evolve_vector_ENO1.m

     文件        672  2005-05-05 11:24  LevelSetMethods\evolve_vector_ENO2.m

     文件        672  2005-05-05 11:25  LevelSetMethods\evolve_vector_ENO3.m

     文件        673  2005-05-05 11:24  LevelSetMethods\evolve_vector_WENO.m

     文件        360  2005-05-05 11:12  LevelSetMethods\get_dt_kappa.m

     文件        356  2005-05-05 11:12  LevelSetMethods\get_dt_normal.m

     文件        399  2005-05-05 11:12  LevelSetMethods\get_dt_normal_kappa.m

............此处省略28个文件信息

评论

共有 条评论