从混在一起CT切片中分离出各套CT

最近遇到这样一个小问题，采集的CT数据，几套数据切片混杂在了一起，切片的名字是机器自动连续编号，无法判断每套数据在那个切片处分开。为此我用matlab写了段小代码，可以实现CT的分开。欢迎大家提意见。

利用matlab自带的函数

info = dicominfo(CT切片名字)

可以读取CT切片的文件头，从中可以看到CT的各个属性，如采集时间，采集者名字，像素等。用户查看info从中找到可以用于区分各套数据属性，根据该属性进行区分。我使用的代码如下：

 %在原始CT，经常会遇到几套数据的CT混在一起，本程序的作用将实现各个CT切边分开。 %读取CT头部的信息，里面ImageComments属性可以实现分开 rootDir=fullfile (pwd,'data'); %寻找data下面的图片，并加载到内存当中 clc; rootDirInfo = dir(rootDir); [rootDirNo,t]=size(rootDirInfo); fileIndex=3; while(fileIndex<=rootDirNo) fileName=rootDirInfo(fileIndex).name; testImagepath = fullfile(rootDir,fileName); info = dicominfo(testImagepath); % ClassifyProperty = info.ImageComments; %每套CT各个属性不同，至于那个属性不同，需要用户根据实际情况自己选定。 ClassifyProperty = info. PatientName.FamilyName; if ( isempty( ClassifyProperty)) sprintf('ClassifyProperty is empty,CT name is %s/n',fileName); fileIndex = fileIndex + 1; continue; end if( 7 == exist(ClassifyProperty,'dir')) %判断目录是否存在 copyfile(testImagepath,ClassifyProperty) else mkdir(ClassifyProperty); copyfile(testImagepath, ClassifyProperty); end fileIndex = fileIndex + 1; end 本文来自CSDN博客，转载请标明出处：http://blog.csdn.net/leidayihao/archive/2011/05/19/6431279.aspx 

 

 

 

混合的CT放在“data”文件下，本段代码与该文件夹同级存放。