[发明专利]自动采样厚度检测的非线性口腔CT全景图像合成方法

权利要求书

1.自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于，步骤包括有：

步骤一、输入口腔CT重建后体数据，得到体数据最大密度投影图像；

步骤二、对步骤一得到的体数据最大密度投影图像进行二维高斯滤波，得到滤波体数据最大密度投影图像；

步骤三、计算步骤二得到的滤波体数据最大密度投影图像的直方图，得到滤波体数据最大密度投影图像分割阈值；

步骤四、用步骤三得到的滤波体数据最大密度投影图像分割阈值对步骤二得到的滤波体数据最大密度投影图像进行分割，得到上下颌中牙齿在序列图像中范围；

步骤五、由步骤四得到的上下颌中牙齿在序列图像中范围生成得到横断面最大密度投影图像；

步骤六、统计步骤五得到的横断面最大密度投影图像的直方图，得到横断面最大密度投影图像分割阈值；

步骤七、根据步骤六得到的横断面最大密度投影图像分割阈值对步骤五得到的横断面最大密度投影图像进行分割，得到横断面最大密度投影图像分割结果；

步骤八、对步骤七得到的横断面最大密度投影图像分割结果进行轮廓统计，然后滤波，得到牙弓轮廓；

步骤九、对步骤八得到的牙弓轮廓进行细化操作，得到初步牙弓曲线；

步骤十、统计步骤九得到的初步牙弓曲线的端点，得到两端点初步牙弓曲线；

步骤十一、对步骤十得到的两端点初步牙弓曲线进行控制点采样操作，得到初步牙弓曲线进控制点；

步骤十二、根据步骤十一得到的初步牙弓曲线进控制点构建B-Cubic样条曲线，得到正式牙弓曲线；

步骤十三、根据步骤八得到的牙弓轮廓和步骤十二得到的正式牙弓曲线计算，得到虚拟全景图像采样厚度；

步骤十四、对步骤十二得到的正式牙弓曲线进行等像素间隔采样，并计算采样所在点的法向量，并根据步骤十三得到的虚拟全景图像采样厚度对步骤一的口腔CT重建后体数据进行切割，得到全景容积图像和虚拟全景容积数据；

步骤十五、沿步骤十四的正式牙弓曲线的法向量的方向对步骤十四的虚拟全景容积数据和全景容积图像进行合成，得到虚拟全景图像；

所述步骤十三具体为，

步骤13.1、对步骤十二得到的正式牙弓曲线进行等间距采集，得到等间距采样点，进入步骤13.2；

步骤13.2、在步骤13.1的等间距采样点处，计算正式牙弓曲线的法向量，进入步骤13.3；

步骤13.3、沿步骤13.3的法向量方向，计算步骤八得到的牙弓轮廓边缘的距离，得到牙弓轮廓边缘距离，进入步骤13.4；

步骤13.4、统计步骤13.3得到牙弓轮廓边缘距离的直方图，得到距离直方图，进入步骤13.5；

步骤13.5、计算步骤13.4得到的距离直方图的最大峰值，并取最大峰值的±20范围的平均值，得到虚拟全景图像采样厚度；

所述步骤十四具体为，

步骤14.1、对步骤十二得到的正式牙弓曲线进行等像素间隔采样，得到等像素间隔采样点，进入步骤14.2；

步骤14.2、并计算等像素间隔采样点所在位置的法向量，沿法向量方向，并根据步骤十三得到的虚拟全景图像采样厚度对步骤一的口腔CT重建后体数据进行切割，进入步骤14.3；

步骤14.3、根据步骤14.1得到的等像素间隔采样点和步骤14.2得到的切割图，依次得到所有等像素间隔采样点的进入切割，每个等像素间隔采样点各得到全景容积图像，进入步骤14.4；

步骤14.4，将所有等像素间隔采样点的全景容积图像依次排列，得到虚拟全景容积数据。

2.根据权利要求1所述的自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于：所述体数据最大密度投影图像为体数据冠状面最大密度投影图像或者体数据矢状面的最大密度投影图像。

3.根据权利要求2所述的自动采样厚度检测的非线性口腔CT全景图像合成方法,其特征在于：所述步骤三具体为，

步骤3.1、计算步骤二得到的滤波体数据最大密度投影图像的直方图，进入步骤3.2；

步骤3.2、计算直方图的峰值，进入步骤3.3；

步骤3.3、从大至小排列直方图的灰度值得到第一个直方图峰值，进入步骤3.4；

步骤3.4、采用高斯函数逼近第一个直方图峰值，取右侧置信区间98％的值，得到滤波体数据最大密度投影图像分割阈值。