Bouhova-Thacker, E., Kostyukhin, V., Koffas, T., Liebig, W., Limper, M., Piacquadio, G. N., et al. (2010). Expected Performance of Vertex Reconstruction in the ATLAS Experiment at the LHC. IEEE Trans. Nucl. Sci., 57(2), 760–767.
Abstract: In the harsh environment of the Large Hadron Collider at CERN (design luminosity of 10(34) cm(-2) s(-1)) efficient reconstruction of vertices is crucial for many physics analyses. Described in this paper is the expected performance of the vertex reconstruction used in the ATLAS experiment. The algorithms for the reconstruction of primary and secondary vertices as well as for finding photon conversions and vertex reconstruction in jets are described. The implementation of vertex algorithms which follows a very modular design based on object-oriented C++ is presented. A user-friendly concept allows event reconstruction and physics analyses to compare and optimize their choice among different vertex reconstruction strategies. The performance of implemented algorithms has been studied on a variety of Monte Carlo samples and results are presented.
|
Rodriguez-Alvarez, M. J., Sanchez, F., Soriano, A., & Iborra, A. (2010). Sparse Givens resolution of large system of linear equations: Applications to image reconstruction. Math. Comput. Model., 52(7-8), 1258–1264.
Abstract: In medicine, computed tomographic images are reconstructed from a large number of measurements of X-ray transmission through the patient (projection data). The mathematical model used to describe a computed tomography device is a large system of linear equations of the form AX = B. In this paper we propose the QR decomposition as a direct method to solve the linear system. QR decomposition can be a large computational procedure. However, once it has been calculated for a specific system, matrices Q and R are stored and used for any acquired projection on that system. Implementation of the QR decomposition in order to take more advantage of the sparsity of the system matrix is discussed.
|
Aguiar, P., Rafecas, M., Ortuño, J. E., Kontaxakis, G., Santos, A., Pavia, J., et al. (2010). Geometrical and Monte Carlo projectors in 3D PET reconstruction. Med. Phys., 37(11), 5691–5702.
Abstract: Purpose: In the present work, the authors compare geometrical and Monte Carlo projectors in detail. The geometrical projectors considered were the conventional geometrical Siddon ray-tracer (S-RT) and the orthogonal distance-based ray-tracer (OD-RT), based on computing the orthogonal distance from the center of image voxel to the line-of-response. A comparison of these geometrical projectors was performed using different point spread function (PSF) models. The Monte Carlo-based method under consideration involves an extensive model of the system response matrix based on Monte Carlo simulations and is computed off-line and stored on disk. Methods: Comparisons were performed using simulated and experimental data of the commercial small animal PET scanner rPET. Results: The results demonstrate that the orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions yield better images in terms of contrast and spatial resolution than those obtained after using the conventional method and the multiray-based S-RT. Furthermore, the Monte Carlo-based method yields slight improvements in terms of contrast and spatial resolution with respect to these geometrical projectors. Conclusions: The orthogonal distance-based ray-tracer and Siddon ray-tracer using PSF image-space convolutions represent satisfactory alternatives to factorizing the system matrix or to the conventional on-the-fly ray-tracing methods for list-mode reconstruction, where an extensive modeling based on Monte Carlo simulations is unfeasible.
|
Blume, M., Martinez-Moller, A., Keil, A., Navab, N., & Rafecas, M. (2010). Joint Reconstruction of Image and Motion in Gated Positron Emission Tomography. IEEE Trans. Med. Imaging, 29(11), 1892–1906.
Abstract: We present a novel intrinsic method for joint reconstruction of both image and motion in positron emission tomography (PET). Intrinsic motion compensation methods exclusively work on the measured data, without any external motion measurements. Most of these methods separate image from motion estimation: They use deformable image registration/optical flow techniques in order to estimate the motion from individually reconstructed gates. Then, the image is estimated based on this motion information. With these methods, a main problem lies in the motion estimation step, which is based on the noisy gated frames. The more noise is present, the more inaccurate the image registration becomes. As we show both visually and quantitatively, joint reconstruction using a simple deformation field motion model can compete with state-of-the-art image registration methods which use robust multilevel B-spline motion models.
|
ANTARES Collaboration(Aguilar, J. A. et al), Bigongiari, C., Dornic, D., Emanuele, U., Gomez-Gonzalez, J. P., Hernandez-Rey, J. J., et al. (2011). A fast algorithm for muon track reconstruction and its application to the ANTARES neutrino telescope. Astropart Phys., 34(9), 652–662.
Abstract: An algorithm is presented, that provides a fast and robust reconstruction of neutrino induced upward-going muons and a discrimination of these events from downward-going atmospheric muon background in data collected by the ANTARES neutrino telescope. The algorithm consists of a hit merging and hit selection procedure followed by fitting steps for a track hypothesis and a point-like light source. It is particularly well-suited for real time applications such as online monitoring and fast triggering of optical follow-up observations for multi-messenger studies. The performance of the algorithm is evaluated with Monte Carlo simulations and various distributions are compared with that obtained in ANTARES data.
|
ANTARES Collaboration(Aguilar, J. A. et al), Bigongiari, C., Dornic, D., Emanuele, U., Gomez-Gonzalez, J. P., Hernandez-Rey, J. J., et al. (2012). A method for detection of muon induced electromagnetic showers with the ANTARES detector. Nucl. Instrum. Methods Phys. Res. A, 675, 56–62.
Abstract: The primary aim of ANTARES is neutrino astronomy with upward going muons created in charged current muon neutrino interactions in the detector and its surroundings. Downward going muons are background for neutrino searches. These muons are the decay products of cosmic-ray collisions in the Earth's atmosphere far above the detector. This paper presents a method to identify and count electromagnetic showers induced along atmospheric muon tracks with the ANTARES detector. The method is applied to both cosmic muon data and simulations and its applicability to the reconstruction of muon event energies is demonstrated.
|
Grkovski, M., Brzezinski, K., Cindro, V., Clinthorne, N. H., Kagan, H., Lacasta, C., et al. (2015). Evaluation of a high resolution silicon PET insert module. Nucl. Instrum. Methods Phys. Res. A, 788, 86–94.
Abstract: Conventional PET systems can be augmented with additional detectors placed in close proximity of the region of interest. We developed a high resolution PET insert module to evaluate the added benefit of such a combination. The insert module consists of two back-to-back 1 mm thick silicon sensors, each segmented into 1040 1 mm(2) pads arranged in a 40 by 26 array. A set of 16 VATAGP7.1 ASICs and a custom assembled data acquisition board were used to read out the signal from the insert module. Data were acquired in slice (20) geometry with a Jaszczak phantom (rod diameters of 12-4.8 mm) Filled with F-18-FDG and the images were reconstructed with ML-EM method. Both data with full and limited angular coverage from the insert module were considered and three types of coincidence events were combined. The ratio of high-resolution data that substantially improves quality of the reconstructed image for the region near the surface of the insert module was estimated to be about 4%. Results from our previous studies suggest that such ratio could be achieved at a moderate technological expense by using an equivalent of two insert modules (an effective sensor thickness of 4 mm).
|
Boronat, M., Fuster, J., Garcia, I., Ros, E., & Vos, M. (2015). A robust jet reconstruction algorithm for high-energy lepton colliders. Phys. Lett. B, 750, 95–99.
Abstract: We propose a new sequential jet reconstruction algorithm for future lepton colliders at the energy frontier. The Valencia algorithm combines the natural distance criterion for lepton colliders with the greater robustness against backgrounds of algorithms adapted to hadron colliders. Results on a detailed Monte Carlo simulation of t (t) over tilde and ZZ production at future linear e(+)e(-) colliders (ILC and CLIC) with a realistic level of background overlaid, show that it achieves better performance in the presence of background than the classical algorithms used at previous e(+)e(-) colliders.
|
Richart, J., Otal, A., Rodriguez, S., Nicolas, A. I., DePiaggio, M., Santos, M., et al. (2015). A practical MRI-based reconstruction method for a new endocavitary and interstitial gynaecological template. J. Contemp. Brachytherapy, 7(5), 407–414.
Abstract: Purpose: There are perineal templates for interstitial implants such as MUPIT and Syed applicators. Their limitations are the intracavitary component deficit and the necessity to use computed tomography (CT) for treatment planning since both applicators are non-magnetic resonance imaging (MRI) compatibles. To overcome these problems, a new template named Template Benidorm (TB) has been recently developed. Titanium needles are usually reconstructed based on their own artifacts, mainly in T1-weighted sequence, using the void on the tip as the needle tip position. Nevertheless, patient tissues surrounding the needles present heterogeneities that complicate the accurate identification of these artifact patterns. The purpose of this work is to improve the titanium needle reconstruction uncertainty for the TB case using a simple method based on the free needle lengths and typical MRI pellets markers. Material and methods: The proposed procedure consists on the inclusion of three small A-vitamin pellets (hyper-intense on MRI images) compressed by both applicator plates defining the central plane of the plate's arrangement. The needles used are typically 20 cm in length. For each needle, two points are selected defining the straight line. From such line and the plane equations, the intersection can be obtained, and using the free length (knowing the offset distance), the coordinates of the needle tip can be obtained. The method is applied in both T1W and T2W acquisition sequences. To evaluate the inter-observer variation of the method, three implants of T1W and another three of T2W have been reconstructed by two different medical physicists with experience on these reconstructions. Results and conclusions: The differences observed in the positioning were significantly smaller than 1 mm in all cases. The presented algorithm also allows the use of only T2W sequence either for contouring or reconstruction purposes. The proposed method is robust and independent of the visibility of the artifact at the tip of the needle.
|
Albiol, F., Corbi, A., & Albiol, A. (2016). Geometrical Calibration of X-Ray Imaging With RGB Cameras for 3D Reconstruction. IEEE Trans. Med. Imaging, 35(8), 1952–1961.
Abstract: We present a methodology to recover the geometrical calibration of conventional X-ray settings with the help of an ordinary video camera and visible fiducials that are present in the scene. After calibration, equivalent points of interest can be easily identifiable with the help of the epipolar geometry. The same procedure also allows the measurement of real anatomic lengths and angles and obtains accurate 3D locations from image points. Our approach completely eliminates the need for X-ray-opaque reference marks (and necessary supporting frames) which can sometimes be invasive for the patient, occlude the radiographic picture, and end up projected outside the imaging sensor area in oblique protocols. Two possible frameworks are envisioned: a spatially shifting X-ray anode around the patient/object and a moving patient that moves/rotates while the imaging system remains fixed. As a proof of concept, experiences with a device under test (DUT), an anthropomorphic phantom and a real brachytherapy session have been carried out. The results show that it is possible to identify common points with a proper level of accuracy and retrieve three-dimensional locations, lengths and shapes with a millimetric level of precision. The presented approach is simple and compatible with both current and legacy widespread diagnostic X-ray imaging deployments and it can represent a good and inexpensive alternative to other radiological modalities like CT.
|
