Monte-Carlo Ray-Tracing for Realistic Interactive Ultrasound Simulation

Ray-based simulations have been shown to generate impressively realistic ultrasound images in interactive frame rates. Recent efforts used GPU-based surface ray-tracing to simulate complex ultrasound interactions such as multiple reflections and refractions. These methods are restricted to perfectly specular reflections (i.e., following only a single reflective/refractive ray), whereas real tissue exhibits roughness of varying degree at tissue interfaces, causing partly diffuse reflections and refractions. Such surface interactions are significantly more complex and can in general not be handled by such deterministic ray-tracing approaches. However, they can be efficiently computed by Monte-Carlo sampling techniques, where many ray paths are generated with respect to a probability distribution. In this paper we introduce Monte-Carlo ray-tracing for ultrasound. This enables the realistic simulation of ultrasound interactions such as soft shadows and fuzzy reflections. We discuss how to properly weight the contribution of each ray path in order to simulate the behavior of a beamformed ultrasound signal. Tracing many individual rays per transducer element is easily parallelizable on modern GPUs, as opposed to previous approaches based on recursive binary ray-tracing.

[1]  T. Varslot,et al.  Computer simulation of forward wave propagation in soft tissue , 2005, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[2]  Michel Bertrand,et al.  Ultrasonic texture motion analysis: theory and simulation , 1995, IEEE Trans. Medical Imaging.

[3]  Hans Hagen,et al.  Towards a Virtual Echocardiographic Tutoring System , 2008, Visualization in Medicine and Life Sciences.

[4]  Nassir Navab,et al.  Automatic CT-ultrasound registration for diagnostic imaging and image-guided intervention , 2008, Medical Image Anal..

[5]  T. Varslot,et al.  Forward propagation of acoustic pressure pulses in 3D soft biological tissue , 2006 .

[6]  Andrew S. Glassner,et al.  An introduction to ray tracing , 1989 .

[7]  J M Thijssen,et al.  Texture of B-mode echograms: 3-D simulations and experiments of the effects of diffraction and scatterer density. , 1985, Ultrasonic imaging.

[8]  J. D'hooge,et al.  A fast convolution-based methodology to simulate 2-Dd/3-D cardiac ultrasound images , 2009, IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.

[9]  Benny Bürger,et al.  Real-Time GPU-Based Ultrasound Simulation Using Deformable Mesh Models , 2013, IEEE Transactions on Medical Imaging.

[10]  Orcun Goksel,et al.  Scatterer reconstruction and parametrization of homogeneous tissue for ultrasound image simulation , 2015, 2015 37th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC).

[11]  Roman Kuchkuda,et al.  An introduction to ray tracing , 1993, Comput. Graph..

[12]  J C Bamber,et al.  Ultrasonic B-scanning: a computer simulation , 1980, Physics in medicine and biology.

[13]  Mario Ries,et al.  Stochastic ray tracing for simulation of high intensity focal ultrasound therapy. , 2014, The Journal of the Acoustical Society of America.

[14]  Arthur Appel,et al.  Some techniques for shading machine renderings of solids , 1968, AFIPS Spring Joint Computing Conference.

[15]  Nassir Navab,et al.  Fast Ultrasound Image Simulation Using the Westervelt Equation , 2010, MICCAI.

[16]  Nassir Navab,et al.  Patient-specific 3D Ultrasound Simulation Based on Convolutional Ray-tracing and Appearance Optimization , 2015, MICCAI.

[17]  Olivier Salvado,et al.  Ultrasound goes GPU: real-time simulation using CUDA , 2009, Medical Imaging.

[18]  Orcun Goksel,et al.  Image-Based PSF Estimation for Ultrasound Training Simulation , 2016, SASHIMI@MICCAI.

[20]  Leonidas J. Guibas,et al.  Metropolis light transport , 1997, SIGGRAPH.

[21]  Turner Whitted,et al.  An improved illumination model for shaded display , 1979, CACM.

[22]  Eric P. Lafortune,et al.  Mathematical Models and Monte Carlo Algorithms for Physically Based Rendering , 1995 .

[23]  Jørgen Jensen,et al.  Simulation of advanced ultrasound systems using Field II , 2004, 2004 2nd IEEE International Symposium on Biomedical Imaging: Nano to Macro (IEEE Cat No. 04EX821).

[24]  Jean Meunier,et al.  A Model For Dynamic Texture Analysis In Two-Dimensional Echocardiograms Of The Myocardium , 1987, Other Conferences.