[发明专利]一种基于奇异排除的滑膜炎压痛点具体位置定量分析方法

权利要求书

1.一种基于奇异排除的滑膜炎压痛点具体位置定量分析方法，其特征在于：包括如下步骤：

S1、CT点云数据的三维重构；

S2、粗选压痛点位置中奇异区域的排除；

S3、竖直方向精选应力集中区域的确定；

S4、融合水平、竖直两向精选压痛点的水平区域确定；

S5、精选穿刺深度的确定：在S4求解的压痛点综合应力λ_j基础上，获得深度修正系数H_j，

其中，V为腔内容积，λ_j为压痛点综合应力，C(λ_j)_max为所有压痛点综合应力中的最大值；j为编号，j＝1，2，3......M，M为精选压痛点的总数量，

然后根据式(4)对精选穿刺深度D_j进行求解，

D_j＝3H_jλ_j (5)

其中，H_j为深度修正系数，λ_j为压痛点综合应力，j为编号，j＝1，2，3......M，M为精选压痛点的总数量，取精选穿刺深度D_j的最大值作为滑膜炎治疗的最适穿刺深度位置。

2.根据权利要求1所述的基于奇异排除的滑膜炎压痛点具体位置定量分析方法，其特征在于：步骤S1三维重构过程如下：通过CT检查方式对滑膜炎症区域进行扫描，将CT切片数据用SOLIDWORKS软件重构，转化为三维.Prt数据，用ANSYS中的网格划分器对导入的.Prt模型进行网格划分，在前处理模块设置好它的约束条件并施加载荷，对后处理结果求解，得到具体参数，包括：网格修正系数δ，腔内容积V，水平单向应力σ_i，水平粗选压痛点i，每个粗选压痛点的应力集中系数K_i，并确定它的量化区域，在ANSYS软件中对i个粗选压痛点附近的应力状态进行分析，取每个应力集中点的区域位置作为压痛点的具体位置，其中，i＝1，2，3......N，N为粗选压痛点的总数量。

3.根据权利要求2所述的基于奇异排除的滑膜炎压痛点具体位置定量分析方法，其特征在于：步骤S2排除过程如下：在S1中重构的三维.Prt数据基础上回归处理，得到新的CT数据，将其与原始CT数据进行对比，确定并筛选出模型重构过程中存在差异的点，在粗选压痛点中排除存在差异的点，对精选之后的压痛点重新编号为j，并确定精选后的水平单向应力σ_j及压痛点的应力集中系数K_j，其中，j＝1，2，3......M，M为精选压痛点的总数量。

4.根据权利要求3所述的基于奇异排除的滑膜炎压痛点具体位置定量分析方法，其特征在于：步骤S3确定过程如下：在S2分析的精选压痛点j的基础上，在前处理模块设置边界条件，施加均布拉力，执行求解，得到竖直单向应力ε_j以及竖直单向应力集中系数G_j，并确定其量化区域，其中，j＝1，2，3......M，M为精选压痛点的总数量。

5.根据权利要求4所述的基于奇异排除的滑膜炎压痛点具体位置定量分析方法，其特征在于：步骤S4确定过程如下：在S3求解的竖直单向应力ε_j以及应力集中系数G_j基础上，对水平方向应力比例系数α、竖直方向应力比例系数β进行求解，

其中，σ_j为精选后的水平单向应力，ε_j为竖直单向应力，j为编号，j＝1，2，3......M，M为精选压痛点的总数量，

根据式(1)、(2)求解以单向应力为主的每个压痛点综合应力λ_j，

其中，α为水平方向应力比例系数，β为竖直方向应力比例系数；K_j为压痛点的应力集中系数，G_j为竖直单向应力集中系数；δ为网格修正系数，γ为材质修正系数，对于两个修正系数，0-20岁取值0.97，20-50岁取值1.06，50岁以上取值1.23；μ为环境修正系数，春季取值1.09，夏季取值1.2，秋季取值1.03，冬季取值0.8，j为编号，j＝1，2，3......M，M为精选压痛点的总数量，

当压痛点综合应力λ_j为最大值时，取该压痛点的水平位置为滑膜炎治疗的最适穿刺水平位置。