医学图像转换的方法和设备的制造方法
【专利摘要】本公开涉及医学图像转换的方法，包括基于查找表通过对原始坐标系中的原始像素值进行插值来得到目标坐标系中的每个目标像素的像素值，其中所述查找表记录目标坐标系中的每个目标像素的边界在原始坐标系中的位置。本公开还涉及相应的医学图像转换的设备。
【专利说明】
医学图像转换的方法和设备
技术领域
[0001] 本公开涉及不同坐标系下的图像转换，尤其涉及可适用于非均匀分布的坐标系之 间的图像转换。
【背景技术】
[0002] 在PET(positron emission tomography，正电子发射断层成像）系统中，对应的每 一对检测器元件之间的连线可称为响应线--line-of-response(LOR)。在PET系统尤其是 环型PET系统中，真实响应线的分布是不均匀的。例如，在如图1中所示的示例性环形PET系 统中，靠近系统中心的响应线较宽，而远离系统中心的响应线较窄。这是系统的几何所决定 的。例如，如图1中所示，环形PET系统的响应线宽度可为d ? sin0，其中d是中心响应线的宽 度，而9则是当前响应线与切线之间的夹角。
[0003] 现有技术中较为常用的PET图像重建法包括迭代式统计重建法。以往，最常用的方 式通常包括先对数据通过例如插值等方式进行预处理以将其转换成均匀分布的虚拟响应 线，然后再进一步处理。然而，此类预处理会破坏数据的统计特性，诸如泊松分布。
[0004] 更新近的技术避免进行此类破坏数据统计特性的预处理，而是可以按照真实响应 线的位置信息来进行PET图像重建。例如，在Dan J Kadramas的研究论文《L0R-0SEM: statistical PET reconstruction from raw line-〇f-response histograms〉〉有详细描 述。由于不存在插值引起的误差和对数据统计特性的破坏，保持了真实响应线的信息，故而 此类方法在计算上会更加准确。
[0005] 然而，按照真实响应线来进行PET图像重建必须解决的问题之一是响应线的非均 匀分布所带来的坐标系间转换的困难。例如，当原始坐标系与目标坐标系之间至少有一者、 或甚至其两者均为非均匀分布时，在两者之间进行转换需要牵涉大量的运算和存储空间， 其处理效率和存储效率均有待提高。
【发明内容】
[0006] 本公开的一方面涉及图像转换的方法，包括基于查找表通过对原始坐标系中的原 始像素值进行插值来得到目标坐标系中的每个目标像素的像素值，其中所述查找表记录目 标坐标系中的每个目标像素的边界在原始坐标系中的位置。
[0007]在该方法方面的一示例实现中，该方法进一步包括基于所述原始坐标系和所述目 标坐标系，建立所述查找表，其中所述建立不依赖于原始图像的具体像素值，其中所述查找 表中的每个条目与目标坐标系中的目标像素对应并且包括整数部分和小数部分。
[0008]在该方法方面的另一示例实现中，所述查找表中的条目的整数部分包括相应目标 像素的左边界所落入的原始像素的索引。
[0009]在该方法方面的进一步示例实现中，所述查找表中的条目的小数部分包括相应目 标像素的左边界到所落入的原始像素的左边界的距离与该原始像素的左右边界之差的比 值。
[0010]在该方法方面的又一示例实现中，该方法还可包括计算原始坐标系中每个原始像 素与在其之前的所有原始像素各自的计数的累加，其中原始像素的计数包括该原始像素的 像素值与该原始像素的大小的乘积。
[0011]在该方法方面的一示例实现中，基于所述查找表通过对原始坐标系中的原始像素 值进行插值来得到目标坐标系中的每个目标像素的像素值进一步包括利用所计算出的原 始像素计数累加来进行所述插值。
[0012]在该方法方面的又一示例实现中，利用所计算出的原始像素计数累加来进行所述 插值进一步包括计算目标像素的左边界所落入的原始像素的左边界以及目标像素的右边 界所落入的原始像素的左边界，其中目标像素的左边界所落入的原始像素的左边界对应于 该目标像素在所述查找表中的对应条目的整数部分，并且目标像素的右边界所落入的原始 像素的左边界对应于该目标像素所述查找表中的对应条目的下一条目的整数部分；以及根 据下式计算该目标像素匕的像素值f[ i ]:
[0014] 其中Xl是目标像素5的左边界，以是原始像素幻的左边界，g[j]是原始像素值数 组，xLut [ i ]是查找表中与目标像素 f i对应的项，sumg [ j ]是对原始像素值数组g [ j ]的前j个 像素的计数进行累加所得的总和，并且id = LxLutw」.名=LxLut[i+1]」 '， 〇
[0015] 在该方法方面的又一示例实现中，该方法还包括对每个维度重复所述方法。
[0016] 在该方法方面的再一示例实现中，所述原始坐标系是S-T坐标系，并且所述目标坐 标系是X-Y坐标系，其中所述目标坐标系的X轴平行于所述原始坐标系的S轴，并且所述目标 坐标系的Y轴平行于所述原始坐标系的T轴。
[0017] 本公开的另一方面涉及用于图像转换的设备，包括配置成执行如上述的方法的处 理器；以及耦合到所述处理器用于存储数据的存储器。
【附图说明】
[0018] 图1示出了示例性环形PET系统的响应线分布。
[0019] 图2示出了根据本公开一示例性实施例的在两个不同坐标系之间进行图像转换的 示例。
[0020] 图3示出了根据本公开一示例性实施例的查找表的示图。
[0021 ]图4示出了根据本公开一示例性实施例的一维数组的坐标系间转换。
[0022] 图5示出了根据本公开一示例性实施例的用于一维数组的坐标系间转换的查找表 计算方法的流程图。
[0023] 图6示出了根据本公开一示例性实施例的二维数组的坐标系间转换。
【具体实施方式】
[0024] 将参照附图详细描述各实施例。在可能之处，相同附图标记将贯穿附图用于指代 相同或类似部分。对特定示例和附图所作的引用是用于解说性目的，而无意限定本发明或 权利要求的范围。措辞"示例性"在本文中用于表示"用作示例、实例或解说"。本文中描述为 "示例性"的任何实现不必然被解释为优于或胜过其他实现。
[0025] 本申请中以PET为例来解说实施例，但是本领域普通技术人员可以明了，本公开的 技术实质可以被应用于医学图像重建领域中的各种图像转换，而不限于PET领域。
[0026] 图2给出了根据本公开一示例性实施例的在两个不同坐标系之间进行图像转换的 示例。如图中所示，X-Y和S-T是两个不同的坐标系，并且这两个坐标系可以各自分别为均匀 分布的或非均匀分布的。另外应注意，为便于说明，图2中示出的像素尺寸仅是示意性的解 说，而并非对本公开构成任何限定。在每一坐标系中的每个像素都有一定的像素大小。例 如，在一维情况下，像素大小可对应于相应像素的长度。又如，在二维情况下，像素大小可对 应于相应像素的面积。在高维情况下的像素大小可依此类推。在一个示例性实施例中，像素 值可以是当前像素大小内的平均值，从而例如一个像素的计数可以是该像素的像素值与该 像素大小的乘积。
[0027] 同一个图像在不同坐标系下表达形式不同，即图像中每个像素值会发生变化。以 图2为例，假定原始坐标系S-T下的像素为g，目标坐标系X-Y下的像素为f。根据一示例性实 施例，要将原始坐标系S-T下沿S轴的一组像素 gi，g2，g3…转换成目标坐标系X-Y下沿X轴的 一组像素…。例如，对于X-Y坐标系下的像素 f2,它的像素值的计算公式可以表示为 例如：
(1)
[0029] 其中，8[1]4[2]4[3]，玎1]，打2]，打3]为各个单位像素即8142,8 3彳1，5彳3的像 素值，11，12,81，82分别表示像素€1，€2^1^2的右边界。在该示例中，乂-￥坐标系下的像素€2的 值的计算涉及S-T坐标系下的三个像素81^2^3。本领域普通技术人员应理解，这仅是示例 性实施例。取决于具体的坐标系的分布以及具体目标像素，计算中所涉及的原始像素的数 目将可能有所不同。
[0030] 现有技术算法中，解决该问题的直接方式是首先找出与目标像素(例如f2)相关的 所有像素(例如，在图2的示例中为 gl，g2，g3)，然后按照上述公式进行计算。然而，这种方式 的缺点是:计算过程中需要的存储空间和运算效率与原始坐标系及目标坐标系中每个像素 的尺寸都相关，当S-T坐标系和X-Y坐标系中像素尺寸差距较大时，会增大存储空间开销，同 时运算速度显著下降。
[0031] 为了解决现有技术中的上述及其他问题，本公开提出了改进方案。根据本方案，当 图像在两个任意给定的平行坐标系间进行转换时，运算效率仅与目标坐标系相关，而计算 过程中额外需要的存储空间仅与原始坐标系相关。本公开的坐标系转换方案适用于一维、 二维、三维和更高维坐标系，因为根据至少一个示例性实施例，图像转换过程是按照维度依 次进行的，每次可以只对一个维度下的图像进行转换。
[0032] 下面以二维坐标系中的一维数组插值为例，原始坐标系S-T下，沿S轴的一维原始 像素值数组记为8[」]，」=〇，1，'"，对应的像素」的右边界坐标记为8 1，幻，'"，以+1;在目标坐 标系X-Y下，沿X轴的一维目标像素值数组记为f [i ]，i = 0，1，…，对应的像素 f:的右边界坐 标记为X1，X2,…，Xl+1。在一示例性而非限定性的实施例中，从原始坐标系S-T下沿S轴的原始 像素值数组g[ j]向目标坐标系X-Y下沿X轴的目标像素6的像素值f[i]的示例性转换过程 可以包括如下所述的处理中的至少一部分。
[0033] 首先，计算查找表xLut。查找表xLut用于记录目标坐标系X-Y中沿X轴的每一个目 标像素6的右边界在原始坐标系S-T中沿S轴的位置。对于目标坐标系X-Y中沿X轴的第i个 目标像素6，它在查找表xLut中对应的值可以包括以下两个部分：（i)整数部分，其可记录 在原始坐标系S-T中沿S轴的符合下述条件的一个原始像素 gj:该原始像素的左边界小于 等于目标像素匕的左边界且该原始像素幻的右边界大于目标像素匕的左边界，换言之，原始 像素^所索引的是覆盖目标像素 h的左边界的那一个原始像素(例如，在图2的示例中，对于 目标像素 f2,所找到的符合条件的原始像素是gl);(ii)小数部分，其可记录目标像素 h的左 边界到原始像素^的左边界的距离与原始像素^的宽度（即，该原始像素左右边界之差）的 比值。本领域普通技术人员将理解，对"左"/"右"边界的命名仅是为了方便起见而非构成限 定。例如，对于纵轴方向而言，所称的"左"/"右"边界分别对应于上/下边界。不失一般性， "左"边界可指像素在所议维度的轴上较小的边界，而"右"边界可指像素在所议维度的轴上 较大的边界。但本公开并不受此限定。
[0034]从上述查找表xLut的定义可见，它只与原始和目标坐标系分别的定义有关，而与 具体的图像无关。所以，只要原始和目标坐标系确定了，便可以事先计算出xLut保存于硬盘 或者内存中，需要的时候直接调用即可。因此，只需事先计算一次xLut，此后便可将其重用 于具体的图像。另外，从上述查找表xLut的定义可见，它所占用的存储空间只与目标坐标系 中的像素个数有关，而与原始坐标系无关，也与原始坐标系及目标坐标系中每个像素的尺 寸无关，并且即便原始坐标系S-T与目标坐标系X-Y中像素尺寸差距较大也不会受到影响。
[0035] 图3示出了根据本公开一示例性实施例的查找表的示图。
[0036] 图4示出了根据本公开一示例性实施例的一维数组的坐标系间转换，并且图5示出 了根据本公开一示例性实施例的用于一维数组的坐标系间转换的查找表计算方法100的流 程图。该示例性实施例以原始坐标系S向目标坐标系X的转换为例来解说，但本领域普通技 术人员应理解，该示例性实施例同样能被现成地应用于其他维度（例如，一维、三维、更高 维）的坐标系之间的相应轴之间的转换。另外，在该示例性实施例中，为简化解说起见，使原 始坐标系中的首个原始像素与目标坐标系中的首个目标像素的左边界对齐(即， XQ = S())。但 本领域普通技术人员能无需任何创造性劳动地将本公开的算法稍加改动即应用于原始坐 标系中的首个原始像素与目标坐标系中的首个目标像素的左边界不对齐的情形。
[0037] 根据该示例性而非限定性的实施例，查找表xLut的计算过程100可包括如下步骤。
[0038] 框 101:初始化 i = 0, j = 0，xo = so，xLut[i] =0,0<;1〈;[111￡?，其中：[111￡?-1是目标像素值 数组f[ i ]中最后一个像素的序号；
[0039] 框102:根据下式判断是否已经遍历所有像素：i〈imax，当该式成立时，表明尚未遍 历所有像素，该流程转至框103，而当该式不成立时，表明已遍历完所有像素，该流程转至框 106；
[0040]框103:判断目标坐标系X中的当前目标像素 fi的右边界xi+1与原始坐标系S中的像 素 gj的左边界Sj和右边界8」+ 1的关系是否满足:Sj<Xi+1〈Sj+1，如果是，则该流程转到框104， 如果否，则该流程转到框105;
[0041] 框104
[0042] 框 i〇5:j = j+l;
[0043] 框106:当已遍历完所有像素时，计算结束，输出查找表xLut。
[0044]然后，可利用事先计算并存储的查找表xLut来进行图像重建。例如，在用迭代算法 进行图像重建的过程中，在对当前图像进行一次计算和更新后，还需要将更新后的图像作 为下一次运算的初始图像，这样反复进行多次才能获得最终图像。相应地，对于图像在不同 坐标系下的转换也要反复进行。每次图像转换包括例如以下两个步骤 [0045]第一，将原始坐标系S-T下的沿S轴的一维原始像素值数组g[k](k = 0，l，-_)的前j 个像素的计数进行累加，得到数组sumg[ j ]，j = 0，1，…，即
(2)
[0047]第二，通过插值得到目标坐标系X-Y下的沿X轴的的每个像素匕的像素值f [i]，即：
(3)
[0049] 其中
[0050] A = LxLut[i]」；乂 =|_xLut[i+1]」，其中L」表示向下取整。
[0051] 本公开的方法不仅适用于一维数组的插值，也适用于更高维度数组。对于一个n维 数组，对第1，2, ..，n维依次进行上述过程即可。图6示出了根据本公开一示例性实施例的二 维数组的坐标系间转换。
[0052]以二维数组为例，根据一个示例性实施例从S-T坐标系转换到X-Y坐标系的过程如 下：
[0053] (1)建立查找表xLut，并且对S-T坐标系中的S坐标进行插值，将二维数组g[s，t]由 S-T坐标系转换得到X-T坐标系，产生中间结果h[x，t];
[0054] (2)建立查找表yLut，并且对X-T坐标系中的T坐标进行插值，将二维数组h[x，t]由 X-T坐标系转换得到X-Y坐标系，产生最终结果f[x，y]。
[0055]本领域技术人员可知，对多维坐标系的转换而言，所转换的维度的次序并不受限 制。例如，在上述二维坐标系的转换中，也可以先对T坐标进行插值转换到中间的S-Y坐标 系，再对S坐标进行插值转换到X-Y坐标系。
[0056]本领域普通技术人员应理解，本公开的有益效果并非由任何单个实施例来全部实 现。各种组合、修改和替换均为本领域普通技术人员在本公开的基础上所明了。
[0057]此外，术语"或"旨在表示包含性"或"而非排他性"或"。即，除非另外指明或从上下 文能清楚地看出，否则短语"X采用A或B"旨在表示任何自然的可兼排列。即，短语"X采用A或 B"箱由以下实例中任何实例得到满足:X米用A;X米用B;或X米用A和B两者。另外，本申请和 所附权利要求书中所用的冠词"一"和"某"一般应当被理解成表示"一个或多个"，除非另外 声明或者可从上下文中清楚看出是指单数形式。
[0058] 各个方面或特征将以可包括数个设备、组件、模块、及类似物的系统的形式来呈 现。应理解和领会，各种系统可包括附加设备、组件、模块等，和/或可以并不包括结合附图 所讨论的全部设备、组件、模块等。也可以使用这些办法的组合。
[0059] 结合本文所公开的实施例描述的各种说明性逻辑、逻辑块、模块、和电路可用通用 处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或其它可 编程逻辑器件、分立的门或晶体管逻辑、分立的硬件组件、或其设计成执行本文所描述功能 的任何组合来实现或执行。通用处理器可以是微处理器，但在替换方案中，处理器可以是任 何常规的处理器、控制器、微控制器、或状态机。处理器还可以被实现为计算设备的组合，例 如DSP与微处理器的组合、多个微处理器、与DSP核心协同的一个或多个微处理器、或任何其 它此类配置。此外，至少一个处理器可包括可作用于执行以上描述的一个或多个步骤和/或 动作的一个或多个模块。例如，在一示例性实施例中，一种医学图像转换的设备可包括配置 成执行如本文中所述的任何方法的处理器；以及耦合到所述处理器用于存储数据的存储 器。
[0060] 此外，结合本文中所公开的方面描述的方法或算法的步骤和/或动作可直接在硬 件中、在由处理器执行的软件模块中、或在这两者的组合中实施。
[0061] 本公开中通篇描述的各种方面的要素为本领域普通技术人员当前或今后所知的 所有结构上和功能上的等效方案通过引述被明确纳入于此，且意在被权利要求书所涵盖。 此外，本文所公开的任何内容都并非旨在贡献给公众一一无论这样的公开是否在权利要求 书中被显式地叙述。
【主权项】
1. 一种图像转换的方法，其特征在于，包括： 基于查找表通过对原始坐标系中的原始像素值进行插值来得到目标坐标系中的每个 目标像素的像素值，其中所述查找表记录目标坐标系中的每个目标像素的边界在原始坐标 系中的位置。2. 如权利要求1所述的方法，其特征在于，进一步包括： 基于所述原始坐标系和所述目标坐标系，建立所述查找表，其中所述建立不依赖于原 始图像的具体像素值，其中 所述查找表中的每个条目与目标坐标系中的目标像素对应并且包括整数部分和小数 部分。3. 如权利要求2所述的方法，其特征在于，所述查找表中的条目的整数部分包括相应目 标像素的左边界所落入的原始像素的索引。4. 如权利要求2所述的方法，其特征在于，所述查找表中的条目的小数部分包括相应目 标像素的左边界到所落入的原始像素的左边界的距离与该原始像素的左右边界之差的比 值。5. 如权利要求1所述的方法，其特征在于，进一步包括： 计算原始坐标系中每个原始像素与在其之前的所有原始像素各自的计数的累加，其中 原始像素的计数包括该原始像素的像素值与该原始像素的大小的乘积。6. 如权利要求5所述的方法，其特征在于，基于所述查找表通过对原始坐标系中的原始 像素值进行插值来得到目标坐标系中的每个目标像素的像素值进一步包括： 利用所计算出的原始像素计数累加来进行所述插值。7. 如权利要求6所述的方法，其特征在于，利用所计算出的原始像素计数累加来进行所 述插值进一步包括： 计算目标像素的左边界所落入的原始像素的左边界W及目标像素的右边界所落入的 原始像素的左边界，其中目标像素的左边界所落入的原始像素的左边界对应于该目标像素 在所述查找表中的对应条目的整数部分，并且目标像素的右边界所落入的原始像素的左边 界对应于该目标像素所述查找表中的对应条目的下一条目的整数部分；W及 根据下式计算该目标像素的像素值f[ i ]其中Xi是目标像素fi的左边界，Sj是原始像素gj的左边界，g[j]是原始像素值数组，xLut [i]是查找表中与目标像素fi对应的项，sumg[j]是对原始像素值数组g[j]的前j个像素的 计数进行累加所得的总和，并i8. 如权利要求1所述的方法，其特祉在十，进一步包巧对每个维度重复所述方法。9. 如权利要求1所述的方法，其特征在于，所述原始坐标系是S-T坐标系，并且所述目标 坐标系是X-Y坐标系，其中所述目标坐标系的X轴平行于所述原始坐标系的巧由，并且所述目 标坐标系的Y轴平行于所述原始坐标系的T轴。10.-种用于图像转换的设备，其特征在于，包括： 配置成执行如权利要求1-9中任一项所述的方法的处理器；W及 禪合到所述处理器用于存储数据的存储器。
【文档编号】G06T3/60GK105913380SQ201610200285
【公开日】2016年8月31日
【申请日】2016年3月31日
【发明人】吕杨
【申请人】上海联影医疗科技有限公司