matlab图像区域提取,matlab – 提取边界内的图像区域

我必须使用2D CT图像进行一个项目,并使用Matlab(仅)对肝脏和肿瘤进行分割.最初我必须单独分割肝脏区域.我使用区域生长进行肝脏分割.它获得种子点作为输入.

输出是具有肝区边界的图像.现在我需要仅由边界包围的区域.

我的程序有一个主程序和一个regionGrowing.m函数.因为我是新用户,所以不允许发布图片.如果你确实需要图像,我会给你发邮件.请帮助我.

% mainreg.m

IR=imread('nfliver5.jpg');

figure, imshow(IR), hold all

poly = regionGrowing(IR,[],15,1200); % click somewhere inside the liver

plot(poly(:,1), poly(:,2), 'LineWidth', 2, 'Color', [1 1 1])

%regionGrowing.m

function [P, J] = regionGrowing(cIM, initPos, thresVal, maxDist, tfMean, tfFillHoles, tfSimplify)

% REGIONGROWING Region growing algorithm for 2D/3D grayscale images

%

% Syntax:

% P = regionGrowing();

% P = regionGrowing(cIM);

% P = regionGrowing(cIM, initPos)

% P = regionGrowing(..., thresVal, maxDist, tfMean, tfFillHoles, tfSimpl)

% [P, J] = regionGrowing(...);

%

% Inputs:

% cIM: 2D/3D grayscale matrix {current image}

% initPos: Coordinates for initial seed position {ginput position}

% thresVal: Absolute threshold level to be included {5% of max-min}

% maxDist: Maximum distance to the initial position in [px] {Inf}

% tfMean: Updates the initial value to the region mean (slow) {false}

% tfFillHoles: Fills enclosed holes in the binary mask {true}

% tfSimplify: Reduces the number of vertices {true, if dpsimplify exists}

%

% Outputs:

% P: VxN array (with V number of vertices, N number of dimensions)

% P is the enclosing polygon for all associated pixel/voxel

% J: Binary mask (with the same size as the input image) indicating

% 1 (true) for associated pixel/voxel and 0 (false) for outside

%

% Examples:

% % 2D Example

% load example

% figure, imshow(cIM, [0 1500]), hold all