createimage.m

来自「这是一个利用水平集方法进行分割的方法」· M 代码 · 共 21 行

function newim = createimage( im, phi )
% CREATEIMAGE Return segmented image as 3 channel uint8 image
%    CREATEIMAGE( im, phi ) Draws the segmented region in green
%    overtop of the original image im. The region is determined by
%    finding the front points from phi

% copy overpixels to red and blue channels of new image
newim( :, :, 1 ) = im;
newim( :, :, 3 ) = im;

% copy over pixels to a temporary image and outline pixels of the
% segmented region
tempim = im;
tempim( find( isfront( phi ) ) ) = 255;

% copy over pixels into green channel
newim( :, :, 2 ) = tempim;

% convert to uint8
newim = uint8( newim );