In our approach, we first specify a 3D model of the image structure by segmenting a CT dataset into the respective tissues followed by assigning the acoustic properties (velocity, impedance, scattering mean and standard deviation, damping factor and packing factor). Given that model, we simulate the ray propagation, beam forming, and finally the backscattering. Due to the inhomogenities of tissue, different physical models for ultrasound simulation are required: Rayleigh scattering is applied for homogenous regions and ray tracing techniques handle abrupt changes in acoustic impedance on tissue boundaries. The latter leads to different phenomena like refraction (Snell's law), reflection and transmission (Fresnel equation). The gradients needed for these methods are precomputed for each model using a central-difference method with multiple neighbours. Absorption is calculated by the Beer-Lambert law.
[1]
J. Jensen,et al.
Calculation of pressure fields from arbitrarily shaped, apodized, and excited ultrasound transducers
,
1992,
IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control.
[2]
Daniel Cohen-Or,et al.
Real-Time Ultrasound Imaging Simulation
,
1998,
Real Time Imaging.
[3]
Thomas W. Sederberg,et al.
Free-form deformation of solid geometric models
,
1986,
SIGGRAPH.
[4]
Guy Cloutier,et al.
Simulation of ultrasound backscattering by red cell aggregates: effect of shear rate and anisotropy.
,
2002,
Biophysical journal.
[5]
L. Rayleigh,et al.
The theory of sound
,
1894
.