一种基于SPECT成像的人体部位切分方法及系统

摘要

本发明涉及一种基于SPECT成像的人体部位切分方法及系统，基于骨显像数据矩阵，在骨显像图像中提取人体区域，以获取骨显像数据矩阵的有效信息，不仅减少了计算量，而且提高了切分精度。然后基于人体区域对应的数据矩阵和人体形态学特征，将人体区域划分为头部区域、胸腔区域、盆骨区域和腿部区域，再基于胸腔区域对应的数据矩阵，将胸腔区域划分为躯干区域、左上肢和右上肢，最后基于腿部区域对应的数据矩阵，将腿部区域划分为左下肢和右下肢，从而将人体划分为头部、躯干、四肢、盆骨等部位，实现了人体部位的精准切分，能够有效针对任意一个SPCET骨显像图像进行头部、躯干、四肢、盆骨等部位的精准切分。

权利要求

1.一种基于SPECT成像的人体部位切分方法，其特征在于，所述切分方法包括如下步骤：

利用SPECT成像方法进行骨扫描，得到骨显像图像，并获取所述骨显像图像中的前位图或后位图所对应的骨显像数据矩阵；

基于所述骨显像数据矩阵，在所述骨显像图像中提取人体区域；

基于所述人体区域对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特征点，将所述人体区域划分为头部区域、胸腔区域、盆骨区域和腿部区域；

基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢；

基于所述腿部区域对应的数据矩阵，选取双腿划分特征点，将所述腿部区域划分为左下肢和右下肢；

其中，所述基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢具体包括：将所述胸腔区域对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘制第二统计图；对所述第二统计图进行曲线拟合，得到第二拟合曲线；

选取所述第二拟合曲线的最大值点作为脊柱特征点；根据所述脊柱特征点将所述胸腔区域对应的数据矩阵划分为左胸腔对应的数据矩阵和右胸腔对应的数据矩阵；

将所述左胸腔对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘制第三统计图；对所述第三统计图进行曲线拟合，得到第三拟合曲线；

根据所述第三拟合曲线确定左胳膊切分特征点，将所述左胸腔对应的数据矩阵划分为左躯干对应的数据矩阵和左上肢对应的数据矩阵；

将所述右胸腔对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘制第四统计图；对所述第四统计图进行曲线拟合，得到第四拟合曲线；

根据所述第四拟合曲线确定右胳膊切分特征点，将所述右胸腔对应的数据矩阵划分为右躯干对应的数据矩阵和右上肢对应的数据矩阵；

合并所述左躯干对应的数据矩阵和所述右躯干对应的数据矩阵，得到躯干区域对应的数据矩阵；

将所述躯干区域对应的数据矩阵、所述左上肢对应的数据矩阵和所述右上肢对应的数据矩阵分别在显示界面上进行可视化，得到躯干区域、左上肢和右上肢。

2.如权利要求1所述的一种基于SPECT成像的人体部位切分方法，其特征在于，在所述骨显像图像中提取人体区域之前，还包括对所述骨显像数据矩阵进行固定阈值去噪，具体包括：

获取历史骨显像图像对应的所有历史数据矩阵；

基于所有所述历史数据矩阵的放射值，以放射值作为横坐标、放射值的个数作为纵坐标绘制直方图；

根据所述直方图确定阈值；

根据所述阈值对所述骨显像数据矩阵进行去噪，将小于所述阈值的放射值归为背景值，得到去噪后的骨显像数据矩阵。

3.如权利要求2所述的一种基于SPECT成像的人体部位切分方法，其特征在于，所述根据所述直方图确定阈值具体包括：在所述直方图中选取与背景值相差最小的放射值作为阈值。

4.如权利要求1或2所述的一种基于SPECT成像的人体部位切分方法，其特征在于，所述基于所述骨显像数据矩阵，在所述骨显像图像中提取人体区域具体包括：基于所述骨显像数据矩阵，从上到下按行进行遍历，选取第一个包含非背景值的行作为起始行；所述非背景值为放射值不为背景值的放射值；

基于所述骨显像数据矩阵，从下到上按行进行遍历，选取第一个包含非背景值的行作为终止行；

选取所述起始行与所述终止行以内的所有行，构成人体区域对应的数据矩阵；

将所述人体区域对应的数据矩阵在显示界面上进行可视化，得到人体区域。

5.如权利要求1所述的一种基于SPECT成像的人体部位切分方法，其特征在于，所述基于所述人体区域对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特征点，将所述人体区域划分为头部区域、胸腔区域、盆骨区域和腿部区域具体包括：

对于所述人体区域对应的数据矩阵，按行进行非背景值计数，得到每一行所包含的非背景值的个数；

以行数作为横坐标，非背景值的个数作为纵坐标，绘制第一统计图；

对所述第一统计图进行曲线拟合，得到第一拟合曲线；

根据人体形态学特征和所述第一拟合曲线，选取第一切分特征点、第二切分特征点和第三切分特征点；

根据所述第一切分特征点、所述第二切分特征点和所述第三切分特征点，将所述人体区域对应的数据矩阵划分为头部区域对应的数据矩阵、胸腔区域对应的数据矩阵、盆骨区域对应的数据矩阵和腿部区域对应的数据矩阵；

将所述头部区域对应的数据矩阵、所述胸腔区域对应的数据矩阵、所述盆骨区域对应的数据矩阵和所述腿部区域对应的数据矩阵分别在显示界面上进行可视化，得到头部区域、胸腔区域、盆骨区域和腿部区域。

6.如权利要求5所述的一种基于SPECT成像的人体部位切分方法，其特征在于，所述根据人体形态学特征和所述第一拟合曲线，选取第一切分特征点、第二切分特征点和第三切分特征点具体包括：

根据人体形态学特征，在所述第一拟合曲线中选取处于第一预设范围内的极小值或最小值作为第一切分特征点，选取处于第二预设范围内的极小值或最小值作为第二切分特征点，选取处于第三预设范围内的极小值或最小值作为第三切分特征点。

7.如权利要求6所述的一种基于SPECT成像的人体部位切分方法，其特征在于，所述根据所述第一切分特征点、所述第二切分特征点和所述第三切分特征点，将所述人体区域对应的数据矩阵划分为头部区域对应的数据矩阵、胸腔区域对应的数据矩阵、盆骨区域对应的数据矩阵和腿部区域对应的数据矩阵具体包括：所述人体区域对应的数据矩阵的起始行和所述第一切分特征点所在行以内的所有行构成头部区域对应的数据矩阵；

所述第一切分特征点所在行和所述第二切分特征点所在行以内的所有行构成胸腔区域对应的数据矩阵；

所述第二切分特征点所在行和所述第三切分特征点所在行以内的所有行构成盆骨区域对应的数据矩阵；

所述第三切分特征点所在行和所述人体区域对应的数据矩阵的终止行以内的所有行构成腿部区域对应的数据矩阵；

或者，

所述人体区域对应的数据矩阵的起始行和所述第一切分特征点所在行以内的所有行构成头部区域对应的数据矩阵；

将所述第二切分特征点所在行上移β行，得到第一实际切分行；所述第一切分特征点所在行和所述第一实际切分行以内的所有行构成胸腔区域对应的数据矩阵；

将所述第二切分特征点所在行下移β行，得到第二实际切分行；将所述第三切分特征点所在行上移r行，得到第三实际切分行；所述第二实际切分行和所述第三实际切分行以内的所有行构成盆骨区域对应的数据矩阵；

所述第三切分特征点所在行和所述人体区域对应的数据矩阵的终止行以内的所有行构成腿部区域对应的数据矩阵。

8.如权利要求1所述的一种基于SPECT成像的人体部位切分方法，其特征在于，所述基于所述腿部区域对应的数据矩阵，选取双腿划分特征点，将所述腿部区域划分为左下肢和右下肢具体包括：

选取所述腿部区域对应的数据矩阵的对称列作为双腿划分特征点所在列；

根据所述双腿划分特征点所在列将所述腿部区域对应的数据矩阵划分为左下肢对应的数据矩阵和右下肢对应的数据矩阵；

将所述左下肢对应的数据矩阵和所述右下肢对应的数据矩阵分别在显示界面上进行可视化，得到左上肢和右下肢。

9.一种基于SPECT成像的人体部位切分系统，其特征在于，所述切分系统包括：骨显像数据矩阵获取模块，用于利用SPECT成像方法进行骨扫描，得到骨显像图像，并获取所述骨显像图像中的前位图或后位图所对应的骨显像数据矩阵；

人体区域提取模块，用于基于所述骨显像数据矩阵，在所述骨显像图像中提取人体区域；

人体区域划分模块，用于基于所述人体区域对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特征点，将所述人体区域划分为头部区域、胸腔区域、盆骨区域和腿部区域；

胸腔区域划分模块，用于基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢；

所述基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢具体包括：将所述胸腔区域对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘制第二统计图；对所述第二统计图进行曲线拟合，得到第二拟合曲线；

选取所述第二拟合曲线的最大值点作为脊柱特征点；根据所述脊柱特征点将所述胸腔区域对应的数据矩阵划分为左胸腔对应的数据矩阵和右胸腔对应的数据矩阵；

将所述左胸腔对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘制第三统计图；对所述第三统计图进行曲线拟合，得到第三拟合曲线；

根据所述第三拟合曲线确定左胳膊切分特征点，将所述左胸腔对应的数据矩阵划分为左躯干对应的数据矩阵和左上肢对应的数据矩阵；

将所述右胸腔对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘制第四统计图；对所述第四统计图进行曲线拟合，得到第四拟合曲线；

根据所述第四拟合曲线确定右胳膊切分特征点，将所述右胸腔对应的数据矩阵划分为右躯干对应的数据矩阵和右上肢对应的数据矩阵；

合并所述左躯干对应的数据矩阵和所述右躯干对应的数据矩阵，得到躯干区域对应的数据矩阵；

将所述躯干区域对应的数据矩阵、所述左上肢对应的数据矩阵和所述右上肢对应的数据矩阵分别在显示界面上进行可视化，得到躯干区域、左上肢和右上肢；

腿部区域划分模块，用于基于所述腿部区域对应的数据矩阵，选取双腿划分特征点，将所述腿部区域划分为左下肢和右下肢。

说明书

一种基于SPECT成像的人体部位切分方法及系统

技术领域

[0001] 本发明涉及医学影像技术领域，特别是涉及一种基于SPECT成像的人体部位切分方法及系统。

背景技术

[0002] 医学影像是现代临床医学疾病诊治的重要技术手段，它能够以非侵入方式捕获和呈现机体的内部结构和功能状态。医学影像作为一般性术语，通常包含医学成像(Medical 
imaging)和医学图像处理(Medical image processing)两个相对独立的构成部分。

[0003] 在医学影像领域，核医学影像被视为一种诊断、治疗、评价和预防各种严重疾病的有效医疗方式。自1895年英国物理学家伦琴发现X射线以来，医学成像从计算机断层扫描
(Computerized Tomography，CT)、核磁共振成像(Magnetic Resonance Imaging，MRI)、超
声波成像(Ultrasound)等传统结构成像模态发展到目前的单光子发射计算机断层成像
(Single‑Photon Emission Computed Tomography，SPECT)和正电子发射断层成像
(Positron Emission Tomography，PET)等功能成像模态，进而催生了SPECT/CT、SPECT/
MRI、PET/CT和PET/MRI等混合成像模态。

[0004] SPECT成像不同于传统X射线成像，传统X射线成像只提供一个器官或身体局部区域的解剖信息，而SPECT成像可以显示在整个身体的器官和组织中的结构和功能变异。
SPECT成像需要预先在患者体内注入放射性药物(如99Tmc‑MIBI)，经过一定时长的代谢后，
用体外探头捕获体内特定部位或全身的药物残留分布，即放射量。大量临床验证表明，病变
部位的放射量通常高于正常部位，这种依据放射量而非器官形态差异实现疾病检测的方
法，有效弥补了传统结构成像的不足。这是因为，有些疾病引发的器官或组织形态改变可能
显著滞后于疾病的真正发作。目前，SPECT核医学成像已经应用于帕金森等精神疾病、甲状
腺、脑部疾病等的计算机自动辅助诊断领域。

[0005] 目前，全身骨扫描显像的诊断主要依靠医师的诊断，这要求临床医生不仅需要掌握大量的肿瘤影像学知识，更需要深入分析SPECT成像原理及其病灶，临床医生不仅需要掌
握肿瘤的生物学行为规律，还要分析病灶的产生并且及时做出诊断。近几年，基于全身骨扫
描的计算机辅助诊断技术以图像预处理、病灶点分割、多视图特征融合分类、假阳性案例排
查等研究为主。SPECT显像原理限制了传统基于分类、分割方法在全身骨扫描显像中的处理
分析，在SPECT成像中，人体骨骼的放射值和背景值相差较小，导致极易发生显像时骨骼和
背景的粘连现象，并且存在体外放射污染，这使得人体骨骼在SPECT显像中显像不清楚，骨
骼与背景的信噪比低、骨骼和软组织的边界不明显、不同年龄的不同患者之间的SPECT显像
各异，而且患者在SPECT扫描时形态各异，这使得基于计算机辅助诊断技术进行SPECT显像
分析非常困难。目前，尚未发现可有效针对每一个SPCET全身骨显像数据矩阵进行头部、躯
干、四肢、盆骨等的精准切分的方法。

发明内容

[0006] 本发明的目的是提供一种基于SPECT成像的人体部位切分方法及系统，在不改变骨显像数据矩阵放射值的基础上，实现了对任意一个骨显像数据矩阵的精准切分。

[0007] 为实现上述目的，本发明提供了如下方案：

[0008] 一种基于SPECT成像的人体部位切分方法，所述切分方法包括如下步骤：

[0009] 利用SPECT成像方法进行骨扫描，得到骨显像图像，并获取所述骨显像图像中的前位图或后位图所对应的骨显像数据矩阵；

[0010] 基于所述骨显像数据矩阵，在所述骨显像图像中提取人体区域；

[0011] 基于所述人体区域对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特征点，将所述人体区域划分为头部区域、胸腔区域、盆骨区域和
腿部区域；

[0012] 基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢；

[0013] 基于所述腿部区域对应的数据矩阵，选取双腿划分特征点，将所述腿部区域划分为左下肢和右下肢。

[0014] 一种基于SPECT成像的人体部位切分系统，所述切分系统包括：

[0015] 骨显像数据矩阵获取模块，用于利用SPECT成像方法进行骨扫描，得到骨显像图像，并获取所述骨显像图像中的前位图或后位图所对应的骨显像数据矩阵；

[0016] 人体区域提取模块，用于基于所述骨显像数据矩阵，在所述骨显像图像中提取人体区域；

[0017] 人体区域划分模块，用于基于所述人体区域对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特征点，将所述人体区域划分为头部区
域、胸腔区域、盆骨区域和腿部区域；

[0018] 胸腔区域划分模块，用于基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢；

[0019] 腿部区域划分模块，用于基于所述腿部区域对应的数据矩阵，选取双腿划分特征点，将所述腿部区域划分为左下肢和右下肢。

[0020] 根据本发明提供的具体实施例，本发明公开了以下技术效果：

[0021] 本发明所提供的一种基于SPECT成像的人体部位切分方法及系统，先利用SPECT成像方法进行骨扫描，得到全身骨显像可视化图像，并获取骨显像图像中的前位图或后位图
所对应的骨显像数据矩阵，并基于骨显像数据矩阵，在骨显像图像中提取人体区域，以获取
骨显像数据矩阵的有效信息，不仅减少了计算量，而且提高了切分精度。然后基于人体区域
对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特
征点，将人体区域划分为头部区域、胸腔区域、盆骨区域和腿部区域，再基于胸腔区域对应
的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将胸腔区域划分为躯干区域、左
上肢和右上肢，最后基于腿部区域对应的数据矩阵，选取双腿划分特征点，将腿部区域划分
为左下肢和右下肢，从而将人体划分为头部、躯干、四肢、盆骨等部位，实现了人体部位的精
准切分，能够有效针对任意一个形态、姿势正常的SPCET骨显像图像进行头部、躯干、四肢、
盆骨等部位的精准切分。

附图说明

[0022] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施
例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获
得其他的附图。

[0023] 图1为本发明实施例1所提供的切分方法的方法流程图。

[0024] 图2为本发明实施例1所提供的数据矩阵的三维示意图；(a)为包括病灶的数据矩阵的三维示意图，其中病灶点为骨转移；(b)为包含注射点的数据矩阵的三维示意图。

[0025] 图3为本发明实施例1所提供的人体区域提取的方法流程图。

[0026] 图4为本发明实施例1所提供的人体区域提取原理图。

[0027] 图5为本发明实施例1所提供的人体区域划分的方法流程图。

[0028] 图6为本发明实施例1所提供的人体区域划分原理图。

[0029] 图7为本发明实施例1所提供的人体区域划分的结果示意图。

[0030] 图8为本发明实施例1所提供的胸腔区域划分的方法流程图。

[0031] 图9为本发明实施例1所提供的胸腔区域划分原理图。

[0032] 图10为本发明实施例1所提供的左胸腔划分原理图。

[0033] 图11为本发明实施例1所提供的右胸腔划分原理图。

[0034] 图12为本发明实施例1所提供的腿部区域划分的方法流程图。

[0035] 图13为本发明实施例1所提供的腿部区域划分原理图。

[0036] 图14为本发明实施例1所提供的人体部位切分的结果示意图。

[0037] 图15为本发明实施例2所提供的切分系统的结构框图。

具体实施方式

[0038] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于
本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他
实施例，都属于本发明保护的范围。

[0039] 本发明的目的是提供一种基于SPECT成像的人体部位切分方法及系统，在不改变骨显像数据矩阵放射值的基础上，实现了对任意一个骨显像数据矩阵的精准切分。

[0040] 为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

[0041] 实施例1：

[0042] 现有技术主要针对传统自然图像的人体分割，不仅算法复杂性高，而且无法直接应用于骨显像图像。在面向全身骨显像图像的切分中，尚未发现一款可有效针每一个全身
骨显像图像的精准人体切分方法。本实施例用于提供一种基于SPECT成像的人体部位切分
方法，能够适用于每一个骨显像图像的人体切分，如图1所示，所述切分方法包括如下步骤：

[0043] S1：利用SPECT成像方法进行骨扫描，得到骨显像图像，并获取所述骨显像图像中的前位图或后位图所对应的骨显像数据矩阵；

[0044] 具体的，S1包括：通过核素探测器捕获核素量，基于计算机辅助设备显示为骨显像图像，并获得骨显像图像的DICOM文件。由于DICOM文件中信息量大，故再通过专业的计算机
技术在DICOM文件中获取骨显像图像对应的数据矩阵，以得到前位图或后位图所对应的骨
显像数据矩阵。如图2所示，其给出了数据矩阵的三维示意图，纵轴代表数据矩阵中的放射
值。本实施例利用后续步骤对前位图对应的骨显像数据矩阵或者后位图对应的骨显像数据
矩阵进行处理，均可以实现对人体的精确切分。

[0045] 由于利用SPECT成像方法所得到的骨显像数据矩阵，会存在由于体外药物残留所导致的噪声数据。作为一种可选的实施方式，在骨显像图像中提取人体区域之前，本实施例
所提供的切分方法还包括对骨显像数据矩阵进行固定阈值去噪，需要说明的是，骨显像数
据矩阵的矩阵大小为1024×256，骨显像数据矩阵中的每个元素的值均是一个放射值。去噪
具体包括：获取历史骨显像图像对应的所有历史数据矩阵，再基于所有历史数据矩阵的放
射值，以放射值作为横坐标、放射值的个数作为纵坐标绘制直方图。根据直方图确定阈值，
具体的，绘制直方图后，在直方图中选取与背景值相差最小的放射值作为阈值，背景值为0，
即选取最接近0的放射值作为阈值。当然，本实施例也可以不绘制直方图，直接基于所有历
史数据矩阵的放射值选取与背景值相差最小的放射值作为阈值即可。在得到阈值后，再根
据阈值对骨显像数据矩阵进行去噪，将小于阈值的放射值归为背景值，得到去噪后的骨显
像数据矩阵。阈值可为5，对于骨显像数据矩阵中的每一个元素，若这个元素的放射值小于
5，则将其放射值归为背景值，即归为0，否则，则不改变该放射值。

[0046] 通过对骨显像数据矩阵进行固定阈值去噪，能够消除显像过程中，体外残留的少量显像剂颗粒被SPECT捕获而形成体外放射值污染所产生的噪声。

[0047] S2：基于所述骨显像数据矩阵，在骨显像图像中提取人体区域；

[0048] 当然，本实施例也可以基于去噪后的骨显像数据矩阵，在骨显像图像中提取人体区域。

[0049] 具体的，如图3和图4所示，S2可以包括：

[0050] S21：基于所述骨显像数据矩阵，从上到下按行进行遍历，选取第一个包含非背景值的行作为起始行；所述非背景值为放射值不为背景值的放射值，即选取第一个包含不为0
的值的行作为起始行。

[0051] S22：基于所述骨显像数据矩阵，从下到上按行进行遍历，选取第一个包含非背景值的行作为终止行；

[0052] S23：选取所述起始行与所述终止行以内的所有行，构成人体区域对应的数据矩阵，进而能够获取骨显像数据矩阵中的有效信息，不仅减少了参数量和计算量，而且提高了
切分精度，高效辅助医生诊断。

[0053] S24：将所述人体区域对应的数据矩阵在显示界面上进行可视化，得到人体区域。

[0054] 本实施例均是对数据矩阵进行处理，而数据矩阵可视化后即为人体图像。进而在将骨显像数据矩阵进行处理得到人体区域对应的数据矩阵后，再将人体区域对应的数据矩
阵进行可视化，即可得到人体区域的图像，这一过程即为在骨显像图像中提取人体区域的
过程。

[0055] S3：基于所述人体区域对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特征点，将所述人体区域划分为头部区域、胸腔区域、盆骨区域
和腿部区域；

[0056] 具体的，如图5和图6所示，S3可以包括：

[0057] S31：对于所述人体区域对应的数据矩阵，按行进行非背景值计数，得到每一行所包含的非背景值的个数；

[0058] S32：以行数作为横坐标，非背景值的个数作为纵坐标，绘制第一统计图；第一统计图参见图6。

[0059] S33：对所述第一统计图进行曲线拟合，得到第一拟合曲线；具体的，将第一统计图中的统计点基于numpy库进行多项式拟合，并绘制曲线，得到第一拟合曲线。

[0060] S34：根据人体形态学特征和所述第一拟合曲线，选取第一切分特征点、第二切分特征点和第三切分特征点；

[0061] 参见图6，a点即为第一切分特征点，b点即为第二切分特征点，c点即为第三切分特征点。由于第一统计图是以行数作为横坐标，故第一拟合曲线上的特征点即对应于数据矩
阵中的一行。为了便于说明，本实施例所提供的图6将特征点所在行指明为特征点。

[0062] 具体的，S34包括：根据人体形态学特征，在第一拟合曲线中选取处于第一预设范围内的极小值或最小值作为第一切分特征点，第一预设范围为人体区域对应的数据矩阵的
第10％‑20％行，数据矩阵可视化即为人体图像，即相当于人体长度的10％‑20％；选取处于
第二预设范围内的极小值或最小值作为第二切分特征点，第二预设范围为人体区域对应的
数据矩阵的第25％‑45％行；选取处于第三预设范围内的极小值或最小值作为第三切分特
征点，第三预设范围为第二切分特征点所在行至人体区域对应的数据矩阵的第70％行，从
而根据人体的形态学特征，确定人体区域划分的特征点，即确定人体区域对应的数据矩阵
的特征行。

[0063] S35：根据所述第一切分特征点、所述第二切分特征点和所述第三切分特征点，将所述人体区域对应的数据矩阵划分为头部区域对应的数据矩阵、胸腔区域对应的数据矩
阵、盆骨区域对应的数据矩阵和腿部区域对应的数据矩阵；

[0064] 作为一种可选的实施方式，S35可以包括：人体区域对应的数据矩阵的起始行和第一切分特征点所在行以内的所有行构成头部区域对应的数据矩阵；第一切分特征点所在行
和第二切分特征点所在行以内的所有行构成胸腔区域对应的数据矩阵；第二切分特征点所
在行和第三切分特征点所在行以内的所有行构成盆骨区域对应的数据矩阵；第三切分特征
点所在行和人体区域对应的数据矩阵的终止行以内的所有行构成腿部区域对应的数据矩
阵，进而将人体区域对应的数据矩阵划分为头部区域对应的数据矩阵、胸腔区域对应的数
据矩阵、盆骨区域对应的数据矩阵和腿部区域对应的数据矩阵。

[0065] 考虑到骨显像图像的特点以及人体区域划分过程中曲线拟合造成的误差，本实施例还提供了另外一种对人体区域对应的数据矩阵进行划分的实施方式，人体区域对应的数
据矩阵的起始行和第一切分特征点所在行以内的所有行构成头部区域对应的数据矩阵；考
虑到腰椎长度，将第二切分特征点所在行上移β行，得到第一实际切分行；第一切分特征点
所在行和第一实际切分行以内的所有行构成胸腔区域对应的数据矩阵；考虑到手的影响，
将第二切分特征点所在行下移β行，得到第二实际切分行；将第三切分特征点所在行上移r
行，得到第三实际切分行；第二实际切分行和第三实际切分行以内的所有行构成盆骨区域
对应的数据矩阵；第三切分特征点所在行和人体区域对应的数据矩阵的终止行以内的所有
行构成腿部区域对应的数据矩阵，进而通过改变特征点所在行，能够减小误差。具体的，β的
取值范围为(5，10)，r的取值范围为(10，20)。

[0066] S36：将所述头部区域对应的数据矩阵、所述胸腔区域对应的数据矩阵、所述盆骨区域对应的数据矩阵和所述腿部区域对应的数据矩阵分别在显示界面上进行可视化，得到
头部区域、胸腔区域、盆骨区域和腿部区域。

[0067] 上述将相应数据矩阵进行可视化的过程，即为将人体区域划分为头部区域、胸腔区域、盆骨区域和腿部区域的过程。参见图7，其给出了人体区域划分的结果示意图，(a)为
头部区域，(b)为胸腔区域，(c)为盆骨区域，(d)为腿部区域。

[0068] S4：基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢；

[0069] 具体的，参见图8，S4可以包括：

[0070] S41：将所述胸腔区域对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐
标，绘制第二统计图；对所述第二统计图进行曲线拟合，得到第二拟合曲线；

[0071] 图9给出了第二统计图及第二拟合曲线的示意图。

[0072] S42：选取所述第二拟合曲线的最大值点作为脊柱特征点；根据所述脊柱特征点将所述胸腔区域对应的数据矩阵划分为左胸腔对应的数据矩阵和右胸腔对应的数据矩阵；

[0073] 在得到脊柱特征点后，由于在确定脊柱特征点时对胸腔区域对应的数据矩阵进行逆时针旋转，故对于胸腔区域对应的原始数据矩阵，以脊柱特征点所在列作为分割列，对胸
腔区域对应的原始数据矩阵进行划分，得到左胸腔对应的数据矩阵和右胸腔对应的数据矩
阵。

[0074] S43：将所述左胸腔对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘
制第三统计图；对所述第三统计图进行曲线拟合，得到第三拟合曲线；

[0075] 图10给出了第三统计图和第三拟合曲线的示意图。

[0076] S44：根据所述第三拟合曲线确定左胳膊切分特征点，将所述左胸腔对应的数据矩阵划分为左躯干对应的数据矩阵和左上肢对应的数据矩阵；

[0077] 以旋转之后的左胸腔对应的数据矩阵的行数作为高度Y，左胳膊切分特征点的取值范围位于20％Y‑70％Y，选取这一范围内的最小值点作为左胳膊切分特征点。在得到左胳
膊切分特征点后，由于在确定左胳膊切分特征点时对左胸腔对应的数据矩阵进行逆时针旋
转，故对于左胸腔对应的原始数据矩阵，以左胳膊切分特征点所在列作为分割列，对左胸腔
对应的原始数据矩阵进行划分，得到左躯干对应的数据矩阵和左上肢对应的数据矩阵。

[0078] 另外，为了进一步减小由于曲线拟合带来的误差，也为了使左胳膊包含更多信息，本实施例在进行左胸腔对应的数据矩阵的划分时，对于左胸腔对应的原始数据矩阵，将左
胳膊切分特征点所在列右移α行作为分割列，对左胸腔对应的原始数据矩阵进行划分，得到
左躯干对应的数据矩阵和左上肢对应的数据矩阵。其中，α的取值范围为(0，5)。

[0079] S45：将所述右胸腔对应的数据矩阵逆时针旋转90度，并按行进行非背景值计数，得到每一行所包含的非背景值的个数；以行数作为横坐标，非背景值的个数作为纵坐标，绘
制第四统计图；对所述第四统计图进行曲线拟合，得到第四拟合曲线；

[0080] 图11给出了第四统计图和第四拟合曲线的示意图。

[0081] S46：根据所述第四拟合曲线确定右胳膊切分特征点，将所述右胸腔对应的数据矩阵划分为右躯干对应的数据矩阵和右上肢对应的数据矩阵；

[0082] 以旋转之后的右胸腔对应的数据矩阵的行数作为高度Y，右胳膊切分特征点的取值范围位于20％Y‑70％Y，选取这一范围内的最小值点作为右胳膊切分特征点。在得到右胳
膊切分特征点后，由于在确定右胳膊切分特征点时对右胸腔对应的数据矩阵进行逆时针旋
转，故对于右胸腔对应的原始数据矩阵，以右胳膊切分特征点所在列作为分割列，对右胸腔
对应的原始数据矩阵进行划分，得到右躯干对应的数据矩阵和右上肢对应的数据矩阵。

[0083] 另外，为了进一步减小由于曲线拟合带来的误差，也为了使右胳膊包含更多信息，本实施例在进行右胸腔对应的数据矩阵的划分时，对于右胸腔对应的原始数据矩阵，将右
胳膊切分特征点所在列左移α行作为分割列，对右胸腔对应的原始数据矩阵进行划分，得到
右躯干对应的数据矩阵和右上肢对应的数据矩阵。其中，α的取值范围为(0，5)。

[0084] S47：合并所述左躯干对应的数据矩阵和所述右躯干对应的数据矩阵，得到躯干区域对应的数据矩阵；

[0085] S48：将所述躯干区域对应的数据矩阵、所述左上肢对应的数据矩阵和所述右上肢对应的数据矩阵分别在显示界面上进行可视化，得到躯干区域、左上肢和右上肢。

[0086] 上述将相应数据矩阵进行可视化的过程，即为将胸腔区域划分为躯干区域、左上肢和右上肢的过程，进而对胸腔区域进行进一步的精细划分。

[0087] S5：基于所述腿部区域对应的数据矩阵，选取双腿划分特征点，将所述腿部区域划分为左下肢和右下肢。

[0088] 具体的，如图12所示，S5可以包括：

[0089] S51：选取所述腿部区域对应的数据矩阵的对称列作为双腿划分特征点所在列；

[0090] 考虑到腿部显像较为清晰、对称，因此对腿部区域进行对称切分，将腿部区域对应的数据矩阵的对称列作为双腿划分特征点所在列。

[0091] S52：根据所述双腿划分特征点所在列将所述腿部区域对应的数据矩阵划分为左下肢对应的数据矩阵和右下肢对应的数据矩阵；

[0092] S53：将所述左下肢对应的数据矩阵和所述右下肢对应的数据矩阵分别在显示界面上进行可视化，得到左上肢和右下肢。

[0093] 另外，本实施例在进行腿部区域划分时，还可以先将腿部区域对应的数据矩阵逆时针旋转90度，然后选取旋转之后的腿部区域对应的数据矩阵的对称行作为双腿划分特征
点所在行，根据双腿划分特征点所在行将旋转之后的腿部区域对应的数据矩阵划分为旋转
之后的左下肢对应的数据矩阵和旋转之后的右下肢对应的数据矩阵，再分别将旋转之后的
左下肢对应的数据矩阵和旋转之后的右下肢对应的数据矩阵顺时针旋转90度，得到左下肢
对应的数据矩阵和右下肢对应的数据矩阵，将左下肢对应的数据矩阵和右下肢对应的数据
矩阵分别在显示界面上进行可视化，得到左上肢和右下肢。

[0094] 如图13所示，为了进一步提高左下肢和右下肢提取的准确性，本实施例在进行腿部区域划分时，还可以先将腿部区域对应的数据矩阵逆时针旋转90度，然后基于旋转后的
腿部区域对应的数据矩阵，从上到下按行进行遍历，选取第一个包含非背景值的行作为起
始行，从下到上按行进行遍历，选取第一个包含非背景值的行作为终止行，然后选取起始行
与终止行以内的所有行，构成旋转之后的腿部区域对应的数据矩阵，然后再利用上述方法
对腿部区域对应的数据矩阵进行划分。

[0095] 图14给出了利用本实施例所提供的切分方法进行人体切分后的结果示意图。由此，本实施例所提供的切分方法，基于人体形态学、统计学及多项式曲线拟合进行人体切分
特征点的寻找并进行切分，使得对于每一个全身骨显像矩阵，既不改变原始数据矩阵信息，
又精确切分出了人体部位，该方法解决了现有人工识别切分SPECT全身骨显像图像造成的
低效性及主观性，此方法不仅能有效针对每一个全身骨显像矩阵进行头部、躯干、盆骨、四
肢的精准切分，而且用时短，耗费空间小，具有高效性和高精度的优点。

[0096] 实施例2：

[0097] 本实施例用于提供一种基于SPECT成像的人体部位切分系统，如图15所示，所述切分系统包括：

[0098] 骨显像数据矩阵获取模块M1，用于利用SPECT成像方法进行骨扫描，得到骨显像图像，并获取所述骨显像图像中的前位图或后位图所对应的骨显像数据矩阵；

[0099] 人体区域提取模块M2，用于基于所述骨显像数据矩阵，在所述骨显像图像中提取人体区域；

[0100] 人体区域划分模块M3，用于基于所述人体区域对应的数据矩阵和人体形态学特征，选取第一切分特征点、第二切分特征点和第三切分特征点，将所述人体区域划分为头部
区域、胸腔区域、盆骨区域和腿部区域；

[0101] 胸腔区域划分模块M4，用于基于所述胸腔区域对应的数据矩阵，选取左胳膊切分特征点和右胳膊切分特征点，将所述胸腔区域划分为躯干区域、左上肢和右上肢；

[0102] 腿部区域划分模块M5，用于基于所述腿部区域对应的数据矩阵，选取双腿划分特征点，将所述腿部区域划分为左下肢和右下肢。

[0103] 本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统
而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说
明即可。

[0104] 本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据
本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不
应理解为对本发明的限制。