求教:三维实体的断层扫描后任意方向截面的平面图形重建问题

[FmzHxj]

问题描述如下:
小弟有一组人体某器官CT图(50张左右,已经把DICOM3转化成BMP),根据每个CT图上器官的轮廓线,小弟已经用蒙皮算法将该器官进行了三维重建.在重建的三维模型里,客户提出需要任意方向截面的纹理,该纹理数据则从50张CT图上的像素数据中通过插值计算而获得.
小弟经过几天思考,设计如下方案:
先将50张CT图中的像素保存到一个矩阵(数组)中,在客户旋转三维模型到某个需要的角度时,获得视点(eye)和up矢量,然后再记录用户截面选取的直线段.
用所获得视点(eye)和up矢量与初始值的相对差异来旋转像素矩阵(数组),然后在矩阵上找到直线段所在的平面的所有像素.
这个方案是不可行的.
第一,目前计算机的情况不允许有这么大的矩阵(数组)存在(图片的分辨率为1280 * 1280,50张啊,可想而知了),但可以动态地按每张CT图与硬盘进行导入导出.
第二,每个图片之间虽然间隔是均匀的,但像素不是连续的,间隔像素有几甚至十几之多,这样导致插值时需要有一个类似中点线之类的算法进行像素衍生,涉及到不止一张这样的CT图,又需要频繁地读硬盘数据,使得在速度上出现严重问题.而且插值算法过于复杂.
希望哪位高手或者有过类似项目设计经验的同志指点迷津,给点提示,谢谢

[asett1]

你用透视算法，再投影到任意平面就ok了

[FmzHxj]

透视？跟这个没什么关系吧，所有投影都是平行投影。

[孤必有邻]

你的问题比较专业，是不是类似这样一个问题：
把一段木头沿某一方向均匀切成小段，记录各段的平面纹理，
然后依据该组纹理图，推算原木在任意切面上的纹理？

[FmzHxj]

楼上所言正是！

[FmzHxj]

小弟的方案目前作了改进，在视觉允许的限度内把图片缩小到了512 * 512，似乎有成功的可能，但是作了很大的牺牲。

[crazy_lazy_pig]

学习一下matlab吧,我曾经用它处理512 * 512 的二值图100张(是血管切片图),做其差值运算也就是一个小时而已,挺管用的.