boyexex

parameters of this function is simple.

im_out = canny3D(im_in, size_of_gaussian_filter, SD_of_gaussian_filter, upper_threshold, lower_threshold)

im_out is a logical 3D array of surfaces; 1 for surface, 0 otherwise.

For the definition of upper_threshold and lower_threshold please check the original paper on Canny edge detection by John Canny.

Didn't work for me, but looks promising

edgevol = canny3D(double(vol),5,2,.3,.1);
??? Subscript indices must either be real positive integers or logicals.

Error in ==> interp3>linear at 374
F = (( arg4(ndx).*(1-t) + arg4(ndx+1).*t ).*(1-s) + ...

Error in ==> interp3 at 223
vi = linear(x,y,z,v,xi,yi,zi,ExtrapVal);

Error in ==> nonmax_sup at 21
imtemp = interp3(imfil_mag,xi,yi,zi);

Error in ==> canny3D at 28
im_sub = nonmax_sup(imfil_x,imfil_y,imfil_z,th_up, th_low);