心肌聚类与取向的制作方法
【专利摘要】一种心脏成像方法，包括采集包括心肌的投影图像表示（SI00）。分割所述心肌并生成掩模（SI02）。优化所述掩模（SI04）。从经优化的掩模确定血池（SI06），并且基于心肌切片的聚类使所述掩模骨架化（S108）。从所述血池和经骨架化的掩模来确定质心（SI10）。从所述经骨架化的掩模来确定心肌参数（SI12）。
【专利说明】心肌聚类与取向
【技术领域】
[0001]本申请涉及医学成像技术。其尤其应用于诊断核成像中对心肌组织的感兴趣区域(ROI)的限定。
【背景技术】
[0002]在诊断核成像中，放射性核素分布因其穿过患者的血流而得到研究，用于对循环系统成像，或用于对积累注射的放射性药物的特定器官成像。在单光子发射计算机断层摄影(SPECT)中，例如，通常被称作伽马照相机的一个或多个辐射探测器用于经由通过放射性衰变事件造成的辐射发射来探测所述放射性药物。典型地，每个伽马照相机包括辐射探测器阵列和设置在所述辐射探测器阵列前方的准直器。所述准直器限定线性或小角度锥形瞄准线，从而探测到的辐射包括投影数据。如果在角度视野的范围上，例如在180°或360°的角度范围上移动所述伽马照相机，则能够使用滤波反投影、期望值最大化或另一图像重建技术，将得到的投影数据重建成所述患者中所述放射性药物分布的图像。有利地，所述放射性药物能够被设计为在选定的组织中进行积累，用于提供那些选定的组织的优先成像，诸如，心脏组织，以达到心脏成像的目的。
[0003]在许多心脏SPECT研究中，最广泛使用的诊断性应用之一包括心肌灌注成像，其中，左心室是尤其感兴趣的。左心室(LV)区在SPECT图像中的低强度与典型地因冠状动脉疾病造成的灌注缺陷有关。LV区中的放射性药物的活性能够用于估计以下参数，诸如，血流率、血流储备、射血分数，或与诊断和处置有关的其他参数。为了估计这些以及其他参数，人们想要在至少所述左心室的感兴趣区域处进行定位，并且确定所述感兴趣区域的几何结构和取向，以准确地估计指示患者的心肌健康的参数。因为对所述心肌参数的准确量化评估部分地取决于准确且可重复地确定所述心肌几何结构和定向，更具体为左心室姿态估计。
【发明内容】
[0004]本发明提供一种克服了上文提及的问题以及其他问题的新型改进的方法和系统。
[0005]根据一个方面，提供一种用于心脏成像的方法。采集对象的功能图像数据，所述功能图像数据至少包括心肌。确定所述心肌的姿态和几何结构。基于所确定的姿态和几何结构，估计所述心肌的诊断参数。
[0006]根据另一方面，所述方法还包括分割所述左心室、生成左心室心肌掩模、计算所述左心室的质心、使所述左心室心肌掩模骨架化以生成心肌骨架以及修剪所述心肌骨架。
[0007]根据另一方面，提供一种心脏成像装置，其包括诊断扫描器，用于从所述诊断扫描器生成对象的功能成像数据(所述功能成像数据包括所述心肌)、确定所述心肌的姿态和几何结构以及基于所确定的姿态和几何结构来估计所述心肌的诊断参数。显示设备至少显示所述诊断参数。
[0008]根据另一方面，提供一种心脏成像装置，其包括采集单元，所述采集单元采集对象的功能图像数据，所述功能图像数据至少包括心肌。所述确定单元确定所述心肌的姿态和几何结构。估计单元基于所确定的姿态和几何结构来估计所述心肌的诊断参数。
[0009]一个优点在于改善心肌姿态计算的准确性和鲁棒性。
[0010]另一优点在于基于所计算的姿态改善参数的准确性。
[0011]本领域技术人员在阅读和理解以下详细描述时，将认识到本发明仍具有的其他优点。
【专利附图】
【附图说明】
[0012]本发明可以采取各种部件和部件的布置，以及各种步骤和步骤的安排的形式。附图仅出于图示优选实施例的目的，并且不应被解释为限制本发明。
[0013]图1为具有处理单元的诊断成像系统的图解图示，所述处理单元被编程为执行基于心肌聚类的骨架化与姿态估计；
[0014]图2为表示用于确定心肌参数的方法的流程图；
[0015]图3为表示用于心肌遮蔽和掩模优化的方法的流程图；
[0016]图4为表示用于心肌掩模优化的方法的另一流程图；
[0017]图5为表示用于基于心肌聚类的骨架化与姿态计算的方法的流程图；
[0018]图6图示了心肌组织的各种解剖结构，更具体为左心室；
[0019]图7A至图7D图示了各种类型的心肌形状簇；
[0020]图8A至图8B、图9A至图9B、图1OA至图10B，以及图1lB至图1lD图示了针对类型I心肌簇的骨架化方法，同时图1lA为表示用于针对类型I心肌簇的骨架化方法的方法的流程图；
[0021]图12A和图12B为表示用于针对类型2心肌簇的骨架化方法的方法的流程图；
[0022]图13A至图13D、图14A至图14E，以及图15A至图15D图示了针对类型2心肌簇的骨架化方法；
[0023]图16为表示用于针对类型3心肌簇的骨架化方法的方法的流程图；
[0024]图17A至图17D以及图18A至图18B图示了针对类型3心肌簇的骨架化方法；
[0025]图19A至图19E图示了针对类型4簇的骨架化方法；
[0026]图20为表示用于心肌血池分割的方法的流程图；
[0027]图21A至图21F图示了血池分割方法；
[0028]图22A为表示用于血池分割的另一种方法的流程图，同时图22B图示了所述血池分割方法；
[0029]图23A为表示用于血池分割的另一种方法的流程图，同时图23B图示了所述血池分割方法；
[0030]图24图示了来自各种形状簇类型的二元掩模及相应的血池的屏幕快照；
[0031]图25为表示骨架修剪方法的流程图；以及
[0032]图26图示了用于骨架修剪方法的直方图。
【具体实施方式】
[0033]参考图1，诊断成像系统10采集检查区域14内的对象12的功能成像数据。尽管描述了 SPECT系统，应认识到还预期其他成像模态，诸如，正电子发射断层摄影(PET)等。诊断成像系统10包括患者支撑体16，其可进行选择性地平移，以便于将正被成像或检查的对象12定位在期望位置处，例如，由此感兴趣区域关于纵轴18居中。任选地，外机架20被可移动地安装在轨道22上，以实现对象12沿纵轴18相对于成像系统10的期望定位。
[0034]内机架24被可旋转地安装在外机架20上，用于步进或连续运动。旋转的内机架24限定对象接收检查区域14。一个或多个探测头26、28、30能够单独定位在可旋转内机架24上。随着可旋转内机架24的旋转，探测头26、28、30作为整体关于检查区域14和对象12旋转。探测头26、28、30为径向、周向及横向可调节的，以变化它们与对象12的距离以及在旋转机架24上的间隔，用于以关于中心轴的多种角度取向中的任何一种来定位所述探测头。
[0035]探测头26、28、30的每个均包括辐射探测器阵列，诸如，响应于来自放射性药物的入射辐射事件而发射闪光或光子的多个闪烁体之一。通过光探测器阵列来观察所述(一个或多个)闪烁体，所述光探测器接收闪光并将其转化成电信号。任选地，也预期直接将辐射转化为电脉冲的探测器阵列。提供适当的准直以限定投影数据，例如，设置在所述探测器阵列前方的辐射吸收蜂窝准直器。解析器电路解析每个接收的辐射事件的X、y坐标，以及入射辐射的能量。以常规方式处理并校正光探测器的相对输出，以生成输出信号，所述输出信号指示(i)在探测头上接收每个辐射事件处的位置坐标、(ii)每个事件的能量以及(iii)所述探测头的角位置。所述能量用于在诸如多重发射辐射源、杂散及二次发射辐射、散射辐射、透射辐射的各种类型的辐射之间进行区分，并且用于消除噪声。
[0036]在SPECT成像中，由在所述探测头上的每个坐标处接收的所述辐射数据限定投影图像表示。在SPECT成像中，准直器限定沿其接收辐射的射线。应认识到，尽管所图示的实施例是关于SPECT成像进行描述的，也预期其他核成像模态，诸如，正电子发射断层摄影(PET)成像系统。
[0037]在PET成像中，针对两个头上的符合辐射事件监测探测头输出。从所述头的位置和取向以及表面上接收所述符合辐射处的位置，计算所述符合事件探测点之间的射线。该射线限定所述辐射事件沿其发生的线。在PET和SPECT中，来自多重角取向的投影数据被存储在投影数据缓冲器36中，并且然后由重建处理器38重建成所述感兴趣区域的横向体积图像表示和投影图像表示，它们被分别存储在体积图像存储器40和投影图像存储器42中。所述投影图像表示能够是体积的二维(2D)轴向表示，在所述体积中，选择沿通过所述体积数据集投影的线的最高衰减体素。与投影图像表示相比，利用体积图像表示，以计算时间和较小特征可视化为代价，保留了 3D空间关系。两个图像表示均可以包括沿纵轴18的图像表示的多个横向切片。功能成像系统10由控制器44操作，以执行所述对象的选定的目标区域的选定的成像序列。控制台46包括显示单元48，其显示图形用户界面(⑶I)，临床医师能够连同用户输入设备50 —起使用所述图形用户界面，用于控制扫描器控制器44，以选择扫描序列或协议。
[0038]参考图1和图2，系统10包括处理单元60，其包括执行用于心肌分割62、心肌掩模优化64、心肌掩模骨架化及修剪66、血池确定68、质心(COM)确定70以及参数化72的算法的若干个处理器。处理单元60被耦合到计算机可读存储单元74，计算机可读存储单元74存储接收到的图像数据；经处理的图像数据；用于确定、处理、生成、重建等的算法；用于确定心肌分割、心肌掩模校正、COM确定、心肌掩模骨架化、骨架修剪等的算法。应认识到，投影数据缓冲器36、体积图像存储器40、投影图像存储器42以及存储介质74可以为单一计算机可读存储器模块的部分，或作为单独的模块来实现。
[0039]具有各种处理器62、64、66、68、70、72的处理单元60执行姿态估计算法78，其解释心肌形状、几何结构、心内膜、心外膜以及对其他器官的邻近度。骨架化算法基于横向心肌切片的聚类，利用厚度切片，并且然后将每个2D骨架合并到3D骨架中。质心计算使血池分割适合于确定针对每个切片的X、y_方向COM,并且使3D骨架适合于确定z-方向COM。
[0040]在利用系统10采集SlOO经重建的投影体积之后，心肌分割处理器62接收存储在体积图像存储器40中的经重建的体积表示(和/或投影图像表示)，并且分析所述体积的每个切片，以限定所述心肌的三维(3D)感兴趣区域(ROI)的位置和边界，而无需诸如在美国申请61/311,406中描述的解剖图像或人类操作者的帮助。一旦已确定了所述R0I，生成相应的ROI的3D 二元掩模S102。能够在投影和/或体积图像表示上重叠所述二元掩模，并且能够将其提供到显示单元48，以便由临床医师进行查看。
[0041]心肌掩模优化处理器64从分割处理器62接收所确定的心肌掩模，并针对缺陷或偏差来校正所述掩模S104。由过度灌注或灌注不足造成的所述心肌掩模中的偏差，特别是在横向体积中的偏差，能够导致由参数化处理器72确定的参数中的显著误差。
[0042]骨架化处理器66从优化处理器64接收经优化的心肌掩模，并且以基于聚类的骨架化算法为基础，确定针对每个经优化的掩模的心肌组织的骨架。所述基于聚类的骨架化算法将经优化的掩模归类成各种形状簇。然后执行适于每种簇类型的2D骨架化方法。将所述2D骨架合成为3D堆叠，然后通过距离分箱和修剪算法来优化所述3D堆叠，所述距离分箱和修剪算法使用从血池分割处理器68确定的COM。血池确定与分割处理器68从优化处理器64接收经优化的心肌掩模，并将心肌空腔的心内膜以内的像素确定指定为血池像素。
[0043]COM处理器70从血池处理器68接收经分割的血池掩模，并基于每个血池的形心(centroid)来确定所述COM的x和y坐标。基于所述心肌掩模在z方向上的范围，计算所述COM的z坐标。骨架化处理器66将所述血池的COM用于修剪算法，以修剪或去除能够影响后续椭球拟合的附加骨架像素。COM处理器68将椭球拟合到经优化的骨架，从所述经优化的骨架能够估计姿态，例如，方位角、倾斜角等。
[0044]几何结构与参数化处理器72从拟合的椭球确定心肌几何结构估计，诸如，心脏方位角、倾斜角等，并且基于所述R01、掩模、血池、骨架和/或椭球以及所确定的心肌几何结构估计来确定可量化参数。在心脏评估的背景下，可量化参数包括心肌血流量、区域性心肌血流量、血流储备、射血分数等。应认识到，也预期调整适用于特定成像研究，诸如，心脏成像、肺部成像、脑成像等的参数。
[0045]在心脏成像与参数化中，左心室是尤其感兴趣的。以下算法将参考左心室(L.V)成像进行描述。例如，处理单元60的处理器62、64、66、68、70能够分别执行L.V.分割、L.V.掩模优化、L.V.掩模骨架化、L.V.血池确定以及L.V质心确定。此外，参数化处理器72确定若干个参数，例如，L.V.的血流量、血流储备、射血分数等。应认识到,处理单元60能够对心脏区域的各种结构执行以下算法。
[0046]掩模优化处理器64针对两种极端情况优化所述心肌掩模:过度分割，其中，心肌组织为严重过度灌注的；以及分割不足，其中，所述心肌组织为灌注严重不足的。L.V.横向体积的顶端和底端对于椭球拟合和L.V.角度确定而言是重要的，这将在后文中进行讨论。L.V.角度的不准确能够显著影响由参数化处理器60确定的参数。
[0047]参考图3，提出用于校正受过度分割影响的二元掩模的算法80，以修剪所述L.V.横向体积的顶部和底部的不必要厚度，并且然后恰当地形成顶部和底部切片。在所述L.V.体积被分割并且生成相应的二元掩模S102之后，然后连同不具有空腔的第一切片和最末切片一起S202，定位具有封闭或半封闭的空腔的第一切片和最末切片S200，所述不具有空腔的第一切片和最末切片将毗邻来自步骤S200的所述切片。一旦确定这四个切片的位置，确定在所述L.V.掩模的顶部和底部的深度S204。在任一端的深度是指除来自S202的切片以外，不具有空腔的切片的堆叠。然后通过去除在S200和S202中确定的四个切片以外的不必要的切片，优化或净化所述L.V.空腔的切片S206。根据预选的阈值，确定经“净化的”切片的数目。当所述L.V.为过度灌注时(导致过度分割)，第一切片(不具有空腔)可以具有比下一切片(具有空腔)更多或附加的像素。最末切片和倒数第二切片会出现相同的情形。在这些情形中，优化所述第一切片和最末切片的心肌掩模，以改变所述掩模的形状，用于分别密切匹配所述第二切片和倒数第二切片的形状。在选择顶部两个切片和最末两个切片S208之后，然后利用所述第二个(具有空腔的切片)遮蔽所述第一切片(不具有空腔)，并且类似地，利用倒数第二切片(具有空腔)遮蔽所述最末切片(不具有空腔)S210。在两个集合上执行2D连通分量分析S212，以在所述L.V.体积的顶部和底部定位间断分量(brokencomponents)ο如果定位了间断分量，则对第一掩模和最末掩模的优化是不正确的。在所述第一切片和最末切片上执行形态学操作S214，诸如，修剪，使得所述第一切片密切匹配所述第二切片的形状，并且所述最末切片密切匹配倒数第二切片的形状。
[0048]参考图4，提出用于校正受分割不足影响的二元掩模的算法82，以填充因灌注不足而缺失所述L.V.的心肌的部分处的间隙。定位所述L.V.横向体积的顶部切片和底部切片S300，并且在每个切片上向所定位的间隙执行二元射线投射S302。如果间隙位于选定的切片中，则执行切片填充算法S304，以填充所定位的心肌间隙。针对所述L.V.体积的接下来的两个中间切片重复步骤302和S304，直到到达公共点。例如，如果在L.V.横向体积中存在25个切片，则将步骤302和S304重复13次。
[0049]在完成过度分割算法80和分割不足算法82之后，结果是在顶部和底部的平衡的心肌掩模，同时解决了诸如针对灌注不足的情况填充心肌间隙的问题。
[0050]参考图5和图6,骨架化处理器66利用基于簇的骨架化算法84,针对L.V.横向体积的每个经优化的二元掩模，生成骨架。基于对所述骨架的椭球拟合，确定心脏角，诸如，方位角和倾斜角。骨架化处理及其准确度对于角度生成是关键的。骨架化算法84基于2D中的所述心肌掩模切片的聚类。在生成所述骨架中，算法84利用了心肌86以及邻接器官的几何结构、心外膜88和心内膜90。基于每个切片中的水平射线(HRj)的数目和垂直射线(VRi)的数目S400，为所述切片分配合适的簇类型S402，并且之后相应地将每个簇类型进行归类，然后进行骨架化S404、S406、S408、S410。堆叠所述2D骨架S412，以形成3D骨架，然后修剪所述3D骨架S414，以去除边界处的任意附加像素。然后将椭球拟合S416到经修剪的骨架，从其确定L.V.角度(方位和倾斜)以及在z方向上的COM S418。
[0051]基于在它们各自方向上的射线命中(ray hits)，限定所述VRi和HR」。根据方程1，如果像素(i，j )为前景像素(F)，并且其连接像素之一为背景像素(B)，则像素(i，j )为射线命中RHi, j。[0052]
[典，U κ^Λ)€β]
履⑴
I假，其他
[0053]方程2和方程3将VRi和HRj分别限定为在一行中与至少4个非零像素交叉的垂
直射线和在一列中与至少4个非零像素交叉的水平射线。
[0054]
窗 A <r, V* /?//
jrn一、，— / -> /.1/ \
VRi = i;-0.\Z )
[fi,其他
[0055]
?假，ΛΛ
[0056]参考图7Α至图7D，限定所述心肌形状簇并且将其归类为类型1(图7Α)、类型2(图7Β)、类型3 (图7C)和类型4 (图7D)，从而将所述骨架化过程调谐到心肌掩模形状的每个引入的类型。以此方式，连同形状信息一起的切片的聚类确保了准确地处理具有或不具有形状缺陷的情况。
[0057]参考图8Α和图8Β，类型I簇被限定为具有VRi和HRj两者，即，VRiX)并且HRjXL为了确定针对类型I簇的骨架S404，执行垂直骨架化和水平骨架化。对于垂直骨架化，投射垂直射线VRi,并且根据方程4记录射线命中RHy:
[0058]
/, I V:
Ρν^ = {}^ΨΕ￥Η^(4)
[0059]其中，PVRi为射线VRi上的射线命中的集合。从所述射线命中PVRi,根据方程5确定骨架像素SVRi:
[0060]
'JPVRlllt^PVRnin, PVR1.,,+PVRljin I…
SVR = N——^——^-^^eF}(5)
* Hl 2 1 2I
[0061]并且SVRi表示所述骨架沿射线VRi的像素位置。每个骨架像素与垂直射线上的两个相应的命中之间是等距的。一旦获得所述骨架像素，则根据方程6，将起始垂直射线SV1和终止垂直射线EV1限定为第一和最末VRi:
[0062]
mm {i|^},m?{/|^}}(6)
[0063]接下来，确定SV1和EV1与所述骨架像素SVRi的交叉，并将其限定为SV水平界限起始(Horizontal Limit Start)svHSnEV水平界限起始evHSpSV水平界限终止(HorizontalLimit End) SvHE1, EV水平界限终止evHE^ SvHS1与之间的骨架像素限定所述骨架的一个水平部分，并且SvHE1与^HE1之间的骨架像素限定所述骨架的第二水平部分。
[0064]参考图9A和图9B，利用水平骨架化来确定类型I簇的骨架的垂直部分。对于水平骨架化，投射水平射线HRp并根据方程7记录射线命中RHiij:
【权利要求】
1.一种用于心脏成像的方法，包括: 采集对象的功能成像数据，所述功能成像数据至少包括心肌； 确定所述心肌的姿态和几何结构； 基于所确定的姿态和几何结构，估计所述心肌的诊断参数。
2.根据权利要求1所述的方法，其中，确定所述心肌的姿态和几何结构的所述步骤包括: 分割左心室； 生成左心室心肌掩模； 计算所述左心室的质心； 使所述左心室心肌掩模骨架化，以生成心肌骨架； 修剪所述心肌骨架。
3.根据权利要求2所述的方法，其中，修剪所述心肌骨架包括: 计算从每个非零骨架像素到所述左心室的所述质心的距离； 将所计算的距离分箱到距离分箱中； 从每个所述距离分箱中的多个 像素来生成直方图； 从所述直方图确定下部分箱阈值、上部分箱阈值、下部分箱频率阈值，以及上部分箱频率阈值； 去除在所述下部分箱阈值与所述上部分箱阈值以及所述下部分箱频率阈值与所述上部分箱频率阈值之外的骨架像素。
4.根据权利要求2和3中的任一项所述的方法，还包括: 针对边界增厚和心肌间隙，校正所述左心室掩模。
5.根据权利要求2至4中的任一项所述的方法，还包括: 分割所述左心室掩模，以生成左心室掩模切片的堆叠； 定位所述堆叠的顶部和底部； 定位具有封闭或半封闭空腔的所述堆叠中的最上部掩模切片和所述堆叠中的最下部掩模切片； 利用所述堆叠中的第二最上部掩模切片来遮蔽所述堆叠中的所述最上部掩模切片； 利用所述堆叠中的第二最下部掩模切片来遮蔽所述堆叠中的所述最下部掩模切片；以及 在最上部切片与最下部切片之间进行填充。
6.根据权利要求2至5中的任一项所述的方法，针对严重灌注情况还包括: 将所述左心室掩模分割成掩模切片； 使用射线投射来定位具有心肌间隙的最上部掩模切片和最下部掩模切片，以定位空腔； 应用切片填充算法，以在所述最上部掩模切片与所述最下部掩模切片之间进行填充。
7.根据权利要求2至6中的任一项所述的方法，其中，骨架化包括: 将所述左心室掩模分割成掩模切片； 使用垂直射线投射和水平射线投射来对所述掩模切片进行聚类，以形成针对每个切片的多个簇；从所述簇限定血池。
8.根据权利要求7所述的方法，其中，聚类所述掩模切片包括以下的一个或多个: 对于每个切片，投射垂直射线以限定所述心肌的射线命中和水平界限，并且投射水平射线以限定所述心肌空腔的射线命中和垂直界限，以及填充所述心肌中超出所述水平界限和垂直界限的部分； 遮蔽每个掩模切片，以去除在所述垂直界限之间的部分；基于紧邻所述垂直界限的像素的强度，对上部垂直界限以上和下部垂直界限以下的像素的强度进行分类；将经分类的所述上部垂直界限以上和所述下部垂直界限以下的像素与在所述垂直界限之间的像素合并，以限定净化的骨架； 遮蔽每个掩模切片，以去除在所述水平界限之间的部分；基于紧邻所述水平界限的像素的强度，对左侧水平界限以上和右侧水平界限以下的像素的强度进行分类；将经分类的所述左侧水平界限以上和所述右侧水平界限以下的像素与在所述水平界限之间的像素合并，以限定净化的骨架。
9.根据权利要求2至8中的任一项所述的方法，还包括: 识别所述心肌掩模中具有至少一个水平射线或垂直射线的切片； 识别接近主动脉瓣的每个水平射线的边缘和点； 识别接近二尖瓣的每个垂直射线的边缘和点，所述水平射线和所述垂直射线的所述边缘端点限定一条线； 将距所述线小于例如7mm的预选距离的像素标记为血池像素。
10.根据权利要求2至9中`的任一项所述的方法，还包括: 识别所述心肌掩模中具有至少一个水平射线或垂直射线的切片； 识别所述心肌掩模上与至少一个所述射线交叉的点； 填充所识别的点之间的像素作为血池像素。
11.被编程为执行根据权利要求1至10中的任一项所述的方法的一个或多个计算机处理器。
12.—种非暂态计算机可读介质，其承载控制一个或多个处理器执行根据权利要求1至10中的任一项所述的方法的软件。
13.—种心脏成像系统,包括: 诊断扫描器，其生成对象的功能成像数据； 一个或多个处理器，其执行根据权利要求2至10中任一项所述的方法；以及 显示器，在其上显示所述心肌的所述诊断参数。
14.一种心脏成像装置，包括: 功能成像数据采集单元，其采集对象的功能图像数据，所述功能成像数据至少包括心肌； 姿态和几何结构确定单元，其确定所述心肌的姿态和几何结构；以及 诊断参数估计单元，其基于所确定的姿态和几何结构，估计所述心肌的诊断参数。
15.根据权利要求14所述的装置，其中，所述姿态和几何结构确定单元包括: 分割单元，其分割所述对象的左心室； 掩模生成单元，其生成左心室心肌掩模；质心单元，其确定所述左心室的质心； 骨架化单元，其使所述左心室心肌掩模骨架化，以生成心肌骨架；以及 修剪单元，其修剪所述心肌骨架。
16.根据权利要求1和2中的任一项所述的装置，还包括以下中的一个或多个:分割单元(62)、掩模优化单元(64)、骨架化单元(66)、血池处理器(68)、质心处理器(70)以及心肌运动图几何结构与参数化处理器(72 )。
17.根据权利要求15所述的装置，其中，所述修剪单元包括: 距离确定模块，其确定从每个非零骨架像素到所述左心室的所述质心的距离； 分箱单元，其将所计算的距离分箱到距离分箱中； 直方图单元，其从每个所述距离分箱中的多个像素生成直方图； 阈值确定单元，其从所述直方图确定下部分箱阈值、上部分箱阈值、下部分箱频率阈值以及上部分箱频率阈值； 去除单元，其去除在所述下部分箱阈值与上部分箱阈值以及所述下部分箱频率阈值与上部分箱频率阈值之外的骨架像素。
18.根据权利要求15至17的中任一项所述的装置，还包括: 分割单元，其分割所述左心室掩模，以生成掩模切片堆叠； 定位单元，其定位所述切片堆叠的顶部和底部、具有封闭或半封闭空腔的所述堆叠中的最上部掩模切片和所述堆叠中的最下部掩模切片；` 遮蔽单元，其利用第二最上部掩模切片来遮蔽所述堆叠中的最上部切片，并且利用所述堆叠中的第二最下部掩模切片来遮蔽所述堆叠的最下部切片；以及填充单元，其在所述最上部切片与所述最下部切片之间进行填充。
19.根据权利要求15至18中的任一项所述的装置，其中，所述骨架化单元还包括: 分割单元，其将所述左心室掩模分割成掩模切片； 聚类单元，其使用垂直射线和水平射线来对所述掩模切片进行聚类，以针对每个切片投射多个簇；以及 血池单元，其从所述簇限定血池。
20.根据权利要求19所述的装置，其中，所述聚类单元包括以下中的一个或多个: 投射单元，其针对每个切片，投射垂直射线，以限定所述心肌空腔的射线命中和水平界限，并且投射水平射线，以限定所述心肌空腔的射线命中和垂直界限； 填充单元，其填充所述心肌超出所述水平界限和所述垂直界限的部分； 遮蔽单元，其遮蔽每个掩模切片，以去除所述垂直界限之间的部分；基于紧邻所述垂直界限的像素的强度，对上部垂直界限以上和下部垂直界限以下的像素的强度进行分类；以及将经分类的所述上部垂直界限以上和所述下部垂直界限以下的像素与在所述垂直界限之间的像素合并，以限定净化的骨架； 遮蔽单元，其遮蔽每个掩模切片，以去除在所述水平界限之间的部分；基于紧邻所述水平界限的所述像素的强度，对左侧垂直界限以左和右侧垂直界限以右的像素的强度进行分类，将经分类的所述水平垂直界限以左和右的像素与在所述水平界限之间的像素合并，以限定净化的骨架。
【文档编号】G06T7/00GK103518225SQ201280022233
【公开日】2014年1月15日 申请日期:2012年5月2日 优先权日:2011年5月10日
【发明者】S·德维韦迪 申请人:皇家飞利浦有限公司