[1]沈傲东,罗立民.基于蒙特卡罗积分的数字影像重建方法[J].东南大学学报(自然科学版),2008,38(6):992-995.[doi:10.3969/j.issn.1001-0505.2008.06.011]
 Shen Aodong,Luo Limin.Digitally reconstructed radiography based on Monte Carlo integral[J].Journal of Southeast University (Natural Science Edition),2008,38(6):992-995.[doi:10.3969/j.issn.1001-0505.2008.06.011]
点击复制

基于蒙特卡罗积分的数字影像重建方法()
分享到:

《东南大学学报(自然科学版)》[ISSN:1001-0505/CN:32-1178/N]

卷:
38
期数:
2008年第6期
页码:
992-995
栏目:
计算机科学与工程
出版日期:
2008-11-20

文章信息/Info

Title:
Digitally reconstructed radiography based on Monte Carlo integral
作者:
沈傲东 罗立民
东南大学影像科学与技术实验室, 南京 210096
Author(s):
Shen Aodong Luo Limin
Laboratory of Image Science and Technology, Southeast University, Nanjing 210096, China
关键词:
蒙特卡罗积分 数字影像重建 重要性采样
Keywords:
Monte Carlo integral digitally reconstructed radiography importance sample
分类号:
TP391.41
DOI:
10.3969/j.issn.1001-0505.2008.06.011
摘要:
提出了一种基于蒙特卡罗积分的数字影像重建方法,首先根据概率密度函数采样形成点云集合,并将空间采样转换为强度采样,采样速度获得了较大提升; 然后利用空域滤波器平滑点云集,避免投影出现孔洞; 最后根据不同的视角对样本投影,统计出像平面上每个区域内的投影数目.为降低图像估计方差,提出结合Russian roulette技术的混合采样方法和Metropolis采样方法.相比较于复杂度为O(N3)3的确定性投影算法,所提出的算法复杂度降为O(N2)2.实验结果表明,该方法能以较快的帧速给出大数据集的数字影像重建,且无需多次调整转换函数即可生成类X射线影像.
Abstract:
A digitally reconstructed radiography(DRR)method based on Monte Carlo integral is presented. The point cloud set is formed firstly according to the probability density function by sampling the intensity instead of the position and the speed of sampling has been improved much more. Then the space filter is applied to smooth the point cloud to avoid the projecting holes. Finally the number of the projecting points is counted from the different views on each area of the image plane. Russian roulette and Metropolis techniques are utilized to decrease the variance of generated image. Compared with the deterministic methods, the algorithm complexity is decreased from O(N3)3 to O(N2)2 for the present method. The experiments demonstrate that the proposed method can obtain the DRR with a higher frame speed for the large data set and without more adjustment for the transfer function.

参考文献/References:

[1] Freud N,Letang J M,Babot D.A hybrid approach to simulate X-ray imaging techniques,combining Monte Carlo and deterministic algorithms[C] //IEEE Nuclear Science Symposium Conference Record.Rome,2004,5:3075-3079.
[2] Rohlfing T,Russakoff D B,Denzler J,et al.Progressive attenuation fields:fast 2D-3D image registration without precomputation[J].Medical Physics,2005,32(9):2870-2880.
[3] Ahn S K,Cho G,Chi Y K,et al.A computer code for the simulation of X-ray imaging systems[C] //IEEE Nuclear Science Symposium Conference Record.Portland,2003,2:838-842.
[4] Fox C,Romeijn H E,Dempsey J F.Fast voxel and polygon ray-tracing algorithms in intensity modulated radiation therapy treatment planning[J].Medical Physics,2006,33(5):1364-1371.
[5] Birkfellner W,Seemann R,Figl M,et al.Wobbled splatting-a fast perspective volume rendering method for simulation of X-ray images from CT[J]. Physics in Medicine and Biology,2005,50(9):N73-N84.
[6] Cai W,Sakas G.Drr volume rendering using splatting in shear warp context[C] //IEEE Nuclear Science Symposium Conference Record.Lyon,2000,3:19/12-19/17.
[7] Zhao H,Reader A.Fast ray-tracing technique to calculate line integral paths in voxel arrays[C] //IEEE Nuclear Science Symposium Conference Record.Portland,2003,4:2808-2812.
[8] Ntasis E,Maniatis T A,Nikita K S.Fourier volume rendering for real time preview of digital reconstructed radiographs:a web-based implementation[J]. Comput Med Imaging Graph,2002,26(1):1-8.
[9] Rogers D W O.Fifty years of Monte Carlo simulations for medical physics[J]. Physics in Medicine and Biology,2006,51(13):R287-R301.
[10] Sakamoto N,Nonaka J,Koyamada K,et al.Particle based volume rendering[C] //Proc of Asia-Pacific Symposium on Visualization.Sydney,2007:129-132.
[11] Csebfalvi B,Szirmay-Kalos S K.Monte carlo volume rendering[C] //Proceedings of the 14th IEEE Visualization Conference.Seattle,2003:59.
[12] Csebfalvi B.Interactive transfer function control for Monte Carlo volume rendering[C] //IEEE Symposium on Volume Visualization and Graphics(VV’04).Austin,2004:33-38.
[13] Li X,Yang J,Zhu Y.Digitally reconstructed radiograph generation by an adaptive Monte Carlo method[J].Physics in Medicine and Biology,2006,51(11):2745-2752.
[14] Horvitz D G,Thompson D J.A generalization of sampling without replacement from a finite universe[J]. Journal of the American Statistical Association,1952,47(260):663-685.
[15] Metropolis N,Rosenbluth A W,Rosenbluth M N,et al.Equation of state calculations by fast computing machines[J]. Journal of Chemical Physics,1953, 21(6):1087-1092.
[16] Lorensen B.The visible human project CT datasets[EB/OL].(2007-01-01)[2008-03-25].https://mri.radiology.uiowa.edu/VHDicom/index.html.

备注/Memo

备注/Memo:
作者简介: 沈傲东(1978—),男,博士生; 罗立民(联系人),男,博士,教授,博士生导师,luo.list@seu.edu.cn.
基金项目: 国家重点基础研究发展计划(973计划)资助项目(2003CB716102).
引文格式: 沈傲东,罗立民.基于蒙特卡罗积分的数字影像重建方法[J].东南大学学报:自然科学版,2008,38(6):992-995.
更新日期/Last Update: 2008-11-20