基于标签化医学图像的多阈值三维重建算法
|
摘要
针对传统面绘制重建方法 MC无法一次性提取多个阈值器官的问题,提出一种将MRI图像中靶区及其周边组织进行标签化分割的多阈值三维重建算法。该算法通过将二维图像中的多阈值器官标签化为简单的整数,降低提取等值面时的数据存储量,进而提高等值面的绘制速度。同时,定义了多阈值三维重建时体素顶点索引方式和等值面相交形态,减少了传统MC算法提取多个阈值器官时存在三角面片与顶点复用的情况,且仅需一次性扫描即可重建多个器官等值面。实验结果表明,本文算法较传统MC算法,重建器官数量越多,器官彼此结构越紧密,三角面片与顶点复用情况减少越明显,在保证重建效果的同时,绘制速度可提高30%。
As a traditional surface reconstruction algorithm the marching cubes(MC)is unable to extract multi-threshold organs at one time.We propose a multi-threshold reconstruction algorithm to divide the target area and surrounding tissues in the MRI image into several label values.Changing the lables of multi-threshold organs into brief numbers can decrease the memory space and accelerate the speed of drawing contour surfaces.We also define the index way of voxel vertexes and isosurface intersecting morphology of multi-threshold three-dimensional reconstruction,which can avoid the situation of triangle facets and vertexes reuse of the MC.Furthermore,the contour surface reconstruction of the multi-organ requires only one time scanning of the MRI image.Experimental results show that the number of facets and vertices is decreased more significantly when the number of organs is larger and their relationship is closer in comparison with the traditional MC algorithm,and the reconstruction speed is improved by 30% while the reconstruction effect is guaranteed.
As a traditional surface reconstruction algorithm the marching cubes(MC)is unable to extract multi-threshold organs at one time.We propose a multi-threshold reconstruction algorithm to divide the target area and surrounding tissues in the MRI image into several label values.Changing the lables of multi-threshold organs into brief numbers can decrease the memory space and accelerate the speed of drawing contour surfaces.We also define the index way of voxel vertexes and isosurface intersecting morphology of multi-threshold three-dimensional reconstruction,which can avoid the situation of triangle facets and vertexes reuse of the MC.Furthermore,the contour surface reconstruction of the multi-organ requires only one time scanning of the MRI image.Experimental results show that the number of facets and vertices is decreased more significantly when the number of organs is larger and their relationship is closer in comparison with the traditional MC algorithm,and the reconstruction speed is improved by 30% while the reconstruction effect is guaranteed.
引文
[1]Wang Zhi-yuan,Li Shu-xiang.A 3Dbrachytherapy treatment planning system[J].Chinese Journal of Medical Instrumentation,2000,21(1):7-10.(in Chinese)
[2]Sutherland J G H,Furutani K M,Thomson R M.A Monte Carlo investigation of lung brachytherapy treatment planning[J].Physics in Medicine and Biology,2013,58(14):4763.
[3]Chen W C,Katcher J.Radioactive seed migration after transperineal interstitial prostate brachytherapy and associated development of small-cell lung cancer[J],Brachytherapy,2012,11(5):354-358.
[4]Lorensen W E,Cline H E.Marching cubes:A high resolution3Dsurface construction algorithm[J].ACM Siggraph Computer Graphics,1987,21(4):163-169.
[5]Zhang Y,Qian J.Dual contouring for domains with topology ambiguity[J].Computer Methods in Applied Mechanics and Engineering,2012,217:34-45.
[6]Jia T,Feng X,Guo X.Improvement of marching cubes algorithm based on sign determination[C]∥Proc of 2010 2nd IEEE International Workshop on Education Technology and Computer Science(ETCS),2010:194-197.
[7]Raman S,Wenger R.Quality isosurface mesh generation using an extended marching cubes lookup table[J].Computer Graphics Forum,2010,27(3):791-798.
[8]Jiang H,Zhang R,Liu H,et al.An improved 3Dreconstruction method[C]∥Proc of IEEE International Conference on Multimedia Information Networking and Security(MINES’09),2009:209-212.
[9]Liu Y,Xing H L.A boundary focused quadrilateral mesh generation algorithm for multi-material structures[J].Journal of Computational Physics,2013,232(1):516-528.
[10]Lederman C,Joshi A,Dinov I,et al.The generation of tetrahedral mesh models for neuroanatomical MRI[J].NeuroImage,2011,55(1):153-164.
[11]Tian S,Liu G,Weng Z,et al.Isosurface generation algorithm based on spatial discrete points with attribute[J].Journal of Convergence Information Technology,2012,7(4):86-96.
[12]Li Zhong-hua,Zhou Fu-gen,Bai Xiang-zhi.A GPU-based fast volume CT reconstructive algorithm method[J].Journal of Biomedical Engineering,2011,28(2):238-242.(in Chinese)
[13]Li Jia,Long Peng-cheng,Luo Yue-tong,et al.Research and application of a fast 3D-reconstruction method in accurate radiotherapy treatment planning system[J].Nuclear Techniques,2010,33(1):69-74.(in Chinese)
[14]Chen Chun-xiao,Zhang Juan.3Dreconstruction of multiple views based on trifocal tensor[J].Journal of Biomedical Engineering,2012,29(4):769-774.(in Chinese)
[15]Lehmann G.Label objects representation and manipulation with ITK[J].The Insight Journal,2008:1-34.
[16]Xu Xiao-ling,Li Xian-min,Li Gui-qing,et al.A fast algorithm of marching cubes in volume reconstruction[J].Journal of System Simulation,2002,14(4):509-513.(in Chinese)
[17]Vollmer J,Mencl R,Mueller H.Improved laplacian smoothing of noisy surface meshes[J].Computer Graphics Forum,1999,18(3):131-138.
[18]Wang H,Chen H,Su Z,et al.Versatile surface detail editing via Laplacian coordinates[J].The Visual Computer,2011,27(5):401-411.
[1]王志远,李树祥.三维近距离放射治疗计划系统的研究[J].医疗卫生装备,2000,21(1):7-10.
[12]李忠华,周付根,白相志.一种基于GPU的体积CT快速重建算法[J].生物医学工程学杂志,2011,28(2):238-242.
[13]李佳,龙鹏程,罗月童,等.精确放射治疗计划系统中快速三维重建方法的研究与应用[J].核技术,2010,33(1):69-74.
[14]陈春晓,张娟.基于三焦点张量的多视图三维重构[J].生物医学工程学杂志,2012,29(4):769-774.
[16]徐晓玲,李现民,李桂清,等.体素重建中的快速移动立方体方法[J].系统仿真学报,2002,14(4):509-513.
[2]Sutherland J G H,Furutani K M,Thomson R M.A Monte Carlo investigation of lung brachytherapy treatment planning[J].Physics in Medicine and Biology,2013,58(14):4763.
[3]Chen W C,Katcher J.Radioactive seed migration after transperineal interstitial prostate brachytherapy and associated development of small-cell lung cancer[J],Brachytherapy,2012,11(5):354-358.
[4]Lorensen W E,Cline H E.Marching cubes:A high resolution3Dsurface construction algorithm[J].ACM Siggraph Computer Graphics,1987,21(4):163-169.
[5]Zhang Y,Qian J.Dual contouring for domains with topology ambiguity[J].Computer Methods in Applied Mechanics and Engineering,2012,217:34-45.
[6]Jia T,Feng X,Guo X.Improvement of marching cubes algorithm based on sign determination[C]∥Proc of 2010 2nd IEEE International Workshop on Education Technology and Computer Science(ETCS),2010:194-197.
[7]Raman S,Wenger R.Quality isosurface mesh generation using an extended marching cubes lookup table[J].Computer Graphics Forum,2010,27(3):791-798.
[8]Jiang H,Zhang R,Liu H,et al.An improved 3Dreconstruction method[C]∥Proc of IEEE International Conference on Multimedia Information Networking and Security(MINES’09),2009:209-212.
[9]Liu Y,Xing H L.A boundary focused quadrilateral mesh generation algorithm for multi-material structures[J].Journal of Computational Physics,2013,232(1):516-528.
[10]Lederman C,Joshi A,Dinov I,et al.The generation of tetrahedral mesh models for neuroanatomical MRI[J].NeuroImage,2011,55(1):153-164.
[11]Tian S,Liu G,Weng Z,et al.Isosurface generation algorithm based on spatial discrete points with attribute[J].Journal of Convergence Information Technology,2012,7(4):86-96.
[12]Li Zhong-hua,Zhou Fu-gen,Bai Xiang-zhi.A GPU-based fast volume CT reconstructive algorithm method[J].Journal of Biomedical Engineering,2011,28(2):238-242.(in Chinese)
[13]Li Jia,Long Peng-cheng,Luo Yue-tong,et al.Research and application of a fast 3D-reconstruction method in accurate radiotherapy treatment planning system[J].Nuclear Techniques,2010,33(1):69-74.(in Chinese)
[14]Chen Chun-xiao,Zhang Juan.3Dreconstruction of multiple views based on trifocal tensor[J].Journal of Biomedical Engineering,2012,29(4):769-774.(in Chinese)
[15]Lehmann G.Label objects representation and manipulation with ITK[J].The Insight Journal,2008:1-34.
[16]Xu Xiao-ling,Li Xian-min,Li Gui-qing,et al.A fast algorithm of marching cubes in volume reconstruction[J].Journal of System Simulation,2002,14(4):509-513.(in Chinese)
[17]Vollmer J,Mencl R,Mueller H.Improved laplacian smoothing of noisy surface meshes[J].Computer Graphics Forum,1999,18(3):131-138.
[18]Wang H,Chen H,Su Z,et al.Versatile surface detail editing via Laplacian coordinates[J].The Visual Computer,2011,27(5):401-411.
[1]王志远,李树祥.三维近距离放射治疗计划系统的研究[J].医疗卫生装备,2000,21(1):7-10.
[12]李忠华,周付根,白相志.一种基于GPU的体积CT快速重建算法[J].生物医学工程学杂志,2011,28(2):238-242.
[13]李佳,龙鹏程,罗月童,等.精确放射治疗计划系统中快速三维重建方法的研究与应用[J].核技术,2010,33(1):69-74.
[14]陈春晓,张娟.基于三焦点张量的多视图三维重构[J].生物医学工程学杂志,2012,29(4):769-774.
[16]徐晓玲,李现民,李桂清,等.体素重建中的快速移动立方体方法[J].系统仿真学报,2002,14(4):509-513.