189 research outputs found
Advanced capabilities for planar X-ray systems
Mención Internacional en el título de doctorThe past decades have seen a rapid evolution towards the use of digital detectors
in radiology and a more flexible robotized movement of the system components,
X-ray tube and detector. This evolution opened the possibility for incorporating
advanced capabilities in these planar X-ray systems, and for providing new valuable
diagnostic information compared to the previous technology. Some of the current
challenges for radiography are to obtain more quantitative images and to reduce the
inherent superposition of tissues because of the 2D nature of the technique.
Dual energy radiography, based on the acquisition of two images at different
source voltages, enables a separate characterization of soft tissue and bone structures.
Its benefits over conventional radiography have been proven in different applications,
since it improves information content without adding significant extra
acquisition time or radiation dose.
In a different direction, a really disruptive advance would be to obtain 3D imaging
with systems designed just for planar images. The incorporation of tomographic
capabilities into these systems would have to deal with the acquisition of a limited
number of projections, with non-standard geometrical configurations.
This thesis presents original contributions in these two directions: dual energy
radiography and 3D imaging with X-ray systems designed for planar imaging. The
work is framed in a line of research of the Biomedical Imaging and Instrumentation
Group from the Bioengineering and Aerospace Department of University Carlos III
de Madrid working jointly with the University Hospital Gregorio Marañón, focused
on the advance of radiology systems. This research line is carried out in collaboration
with the group of Computer Architecture, Communications and Systems (ARCOS),
from the same university, the Imaging Research Laboratory (IRL) of the University
of Washington and the research center CREATIS, France. The research has a clear
focus on technology transfer to the industry through the company Sedecal, a Spanish
multinational among the 10 best world companies in the medical imaging field.
The first contribution of this thesis is a complete novel protocol to incorporate
dual energy capabilities that enable quantitative planar studies. The proposal is
based on the use of a preliminary calibration with a very simple and low-cost phantom
formed by two parts that represent soft tissue and bone equivalent materials.
This calibration is performed automatically with no strict placement requirements.
Compared to current Dual-energy X-ray Absorptiometry (DXA) systems, 1) it provides
real mass-thickness values directly, enabling quantitative planar studies instead
of relative comparisons, and 2) it is based on an automatic preliminary calibration without the need of interaction of an experienced technician.
The second contribution is a novel protocol for the incorporation of tomographic
capabilities into X-ray systems originally intended for planar imaging. For this purpose,
we faced three main challenges.
First, the geometrical trajectory of equipment follows non-standard circular orbits,
thus posing severe difficulties for reconstruction. To handle this, the proposed
protocol comprises a new geometrical calibration procedure that estimates all the
system parameters per-projection.
Second, the reconstruction of a limited number of projections from a reduced angular
span leads to severe artifacts when using conventional reconstruction methods.
To deal with these limited-view data, the protocol includes a novel advanced reconstruction
method that incorporates the surface information of the sample, which
can be extracted with a 3D light surface scanner. These data are introduced as an
imposed constraint following the Split Bregman formulation. The restriction of the
search space by exploiting the surface-based support becomes crucial for a complete
recovery of the external contour of the sample and surroundings when the angular
span is extremely reduced. The modular, efficient and flexible design followed for its
implementation allows for the reconstruction of limited-view data with non-standard
trajectories.
Third, the optimization of the acquisition protocols has not yet explored with
these systems. This thesis includes a study of the optimum acquisition protocols
that allowed us to identify the possibilities and limitations of these planar systems.
Using the surface-constrained method, it is possible to reduce the total number of
projections up to 33% and the angular span down to 60 degrees.
The contributions of this thesis open the way to provide depth and quantitative
information very valuable for the improvement of radiological diagnosis. This could
impact considerably the clinical practice, where conventional radiology is still the
imaging modality most used, accounting for 80-90% of the total medical imaging
exams. These advances open the possibility of new clinical applications in scenarios
where 1) the reduction of the radiation dose is key, such as lung cancer screening or
Pediatrics, according to the ALARA criteria (As Low As Reasonably Achievable),
2) a CT system is not usable due to movement limitations, such as during surgery
or in an ICU and 3) where costs issues complicate the availability of CT systems,
such as rural areas or underdeveloped countries.
The results of this thesis has a clear application in the industry, since it is part
of a proof of concept of the new generation of planar X-ray systems that will be
commercialized worldwide by the company SEDECAL (Madrid, Spain).Los últimos años están viendo un rápido avance de los sistemas de radiología hacia el
uso de detectores digitales y a una mayor flexibilidad de movimientos de los principales
componentes del sistema, el tubo de rayos X y el detector. Esta evolución abre
la posibilidad de incorporar capacidades avanzadas en sistemas de imagen plana por
rayos X proporcionando nueva información valiosa para el diagnóstico. Dos retos en
radiografía son obtener imágenes cuantitativas y reducir la superposición de tejidos
debida a la naturaleza proyectiva de la técnica.
La radiografía de energía dual, basada en la adquisición de dos imágenes a diferente
kilovoltaje, permite obtener imágenes de tejido blando y hueso por separado.
Los beneficios de esta técnica que aumenta la cantidad de información sin añadir
un tiempo de adquisición o de dosis de radiación extra significativos frente al uso de
radiografía convencional, han sido demostrados en diferentes aplicaciones.
En otra dirección, un avance realmente disruptivo sería la obtención de imagen
3D con sistemas diseñados únicamente para imagen plana. La incorporación de capacidades
tomográficas en estos sistemas tendría que lidiar con la adquisición de un
número limitado de proyecciones siguiendo trayectorias no estándar.
Esta tesis presenta contribuciones originales en esas dos direcciones: radiografía
de energía dual e imagen 3D con sistemas de rayos X diseñados para imagen plana.
El trabajo se encuadra en una línea de investigación del grupo de Imagen Biomédica
e Instrumentación del Departamento de Bioingeniería e Ingeniería Aerospacial de
la Universidad Carlos III de Madrid junto con el Hospital Universitario Gregorio
Marañon, centrada en el avance de sistemas de radiología. Esta línea de investigación
se desarollada en colaboración con el grupo Computer Architecture, Communications
and Systems (ARCOS), de la misma universidad, el grupo Imaging Research Laboratory
(IRL) de la Universidad de Washington y el centro de investigación CREATIS,
de Francia. Se trata de una línea de investigación con un claro enfoque de transferencia
tecnológica a la industria a través de la compañía SEDECAL, una multinacional
española de entre las 10 líderes del mundo en el campo de la radiología.
La primera contribución de esta tesis es un protocolo completo para incorporar
capacidades de energía dual que permitan estudios cuantitativos de imagen plana.
La propuesta se basa en una calibración previa con un maniquí simple y de bajo coste
formado por dos materiales equivalentes de tejido blando y hueso respectivamente.
Comparado con los sistemas actuales DXA (Dual-energy X-ray Absorptiometry),
1) proporciona valores reales de tejido atravesado, 2) se basa en una calibración
automática que no requiere la interacción de un técnico con gran experiencia. La segunda contribución es un protocolo nuevo para la incorporación de capacidades
tomográficas en sistemas de rayos X originariamente diseñados para imagen
plana. Para ello, nos enfrentamos a tres principales dificultades.
En primer lugar, las trayectorias que pueden seguir la fuente y el detector en
estos sistemas no constituyen órbitas circulares estándares, lo que plantea retos importantes
en la caracterización geométrica. Para solventarlo, el protocolo propuesto
incluye una calibración geométrica que estima todos los parámetros geométricos del
sistema para cada proyección.
En segundo lugar, la reconstrucción de un número limitado de proyecciones
adquiridas en un rango angular reducido da lugar a artefactos graves cuando se
reconstruye con algoritmos convencionales. Para lidiar con estos datos de ángulo
limitado, el protocolo incluye un nuevo método avanzado de reconstrucción que incorpora
la información de superficie de la muestra, que se puede se obtener con un
escáner 3D. Esta información se impone como una restricción siguiendo la formulación
de Split Bregman, para compensar la falta de datos. La restricción del espacio
de búsqueda a través de la explotación del soporte basado en superficie, es crucial
para una recuperación completa del contorno externo de la muestra cuando el rango
angular es extremadamente pequeño. El diseño modular, eficiente y flexible de la
implementación propuesta permite reconstruir datos de ángulo limitado obtenidos
con posiciones de fuente y detector no estándar.
En tercer lugar, hasta la fecha, no se ha explorado la optimización del protocolo
de adquisición con estos sistemas. Esta tesis incluye un estudio de los protocolos
óptimos de adquisición que permitió identificar las posibilidades y limitaciones de
estos sistemas de imagen plana. Gracias al método de reconstrucción basado en
superficie, es posible reducir el número total de proyecciones hasta el 33% y el rango
angular hasta 60 grados.
Las contribuciones de esta tesis abren la posibilidad de proporcionar información
de profundidad y cuantitativa muy valiosa para la mejora del diagnóstico radiológico.
Esto podría impactar considerablemente en la práctica clínica, donde la radiología
convencional es todavía la modalidad de imagen más utilizada, abarcando el 80-
90% del total de los exámenes de imagen médica. Estos avances abren la posibilidad
de nuevas aplicaciones clínicas en escenarios donde 1) la reducción de la dosis de
radiación es clave, como en screening de cáncer de pulmón, de acuerdo con el criterio
ALARA (As Low As Reasonably Achievable), 2) no se puede usar un sistema
TAC por limitaciones de movimiento como en cirugía o UCI, o 3) el coste limita la
disponibilidad de sistemas TAC, como en zonas rurales o en países subdesarrollados.
Los resultados de esta tesis presentan una clara aplicación industrial, ya que
son parte de un prototipo de la nueva generación de sistemas planos de rayos X que
serán distribuidos mundialmente por la compañía SEDECAL.This thesis has been developed as part of several research projects with public funding:
- DPI2016-79075-R. ”Nuevos escenarios de tomografía por rayos X”, IP: Mónica
Abella García, Ministerio de Economía y Competitividad, 01/01/2017-31/12/2019,
147.620 e.
- ”Nuevos escenarios de tomografía por rayos X (NEXT) DPI2016-79075-R.
Ministerio de Economía”, Industria y Competitividad. (Universidad Carlos
III de Madrid). 30/12/2016-29/12/2019. 147.620 e.
(…)
- FP7-IMI-2012 (GA-115337), ”PreDict-TB: Model-based preclinical development
of anti-tuberculosis drug combinations”. FP7-IMI - Seventh Framework
Programme (EC-EFPIA). Unión Europea. (Universidad Carlos III de Madrid).
01/05/2012-31/10/2017.
(…)
- TEC2013-47270-R, ”Avances en Imagen Radiológica (AIR)”, Ministerio de
Economía y Competitividad”, 01/01/2014-31/12/2016. IP: Mónica Abella Garcia
and Manuel Desco Menéndez. 160.204 e
(…)
- RTC-2014-3028-1, ”Nuevos Escenarios Clínicos con Radiología Avanzada (NECRA)”,
Ministerio de Economía y Competitividad, 01/06/2014-31/12/2016 IP: Mónica
Abella García. 2014-2016. 219.458,96 e
- IDI-20130301, ”Nuevo sistema integral de radiografía (INNPROVE: INNovative
image PROcessing in medicine and VEterinary)”, IP: Mónica Abella García
and Manuel Desco Menéndez. Ministerio de Economía y Competitividad.
Subcontratación CDTI, 14/01/2013-31/03/2015. Total: 1.860.629e (UC3M:
325.000e). (Art. 83)
- IPT-2012-0401-300000 INNPACTO 2012, ”Tecnologías para Procedimientos
Intraoperatorios Seguros y Precisos. XIORT. MINECO. (Universidad Carlos
III de Madrid). 01/01/2013-31/12/2015.Programa Oficial de Doctorado en Ingeniería MatemáticaPresidente: Doménec Ros Puig.- Secretario: Cyril Riddell.- Vocal: Yannick Boursie
Enhancing Compressed Sensing 4D Photoacoustic Tomography by Simultaneous Motion Estimation
A crucial limitation of current high-resolution 3D photoacoustic tomography
(PAT) devices that employ sequential scanning is their long acquisition time.
In previous work, we demonstrated how to use compressed sensing techniques to
improve upon this: images with good spatial resolution and contrast can be
obtained from suitably sub-sampled PAT data acquired by novel acoustic scanning
systems if sparsity-constrained image reconstruction techniques such as total
variation regularization are used. Now, we show how a further increase of image
quality can be achieved for imaging dynamic processes in living tissue (4D
PAT). The key idea is to exploit the additional temporal redundancy of the data
by coupling the previously used spatial image reconstruction models with
sparsity-constrained motion estimation models. While simulated data from a
two-dimensional numerical phantom will be used to illustrate the main
properties of this recently developed
joint-image-reconstruction-and-motion-estimation framework, measured data from
a dynamic experimental phantom will also be used to demonstrate their potential
for challenging, large-scale, real-world, three-dimensional scenarios. The
latter only becomes feasible if a carefully designed combination of tailored
optimization schemes is employed, which we describe and examine in more detail
Compressed Sensing Based Reconstruction Algorithm for X-ray Dose Reduction in Synchrotron Source Micro Computed Tomography
Synchrotron computed tomography requires a large number of angular projections to reconstruct tomographic images with high resolution for detailed and accurate diagnosis. However, this exposes the specimen to a large amount of x-ray radiation. Furthermore, this increases scan time and, consequently, the likelihood of involuntary specimen movements. One approach for decreasing the total scan time and radiation dose is to reduce the number of projection views needed to reconstruct the images. However, the aliasing artifacts appearing in the image due to the reduced number of projection data, visibly degrade the image quality. According to the compressed sensing theory, a signal can be accurately reconstructed from highly undersampled data by solving an optimization problem, provided that the signal can be sparsely represented in a predefined transform domain. Therefore, this thesis is mainly concerned with designing compressed sensing-based reconstruction algorithms to suppress aliasing artifacts while preserving spatial resolution in the resulting reconstructed image. First, the reduced-view synchrotron computed tomography reconstruction is formulated as a total variation regularized compressed sensing problem. The Douglas-Rachford Splitting and the randomized Kaczmarz methods are utilized to solve the optimization problem of the compressed sensing formulation.
In contrast with the first part, where consistent simulated projection data are generated for image reconstruction, the reduced-view inconsistent real ex-vivo synchrotron absorption contrast micro computed tomography bone data are used in the second part. A gradient regularized compressed sensing problem is formulated, and the Douglas-Rachford Splitting and the preconditioned conjugate gradient methods are utilized to solve the optimization problem of the compressed sensing formulation. The wavelet image denoising algorithm is used as the post-processing algorithm to attenuate the unwanted staircase artifact generated by the reconstruction algorithm.
Finally, a noisy and highly reduced-view inconsistent real in-vivo synchrotron phase-contrast computed tomography bone data are used for image reconstruction. A combination of prior image constrained compressed sensing framework, and the wavelet regularization is formulated, and the Douglas-Rachford Splitting and the preconditioned conjugate gradient methods are utilized to solve the optimization problem of the compressed sensing formulation. The prior image constrained compressed sensing framework takes advantage of the prior image to promote the sparsity of the target image. It may lead to an unwanted staircase artifact when applied to noisy and texture images, so the wavelet regularization is used to attenuate the unwanted staircase artifact generated by the prior image constrained compressed sensing reconstruction algorithm.
The visual and quantitative performance assessments with the reduced-view simulated and real computed tomography data from canine prostate tissue, rat forelimb, and femoral cortical bone samples, show that the proposed algorithms have fewer artifacts and reconstruction errors than other conventional reconstruction algorithms at the same x-ray dose
Acceleration of Subtractive Non-contrast-enhanced Magnetic Resonance Angiography
Although contrast-enhanced magnetic resonance angiography (CE-MRA) is widely established as a clinical examination for the diagnosis of human vascular diseases, non-contrast-enhanced MRA (NCE-MRA) techniques have drawn increasing attention in recent years. NCE-MRA is based on the intrinsic physical properties of blood and does not require the injection of any exogenous contrast agents. Subtractive NCE-MRA is a class of techniques that acquires two image sets with different vascular signal intensity, which are later subtracted to generate angiograms.
The long acquisition time is an important drawback of NCE-MRA techniques, which not only limits the clinical acceptance of these techniques but also renders them sensitive to artefacts from patient motion. Another problem for subtractive NCE-MRA is the unwanted residual background signal caused by different static background signal levels on the two raw image sets. This thesis aims at improving subtractive NCE-MRA techniques by addressing both these limitations, with a particular focus on three-dimensional (3D) femoral artery fresh blood imaging (FBI).
The structure of the thesis is as follows:
Chapter 1 describes the anatomy and physiology of the vascular system, including the characteristics of arteries and veins, and the MR properties and flow characteristics of blood. These characteristics are the foundation of NCE-MRA technique development.
Chapter 2 introduces commonly used diagnostic angiographic methods, particularly CE-MRA and NCE-MRA. Current NCE-MRA techniques are reviewed and categorised into different types. Their principles, implementations and limitations are summarised.
Chapter 3 describes imaging acceleration theories including compressed sensing (CS), parallel imaging (PI) and partial Fourier (PF). The Split Bregman algorithm is described as an efficient CS reconstruction method. The SPIRiT reconstruction for PI and homodyne detection for PF are also introduced and combined with Split Bregman to form the basis of the reconstruction strategy for undersampled MR datasets. Four image quality metrics are presented for evaluating the quality of reconstructed images.
In Chapter 4, an intensity correction method is proposed to improve background suppression for subtractive NCE-MRA techniques. Residual signals of background tissues are removed by performing a weighted subtraction, in which the weighting factor is obtained by a robust regression method. Image sparsity can also be increased and thereby potentially benefit CS reconstruction in the following chapters.
Chapter 5 investigates the optimal k-space sampling patterns for the 3D accelerated femoral artery FBI sequence. A variable density Poisson-disk with a fully sampled centre region and missing partial Fourier fractions is employed for k-space undersampling in the ky-kz plane. Several key parameters in sampling pattern design, such as partial Fourier sampling ratios, fully sampled centre region size and density decay factor, are evaluated and optimised.
Chapter 6 introduces several reconstruction strategies for accelerated subtractive NCE-MRA. A new reconstruction method, k-space subtraction with phase and intensity correction (KSPIC), is developed. By performing subtraction in k-space, KSPIC can exploit the sparsity of subtracted angiogram data and potentially improve the reconstruction performance. A phase correction procedure is used to restore the polarity of negative signals caused by subtraction. The intensity correction method proposed in Chapter 4 is also incorporated in KSPIC as it improves background suppression and thereby sparsity.
The highly accelerated technique can be used not only to reduce the acquisition time, but also to enable imaging with increased resolution with no time penalty. A time-efficient high-resolution FBI technique is proposed in Chapter 7. By employing KSPIC and modifying the flow-compensation/spoiled gradients, the image matrix size can be increased from 256×256 to up to 512×512 without prolonging the acquisition time.
Chapter 8 summarises the overall achievements and limitations of this thesis, as well as outlines potential future research directions.Cambridge Trust
China Scholarship Council
Addenbrooke’s Charitable Trust
National Institute of Health Research, Cambridge Biomedical Research Cente
- …