matlab开发-FloodFill3D

matlab开发-FloodFill3D。该程序在二进制矩阵中填充三维区域。

function [B] = FloodFill3D（A slice）;
% [B] = FloodFill3D（A slice）;
% This program flood fills a 6-connected 3D region. The input matrix MUST
% be a binary image. The user will select a seed （point） in the matrix to
% initiate the flood fill. You must specify the matrix slice in which you
% wish to place the seed.
% 
% A = binary matrix
% slice = a chosen slice in the matrix where you wish to place the seed.
%
% Enjoy
% F. Dinath

A = single（A）;      % In case a logical matrix comes in.

A（1::） = NaN;     % Pad the border of the matrix
A（end::） = NaN;   % so the program doesn‘t attempt 
A（:1:） = NaN;     % to seek voxels outside the matrix
A（:end:） = NaN;   % boundry during the for loop below.
A（::1） = NaN;     %
A（::end） = NaN;   %

imagesc（A（::slice））;
title（‘select seed on figure‘）;

k = waitforbuttonpress;
point = get（gca‘CurrentPoint‘）; % button down detected
point = [fliplr（round（point（21:2））） slice];

if A（point（1） point（2） point（3））;
    A（point（1） point（2） point（3）） = NaN;
    a{1} = sub2ind（size（A） point（1） point（2） point（3））;

    i = 1;

    while 1

        i = i+1;
        a{i} = [];

        [x y z] = ind2sub（size（A） a{i-1}）;

        ob = nonzeros（（A（sub2ind（size（A） x y z-1）） == 1）.*sub2ind（size（A） x y z-1））;
        A（ob） = NaN;
        a{i} = [a{i} ob‘];

        ob = nonzeros（（A（sub2ind（size（A） x y z+1）） == 1）.*sub2ind（size（A） x y z+1））;
        A（ob） = NaN;
        a{i} = [a{i} ob‘];

        ob = nonzeros（（A（sub2ind（size（A） x-1 y z）） == 1）.*sub2ind（size（A） x-1 y z））;
        A（ob） = NaN;
        a{i} = [a{i} ob‘];

        ob = nonzeros（（A（sub2ind（size（A） x+1 y z）） == 1）.*sub2ind（size（A） x+1 y z））;
        A（ob） = NaN;
        a{i} = [a{i} ob‘];

        ob = nonzeros（（A（sub2ind（size（A） x y-1 z）） == 1）.*sub2ind（size（A） x y-1 z））;
        A（ob） = NaN;
        a{i} = [a{i} ob‘];

        ob = nonzeros（（A（sub2ind（size（A） x y+1 z）） == 1）.*sub2ind（size（A） x y+1 z））;
        A（ob） = NaN;
        a{i} = [a{i} ob‘];

        if isempty（a{i}）;
            break;
        end
%         imagesc（A（::slice））;
%         drawnow;
    end
end

b = cell2mat（a）;
b = sort（b2）;

B = logical（zeros（size（A）））;
B（b） = 1;

imagesc（B（::slice））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2245  2006-09-05 22:55  FloodFill3D.m
     文件         359  2006-09-05 22:56  FloodFill3Dexample.m
     文件      979915  2006-09-05 22:37  A.mat