首页 | 本学科首页   官方微博 | 高级检索  
     

医学图像三维重建中的快速表面绘制
引用本文:李飞鹏, 廖孟扬, 徐振勤, 王思贤. 医学图像三维重建中的快速表面绘制[J]. 武汉大学学报 ( 信息科学版), 2000, 25(2): 153-157.
作者姓名:李飞鹏  廖孟扬  徐振勤  王思贤
作者单位:1 武汉测绘科技大学测绘遥感信息工程国家重点实验室, 武汉市珞喻路129号, 430079;;2 武汉大学图像及信息研究所, 武汉市珞珈山, 430072
摘    要:简述了医学图像三维可视化表面绘制的基本原理。由于三维体数据中的等值面直接对应着组织器官的表面,而自然器官的表面应处处连续,为进一步加快绘制速度,假定待提取等值面连续,对三维体数据先进行 8像素等间隔采样和阈值判断,分离出边界点,然后只对边界点邻域做阈值判断。与DividingCubes算法相比,优化后的算法计算量不到原来的 1/4,绘制结果图像也比较逼真。

关 键 词:三维可视化  表面绘制  图像分割  等值面
收稿时间:1999-09-09

Fast Surface Rendering for Three Dimensional Medical Image Reconstruction
LI Feipeng, LIAO Mengyang, XU Zhenqin, WANG Sixian. Fast Surface Rendering for Three Dimensional Medical Image Reconstruction[J]. Geomatics and Information Science of Wuhan University, 2000, 25(2): 153-157.
Authors:LI Feipeng  LIAO Mengyang  XU Zhenqin  WANG Sixian
Affiliation:1 National Laboratory for Information Engineering in Surveying, Mapping and Remote Sensing, WTUSM, 129 Luoyu Road, Wuhan, China, 430079;;2 Image Processing Lab, Wuhan University, Luojia Hill, Wuhan, China, 430072
Abstract:Surface rendering, together with volume rendering, is the most popular sort of algorithm used in scientific visualization. Though developed later than surface rendering and having a wider adaptability to object data sets, the algorithm of volume rendering still can not be put into use in practical medical imaging application,simply because it could not meet the speed requirement for interactive operation due to the large amount of computation involved.Surface rendering thus becomes the general method for visualization in most medical imaging, for it is at least more computation-economic than its counterpart. The common algorithm of surface rendering can be generalized by two main steps - isosurface extracting and isosurface shading. Of all methods proposed for isosurface extracting, Marching Cubes (or Dividing Cubes) is the most widely used one. To pick up an isosurface for display, it traverses all the cells in the data set, and for each cell it compares the included eight vertexes to a certain threshold to decide whether it is a boundary cell or not. However, extensive investigation shows that, of all the units checked out, there are usually only less than 10% are useful. Considering the inefficiency in the search for boundary cells the whole run time can be further reduced by applying an optimized search order. Because each isosurface we are extracting from the volume data set actually corresponds to the epidermis of a human organ, we assume that the isosurfaces are continual as well. With the continuity assumption, the isosurface extracting can be accomplished by two steps. First, the volume data set are subsampled by 1/8 voxel.After comparing with the threshold, all the voxels in the subset are classified into boundary points and non-boundary points. After the search for isosurface is carried out, but only in the neighboring space of boundary points, those voxels labeled as non-boundary points will be completely neglected for further processing. As a result, the optimized algorithm is at least 4 times more efficient, with the possibility that the rendered image be mildly distorted by an average error of less than 1 voxel in all directions.When applied to practical use - usually large image series, it is foreseeable that the distortion will be so slight as nearly undetectable. Since surface fitting alone can not produce any image usable, the segmentation of volume data has to be discussed. For many CT image series, the separation of different tissues such as bones, teeth from brain could be accomplished conveniently by setting an apposite threshold, but so far as the MRI and B scan image series are concerned, other more elaborate methods should be applied. This paper discusses the improvement of surface rendering for practical use. In section 1, the basic principle and rendering scheme are presented. Section 2 deals with the volume data segmentation. Then in section 3, an optimized isosurface searching order is proposed. Section 4 discusses isosurface shading. In section 5 the rendering speed of the new method is provided with generated image in contrast to the original one.
Keywords:D visualization  surface rendering  image segmentation  isosurface
本文献已被 CNKI 万方数据 等数据库收录!
点击此处可从《武汉大学学报(信息科学版)》浏览原始摘要信息
点击此处可从《武汉大学学报(信息科学版)》下载全文
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号