功能磁共振数据处理中的功率谱定量分析方法
专利汇可以提供功能磁共振数据处理中的功率谱定量分析方法专利检索，专利查询，专利分析的服务。并且本 发明 涉及 磁共振成像 技术领域，一种利用 功率谱 定量分析方法，对脑的功能磁共振数据进行处理。在采集到功能磁共振数据并经过必要的预处理后，对时间序列进行 傅立叶变换 ，然后对各频段成分的功率进行定量分析，对全脑数据逐个 体素 进行，从而获得全脑各部位涨落幅度的定量描述。步骤为：S1、脑功能磁共振数据的获取；S2、预处理；S3、傅立叶变换；S4、功率谱定量分析。该方法主要用于各种相对稳定的状态(如静息状态)以及其它非线性状态(如药物应用对脑血流的影响)等情况下的脑功能磁共振数据分析。该方法计算过程快速，在普通微机上即可完成。该方法可广泛应用于脑功能磁共振的临床与 基础 研究中。，下面是功能磁共振数据处理中的功率谱定量分析方法专利的具体信息内容。
1、一种对脑的功能磁共振数据进行处理的功率谱定量分析方法，其 特征是，利用计算机设备，采用功率谱对功能磁共振时间序列中不同频率 信号的涨落幅度进行定量分析，以达到检测脑功能的目的，功率谱定量分 析方法，是对原始数据进行快速傅立叶变换，然后对各频率的功率谱进 行定量分析，对全脑数据逐个体素进行，因而获得全脑各部位低频涨落 幅度的定量描述，以此对脑功能进行研究。
2、根据权利要求1所述的对脑的功能磁共振数据进行处理的功率谱 定量分析方法，其特征在于，对静息状态下脑功能磁共振数据的自发性 低频涨落进行定量分析，以了解静息状态下的脑功能状态。
3.根据权利要求1所述的对脑的功能磁共振数据进行处理的功率谱 定量分析方法，具体步骤如下：
步骤S1、脑功能磁共振数据的获取：脑功能磁共振数据的采集在具 备平面回波成像(EPI)序列的磁共振扫描仪上完成；
步骤S2、预处理：数据采集完毕后一般需要进行常规预处理，包括 头动校正、层面间的时间校正、去线性漂移、空间滤波、空间标准化；
步骤S3、傅立叶变换：预处理完成后，通过傅立叶变换，将时间序 列转换为功率谱，即将信号从时域转换到频域；
步骤S4、功率谱定量分析：根据信号分析理论知识可知，信号从时 域转换到频域后，对应某一频率的功率谱强度为原信号中该频率成分时 域信号振幅的平方。
4.根据权利要求3所述的对脑的功能磁共振数据进行处理的功率谱 定量分析方法，其特征在于，步骤S1、脑功能磁共振数据的获取，成像 的具体参数无特殊要求，采样时间最好在7分钟左右，重复时间选用2 秒或3秒，空间分辨率为数毫米，如3×3mm2，如果为静息状态，对受 试者的配合无特殊要求。
5.根据权利要求3所述的对脑的功能磁共振数据进行处理的功率谱 定量分析方法，其特征在于，步骤S2、预处理，基本过程完成后，对时 间序列进行带通滤波，如果准备分析低频涨落(LFF)，选择带宽为0.01- 0.08Hz之间。
6.根据权利要求3所述的对脑的功能磁共振数据进行处理的功率谱 定量分析方法，其特征在于，步骤S4、功率谱定量分析，不同受试者之 间的比较，需要进行开方，另外，为了进一步减小由于基础水平造成的 不同受试者之间的差异，对功率谱的值开方后，还需要根据全脑的均值 进行标准化。
技术领域
本发明涉及磁共振成像技术领域，一种利用功率谱定量分析方法， 对脑的功能磁共振数据进行处理。
背景技术
近年来有研究发现(Biswal等，1995)，即使在休息状态下，活体脑 内血液动力学也会出现一种低频涨落(low frequency fluctuation， LFF)，这种低频涨落反映了脑的自组织性，该研究发现这种低频涨落在 不同的脑运动皮层之间有高度的一致性。对这种慢波涨落的分析有非常 重要的意义，例如，病人无须完成复杂的任务，只需要安静地躺在磁共 振检查仪内数分钟，即有可能检查出其脑内某些功能的变化情况。但由 于这种慢波涨落不符合前面提到的多次重复事件的要求，对这种慢波涨 落的分析尚无有效的办法，目前通常采用的方法是：先让受试者完成一 定的任务，反过来再根据任务中激活的脑区定位。然后，将任务状态下 得到的脑区为依据，再研究休息状态下的慢波涨落。显然这种采用任务 状态下定位的方法并不能反映休息状态下的脑活动情况。并且，在现实 中所需要完成的任务可以说是无限多，标准也无法统一。还有研究者(Li 等，2000)采用人工方法，先在结构磁共振图像上画出海马的形状，然后 分析海马内部所有部位静息状态下功能磁共振低频涨落的一致性，发现 临床前期Alzheimer氏病患者海马内部一致性较正常对照组降低，并认 为该指标可能成为Alzheimer氏病的早期诊断指标。目前，我们尚未检 索到对低频涨落的幅度进行定量分析的报告。为此，我们提出了功率谱 分析方法，以定量描述低频涨落的幅度。
参考文献：
Biswal BB，Yetkin FZ，Haughton VM，Hyde JS.Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med 1995：34：537-541.
Li SJ，Li Z，Wu GH，Zhang MJ，Franczak M，Antuono PG.Alzheimer disease：evaluation of a functional MR imaging index as a marker. Radiology 2002；225：253-259.
发明内容
本发明的核心部分在于，将功能磁共振原始时间序列进行傅立叶变 换后，定量分析功率谱各频段的能量。对全脑数据逐个体素进行，因而 获得全脑各部位涨落幅度的定量描述，从而对脑功能进行评价。这种方 法主要用于各种相对稳定的状态(如静息状态)的数据分析以及其它非 线性状态(如药物应用对脑血流的影响)等情况下的脑功能磁共振数据分 析。
对脑的功能磁共振数据进行处理的功率谱定量分析方法，利用计算机 设备，采用功率谱对功能磁共振时间序列中不同频率信号的涨落幅度进 行定量分析，以达到检测脑功能的目的，功率谱定量分析方法，是对原 始数据进行快速傅立叶变换，然后对各频率的功率谱进行定量分析，对 全脑数据逐个体素进行，因而获得全脑各部位低频涨落幅度的定量描述， 以此对脑功能进行研究。
功率谱定量分析方法，对静息状态下脑功能磁共振数据的自发性低 频涨落进行定量分析，以了解静息状态下的脑功能状态。
该方法计算过程快速，在普通的微机上即可完成。该方法可广泛应 用于脑功能磁共振的临床与基础研究中。其实现过程可分为4个步骤， 如图1所示：
附图说明
图1是本发明的对脑的功能磁共振数据进行处理的功率谱定量分析 方法流程图。
图2是本发明的将时间序列从时域(a)转换到频域(b)图。
图3是青年组与老年组低频涨落幅度统计差异图。
具体实施方式
步骤S1、脑功能磁共振数据的获取：脑功能磁共振数据的采集在具 备平面回波成像(EPI)序列的磁共振扫描仪上完成。成像的具体参数无特 殊要求，采样时间最好在7分钟左右，重复时间(repetition time，TR) 一般选用2秒或3秒，空间分辨率一般为数毫米，如3×3mm2。如果为 静息状态，对受试者的配合无特殊要求，只需要告诉受试者安静闭目， 尽量保持头部不动，并且尽量不要思考问题。
步骤S2、预处理：数据采集完毕后一般需要进行常规预处理，包括 头动校正、层面间的时间校正、去线性漂移、空间滤波(也称为空间平滑)、 空间标准化等，但这些过程要根据使用目的来决定，并非必须的过程。 这些基本过程完成后，对时间序列进行带通滤波，如果准备分析低频涨 落(LFF)，一般选择带宽为0.01-0.08Hz之间。
步骤S3、傅立叶变换：预处理完成后，通过傅立叶变换，将时间序 列转换为功率谱，即将信号从时域转换到频域。
步骤S4、功率谱定量分析：根据信号分析理论知识可知，信号从时 域转换到频域后，对应某一频率的功率谱强度为原信号中该频率成分时 域信号振幅的平方。为便于不同受试者之间的比较，需要进行开方。另 外，为了进一步减小由于基础水平造成的不同受试者之间的差异，对功 率谱的值开方后，还需要根据全脑的均值进行标准化。至此，得到的某 一频段的值即可用于不同受试者(如病人与正常人)之间的比较。
下面举例说明本发明的实现过程：
步骤Sl、数据获取：选择16个正常青年受试者(男女各8人，年龄22-32 岁)和16个正常老年人(男女各8人，年龄60～78岁，无高血压和脑疾病)。 本研究的目的是比较正常青年人和正常老年人脑低频涨落幅度的差异， 以了解老化的生理过程。所有受试者均进行脑功能磁共振扫描，设备为 GE 1.5T，主要参数如下：TR＝2s，分辨率3.75×3.75mm，层厚6mm，扫 描时间为6分40秒。扫描时嘱受试者头部尽量保持不动，闭目，尽量不 要思考问题。
步骤S2、预处理：获取数据后，丢弃前10个扫描点(20秒)，以减 少受试者在适应扫描过程中的差异。预处理过程包括：层面间的时间校 正、头动校正、空间标准化、去线性漂移、空间平滑以及带通滤波 (0.01-0.08Hz)。
步骤S3、傅立叶变换：预处理完成后，通过傅立叶变换，将时间序 列转换为功率谱，即将信号从时域转换到频域(图2)。
步骤S4、功率谱定量分析：信号从时域转换到频域后，将频率为0.01-0.08 Hz之间的功率值开方，然后对这段开方后的数据进行平均，这样就得到 了每个体素在0.01-0.08Hz的涨落幅度。为了进一步减小由于基础水平 造成的不同受试者之间的差异，计算该频率段涨落幅度的均值，然后每 个体素的涨落幅度除以全脑均值，得到了标准化的涨落幅度。青年组与 老年组进行双样本t统计检验后，发现老年人在扣带回后部、顶叶内侧 面、辅助运动区等部位的低频涨落幅度明显减少；相反，在额顶交界处 等白质部位，老年组的低频涨落幅度却增强(图3)。
图3、青年组与老年组低频涨落幅度统计差异图。标号为1的黑色区 域(扣带回后部、顶叶内侧面、辅助运动区等部位)代表青年组高于老年 组，标号为2的白色区域(颞叶底部、额顶交界白质等部位)则代表老年 组高于青年组。统计检验的阈值为|t|＞3.646(p＜0.001)，最小团块体 积＞648mm3。
总之，休息状态下功能磁共振扫描有非常广泛的应用前景。与执行 任务的情况相比较，休息状态的扫描有如下主要优点：(1)不同检查地 点检查条件容易统一，这样才可能取得大样本数据，从而对临床诊断提 供帮助，而目前功能磁共振由于采用的任务繁多，各单位无统一标准， 不同单位间数据可比性差；(2)病人容易配合检查：执行任务情况下的 功能磁共振要求受试者完成不同难度的任务，对某些病人来说配合难度 较大，而休息状态下扫描则只要求受试者与一般的磁共振检查一样，躺 在检查仪里不动即可；(3)与休息状态相类似的状态还有睡眠、昏迷或 口服药物后(后者在持续数分钟的检查时间内是一种稳定状态)，对这些 状态的研究也有非常重要的生理及临床意义。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种行人导航定位中的自适应零速修正方法 | 2020-05-08 | 167 |
| 一种机载分布式POS的实时传递对准的方法及装置 | 2020-05-12 | 488 |
| 一种基于滤波最优平滑确定故障首达时刻的故障诊断方法 | 2020-05-16 | 371 |
| 一种基于遗传算法的视觉惯性组合的SLAM方法 | 2020-05-16 | 526 |
| 一种结合小波变换和图像分割网络的QRS波识别方法 | 2020-05-12 | 173 |
| 一种水下航行器的水下信标定位方法 | 2020-05-14 | 82 |
| 一种传感器融合采集切削信号的方法 | 2020-05-15 | 750 |
| 一种改进SPC的小区供水漏损监测预警方法 | 2020-05-15 | 678 |
| 一种柴油在线调和方法 | 2020-05-11 | 453 |
| 基于互联网的图片自动合成系统 | 2020-05-08 | 1021 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
