CN112017147B - 一种骨骼图像的提取方法及装置 - Google Patents

Abstract

本申请提供一种骨骼图像的提取方法及装置；该提取方法可以包括：提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像；沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度；从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度；将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像。

Description

一种骨骼图像的提取方法及装置

技术领域

本申请涉及图像处理技术领域，尤其涉及一种骨骼图像的提取方法及装置。

背景技术

X光图像是骨骼手术中常用的成像方式，随着计算机辅助治疗在临床手术中的应用，如何将X光图像有效利用，提高脊骨骼手术的精确性和安全性，成为热门的研究趋势。

在研究和应用X光图像时，往往需要对成像目标进行目的区域的提取。例如，在脊柱手术中，通常需要提取出X光图像中的脊柱区域，一方面可降低X光图像的噪声，另一方面还可为脊柱分节提供更精确的输入数据，从而提高后续针对该脊柱的X光图像与CT图像进行配准的精度。

发明内容

有鉴于此，本申请提供一种骨骼图像的提取方法及装置。

为实现上述目的，本申请提供技术方案如下：

根据本申请的第一方面，提出了一种骨骼图像的提取方法，包括：

提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像；

沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度；

从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度；

将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像。

可选的，所述目标图像为针对所述目标骨骼进行X光成像得到的X光图像；所述提取目标图像的前景部分，包括：

对所述X光图像进行二值化处理，以提取出所述X光图像的前景部分。

可选的，所述沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度，包括：

确定所述前景部分在各投影方向下像素点数量超过预设数量阈值的投影部分；

将所述投影部分的像素点分布宽度作为所述前景部分在各投影方向下的像素点分布宽度。

可选的，所述从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度，包括：

从各投影方向下的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，

所述沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度，包括：按照预设角度间隔对所述前景部分环绕投影360°，并将与各投影方向下的像素点分布宽度中最短的像素点分布宽度对应的投影方向作为初步投影方向θ；以上一次得到的初步投影方向θ为中心，在[θ-Δθ，θ+Δθ]范围内按照预设角度间隔对所述前景部分进行投影，直到针对所述前景部分进行投影的循环次数达到预设循环阈值为止；其中，每次循环过程中投影的角度间隔小于上一次循环过程中投影的角度间隔，每次循环过程中Δθ的取值小于上一次循环过程中Δθ的取值；

所述从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度，包括：获取最后一次循环过程中各投影方向下的像素点分布宽度，从获取到的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，每次循环过程中Δθ的取值为上一次循环过程中投影的角度间隔。

根据本申请的第二方面，提出了一种骨骼图像的提取装置，包括：

前景提取单元，提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像；

投影单元，沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度；

选取单元，从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度；

骨骼提取单元，将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像。

可选的，所述目标图像为针对所述目标骨骼进行X光成像得到的X光图像；所述前景提取单元具体用于：

对所述X光图像进行二值化处理，以提取出所述X光图像的前景部分。

可选的，所述投影单元具体用于：

确定所述前景部分在各投影方向下像素点数量超过预设数量阈值的投影部分；

将所述投影部分的像素点分布宽度作为所述前景部分在各投影方向下的像素点分布宽度。

可选的，所述选取单元具体用于：

从各投影方向下的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，

投影单元具体用于：按照预设角度间隔对所述前景部分环绕投影360°，并将与各投影方向下的像素点分布宽度中最短的像素点分布宽度对应的投影方向作为初步投影方向θ；以上一次得到的初步投影方向θ为中心，在[θ-Δθ，θ+Δθ]范围内按照预设角度间隔对所述前景部分进行投影，直到针对所述前景部分进行投影的循环次数达到预设循环阈值为止；其中，每次循环过程中投影的角度间隔小于上一次循环过程中投影的角度间隔，每次循环过程中Δθ的取值小于上一次循环过程中Δθ的取值；

所述选取单元具体用于：获取最后一次循环过程中各投影方向下的像素点分布宽度，从获取到的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，每次循环过程中Δθ的取值为上一次循环过程中投影的角度间隔。

根据本申请的第三方面，提供一种电子设备，包括：

处理器；

用于存储处理器可执行指令的存储器；

其中，所述处理器通过运行所述可执行指令以实现如上述实施例中任一所述的骨骼图像的提取方法。

根据本公开实施例的第四方面，提供一种计算机可读存储介质，其上存储有计算机指令，该指令被处理器执行时实现如上述实施例中任一所述骨骼图像的提取方法的步骤。

由以上技术方案可见，本申请基于同一目标骨骼的延伸方向和宽度较为固定的特点，通过在多个不同投影方向上对目标图像的前景部分进行投影，可以根据投影所得图像中的像素点分布宽度来确定出目标骨骼的骨骼宽度和延伸方向，从而根据确定出的骨骼宽度和延伸方向来提取出目标骨骼的骨骼图像。

附图说明

图1是本申请一示例性实施例示出的一种骨骼图像的提取方法的流程图。

图2是本申请一示例性实施例示出的另一种骨骼图像的提取方法的流程图。

图3是本申请一示例性实施例示出的提取前景部分的示意图。

图4A-4B是本申请一示例性实施例示出的投影示意图。

图5是本申请一示例性实施例示出的提取脊柱区域的示意图。

图6是本申请一示例性实施例示出的一种电子设备的结构示意图。

图7是本申请一示例性实施例示出的一种骨骼图像的提取装置的框图。

具体实施方式

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本申请相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本申请的一些方面相一致的装置和方法的例子。

在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。

应当理解，尽管在本申请可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本申请范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。

请参见图1，图1是本申请一示例性实施例示出的一种骨骼图像的提取方法的流程图。如图1所示，该方法可以包括以下步骤：

步骤102，提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像。

在本实施例中，对目标骨骼进行成像得到的目标图像包括前景部分和背景部分，其中前景部分包含目标骨骼的骨骼图像，而本申请正是提供一种骨骼图像的提取方案，以从目标图像中提取出包含目标骨骼的骨骼区域(即骨骼图像)。

进一步的，以目标图像为针对目标骨骼进行X光成像得到的X光图像为例，可对该X光图像进行二值化处理，以提取出该X光图像中的前景部分。

步骤104，沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度。

在本实施例中，在沿多个不同投影方向对前景部分进行投影得到投影图像后，理想情况下投影图像中只有对应于目标骨骼的骨骼区域才会产生投影值(包括横坐标像素点分布宽度，以及纵坐标像素点数量)，也即只有对应于骨骼区域的投影值非零。但是，在实际情况中，除目标骨骼的骨骼区域以外的其他部分的投影值也可能非零。以提取出脊柱的前景部分为例，该前景部分可能还包含靠近脊柱的其他骨骼(或者其他人体组织)，或者一些在成像过程中产生的噪点，导致前景部分也不可避免地包含了部分非脊柱区域。

对此，本申请针对上述误差设定一预设数量阈值th，只保留各投影方向下像素点数量大于数量阈值th的部分。作为一示例性实施例，可先确定前景部分在各投影方向下像素点数量超过预设数量阈值的投影部分，再将该投影部分的像素点分布宽度作为前景部分在各投影方向下的像素点分布宽度。可见，本申请通过设定数量阈值th，可有效排除不属于目标骨骼区域的其他区域造成的干扰，从而提高后续测量目标骨骼的骨骼宽度的精度。当然，数量阈值th的具体取值可由开发人员根据实际情况灵活设定，本申请并不对此进行限制。

步骤106，从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度。

在一实施例中，可从各投影方向下的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为目标骨骼的骨骼宽度。其中，可通过缩小各个投影方向之间角度的间隔来提高测量骨骼宽度的精度。

在另一实施例中，可先以每次投影的角度间隔较大的投影方式，对前景部分进行较为粗略的环绕投影，然后再以环绕投影得到的像素点分布宽度为依据，进一步确定出范围更小的投影方向范围，并以角度间隔更小的投影方式，在该范围更小的投影方向范围中进行投影，进而得到更为精确的骨骼宽度。

具体而言，按照预设角度间隔对前景部分环绕投影360°，并将与各投影方向下的像素点分布宽度中最短的像素点分布宽度对应的投影方向作为初步投影方向θ；以上一次得到的初步投影方向θ为中心，在[θ-Δθ，θ+Δθ]范围内按照预设角度间隔对前景部分进行投影，直到针对前景部分进行投影的循环次数达到预设循环阈值为止；其中，每次循环过程中投影的角度间隔小于上一次循环过程中投影的角度间隔，每次循环过程中Δθ的取值小于上一次循环过程中Δθ的取值。基于上述投影得到的像素点分布宽度，可获取最后一次循环过程中各投影方向下的像素点分布宽度，从获取到的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为目标骨骼的骨骼宽度。

进一步的，Δθ的初始值可设定为在环绕投影阶段采用的角度间隔，而每次循环过程中Δθ的取值可设定为上一次循环过程中投影的角度间隔。通过上述设定Δθ取值的方式，可在保证骨骼宽度的测量精度的情况下，有效减少投影次数和计算量，从而提升提取骨骼图像的效率。当然，上述Δθ和循环阈值的具体取值可由开发人员根据实际情况灵活设定，只要符合上述实施例中的设定规则即可，本申请并不对此进行限制。

步骤108，将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像。

需要说明的是，本申请的骨骼图像的提取方案，可应用于任何用于提取骨骼图像的电子设备中，例如电脑、X光机、CT设备等，本申请并不对此进行限制。

由以上技术方案可见，本申请基于同一目标骨骼的延伸方向和宽度较为固定的特点，通过在多个不同投影方向上对目标图像的前景部分进行投影，可以根据投影所得图像中的像素点分布宽度来确定出目标骨骼的骨骼宽度和延伸方向，从而根据确定出的骨骼宽度和延伸方向来提取出目标骨骼的骨骼图像。

为了便于理解，下面以提取X光图像中的脊柱区域为例，结合附图对本申请的骨骼图像请提取方案进行详细描述。

请参见图2，图2是本申请一示例性实施例示出的另一种骨骼图像的提取方法的流程图。如图2所示，该方法可以包括以下步骤：

步骤202，提取X光图像的二值图像。

在本实施例中，可对X光图像(对脊柱进行X光成像得到的X光图像)进行二值化处理，从而提取出该X光图像的前景部分(即步骤202中的二值图像)。其中，所采用的二值化算法可参考相关技术中的记载，在此不再赘述。例如，可采用大金法、双峰法、P参数法、迭代法和OTSU法等；当然，本申请并不对此进行限制。

举例而言，如图3所示，对脊柱的X光图像30进行二值化处理后，可得到二值图像31。通过对X光图像进行二值化处理，可得到包含脊柱区域的前景部分，并且前景部分(二值图像31中的黑色部分)与背景部分(二值图像31中的白色部分)呈现出明显的黑白效果，从而使得图像中的数据量大大减少，进而可凸显出脊柱的轮廓、纹理等特征。

步骤204，对二值图像进行环绕投影。

在本实施例中，可从任意投影方向(即任意角度)开始，按照预设角度间隔对二值图像进行360°的环绕投影(也即投影一周)，那么对应于每个投影方向，均可生成一投影图像。

举例而言，如图4A所示，以θ1的投影方向(即与水平方向呈θ1度)对二值图像31进行投影，可得到图4A所示的投影图像(投影分布坐标图中的图像)。其中，图中投影分布坐标图的横坐标表示前景部分的像素点分布宽度(比如，以一个像素点的宽度为单位来表示)，纵坐标表示前景部分在投影方向上像素点的数量。同理，如图4B所示，以θ2的投影方向(即与水平方向呈θ2度)对二值图像31进行投影，可得到图4B所示的投影图像。

由上述针对二值图像31的投影方式可知：投影图像的横坐标范围可表示二值图像31中前景部分(即脊柱区域)在相应投影方向下所得投影的宽度。因此，在所有投影方向下所得投影的宽度之中，数值最小的宽度可用于表示脊柱的宽度。其中，脊柱的宽度可理解为“当把脊柱抽象为圆柱体时，该圆柱体的直径便为脊柱的宽度”；当然，也可参考相关技术中其他对脊柱宽度的定义，本申请并不对此进行限制。

步骤206，测量脊柱投影图像的宽度。

步骤208，计算最短宽度和初步投影方向θ。

在本实施例中，在沿各个投影方向对二值图像31进行投影得到投影图像后，理想情况下投影图像中只有脊柱区域才会产生投影值(包括横坐标像素点分布宽度，以及纵坐标像素点数量)，也即只有对应于脊柱区域的投影值非零。但是，在实际情况中，除脊柱区域以外的其他部分的投影值也可能非零。比如，二值图像31还可能包含靠近脊柱的其他骨骼(或者其他人体组织)，或者一些在成像过程中产生的噪点，导致前景部分也不可避免地包含了部分非脊柱区域。例如，图3中的噪点区域310；显然，该部分区域的图像并非脊柱。

对此，本申请针对上述误差设定一预设数量阈值th，只保留各投影方向下像素点数量大于数量阈值th的部分。承接于上述举例，如图4A-4B所示，设定有数量阈值th(即图中与横坐标平行的直线th)，那么，只有投影图像中超过该阈值th的部分才可进入统计范围，而位于该阈值th以下的部分，则不计入统计范围(即通过阈值th来筛选出存在干扰的噪点)。例如，图4A投影图像中超过阈值th的部分的宽度为w1，即脊柱图像在投影方向θ1下的投影图像的宽度为w1；同理，图4B投影图像中超过阈值th的部分的宽度为w2，即脊柱图像在投影方向θ2下的投影图像的宽度为w2。

可见，本申请通过设定数量阈值th，可有效排除不属于脊柱区域的其他区域造成的干扰，从而提高后续测量脊柱的骨骼宽度的精度。当然，数量阈值th的具体取值可由开发人员根据实际情况灵活设定，本申请并不对此进行限制。

步骤210，在[θ-Δθ，θ+Δθ]范围内重新投影。

步骤212，重新计算最短宽度和初步投影方向θ。

步骤214，若投影的循环次数达到循环阈值，则转入步骤216；否则，返回步骤210。

承接于上述举例，在通过步骤208首次计算得到初步投影方向θ后，可进一步基于该初步投影方向θ确定出范围更小的投影方向范围，并以角度间隔更小的投影方式，在该范围更小的投影方向范围中进行投影，进而得到更为精确的骨骼宽度。

举例而言，在步骤204环绕投影时，可每隔5°对二值图像31环绕投影360°，得到对应于各个投影方向的投影宽度(投影图像的宽度)w1、w2、……wi，再计算出最短宽度w＝min{w1,w2,…,wi,…}，并记录下对应于最短宽度w的投影方向(即初步投影方向θ)。再以初步投影方向θ为中心，在[θ-Δθ，θ+Δθ]范围(比如，每次循环过程中Δθ的取值可设定为上一次循环过程中投影的角度间隔，即此时Δθ＝5°)内，每隔1°(或者其他小于5°的角度间隔)重新进行投影得到各个投影方向的投影宽度，同样计算出最短宽度w以及对应于最短宽度w的投影方向(即重复执行步骤210-212)。当循环执行步骤210-212的次数达到预设循环阈值时，跳出循环，将最后一次循环过程中得到的最短宽度作为脊柱宽度，以及将对应于该最短宽度的投影方向作为脊柱的延伸方向。

步骤216，确定脊柱宽度和脊柱延伸方向。

步骤218，根据脊柱宽度和脊柱延伸方向提取脊柱区域。

承接于上述举例，假定通过上述过程最终得到脊柱宽度为d，脊柱延伸方向为a；那么，可构建如图5所示的截取区域320(32中的空白区域)。进一步的，通过截取区域320去截取X光图像30，便可提取出X光图像30中的脊柱区域33。

由以上技术方案可见，本申请基于同一脊柱的延伸方向和宽度较为固定的特点，通过在多个不同投影方向上对脊柱X光图像的前景部分进行投影，可以根据投影所得图像中的像素点分布宽度来确定出脊柱的骨骼宽度和延伸方向，从而根据确定出的骨骼宽度和延伸方向来提取出脊柱区域。

图6示出了根据本申请的一示例性实施例的一种电子设备的示意结构图。请参考图6，在硬件层面，该电子设备包括处理器602、内部总线604、网络接口606、内存608以及非易失性存储器610，当然还可能包括其他业务所需要的硬件。处理器602从非易失性存储器610中读取对应的计算机程序到内存608中然后运行，在逻辑层面上形成骨骼图像的提取装置。当然，除了软件实现方式之外，本申请并不排除其他实现方式，比如逻辑器件抑或软硬件结合的方式等等，也就是说以下处理流程的执行主体并不限定于各个逻辑单元，也可以是硬件或逻辑器件。

请参考图7，在软件实施方式中，该骨骼图像的提取装置可以包括：

前景提取单元71，提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像；

投影单元72，沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度；

选取单元73，从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度；

骨骼提取单元74，将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像。

可选的，所述目标图像为针对所述目标骨骼进行X光成像得到的X光图像；所述前景提取单元71具体用于：

对所述X光图像进行二值化处理，以提取出所述X光图像的前景部分。

可选的，所述投影单元72具体用于：

确定所述前景部分在各投影方向下像素点数量超过预设数量阈值的投影部分；

将所述投影部分的像素点分布宽度作为所述前景部分在各投影方向下的像素点分布宽度。

可选的，所述选取单元73具体用于：

从各投影方向下的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，

投影单元72具体用于：按照预设角度间隔对所述前景部分环绕投影360°，并将与各投影方向下的像素点分布宽度中最短的像素点分布宽度对应的投影方向作为初步投影方向θ；以上一次得到的初步投影方向θ为中心，在范围内按照预设角度间隔对所述前景部分进行投影，直到针对所述前景部分进行投影的循环次数达到预设循环阈值为止；其中，每次循环过程中投影的角度间隔小于上一次循环过程中投影的角度间隔，每次循环过程中的取值小于上一次循环过程中的取值；

所述选取单元73具体用于：获取最后一次循环过程中各投影方向下的像素点分布宽度，从获取到的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，每次循环过程中的取值为上一次循环过程中投影的角度间隔。

上述装置中各个单元的功能和作用的实现过程具体详见上述方法中对应步骤的实现过程，在此不再赘述。

对于装置实施例而言，由于其基本对应于方法实施例，所以相关之处参见方法实施例的部分说明即可。以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本申请方案的目的。本领域普通技术人员在不付出创造性劳动的情况下，即可以理解并实施。

在示例性实施例中，还提供了一种包括指令的非临时性计算机可读存储介质，例如包括指令的存储器，上述指令可由上述骨骼图像的提取装置的处理器执行以实现如上述实施例中任一所述的方法，比如该方法可以包括：

提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像；

沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度；

从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度；

将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像。

可选的，所述目标图像为针对所述目标骨骼进行X光成像得到的X光图像；所述提取目标图像的前景部分，包括：

对所述X光图像进行二值化处理，以提取出所述X光图像的前景部分。

可选的，所述沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度，包括：

确定所述前景部分在各投影方向下像素点数量超过预设数量阈值的投影部分；

将所述投影部分的像素点分布宽度作为所述前景部分在各投影方向下的像素点分布宽度。

可选的，所述从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度，包括：

从各投影方向下的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，

所述沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度，包括：按照预设角度间隔对所述前景部分环绕投影360°，并将与各投影方向下的像素点分布宽度中最短的像素点分布宽度对应的投影方向作为初步投影方向θ；以上一次得到的初步投影方向θ为中心，在[θ-Δθ，θ+Δθ]范围内按照预设角度间隔对所述前景部分进行投影，直到针对所述前景部分进行投影的循环次数达到预设循环阈值为止；其中，每次循环过程中投影的角度间隔小于上一次循环过程中投影的角度间隔，每次循环过程中Δθ的取值小于上一次循环过程中Δθ的取值；

所述从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度，包括：获取最后一次循环过程中各投影方向下的像素点分布宽度，从获取到的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

可选的，每次循环过程中Δθ的取值为上一次循环过程中投影的角度间隔。

其中，所述非临时性计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等，本申请并不对此进行限制。

以上所述仅为本申请的较佳实施例而已，并不用以限制本申请，凡在本申请的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本申请保护的范围之内。

Claims (10)

1.一种骨骼图像的提取方法，其特征在于，包括：

提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像；

沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度；

从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度；

将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像；

所述从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度，包括：

从各投影方向下的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度；

所述沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度，包括：按照预设角度间隔对所述前景部分环绕投影360°，并将与各投影方向下的像素点分布宽度中最短的像素点分布宽度对应的投影方向作为初步投影方向θ；以上一次得到的初步投影方向θ为中心，在范围内按照预设角度间隔对所述前景部分进行投影，直到针对所述前景部分进行投影的循环次数达到预设循环阈值为止；其中，每次循环过程中投影的角度间隔小于上一次循环过程中投影的角度间隔，每次循环过程中/>的取值小于上一次循环过程中/>的取值；

所述从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度，包括：获取最后一次循环过程中各投影方向下的像素点分布宽度，从获取到的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

2.根据权利要求1所述的方法，其特征在于，所述目标图像为针对所述目标骨骼进行X光成像得到的X光图像；所述提取目标图像的前景部分，包括：

对所述X光图像进行二值化处理，以提取出所述X光图像的前景部分。

3.根据权利要求1所述的方法，其特征在于，所述沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度，包括：

确定所述前景部分在各投影方向下像素点数量超过预设数量阈值的投影部分；

将所述投影部分的像素点分布宽度作为所述前景部分在各投影方向下的像素点分布宽度。

4.根据权利要求1所述的方法，其特征在于，每次循环过程中的取值为上一次循环过程中投影的角度间隔。

5.一种骨骼图像的提取装置，其特征在于，包括：

前景提取单元，提取目标图像的前景部分，所述前景部分包含待提取的目标骨骼的骨骼图像；

投影单元，沿多个不同投影方向对所述前景部分进行投影，以得到所述前景部分在各投影方向下的像素点分布宽度；

选取单元，从各投影方向下的像素点分布宽度中选取特定像素点分布宽度，以用于表示所述目标骨骼的骨骼宽度；

骨骼提取单元，将对应于所述特定像素点分布宽度的投影方向作为所述目标骨骼的延伸方向，并根据所述延伸方向与所述骨骼宽度提取所述骨骼图像；

所述选取单元具体用于：从各投影方向下的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度；

所述投影单元具体用于：按照预设角度间隔对所述前景部分环绕投影360°，并将与各投影方向下的像素点分布宽度中最短的像素点分布宽度对应的投影方向作为初步投影方向θ；以上一次得到的初步投影方向θ为中心，在范围内按照预设角度间隔对所述前景部分进行投影，直到针对所述前景部分进行投影的循环次数达到预设循环阈值为止；其中，每次循环过程中投影的角度间隔小于上一次循环过程中投影的角度间隔，每次循环过程中/>的取值小于上一次循环过程中/>的取值；

所述选取单元具体用于：获取最后一次循环过程中各投影方向下的像素点分布宽度，从获取到的像素点分布宽度中选取出最短的像素点分布宽度，并将所选取的像素点分布宽度作为所述目标骨骼的骨骼宽度。

6.根据权利要求5所述的装置，其特征在于，所述目标图像为针对所述目标骨骼进行X光成像得到的X光图像；所述前景提取单元具体用于：

对所述X光图像进行二值化处理，以提取出所述X光图像的前景部分。

7.根据权利要求5所述的装置，其特征在于，所述投影单元具体用于：

确定所述前景部分在各投影方向下像素点数量超过预设数量阈值的投影部分；

将所述投影部分的像素点分布宽度作为所述前景部分在各投影方向下的像素点分布宽度。

8.根据权利要求5所述的装置，其特征在于，每次循环过程中的取值为上一次循环过程中投影的角度间隔。

9.一种电子设备，其特征在于，包括：

处理器；

用于存储处理器可执行指令的存储器；

其中，所述处理器通过运行所述可执行指令以实现如权利要求1-4中任一项所述的方法。

10.一种计算机可读存储介质，其上存储有计算机指令，其特征在于，该指令被处理器执行时实现如权利要求1-4中任一项所述方法的步骤。