Book of Abstracts 15th International Symposium on Computer Methods in Biomechanics and Biomedical Engineering and 3rd Conference on Imaging and Visualization

In this edition, the two events will run together as a single conference, highlighting the strong connection with the Taylor & Francis journals: Computer Methods in Biomechanics and Biomedical Engineering (John Middleton and Christopher Jacobs, Eds.) and Computer Methods in Biomechanics and Biomedical Engineering: Imaging and Visualization (JoãoManuel R.S. Tavares, Ed.). The conference has become a major international meeting on computational biomechanics, imaging andvisualization. In this edition, the main program includes 212 presentations. In addition, sixteen renowned researchers will give plenary keynotes, addressing current challenges in computational biomechanics and biomedical imaging. In Lisbon, for the first time, a session dedicated to award the winner of the Best Paper in CMBBE Journal will take place. We believe that CMBBE2018 will have a strong impact on the development of computational biomechanics and biomedical imaging and visualization, identifying emerging areas of research and promoting the collaboration and networking between participants. This impact is evidenced through the well-known research groups, commercial companies and scientific organizations, who continue to support and sponsor the CMBBE meeting series. In fact, the conference is enriched with five workshops on specific scientific topics and commercial software.info:eu-repo/semantics/draf

Is the Rayleigh-Sommerfeld diffraction always an exact reference for high speed diffraction algorithms?

In several areas of optics and photonics like wave propagation, digital holography, holographic microscopy, diffraction imaging, biomedical imaging and diffractive optics, the behavior of the electromagnetic waves has to be calculated with the scalar theory of diffraction by computational methods. Many of these high speed diffraction algorithms based on a fast Fourier transformation are in principle approximations of the Rayleigh-Sommerfeld Diffraction (RSD) theory. However, to investigate their numerical accuracy, they should be compared with and verified by RSD. All numerical simulations are in principle based on a sampling of the analogue continuous field. In this article we demonstrate a novel validity condition for the well-sampling in RSD, which makes a systematic treatment of sampling in RSD possible. We show the fundamental restrictions due to this condition and the anomalies caused by its violation. We also demonstrate that the restrictions are completely removed by a sampling below the Abbe resolution limit. Furthermore, we present a very general unified approach for applying the RSD outside its validity domain by the combination of a forward and reverse calculation

Nanoinformatics: developing new computing applications for nanomedicine

Nanoinformatics has recently emerged to address the need of computing applications at the nano level. In this regard, the authors have participated in various initiatives to identify its concepts, foundations and challenges. While nanomaterials open up the possibility for developing new devices in many industrial and scientific areas, they also offer breakthrough perspectives for the prevention, diagnosis and treatment of diseases. In this paper, we analyze the different aspects of nanoinformatics and suggest five research topics to help catalyze new research and development in the area, particularly focused on nanomedicine. We also encompass the use of informatics to further the biological and clinical applications of basic research in nanoscience and nanotechnology, and the related concept of an extended ?nanotype? to coalesce information related to nanoparticles. We suggest how nanoinformatics could accelerate developments in nanomedicine, similarly to what happened with the Human Genome and other -omics projects, on issues like exchanging modeling and simulation methods and tools, linking toxicity information to clinical and personal databases or developing new approaches for scientific ontologies, among many others

Eigenspace-Based Minimum Variance Combined with Delay Multiply and Sum Beamformer: Application to Linear-Array Photoacoustic Imaging

In Photoacoustic imaging, Delay-and-Sum (DAS) algorithm is the most commonly used beamformer. However, it leads to a low resolution and high level of sidelobes. Delay-Multiply-and-Sum (DMAS) was introduced to provide lower sidelobes compared to DAS. In this paper, to improve the resolution and sidelobes of DMAS, a novel beamformer is introduced using Eigenspace-Based Minimum Variance (EIBMV) method combined with DMAS, namely EIBMV-DMAS. It is shown that expanding the DMAS algebra leads to several terms which can be interpreted as DAS. Using the EIBMV adaptive beamforming instead of the existing DAS (inside the DMAS algebra expansion) is proposed to improve the image quality. EIBMV-DMAS is evaluated numerically and experimentally. It is shown that EIBMV-DMAS outperforms DAS, DMAS and EIBMV in terms of resolution and sidelobes. In particular, at the depth of 11 mm of the experimental images, EIBMV-DMAS results in about 113 dB and 50 dB sidelobe reduction, compared to DMAS and EIBMV, respectively. At the depth of 7 mm, for the experimental images, the quantitative results indicate that EIBMV-DMAS leads to improvement in Signal-to-Noise Ratio (SNR) of about 75% and 34%, compared to DMAS and EIBMV, respectively.Comment: arXiv admin note: substantial text overlap with arXiv:1709.0796

Quasi-static imaged-based immersed boundary-finite element model of human left ventricle in diastole

SUMMARY: Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid–structure interaction simulations of the heart

Exponential Krylov time integration for modeling multi-frequency optical response with monochromatic sources

Light incident on a layer of scattering material such as a piece of sugar or white paper forms a characteristic speckle pattern in transmission and reflection. The information hidden in the correlations of the speckle pattern with varying frequency, polarization and angle of the incident light can be exploited for applications such as biomedical imaging and high-resolution microscopy. Conventional computational models for multi-frequency optical response involve multiple solution runs of Maxwell's equations with monochromatic sources. Exponential Krylov subspace time solvers are promising candidates for improving efficiency of such models, as single monochromatic solution can be reused for the other frequencies without performing full time-domain computations at each frequency. However, we show that the straightforward implementation appears to have serious limitations. We further propose alternative ways for efficient solution through Krylov subspace methods. Our methods are based on two different splittings of the unknown solution into different parts, each of which can be computed efficiently. Experiments demonstrate a significant gain in computation time with respect to the standard solvers.Comment: 22 pages, 4 figure

Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models

Significance: Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop an approach to reduce this bottleneck, which has significant implications for quantitative tomographic imaging in a variety of medical and industrial applications.Aim: Our aim is to enable computationally efficient image reconstruction in (hybrid) diffuse optical modalities using stochastic forward models.Approach: Using Monte Carlo, we compute a fully stochastic gradient of an objective function for a given imaging problem. Leveraging techniques from the machine learning community, we then adaptively control the accuracy of this gradient throughout the iterative inversion scheme to substantially reduce computational resources at each step.Results: For example problems of quantitative photoacoustic tomography and ultrasound-modulated optical tomography, we demonstrate that solutions are attainable using a total computational expense that is comparable to (or less than) that which is required for a single high-accuracy forward run of the same Monte Carlo model.Conclusions: This approach demonstrates significant computational savings when approaching the full nonlinear inverse problem of optical property estimation using stochastic methods