| Home | << 1 2 >> |
|
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.
|
|
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.
|
|
Albiol, F., Corbi, A., & Albiol, A. (2017). 3D measurements in conventional X-ray imaging with RGB-D sensors. Med. Eng. Phys., 42, 73–79.
Abstract: A method for deriving 3D internal information in conventional X-ray settings is presented. It is based on the combination of a pair of radiographs from a patient and it avoids the use of X-ray-opaque fiducials and external reference structures. To achieve this goal, we augment an ordinary X-ray device with a consumer RGB-D camera. The patient' s rotation around the craniocaudal axis is tracked relative to this camera thanks to the depth information provided and the application of a modern surface-mapping algorithm. The measured spatial information is then translated to the reference frame of the X-ray imaging system. By using the intrinsic parameters of the diagnostic equipment, epipolar geometry, and X-ray images of the patient at different angles, 3D internal positions can be obtained. Both the RGB-D and Xray instruments are first geometrically calibrated to find their joint spatial transformation. The proposed method is applied to three rotating phantoms. The first two consist of an anthropomorphic head and a torso, which are filled with spherical lead bearings at precise locations. The third one is made of simple foam and has metal needles of several known lengths embedded in it. The results show that it is possible to resolve anatomical positions and lengths with a millimetric level of precision. With the proposed approach, internal 3D reconstructed coordinates and distances can be provided to the physician. It also contributes to reducing the invasiveness of ordinary X-ray environments and can replace other types of clinical explorations that are mainly aimed at measuring or geometrically relating elements that are present inside the patient's body.
|
|
Muñoz, E., Barrientos, L., Bernabeu, J., Borja-Lloret, M., Llosa, G., Ros, A., et al. (2020). A spectral reconstruction algorithm for two-plane Compton cameras. Phys. Med. Biol., 65(2), 025011–17pp.
Abstract: One factor limiting the current applicability extent of hadron therapy is the lack of a reliable method for real time treatment monitoring. The use of Compton imaging systems as monitors requires the correct reconstruction of the distribution of prompt gamma productions during patient irradiation. In order to extract the maximum information from all the measurable events, we implemented a spectral reconstruction method that assigns to all events a probability of being either partial or total energy depositions. The method, implemented in a list-mode maximum likelihood expectation maximization algorithm, generates a four dimensional image in the joint spatial-spectral domain, in which the voxels containing the emission positions and energies are obtained. The analytical model used for the system response function is also employed to derive an analytical expression for the sensitivity, which is calculated via Monte Carlo integration. The performance of the method is evaluated through reconstruction of various experimental and simulated sources with different spatial and energy distributions. The results show that the proposed method can recover the spectral and spatial information simultaneously, but only under the assumption of ideal measurements. The analysis of the Monte Carlo simulations has led to the identification of two important degradation sources: the mispositioning of the gamma interaction point and the missing energy recorded in the interaction. Both factors are related to the high energy transferred to the recoil electrons, which can travel far from the interaction point and even escape the detector. These effects prevent the direct application of the current method in more realistic scenarios. Nevertheless, experimental point-like sources have been accurately reconstructed and the spatial distributions and spectral emission of complex simulated phantoms can be identified.
Keywords: Compton imaging; Compton camera; hadron therapy; image reconstruction
|
|
Roser, J., Muñoz, E., Barrientos, L., Barrio, J., Bernabeu, J., Borja-Lloret, M., et al. (2020). Image reconstruction for a multi-layer Compton telescope: an analytical model for three interaction events. Phys. Med. Biol., 65(14), 145005–17pp.
Abstract: Compton Cameras are electronically collimated photon imagers suitable for sub-MeV to few MeV gamma-ray detection. Such features are desirable to enablein vivorange verification in hadron therapy, through the detection of secondary Prompt Gammas. A major concern with this technique is the poor image quality obtained when the incoming gamma-ray energy is unknown. Compton Cameras with more than two detector planes (multi-layer Compton Cameras) have been proposed as a solution, given that these devices incorporate more signal sequences of interactions to the conventional two interaction events. In particular, three interaction events convey more spectral information as they allow inferring directly the incident gamma-ray energy. A three-layer Compton Telescope based on continuous Lanthanum (III) Bromide crystals coupled to Silicon Photomultipliers is being developed at the IRIS group of IFIC-Valencia. In a previous work we proposed a spectral reconstruction algorithm for two interaction events based on an analytical model for the formation of the signal. To fully exploit the capabilities of our prototype, we present here an extension of the model for three interaction events. Analytical expressions of the sensitivity and the System Matrix are derived and validated against Monte Carlo simulations. Implemented in a List Mode Maximum Likelihood Expectation Maximization algorithm, the proposed model allows us to obtain four-dimensional (energy and position) images by using exclusively three interaction events. We are able to recover the correct spectrum and spatial distribution of gamma-ray sources when ideal data are employed. However, the uncertainties associated to experimental measurements result in a degradation when real data from complex structures are employed. Incorrect estimation of the incident gamma-ray interaction positions, and missing deposited energy associated with escaping secondaries, have been identified as the causes of such degradation by means of a detailed Monte Carlo study. As expected, our current experimental resolution and efficiency to three interaction events prevents us from correctly recovering complex structures of radioactive sources. However, given the better spectral information conveyed by three interaction events, we expect an improvement of the image quality of conventional Compton imaging when including such events. In this regard, future development includes the incorporation of the model assessed in this work to the two interaction events model in order to allow using simultaneously two and three interaction events in the image reconstruction.
|
|
Ros Garcia, A., Barrio, J., Etxebeste, A., Garcia-Lopez, J., Jimenez-Ramos, M. C., Lacasta, C., et al. (2020). MACACO II test-beam with high energy photons. Phys. Med. Biol., 65(24), 245027–12pp.
Abstract: The IRIS group at IFIC Valencia is developing a three-layer Compton camera for treatment monitoring in proton therapy. The system is composed of three detector planes, each made of a LaBr3<i monolithic crystal coupled to a SiPM array. Having obtained successful results with the first prototype (MACACO) that demonstrated the feasibility of the proposed technology, a second prototype (MACACO II) with improved performance has been developed, and is the subject of this work. The new system has an enhanced detector energy resolution which translates into a higher spatial resolution of the telescope. The image reconstruction method has also been improved with an accurate model of the sensitivity matrix. The device has been tested with high energy photons at the National Accelerator Centre (CNA, Seville). The tests involved a proton beam of 18 MeV impinging on a graphite target, to produce 4.4 MeV photons. Data were taken at different system positions of the telescope with the first detector at 65 and 160 mm from the target, and at different beam intensities. The measurements allowed successful reconstruction of the photon emission distribution at two target positions separated by 5 mm in different telescope configurations. This result was obtained both with data recorded in the first and second telescope planes (two interaction events) and, for the first time in beam experiments, with data recorded in the three planes (three interaction events).
Keywords: Compton imaging; Compton camera; proton therapy; LaBr3; test-beam; image reconstruction
|
|
Roser, J., Barrientos, L., Bernabeu, J., Borja-Lloret, M., Muñoz, E., Ros, A., et al. (2022). Joint image reconstruction algorithm in Compton cameras. Phys. Med. Biol., 67(15), 155009–15pp.
Abstract: Objective. To demonstrate the benefits of using an joint image reconstruction algorithm based on the List Mode Maximum Likelihood Expectation Maximization that combines events measured in different channels of information of a Compton camera. Approach. Both simulations and experimental data are employed to show the algorithm performance. Main results. The obtained joint images present improved image quality and yield better estimates of displacements of high-energy gamma-ray emitting sources. The algorithm also provides images that are more stable than any individual channel against the noisy convergence that characterizes Maximum Likelihood based algorithms. Significance. The joint reconstruction algorithm can improve the quality and robustness of Compton camera images. It also has high versatility, as it can be easily adapted to any Compton camera geometry. It is thus expected to represent an important step in the optimization of Compton camera imaging.
|