matlab实例之3维图像连通性分析

bwlabeln函数

bwlabeln - Label connected components in binary image

This MATLAB function returns a label matrix, L, containing labels for the
connected components in BW.

L = bwlabeln(BW)
L = bwlabeln(BW,conn)
[L,n] = bwlabeln(___)

一个matlab 实例

clear all; clc; close all;

addpath /home/dzhang/Matlab_toolbox/matlab_nifti/;
addpath /nasdisk/disk2/SleepDeprivation/voxelwise_analysis/;

CutNum = 20;
corr_type = 'pearson';

%% local FCD part
%% with 8AM
tmp = load_untouch_nii(['./RES/alphasim/lfcd_alphasim_NonVul_8am.nii']);
img = tmp.img;

BW = double(abs(img)>0);
[L,NUM] = bwlabeln(BW,6);

howmany = [];
Lcut = [];
count = 1;
for p = 1:NUM
   howmany = length(find(L(:)==p)); 
   if howmany < CutNum
      img(find(L==p)) = 0; 
   end
end

temp = load_untouch_nii('temp4.nii');
temp.hdr.dime.dim(6) = 1;
temp.img = zeros(61,73,61);
temp.img(:) = img;   
save_untouch_nii(temp,['./RES/alphasim/lfcd_alphasim_NonVul_8am_cut.nii']);