This paper presents a multi-dimensional point scatterer distribution model for the context of ultrasound image simulation. The model has a simple parameterisation, has low computational requirements and is flexible enough to model spatial organisation of scatterers ranging from highly clustered to nearly regular. The model extends an existing 1D model by mapping 1D scatterer positions to a Hilbert space-filling curve. The flexibility of the heuristic model is illustrated through experiments where common statistical models of ultrasonic speckle are fitted to simulated data. The results agree with theoretical predictions. Index terms: acoustic imaging, simulation, fractals, statistics, speckle
BACKGROUNDAn ultrasound (US) image simulator is a useful validation platform for image processing applications where ground truth pertaining to image content or acquisition conditions is needed but not available from US data without additional equipment. It also serves as a test platform during application software development when data acquisition is impossible or awkward. Typical simulators take a list of point scatterers with their strength and position as input along with US transducer specifications and simulate the resulting backscattered signal.The density and spatial organisation of scatterers are two of many factors which determine the appearance and statistics of US speckle. In order to validate image processing applications or study image properties across different tissue types, it is useful to be able to generate random lists of point scatterers with varying density and spatial organisation, ranging from highly clustered to random to nearly regular. There is currently a lack of models displaying such flexibility along with computational efficiency and simple parameterisation.The Neyman-Scott point process model has been used for studying the effect of blood cell aggregation on ultrasonic signals [1]. The model creates random cluster centers and spawns daughter points surrounding them. It is not appropriate for generating quasi-periodic patterns. Such patterns can Catherine Laporte is funded by a doctoral scholarship from the Natural Science and Engineering Research Council of Canada.be generated with varying regularity by perturbing a regular point lattice [2], but this approach cannot generate clustered point patterns. An alternative is a Gibbs-Markov area interaction process which imposes pairwise repulsive and attractive constraints between points. This model was used to study US backscattering from aggregates of non-overlapping blood cells [3]. This type of model is computationally demanding and fails to produce a broad enough variety of clustered patterns when no repulsive constraints are used to maintain a minimal distance between points [4]. In short, previous attempts to parameterise variation in the spatial organisation of scatterers have been geared towards specific applications and have not modeled the full continuum of spatial organisations ranging from clustered to regular. One notable excep...