CN103565449A - 一种脊柱x线图像的处理方法和系统 - Google Patents
一种脊柱x线图像的处理方法和系统 Download PDFInfo
- Publication number
- CN103565449A CN103565449A CN201210271466.6A CN201210271466A CN103565449A CN 103565449 A CN103565449 A CN 103565449A CN 201210271466 A CN201210271466 A CN 201210271466A CN 103565449 A CN103565449 A CN 103565449A
- Authority
- CN
- China
- Prior art keywords
- vertebral body
- spinal column
- ray image
- image
- coordinate information
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了脊柱X线图像的处理方法和系统。该方法包括:利用第二脊柱X线图像对预定的第一脊柱X线图像进行配准得到所述第一脊柱X线图像的空间变换信息;根据所述空间变换信息和所述第一脊柱X线图像中位于至少一对椎体边缘的预定点的第一坐标信息获得配准后所述预定点的第二坐标信息;利用所述第二坐标信息计算脊柱形变量并输出。本发明提供的方法和装置利用预先设定的脊柱X线图像及其中若干点的坐标信息可以自动计算脊柱形变量,能够减少或者避免人工操作对计算结果的影响,使计算结果更准确。
Description
技术领域
本发明涉及医学图像处理技术,特别涉及一种脊柱X线图像的处理方法和系统。
背景技术
脊柱形变,如后凸、侧弯等,早期发现可以防止形变的进一步发展。从脊柱的X线图像中可以看到脊柱的形变,并且通过对脊柱X线图像的定量测量可以确定脊柱的形变是否超过正常范围。
目前已出现一些手动和半自动的技术用于定量测量X射线图像中的脊柱形变,用户需要手动标记出X线图像中脊柱曲线两端最倾斜的椎体的边缘,甚至需要手动确定椎体边缘曲线的切线来进行形变量的计算。这些手动和半自动的方法需要用户的大量手动操作,导致测量结果的准确度和可靠性降低,针对同一图像,不同测量者的测量结果、同一测量者的两次测量结果之间均会存在差异。
发明内容
鉴于现有技术中存在的脊柱形变测量结果受手动操作的影响不够准确和可靠的问题,本发明提供了一种脊柱X线图像的处理方法和系统,能够减少或者避免人工操作对计算结果的影响,使计算结果更准确。
本发明实施例提供的一种脊柱X线图像的处理方法包括:
利用第二脊柱X线图像对预定的第一脊柱X线图像进行配准得到所述第一脊柱X线图像的空间变换信息;
根据所述空间变换信息和所述第一脊柱X线图像中位于至少一对椎体边缘的预定点的第一坐标信息获得配准后所述预定点的第二坐标信息;
利用所述第二坐标信息计算脊柱形变量。
本发明实施例提供的一种脊柱X线图像的处理系统包括:
图像配准模块,用于利用第二脊柱X线图像对预定的第一脊柱X线图像进行配准得到所述第一脊柱X线图像的空间变换信息;
坐标变换模块,用于根据所述图像配准模块获得的所述空间变换信息和所述第一脊柱X线图像中位于至少一对椎体边缘的预定点的第一坐标信息获得配准后所述预定点的第二坐标信息;
形变计算模块,用于利用所述第二坐标信息计算脊柱形变量。
本发明实施例提供的方法和系统利用预先设定的脊柱X线图像及其中若干点的坐标信息可以自动计算脊柱形变量,能够减少或者避免人工操作对计算结果的影响,使计算结果更准确。
附图说明
图1为本发明实施例提供的一种脊柱X线图像的处理方法的流程图。
图2为本发明实施例的一种脊柱X线图像的处理系统的结构示意图。
图3为本发明实施例的一种脊柱X线图像的处理系统的组成示意图。
图4为Cobb角的计算原理示意图。
图5为本发明实施例提供的一种脊柱X线图像的处理方法的流程图。
具体实施方式
为使本发明的实施例的目的、技术方案及优点更加清楚明白,以下参照附图并举实施例,对本发明实施例进一步详细说明。
图1为本发明实施例提供的一种脊柱X线图像的处理方法的流程图。如图1所示,该方法主要包括:
步骤101,利用第二脊柱X线图像对预定的第一脊柱X线图像进行配准得到所述第一脊柱X线图像的空间变换信息;
步骤102,根据所述空间变换信息和、所述第一脊柱X线图像中位于至少一对椎体边缘的预定点的第一坐标信息获得配准后所述预定点的第二坐标信息;
步骤103,利用所述第二坐标信息计算脊柱形变量并输出。
其中,第二脊柱X线图像可以是被测量者的脊柱X线图像。
第一脊柱X线图像是作为模板的标准脊柱X线图像,可以预先存储在图像处理系统中,也可以是该设备通过U盘或网络等其它途径获取得到的。
配准是指在医学诊断和治疗过程中,需要对比分析多幅图像时,对于一幅医学图像寻求一种或一系列空间变换,使它与另一幅医学图像上的对应点达到空间上的一致。这种一致是指人体上的同一解剖点在两张匹配图像上有相同的空间位置。配准的结果应使两幅图像上所有的解剖点,或至少是所有具有诊断意义的点及手术感兴趣的点都达到匹配。
医学图像配准可分为刚性配准和非刚性配准两类。刚性配准是指将图像作为刚体,仅进行线性的变换;而非刚性配准则将图像作为非刚体,进行非线性的变换,例如伸缩、旋转、扭曲等。由于脊柱的形变通常包括椎体的错位以及相对其它椎体的弯曲或扭曲等,例如脊柱侧弯、后凸时表现出的脊柱曲线的变化,本发明实施例优选使用非刚性配准。
步骤101中的配准过程可以使用现有的各种配准算法,如灰度驱动配准算法、模型驱动配准算法及混合驱动配准算法等,也可以使用专门设计的配准算法,本发明的实施例对此不加以限定。
通过步骤101的配准,对第一脊柱X线图像中的各个点进行位置变换,使第一脊柱X线图像与第二脊柱X线图像中的同一解剖点的空间位置一致。根据本发明的一个实施例,通过配准获得的所述空间变换信息可以包括配准后第一脊柱X线图像中若干个点的位移矢量,在步骤102中,利用该位移矢量和这若干个点在配准前的第一坐标信息,就可以对所述若干个点的第一坐标进行坐标变换得到配准后的第二坐标信息。
所述若干个点是预先在第一脊柱X线图像中选定的,用于在后面的步骤中计算脊柱的形变量。根据本发明的一个实施例,所述若干个点可以分别位于第一脊柱X线图像中两个椎体的边缘,这两个椎体以下称为第一椎体和第二椎体。这些点可以事先由具有丰富经验的医师或专业人员或利用专门的设备在第一脊柱X线图像中标记出来,使其尽量精确地体现椎体边缘的位置。这些点的坐标信息可以根据标记从第一脊柱X线图像数据中获得,也可以单独存储。
由于待测量者的脊柱X线图像是在进行检查时获得的,如果没有条件对其中的椎骨边缘进行精确的标记,例如由于时间的限制和医师经验的限制,直接从其中获得准确的形变测量结果是困难的。利用图像识别技术进行椎体边缘自动识别则受到图像识别算法自身性能的限制,准确性和可靠性都不能保证。而在上面的实施例中,通过利用事先进行了精确标记的第一脊柱X线图像,在配准后只需要对预先标记的点进行坐标变换即可自动得出形变测量结果,不需要医师在检查过程中进行手动操作,只要事先获得的模板图像足够准确,就能获得相对更为准确可靠的测量结果。
根据一个实施例,步骤103中可以根据预设的算法计算脊柱形变量。例如,可以根据所述若干个点的坐标信息确定所述第一椎体的第一边缘的切线和第二椎体的第二边缘的切线,计算所述切线之间的夹角;或者利用所述第二坐标信息计算所述第一椎体和第二椎体之间的Cobb角,等。
所述第一椎体和第二椎体可以是预先确定的用于进行脊柱形变测量的椎体。例如,可以由经验丰富的医师根据经验选择具有典型意义的一对或多对椎体作为测量使用的椎体。根据一个实施例,所述第一椎体和第二椎体可以是以下组合中的一个或多个:脊柱的顶部椎体和底部椎体,颈椎中的顶部椎体和底部椎体,胸椎中的顶部椎体和底部椎体,腰椎中的顶部椎体和底部椎体。当预设了多对需要测量的椎体时,可以针对这多对椎体重复执行上述步骤102、103来获得每对椎体的测量结果。
所述第一椎体和第二椎体还可以是用户选择的椎体,例如用户可以通过输入装置(如鼠标、键盘、触控装置等)来选择需要测量的椎体。根据一个实施例,在计算得出预先确定的第一椎体和第二椎体的Cobb角后,医师需要了解位于第一椎体和第二椎体之间的第三椎体到第四椎体这一段脊柱户的鼠标操作区域确定用户选择了第三椎体和第四椎体,然后将第三椎体和第四椎体作为第一椎体和第二椎体,执行上述过程以计算第三椎体和第四椎体之间的Cobb角。
由于脊柱的情况与性别、年龄、人种等有关,同时,脊柱的形变还包括侧弯、后凸等不同情况,不同情况采集的X线图也有所不同,因此,根据一个实施例,该方法还可以包括:
从预先存储的至少两个脊柱X线图像中选择一个作为所述第一脊柱X线图像。
预先存储的至少两个脊柱X线图像可以分别对应不同年龄、不同性别、不同的人种、不同的图像采集角度等等,这样,就可以根据年龄、性别、人种、图像采集角度中的一个或多个从所述至少两个脊柱X线图像中选择一个作为所述第一脊柱X线图像,从而得到更加科学和准确的测量结果。
根据一个实施例,可以预先对第一脊柱X线图像中的脊柱图像按照椎体进行分段并标记每个椎体的信息,如编号、名称和边缘等,输出脊柱形变量的时候,还可以同时输出所述第一椎体和第二椎体的信息,这样,可以让用户更直观地了解测量结果。例如,可以输出“椎体L2与L5之间角度为10度”等。
本发明实施例还提供一种脊柱X线图像的处理系统。图2为本发明实施例的一种脊柱X线图像的处理系统的结构示意图。如图2所示,该系统可包括:图像配准模块201、坐标变换模块202、形变计算模块203、和输出模块204。
图像配准模块201主要用于利用第二脊柱X线图像对预定的第一脊柱X线图像进行配准得到所述第一脊柱X线图像的空间变换信息,将所述空间变换信息提供给坐标变换模块202。
坐标变换模块202主要用于根据所述图像配准模块201获得的所述空间变换信息和所述第一脊柱X线图像中位于至少一对椎体边缘的预定点的第一坐标信息获得配准后所述若干个点的第二坐标信息。
形变计算模块203主要用于利用所述坐标变换模块202提供的第二坐标信息计算脊柱形变量,并将得到的所述脊柱形变量提供给输出模块204用于输出。
其中,第一脊柱X线图像可以预先存储在图像处理系统的存储模块中,该存储模块可以是一个单独的存储装置,也可以是分布在多个设备中的多个存储装置,还可以是内嵌在各模块中的存储单元的统称。可选的,第一脊柱X线图像也可以通过U盘或网络等其它途径获取。
根据一个实施例,图像配准模块201可以获得配准后所述若干个点的位移矢量作为所述空间变换信息;
坐标变换模块202可以利用所述图像配准模块201获得的所述位移矢量和所述第一坐标信息对所述若干个点的第一坐标进行坐标变换,得到所述第二坐标信息。
根据一个实施例,形变计算模块203可以根据所述若干个点的坐标信息确定第一脊柱X线图像中第一椎体的第一边缘的切线和第二椎体的第二边缘的切线,计算所述切线之间的夹角作为所述脊柱形变量。
根据一个实施例,所述形变计算模块203用于计算第一脊柱X线图像中第一椎体和第二椎体之间的Cobb角作为所述脊柱形变量。
根据一个实施例,坐标变换模块202可以根据预设的策略将脊柱的顶部椎体和底部椎体、或颈椎中的顶部椎体和底部椎体、或胸椎中的顶部椎体和底部椎体、或腰椎中的顶部椎体和底部椎体、或用户选择的两个椎体确定为所述第一椎体和第二椎体,从所述存储模块获取所述第一脊柱X线图像中的所述第一椎体和第二椎体边缘上的若干个点的坐标信息作为所述第一坐标信息。
下面再举一个例子。图3为本发明实施例的一种脊柱X线图像的处理系统的结构示意图。如图3所示,该系统可以包括:图像配准模块301、坐标变换模块302、形变计算模块303、输出模块304、图像选择模块305、目标图像数据库306和模板图像数据库307、输入模块308。
目标图像数据库306中存储有被测量者的脊柱X射线图像。
模板图像数据库307中预先存储有多个模板脊柱X射线图像。这些模板脊柱X线图像分别对应不同的年龄、性别、人种、图像采集角度等中的一种或多种。在这些模板图像中,所有的椎体已经由专家仔细分段并标记出其边缘。
图像选择模块305用于从目标图像数据库306选择目标脊柱X线图像,从模板图像数据库307选择模板脊柱X线图像,将所述目标脊柱X线图像的信息和所述模板脊柱X线图像的信息提供给所述图像配准模块301。
所述图像选择模块305可以根据被测量者的年龄、性别、人种、所述目标脊柱X线图像对应的图像采集角度等信息中的一个或多个来选择所述模板脊柱X线图像。
图像配准模块301主要用于利用所述目标图像数据库306中的所述目标脊柱X线图像对模板图像数据库307中预先存储的所述模板脊柱X线图像进行非刚性图像配准,将模板脊柱X线图像中已分段并标记的椎体进行扭曲使之与目标脊柱X线图像中的各椎体相对应,得到所述模板脊柱X线图像的空间变换信息,将所述空间变换信息提供给坐标变换模块302。根据一个实施例,图像配准模块301还可以在完成图像配准后根据模板脊柱X线图像中椎体的分段信息和标记信息对所述目标脊柱X线图像中的脊柱进行分段和标记。
坐标变换模块302主要用于确定需要测量的第一椎体和第二椎体,从模板图像数据库307中获取预先存储的、模板脊柱X线图像中所述第一椎体和所述第二椎体边缘的若干个点的第一坐标信息,根据所述图像配准模块301获得的所述空间变换信息对所述第一坐标进行坐标变换,获得配准后所述若干个点的第二坐标信息,将第二坐标信息提供给形变计算模块303。
坐标变换模块302可以根据预设的策略确定需要测量的第一椎体和第二椎体,例如,预设的第一椎体和第二椎体可以为脊柱的顶部椎体和底部椎体、或颈椎中的顶部椎体和底部椎体、或胸椎中的顶部椎体和底部椎体、或腰椎中的顶部椎体和底部椎体中的一对或多对椎体;也可以根据用户通过输入模块308输入的选择信息确定一对或多对椎体。当根据预设的策略需要分别针对多对椎体计算测量值时,坐标变换模块302可以依次获取多对椎体的第一坐标信息,进行坐标变换得到第二坐标信息,再将各对椎体的第二坐标信息提供给形变计算模块303,触发形变计算模块303分别针对每对椎体的第二坐标信息计算脊柱形变量。
形变计算模块303主要用于利用所述坐标变换模块302提供的第一椎体和第二椎体的第二坐标信息计算第一椎体和第二椎体之间的Cobb角,并将得到的所述Cobb角提供给输出模块304。
图4为Cobb角的计算原理示意图。如图4所示,形变计算模块303利用配准后的模板脊柱X线图像中的第一椎体401上边缘的若干个点的第二坐标值确定第一椎体上边缘的切线,利用第二椎体402下边缘的若干个点的第二坐标值确定第二椎体下边缘的切线,计算这两条切线之间的夹角403作为Cobb角。
根据一个实施例,坐标变换模块302还可以将第一椎体和第二椎体的信息,如标识、序号、或名称等,提供给形变计算模块303,形变计算模块303将第一椎体和第二椎体的信息连同计算得到的形变量一起提供给所述输出模块。
输出模块304输出所述第一椎体和第二椎体的信息以及所述脊柱形变量。输出模块304可以包括用户界面,用于显示第一椎体和第二椎体的信息以及所述脊柱形变量,还可以显示原始的、或进行了分段和标记后的目标脊柱X线图像。
上述各模块仅是根据功能进行划分得到的逻辑模块,实现时可以由同一模块实现,也可以分由多个子模块实现。以上仅列出与本发明方案有关的功能模块,实际上,该系统还可以包括其它的模块。例如,该装置还可以包括:中央处理单元(CPU)、存储模块、内部总线等。如前所述,上述模块的功能在实现时可能需要其它模块协助,例如利用CPU的计算功能,利用存储模块存储计算的中间量等。简洁起见,均描述为由上述模块完成。总线则提供各模块间的信息交互通道,可以是一根连接各模块的总线,也可以由上述各个模块间的多条线路组成。当某些模块,例如目标图像数据库306、模板图像数据库307等与其它模块位于不同的设备中时,总线与这些模块之间可能需要通过有线或无线的链路进行连接,还有可能通过本地网络或互联网络相连。例如目标图像数据库306、模板图像数据库307可以由一个服务器实现,当其它模块需要读取其中的数据时,可以通过某种通信协议通过网络连接到服务器来获取该数据。
图5为本发明实施例提供的一种脊柱X线图像的处理方法的流程图。如图5所示,该方法可以包括以下步骤。
步骤501,从目标图像数据库选择目标脊柱X线图像。
步骤502,从模板图像数据库选择模板脊柱X线图像。
目标图像数据库中可以存储有目标脊柱X线图像对应的被测量者的信息,如年龄、性别、人种等,以及图像的采集角度,如侧面图、正面图等。选定了目标脊柱X线图像后,可以获取相应的被测量者的信息和图像采集角度等。
模板图像数据库中可以存储多个模板图像,同时也存储各模板图像对应的信息,如适用年龄、性别、人种及图像采集角度等。
可以按照预定的策略,根据被测量者的年龄、性别、人种、所述目标脊柱X线图像对应的图像采集角度等信息中的一个或多个来选择所述模板脊柱X线图像。例如,可以仅考虑年龄和图像采集角度时,如果待测量者的年龄为12岁,检查的项目为脊柱侧弯,则目标脊柱X线图像为被测量者的正面或背面X线图像,并可以从模板图像数据库中选择对应年龄10-15岁的背面脊柱X线图像。
步骤503,根据目标脊柱X线图像对模板脊柱X线图像进行非刚性配准,得到空间变换信息。
空间变换信息可以是模板脊柱X线图像在配准前后多个点的位移矢量,例如可以包括空间位移和旋转角度等。
步骤504,确定待测量的第一椎体和第二椎体。可以直接读取预定的待测量椎体的信息,也可以根据用户输入的选择信息确定待测量的第一椎体和第二椎体。
步骤505,获得预先存储的第一椎体和第二椎体在模板脊柱X线图像中的坐标信息。其中,坐标信息可以为原始模板脊柱X线图像中第一椎体的第一边缘和第二椎体的第二边缘上的若干点的坐标。
步骤506,根据空间变换信息对坐标信息进行坐标变换,得到配准后的所述若干点的坐标信息。根据一个实施例,例如采用直角坐标系,第一椎体的第一边缘上的一个点的坐标是(2,2),经过配准,得到该点的空间位移是5,旋转角度是30度,则可以计算得到配准后,该点的坐标变为(5,6)。
步骤507,根据配准后的坐标信息计算第一椎体和第二椎体之间的Cobb角。
步骤508,输出计算得到的Cobb角和第一椎体、第二椎体的信息。
步骤509,判断用户是否输入待测量的第一椎体和第二椎体,如果是,则执行步骤505,否则结束计算。
通过上述实施例中的流程,就可以得出被测量者的第一椎体和第二椎体之间的Cobb角,从而可以根据Cobb角来对被测量者的脊柱形变情况进行定量地评估。如果需要进一步获得更详细的测量结果,可以重复步骤504到508,在步骤504中输入需要进一步测量的两个椎体,在步骤508得到输出的这两个椎体间的Cobb角。
需要说明的是,上述各流程和各系统结构图中不是所有的步骤和模块都是必须的,可以根据实际的需要忽略某些步骤或模块。各步骤的执行顺序不是固定的,可以根据需要进行调整。上述各实施例中描述的系统结构可以是物理结构,也可以是逻辑结构,即,有些模块可能由同一物理实体实现,或者,有些模块可能分由多个物理实体实现,或者,可以由多个独立设备中的某些部件共同实现。
另外,上面描述中采用“第一”、“第二”仅仅为了方便区分具有同一含义的两个对象,并不表示其有实质的区别。
以上各实施例中,硬件单元可以通过机械方式或电气方式实现。例如,一个硬件单元可以包括永久性专用的电路或逻辑(如专门的处理器,FPGA或ASIC)来完成相应操作。硬件单元还可以包括可编程逻辑或电路(如通用处理器或其它可编程处理器),可以由软件进行临时的设置以完成相应操作。具体的实现方式(机械方式、或专用的永久性电路、或者临时设置的电路)可以基于成本和时间上的考虑来确定。
本发明还提供了一种机器可读介质,存储用于使一机器执行如本文所述的图像处理方法的指令。具体地,可以提供配有存储介质的系统或者装置,在该存储介质上存储着实现上述实施例中任一实施例的功能的软件程序代码,且使该系统或者装置的计算机(或CPU或MPU)读出并执行存储在存储介质中的程序代码。
在这种情况下,从存储介质读取的程序代码本身可实现上述实施例中任何一项实施例的功能,因此程序代码和存储程序代码的存储介质构成了本发明的一部分。
用于提供程序代码的存储介质实施例包括软盘、硬盘、磁光盘、光盘(如CD-ROM、CD-R、CD-RW、DVD-ROM、DVD-RAM、DVD-RW、DVD+RW)、磁带、非易失性存储卡和ROM。可选择地,可以由通信网络从服务器计算机上下载程序代码。
此外,应该清楚的是,不仅可以通过执行计算机所读出的程序代码,而且可以通过基于程序代码的指令使计算机上操作的操作系统等来完成部分或者全部的实际操作,从而实现上述实施例中任意一项实施例的功能。
此外,可以理解的是,将由存储介质读出的程序代码写到插入计算机内的扩展板中所设置的存储器中或者写到与计算机相连接的扩展单元中设置的存储器中,随后基于程序代码的指令使安装在扩展板或者扩展单元上的CPU等来执行部分和全部实际操作,从而实现上述实施例中任一实施例的功能。
Claims (17)
1.一种脊柱X线图像的处理方法,包括:
利用第二脊柱X线图像对预定的第一脊柱X线图像进行配准,得到所述第一脊柱X线图像的空间变换信息;
根据所述空间变换信息和所述第一脊柱X线图像中位于至少一对椎体边缘的预定点的第一坐标信息,获得配准后所述预定点的第二坐标信息;
利用所述第二坐标信息计算脊柱形变量。
2.根据权利要求1所述的方法,其特征在于,
所述空间变换信息包括配准后所述预定点的位移矢量;
所述根据所述空间变换信息和所述第一脊柱X线图像中预定点的第一坐标信息获得配准后所述预定点的第二坐标信息包括:利用所述位移矢量和所述第一坐标信息对所述预定点的第一坐标进行坐标变换,得到所述第二坐标信息。
3.根据权利要求1所述的方法,其特征在于,
所述预定点分别位于所述第一脊柱X线图像中第一椎体和第二椎体的边缘;
所述利用所述第二坐标信息计算脊柱形变量包括:
根据所述预定点的第二坐标信息确定所述第一椎体的第一边缘的切线和所述第二椎体的第二边缘的切线,计算所述切线之间的夹角。
4.根据权利要求1所述的方法,其特征在于,
所述预定点分别位于第一脊柱X线图像中第一椎体和第二椎体的边缘;
所述利用所述第二坐标信息计算脊柱形变量包括:
根据所述预定点的第二坐标信息计算所述第一椎体和第二椎体之间的Cobb角。
5.根据权利要求3或4所述的方法,其特征在于,所述第一椎体和第二椎体为以下组合中的任一个:
脊柱的顶部椎体和底部椎体,颈椎中的顶部椎体和底部椎体,胸椎中的顶部椎体和底部椎体,腰椎中的顶部椎体和底部椎体,用户选择的一对椎体。
6.根据权利要求1所述的方法,其特征在于,进一步包括:
从预先存储的至少两个脊柱X线图像中选择一个作为所述第一脊柱X线图像。
7.根据权利要求6所述的方法,其特征在于,所述从预先存储的至少两个脊柱X线图像中选择一个作为所述第一脊柱X线图像包括:
根据年龄、性别、人种、图像采集角度中的一个或多个从所述至少两个脊柱X线图像中选择一个作为所述第一脊柱X线图像。
8.根据权利要求1所述的方法,其特征在于,进一步包括:
输出所述脊柱形变量和所述至少一对椎体的信息。
9.一种脊柱X线图像的处理系统,包括:
图像配准模块,用于利用第二脊柱X线图像对预定的第一脊柱X线图像进行配准得到所述第一脊柱X线图像的空间变换信息;
坐标变换模块,用于根据所述图像配准模块获得的所述空间变换信息和所述第一脊柱X线图像中位于至少一对椎体边缘的预定点的第一坐标信息获得配准后所述预定点的第二坐标信息;
形变计算模块,用于利用所述第二坐标信息计算脊柱形变量。
10.根据权利要求9所述的系统,其特征在于,
所述图像配准模块用于获得配准后所述预定点的位移矢量作为所述空间变换信息;
所述坐标变换模块用于利用所述图像配准模块获得的所述位移矢量和所述第一坐标信息对所述预定点的第一坐标进行坐标变换,得到所述第二坐标信息。
11.根据权利要求9所述的系统,其特征在于,
所述预定点分别位于所述第一脊柱X线图像中第一椎体和第二椎体的边缘;
所述形变计算模块用于根据所述预定点的第二坐标信息确定所述第一椎体的第一边缘的切线和所述第二椎体的第二边缘的切线,计算所述切线之间的夹角作为所述脊柱形变量。
12.根据权利要求9所述的系统,其特征在于,
所述预定点分别位于所述第一脊柱X线图像中第一椎体和第二椎体的边缘;
所述形变计算模块用于根据所述预定点的第二坐标信息计算第一脊柱X线图像中第一椎体和第二椎体之间的Cobb角作为所述脊柱形变量。
13.根据权利要求11或12所述的系统,其特征在于,所述第一椎体和第二椎体为以下组合中的任一个:
脊柱的顶部椎体和底部椎体,颈椎中的顶部椎体和底部椎体,胸椎中的顶部椎体和底部椎体,腰椎中的顶部椎体和底部椎体,用户选择的一对椎体。
14.根据权利要求9所述的系统,其特征在于,进一步包括:
图像选择模块,用于从存储模块中预先存储的至少两个脊柱X线图像中选择一个作为所述第一脊柱X线图像提供给所述图像配准模块。
15.根据权利要求14所述的系统,其特征在于,
所述图像选择模块用于根据年龄、性别、人种、图像采集角度中的一个或多个从所述至少两个脊柱X线图像中选择一个作为所述第一脊柱X线图像提供给所述图像配准模块。
16.根据权利要求9所述的系统,其特征在于,进一步包括:
输出模块,用于输出所述脊柱形变量和所述至少一对椎体的信息。
17.一种机器可读介质,在其上存储有机器可执行指令,其中,当所述机器可执行指令被执行时,使得机器实施权利要求1-8中的任意一个所包括的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210271466.6A CN103565449A (zh) | 2012-07-31 | 2012-07-31 | 一种脊柱x线图像的处理方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210271466.6A CN103565449A (zh) | 2012-07-31 | 2012-07-31 | 一种脊柱x线图像的处理方法和系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103565449A true CN103565449A (zh) | 2014-02-12 |
Family
ID=50038775
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210271466.6A Pending CN103565449A (zh) | 2012-07-31 | 2012-07-31 | 一种脊柱x线图像的处理方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103565449A (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105719273A (zh) * | 2014-12-05 | 2016-06-29 | Ge医疗系统环球技术有限公司 | 一种在医学影像上测量脊椎的旋转参数的方法及装置 |
CN107361792A (zh) * | 2016-05-12 | 2017-11-21 | 宗晓宇 | 一种中医诊疗设备及其使用方法 |
CN108030495A (zh) * | 2017-11-10 | 2018-05-15 | 深圳纳富特科技有限公司 | 身体姿势的测量方法、装置、保健器械及计算机可读介质 |
CN108309334A (zh) * | 2017-12-08 | 2018-07-24 | 李书纲 | 一种脊柱x线影像的数据处理方法 |
CN108573502A (zh) * | 2018-03-06 | 2018-09-25 | 安徽大学 | 一种自动测量Cobb角的方法 |
CN109493334A (zh) * | 2018-11-12 | 2019-03-19 | 深圳码隆科技有限公司 | 测量脊柱弯曲的方法及装置 |
CN110037699A (zh) * | 2019-04-23 | 2019-07-23 | 上海联影医疗科技有限公司 | 用于磁共振扫描的病床组件及应用磁共振系统扫描的方法 |
CN110072436A (zh) * | 2016-12-08 | 2019-07-30 | 皇家飞利浦有限公司 | 脊柱医学成像数据的简化导航 |
CN112472111A (zh) * | 2020-12-14 | 2021-03-12 | 上海交通大学医学院附属新华医院 | 一种脊柱侧凸早期自检与康复云平台 |
US11145060B1 (en) | 2020-07-20 | 2021-10-12 | International Business Machines Corporation | Automatic detection of vertebral dislocations |
CN115187606A (zh) * | 2022-09-14 | 2022-10-14 | 中国医学科学院北京协和医院 | 一种青少年特发性脊柱侧凸pumc分型方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1449721A (zh) * | 2002-04-03 | 2003-10-22 | 株式会社东芝 | 具有医学图像测量功能的医学图像处理装置 |
CN1969756A (zh) * | 2005-09-07 | 2007-05-30 | 美国西门子医疗解决公司 | 用图像和测角来计算机辅助检测脊柱弯曲度的系统和方法 |
EP1968015A1 (en) * | 2004-11-10 | 2008-09-10 | Agfa Healthcare | Method of performing measurements on digital images |
WO2011021181A1 (en) * | 2009-08-16 | 2011-02-24 | Ori Hay | Assessment of spinal anatomy |
-
2012
- 2012-07-31 CN CN201210271466.6A patent/CN103565449A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1449721A (zh) * | 2002-04-03 | 2003-10-22 | 株式会社东芝 | 具有医学图像测量功能的医学图像处理装置 |
EP1968015A1 (en) * | 2004-11-10 | 2008-09-10 | Agfa Healthcare | Method of performing measurements on digital images |
CN1969756A (zh) * | 2005-09-07 | 2007-05-30 | 美国西门子医疗解决公司 | 用图像和测角来计算机辅助检测脊柱弯曲度的系统和方法 |
WO2011021181A1 (en) * | 2009-08-16 | 2011-02-24 | Ori Hay | Assessment of spinal anatomy |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105719273A (zh) * | 2014-12-05 | 2016-06-29 | Ge医疗系统环球技术有限公司 | 一种在医学影像上测量脊椎的旋转参数的方法及装置 |
CN107361792B (zh) * | 2016-05-12 | 2020-05-19 | 宗晓宇 | 一种中医诊疗设备及其使用方法 |
CN107361792A (zh) * | 2016-05-12 | 2017-11-21 | 宗晓宇 | 一种中医诊疗设备及其使用方法 |
CN110072436B (zh) * | 2016-12-08 | 2022-12-13 | 皇家飞利浦有限公司 | 脊柱医学成像数据的简化导航 |
CN110072436A (zh) * | 2016-12-08 | 2019-07-30 | 皇家飞利浦有限公司 | 脊柱医学成像数据的简化导航 |
CN108030495A (zh) * | 2017-11-10 | 2018-05-15 | 深圳纳富特科技有限公司 | 身体姿势的测量方法、装置、保健器械及计算机可读介质 |
CN108030495B (zh) * | 2017-11-10 | 2020-09-04 | 深圳市瑞贝特科技有限公司 | 身体姿势的测量方法、装置、保健器械及计算机可读介质 |
CN108309334A (zh) * | 2017-12-08 | 2018-07-24 | 李书纲 | 一种脊柱x线影像的数据处理方法 |
CN108573502A (zh) * | 2018-03-06 | 2018-09-25 | 安徽大学 | 一种自动测量Cobb角的方法 |
CN108573502B (zh) * | 2018-03-06 | 2021-07-06 | 安徽大学 | 一种自动测量Cobb角的方法 |
CN109493334B (zh) * | 2018-11-12 | 2020-12-29 | 深圳码隆科技有限公司 | 测量脊柱弯曲的方法及装置 |
CN109493334A (zh) * | 2018-11-12 | 2019-03-19 | 深圳码隆科技有限公司 | 测量脊柱弯曲的方法及装置 |
CN110037699A (zh) * | 2019-04-23 | 2019-07-23 | 上海联影医疗科技有限公司 | 用于磁共振扫描的病床组件及应用磁共振系统扫描的方法 |
CN110037699B (zh) * | 2019-04-23 | 2024-01-19 | 上海联影医疗科技股份有限公司 | 应用磁共振系统扫描的方法 |
US11145060B1 (en) | 2020-07-20 | 2021-10-12 | International Business Machines Corporation | Automatic detection of vertebral dislocations |
CN112472111A (zh) * | 2020-12-14 | 2021-03-12 | 上海交通大学医学院附属新华医院 | 一种脊柱侧凸早期自检与康复云平台 |
CN112472111B (zh) * | 2020-12-14 | 2023-08-29 | 上海交通大学医学院附属新华医院 | 一种脊柱侧凸早期自检与康复云平台 |
CN115187606A (zh) * | 2022-09-14 | 2022-10-14 | 中国医学科学院北京协和医院 | 一种青少年特发性脊柱侧凸pumc分型方法 |
CN115187606B (zh) * | 2022-09-14 | 2022-12-13 | 中国医学科学院北京协和医院 | 一种青少年特发性脊柱侧凸pumc分型方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103565449A (zh) | 一种脊柱x线图像的处理方法和系统 | |
Wu et al. | Automated comprehensive adolescent idiopathic scoliosis assessment using MVC-Net | |
JP4744883B2 (ja) | 画像位置合わせ方法及び医用画像データ処理装置 | |
US8509502B2 (en) | Spine labeling | |
US8498459B2 (en) | System and method for verifying registration accuracy in digital medical images | |
US8965083B2 (en) | Automatic detection of vertebrae boundaries in spine images | |
CN110175502A (zh) | 一种脊柱Cobb角测量方法、装置、可读存储介质及终端设备 | |
JP2018530401A5 (zh) | ||
US6668083B1 (en) | Deriving geometrical data of a structure from an image | |
US8897517B2 (en) | Imaging system with reporting function and method of operation thereof | |
US9576392B2 (en) | Image display device and method, and medium containing program | |
US9881125B2 (en) | Ultrasound measurement of biometrics of fetus | |
CN111311655B (zh) | 多模态图像配准方法、装置、电子设备、存储介质 | |
Joshi et al. | Automatic intra‐subject registration‐based segmentation of abdominal fat from water–fat MRI | |
EP3490482B1 (en) | System and method for verification of fiducial correspondence during image-guided surgical procedures | |
US20120027277A1 (en) | Interactive iterative closest point algorithm for organ segmentation | |
CN106295110B (zh) | 医学图像处理方法和图像处理系统 | |
CN103838951B (zh) | 标注放射学图像的方法和设备 | |
CN112529860A (zh) | 脊柱图像处理方法、装置、电子设备及存储介质 | |
JP2019520954A (ja) | 脊柱の弯曲角度を測定する方法及びデバイス | |
EP3089107B1 (en) | Computer program product and method for determining lesion similarity of medical image | |
Zhang et al. | An automated cobb angle estimation method using convolutional neural network with area limitation | |
Alukaev et al. | A deep learning framework for vertebral morphometry and Cobb angle measurement with external validation | |
CN111782529A (zh) | 辅助诊断系统的测试方法、装置、计算机设备及存储介质 | |
CN112535489B (zh) | 一种利用波形图测量脊柱侧弯Cobb角的方法及其系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140212 |