。. . .dbS9是相应的BOLD 时间序列，计算dbl。. . .dbS9与cU. . .cU之间的相关系数；预定时间差为0时，dbll. . .db9。是 相应的BOLD时间序列，计算dbll. ..db9。与cU. ..cU之间的相关系数；预定时间差为+1时， dbl2. . .db91是相应的BOLD时间序列，计算dbl2. . .db91与cU. . .Clni9。之间的相关系数；预定时 间差为+2时，dbl3. . .db92是相应的BOLD时间序列，计算dbl3. . .db92与cU. . .cU之间的相关 系数；预定时间差为+3时，dbl4. . .db93是相应的BOLD时间序列，计算dbl4. . .db93与dmll. . .dm9。 之间的相关系数；预定时间差为+4时，dbl5. . .db94是相应的BOLD时间序列，计算dbl5. . .db94 与cU. . .cU之间的相关系数；预定时间差为+5时，dbl6. . .db95是相应的BOLD时间序列，计 算dbl6. . .db95与dmll. . .dm9。之间的相关系数；预定时间差为+6时，dbl7. . .db96是相应的BOLD 时间序列，计算dbl7. . .db96与Clnill. . .cU。之间的相关系数；预定时间差为+7时，dbls. . .db97 是相应的BOLD时间序列，计算dbls. . .db97与cU. . .Clni9。之间的相关系数；预定时间差为+8 时，dbl9. . .db9S是相应的BOLD时间序列，计算dbl9. . .db9S与cU. . .cU之间的相关系数；预 定时间差为+9时，db2Q. . .db99是相应的BOLD时间序列，计算db2Q. . .db99与dmll. . .dm9。之间 的相关系数；预定时间差为+10时，db21. . .dbl。。是相应的BOLD时间序列，计算db21. . .dbl。。与 dnll. . .cU之间的相关系数。
[0098] 如，假定当预定时间差为0时，dbll. . .db9。与Clnill. . .cU。之间的相关系数在这些 相关系数中具有最大值，因此在时间延迟图中将体素2的选中时间差设置为0。
[0099] 同理，对体素3,计算预定时间差范围为-IOs到+IOs的时间序列与待测部位平均 时间序列（(!_...(!_)之间的相关系数。然后，从所计算出的相关系数中确定出最大值，并 将该最大值相关系数所对应的预定时间差设置为时间延迟图中体素3的选中时间差。
[0100] 具体地，对体素3,分别计算：
[0101] 预定时间差为-10时，del. ..cU是相应的BOLD时间序列，计算del. ..cU与CU. ..Clni9。之间的相关系数；预定时间差为-9时，CU. ..CU是相应的BOLD时间序列，计算 cU. . .cU与cU. . .cU之间的相关系数；预定时间差为-8时，cU. .cU是相应的BOLD时 间序列，计算cU..cU与cU. ..cU之间的相关系数；预定时间差为-7时，cU..cU是相 应的BOLD时间序列，计算cU. ..cU与cU. ..cU之间的相关系数；预定时间差为-6时， cU. ..cU3是相应的BOLD时间序列，计算cU..cU与cU. ..cU之间的相关系数；预定时 间差为-5时，cU. ..cU是相应的BOLD时间序列，计算cU. ..cU与cU. ..cU之间的相关 系数；预定时间差为-4时，cU..cU是相应的BOLD时间序列，计算cU..cU与cU. ..cU 之间的相关系数；预定时间差为-3时，cU. ..cU是相应的BOLD时间序列，计算cU. ..cU 与cU. ..cU之间的相关系数；预定时间差为-2时，cU. ..cU是相应的BOLD时间序列，计 算(1。9. . . (1。88与dmll. . .dm9。之间的相关系数；预定时间差为-1时，dd。. . . (1。89是相应的BOLD 时间序列，计算cU. ..cU与cU. ..cU之间的相关系数；预定时间差为0时，cU. ..cU是 相应的BOLD时间序列，计算cU. ..cU与cU. ..cU之间的相关系数；预定时间差为+1时， cU. ..cU是相应的BOLD时间序列，计算cU. ..cU与cU. ..cU之间的相关系数；预定时 间差为+2时，cU. ..cU是相应的BOLD时间序列，计算cU. ..cU与cU. ..cU之间的相关 系数；预定时间差为+3时，(1。14. . . (1。93是相应的BOLD时间序列，计算(1。14. . . (1。93与dmll. . .dm9。 之间的相关系数；预定时间差为+4时，cU5. . .cU4是相应的BOLD时间序列，计算cU5. . .cU4 与cU. . .cU之间的相关系数；预定时间差为+5时，cU. . .cU是相应的BOLD时间序列，计 算(1。16. . . (1。95与dmll. . .dm9。之间的相关系数；预定时间差为+6时，(1。17. . . (1。96是相应的BOLD 时间序列，计算(1。17. . .cU6与Clnill. . .cU。之间的相关系数；预定时间差为+7时，cUs. . .cU7 是相应的BOLD时间序列，计算dels. ..de97与cU. ..Clni9。之间的相关系数；预定时间差为+8 时，cU. ..cU是相应的BOLD时间序列，计算cU. ..cU与cU. ..cU之间的相关系数；预 定时间差为+9时，（1。2。. . . (1。99是相应的BOLD时间序列，计算(1。2。. . . (1。99与dmll. . .dm9。之间 的相关系数；预定时间差为+10时，cU. ..cU。。是相应的BOLD时间序列，计算cU. ..cU。。与 dnll. . .cU之间的相关系数。
[0102] 如，假定当预定时间差为+5时，dcl6. . .cU与cU. . .cU之间的相关系数在这些 相关系数中具有最大值，因此在时间延迟图中将体素3的选中时间差设置为+5。
[0103] 然后，再将时间延迟图中各个体素的选中时间差所对应的BOLD时间序列相加并 平均，以重新计算待测部位的平均BOLD时间序列。
中ta为体素1在时间延迟图中的选中时间差（即-3);其中tb为体素2在时间延迟图中 的选中时间差（即〇);其中tc为体素3在时间延迟图中的选中时间差（即+5)。
[0106] 当重新计算待测部位的平均BOLD时间序列Clmi之后，可以再次计算每个体素预定 时间差为-IOs到+IOs的数据与更新后的待测部位平均BOLD时间序列Clnin之间的相关系 数，而且可以针对每个体素，分别再从所计算出的各个相关系数中确定出最大值，并将该最 大值相关系数所对应的预定时间差设置为时间延迟图中各自体素的选中时间差，从而更新 时间延迟图。
[0107] 本发明实施方式可以通过上述迭代方式多次更新平均BOLD时间序列和时间延迟 图，直到满足预先设定的迭代停止条件即停止迭代计算。
[0108] 在计算出更新的平均BOLD时间序列和时间延迟图之后，可以将该时间延迟图中 各个体素的时间差与该时间延迟图中一预定体素的一时间差进行比较，以获取待测区域内 各个体素的血流灌注延迟时间。
[0109] 在计算出更新的平均BOLD时间序列和时间延迟图之后，可以基于更新的平均 BOLD时间序列以及时间延迟图中各个体素的选中时间差所对应的BOLD时间序列，确定各 个体素的血氧代谢程度。比如，对于每一体素，分别利用更新的平均BOLD时间序列作为回 归因子，计算该更新的平均BOLD时间序列在每个体素对应于时间延迟图中选中时间差的 各自BOLD时间序列中所占的成分，以作为每一体素的血氧代谢强弱特征。
[0110] 上述以静息态功能磁共振为实例对本发明实施方式进行了详细描述。本领域技术 人员可以意识到，本发明实施方式并不局限于以静息态功能磁共振成像系统，还可以适用 于任务态fMRI等其它功能磁共振成像系统。
[0111] 基于上述详细分析，本发明实施方式还提出了一种获取功能磁共振检测值的装 置。
[0112] 图3为根据本发明获取功能磁共振检测值的装置结构图。
[0113] 如图3所示，该装置包括：
[0114] 时间差确定单元301，用于基于一待测区域内各个体素在一预定时间点范围与一 预定时间差范围内的多个预定时间差确定的多个时间顺序所分别相应的各个BOLD时间序 列与该待测区域在所述预定时间点范围内的一平均BOLD时间序列确定各个体素的一选中 时间差，并确定该待测区域的一时间延迟图，所述时间延迟图包括各个体素的所述选中时 间差；
[0115] -检测值确定单元302,用于基于所述时间延迟图确定一功能磁共振检测值。
[0116] 在一个实施方式中，时间差确定单元301，还基于各个体素的所述选中时间差所对 应的BOLD时间序列，以迭代方式更新各个体素的所述选中时间差。而且，迭代的停止条件 包括下列中的至少一个：选中时间差不同的体素数目小于一预定值；选中时间差不同的体 素数目不变化，等等。
[0117] 在一个实施方式中，时间差确定单元301，用于基于各个体素的所述选中时间差所 对应的BOLD时间序列，以迭代方式更新各个体素的所述选中时间差。
[0118] 在一个实施方式中，时间差确定单元302,用于确定待测区域内各个体