USCT Data Challenge PDF
USCT Data Challenge PDF
net/publication/315003782
CITATIONS READS
7 86
5 authors, including:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Hartmut Gemmeke on 20 July 2018.
ABSTRACT
In the past years we have perceived within the USCT research community a demand for freely available USCT data sets.
Inspired by the idea of Open Science, this collection of data sets could stimulate the collaboration and the exchange of
ideas and experiences between USCT researchers. In addition, it may lead to comprehensive comparison of different
reconstruction algorithms and their results. Finally, by collecting feedback from the users about data and system
architecture, valuable information is gathered for further development of measurement setups. For the above reasons, we
have initiated a digital portal with several reference data sets and access scripts under free licenses. To kick off this
initiative, we organized a USCT data challenge event at SPIE Medical Imaging 2017.
Keywords: Ultrasound computer tomography, data challenge, open science, open data
1. INTRODUCTION
Ultrasound Computer Tomography is an exciting new technology mostly aimed at breast cancer imaging. Due to the
complex interaction of ultrasound with human tissue, the large amount of raw data, and the large volumes of interest,
image reconstruction is challenging both from a physical and a computational point of view. Many different image
reconstruction methods can be applied, ranging from simple ray based methods to full inversion with very different
approximations of the linear acoustic wave equation, see e.g. [1]. While the ray based algorithms offer fast
reconstruction, the inversion algorithms promise high image quality and resolution.
In recent years various groups have proposed a number of different approaches for image reconstruction, e.g. [1-23].
Direct comparison and analysis of the algorithms is difficult as they are usually applied to simulated data or – if available
– to real data acquired with very different measurement setups, e.g. [24-28]. The relatively few experimental data sets are
usually not commonly available. Moreover, even if they would be readily available, using the different date sets to test a
new reconstruction method would require accurate (and often not freely available) knowledge about the experimental
setup and data handling.
To mitigate such challenges for the scientific communities, there is an ongoing endeavor for so-called “Open science”.
Several leading European academic institutions signed in 2003 the “Berlin Declaration” on “open access to knowledge”
to enable easier exchange, participation and collaboration within the scientific community [33-35].
Following this idea and to mend gaps within the US(C)T community, the establishment of open and easy to use data and
code interfaces is the long term goal of this initiative. To stimulate the exchange of the available reconstruction
algorithms and raw data sets of different USCT devices, a reference database with freely available and open licensed
USCT data for comparison of reconstruction algorithms will be established, maintained and updated. Additionally, the
feedback about data and system architecture of the scientists working on reconstruction methods will help to drive
further development of the various measurement setups.
2. DESCRIPTION OF CHALLENGE
This challenge aims on applying available image reconstruction algorithms on provided USCT data in order to establish
a first intercommunication and standards for open data interface. In addition, this challenge marks the start of the USCT
reference database. The data sets and software for data access are available via the USCT challenge home page and the
linked data and code repository [29]. The long term goals of this work is to build up a free and open licensed reference
Medical Imaging 2017: Ultrasonic Imaging and Tomography, edited by Neb Duric, Brecht Heyde, Proc. of
SPIE Vol. 10139, 101391N · © 2017 SPIE · CCC code: 1605-7422/17/$18 · doi: 10.1117/12.2272593
The materialss are providedd using a free and open liceense, i.e. the BSDB 3-clause license [30] ffor code and data, allowingg
free use and publication
p off results. Afterr the challengee the data setss may remain with
w the particcipants, and can
c be used forr
additional ressearch and othher purposes. Further publiication and reesults should follow
f good aacademic pracctices and citee
this paper.
3. USC
CT SYSTEM OF
O TU DELF
FT
The Delft Brreast Ultrasouund Scanner (DBUS) is deppicted on the left-hand sidde of Figure 11. The system m consists of a
water tank with
w dimensionns 0.75 m x 0.75 0 m x 0.655 m with a waater level of 0.45 0 m. The ttemperature ofo the water iss
continuously monitored viaa thermocouples and kept constant
c n 1oC using heating
within h mats annd a temperatture controllerr.
On top of thee system, two rotary stages (LG Motion LGR1090-PD D), controlled by motor drivvers (Parker),, are mountedd.
The first rotaary stage rotaates the objectt, the second the receiver (0.5 MHz, Paanametrics V3318). The sou urce, which iss
identical to thhe receiver, is mounted at a fixed positionn in the cornerr of the tank. In
I this way a 22D tomograph hic scan of thee
object is obtaained. For eacch A-scan, ann electric pulsse with a centter frequency of approx. 00.5 MHz is geenerated by ann
arbitrary wavve form generrator (Agilent 33250A), am mplified (Electtronics and Innnovation 2100L 40 dB), and d successivelyy
damped usingg a variable atttenuator (JFW W Industries, 50BR-036).
5 The
T resulting wavew field is m
measured usin
ng an identicaal
transducer coonnected to a 14 bit A/D converter (S Spectrum M3i.4142-exp - PCI) which is set to a sample rate of
400 MHz. Foor each A-scaan, in additionn to the raw (unfiltered)
( R data, information such as temperaturre, source andd
RF
receiver posittions, etc. is stored.
s Togethher with a scann including ann object, a scaan in absence of an object, referred to ass
an empty scann, is made forr reference. An
A example off an empty scaan is shown inn Figure 2. Thhe A-scan, preesented in timee
and frequencyy domain, corrresponds to a measurementt with the sourrce and receivver facing eachh other.
4. USCT
U SYSTE
EM OF KIT
The KIT’s 3DD Ultrasound Computer Toomography sysstem (3D USC CT) is depicteed on the rightt hand side off Figure 1. Thee
device has a semi-ellipsoiddal 3D apertuure. Approx. sppherical wavee fronts are geenerated by eaach emitter att 2.5 MHz andd
-10
0.2
-20
Ñ n4
E -30
.- 0.6
-40
0.8
-50
-60
50 100 150 200 250 300 350
[o]
5 -40
1 ô -50
ó
co
0.5 -60
ii
ii 0 o.
-70
-0.5 -80
o
1 ti -90
-1.5 -100
0.25 0.3 0.35 0 0.2 0.4 0.6 0.8 1
t [ms] f [MHz]
Figure 2: The RF data measured with the DBUS system in absence of an object; (top) B-scan showing the field measured as a
function of time and receiver angle for single source position; (bottom-left and -right) single A-scan in time and frequency domain.
with a bandwidth of 1.5 MHz at -6 dB. The semi-elliptical aperture has a diameter of 26 cm and a height of 16 cm.
Rotational and translational movements, so-called aperture positions, of the complete sensor system create additional
virtual positions of the transducers.
The 2041 individual transducers are either operated as emitter (628) or receiver (1413). The transducers have opening
angles of 38.2° (standard deviation ±1.5°) at -6 dB. Four emitters and nine receivers are grouped together including pre-
amplifier and control electronics in so-called Transducer Array Systems (TAS). Each of the 157 TAS contains a
temperature sensor for tracking the temperature distribution within the water basin and the shift at each TAS position
during measurements. Additionally, two calibrated PT100 temperature sensors are embedded in the TAS holder to
enable increased accuracy.
The data acquisition is carried out with an FPGA based system, which can store up to 80 of A-scans [31]. The
digitalization is performed by 480 parallel channels (12 at 20 ), enabling data acquisition at one aperture
position in approx. ten seconds. After digitization, the parallel data streams are processed as follows: First, the data
streams are bandpass filtered (1.67 to 3.33 at −60 ). Next, the data rate is reduced by a factor of six by
performing bandpass undersampling. Finally, the reduced data is stored in the internal memory buffer. Using this
approach up to 47 data sets at different aperture positions can be stored in one data acquisition step. A detailed
description of the 3D USCT system can be found in [32].
The emitters are excited with a coded excitation signal, e.g. frequency coded chirps can be applied to increase the signal-
to-noise ratio of the data. Also the gain of the receiving channels is set individually based on an initial measurement. The
applied coded excitation, the individual gain, the temperature data and the spatial positions of the aperture are stored along
with the A-scans for each measurement and can be used for signal (pre-) processing and image reconstruction. Empty
measurements are also provided and exemplary data is depicted in Figure 3.
0.5
time/s
1.5
2.5
0.8
1
0.6
0.4
0.5
0.2
Amplitude
p(t)
0 0
-0.2
-0.5 -0.4
-0.6
-1
-0.8
Figure 3: RF data measured with the 3D USCT system in absence of an object; (top) B-scan of the raw data of one emitter in the first
TAS of top transducer row of aperture as a function of time and one receiver of each TAS counter-clockwise from top to bottom;
bottom-left shows a single A-scan (receiver opposite of emitter element) in time and bottom-right a frequency coded chirp applied
for coded excitation.
5. DATASETS
In total eight data sets have been made available; two sets are obtained with the system from the TU Delft, and six with
the system from KIT. The data is available under “Open Data Commons Attribution License” and the access scripts
under the 3-clause BSD license.
A compact overview of the available data sets is presented in Table 1 and Figure 4 shows photos of the phantoms and
exemplary reconstructed slices.
Figure 4: Phanntoms scanned; (left) photo, (rright) reconstruuction using different reconstrruction algorithhms. From top left to bottom
right: TU Delfft: agar based phantom
p reconsstructed with Delay
D and Sum;; KIT: gelatin 3 phantom recoonstructed with h SAFT in low
resolution and displayed as maximum
m intenssity projection to
t enhance the visibility of the bottom of thee plastic cup, tu
urkey phantom
reconstructed with
w SAFT, andd finally nylon thread
t phantomm with speed of sound correctedd SAFT reconstruction.
Gelatin 3: Thhe phantom coonsists of a geelatin phantomm with diameeter of approx. 0.07 m at thhe bottom and d 0.10 m at thee
top, a height of approx. 0.113 m of whichh approx. 0.100 m where immmersed into thhe USCT aperrture. The speeed of sound of
the gelatin was
w approx. 15151 m/s. Thee gelatin wass embedded in i a plastic cup.
c Two inclusions were generated byy
embedding drrinking strawss with diameteer 5 mm in thee gelatin durin
ng curing. Prioor to scanningg the object, th
he straws weree
removed andd the inclusionns were filledd with water. The phantom m was positiooned approx. centrally in the 3D USCT T
aperture.
6. FUTURE
E PROSPECTIIVE FOR USC
CT REFERE
ENCE DATA BASE
The aim of thhis work is to start a freely available USC CT data base enabling to teest and analyzze reconstructiion algorithmss
or to comparre the perform mance of com mpeting algorithms. In a first step inittial data basee entries are provided andd
participants of
o this challennge are asked tot report theirr experiences in using the data.
d Additionnally a panel discussion
d willl
After the kick off of opening the data base for the public and establishing easy to use data access further challenges are
planned, e.g. challenging the different algorithms by comparison of obtained image quality or computational
performance. The data sets are freely available and open licensed, so they can also be used outside of challenges for
evaluation of advanced reconstruction techniques on real data, further development of algorithms for image
reconstruction and signal processing.
The feedback about data and USCT systems resulting from the use with different algorithms can also lead to drive
further development of the system architecture, e.g. to research the optimization of transducer positioning or limits on
signal-to-noise levels.
We hope also that these challenges will lead to a growth of the data base by increasing the number of imaged objects
and/or data acquired with different USCT architectures.
Finally, the overall motivation for all of our work is to support early breast cancer diagnosis. Thus, we also need to learn
more about the properties of breast cancer tissues and how to image them to make the most distinguishing properties
available for diagnosis. This upcoming knowledge needs to be included into the data base by imaging and if possible
quantify the properties of appropriate phantoms and tissue structures.
The authors would like to acknowledge the following people: Bastiaan Dekker and Elmer Jansen for building the 2D
setup and performing the measurements; Clemens Kaiser, Wei Yap Tan and Carina Wittenbeck for building and
measuring phantoms for the 3D USCT and the many colleagues in the institutes who made the built up of USCT systems
possible.
REFERENCES
[1] Ozmen, N.; Dapp, R.; Zapf, M. et al.: “Comparing different ultrasound imaging methods for breast cancer
detection,” IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, 62(4), pp. 637–46 (2015).
[2] Natterer, F.: “Acoustic mammography in the time domain,” Technical report, University Muenster, (2008).
[3] Zhang, X.; Broschat, S.L.; and Flynn, P.J.: “A numerical study of conjugate gradient directions for an ultrasound
inverse problem,” Journal of Computational Acoustics, 12(4), pp. 587-604 (2004).
[4] Goncharsky, A.V.; and Romanov, S.Y.: “Supercomputer technologies in inverse problems of ultrasound
tomography,” Inverse Problems, 29:075004, (2013).
[5] Sandhu, G.Y.; Li; Roy, O.; Schmidt, S.; and Duric, N.: “Frequency domain ultrasound waveform tomography:
Breast imaging using a ring transducer,” Physics in Medicine and Biology, 60(14), pp. 5381-98 (2015).
[6] Wiskin, J.; Borup, D.; Johnson, S.; et al.: “Three dimensional nonlinear inverse scattering: Quantitative transmission
algorithms, refraction corrected reflection, scanner design and clinical results,” Proc. Meetings on Acoustics, 19,
(2013).
[7] Chenevert, T.L.; Bylski, D.I.; Carson, P.L.; et al.: “Ultrasonic computed tomography of the breast. Improvement of
image quality by use of cross-correlation time-of-flight and phase-insensitive attenuation measurements,”
Radiology, 152 (1), pp. 155-9 (1984).
[8] Kaveh, M.; Mueller, R.K.; Rylander, R.; et al.: “Experimental results in ultrasonic diffraction tomography,”
Acoustical Imaging (K.Y. Wang, Ed.), 9, pp. 433-50 (1980).
[9] Simonetti, F.; Huang, L.; and Duric, N.: “A multiscale approach to diffraction tomography of complex three-
dimensional objects,” Applied Physics Letters, 95, (2009).
[10] Sponheim, N.; Gelius, L.-J.; Johansen, I.: and Stamnes, J.J.: “Quantitative results in ultrasonic tomography of large
objects using line sources and curved detector arrays,” IEEE Transactions on Ultrasonics, Ferroelectrics, and
Frequency Control, 38(4), pp. 370-9 (1991).
[11] Greenleaf, J.F.; Ylitalo, J.; and Gisvold, J.J.: “Ultrasonic computed tomography for breast examination,” IEEE
Engineering in Medicine and Biology Magazine: The Quarterly Magazine of the Engineering in Medicine &
Biology Society, 6(4), pp. 27-32 (1987).
[12] Andre, M.P.; Janee, H.S.; Martin, P.J.; et al.: “High speed data acquisition in a diffraction tomography system
employing large-scale toroidal arrays,” International Journal of Imaging Systems and Technology, 8(1), pp. 137-47
(1997).
[13] Lavarello, R.J.; and Oelze, M.L.: “Tomographic reconstruction of three-dimensional volumes using the distorted
Born iterative method,” IEEE Transactions on Medical Imaging, 28(10), pp. 1643-53 (2009).
[14] Andersen, A.H.: “A ray tracing approach to restoration and resolution enhancement in experimental ultrasound
tomography,” Ultrasonic Imaging, 12(4), pp. 268-91 (1990).
[15] Hormati, A.; Jovanovic, I.; Roy, O.; and Vetterli, M.: “Robust ultrasound travel-time tomography using the bent ray
model” Proc. SPIE 7629, Medical Imaging 2010: Ultrasonic Imaging, Tomography, and Therapy, 7629, (2010).
[16] Li, C.; Duric, N.; and Huang, L.: “Breast imaging using transmission ultrasound: Reconstructing tissue parameters
of sound speed and attenuation,” BioMedical Engineering and Informatics: New Development and the Future -
Proceedings of the 1st International Conference on BioMedical Engineering and Informatics, 2, pp. 708-12 (2008).
[17] Dapp; R.: “Abbildungsmethoden für die Brust mit einem 3D-Ultraschall-Computertomographen,” PhD thesis,
Karlsruhe Institute of Technology, (2013).
[18] Kretzek, E.; and Ruiter, N.V.: “GPU based 3D SAFT reconstruction including phase aberration,” Proc. SPIE 9040,
Medical Imaging 2014: Ultrasonic Imaging and Tomography, p. 90400W (2014).
[19] Kretzek, E.; Hopp, T.; and Ruiter, N.V.: “GPU-based 3D SAFT reconstruction including attenuation correction,”
Proc. SPIE 9419, Medical Imaging 2015: Ultrasonic Imaging and Tomography, p. 94190E, (2015).