[发明专利]一种基于流体计算的非先验性骨折识别方法

权利要求书

1.一种基于流体计算的非先验性骨折识别方法，其特征在于：

S1、将待检测骨折CT数据集的几何外边界构成的六面体转换为流体计算域；将CT各像素点转换为对应的计算域空间网格节点；

S2、将CT各像素点的HU值信息，转换成对应网格单元的物性参数，并识别皮质骨区域；

S3、在计算域外边界上选择一区域作为流体入口；给定流体入口的气体参数；

S4、基于上述条件，使用计算流体力学技术求解，直到一定步数；

S5、基于流体计算的可视化计算结果，判定是否存在皮质骨破损。

2.根据权利要求1所述的一种基于流体计算的非先验性骨折识别方法，其特征在于，所述步骤S1中，CT数据点集的空间几何外边界构成了三维长方体，此长方体为空间计算域；规则排列的每八个数据点构成一个微长方体，此微长方体为流体计算单元；流体计算单元的长宽高对应CT数据集的三方向空间单位间隔；在进行下述流体力学计算时直接采用此空间计算域和空间计算单元，无须另行划分网格。

3.根据权利要求1所述的一种基于流体计算的非先验性骨折识别方法，其特征在于，所述步骤S2中，将全部流体计算单元的对应CT数据集的微长方体的HU值转换为物性数据；HU值为皮质骨的空间计算单元标记为壁面单元；松质骨和软组织（例如：皮肤、肌肉、肌腱、血管、体液、神经和内脏）对应的空间计算单元标记为空气。

4.根据权利要求1所述的一种基于流体计算的非先验性骨折识别方法，其特征在于，所述步骤S3中，选择空间计算域外部边界的一个长方形区域作为流体进口；该区域计算单元的物性必须为空气；流动方向朝向空间计算域内侧。

5.根据权利要求1所述的一种基于流体计算的非先验性骨折识别方法，其特征在于，所述步骤S4中，选择使用非定常三维无粘可压缩流体计算方法，在上述计算域上求解非定常三维无粘可压缩方程，直到计算收敛。

6.根据权利要求1所述的一种基于流体计算的非先验性骨折识别方法，其特征在于，所述步骤S5中，基于流体计算的可视化计算结果，在连续的皮质骨内侧区域查看气体压力变化，若气体压力变化为零，则不存在皮质骨破损；若皮质骨内侧区域压力变化不为零，则皮质骨区域存在骨折；观察皮质骨区域的流动速度不为零的区域，该区域即可识别为骨折区域。