本发明涉及磁共振指纹成像技术领域，尤其涉及一种用于磁共振指纹成像的快速字典搜索方法。
背景技术：
磁共振指纹法成像(magneticresonancefingerprinting，mrf)是一种磁共振成像的新方法。不同于传统方法利用一次磁共振扫描只获取一种物理特征图像，该方法可以在一次扫描后同时获取多个物理特征图像，从而提高磁共振扫描的效率。在mrf方法的数据采集序列中，激发角度(flipangles，fa)，回波时间(echotimes，te)和重复时间(repetitiontimes，tr)这些序列参数的设置都是随机变化的，这使得所产生的磁共振信号对不同的序列参数的依赖互不相干，从而保证了不同参数组合所产生的磁共振信号的唯一性。
mrf图像重构的方法不涉及传统mr图像重构方法中的矩阵求逆的过程，因此避免了求逆过程中带来的噪音放大问题，这使得mrf方法具有扩展性，多参数同时成像，层内移动不敏感性等优点。然而在实际应用中，mrf方法面临的主要问题之一是其计算过程中消耗大量计算内存和硬盘存储空间及计算量巨大。早期mrf使用遍历搜寻的方法对信号字典进行匹配，这需要大量的内存或硬盘空间和计算时间来产生最终的图像，从而制约了mrf方法向临床医学应用的发展。
技术实现要素：
本发明的目的是克服现有磁共振指纹成像的字典搜索方法计算量巨大，计算时间长，搜索效率低的技术问题，提供了一种用于磁共振指纹成像的快速字典搜索方法，其能够有效降低搜索计算复杂度，减少计算时间，提高搜索效率。
本发明的一种用于磁共振指纹成像的快速字典搜索方法，包括以下步骤：
s1：对任意给定的mrf扫描序列及其扫描参数，通过磁共振扫描设备运行该mrf扫描序列得到每个像素位置处的实际采集到的磁共振指纹s；
s2：搜索每个像素相应的物理特征参数，搜索某个像素相应的物理特征参数包括以下步骤：将该像素的磁共振指纹信号s与其对应的信号字典d进行匹配，找到与该像素的磁共振指纹s最接近的匹配项，该最接近的匹配项就是该像素相应的物理特征参数；
s3：将所有像素的物理特征参数分别集成后即得到所需的完整的物理特征图像。
在本技术方案中，扫描序列可以是ir-bssfp序列、fisp序列或其它mrf序列。扫描参数可以包括射频激发翻转角度、tr、te、ti等。磁共振指纹s为磁共振信号的时间序列。信号字典d可预先生成或动态生成。物理特征参数可以包括t1，t2，δb等。本方法不需要预先生成完整的磁共振数据字典，具体搜索哪个像素相应的物理特征参数时，可以动态生成其对应的信号字典d。
例如：使用ir-bssfp序列(翻转恢复式的平衡稳态自由进动序列)作为扫描序列。多个单独的扫描序列片段连接在一起顺序运行，从而产生多点的磁共振指纹s。信号字典d可以在任何计算机上产生不依赖于任何磁共振设备。
在磁共振扫描设备上进行实际数据采集时所用的mrf扫描序列参数(fa，te，tr)应与信号字典d产生过程中所用的mrf扫描序列参数(fa，te，tr)完全一致。这样，在图像中每个象素对应的实际采集到的磁共振指纹s，我们都可以将其与信号字典d进行匹配，以找到最接近的匹配项，从而得到这个像素相应的物理特征参数。
作为优选，每个像素的磁共振指纹信号s对应的信号字典d是根据磁共振扫描设备运行的mrf扫描序列、扫描参数以及待量化的物理特征参数的当前匹配位置的理想值，通过数值计算方法求解bloch方程得到理想磁共振信号
作为优选，所述步骤s2中将某个像素的磁共振指纹信号s与其对应的信号字典d进行匹配之前先执行以下步骤：动态生成该像素的磁共振指纹信号s对应的信号字典d。
作为优选，所述bloch方程为：
作为优选，物理特征参数包括t1，t2，δb，信号字典d的每个字典条目包括t1，t2，δb，将某个象素对应的磁共振指纹s与其对应的信号字典d进行匹配，找到最接近的匹配项，从而得到这个像素相应的物理特征参数的方法包括以下步骤：
m1：设定t1和t2的初始值
m2：根据当前查找到的最新t1值和最新t2值，在指定δb范围根据周期tr以最低边界为起点以周期tr为间隔，将磁共振指纹s与信号字典d的字典条目依次匹配，找到当前最佳匹配点(t1，t2，δb)，以找到的当前最佳匹配点(t1，t2，δb)为中心左右距离半个周期为范围，步长为最小δb搜索分辨率的5倍再次搜索找到最新的当前最佳匹配点(t1，t2，δb)；在指定的最大搜索范围内从当前最佳匹配点(t1，t2，δb)的δb为种子点向左和向右以周期为步长遍历整个范围找到最新的当前最佳匹配点(t1，t2，δb)；以当前最佳匹配点(t1，t2，δb)为中心，左右各半个周期为范围，最小分辨率为步长进行字典条目匹配来获取或者实时生成，最终确定最佳匹配点(t1，t2，δb)，将相对应的δb值设定为最新δb值；
m3：根据当前查找到的最新δb值和最新t2值，采用聚焦变焦技术快速查找t1，包括以下步骤：
将当前t1搜索视野窗的中心设为当前的t1值，t1搜索视野窗大小为t1搜索最大范围的1/4，比较t1搜索视野窗的中心、t1搜索视野窗两个边缘点、磁共振指纹之间匹配程度大小，并将t1搜索视野窗平移使得它的中心位于这3个位置中匹配程度最高的地方，重复这个比较和平移直至t1搜索视野窗在搜索范围内不再移动，然后将t1搜索视野窗的视野减少到当前大小的1/4，重复以上搜索；
多次重复以上的缩小视野平移t1搜索视野窗的步骤直至视野大小与指定分辨率相同为止，最终t1搜索视野窗的中心点的t1值设定为最新t1值；
m4：根据当前查找到的最新δb值和最新t1值，采用步骤m3中聚焦变焦技术快速查找t2，最终t2搜索视野窗的中心点的t2值设定为最新t2值；
m5：重复执行一次步骤m3和步骤m4；
m6：最终得到的t1值、t2值、δb值为这个像素相应的物理特征参数(t1，t2，δb)。
字典搜索和匹配所需的总体数目很小，所以不需要预先将所有字典条目都产生出来，可以做到完全动态生成，需要匹配一项就产生一项。在具体搜索δb时，由于δb的搜索几乎与其它参数的搜索是独立的，所以单独对于δb可以产生一个预先生成的字典，这样匹配的时候不用动态产生字典，对于很多体素而言，使用一个小的δb字典可以进一步搜索匹配节省时间。
作为优选，物理特征参数包括t1，t2，δb，信号字典d的每个字典子集包括t1，t2，δb，将某个象素对应的磁共振指纹s与其对应的信号字典d进行匹配，找到最接近的匹配项，从而得到这个像素相应的物理特征参数的方法包括以下步骤：
m1：设定t1和t2的初始值；
m2：根据当前查找到的最新t1值和最新t2值，从信号字典d中提取一个信号字典子集d1＝d(t1，t2，δb)，这个信号字典子集d1只包含对应的t1值和t2值为当前最新值且δb可以为任意值的项，即该信号字典子集d1中的项随δb值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d1中的每一项进行比较，从信号字典子集d1中找到与磁共振指纹s匹配度最高的信号字典子集项，从而确定与该信号字典子集项相对应的δb值，将该δb值设定为最新δb值；
m3：根据当前查找到的最新δb值和最新t2值，从信号字典d中提取一个新的信号字典子集d2＝d(t1，t2，δb)，这个信号字典子集d2只包含对应的δb值和t2值为当前最新值且t1可以为任意值的项，即该信号字典子集d2中的项随t1值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d2中的每一项进行比较，从信号字典子集d2中找到与磁共振指纹s匹配度最高的信号字典子集项，从而确定与该信号字典子集项相对应的t1值，将该t1值设定为最新t1值；
m4：根据当前查找到的最新δb值和最新t1值，从信号字典d中提取一个新的信号字典子集d3＝d(t1，t2，δb)，这个信号字典子集d3只包含对应的δb值和t1值为当前最新值且t2可以为任意值的项，即该信号字典子集d3中的项随t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d3中的每一项进行比较，从信号字典子集d3中找到与磁共振指纹s匹配度最高的信号字典子集项，从而确定与该信号字典子集项相对应的t2值，将该t2值设定为最新t2值；
m5：重复执行一次步骤m3和步骤m4；
m6：重复执行一次步骤m2；
m7：根据当前查找到的最新δb值，从信号字典d中提取一个新的信号字典子集d4＝d(t1，t2，δb)，这个信号字典子集d4只包含对应的δb值为当前最新值且t1、t2可以为任意值的项，即该信号字典子集d4中的项随t1值及t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d4中的每一项进行比较，从信号字典子集d4中找到与磁共振指纹s匹配度最高的信号字典子集项，该信号字典子集项包含的t1值、t2值、δb值为这个像素相应的物理特征参数(t1，t2，δb)。
衡量实际采集到的磁共振指纹s与信号字典子集d4中的每一项的匹配度可采用相关系数cc、距离或点乘内积。
作为优选，物理特征参数包括t1，t2，δb，信号字典d的每个字典子集包括t1，t2，δb，将某个象素对应的磁共振指纹s与其对应的信号字典d进行匹配，找到最接近的匹配项，从而得到这个像素相应的物理特征参数的方法包括以下步骤：
m1：设定t1和t2的初始值；
m2：根据当前查找到的最新t1值和最新t2值，从信号字典d(t1，t2，δb)中提取一个信号字典子集d1＝d(t1，t2，δb)，这个信号字典子集d1只包含对应的t1值和t2值为当前最新值且δb可以为任意值的项，即该信号字典子集d1中的项随δb值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d1中的每一项进行相关系数计算，得到若干个相关系数cc，
相关系数cc的计算公式为：
xj代表采集到的磁共振指纹s中时间序列的某个时间点的值，yj代表信号字典子集d1中时间序列的某个时间点的值，n为时间点个数，
比较计算得到的若干个相关系数cc值，从中找出最大的相关系数cc值以及其对应的信号字典子集项，从而确定与该信号字典子集项相对应的δb值，将该δb值设定为最新δb值；
m3：根据当前查找到的最新δb值和最新t2值，从信号字典中提取一个新的信号字典子集d2＝d(t1，t2，δb)，这个信号字典子集d2只包含对应的δb值和t2值为当前最新值且t1可以为任意值的项，即该信号字典子集d2中的项随t1值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d2中的每一项进行相关系数计算，得到若干个相关系数cc，并找出最大的相关系数cc值以及其对应的信号字典子集项，从而确定与该信号字典子集项相对应的t1值，将该t1值设定为最新t1值；
m4：根据当前查找到的最新δb值和最新t1值，从信号字典中提取一个新的信号字典子集d3＝d(t1，t2，δb)，这个信号字典子集d3只包含对应的δb值和t1值为当前最新值且t2可以为任意值的项，即该信号字典子集d3中的项随t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d3中的每一项进行相关系数计算，得到若干个相关系数cc，并找出最大的相关系数cc值以及其对应的信号字典子集项，从而确定与该信号字典子集项相对应的t2值，将该t2值设定为最新t2值；
m5：重复执行一次步骤m3和步骤m4；
m6：重复执行一次步骤m2；
m7：根据当前查找到的最新δb值，从信号字典中提取一个新的信号字典子集d4＝d(t1，t2，δb)，这个信号字典子集d4只包含对应的δb值为当前最新值且t1、t2可以为任意值的项，即该信号字典子集d4中的项随t1值及t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d4中的每一项进行相关系数计算，得到若干个相关系数cc，并找出最大的相关系数cc值以及其对应的信号字典子集项，该信号字典子集项为最接近的匹配项，该信号字典子集项包含的t1值、t2值、δb值为这个像素相应的物理特征参数(t1，t2，δb)。
引入了一维及两维的快速递归搜索方法，利用搜索空间的凸性(convexity)特点来加快mrf成像方法中图像重构时的搜索算法速度。
作为优选，所述步骤m3、m4采用一维快速递归搜索方法找出最大的相关系数cc值，所述用一维快速递归搜索方法找出最大的相关系数cc值的方法包括以下步骤：
a1：在给定的搜索范围及步长的条件下，比较当前位置的cc值与最邻近左右两个位置的cc值；
a2：如果当前位置的cc值小于最邻近的某个位置的cc值，则将当前位置移到cc值最大的那个邻近位置，并重复执行步骤a2；
a3：如果当前位置的cc值不小于最邻近左右两个位置的cc值，则保持当前位置不变，并减小搜索步长；
a4：如果搜索步长大于或等于指定的最小步长，则重复执行a1至a3，如果搜索步长小于指定的最小步长，则结束步骤，所得到的当前位置所对应的变量值即为最大的相关系数cc值。
作为优选，所述步骤m7采用二维快速递归搜索方法找出最大的相关系数cc值，所述用二维快速递归搜索方法找出最大的相关系数cc值的方法包括以下步骤：
b1：在给定的搜索范围及步长的条件下，比较当前位置的cc值与最邻近的左边位置、上方位置及左上方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b2：比较当前位置的cc值与最邻近的左边位置、下方位置及左下方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b3：比较当前位置的cc值与最邻近的右边位置、上方位置及右上方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b4：比较当前位置的cc值与最邻近的右边位置、下方位置及右下方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b5：如果所有的邻近位置的cc值都不大于当前位置的cc值，则保持当前位置不变，并减小搜索步长；
b6：如果搜索步长大于或等于指定的最小步长，则重复执行b1至b5，如果搜索步长小于指定的最小步长，则结束步骤，所得到的当前位置所对应的变量值即为最大的相关系数cc值。
本发明的有益效果是：(1)基于被处理的mr参数相对独立的假设，将一个高维空间中的搜索问题简化为较低维度空间的搜索，这有效降低了mrf图像重构计算的复杂度，显著减少计算时间及占用的计算内存需求，极大地提高了计算效率，引入了一维及两维的快速递归搜索方法，利用搜索空间的凸性(convexity)特点来加快mrf成像方法中图像重构时的搜索算法速度。(2)本方法的优点是快速搜索，每个像素搜索时所需要的匹配次数非常小，所以可以不需要预先生成一个完整包含所有可能的字典，具体搜索哪个像素相应的物理特征参数时，可以动态生成其对应的信号字典d。
附图说明
图1是本发明的一种工作流程图。
具体实施方式
下面通过实施例，并结合附图，对本发明的技术方案作进一步具体的说明。
实施例1：本实施例的一种用于磁共振指纹成像的快速字典搜索方法，如图1所示，包括以下步骤：
s1：对任意给定的mrf扫描序列及其扫描参数，通过磁共振扫描设备运行该mrf扫描序列并在数据重建以后，得到每个像素位置处的实际采集到的磁共振指纹s；
s2：搜索每个像素相应的物理特征参数，搜索某个像素相应的物理特征参数包括以下步骤：动态生成该像素的磁共振指纹信号s对应的信号字典d，将该像素的磁共振指纹信号s与其对应的信号字典d进行匹配，找到与该像素的磁共振指纹s最接近的匹配项，该最接近的匹配项就是该像素相应的物理特征参数；
s3：将所有像素的物理特征参数分别集成后即得到所需的完整的物理特征图像。
每个像素的磁共振指纹信号s对应的信号字典d是根据磁共振扫描设备运行的mrf扫描序列、扫描参数以及待量化的物理特征参数的当前匹配位置的理想值，通过数值计算方法求解bloch方程得到理想磁共振信号
所述bloch方程为：
扫描序列可以是ir-bssfp序列、fisp序列或其它mrf序列。扫描参数可以包括射频激发翻转角度、tr、te、ti等。磁共振指纹s为磁共振信号的时间序列。信号字典d可预先生成或动态生成。
物理特征参数包括t1，t2，δb，信号字典d的每个字典条目包括t1，t2，δb，将某个象素对应的磁共振指纹s与其对应的信号字典d进行匹配，找到最接近的匹配项，从而得到这个像素相应的物理特征参数的方法包括以下步骤：
m1：设定t1和t2的初始值
m2：根据当前查找到的最新t1值和最新t2值，在指定δb范围(比如-300hz到300hz)根据周期tr(1/min)以最低边界(比如-300hz)为起点以周期tr为间隔，将磁共振指纹s与信号字典d的字典条目(该条目可以预先生成或动态生成)依次匹配，找到当前最佳匹配点(t1，t2，δb)，以找到的当前最佳匹配点(t1，t2，δb)为中心左右距离半个周期为范围，步长为最小δb搜索分辨率(比如1hz)的5倍再次搜索找到最新的当前最佳匹配点(t1，t2，δb)；在指定的最大搜索范围内从当前最佳匹配点(t1，t2，δb)的δb为种子点向左和向右以周期为步长遍历整个范围找到最新的当前最佳匹配点(t1，t2，δb)；以当前最佳匹配点(t1，t2，δb)为中心，左右各半个周期为范围，最小分辨率为步长进行字典条目匹配(字典条目根据t1、t2、δb)来获取或者实时生成，最终确定最佳匹配点(t1，t2，δb)，将相对应的δb值设定为最新δb值；
m3：根据当前查找到的最新δb值和最新t2值，采用聚焦变焦技术快速查找t1，包括以下步骤：
将当前t1搜索视野窗的中心设为当前的t1值，t1搜索视野窗大小为t1搜索最大范围的1/4(针对具体应用比如大脑成像，该初始值可以根据具体物体的t1作相应调整，比如设为1000msec)，比较t1搜索视野窗的中心、t1搜索视野窗两个边缘点、磁共振指纹之间匹配程度大小，并将t1搜索视野窗平移使得它的中心位于这3个位置中匹配程度最高的地方，重复这个比较和平移直至t1搜索视野窗在搜索范围内不再移动，然后将t1搜索视野窗的视野大小减少到当前大小的1/4，重复以上搜索；
多次重复以上的缩小视野平移t1搜索视野窗的步骤直至视野大小与指定分辨率相同为止，最终t1搜索视野窗的中心点的t1值设定为最新t1值；
m4：根据当前查找到的最新δb值和最新t1值，采用步骤m3中聚焦变焦技术快速查找t2，最终t2搜索视野窗的中心点的t2值设定为最新t2值；
m5：重复执行一次步骤m3和步骤m4；
m6：最终得到的t1值、t2值、δb值为这个像素相应的物理特征参数(t1，t2，δb)。
字典搜索和匹配所需的总体数目很小，所以不需要预先将所有字典条目都产生出来，可以做到完全动态生成，需要匹配一项就产生一项。在具体搜索δb时，由于δb的搜索几乎与其它参数的搜索是独立的，所以单独对于δb可以产生一个预先生成的字典，这样匹配的时候不用动态产生字典，对于很多体素而言，使用一个小的δb字典可以进一步搜索匹配节省时间。
例如：使用ir-bssfp序列(翻转恢复式的平衡稳态自由进动序列)作为扫描序列。多个单独的扫描序列片段连接在一起顺序运行，从而产生多点的磁共振指纹s。信号字典d可以在任何计算机上产生不依赖于任何磁共振设备。
在磁共振扫描设备上进行实际数据采集时所用的mrf扫描序列参数(fa，te，tr)应与信号字典d产生过程中所用的mrf扫描序列参数(fa，te，tr)完全一致。这样，在图像中每个象素对应的实际采集到的磁共振指纹s，我们都可以将其与信号字典d进行匹配，以找到最接近的匹配项，从而得到这个像素相应的物理特征参数。
本方法根据具体的快速多次扫描的序列参数产生一个磁共振信号字典d(字典容量如果超出计算机存储空间的限度时可以适当降低分辨率和缩小搜索范围)。任意抽取几个字典项目作为采集到的磁共振信号(也可以称为磁共振指纹)的理想替代品。对每个理想磁共振指纹分别计算它和所有字典项目匹配时的相似程度或匹配误差，该匹配误差相对于各字典项目所对应的生产物理参数比如t1/t2/δb/perfusion等形成一个多维函数。分析该多维匹配误差函数图像，寻找它对每一维度的近似性质并设计相应的快速搜索算法。由于该算法导致了需要匹配的字典数目急剧下降，因而在实际mrf参数估计时只需要产生非常小量的指纹字典甚至完全不需要预先生成字典。根据该技术可以预先指定一个期望获得的匹配误差函数图(比如平滑和只有一个全局最高峰等性质)反向设计扫描序列参数使得真正的匹配函数图与理想的非常接近，从而可以使用本技术方案作快速搜索。
本方法的优点是快速搜索，每个像素搜索时所需要的匹配次数非常小，所以可以不需要预先生成一个完整包含所有可能的字典，具体搜索到哪时所需要的字典项可以动态生成。当然如果预先已经生成了完备的字典，该方法照样可以使用。前者可以针对任意数目待求参数的磁共振指纹技术，也可以任意增加精度。后者(有预先生成字典的情况)在待求参数数目很大比如超过3个，精度很高时无法实现因为完备的字典没有硬盘可以存放，数据读取也极其缓慢。
实施例2：本实施例的一种用于磁共振指纹成像的快速字典搜索方法，包括以下步骤：
s1：对任意给定的mrf扫描序列及其扫描参数，通过磁共振扫描设备运行该mrf扫描序列并在数据重建以后，得到每个像素位置处的实际采集到的磁共振指纹s；
s2：搜索每个像素相应的物理特征参数(t1，t2，δb)，搜索某个像素相应的物理特征参数(t1，t2，δb)包括以下步骤：动态生成该像素的磁共振指纹信号s对应的信号字典d(t1，t2，δb)，将该像素的磁共振指纹信号s与其对应的信号字典d(t1，t2，δb)进行匹配，找到与该像素的磁共振指纹s最接近的匹配项，该最接近的匹配项就是该像素相应的物理特征参数(t1，t2，δb)；
s3：将所有像素的物理特征参数(t1，t2，δb)分别集成后即得到所需的完整的物理特征图像。
每个像素的磁共振指纹信号s对应的信号字典d(t1，t2，δb)是根据磁共振扫描设备运行的mrf扫描序列、扫描参数以及待量化的物理特征参数的当前匹配位置的理想值，通过数值计算方法求解bloch方程得到理想磁共振信号
所述bloch方程为：
将某个象素对应的磁共振指纹s与其对应的信号字典d进行匹配，找到最接近的匹配项，从而得到这个像素相应的物理特征参数的方法包括以下步骤：
m1：设定t1和t2的初始值；
m2：根据当前查找到的最新t1值和最新t2值，从信号字典d中提取一个信号字典子集d1＝d(t1，t2，δb)，这个信号字典子集d1只包含对应的t1值和t2值为当前最新值且δb可以为任意值的项，即该信号字典子集d1中的项随δb值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d1中的每一项进行比较，从信号字典子集d1中找到与磁共振指纹s匹配度最高的信号字典子集项，从而确定与该信号字典子集项相对应的δb值，将该δb值设定为最新δb值；
m3：根据当前查找到的最新δb值和最新t2值，从信号字典d中提取一个新的信号字典子集d2＝d(t1，t2，δb)，这个信号字典子集d2只包含对应的δb值和t2值为当前最新值且t1可以为任意值的项，即该信号字典子集d2中的项随t1值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d2中的每一项进行比较，从信号字典子集d2中找到与磁共振指纹s匹配度最高的信号字典子集项，从而确定与该信号字典子集项相对应的t1值，将该t1值设定为最新t1值；
m4：根据当前查找到的最新δb值和最新t1值，从信号字典d中提取一个新的信号字典子集d3＝d(t1，t2，δb)，这个信号字典子集d3只包含对应的δb值和t1值为当前最新值且t2可以为任意值的项，即该信号字典子集d3中的项随t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d3中的每一项进行比较，从信号字典子集d3中找到与磁共振指纹s匹配度最高的信号字典子集项，从而确定与该信号字典子集项相对应的t2值，将该t2值设定为最新t2值；
m5：重复执行一次步骤m3和步骤m4；
m6：重复执行一次步骤m2；
m7：根据当前查找到的最新δb值，从信号字典d中提取一个新的信号字典子集d4＝d(t1，t2，δb)，这个信号字典子集d4只包含对应的δb值为当前最新值且t1、t2可以为任意值的项，即该信号字典子集d4中的项随t1值及t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d4中的每一项进行比较，从信号字典子集d4中找到与磁共振指纹s匹配度最高的信号字典子集项，该信号字典子集项包含的t1值、t2值、δb值为这个像素相应的物理特征参数(t1，t2，δb)。
衡量实际采集到的磁共振指纹s与信号字典子集d4中的每一项的匹配度可采用相关系数cc、距离或点乘内积。
本方法的优点是快速搜索，每个像素或体素搜索时所需要的匹配次数非常小，所以可以不需要预先生成一个完整包含所有可能的字典，具体搜索到哪时所需要的字典项可以动态生成。当然如果预先已经生成了完备的字典，该方法照样可以使用。前者可以针对任意数目待求参数的磁共振指纹技术，也可以任意增加精度。后者(有预先生成字典的情况)在待求参数数目很大比如超过3个，精度很高时无法实现因为完备的字典没有硬盘可以存放，数据读取也极其缓慢。
实施例3：本实施例的一种用于磁共振指纹成像的快速字典搜索方法，包括以下步骤：
s1：对任意给定的mrf扫描序列及其扫描参数，通过磁共振扫描设备运行该mrf扫描序列并在数据重建以后，得到每个像素位置处的实际采集到的磁共振指纹s；
s2：搜索每个像素相应的物理特征参数(t1，t2，δb)，搜索某个像素相应的物理特征参数(t1，t2，δb)包括以下步骤：动态生成该像素的磁共振指纹信号s对应的信号字典d(t1，t2，δb)，将该像素的磁共振指纹信号s与其对应的信号字典d(t1，t2，δb)进行匹配，找到与该像素的磁共振指纹s最接近的匹配项，该最接近的匹配项就是该像素相应的物理特征参数(t1，t2，δb)；
s3：将所有像素的物理特征参数(t1，t2，δb)分别集成后即得到所需的完整的物理特征图像。
每个像素的磁共振指纹信号s对应的信号字典d(t1，t2，δb)是根据磁共振扫描设备运行的mrf扫描序列、扫描参数以及待量化的物理特征参数的当前匹配位置的理想值，通过数值计算方法求解bloch方程得到理想磁共振信号
所述bloch方程为：
将某个象素对应的磁共振指纹s与其对应的信号字典d进行匹配，找到最接近的匹配项，从而得到这个像素相应的物理特征参数的方法包括以下步骤：
m1：设定t1和t2的初始值；
m2：根据当前查找到的最新t1值和最新t2值，从信号字典d(t1，t2，δb)中提取一个信号字典子集d1＝d(t1，t2，δb)，这个信号字典子集d1只包含对应的t1值和t2值为当前最新值且δb可以为任意值的项，即该信号字典子集d1中的项随δb值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d1中的每一项进行相关系数计算，得到若干个相关系数cc，
相关系数cc的计算公式为：
xj代表采集到的磁共振指纹s中时间序列的某个时间点的值，yj代表信号字典子集d1中时间序列的某个时间点的值，n为时间点个数，
比较计算得到的若干个相关系数cc值，从中找出最大的相关系数cc值以及其对应的信号字典子集项，从而确定与该信号字典子集项相对应的δb值，将该δb值设定为最新δb值；
m3：根据当前查找到的最新δb值和最新t2值，从信号字典中提取一个新的信号字典子集d2＝d(t1，t2，δb)，这个信号字典子集d2只包含对应的δb值和t2值为当前最新值且t1可以为任意值的项，即该信号字典子集d2中的项随t1值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d2中的每一项进行相关系数计算，得到若干个相关系数cc，并找出最大的相关系数cc值以及其对应的信号字典子集项，从而确定与该信号字典子集项相对应的t1值，将该t1值设定为最新t1值；
m4：根据当前查找到的最新δb值和最新t1值，从信号字典中提取一个新的信号字典子集d3＝d(t1，t2，δb)，这个信号字典子集d3只包含对应的δb值和t1值为当前最新值且t2可以为任意值的项，即该信号字典子集d3中的项随t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d3中的每一项进行相关系数计算，得到若干个相关系数cc，并找出最大的相关系数cc值以及其对应的信号字典子集项，从而确定与该信号字典子集项相对应的t2值，将该t2值设定为最新t2值；
m5：重复执行一次步骤m3和步骤m4；
m6：重复执行一次步骤m2；
m7：根据当前查找到的最新δb值，从信号字典中提取一个新的信号字典子集d4＝d(t1，t2，δb)，这个信号字典子集d4只包含对应的δb值为当前最新值且t1、t2可以为任意值的项，即该信号字典子集d4中的项随t1值及t2值变化而变化，然后将实际采集到的磁共振指纹s与信号字典子集d4中的每一项进行相关系数计算，得到若干个相关系数cc，并找出最大的相关系数cc值以及其对应的信号字典子集项，该信号字典子集项为最接近的匹配项，该信号字典子集项包含的t1值、t2值、δb值为这个像素相应的物理特征参数(t1，t2，δb)。
引入了一维及两维的快速递归搜索方法，利用搜索空间的凸性(convexity)特点来加快mrf成像方法中图像重构时的搜索算法速度。
步骤m3、m4采用一维快速递归搜索方法找出最大的相关系数cc值，用一维快速递归搜索方法找出最大的相关系数cc值的方法包括以下步骤：
a1：在给定的搜索范围及步长的条件下，比较当前位置的cc值与最邻近左右两个位置的cc值；
a2：如果当前位置的cc值小于最邻近的某个位置的cc值，则将当前位置移到cc值最大的那个邻近位置，并重复执行步骤a2；
a3：如果当前位置的cc值不小于最邻近左右两个位置的cc值，则保持当前位置不变，并减小搜索步长；
a4：如果搜索步长大于或等于指定的最小步长，则重复执行a1至a3，如果搜索步长小于指定的最小步长，则结束步骤，所得到的当前位置所对应的变量值即为最大的相关系数cc值。
步骤m7采用二维快速递归搜索方法找出最大的相关系数cc值，用二维快速递归搜索方法找出最大的相关系数cc值的方法包括以下步骤：
b1：在给定的搜索范围及步长的条件下，比较当前位置的cc值与最邻近的左边位置、上方位置及左上方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b2：比较当前位置的cc值与最邻近的左边位置、下方位置及左下方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b3：比较当前位置的cc值与最邻近的右边位置、上方位置及右上方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b4：比较当前位置的cc值与最邻近的右边位置、下方位置及右下方位置的cc值，如果这些邻近位置cc值有至少一个大于当前位置的cc值，则找出其中最大cc值所在的位置，并将当前位置移到该位置，接着重复执行步骤b1；
b5：如果所有的邻近位置的cc值都不大于当前位置的cc值，则保持当前位置不变，并减小搜索步长；
b6：如果搜索步长大于或等于指定的最小步长，则重复执行b1至b5，如果搜索步长小于指定的最小步长，则结束步骤，所得到的当前位置所对应的变量值即为最大的相关系数cc值。
本方法的优点是快速搜索，每个像素搜索时所需要的匹配次数非常小，所以可以不需要预先生成一个完整包含所有可能的字典，具体搜索到哪时所需要的字典项可以动态生成。当然如果预先已经生成了完备的字典，该方法照样可以使用。前者可以针对任意数目待求参数的磁共振指纹技术，也可以任意增加精度。后者(有预先生成字典的情况)在待求参数数目很大比如超过3个，精度很高时无法实现因为完备的字典没有硬盘可以存放，数据读取也极其缓慢。