一种脊柱侧弯椎间刚度个性化计算方法技术

本发明专利技术提供了一种脊柱侧弯椎间刚度个性化计算方法，包括：步骤1、根据患者卧位CT图像提取几何、材料信息并初步计算椎间刚度；步骤2、构建患者整个脊柱的多体模型，计算站位的脊柱长度；步骤3、用临床测得的站位脊柱长度校正椎间刚度矩阵。本发明专利技术可以较为准确地估测脊柱侧弯患者地椎间刚度，为脊柱侧弯患者的手术多体仿真提供个性化的刚度矩阵设置，提高仿真的准确性，为临床手术策略的制订提供高效的指导。导。导。

【技术实现步骤摘要】
一种脊柱侧弯椎间刚度个性化计算方法


[0001]本专利技术涉及脊柱力学仿真的
，尤其涉及脊柱侧弯椎间刚度个性化计算方法。

技术介绍

[0002]脊柱侧弯是指脊柱的一个或数个节段向侧方弯曲并伴有椎体旋转和矢状面上后凸或前凸的增加或减少的脊柱畸形,是一种三维畸形。严重的脊柱侧弯不但会影响患者外观，还会导致脏器、神经等压迫，导致出现心肺功能衰竭、截瘫等并发症，给患者带来巨大的痛苦。对于严重的脊柱侧弯，临床上一般采用手术的方法，将患者的脊柱与两根矫形棒固定，来达到矫形的效果。但是目前脊柱手术的术前规划主要依赖医生的主观经验，难以准确预测矫形手术的效果。借助多体模型，对患者的手术过程进行个性化建模，可以快速求解术后Cobb角的变化以及椎体和钉棒受力，从而对手术效果和风险进行全面准确的评估。确定椎体间的刚度，是脊柱矫形手术多体仿真的重要环节。
[0003]临床上医生通过让患者弯腰对脊柱整体的柔韧度进行估测，但是脊柱的多体仿真需要精准的椎间刚度。目前没有可行的方法通过在体实验的方式测量椎间刚度。
[0004]鉴于这种现状，目前的脊柱多体仿真，大多采用固定的椎间刚度，没有考虑患者间椎间刚度的差异。而不同的脊柱侧弯患者脊柱几何形态、组织材料特性差异非常大，因此，采用固定的刚度设置无法准确完成脊柱侧弯手术的多体仿真。

技术实现思路

[0005]有鉴于此，本专利技术提供了一种脊柱侧弯椎间刚度个性化计算方法，可以为脊柱侧弯患者的手术多体仿真提供个性化的刚度矩阵设置，提高仿真的准确性。
[0006]步骤1、根据患者卧位CT图像提取几何、材料信息并初步计算椎间刚度。几何信息为相邻椎体的尺寸、相邻椎体的角度，材料信息为椎间盘区域的CT灰度值，初步计算采用的刚度计算公式为：
[0007]其中M0为标准椎体的刚度矩阵，β为相邻椎体终板平面的角度，K为角度校正矩阵，α为椎体尺寸，α0为标准椎体尺寸，γ为椎间盘CT灰度值。
[0008]步骤2、构建患者整个脊柱的多体模型，计算站位的脊柱长度。采用步骤1计算所得的刚度矩阵，每节椎体的载荷为所承受的重力，分析的类型为静力学分析，脊柱长度为T1节段到L5节段在竖直方向的距离。
[0009]步骤3、用临床测得的站位脊柱长度校正椎间刚度矩阵。
[0010]优选的，步骤一中患者椎体几何信息的提取方法为：在相邻两椎体的两个终板分别选取8个点，用最小二乘法分别拟合两个平面，计算两个平面之间的角度，椎体尺寸为两椎体宽度的平均值。
[0011]优选的，步骤三中校正椎间刚度矩阵的方法为：整个脊柱每个椎间刚度矩阵乘以
一个相同的系数。
[0012]优选的，步骤3中采用的迭代方法为二分法，如果初始仿真刚度矩阵偏小，则系数乘以2倍，直到刚度矩阵过大。
[0013]优选的，所述步骤3中的收敛条件为：仿真值与真实值差别小于3mm。
附图说明：
[0015]图1为本专利技术实施例的整体流程图；
[0016]图2为本专利技术实施例的选取标记点示意图；
[0017]图3本专利技术实施例的刚度矩阵初步计算流程图；
[0018]图4本专利技术实施例的椎体标记点示意图；
[0019]图5本专利技术实施例的椎体自身坐标系；
[0020]图6本专利技术实施例的模型重建流程图；
[0021]图7本专利技术实施例的模型重建效果图；
[0022]图8本专利技术实施例的站立态仿真对比示意；
[0023]图9本专利技术实施例的刚度矩阵校正流程图；
[0024]图10本专利技术实施例的术前术后对比示意图。
具体实施方式：
[0026]为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本专利技术进一步详细说明。
[0027]实施例。一种脊柱侧弯椎间刚度个性化计算方法，如图1所示，包括步骤1、通过卧位CT图像对患者椎体的几何尺寸和椎间盘CT值进行测量，基于椎间刚度矩阵与椎间几何尺寸、材料的关系对患者的椎间刚度进行初步的估算；步骤2、根据初步估算的脊柱椎间刚度建立患者的脊柱多体动力学模型，对患者的站位和卧位姿态进行模拟，计算站位的脊柱长度；步骤3、基于医学影像的真实脊柱长度对比，通过仿真与实际的差异整体调整患者的椎间刚度，使仿真与实际脊柱长度相一致，实现对椎间刚度的校准。
[0028]所述步骤1中几何信息为相邻椎体的尺寸、相邻椎体的角度，具体方法为：在相邻两椎体的两个终板分别选取8个点，这8个点要求均匀分散在整个终板平面，以减小拟合平面的误差，用最小二乘法分别拟合两个平面，基于两个平面的法向量计算两个平面之间的角度，椎体尺寸为两椎体宽度的平均值。如图2所示，为本专利技术一个实施例的选取标记点示意。如图3所示，为本专利技术一个实施例的刚度矩阵初步计算流程。根据患者椎间盘的CT灰度值确定椎间盘材料的弹性模量，椎间盘弹性模量与椎间盘CT灰度值呈线性关系。根据如下公式初步计算椎体间刚度矩阵：
[0029]其中M0为标准椎体的刚度矩阵，为6*6的矩阵，来源于实验测量的数据；β为相邻椎体终板平面在冠状面的角度；K为角度校正矩阵，同样为6*6的矩阵，体现椎体冠状面夹角对刚度矩阵的影响；α为椎体尺寸；α0为标准椎体尺寸，即椎体间刚度矩阵与椎体的尺寸为线性关系；γ为椎间盘CT灰度值。
[0030]该公式的特点在于：考虑了患者脊柱几何形态和材料属性对椎体间刚度矩阵的影响。患者的12节胸椎和5节腰椎均可通过以上公式进行初步计算。刚度矩阵表示的是相邻椎体之间位移旋转和力、力矩之间的关系：F＝MX，
[0031]其中F代表椎体所受三个方向的力和力矩，X代表椎体在三个方向上的位移和旋转，二者在一定范围内可看作是线性关系。
[0032]所述步骤2根据初步估算的脊柱椎间刚度建立患者的脊柱多体动力学模型，对患者的站位和卧位姿态进行模拟，计算站位的脊柱长度，其模型重建流程如图6所示。具体仿真步骤为：从患者CT图像提取患者12节胸椎、5节腰椎和骶骨的特征标记点，每个椎体提取4个标记点，分别为两个椎弓根的中点以及椎体上下终板的中心点，骶骨上终板提取3个标记点，这些标记点足以准确计算椎体的空间位置和姿态。如图4所示，为本专利技术一个实施例的椎体标记点示意图。
[0034]图5所示，本专利技术实施例的椎体自身坐标系。将椎体上下终板中点坐标求平均值，可求椎体中心点的坐标；同理可以求出骶骨上终板中心点坐标；将椎体中心点的坐标减去骶骨上终板中心点坐标，可求得每个椎体的相对位置。
[0035]通过标记点的坐标，可以计算椎体在各个方向上的倾斜角，由椎体上下终板的中点可以计算椎体在冠状面和矢状面的倾斜角；由椎体椎弓根的两个标记点可以计算椎体在横断面上的倾斜角；将这三个倾斜角进行组合，可以计算椎体自身坐标轴的方向，然后通过自身坐标轴的方向分步计算椎体的旋转矩阵，最后相乘得到椎体总体的旋转矩阵。V
rot
＝cosθν+(1
‑
cosθ)(v
·
k)k+sinθk
本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种脊柱侧弯椎间刚度个性化计算方法，其特征在于，具体包括：步骤1：根据患者卧位CT图像提取几何、材料信息并初步计算椎间刚度；步骤2：构建患者整个脊柱的多体模型，计算站位的脊柱长度；步骤3：用临床测得的站位脊柱长度校正椎间刚度矩阵。2.根据权利要求1所述的脊柱侧弯椎间刚度个性化计算方法，其特征在于，所述步骤1中几何信息为相邻椎体的尺寸、相邻椎体的角度，具体方法为：在相邻两椎体的两个终板分别选取8个点，用最小二乘法分别拟合两个平面，计算两个平面之间的角度，椎体尺寸为两椎体宽度的平均值。3.根据权利要求1所述的脊柱侧弯椎间刚度个性化计算方法，其特征在于，所述步骤1中材料信息为椎间盘区域的CT值。4.根据权利要求1所述的脊柱侧弯椎间刚度个性化计算方法，其特征在于，所述步骤1采用的刚度计算公式为：其中M0为标准椎体的刚度矩阵，β为相邻椎体终板平面的角度，K为角度校正矩阵，α为椎体尺寸，α0...

【专利技术属性】
技术研发人员：陈文轩，张明峥，蒲放，弓明侠，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：