Codebase list octave-iso2mesh / debian/1.9.6+ds-3 imedge3d.m
debian/1.9.6+ds-3

Tree @debian/1.9.6+ds-3 (Download .tar.gz)

imedge3d.m @debian/1.9.6+ds-3 — raw · history · blame 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
function imgdiff=imedge3d(binimg,isdiff)
%
% imgdiff=imedge3d(binimg,isdiff)
%
% Extract the boundary voxels from a binary image
% 
% author: Aslak Grinsted <ag at glaciology.net>
% modified by Qianqian Fang <q.fang at neu.edu>
%
% input: 
%   binimg: a 3D binary image
%   isdiff: if isdiff=1, output will be all voxels which 
%         is different from its neighbors; if isdiff=0 or 
%         ignored, output will be the edge voxels of the 
%         non-zero regions in binimg
%
% output:
%   imgdiff: a 3D logical array with the same size as binimg
%            with 1 for voxels on the boundary and 0 otherwise 
% 
% -- this function is part of iso2mesh toolbox (http://iso2mesh.sf.net)
%

invol=1;
if(nargin==2)
	invol=isdiff;
end
binimg=logical(binimg);
imgdiff=xor(binimg,binimg(:,:,[1 1:end-1]));
imgdiff=imgdiff|xor(binimg,binimg(:,:,[2:end end]));
imgdiff=imgdiff|xor(binimg,binimg(:,[1 1:end-1],:));
imgdiff=imgdiff|xor(binimg,binimg(:,[2:end end],:));
imgdiff=imgdiff|xor(binimg,binimg([1 1:end-1],:,:));
imgdiff=imgdiff|xor(binimg,binimg([2:end end],:,:));
if(invol)
	imgdiff=imgdiff&binimg;
end