Search

With the goal of developing an accurate and fast image reconstruction algorithm for ultrasound computed tomography, we combine elements of model- and data-driven approaches and propose a learned method which addresses the disadvantages of both approaches. We design a deep neural network which accounts for a nonlinear forward operator and primal-dual algorithm by its inherent network architecture. The network is trained end-to-end, with ultrasound pressure field data as input to get directly an optimized reconstruction of speed of sound and attenuation images. The training and test data are based on a set of Optical and Acoustic Breast Phantom Database, where we use the image as ground truth and simulate pressure field data according to our forward model. Extensive experiments show that our method achieves significant improvements over state-of-the-art reconstruction methods in this field. Experiments show that the proposed algorithm improves the measures structural similarity measure (SSIM) from 0.74 to 0.95 and root mean squared error (RMSE) from 0.13 to 0.09 on average concerning the speed of sound reconstruction, while it improves the SSIM from 0.60 to 0.94 and RMSE from 0.24 to 0.10 on average in attenuation reconstruction.

Objective. The image reconstruction of ultrasound computed tomography is computationally expensive with conventional iterative methods. The fully learned direct deep learning reconstruction is promising to speed up image reconstruction significantly. However, for direct reconstruction from measurement data, due to the lack of real labeled data, the neural network is usually trained on a simulation dataset and shows poor performance on real data because of the simulation-to-real gap. Approach. To improve the simulation-to-real generalization of neural networks, a series of strategies are developed including a Fourier-transform-integrated neural network, measurement-domain data augmentation methods, and a self-supervised-learning-based patch-wise preprocessing neural network. Our strategies are evaluated on both the simulation dataset and real measurement datasets from two different prototype machines. Main results. The experimental results show that our deep learning methods help to improve the neural networks’ robustness against noise and the generalizability to real measurement data. Significance. Our methods prove that it is possible for neural networks to achieve superior performance to traditional iterative reconstruction algorithms in imaging quality and allow for real-time 2D-image reconstruction. This study helps pave the path for the application of deep learning methods to practical ultrasound tomography image reconstruction based on simulation datasets.

Deep neural networks have proven to excel classical medical image reconstruction techniques. Some networks are based on fully connected (FC) layers to achieve domain transformation such as from the data acquisition domain to the image domain. However, FC layers result in huge numbers of parameters which take a lot of GPU memory. Hence, they do not scale well, and the overall performance is limited. For ultrasound computed tomography (USCT) application, we propose a memory-efficient convolutional network that reconstructs images from the frequency domain to image domain with much less parameters compared with multilayer perceptron, by using data-driven learning. Extensive experiments demonstrate that our method achieves high reconstruction quality. It improves the structural similarity measure (SSIM) from 0.73 to 0.99 when compared with state-of-the-art reconstruction methods in this field while reduces 2/3 parameters when compared with deep neural network with FC layers to reconstruct images from frequency domain to image domain.

Ultrasound transmission tomography promises a high potential and novel imaging method for early breast cancer diagnosis; it can quantitatively characterize tissues or materials by the attenuation and speed of sound (SoS). Reconstruction of ultrasound transmission tomography is an inverse problem that can be solved iteratively based on a paraxial approximation of the Helmholtz equation as forward model, which is highly non-linear and time-consuming. In order to address these problems and reconstruct desired images, we design a dual domain network architecture for ultrasound transmission tomography reconstruction. It can enhance the information of measurement domain and directly reconstruct from pressure field measurements without using any initialization of reconstruction and fully connected layer. We train the network on simulated ImageNet data and transfer it for ultrasound transmission tomography images to avoid overfitting when the amount of ultrasound transmission tomography images is limited. Our experimental results demonstrate that a dual domain network produces significant improvements over state-of-the-art methods. It improves the measured structural similarity measure (SSIM) from 0.54 to 0.90 and normalized root mean squared error (nRMSE) from 0.49 to 0.01 on average concerning the SoS reconstruction, and from 0.46 to 0.98 for SSIM, from 353 to 0.03 for nRMSE on average concerning the attenuation reconstruction.

Breast cancer is the most common cancer type among women. Approximately 40,000 women are expected to die from breast cancer every year. While digital mammography has a central role in the early diagnosis of breast cancer, many cancers are not visible in mammography, for example in women with dense breast tissue. Contrast enhanced magnetic resonance imaging (CE-MRI) of the breast is often used to detect lesions not visible in mammography. Lesions with suspicious characteristics on CE-MRI need to be further assessed with MRI-guided biopsy. However, MRI-guided biopsy is expensive, time consuming, and not widely available. In this paper, a novel method for a matching tool between MRI and spot mammograms is proposed. Our aim is to transfer information that is only visible in MRI onto mammographic spot projections, to enable X-ray guided biopsy even if the lesion is only visible in MRI. Two methods of registration in combination are used; a biomechanical model based registration between MRI and full view X-ray mammograms and a subsequent image based registration between full mammograms and spot mammograms. Preliminary results assessed for one patient from the Medical University of Vienna are presented. The target registration error (TRE) of biomechanical model based registration is 2.4 mm and the TRE of the image based registration is 9.5 mm. The total TRE of the two steps is 7.3 mm.

Ultrasound computer tomography (USCT) is a promising modality for breast cancer diagnosis which images the reflectivity, sound speed and attenuation of tissue. Elastic properties of breast tissue, however, cannot directly be imaged although they have shown to be applicable as a discriminator between different tissue types. In this work we propose a novel approach combining USCT with the principles of strain elastography. Socalled USCT-SE makes use of imaging the breast in two deformation states, estimating the deformation field based on reconstructed images and thereby allows localizing and distinguishing soft and hard masses. We use a biomechanical model of the breast to realistically simulate both deformation states of the breast. The analysis of the strain is performed by estimating the deformation field from the deformed to the undeformed image by a non-rigid registration. In two experiments the non-rigid registration is applied to ground truth sound speed images and simulated SAFT images. Results of the strain analysis show that for both cases soft and hard lesions can be distinguished visually in the elastograms. This paper provides a first approach to obtain mechanical information based on external mechanical excitation of breast tissue in a USCT system.

Focused ultrasound (FUS) is a promising method to open the blood brain barrier (BBB) for treatment of neurodegenerative diseases. Accurate targeting is essential for a successful BBB opening (BBBo). We aim to develop a robust therapy planning for BBBo in mice, which is challenging due to the size of the brain and the influence of the skull on the ultrasound pressure distribution. For enabling mouse individual therapy planning, a simulation tool is proposed, developed and validated. We used the k-Wave toolbox to enable 3D acoustic simulations of the commercial FUS system from Image Guided Therapy (IGT). Micro-CT scans were used to model the geometry of skulls. Simulations using a mouse skull showed an attenuation of approx. 20–24% depending on the position of penetration, which was validated by hydrophone measurements in the same range. Based on these validations we planned BBBo in m ice by placing the transducer at different positions over the mouse brain and varying the excitation amplitude. With different transducer positions, the peak pressure in the brain varied between 0.54 MPa and 0.62 MPa at 11% output level, which is expected to enable safe BBBo. Subsequently, in vivo experiments were conducted using the aforementioned simulation parameters. BBBo was confirmed by contrast enhanced T1 weighted magnetic resonance images immediately after sonication.

Conventional iterative methods for image reconstruction in ultrasound transmission tomography need to perform many iterations where at each iteration one has to compute the complex forward model of ultrasound wave propagation, and hence they are time-consuming. We use a U-net neural network to accelerate the reconstruction, by training the network to map from an initial reconstruction obtained via a few iterations of L-BFGS method to the target ground truth image. Since the computation of a forward pass of the neural network is very fast, we can expect a significant acceleration using the trained network for image reconstruction. Experiments show that our trained network can replace 40 L-BFGS iterations to generate equivalent reconstructions with slightly better quantitative quality in terms of normalized root mean square error and better visual quality due to the network's denoising effect. It can achieve up to 283× speedup compared with L-BFGS method for reconstructing small-size sound speed images with 80×80 pixels. This implies that we can expect even greater acceleration effects when applying such approach to reconstruct large-size 3D images.

Ultrasound transmission tomography is a promising modality for breast cancer diagnosis. For image reconstruc- tion approximations to the acoustic wave equation such as straight or bent rays are commonly used due to their low computational complexity. For sparse apertures the coverage of the volume by rays is very limited, thereby requiring strong regularization in the inversion process. The concept of fat rays reduces the sparseness and includes the contributions to the measured signal originating from the first Fresnel zone. In this work we investi- gate the application of the fat ray concept to ultrasound transmission tomography. We implement a straight ray, bent ray and fat ray forward model. For the inversion process a least squares solver (LSQR), a simultaneous al- gebraic reconstruction technique (SART) and a compressive sensing based total variation minimization (TVAL3) is applied. The combination of forward models and inversion processes has been evaluated by synthetic data. TVAL3 outperforms SART and LSQR, especially for sparse apertures. The fat ray concept is able to decrease the error with respect to the ground truth compared to the bent ray method especially for SART and LSQR inversion, and especially for very sparse apertures.

A 3D ultrasound computer tomography (USCT) device with a nearly isotropic and spatially invariant 3D point spread function has been constructed at Institute for Data Processing and Electronic (IPE), Karlsruhe Institute of Technology (KIT). This device is currently applied in clinical studies for breast cancer screening. In this paper, a new method to develop an automated segmentation algorithm for USCT acquired images is proposed. The method employs distance regularized level set evolutionary (DRLSE) active contours along with surface fitting extrapolation and 3D binary mask generation for fully automatic segmentation outcome. In the first stage of the proposed algorithm, DRLSE is applied to those 3D USCT slice images which contain breast and are less affected by noise and ring artifacts named as Cat2. The DRLSE segmentation results are employed to extrapolate the rest of slice images known as Cat1. To overcome defectively segmented slice images, a 3D binary mask is generated out of USCT attenuation images. The 3D binary mask is multiplied by the DRLSE-based segmentation results to form finally segmented 3D USCT images. The method was tested on 12 clinical dataset images. According to F-measure criterion, the proposed method shows higher performance than the previously proposed semiautomatic segmentation one.

A promising candidate for improved imaging of breast cancer is ultrasound computer tomography (USCT). Current experimental USCT systems are still focused in elevation dimension resulting in a large slice thickness, limited depth of field, loss of out-of-plane reflections, and a large number of movement steps to acquire a stack of images. 3D USCT emitting and receiving spherical wave fronts overcomes these limitations. We built an optimized 3D USCT, realizing for the first time the full benefits of a 3D system. The point spread function could be shown to be nearly isotropic in 3D, to have very low spatial variability and fit the predicted values. The contrast of the phantom images is very satisfactory in spite of imaging with a sparse aperture. The resolution and imaged details of the reflectivity reconstruction are comparable to a 3 T MRI volume. Important for the obtained resolution are the simultaneously obtained results of the transmission tomography. The KIT 3D USCT was then tested in a pilot study on ten patients. The primary goals of the pilot study were to test the USCT device, the data acquisition protocols, the image reconstruction methods and the image fusion techniques in a clinical environment. The study was conducted successfully; the data acquisition could be carried out for all patients with an average imaging time of six minutes per breast. The reconstructions provide promising images. Overlaid volumes of the modalities show qualitative and quantitative information at a glance. This paper gives a summary of the involved techniques, methods, and first results.

Image reconstruction of ultrasound computed tomography based on the wave equation is able to show much more structural details than simpler ray-based image reconstruction methods. However, to invert the wave-based forward model is computationally demanding. To address this problem, we develop an efficient fully learned image reconstruction method based on a convolutional neural network. The image is reconstructed via one forward propagation of the network given input sensor data, which is much faster than the reconstruction using conventional iterative optimization methods. To transform the ultrasound measured data in the sensor domain into the reconstructed image in the image domain, we apply multiple down-scaling and up-scaling convolutional units to efficiently increase the number of hidden layers with a large receptive and projective field that can cover all elements in inputs and outputs, respectively. For dataset generation, a paraxial approximation forward model is used to simulate ultrasound measurement data. The neural network is trained with a dataset derived from natural images in ImageNet and tested with a dataset derived from medical images in OA-Breast Phantom dataset. Test results show the superior efficiency of the proposed neural network to other reconstruction algorithms including popular neural networks. When compared with conventional iterative optimization algorithms, our neural network can reconstruct a 110 × 86 image more than 20 times faster on a CPU and 1000 times faster on a GPU with comparable image quality and is also more robust to noise.

In ultrasound transmission tomography, image reconstruction is an inverse problem which is solved iteratively based on a forward model that simulates the wave propagation of ultrasound. A commonly used forward model is paraxial approximation of the Helmholtz equation, which is time-consuming. Hence developing optimizers that minimize the number of forward solutions is crucial to achieve clinically acceptable reconstruction time, while the state-of-the-art methods in this field such as Gauss-Newton conjugate gradient (CG) and nonlinear CG are not capable of reaching this goal. To that end, we focus on Jacobian-free optimizers or accelerators in this paper, since the computation of the Jacobian is expensive. We investigate the limited memory Broyden–Fletcher–Goldfarb–Shanno (L-BFGS) algorithm as a preconditioning technique due to its ability to efficiently approximate inverse Hessian without performing forward model or its adjoint. We show L-BFGS can reach a speedup of more than one order of magnitude for the noise-free case, while the method still halves the reconstruction time in presence of noise in the data. The performance drop is explained by perturbed gradients due to noise in the data. We also show when used alone as a quasi-Newton method, L-BFGS is competitive with the accelerated CG based methods regarding the number of iterations, and outperforms them regarding reconstruction time.

The purpose of this paper is to investigate Nesterov’s accelerated gradient (NAG) method for the reconstruction of speed of sound and attenuation images in ultrasound computed tomography. The inverse problem of reconstruction is tackled via minimizing the deviation between exact measurements and the predicted measurements based on a paraxial approximation of the Helmholtz equation which simulates the ultrasound wave forward propagation. To solve this optimization problem, NAG is performed and compared with other algorithms. Also, a line search method is used to compute the step size for each iteration since finding proper step sizes is crucial for the convergence of such optimization algorithms. The strong Wolfe conditions are adopted as the termination condition for line search. We have compared five algorithms, namely Gauss-Newton conjugate gradient, gradient descent, NAG, gradient descent with line search, and NAG with line search. On one hand, NAG with line search has the fastest convergence rate in respect to the number of used iterations compared to the other methods. However, due to the increased computational complexity of line search for each iteration, it requires extra computational time. On the other hand, NAG with a fixed step size for all iterations is the fastest method among all the tested methods regarding computational time.

Ultrasound transmission tomography can offer quantitative characterization of breast tissue by reconstructing its sound speed and attenuation images. The image reconstruction process is a nonlinear inverse problem which we iteratively solve using the paraxial approximation of wave equation. The problem is tackled via the Gauss-Newton method yielding a set of linear systems and iteratively solving these linear systems. In this paper, we study multigrid methods for solving these linear systems. We test three multigrid schemes including V-cycle, W-cycle, and full multigrid (FMG), with up to four resolution levels. At each grid level, we use the conjugate gradient (CG) method as a standard solver. Our interest is at first by how far these schemes allow us to reduce the computations alone. For performance evaluation, we compare the multigrid methods with fixed-grid CG method where we directly apply CG on the finest grid. Results show that all tested multigrid methods have an accelerating effect in terms of needing fewer CG iterations on the finest grid. The best-case reduction is 32% for V-cycle, 33% for W-cycle, and 27% for FMG. This means that our multigrid scheme has the potential for significantly reducing reconstruction time for large-scale 2D or 3D images, where the computation cost on the finest grid is very high.

Increasing interest in multimodal breast cancer diagnosis has led to the development of methods for MRI to X-ray mammography registration. The severe breast deformation in X-ray mammography is often tackled by biomechanical models, yet there is no common consensus in literature about the required complexity of the deformation model and the simulation strategy. We present for the first time an automated patient-specific biomechanical model based image registration of MRI to digital breast tomosynthesis (DBT). DBT provides three-dimensional information of the compressed breast and as such drives the registration by a volume similarity metric. We compare different simulation strategies and propose a patient-specific optimization of simulation and model parameters. The average three-dimensional breast overlap measured by Dice coefficient of DBT and registered MRI improves for four analyzed subjects by including the estimation of unloaded state, simulation of gravity, and a concentrated pull force that mimics manual positioning of the breast on the plates from 88.1% for a mere compression simulation to 93.1% when including all our proposed simulation steps, whereas additional parameter optimization further increased the value to 94.4%.

Increasing interest in multimodal breast cancer diagnosis has led to the development of methods for MRI to X-ray mammography registration to provide direct correlation of modalities. The severe breast deformation in X-ray mammography is often tackled by biomechanical models, which however have not yet brought the registration accuracy to a clinically applicable level. We present a novel registration approach of MRI to X-ray tomosynthesis. Tomosynthesis provides three-dimensional information of the compressed breast and as such has the ability to open new possibilities in the registration of MRI and X-ray data. By bundling the 3D information from the tomosynthesis volume with the 2D projection images acquired at different measuring angles, we provide a correlation between the registration error in 3D and 2D and evaluate different 3D- and 2D-based similarity metrics to drive the optimization of the automated patient-specific registration approach. From the preliminary study of four analysed patients we found that the projected registration error is in general larger than the 3D error in case of small registration errors in the cranio-caudal direction. Although both image shape and intensitybased 2D similarity metrics showed a clear correlation with the 2D registration error at different projection angles, metrics that relied on the combined 2D and 3D information yielded in most of the cases the minimal registration error and as such had better performance than similarity metrics that rely only on the shape similarity of volumes.

In Ultrasound computer tomography (USCT) Synthetic aperture focusing technique (SAFT) is often applied for reflectivity image reconstruction. Phase aberration correction is essential to cope with the large sound speed differences in water and the different human tissues. In this paper we compare two approaches for phase aberration correction: a straight ray approximation using the Bresenham algorithm (B-SAFT) and a bent ray approximating using a multi-stencil Fast Marching Method (FMM-SAFT). The analysis is carried out with simulated point scatterers and simulated phantoms to measure the effect on the image resolution and contrast. The method is additionally applied to experimental data. B-SAFT degrades the image resolution and contrast in cases of large sound speed differences of objects and if the reconstructed point is close to a boundary where a change in impedance is present. FMM-SAFT is able to recover the image quality in these cases if the sound speed distribution is known accurately and with high resolution. If these requirements cannot be met, B-SAFT proved to be more robust.

Ultrasound Computer Tomography (USCT) is a promising breast imaging modality under development. Comparison to a standard method like mammography is essential for further development. Due to significant differences in image dimensionality and compression state of the breast, correlating USCT images and X-ray mammograms is challenging. In this paper we present a 2D/3D registration method to improve the spatial correspondence and allow direct comparison of the images. It is based on biomechanical modeling of the breast and simulation of the mammographic compression. We investigate the effect of including patient-specific material parameters estimated automatically from USCT images. The method was systematically evaluated using numerical phantoms and in-vivo data. The average registration accuracy using the automated registration was 11.9mm. Based on the registered images a method for analysis of the diagnostic value of the USCT images was developed and initially applied to analyze sound speed and attenuation images based on X-ray mammograms as ground truth. Combining sound speed and attenuation allows differentiating lesions from surrounding tissue. Overlaying this information on mammograms, combines quantitative and morphological information for multimodal diagnosis.

Ultrasound Computer Tomography (USCT) is a promising new imaging system for breast cancer diagnosis. An essential step before further processing is to remove the water background from the reconstructed images. In this paper we present a fully-automated image segmentation method based on three-dimensional active contours. The active contour method is extended by applying gradient vector flow and encoding the USCT aperture characteristics as additional weighting terms. A surface detection algorithm based on a ray model is developed to initialize the active contour, which is iteratively deformed to capture the breast outline in USCT reflection images. The evaluation with synthetic data showed that the method is able to cope with noisy images, and is not influenced by the position of the breast and the presence of scattering objects within the breast. The proposed method was applied to 14 in-vivo images resulting in an average surface deviation from a manual segmentation of 2.7 mm. We conclude that automated segmentation of USCT reflection images is feasible and produces results comparable to a manual segmentation. By applying the proposed method, reproducible segmentation results can be obtained without manual interaction by an expert.

Synthetic Aperture Focusing Technique (SAFT) allows fast data acquisition and optimally focused images. The computational burden for 3D imaging is large as for each voxel the delay for each acquired A-scan has to be calculated, e.g. O(N5) for N3 voxels and N2 A-scans. For 3D reconstruction of objects which are large in terms of the wavelength, e.g. ≥ (100 λ)3, the computation of one volume takes several days on a current multicore PC. If the 3D distribution of the speed of sound is applied to correct the delays, the computation time increases further. In this work a time of flight interpolation based GPU implementation (TOFI-SAFT) is presented which accelerates our previous GPU implementation of speed of sound corrected SAFT by a factor of 7 to 16 min. with only minor reduction of image quality.

Computer aided diagnosis (CAD) of breast cancer is mainly focused on monomodal applications. Here we present a fully automated multimodal CAD, which uses patient-specific image registration of MRI and two-view X-ray mammography. The image registration estimates the spatial correspondence between each voxel in the MRI and each pixel in cranio-caudal and mediolateral-oblique mammograms. Thereby we can combine features from both modalities. As a proof of concept we classify fixed regions of interest (ROI) into normal and suspect tissue. We investigate the classification performance of the multimodal classification in several setups against a classification with MRI features only. The average sensitivity of detecting suspect ROIs improves by approximately 2% when combining MRI with both mammographic views compared to MRI-only detection, while the specificity stays at a constant level. We conclude that automatically combining MRI and X-ray can enhance the result of a breast CAD system.

For iterative image reconstruction of transmission tomography we apply the paraxial approximation of the Helmholtz equation for a spherical transducer arrangement. We choose this approach due to its three order of magnitude lower complexity than full wave solutions with the same precision for transmission tomography. In homogeneous media we prove that our forward solution is exact. With the help of this forward solution 2D and 3D ultrasound measurements could be simulated for transmission tomography. 2D reconstructions of a breast-like numerical phantom had a deviation in sound speed of 0.14 m/s and a deviation in attenuation of 6.5% from the ground truth. Applications up to now are breast cancer diagnostics and non-destructive testing.

In our first clinical study with a full 3D Ultrasound Computer Tomography (USCT) system patient data was acquired in eight minutes for one breast. In this paper the patient movement during the acquisition was analyzed quantitatively and as far as possible corrected in the resulting images. The movement was tracked in ten successive reflectivity reconstructions of full breast volumes acquired during 10 s intervals at different aperture positions, which were separated by 41 s intervals. The mean distance between initial and final position was 2.2 mm (standard deviation (STD) ± 0.9 mm, max. 4.1 mm, min. 0.8 mm) and the average sum of all moved distances was 4.9 mm (STD ± 1.9 mm, max. 8.8 mm, min. 2.7 mm). The tracked movement was corrected by summing successive images, which were transformed according to the detected movement. The contrast of these images increased and additional image content became visible.

Ultrasound Computer Tomography (USCT) is a promising new imaging method for breast cancer diagnosis. We developed a 3D USCT system and tested it in a pilot study with encouraging results: 3D USCT was able to depict two carcinomas, which were present in contrast enhanced MRI volumes serving as ground truth. To overcome severe differences in the breast shape, an image registration was applied. We analyzed the correlation between average sound speed in the breast and the breast density estimated from segmented MRIs and found a positive correlation with R=0.70. Based on the results of the pilot study we now carry out a successive clinical study with 200 patients. For this we integrated our reconstruction methods and image post-processing into a comprehensive workflow. It includes a dedicated DICOM viewer for interactive assessment of fused USCT images. A new preview mode now allows intuitive and faster patient positioning. We updated the USCT system to decrease the data acquisition time by approximately factor two and to increase the penetration depth of the breast into the USCT aperture by 1 cm . Furthermore the compute-intensive reflectivity reconstruction was considerably accelerated, now allowing a sub-millimeter volume reconstruction in approximately 16 minutes. The updates made it possible to successfully image first patients in our ongoing clinical study.

Ultrasound Computer Tomography (USCT), developed at KIT, is a promising new imaging system for breast cancer diagnosis, and was successfully tested in a pilot study. The 3D USCT II prototype consists of several hundreds of ultrasound (US) transducers on a semi-ellipsoidal aperture. Spherical waves are sequentially emitted by individual transducers and received in parallel by many transducers. Reflectivity volumes are reconstructed by synthetic aperture focusing (SAFT). However, straight forward SAFT imaging leads to blurred images due to system imperfections. We present an extension of a previously proposed approach to enhance the images. This approach includes additional a priori information and system characteristics. Now spatial phase compensation was included. The approach was evaluated with a simulation and clinical data sets. An increase in the image quality was observed and quantitatively measured by SNR and other metrics.

Computer aided detection CADe of breast cancer is mainly focused on monomodal applications. We propose an automated multimodal CADe approach, which uses patient-specific image registration of MRI and X-ray mammography to estimate the spatial correspondence of tissue structures. Then, based on the spatial correspondence, features are extracted from both MRI and X-ray mammography. As proof of principle, distinct regions of interest ROI were classified into normal and suspect tissue. We investigated the performance of different classifiers, compare our combined approach against a classification with MRI features only and evaluate the influence of the registration error. Using the multimodal information, the sensitivity for detecting suspect ROIs improved by 7i¾?% compared to MRI-only detection. The registration error influences the results: using only datasets with a registration error below $$10\,mm$$, the sensitivity for the multimodal detection increases by 10i¾?% to a maximum of 88i¾?%, while the specificity remains constant. We conclude that automatically combining MRI and X-ray can enhance the result of a CADe system.