用动态数字模型检测CT血管造影术重建算法性能的方法技术

一种生物医学成像技术领域的用动态数字模型检测ＣＴ血管造影术重建算法性能的方法。本发明专利技术采用包含运动的造影剂团注和简单的几何结构构成的三维动态数字模型，利用计算机按实际设备中的参数设置对模型进行仿真数据采集，再用得到的仿真数据按算法重建，最后比较重建图像与模型之间的差异，以检测算法的性能。所述三维动态数字仿真模型，包括模拟造影剂团注的动态椭球，其各轴长固定，椭球中心坐标是一个随时间变化的函数，表示为Ｘ＝ｘ（ｔ），Ｙ＝ｙ（ｔ），Ｚ＝ｚ（ｔ）；以及由４个半圆柱和２个圆柱组成的柱状模型。本发明专利技术便于操作和重复定量执行地检测ＣＴ血管造影术重建算法，以判断不同重建算法的优劣或同一种算法在不同的参数条件下的性能差异，进而找出最佳参数设置。

【技术实现步骤摘要】

本专利技术涉及的是一种生物医学成像
的方法，具体地说，是一种用动态数字模型检测CT血管造影术重建算法性能的方法。
技术介绍
CT血管造影术的主要原理是在目标周围静脉内以一定速度注入团注造影剂，经时间延迟，至靶血管内造影剂充盈至高峰期，用CT设备对其进行高速、连续的容积数据采集，然后通过适当的图像重建算法做后处理完成二维及三维图像的重建。最终重建图像的质量受到造影剂的剂量，注射速度，延迟时间，CT设备和重建算法这多方面的影响。经对现有技术的文献检索发现，M.Piotin等在CT angiography，MRangiography and rotational digital subtraction angiography forvolumetric assessment of intracranial aneurysms.An experimental study(CT血管造影术，MR血管造影术，数字减影血管造影术对颅内动脉瘤体成像的初步研究)，Neuroradioragy(神经放射学)，2003，Vol.45，No.6，pp404-409中通过对实物血管模型进行实验，从而对不同的血管造影术重建算法的性能进行了比较。采用实物血管模型的不足指出在于实物模型架设复杂，且对硬件设备要求高；对实物模型的清洁和维护麻烦；每次实验的初始状态和参数不易控制，存在人为误差，最终影响算法性能指标间的可比性。
技术实现思路
本专利技术的目的在于针对现有技术的不足，提供一种用动态数字模型检测CT血管造影术重建算法性能的方法。本专利技术采用三维动态数字模型，便于操作和重复定量执行地检测CT血管造影术重建算法，以判断不同重建算法的优劣或同一种算法在不同的参数条件下的性能差异，进而找出最佳参数设置。本专利技术是通过以下技术方案实现的，本专利技术采用包含运动的造影剂团注和简单的几何结构构成的三维动态数字模型，利用计算机按实际设备中的参数设置对模型进行仿真数据采集，再用得到的仿真数据按算法重建，最后比较重建图像与模型之间的差异，以检测算法的性能。本专利技术包括以下步骤(1)用椭球和半圆柱体组合构成三维动态数字仿真模型。所述三维动态数字仿真模型，包括模拟造影剂团注的动态椭球，其各轴长固定，椭球中心坐标是一个随时间变化的函数，可表示为X＝x(t)，Y＝y(t)，Z＝z(t)；以及由4个半圆柱和2个圆柱组成的柱状模型。所述的椭球、半圆柱体和圆柱体组合构成三维动态数字仿真模型，其具体过程如下设该数字模型整体为一长方体结构，定义其三维直角坐标系的原点位于该长方体的中心，X轴、Y轴、Z轴分别延长方体各棱，并以LX、LY、LZ分别代表长方体X轴、Y轴和Z轴的各棱长。以下所表示的各几何结构的参数均以该坐标系内归一化数值描述。所述归一化数值a是指，a=PL,]]>其中P指是指坐标系中某一分量的实际坐标值，L指该分量P所属的坐标轴的半棱长。●定义空间椭球方程为x2A2+y2B2+z2C2≤1,]]>一般情况下，椭球中心坐标位于(EX，EY，EZ)，A、B、C分别代表椭球各轴长。●定义空间半圆柱方程为x2+y2≤R2|z|≤ZLim0≤y≤R]]>或x2+y2≤R2|z|≤ZLim-R≤y≤0,]]>一般情况下，半圆中心坐标位于(CX，CY，CZ)，R代表圆的半径长。●定义空间圆柱方程为x2+y2≤R2|z|≤ZLim,]]>一般情况下，半圆中心坐标位于(CX，CY，CZ)，R代表圆的半径长。●以方程x2+y2≤R2|z|≤ZLim0≤y≤R]]>建立上半圆柱，半圆柱中心坐标为(0，0，0)；半径R＝0.8；ZLim＝1.0；密度d＝d1。所述的密度d代表各组织的相对平均密度结构，不同的密度对应的不同组织在模型中将反映为不同的灰度级，这里定义人骨的密度为2.00，其它各组织密度参考此密度给出。●以方程x2+y2≤R2|z|≤ZLim-R≤y≤0]]>建立下半圆柱，半圆柱中心坐标为(0，0，0)；半径R＝0.8；ZLim＝1.0；密度d＝d2。●以方程x2+y2≤R2|z|≤ZLim0≤y≤R]]>建立上半圆柱，半圆柱中心坐标为(0.4，0，0)；半径R＝0.4；ZLim＝1.0；密度d＝d3。●以方程x2+y2≤R2|z|≤ZLim-R≤y≤0]]>建立下半圆柱，半圆柱中心坐标为(-0.4，0，0)；半径R＝0.4；ZLim＝1.0；密度d＝d4。●以圆柱体建立血管一，圆柱体的中心坐标为(0.4，0，0)；半径R＝0.1；ZLim＝1.0；密度d＝d5。●以圆柱体建立血管二，圆柱体的中心坐标为(-0.4，0，0)；半径R＝0.1；ZLim＝1.0；密度d＝d6。●以运动的椭球体构成造影剂团注一，其各轴长A＝0.1，B＝0.1，C＝0.3；椭球体的中心坐标为与时间有关的函数，其中X＝0.4，Y＝0，Z＝z(t)；密度d＝d7。●以运动的椭球体构成造影剂团注二，其各轴长A＝0.1，B＝0.1，C＝0.3；椭球体的中心坐标为与时间有关的函数，其中X＝-0.4，Y＝0，Z＝z(t)；密度d＝d8。(2)根据需要检测的CT血管造影术重建算法所要求的数据采集方式(如基于团注追踪的多层螺旋CT扫描方式)，根据实际X射线采集原理，利用公式p＝∫ρ·dl(其中，dl为空间内某一X射线经过的路径，ρ为该路径上所经过的物体密度，p则为采得的数据)，用数字方法计算X射线路经，并用模型密度替换实际物体密度，对公式离散化处理后按扫描模式获得三维动态数字仿真模型在不同时刻下的仿真投影数据。所述的仿真投影是根据实际投影数据采集原理，用解析式或数值积分的方法获取仿真投影数据的过程。(3)用需要检测的CT血管造影术重建算法，对第(2)步中获得的仿真投影数据进行图像重建。(4)利用第(3)步中获得的重建结果与模型原始图像进行比较分析，对需要检测的CT血管造影术重建算法的性能指标给出评定，如用重建图像与模型之间的标准差以检测重建图像的空间分辨率，标准差越小说明图像质量越高。具体所选用的性能指标和其数值随重建算法不同和实际需要而变化。(5)改变模型中的用于描述模型的参数(如造影剂密度和团注运动速度)和重建算法中的参数，重复(1)-(4)步，进一步检测重建算法在不同参数条件下性能变化情况。本专利技术用动态数字模型检测CT血管造影术重建算法性能的方法的原理在于用动态的虚拟数字模型模拟造影剂团注在血管中扩散的过程，并通过计算机模拟实际数据采集过程，获取仿真投影，并对仿真投影数据进行重建，最后通过重建图像与数字模型间的比较和分析，给出对重建算法优劣和适用条件的检测结果。本专利技术的有益效果是(1)采用数字模型检测算法性能，可以放弃实物模型，从而降低成本。(2)数字模型的初始状态和参数设置不会引本文档来自技高网...

【技术保护点】
一种用动态数字模型检测ＣＴ血管造影术重建算法性能的方法，其特征在于：包括以下步骤：　　　　（１）用椭球和半圆柱体组合构成三维动态数字仿真模型，所述三维动态数字仿真模型，包括模拟造影剂团注的动态椭球，其各轴长固定，椭球中心坐标是一个随时间变化的函数，表示为Ｘ＝ｘ（ｔ），Ｙ＝ｙ（ｔ），Ｚ＝ｚ（ｔ）；以及由４个半圆柱和２个圆柱组成的柱状模型；　　　　（２）根据需要检测的ＣＴ血管造影术重建算法所要求的数据采集方式，根据实际Ｘ射线采集原理，利用公式ｐ＝∫ρ.ｄｌ，其中，ｄｌ为空间内某一Ｘ射线经过的路径，ρ为该路径上所经过的物体密度，ｐ则为采得的数据，用数字方法计算Ｘ射线路经，并用模型密度替换实际物体密度，对公式离散化处理后按扫描模式获得三维动态数字仿真模型在不同时刻下的仿真投影数据；　　　　（３）用需要检测的ＣＴ血管造影术重建算法，对第（２）步中获得的仿真投影数据进行图像重建；　　　　（４）利用第（３）步中获得的重建结果与模型原始图像进行比较分析，对需要检测的ＣＴ血管造影术重建算法的性能指标给出评定，具体所选用的性能指标和其数值随重建算法而变化；　　　　（５）改变模型中的用于描述模型的参数和重建算法中的参数，重复（１）－（４）步，进一步检测重建算法在不同参数条件下性能变化情况。...

【技术特征摘要】

【专利技术属性】
技术研发人员：赵俊，刘尊钢，庄天戈，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：31[中国|上海]