基于二维X射线图像序列滤波反投影的分块三维重建方法
专利汇可以提供基于二维X射线图像序列滤波反投影的分块三维重建方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种基于二维 X射线 图像序列 滤波反投影 的分 块 三维重建 方法，包含如下步骤：步骤一，根据初始分块数目，建立立方体分块；步骤二，对各个立方体分块进行局部重建获得分块重建结果；步骤三，组合各个分块的重建结果；步骤四，根据分块的重建结果计算所有分块数目集合，根据步骤一至步骤三进行重建并反馈其时间复杂度，直至分块数目集合中各元素均计算完毕，获得反馈的时间复杂度集合；从反馈的时间复杂度集合中选取具有最小时间复杂度的分块数目，将其作为最优分块数目，重新执行步骤一至步骤三，获得最优分块重建结果，完成三维重建。本发明能够使多线程编程技术在三维重建领域得到充分的运用，从而极大提高了重建速度。，下面是基于二维X射线图像序列滤波反投影的分块三维重建方法专利的具体信息内容。
1.一种基于二维X射线图像序列滤波反投影的分块三维重建方法，其特征在于，包含如下步骤：
步骤一，根据初始分块数目，建立立方体分块：包括建立离散坐标系和确定各分块坐标区；
步骤二，对各个立方体分块进行局部重建获得分块重建结果：包括待重建数据规范化、待重建数据预处理以及基于二维X射线图像的分块反投影重建；
步骤三，组合各个分块的重建结果；
步骤四，根据分块的重建结果计算所有分块数目集合，根据步骤一至步骤三进行重建并反馈其时间复杂度，直至分块数目集合中各元素均计算完毕，获得反馈的时间复杂度集合；
从反馈的时间复杂度集合中选取具有最小时间复杂度的分块数目，将其作为最优分块数目，重新执行步骤一至步骤三，获得最优分块重建结果，完成三维重建。
2.根据权利要求1所述的一种基于二维X射线图像序列滤波反投影的分块三维重建方法，其特征在于，步骤一包括如下步骤：
初始化分块数目N0＝4；
建立立方体分块结构：先构建一个重建立方体，以重建立方体其中一个顶点为坐标原点，与原点相连的三边分别作为x、y、z三个坐标轴建立右手坐标系，用坐标点(x，y，z)表示任一分块中的一个体素，则任意体素坐标点(x，y，z)满足：
，其中参数
XDIM、YDIM、ZDIM分别表示x、y、z轴的最大值，表示自然数集合；
确定各分块坐标区：保持y、z轴不变，对x轴进行分割，通过如下公式调整参数XDIM使其为分块数目N0的整数倍：
把x轴[0，XDIM]的范围平均分割为N0个坐标范围，即：
N0个x轴坐标范围分别结合y、z轴构成N0个立方体分块。
3.根据权利要求2所述的一种基于二维X射线图像序列滤波反投影的分块三维重建方法，其特征在于，步骤二中基于二维X射线图像的分块反投影重建包括如下步骤：
建立用于重建的二维X射线图像序列中每幅图像的像素和重建立方体中体素坐标的关系：图像序号用proj_num表示，θ表示相邻图像之间的拍摄夹角，依次把每幅图像数据存入一维序列，体素坐标为(x，y，z)，当前要进行反投影的二维X射线图像的拍摄角度为φ＝proj_num·θ，所述体素在序号为proj_num的二维X射线图像中对应的坐标(x′p，y′p)为：
其中，D为成像装置中X射线发射源到探测器的距离，d为成像装置中X射线发射源到待重建物体的几何中心的距离；
令P(i′，j′)为一幅二维X射线图像中任意坐标点(i′，j′)的像素值，用差值法获得对应的经滤波后的修正像素值
用三个维度的相对位置(i，j，k)表示一个立方体数据中的体素，以立方体数据的几何中心为原点，分别以i、j、k三个维度为x，y，z轴建立右手直角坐标系；
对 于 每 一 个 基 于 i轴 分 割 开 来 的 分 块 区 域，其i 轴 范 围 分 别 为结合j、k轴表示为N0
个长方体的形式，以三元组的形式表示为：
其中V表示体素值，l为分块序号；
对于每一个分块内体素(i，j，k)，二维X射线图像序列在分块Vl[i][j][k]上的反投影如下式所示：
其中 表示在序号为proj_num的二维X射线图像中第i′行第j′个
像素经滤波后的修正像素值，NUM表示二维X射线图像的最大序号，二维X射线图像总数为NUM+1，对于分块中的每个体素都用进行反投影运算，获得最后的分块重建结果
4.根据权利要求3所述的一种基于二维X射线图像序列滤波反投影的分块三维重建方法，其特征在于，步骤三中组合分块重建结果步骤如下：
利用步骤二中获得的分块重建结果进行拼装组合，描述为：
其中 当(i，j，k)遍历立
方体中所有的体素时，V[i][j][k]即为重建结果。
5.根据权利要求1所述的一种基于二维X射线图像序列滤波反投影的分块三维重建方法，其特征在于，步骤四中根据反馈的时间复杂度集合计算最优分块数目包括如下步骤：
步骤41，定义综合贡献因子：分块数目N0由x轴最大值XDIM约束，令参数α＝XDIM mod N0；
定义综合贡献因子η：
λ为表征精度和性能的参数，λ∈[0，1]；t0和t分别表示分块数目N0取初始值4时和取当前的数值时的重建时间；
步骤42，确定分块数目：先选取所有满足参数α＝0的分块，并统计其数目为m，生成集合 对应的重建时间 记为 然后从集合 中选取重建
时间最少且分块数目也最少的分块数记为 对应重建时间为tk，定义最优分块数目集合
计算最优分块数目集合 中各个分块数目所对应的综合贡献因子ηi：
其中， 和 分别表示分块数目为(i+Nk)和Nk时的重建时间，
求出最优分块数目集合 中每个分块数对应的综合贡献因子ηi后，取该值最小且分块数目最少的分块数为最优分块数目N0*。
基于二维X射线图像序列滤波反投影的分块三维重建方法
技术领域
背景技术
发明内容
贡献因子ηi：
XDIM、YDIM、ZDIM分别表示x、y、z轴的最大值，表示自然数集合。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 锥束CT圆加直线轨迹反投影滤波重建方法 | 2020-05-11 | 95 |
| 针对有限角度断层摄影中的滤波反投影的图像重建的方法、计算设备、装置和可读介质 | 2020-05-12 | 857 |
| C型臂半精确滤波反投影断层成像方法 | 2020-05-11 | 683 |
| 针对有限角度断层摄影中的滤波反投影的图像重建的方法 | 2020-05-11 | 246 |
| 针对半覆盖螺旋锥束CT的单层重排滤波反投影重建方法 | 2020-05-11 | 467 |
| 用于两个倾斜的圆形的计算机断层成像重建 | 2020-05-12 | 316 |
| 多源螺旋CT反投影滤波精确重建系统 | 2020-05-11 | 334 |
| 一种基于C形臂的锥束CT三维数字成像方法 | 2020-05-12 | 456 |
| 锥束CT圆加直线轨迹反投影滤波重建方法 | 2020-05-11 | 166 |
| 一种由聚焦堆栈重建光场的滤波反投影方法和装置 | 2020-05-12 | 830 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
