本发明涉及数据处理领域,具体涉及一种脑功能图的感兴趣区提取方法。
背景技术:
在各种脑功能成像后处理计算中,按照各脑区图谱提取各种参数像素值是医学生以及医学科研工作者经常遇到的问题。很多脑部疾病的发病位置不一定是全脑,有可能是某一些脑区,如颞叶,顶叶等脑区,因此在进行科研时,需要把某一个感兴趣脑区提取出来进行处理和评估等工作。
技术实现要素:
本发明针对以上目的,设计了提取某一特定脑区的方法,并计算出该感兴趣区中的脑信号均值和峰值。为医学科研工作及临床诊断提供准确的脑区信号值提取。
本发明是通过如下技术方案实现的:
一种脑功能图的感兴趣区提取方法,包括如下步骤:
步骤一)根据脑功能图建立全脑信号值矩阵epi(i,j,k),其中i,j,k为脑空间坐标;
步骤二)建立感兴趣脑区的模板矩阵roireslice(i,j,k),模板矩阵大小与全脑信号矩阵相同,roireslice(i,j,k)属于感兴趣区内的位点的值为1,不属于感兴趣区内的位点的值为0,因此模板矩阵标示出了感兴趣区的空间位置。各个脑区的模板矩阵可以从某些脑功能图处理程序中得到,也可以根据需要自行划定。
步骤三)全脑信号值矩阵epi(i,j,k)乘以感兴趣脑区的模板矩阵roireslice(i,j,k),即得到感兴趣区的信号值矩阵epi3(i,j,k);
步骤四)对感兴趣区的信号值矩阵epi3(i,j,k)中的所有不为0的元素进行求平均值,即得到感兴趣区的信号值的平均值;再逐个比较感兴趣区的信号值矩阵epi3(i,j,k)中的所有不为0的元素的大小,得到感兴趣区内的信号峰值。
所述epi(i,j,k)、roireslice(i,j,k)和epi3(i,j,k)均为91*109*91的矩阵。
与现有技术相比,本发明具有以下优点:
本方法对感兴趣区核磁信号参数的提取可以对各个脑疾病的深入研究提供详细数据,在此基础上,可以进行感兴趣区的峰值和均值计算,有助于加深对神经系统疾病的认识和研究。
具体实施方式
实施例1
本方法可以采用matlab编程实现,将从磁共振扫描中获取的各类标准化的脑功能像素值(信号值)与各脑区模板像素值一一对应(本实施例中的脑区模板采用公开软件pick_atlas软件,在mni标准脑空间坐标系中;脑区模板也可以从其他软件或材料中取得,只要能够完成对区域的定位即可),提取所需要的脑区的各类参数,如感兴趣脑区的均值,峰值等。将感兴趣的脑区模板置位1,其它脑区置为0。
epi3(i,j,k)=roireslice(i,j,k)*epi(i,j,k)(1)
其中,i,j,k是mni脑空间坐标,矩阵值为(91,109,91)。epi(i,j,k)是实验组或者正常对照组的脑部图像在配准后的全脑信号值。roireslice(i,j,k)是pick_atlas脑模板软件中感兴趣的脑区的信号值,如假设感兴趣是颞叶脑模板,则roireslice(i,j,k)在颞叶的mni空间坐标点中值为1,在非颞叶的空间坐标中值为0。感兴趣区可以选择一个或多个脑区的模板组合。
把全脑的像素值epi(i,j,k)与模板矩阵值roireslice(i,j,k)相乘(公式1),即提取了感兴趣脑区的像素值epi(i,j,k),屏蔽了非感兴趣的信号值。提取感兴趣脑区的像素值后,可进行感兴趣各类参数的计算,如求感兴趣平均值,感兴趣峰值的等计算。
程序代码如下:
roireslice=zeros(91,109,91);定义一组矩阵表示感兴趣区,并且初始化矩阵为0
epi=zeros(91,109,91);定义一组矩阵表示全脑数据,并且初始化矩阵为0
v=spm_vol_nifti('d:\results\dki-roi\angular_l.nii')
读取模版(前角脑取数据)数据赋给矩阵v
[y,xyz]=spm_read_vols(v);将矩阵v的值赋给y
fork=1:91对三维空间坐标(x,y,z)中的z轴进行循环
forj=1:109对三维空间坐标(x,y,z)中的y轴进行循环
fori=1:91对三维空间坐标(x,y,z)中的x轴进行循环
roireslice(i,j,k)=y(i,j,k);将y(i,j,k)的像素值赋给roireslice(i,j,k)
end
end
end
求感兴趣区峰值
path=strcat('d:\dki-1\dax\result\spmf_0001.img');读取脑数据
s1=analyze75info(path);读取脑数据头文件信息
epi=analyze75read(s1);读取脑数据各像素值
s2=0;
mean=0;
num=0;
max=0;
maxx=0;
maxy=0;
maxz=0;
fork=1:91
forj=1:109
fori=1:91
ifroireslice(i,j,k)>0&&epi(i,j,k)>max求感兴趣内的全脑最大值
max=epi(i,j,k);
maxx=i;保留最大像素值的x轴空间坐标
maxy=j;保留最大像素值的x轴空间坐标
maxz=k;保留最大像素值的x轴空间坐标
end
end
end
end
maxmax为最大值即峰值点的像素值
maxxmaxx为最大值即峰值点的x轴坐标
maxymaxy为最大值即峰值点的y轴坐标
maxzmaxz为最大值即峰值点的z轴坐标
求感兴趣区均值
fork=1:91
forj=1:109
fori=1:91
ifroireslice(i,j,k)>0&&epi(i,j,k)>0
计算在感兴趣区(roireslice(i,j,k)>0)内的脑像素值(epi(i,j,k)>0)的个数和像素值总和
num=num+1;
满足条件的个数加1
s2=s2+epi(i,j,k);
满足条件的像素值相加
end
end
end
end
mean=s2/num求均值mean。
通过以上程序,可以得到脑兴趣区内的均值和峰值,可以对感兴趣的特征进行进一步的分析、评价。