一种脑功能图的感兴趣区提取方法技术

本发明专利技术公开了一种脑功能图的感兴趣区提取方法，把全脑的像素值EPI(i,j,k）与模板矩阵值ROIReslice（i,j,k）相乘（公式1），即提取了感兴趣脑区的像素值EPI (i,j,k)，屏蔽了非感兴趣的信号值。本方法对感兴趣区核磁信号参数的提取可以对各个脑疾病的深入研究提供详细数据，在此基础上，可以进行感兴趣区的峰值和均值计算，有助于加深对神经系统疾病的认识和研究。

A region of interest extraction method for brain function map

The invention discloses a brain sense extraction method of region of interest, the pixel value of whole brain EPI (I, J, K) and ROIReslice (I, J template matrix multiplication, K) (formula 1), which extracts the pixel brain region of interest in the value of EPI (I, J K), signal shielding, non interest value. This method can provide detailed data for the in-depth study of various brain diseases, and can calculate the peak value and mean value of the region of interest, which is helpful to deepen the understanding and research of nervous system diseases.

【技术实现步骤摘要】
一种脑功能图的感兴趣区提取方法
本专利技术涉及数据处理领域，具体涉及一种脑功能图的感兴趣区提取方法。
技术介绍
在各种脑功能成像后处理计算中，按照各脑区图谱提取各种参数像素值是医学生以及医学科研工作者经常遇到的问题。很多脑部疾病的发病位置不一定是全脑，有可能是某一些脑区，如颞叶，顶叶等脑区，因此在进行科研时，需要把某一个感兴趣脑区提取出来进行处理和评估等工作。
技术实现思路
本专利技术针对以上目的，设计了提取某一特定脑区的方法，并计算出该感兴趣区中的脑信号均值和峰值。为医学科研工作及临床诊断提供准确的脑区信号值提取。本专利技术是通过如下技术方案实现的：一种脑功能图的感兴趣区提取方法，包括如下步骤：步骤一）根据脑功能图建立全脑信号值矩阵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);定义一组矩阵表示感兴趣区，并且初始化矩阵为0EPI=zeros(91,109,91);定义一组矩阵表示全脑数据，并且初始化矩阵为0V=spm_vol_nifti('D:\results\DKI-ROI\Angular_L.nii')读取模版（前角脑取数据）数据赋给矩阵V[Y,XYZ]=spm_read_vols(V);将矩阵V的值赋给Yfork=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)endendend求感兴趣区峰值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:91forj=1:109fori=1:91ifROIReslice(i,j,k)>0&&EPI(i,j,k)>MAX求感兴趣内的全脑最大值MAX=EPI(i,j,k);MAXX=i;保留最大像素值的X轴空间坐标MAXY=j;保留最大像素值的X轴空间坐标MAXZ=k;保留最大像素值的X轴空间坐标endendendendMAXMAX为最大值即峰值点的像素值MAXXMAXX为最大值即峰值点的X轴坐标MAXYMAXY为最大值即峰值点的Y轴坐标MAXZMAXZ为最大值即峰值点的Z轴坐标求感兴趣区均值fork=1:91forj=1:109fori=1:91ifROIReslice(i,j,k)>0&&EPI(i,j,k)>0计算在感兴趣区（ROIReslice(i,j,k)>0）内的脑像素值(EPI(i,j,k)>0)的个数和像素值总和num=num+1；满足条件的个数加1S2=S2+EPI(i,j,k)；满足条件的像素值相加endendendendmean=S2/num求均值mean。通过以上程序，可以得到脑兴趣区内的均值和峰值，可以对感兴趣的特征进行进一步的分析、评价。本文档来自技高网...

【技术保护点】
一种脑功能图的感兴趣区提取方法，其特征在于，包括如下步骤： 步骤一）根据脑功能图建立全脑信号值矩阵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的元素的大小，得到感兴趣区内的信号峰值。

【技术特征摘要】
1.一种脑功能图的感兴趣区提取方法，其特征在于，包括如下步骤：步骤一）根据脑功能图建立全脑信号值矩阵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），即得到感兴趣区的...

【专利技术属性】
技术研发人员：王筝，曾祥柱，刘颖，袁慧书，
申请(专利权)人：北京大学第三医院，
类型：发明
国别省市：北京,11