断层数据三维重构的子像素精度方法.pdf

断层数据三维重构的子像素精度方法 山东大学计算机科学与技术学院 杨兴强张彩明刘毅 2006-08-29 主要内容 • 等值面精度 • 等值面拓扑 Marching Cubes 方法 • 基本思想 – 像素值—顶点值， Cube – 线性插值 – 三角形网格 • 二义性问题（立方体内部） – 二维 – 三维 • 精度 – 尖点（增加Cube 内部点） • 效率 – 简化、八叉树方法 1. Marching Cubes中的误差 • 例子 • A、B 、C密度值分 别为700、550、200 • 等值线阈值为680 MC 求出的等值点AX=0.133 ，实际边界AY=1.2 原 因 • 对像素值的理解，MC方法将它看作一个 顶点的值。 • 直观上 v ∫∫E (x, y)dxdy S • 包含感兴趣物体的像素可能被划出等值 面 – 精度的损失 – 拓扑结构的不一致性 • 线性插值 像素数剧场和原始数据场 拓扑结构的不一致性 MC 线性插值 等值线阈值 2. 求边界点的新方法 m mB α(1A +) −α mC m m − B C α m m − A C 边界像素的确定 与MC 方法相比拓扑的变化 拓扑结构完全相同. 问题? 3. 解决整体拓扑的不一致性 假设1：T感兴趣物体内部的密度相对均匀，并且与背景 的密度值有显著的不同。 假设2：感兴趣物体的密度值总和是不变的，即使感兴 趣物体的密度被分布在相邻的几个像素内。 边界内像素： 边界外像素： 过渡像素： 边界像素： 算法1：求过渡像素 (1)求出所有内部像素，记为Sin ，周围的像素记为S0 . i=0; (2) 若 q ∈S ，S ={p | D(p )D(q), q ∈S ,p 与 q相邻} i i+1 i (3) If Si+1≠Φ then let i i+1; repeat step (2); Else step (4). (4) S =S +S +…+S , S 是过渡像素集合. t 1 2 i t 算法2：求边界像素 (1) Suppose work = D , which express the excessive density in a the transition pixels, letj =1. (2) compute D =|S | · (T −D ), where T is the density of the j j b interested object, the meaning of S is same as algorithm 1.