一种医学影像对齐方法、介质及电子设备与流程

1.本发明属于图像处理领域，涉及一种图像对齐方法，特别是涉及一种医学影像对齐方法、介质及电子设备。


背景技术：

2.ct(computed tomography，电子计算机断层扫描)检查是现代一种较先进的医学影像检查技术。在ct检查过程中往往需要拍摄同一病人的平扫期ct影像、静脉期ct影像和动脉期ct影像等。其中，平扫期ct影像是指未注射造影剂时获得的ct影像，静脉期ct影像是指静脉血管充盈显影时刻获得的ct影像，动脉期ct影像是指动脉血管充盈显影时刻获得的ct影像。然而，实际检查过程中，各期ct影像的拍摄时间可能并不相同，病人的呼吸和移动等因素会导致不同拍摄时间下病人在ct扫描机中的位置存在偏差，因而造成不同期ct影像在垂直方向上位置的不对应，例如，可能出现平扫期的第一张ct影像和动脉期的第二张ct影像相对应，或平扫期的第三张ct影像和动脉期的第一张ct影像相对应。现有技术中往往需要人工方式实现不同期ct影像的对齐，效率较低。


技术实现要素：

3.鉴于以上所述现有技术的缺点，本发明的目的在于提供一种医学影像对齐方法、介质及电子设备，用于解决现有技术采用人工对齐方式效率较低的问题。
4.为实现上述目的及其他相关目的，本发明的第一方面提供一种医学影像对齐方法，用于将目标影像与参考影像进行对齐。所述医学影像对齐方法包括：获取所述目标影像中参考对象的特征值，并根据所述目标影像中所述参考对象的特征值获得所述目标影像的特征序列；获取所述参考影像中所述参考对象的特征值，并根据所述参考影像中所述参考对象的特征值获得所述参考影像的特征序列；根据所述目标影像的特征序列和所述参考影像的特征序列获取平移距离；根据所述平移距离对所述目标影像进行平移。
5.于所述第一方面的某些实施例中，获取所述目标影像中参考对象的特征值的实现方法包括：获取所述参考对象的ct值范围；根据所述ct值范围获取所述参考对象的像素点，进而生成所述参考对象的模板；根据所述参考对象的蒙版获取所述参考对象的特征值。
6.于所述第一方面的某些实施例中，所述参考对象的特征值包括：所述参考对象对应凸包的几何特征、所述参考对象的中心点和/或所述参考对象的面积。
7.于所述第一方面的某些实施例中，根据所述目标影像中参考对象的特征值获得所述目标影像的特征序列的实现方法包括：根据所述目标影像中参考对象的特征值获得所述目标影像对应的特征向量；根据所述目标影像对应的特征向量获得所述目标影像的特征序列。
8.于所述第一方面的某些实施例中，根据所述目标影像的特征序列和所述参考影像的特征序列获取平移距离的实现方法包括：获取所述平移距离的取值范围；所述取值范围包括至少两个整数值；根据所述目标影像的特征序列和所述参考影像的特征序列，依次计
算所述取值范围内的各整数值对应的序列方差；根据所述序列方差选取相应的整数值作为所述平移距离。
9.于所述第一方面的某些实施例中，根据所述目标影像的特征序列和所述参考影像的特征序列计算与整数值x对应的序列方差的公式为：
[0010][0011]
其中，n为所述目标影像包含的ct影像数量，m为所述目标影像中各ct影像对应的参考对象的特征值数量；参考对象的特征值数量；表示中第i行第j列的特征值，表示中第i行第j列的特征值。
[0012]
于所述第一方面的某些实施例中，所述参考对象为脊柱。
[0013]
于所述第一方面的某些实施例中，所述目标影像包括：平扫期ct影像、静脉期ct影像和/或动脉期ct影像。
[0014]
本发明的第二方面还提供一种计算机可读存储介质，其上存储有计算机可读程序。所述计算机程序被处理器执行时实现本发明所述医学影像对齐方法。
[0015]
本发明的第三方面还提供一种电子设备，所述电子设备包括：存储器，存储有一计算机程序；处理器，与所述存储器通信相连，调用所述计算机程序时执行本发明所述医学影像对齐方法；显示器，与所述处理器和所述存储器通信相连，用于显示所述医学影像对齐方法的相关gui交互界面。
[0016]
如上所述，本发明所述的医学影像对齐方法、介质及电子设备，具有以下有益效果：
[0017]
所述医学影像对齐方法能够自动获取目标影像和参考影像的特征序列，并依据所述特征序列获取目标影像的平移距离进而对所述目标影像进行平移以实现所述目标影像与所述参考影像相对齐。整个过程基本无需人工参与，操作方便且效率较高。
附图说明
[0018]
图1a显示为本发明所述医学影像对齐方法于一具体实施例中的平扫期ct影像示例图。
[0019]
图1b显示为本发明所述医学影像对齐方法于一具体实施例中的静脉期ct影像示例图。
[0020]
图1c显示为本发明所述医学影像对齐方法于一具体实施例中的动脉期ct影像示例图。
[0021]
图2显示为本发明所述医学影像对齐方法于一具体实施例中的流程图。
[0022]
图3a显示为本发明所述医学影像对齐方法于一具体实施例中目标影像包含的ct序列示例图。
[0023]
图3b显示为本发明所述医学影像对齐方法于一具体实施例中参考影像包含的ct序列示例图。
[0024]
图4a显示为本发明所述医学影像对齐方法于一具体实施例中获取特征值的流程图。
[0025]
图4b显示为本发明所述医学影像对齐方法于一具体实施例中获取参考对象蒙版的流程图。
[0026]
图4c显示为本发明所述医学影像对齐方法于一具体实施例中的ct图像示例图。
[0027]
图4d显示为本发明所述医学影像对齐方法于一具体实施例中的ct图像示例图。
[0028]
图4e显示为本发明所述医学影像对齐方法于一具体实施例中的ct图像示例图。
[0029]
图4f显示为本发明所述医学影像对齐方法于一具体实施例中的ct图像示例图。
[0030]
图5显示为本发明所述医学影像对齐方法于一具体实施例中获取特征序列的流程图。
[0031]
图6显示为本发明所述医学影像对齐方法于一具体实施例中步骤s23的流程图。
[0032]
图7a显示为本发明所述医学影像对齐方法于又一具体实施例中的流程图。
[0033]
图7b显示为本发明所述医学影像对齐方法于一具体实施例中步骤s72的流程图。
[0034]
图7c显示为本发明所述医学影像对齐方法于一具体实施例中步骤s74的流程图。
[0035]
图8显示为本发明所述电子设备于一具体实施例中的结构示意图。
[0036]
元件标号说明
[0037]4ꢀꢀꢀꢀꢀꢀꢀꢀꢀ
选择框
[0038]
411
ꢀꢀꢀꢀꢀꢀꢀ
上边框
[0039]
412
ꢀꢀꢀꢀꢀꢀꢀ
左边框
[0040]
413
ꢀꢀꢀꢀꢀꢀꢀ
右边框
[0041]
414
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
下边框
[0042]4’ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
外边框
[0043]
411
’ꢀꢀꢀꢀꢀꢀꢀꢀ
调整后的上边框
[0044]
412
’ꢀꢀꢀꢀꢀꢀꢀꢀ
调整后的左边框
[0045]
413
’ꢀꢀꢀꢀꢀꢀꢀꢀ
调整后的右边框
[0046]
414
’ꢀꢀꢀꢀꢀꢀꢀꢀ
调整后的下边框
[0047]
51
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
脊柱
[0048]
52
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
肋骨
[0049]
800
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
电子设备
[0050]
810
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
存储器
[0051]
820
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
处理器
[0052]
830
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
显示器
[0053]
s21～s24
ꢀꢀꢀꢀꢀ
步骤
[0054]
s211a～s213b 步骤
[0055]
s211b～s212b 步骤
[0056]
s231～s233
ꢀꢀꢀ
步骤
[0057]
s41～s44
ꢀꢀꢀꢀꢀ
步骤
[0058]
s71～s74
ꢀꢀꢀꢀꢀ
步骤
[0059]
s721～s724
ꢀꢀꢀ
步骤
[0060]
s741～s743
ꢀꢀꢀ
步骤
具体实施方式
[0061]
以下通过特定的具体实例说明本发明的实施方式，本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用，本说明书中的各项细节也可以基于不同观点与应用，在没有背离本发明的精神下进行各种修饰或改变。需说明的是，在不冲突的情况下，以下实施例及实施例中的特征可以相互组合。
[0062]
需要说明的是，以下实施例中所提供的图示仅以示意方式说明本发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其组件布局型态也可能更为复杂。
[0063]
实际ct检查过程中，不同期ct影像的拍摄时间可能并不相同，病人的呼吸和移动等因素会导致不同拍摄时间下病人在ct扫描机中的位置存在偏差，因而造成不同期ct影像在垂直方向上位置的不对应。请参阅图1a～图1c，显示为一实施例中获得的不同期ct影像示例图；其中，图1a显示为平扫期拍摄的第一张ct影像，图1b显示为静脉期拍摄的第一张ct影像，图1c显示为动脉期拍摄的第一张ct影像。可以看出，图1a所示平扫期ct影像中器官和骨骼的位置与形状与图1b所示ct影像中器官和骨骼的位置与形状基本一致，因此图1a所示ct影像与图1b所示ct影像相对齐；图1a所示平扫期ct影像中的器官和骨骼的位置与形状与图1c所示动脉期ct影像中的器官和骨骼的位置与形状相差较大，故图1a所示ct影像与图1c所示ct影像不对齐，此时需要对平扫期ct影像和动脉期ct影像进行对齐处理。现有技术往往采用人工方式实现不同期ct影像的对齐，效率较低。
[0064]
针对这一问题，本发明提供一种医学影像对齐方法，用于将目标影像与参考影像进行对齐。所述医学影像对齐方法包括：获取所述目标影像中参考对象的特征值，并根据所述目标影像中参考对象的特征值获得所述目标影像的特征序列；获取所述参考影像中参考对象的特征值，并根据所述参考影像中参考对象的特征值获得所述参考影像的特征序列；根据所述目标影像的特征序列和所述参考影像的特征序列获取平移距离；根据所述平移距离对所述目标影像进行平移。所述医学影像对齐方法能够自动获取目标影像和参考影像的特征序列进而获得所述目标影像的平移距离，并根据所述平移距离对所述目标影像进行平移从而实现所述目标影像和所述参考影像的对齐。本发明所述医学影像对齐方法基本无需人工参与，提升了影像对齐处理的效率。
[0065]
请参阅图2，于本发明的一实施例中，所述医学影像对齐方法包括：
[0066]
s21，获取所述目标影像中参考对象的特征值，并根据所述目标影像中参考对象的特征值获得所述目标影像的特征序列；其中，所述参考对象为所述目标影像与所述参考影像中的骨骼或器官；所述目标影像包括至少2张ct影像组成的ct序列，例如：图3a显示为本
实施例中所述目标影像包含的ct序列，该ct序列由多张ct影像组成。
[0067]
s22，获取所述参考影像中参考对象的特征值，并根据所述参考影像中参考对象的特征值获得所述参考影像的特征序列；其中，所述参考影像包括至少2张ct影像组成的ct序列，例如，图3b显示为本实施例中所述参考影像包含的ct序列，该ct序列由多张ct影像组成。
[0068]
s23，根据所述目标影像的特征序列和所述参考影像的特征序列获取平移距离；
[0069]
s24，根据所述平移距离对所述目标影像进行平移，以使所述目标影像和所述参考影像相对齐。
[0070]
需要说明的是，图3a和图3b仅示例性的说明所述目标影像和参考影像。实际应用中，通过ct扫描机将获得多期ct影像，选取其中一期ct影像作为所述参考影像，选取所述参考影像之外的某一期ct影像作为所述目标影像，通过本实施例所述医学影像对齐方法即可实现所述目标影像和所述参考影像的对齐。本实施例所述医学影像对齐方法基本无需人工参与，有利于提升对所述医学影像进行对齐处理的效率。
[0071]
于本发明的一实施例中，所述平移距离n为正整数、零或负整数。其中：
[0072]
当n为正整数时，表示所述目标影像的第1张ct影像与所述参考影像的第n 1张ct影像相对齐。若所述目标影像的拍摄顺序是自下往上，所述目标影像的第1张ct影像为最下端的影像，则需要将所述目标影像中的所有ct影像向上平移n张，从而使得所述目标影像的第1张ct影像与所述参考影像的第n 1张ct影像对齐，所述目标影像的第2张ct影像与所述参考影像的第n 2张ct影像对齐，以此类推。若所述目标影像的拍摄顺序是自上往下，所述目标影像的第1张ct影像为最上端的影像，则需要将所述目标影像中的所有ct影像向下平移n张，从而使得所述目标影像的第1张ct影像与所述参考影像的第n 1张ct影像对齐，所述目标影像的第2张ct影像与所述参考影像的第n 2张ct影像对齐，以此类推。例如，若图3a中的第1张ct影像与图3b中的第3张ct影像相对齐，则将图3a所示ct序列整体上移2张从而实现图3a和图3b的对齐，即：图3a的第1张ct影像与图3b的第3张ct影像相对齐，图3a的第2张ct影像与图3b的第4张ct影像相对齐。
[0073]
当n为零时，所述目标影像和所述参考影像已经对齐，此时无需对所述目标影像进行平移处理。
[0074]
当n为负整数时，表示所述目标影像的第-n 1张ct影像与所述参考影像的第1张ct影像相对齐。若所述目标影像的拍摄顺序是自下往上，所述目标影像的第1张ct影像为最下端的影像，则需要将所述目标影像中的所有ct影像向下平移-n张，以使所述目标影像的第-n 1张ct影像与所述参考影像的第1张ct影像相对齐，所述目标影像的第-n 2张ct影像与所述参考影像的第2张ct影像相对齐，以此类推。若所述目标影像的拍摄顺序是自上往下，所述目标影像的第1张ct影像为最上端的影像，则需要将所述目标影像中的所有ct影像向上平移-n张，以使所述目标影像的第-n 1张ct影像与所述参考影像的第1张ct影像相对齐，所述目标影像的第-n 2张ct影像与所述参考影像的第2张ct影像相对齐，以此类推。
[0075]
请参阅图4a，于本发明的一实施例中，所述目标影像包括至少2张ct影像；所述目标影像中各ct影像均与所述参考对象的一个或多个特征值相对应；所述目标影像中参考对象的特征值包括所述目标影像中各ct影像对应的参考对象的特征值。此外，于本实施例中，所述参考对象与周围相邻器官或骨骼的ct值差异明显，因而能够通过ct值将所述参考对象
与周围相邻器官或骨骼进行区分。于本实施例中，获取所述目标影像中任一ct影像对应的参考对象的特征值的实现方法包括：
[0076]
s211a，获取所述参考对象的ct值范围；ct影像的数值单位是亨氏单位(hounsfield unit，hu)，本质上反映了人体器官或组织的密度值。不同参考对象具有不同的ct值范围，例如：对于骨骼，其ct值范围大于1000hu。
[0077]
s212a，根据所述ct值范围获取所述参考对象的像素点，进而生成所述参考对象的蒙版；理想状况下，所述参考对象的蒙版能够全覆盖所述参考对象，并由所述参考对象的全部像素点组成。
[0078]
s213a，根据所述参考对象的蒙版获取所述参考对象的特征值。
[0079]
上述步骤s211a～s213a能够获得所述目标影像中任意ct影像对应的参考对象的特征值；进一步的，将所述目标影像中所有ct影像对应的参考对象的特征值合并在一起，即可获得所述目标影像中参考对象的特征值。
[0080]
请参阅图4b，于本发明的一实施例中，对于任一ct影像，根据所述ct值范围获取所述参考对象的像素点，进而生成所述参考对象的蒙版的一种实现方法包括：
[0081]
s41，在所述ct影像内，删除所有ct值位于所述ct值范围之外的像素点。由于本实施例中所述参考对象的ct值与周围相邻器官或骨骼的ct值相差较大，因而能够通过ct值将所述参考对象与周围相邻器官或骨骼进行区分，通过步骤s41能够删除与所述参考对象相邻的器官和/或骨骼，进而获得一个或多个位置上离散的器官或骨骼。例如，选取脊柱作为所述参考对象时，经过步骤s41后所述ct影像内仅包含脊柱、肋骨等位置上离散的骨骼。
[0082]
s42，根据所述参考对象的位置和形状生成一个选择框；所述选择框包含所述参考对象的部分像素点。所述选择框可以根据先验知识由人工指定，也可以通过采用人工智能技术进行选取，此处不作限制。由于所述参考对象为一封闭且连续的形状，因此，所述选择框的至少一条边框包含有所述参考对象的像素点。优选地，所述选择框为矩形框。
[0083]
s43，对所述选择框的边框进行调整，以使所述选择框的所有边框均不包含所述参考对象的像素点；此时，所述选择框即为脊柱的外边框。具体地，遍历所述选择框的所有边框，并选取包含有所述参考对象像素点的边框作为待调整边框；对所有待调整边框向指定方向调整，直到所有待调整边框均不包含所述参考对象的像素点。
[0084]
例如，对于矩形选择框，若该选择框的下边框包含所述参考对象的像素点，则执行一调整步骤，并判断调整后的下边框是否仍包含所述参考对象的像素点：若所述调整后的下边框不包含所述参考对象的像素点，则对所述下边框的调整结束；否则，重复所述调整步骤直到所述调整后的下边框不包含所述参考对象的像素点。其中，所述调整步骤为将所述下边框向下平移一段距离并相应延伸所述左边框和右边框。
[0085]
s44，所述外边框内所有像素点组成的集合即为所述参考对象的蒙版。
[0086]
于本实施例中，通过参考对象的位置和形状获取所述参考对象的蒙版，并以所述蒙版为依据提取所述参考对象的特征值。由于参考对象为器官或骨骼，其形状多为不规则形，直接获取所述参考对象的特征值难度较大，而通过将所述选择框设置为矩形等常规形状能够简化所述蒙版的获取过程，在此基础上对参考对象特征值的提取相对简单。
[0087]
于本发明的一实施例中，所述医学影像中的一张ct图像如图4c所示，其中，选取脊柱51作为所述参考对象。于本实施例中，步骤s41删除所有ct值小于1000hu的像素点，删除
后所述ct图像仅包含脊柱51和肋骨52。请参阅图4d，本实施例步骤s42中根据先验知识生成一矩形的选择框4，该选择框4具有上边框411、左边框412、右边框413和下边框414。步骤s43中对所述选择框的边框进行调整，直到所述选择框的所有边框均不包含所述脊柱51的像素点。例如：请参阅图4e，首先将右边框413不断向右平移直到所述右边框413上不再包含所述脊柱51的像素点，此时得到调整后的右边框413’，相应延长所述上边框411和所述下边框414使得所述选择框保持为矩形；其后将下边框414不断向下平移直到所述下边框414上不再包含所述脊柱51的像素点，此时得到调整后的下边框414’，相应延长所述左边框412和所述调整后的右边框413’使得所述选择框保持为矩形。接下来，请参阅图4f，将左边框412不断向左平移直到所述左边框412不再包含所述脊柱51的像素点，此时得到调整后的左边框412’；将上边框411向上平移，直到所述右边框411不再包含所述脊柱51的像素点，此时得到调整后的上边框411’，相应延长所述调整后的上边框411’、所述调整后的左边框412’以及所述调整后的右边框413’，使得调整后的选择框保持为矩形，所述调整后的选择框即为脊柱的外边框4’。此时，所述外边框4’内包含了所述脊柱51的所有像素点；在所述外边框4’内，所述脊柱51的所有像素点组成的集合即为所述脊柱51的蒙版。
[0088]
于本发明的一实施例中，对于任一ct影像，根据所述ct值范围获取所述参考对象的像素点，进而生成所述参考对象的蒙版的一种实现方法包括：在所述ct影像内，删除所有ct值位于所述ct值范围之外的像素点；通过该步骤能够删除与所述参考对象相邻的器官和/或骨骼，进而获得一个或多个位置上离散的器官或骨骼。从所述一个或多个位置上离散的器官或骨骼中选取面积最大的闭合区域作为所述参考对象的蒙版。
[0089]
于本发明的一些实施例中，获取所述参考影像中参考对象的特征值的实现方法与获取所述目标影像中参考对象的特征值的实现方法类似，此处不再赘述。
[0090]
于本发明的一实施例中，所述参考对象对应凸包的几何特征、所述参考对象的中心点和/或所述参考对象的面积。其中，所述参考对象对应凸包可以通过所述蒙版计算得到，所述凸包的几何特征可以为所述凸包的面积、边长和/或边数等；所述参考对象的中心点可以通过所述蒙版的外边框的中心点代替；此外，本实施例中可以通过统计所述蒙版中包含的像素点个数作为所述参考对象的面积。
[0091]
请参阅图5，于本发明的一实施例中，根据所述目标影像中参考对象的特征值获得所述目标影像的特征序列的实现方法包括：
[0092]
s211b，根据所述目标影像中参考对象的特征值获得所述目标影像对应的特征向量。其中，所述目标影像中的所有ct影像均对应一个特征向量；某ct影像对应的特征向量由该ct影像对应的参考对象的特征值组成；所述目标影像对应的参考向量包括所述目标影像中所有ct影像对应的特征向量。例如，对于所述目标影像中的第l张ct影像，其对应的特征向量f
l
＝[f
l,1
,f
l,2
...f
l,m
]，其中，m为任一张ct影像对应的参考对象的特征值数量，其取值为正整数；f
l,k
为该第l张ct影像对应的第k个特征值的数值；k为正整数且k≤m；l为正整数且l≤n；
[0093]
s212b，根据所述目标影像对应的特征向量获得所述目标影像的特征序列。所述目标影像的特征序列由所述目标影像中所有ct影像对应的特征向量组成，优选地，所述目标影像的特征序列的计算公式为：其中，n为所述目标影像包含的ct影像数量，其取值为正整数。
[0094]
所述参考影像的特征序列的计算方法与所述目标影像的特征序列的计算方法相似，此处不再赘述。
[0095]
请参阅图6，于本发明的一实施例中，根据所述目标影像的特征序列和所述参考影像的特征序列获取平移距离的实现方法包括：
[0096]
s231，获取所述平移距离的取值范围；所述取值范围包括至少两个整数值；所述整数值可以为正整数、负整数或零；所述取值范围例如为[-3,3]或[-4,4]。
[0097]
s232，根据所述目标影像的特征序列和所述参考影像的特征序列，依次计算所述取值范围内的各整数值对应的序列方差，即：遍历所述取值范围内的所有整数值，并对每一整数值计算各自对应的序列方差；
[0098]
s233，根据所述序列方差选取相应的整数值作为所述平移距离。优选地，选取对应的序列方差最小的整数值作为所述平移距离。
[0099]
本实施例中，通过计算所有整数值对应的序列方差并选取序列方差最小的整数值作为所述平移距离，能够保证平移后的目标影像与所述参考影像之间的差值最小，从而实现所述参考影像和所述目标影像良好的对齐效果。
[0100]
于本发明的一实施例中，根据所述目标影像的特征序列和所述参考影像的特征序列计算与所述取值范围内任一整数值x对应的序列方差的公式为：
[0101][0102]
其中，n为所述目标影像包含的ct影像数量，m为所述目标影像中各ct影像对应的参考对象的特征值数量。
[0103]
差值表示所述目标影像中第i张图像平移x矢量距离后，其第j个特征值与所述参考影像中对应图像的对应特征值的差值，且其中，表示中第i行第j列的特征值，表示中第i行第j列的特征值。当所述目标影像的特征序列正向平移时，即x≥0时，其第i张图像的对应图像为所述参考影像中的第i x张图像；当所述目标影像的特征序列负向平移时，即x<0时，其第i张图像的对应图像为所述参考影像中的第i-x张图像。
[0104]
对于所述差值其在整个特征序列中的期望值表示所述目标影像整体平移x矢量距离后，其第j个特征值与所述参考影像中对应图象的对应特征值的期望。所述期
望值的计算公式为
[0105]
上述的计算公式中选择对m个特征序列差值的方差做整体求和处理，使得在计算过程中能够参考多个特征，从而使得平移后的整体效果最优，以防单个特征的提取不准确造成的误差。
[0106]
于本发明的一实施例中，所述参考对象为脊柱。对于腹部ct影像来说，脊柱的位置和形状相对比较固定，因此选取脊柱作为所述参考对象能够提升所述医学影像的对齐效果。
[0107]
于本发明的一实施例中，所述目标影像包括：平扫期ct影像、静脉期ct影像和/或动脉期ct影像。
[0108]
于本发明的一实施例中，选取脊柱作为所述参考对象。请参阅图7a，所述医学影像对齐方法包括：
[0109]
s71，提取图像脊柱蒙版；
[0110]
s72，根据所述脊柱蒙版提取各ct影像中脊柱的特征值；
[0111]
s73，分别获得平扫期ct影像、静脉期ct影像和动脉期ct影像的特征序列；
[0112]
s74，以所述平扫期ct影像为基准，分别平移所述静脉期ct影像和所述动脉期ct影像以使三期ct影像对齐。
[0113]
本实施例中步骤s71采用阈值法提取脊柱蒙版。ct影像的数值单位是亨氏单位，其反映了人体器官或组织的密度值。通常情况下，人体骨骼的密度大于1000hu值。因此，采用阈值法提取脊柱蒙版时只需将1000hu以下的组织过滤掉即可获得骨骼。为了避免肋骨等非脊柱骨骼对提取结果的影响，步骤s71中进一步采用步骤s41～s44所述蒙版提取方法获取所述脊柱的蒙版。
[0114]
请参阅图7b，本实施例中，步骤s72中提取各ct影像中脊柱的特征值的实现方法包括：
[0115]
s721，计算各ct影像中脊柱的中心点。由于本实施例中所述脊柱蒙版的外边框为矩形，根据该外边框的对称性，可近似认为所述蒙版的中心点即为所述脊柱的中心点；
[0116]
s722，计算所述脊柱的面积。具体地，可以通过统计所述脊柱蒙版所包含的像素点的数量作为所述脊柱面积。
[0117]
s723，获取脊柱的几何特征。具体地，通过对提取的脊柱蒙版进行计算可以得到凸包，并依次计算所述凸包的面积、参数等作为脊柱的几何特征。
[0118]
s724，根据所述脊柱的中心点、面积和几何特征组合成特征向量。
[0119]
通过步骤s721～s724可以获得任一ct影像中脊柱的特征值，通过获取平扫期所有ct影像的特征值即可获得平扫期ct影像的特征值，进而可以获得所述平扫期ct影像的特征序列；通过获取静脉期所有ct影像的特征值即可获得静脉期ct影像的特征值，进而可以获得静脉期ct影像的特征序列；通过获取动脉期所有ct影像的特征值即可获得动脉期ct影像的特征值，进而可以获得动脉期ct影像的特征序列。
[0120]
请参阅图7c，于本实施例中，步骤s74的实现方法包括：
[0121]
s741，依次选取[-m，m]范围内的所有整数值，并根据各选取的整数值依次沿脊柱方向平移静脉期ct影像；其中，m为正整数，优选为3；
[0122]
s742，基于静脉期ct影像的特征序列和动脉期ct影像的特征序列，依次计算得出各选取的整数值对应的序列方差；
[0123]
s743，选取对应的序列方差最小的整数值作为平移距离，进而对所述静脉期ct影像进行平移，从而使所述静脉期ct影像与所述平扫期ct影像进行对齐。
[0124]
对所述动脉期ct影像进行平移以使所述动脉期ct影像与所述平扫期ct影像进行对齐的方式与s741～s743类似，此处不再赘述。
[0125]
基于以上对所述医学影像对齐方法的介绍，本发明还提供一种计算机可读存储介质，其上存储有计算机可读程序。所述计算机程序被处理器执行时实现本发明所述医学影像对齐方法。
[0126]
基于以上对所述医学影像对齐方法的介绍，本发明还提供一种电子设备。请参阅图8，所述电子设备800包括：存储器810，存储有一计算机程序；处理器820，与所述存储器通信相连，调用所述计算机程序时执行本发明所述医学影像对齐方法；显示器830，与所述处理器和所述存储器通信相连，用于显示所述医学影像对齐方法的相关gui交互界面。
[0127]
本发明所述的医学影像对齐方法的保护范围不限于本实施例列举的步骤执行顺序，凡是根据本发明的原理所做的现有技术的步骤增减、步骤替换所实现的方案都包括在本发明的保护范围内。
[0128]
本发明所述医学影像对齐方法能够自动获取目标影像和参考影像的特征序列，并依据所述特征序列获取目标影像的平移距离进而对所述目标影像进行平移以实现所述目标影像与所述参考影像相对齐。整个过程基本无需人工参与，操作方便且效率较高。
[0129]
本发明所述医学影像对齐方法能够根据参考对象的三维生长形状得到一个误差最小的平移距离，进而允许对ct影像进行对齐预处理，从而使得通过人工智能进行ct影像数据处理成为可能。
[0130]
综上所述，本发明有效克服了现有技术中的种种缺点而具高度产业利用价值。
[0131]
上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟悉此技术的人士皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因此，举凡所属技术领域中具有通常知识者在未脱离本发明所揭示的精神与技术思想下所完成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。