[发明专利]一种快速高质量最大密度投影实现方法

权利要求书

1.一种快速高质量最大密度投影实现方法，包括以下步骤：

A.数据预处理，将数据分块整理存储并建立最大值最小值八叉树；

B.光线初始化，计算出光线的进入点并根据进入点位置将所有光线分配给其进入的数据块；

C.根据光线方向将数据块分组，采用并行计算分别处理每组数据块；

D.处理组内每个数据块的所有光线，当光线穿出数据块时，判断其是否已经穿出数据体；

若没有穿出数据体，根据采样位置将光线分配给下一个数据块后进入E步骤；

E.继续处理下一组数据块；

F.所有分组处理完成，

其特征在于：所述B步骤中的光线初始化，先利用移动模板的方法得到每条光线的进入深度后，再通过叠加的方法计算得到每条光线的进入点。

2.根据权利要求1所述的一种快速高质量最大密度投影实现方法，其特征在于：所述A步骤中的数据块大小为2ⁿ*2ⁿ*2ⁿ，n大于1，数据块小于CPU的缓存容量。

3.根据权利要求2所述的一种快速高质量最大密度投影实现方法，其特征在于：所述B步骤的所述移动模板大小为(2ⁿ+m)*(2ⁿ+m)*(2ⁿ+m)，其中m为大等于1的自然数，所述移动模板中心与每个待投影的数据块的中心重叠，所述移动模板记录有模板中光线的数目和模板数据，所述模板数据包括光线相对于模板中心的2D坐标偏移值和光线进入深度。

4.根据权利要求3所述的一种快速高质量最大密度投影实现方法，其特征在于，所述光线进入深度和进入点的计算方法为：

三维坐标表示为vector(x，y，z)，vector.x、vector.y、vector.z分别表示其X、Y、Z分量，

图像表示大小size(cx，cy)，size.cx、size.cy分别表示图像的宽和高，

体数据坐标系约定为Voxel坐标系，投影平面空间坐标系约定为View坐标系；体数据到投影平面空间的三维坐标变换为VoxeltoView，投影平面空间到体数据的三维坐标变换为ViewtoVoxel；

Voxel坐标vecVoxel与对应的View坐标vecView之间的转换关系为：

vecVoxel＝ViewtoVoxel×vecView

vecView＝VoxeltoView×vecVoxel

投影大小为imagesize，根据VoxeltoView计算数据块间在X、Y、Z三个方向的3D View坐标偏移DeltaBX、DeltaBY，DeltaBZ：

DeltaBX=VoxeltoView×Voxel(32,0,0)-VoxeltoView×Voxel(0,0,0)DeltaBY=VoxeltoView×Voxel(0,32,0)-VoxeltoView×Voxel(0,0,0)DeltaBZ=VoxeltoView×Voxel(0,0,32)-VoxeltoView×Voxel(0,0,0)]]>

当前数据块中心在投影空间中的坐标为CurPos，初始值为(0，0，0)，对于顺序为(bx，by，bz)的数据块而言，其对应的位置为：

Pos＝CurPos+DeltaBX×bx+DeltaBY×by+DeltaBZ×bz

对于该数据块内的第i个光线，其进入深度为：

Depth[i]＝Pos.z+EntryDepth[i]

其在投影面上对应的2D坐标索引值：

Index[i]＝Pos.x+Pos.y×imagesize.cx+PixelOffset[i]

以数据块为单位进行循环，计算出数据块内所有光线在的进入深度Depth[i]，根据其在投影面上的索引值Index[i]得到投影面上对应的进入深度Depth[index[i]]，如果Depth[i]＜Depth[index[i]]，则将更新Depth[index[i]]的值为Depth[i]；

得到投影面上所有光线的进入深度后，根据ViewtoVoxel求出光线在X、Y、Z三个方向3D Voxel坐标偏移DeltaX、DeltaY、DeltaZ：

DeltaX=ViewtoVoxel×View(1,0,0)-ViewtoVoxel×View(0,0,0)DeltaY=ViewtoVoxel×View(0,1,0)-ViewtoVoxel×View(0,0,0)DeltaZ=ViewtoVoxel×View(0,0,1)-ViewtoVoxel×View(0,0,0),]]>

同时求出Voxel(0，0，0)对应的View坐标Pos：

Pos＝ViewToVoxel×View(0，0，0)，

则对于投影平面上点(x，y)，令其进入深度为depth，其进入点为：

voxel＝Pos+x×DeltaX+y×DeltaY+depth×DeltaZ。