High-resolution three-dimensional 3 T magnetic resonance angiography for the evaluation of experimental aneurysm in the rabbit

OBJECTIVE: The standard technique of two-dimensional intra-arterial digital subtraction angiography (2D-DSA) for the imaging of experimental rabbit aneurysms is invasive and has considerable surgical risks. Therefore, minimally invasive techniques ideally providing three-dimensional imaging for intervention planning and follow-up are needed. This study evaluates the feasibility and quality of three-dimensional 3-T magnetic resonance angiography (3D-3T-MRA) and compares 3D-3T-MRA with 2D-DSA in experimental aneurysms in the rabbit. METHOD: Three microsurgically created aneurysms in three rabbits were evaluated using 2D-DSA and 3D-3T-MRA. Imaging of the aneurysms was performed 2 weeks after creation using 2D-DSA and contrast-enhanced (CE) MRA. Measurements included aneurysm dome (length and width) and aneurysm neck. Aneurysm volumes were determined using CE-MRA. RESULTS: The measurements of the aneurysms' dimensions and the evaluation of vicinity vessels with both techniques showed a good correlation. The mean aneurysm length, aneurysm width and neck width measured with DSA (6.9, 4.1 and 2.8 mm, respectively) correlated with the measurements performed in 3D-3T-MRA (6.9, 4 and 2.5 mm, respectively). The mean aneurysm volumes measured with CE-MRA was 46.7 mm(3). CONCLUSION: 3D-3T CE-MRA is feasible and less invasive and is a safer imaging alternative to DSA for experimental aneurysm. Additionally, aneurysm technique this precise offers the possibility of repetitive 3D aneurysm volumetry for long-term follow-up studies after endovascular aneurysm occlusion.

Staging of osteoarthritis for clinical trials on femoroacetabular impingement

Future clinical trials investigating the natural history and treatment of femoroacetabular impingement (FAI) will require multimodal staging systems for hip osteoarthritis because the optimal system will differ based on the size of the study population, the specific objective in question, and the time frame in which the investigator expects to see the specified end point. Plain radiographs are readily available, low in cost, and of unquestioned validity, but they are relatively insensitive to early joint damage. MRI allows assessment of both bony and soft-tissue pathology within the joint, and it is much more sensitive for early joint damage because cartilage is visualized directly. Biochemical imaging techniques such as delayed gadolinium-enhanced MRI of cartilage, T2 mapping, and T1rho offer the potential to identify biochemical damage to cartilage before the onset of irreversible tissue loss. In the future, biomarkers may allow earlier detection of osteoarthritis before the development of radiographic evidence of disease.

Inheritance of ectopic ureters in Entlebucher Mountain Dogs

To test the hypothesis of a heritable base of ectopic ureters (EU) in Entlebucher Mountain Dogs (EMD) and to elucidate associated risk factors and mode of inheritance of the disease, 565 EMD were clinically investigated and population genetic analyses performed. Based on the location of the most caudal termination of the ureteral openings, 552 EMD were classified into three phenotype groups trigone, intravesically and extravesically ectopic based on results of abdominal sonography, urethra-cystoscopy and/or contrast-enhanced computed tomography. One-third (32.9%) of the phenotyped animals had normal terminations of both ureters in the bladder trigone, 47.3% had at least one intravesicular ectopic termination and 19.8% had at least one extravesicular ectopic termination. Multivariate mixed logistic regression revealed gender as a risk factor associated with EU as males were more often affected than females. Complex segregation analysis indicated a hereditary basis for EU in EMD and the involvement of a major gene in the occurrence of the extravesicular EU phenotype.

IMPROVING QUANTITATIVE TREATMENT RESPONSE MONITORING WITH DEFORMABLE IMAGE REGISTRATION

Quantitative imaging with 18F-FDG PET/CT has the potential to provide an in vivo assessment of response to radiotherapy (RT). However, comparing tissue tracer uptake in longitudinal studies is often confounded by variations in patient setup and potential treatment induced gross anatomic changes. These variations make true response monitoring for the same anatomic volume a challenge, not only for tumors, but also for normal organs-at-risk (OAR). The central hypothesis of this study is that more accurate image registration will lead to improved quantitation of tissue response to RT with 18F-FDG PET/CT. Employing an in-house developed “demons” based deformable image registration algorithm, pre-RT tumor and parotid gland volumes can be more accurately mapped to serial functional images. To test the hypothesis, specific aim 1 was designed to analyze whether deformably mapping tumor volumes rather than aligning to bony structures leads to superior tumor response assessment. We found that deformable mapping of the most metabolically avid regions improved response prediction (P<0.05). The positive predictive power for residual disease was 63% compared to 50% for contrast enhanced post-RT CT. Specific aim 2 was designed to use parotid gland standardized uptake value (SUV) as an objective imaging biomarker for salivary toxicity. We found that relative change in parotid gland SUV correlated strongly with salivary toxicity as defined by the RTOG/EORTC late effects analytic scale (Spearman’s ρ = -0.96, P<0.01). Finally, the goal of specific aim 3 was to create a phenomenological dose-SUV response model for the human parotid glands. Utilizing only baseline metabolic function and the planned dose distribution, predicting parotid SUV change or salivary toxicity, based upon specific aim 2, became possible. We found that the predicted and observed parotid SUV relative changes were significantly correlated (Spearman’s ρ = 0.94, P<0.01). The application of deformable image registration to quantitative treatment response monitoring with 18F-FDG PET/CT could have a profound impact on patient management. Accurate and early identification of residual disease may allow for more timely intervention, while the ability to quantify and predict toxicity of normal OAR might permit individualized refinement of radiation treatment plan designs.

State-of-the-art aortic imaging: Part II - applications in transcatheter aortic valve replacement and endovascular aortic aneurysm repair.

Transcatheter aortic valve replacement (TAVR) as well as thoracic and abdominal endovascular aortic repair (TEVAR and EVAR) rely on accurate pre- and postprocedural imaging. This review article discusses the application of imaging, including preprocedural assessment and measurements as well as postprocedural imaging of complications. Furthermore, the exciting perspective of computational fluid dynamics (CFD) based on cross-sectional imaging is presented. TAVR is a minimally invasive alternative for treatment of aortic valve stenosis in patients with high age and multiple comorbidities who cannot undergo traditional open surgical repair. Given the lack of direct visualization during the procedure, pre- and peri-procedural imaging forms an essential part of the intervention. Computed tomography angiography (CTA) is the imaging modality of choice for preprocedural planning. Routine postprocedural follow-up is performed by echocardiography to confirm treatment success and detect complications. EVAR and TEVAR are minimally invasive alternatives to open surgical repair of aortic pathologies. CTA constitutes the preferred imaging modality for both preoperative planning and postoperative follow-up including detection of endoleaks. Magnetic resonance imaging is an excellent alternative to CT for postoperative follow-up, and is especially beneficial for younger patients given the lack of radiation. Ultrasound is applied in screening and postoperative follow-up of abdominal aortic aneurysms, but cross-sectional imaging is required once abnormalities are detected. Contrast-enhanced ultrasound may be as sensitive as CTA in detecting endoleaks.

The ovine cerebral venous system: comparative anatomy, visualization, and implications for translational research.

Cerebrovascular diseases are significant causes of death and disability in humans. Improvements in diagnostic and therapeutic approaches strongly rely on adequate gyrencephalic, large animal models being demanded for translational research. Ovine stroke models may represent a promising approach but are currently limited by insufficient knowledge regarding the venous system of the cerebral angioarchitecture. The present study was intended to provide a comprehensive anatomical analysis of the intracranial venous system in sheep as a reliable basis for the interpretation of experimental results in such ovine models. We used corrosion casts as well as contrast-enhanced magnetic resonance venography to scrutinize blood drainage from the brain. This combined approach yielded detailed and, to some extent, novel findings. In particular, we provide evidence for chordae Willisii and lateral venous lacunae, and report on connections between the dorsal and ventral sinuses in this species. For the first time, we also describe venous confluences in the deep cerebral venous system and an 'anterior condylar confluent' as seen in humans. This report provides a detailed reference for the interpretation of venous diagnostic imaging findings in sheep, including an assessment of structure detectability by in vivo (imaging) versus ex vivo (corrosion cast) visualization methods. Moreover, it features a comprehensive interspecies-comparison of the venous cerebral angioarchitecture in man, rodents, canines and sheep as a relevant large animal model species, and describes possible implications for translational cerebrovascular research.

State-of-the-art of imaging detecting endoleaks post-EVAR with special focus on low-flow endoleaks.

Endovascular aortic repair (EVAR) necessitates lifelong surveillance for the patient, in order to detect complications timely. Endoleaks (ELs) are among the most common complications of EVAR. Especially type II ELs can have a very unpredictable clinical course and this can range from spontaneous sealing to aortic rupture. Subgroups of this type of EL need to be identified in order to make a proper risk stratification. Aim of this review is to describe the existing imaging techniques, including their advantages and disadvantages in the context of post-EVAR surveillance with a particular emphasis on low-flow ELs. Low flow ELs cause pressurization of the aortic aneurysm sac with a low velocity filling, leading to difficulty of detection by routine imaging protocols for EVAR surveillance, e.g. bi- or triphasic multislice computed tomographic angiography, magnetic resonance imaging and contrast enhanced ultrasound. In this article, we review the imaging possibilities of ELs and discuss the different imaging strategies available for depicting low flow ELs.

Hepatic steatosis is associated with surgical-site infection after hepatic and colorectal surgery.

BACKGROUND Obesity and increased visceral fat deposits are important risk factors for surgical-site infection (SSI). Interestingly, a potential role of hepatic steatosis on complications after extrahepatic surgery remains unknown. The aim of the present study was to investigate the impact of hepatic steatosis on SSI in patients that underwent open abdominal surgery. METHODS A total of 231 patients that underwent either liver (n = 116) or colorectal (n = 115) resection and received preoperative contrast-enhanced computed tomography scans were retrospectively investigated. Signal attenuation of the liver parenchyma was measured on computed tomography scans to assess hepatic steatosis. RESULTS More SSIs (including types 1, 2, and 3) were found in the group with hepatic steatosis (56/118 [47.5%]) compared with the control group (30/113 [26.6%]; P = .001). Patients with hepatic steatosis showed greater median body mass index than patients without hepatic steatosis (26.6 kg/m(2) [range 16.8-47.0 kg/m(2)] vs 23.2 kg/m(2) [15.9-32.7 kg/m(2)]; P < .001). Patients with hepatic steatosis experienced longer median operation times (297 minutes [52-708 minutes] vs 240 minutes [80-600 minutes]; P = .003). In a multivariate analysis, hepatic steatosis was identified as an independent risk factor for SSI in patients undergoing hepatic (odds ratio 10.33 [95% confidence interval 1.19-89.76]; P = .03) or colorectal (odds ratio 6.67 [95% confidence interval 1.12-39.33]; P = .04) operation. CONCLUSION Hepatic steatosis is associated with SSI after hepatic and colorectal operation.

Feasible Dose Reduction in Routine Chest Computed Tomography Maintaining Constant Image Quality Using the Last Three Scanner Generations: From Filtered Back Projection to Sinogram-affirmed Iterative Reconstruction and Impact of the Novel Fully Integrated Detector Design Minimizing Electronic Noise

OBJECTIVE The aim of the present study was to evaluate a dose reduction in contrast-enhanced chest computed tomography (CT) by comparing the three latest generations of Siemens CT scanners used in clinical practice. We analyzed the amount of radiation used with filtered back projection (FBP) and an iterative reconstruction (IR) algorithm to yield the same image quality. Furthermore, the influence on the radiation dose of the most recent integrated circuit detector (ICD; Stellar detector, Siemens Healthcare, Erlangen, Germany) was investigated. MATERIALS AND METHODS 136 Patients were included. Scan parameters were set to a thorax routine: SOMATOM Sensation 64 (FBP), SOMATOM Definition Flash (IR), and SOMATOM Definition Edge (ICD and IR). Tube current was set constantly to the reference level of 100 mA automated tube current modulation using reference milliamperes. Care kV was used on the Flash and Edge scanner, while tube potential was individually selected between 100 and 140 kVp by the medical technologists at the SOMATOM Sensation. Quality assessment was performed on soft-tissue kernel reconstruction. Dose was represented by the dose length product. RESULTS Dose-length product (DLP) with FBP for the average chest CT was 308 mGy*cm ± 99.6. In contrast, the DLP for the chest CT with IR algorithm was 196.8 mGy*cm ± 68.8 (P = 0.0001). Further decline in dose can be noted with IR and the ICD: DLP: 166.4 mGy*cm ± 54.5 (P = 0.033). The dose reduction compared to FBP was 36.1% with IR and 45.6% with IR/ICD. Signal-to-noise ratio (SNR) was favorable in the aorta, bone, and soft tissue for IR/ICD in combination compared to FBP (the P values ranged from 0.003 to 0.048). Overall contrast-to-noise ratio (CNR) improved with declining DLP. CONCLUSION The most recent technical developments, namely IR in combination with integrated circuit detectors, can significantly lower radiation dose in chest CT examinations.

Dissecting abdominal aortic aneurysm in Ang II-infused mice: suprarenal branch ruptures and apparent luminal dilatation.

AIMS In this work, we provide novel insight into the morphology of dissecting abdominal aortic aneurysms in angiotensin II-infused mice. We demonstrate why they exhibit a large variation in shape and, unlike their human counterparts, are located suprarenally rather than infrarenally. METHODS AND RESULTS We combined synchrotron-based, ultra-high resolution ex vivo imaging (phase contrast X-Ray tomographic microscopy) with in vivo imaging (high-frequency ultrasound and contrast-enhanced micro-CT) and image-guided histology. In all mice, we observed a tear in the tunica media of the abdominal aorta near the ostium of the celiac artery. Independently we found that, unlike the gradual luminal expansion typical for human aneurysms, the outer diameter increase of angiotensin II-induced dissecting aneurysms in mice was related to one or several intramural haematomas. These were caused by ruptures of the tunica media near the ostium of small suprarenal side branches, which had never been detected by the established small animal imaging techniques. The tear near the celiac artery led to apparent luminal dilatation, while the intramural haematoma led to a dissection of the tunica adventitia on the left suprarenal side of the aorta. The number of ruptured branches was higher in those aneurysms that extended into the thoracic aorta, which explained the observed variability in aneurysm shape. CONCLUSION Our results are the first to describe apparent luminal dilatation, suprarenal branch ruptures, and intramural haematoma formation in dissecting abdominal aortic aneurysms in mice. Moreover, we validate and demonstrate the vast potential of phase contrast X-ray tomographic microscopy in cardiovascular small animal applications.

Magnetic resonance imaging features of sinonasal disorders in horses.

Diseases of paranasal sinuses and nasal passages in horses can be a diagnostic challenge because of the complex anatomy of the head and limitations of many diagnostic modalities. Our hypothesis was that magnetic resonance (MR) imaging would provide excellent anatomical detail and soft tissue resolution, and would be accurate in the diagnosis of diseases of the paranasal sinuses and nasal passages in horses. Fourteen horses were imaged. Inclusion criteria were lesions located to the sinuses or nasal passages that underwent MR imaging and subsequent surgical intervention and/or histopathologic examination. A low field, 0.3 tesla open magnet was used. Sequences in the standard protocol were fast spin echo T2 sagittal and transverse, spin echo T1 transverse, short-tau inversion recovery (STIR) dorsal, gradient echo 3D T1 MPR dorsal (plain and contrast enhanced), spin echo T1 fatsat (contrast enhanced). Mean scan time to complete the examination was 53 min (range 39-99 min). Lesions identified were primary or secondary sinusitis (six horses), paranasal sinus cyst (four horses), progressive ethmoid hematoma (two horses), and neoplasia (two horses). The most useful sequences were fast spin echo T2 transverse and sagittal, STIR dorsal and FE3D MPR (survey and contrast enhanced). Fluid accumulation, mucosal thickening, presence of encapsulated contents, bone deformation, and thickening were common findings observed in MR imaging. In selected horses, magnetic resonance imaging is a useful tool in diagnosing lesions of the paranasal sinuses and nasal passages.

Spinal Spot Sign.

INTRODUCTION A marker predictive of hematoma expansion in the central nervous system could aid the selection of patients for hemostatic or surgical treatment. CASE REPORT Here, we present a 83-year-old patient with acute spinal subdural hematoma with paraparesis progressing to paraplegia. A contrast extravasation within the intraspinal hematoma was visualized on spinal MR indicating active bleeding (spinal spot sign). A second acquisition of contrast-enhanced MR images showed progression of contrast extravasation helping to different active bleeding from spinal arteriovenous malformations/fistula. CONCLUSIONS A "spinal spot sign" may be important for treatment decisions, notably in patients with incomplete neurological deficits at the time of imaging.

Real-Time In Vivo Characterization of Primary Liver Tumors With Diffuse Optical Spectroscopy During Percutaneous Needle Interventions: Feasibility Study in Woodchucks

OBJECTIVE This study presents the first in vivo real-time optical tissue characterization during image-guided percutaneous intervention using near-infrared diffuse optical spectroscopy sensing at the tip of a needle. The goal of this study was to indicate transition boundaries from healthy tissue to tumors, namely, hepatic carcinoma, based on the real-time feedback derived from the optical measurements. MATERIALS AND METHODS Five woodchucks with hepatic carcinoma were used for this study. The woodchucks were imaged with contrast-enhanced cone beam computed tomography with a flat panel detector C-arm system to visualize the carcinoma in the liver. In each animal, 3 insertions were performed, starting from the skin surface toward the hepatic carcinoma under image guidance. In 2 woodchucks, each end point of the insertion was confirmed with pathologic examination of a biopsy sample. While advancing the needle in the animals under image guidance such as fluoroscopy overlaid with cone beam computed tomography slice and ultrasound, optical spectra were acquired at the distal end of the needles. Optical tissue characterization was determined by translating the acquired optical spectra into clinical parameters such as blood, water, lipid, and bile fractions; tissue oxygenation levels; and scattering amplitude related to tissue density. The Kruskal-Wallis test was used to study the difference in the derived clinical parameters from the measurements performed within the healthy tissue and the hepatic carcinoma. Kurtoses were calculated to assess the dispersion of these parameters within the healthy and carcinoma tissues. RESULTS Blood and lipid volume fractions as well as tissue oxygenation and reduced scattering amplitude showed to be significantly different between the healthy part of the liver and the hepatic carcinoma (P < 0.05) being higher in normal liver tissue. A decrease in blood and lipid volume fractions and tissue oxygenation as well as an increase in scattering amplitude were observed when the tip of the needle crossed the margin from the healthy liver tissue to the carcinoma. The kurtosis for each derived clinical parameter was high in the hepatic tumor as compared with that in the healthy liver indicating intracarcinoma variability. CONCLUSIONS Tissue blood content, oxygenation level, lipid content, and tissue density all showed significant differences when the needle tip was guided from the healthy tissue to the carcinoma and can therefore be used to identify tissue boundaries during percutaneous image-guided interventions.

Multimodal non-rigid image registration methods for therapy and surgery planning

La planificación pre-operatoria se ha convertido en una tarea esencial en cirugías y terapias de marcada complejidad, especialmente aquellas relacionadas con órgano blando. Un ejemplo donde la planificación preoperatoria tiene gran interés es la cirugía hepática. Dicha planificación comprende la detección e identificación precisa de las lesiones individuales y vasos así como la correcta segmentación y estimación volumétrica del hígado funcional. Este proceso es muy importante porque determina tanto si el paciente es un candidato adecuado para terapia quirúrgica como la definición del abordaje a seguir en el procedimiento. La radioterapia de órgano blando es un segundo ejemplo donde la planificación se requiere tanto para la radioterapia externa convencional como para la radioterapia intraoperatoria. La planificación comprende la segmentación de tumor y órganos vulnerables y la estimación de la dosimetría. La segmentación de hígado funcional y la estimación volumétrica para planificación de la cirugía se estiman habitualmente a partir de imágenes de tomografía computarizada (TC). De igual modo, en la planificación de radioterapia, los objetivos de la radiación se delinean normalmente sobre TC. Sin embargo, los avances en las tecnologías de imagen de resonancia magnética (RM) están ofreciendo progresivamente ventajas adicionales. Por ejemplo, se ha visto que el ratio de detección de metástasis hepáticas es significativamente superior en RM con contraste Gd–EOB–DTPA que en TC. Por tanto, recientes estudios han destacado la importancia de combinar la información de TC y RM para conseguir el mayor nivel posible de precisión en radioterapia y para facilitar una descripción precisa de las lesiones del hígado. Con el objetivo de mejorar la planificación preoperatoria en ambos escenarios se precisa claramente de un algoritmo de registro no rígido de imagen. Sin embargo, la gran mayoría de sistemas comerciales solo proporcionan métodos de registro rígido. Las medidas de intensidad de voxel han demostrado ser criterios de similitud de imágenes robustos, y, entre ellas, la Información Mutua (IM) es siempre la primera elegida en registros multimodales. Sin embargo, uno de los principales problemas de la IM es la ausencia de información espacial y la asunción de que las relaciones estadísticas entre las imágenes son homogéneas a lo largo de su domino completo. La hipótesis de esta tesis es que la incorporación de información espacial de órganos al proceso de registro puede mejorar la robustez y calidad del mismo, beneficiándose de la disponibilidad de las segmentaciones clínicas. En este trabajo, se propone y valida un esquema de registro multimodal no rígido 3D usando una nueva métrica llamada Información Mutua Centrada en el Órgano (Organ-Focused Mutual Information metric (OF-MI)) y se compara con la formulación clásica de la Información Mutua. Esto permite mejorar los resultados del registro en áreas problemáticas incorporando información regional al criterio de similitud, beneficiándose de la disponibilidad real de segmentaciones en protocolos estándares clínicos, y permitiendo que la dependencia estadística entre las dos modalidades de imagen difiera entre órganos o regiones. El método propuesto se ha aplicado al registro de TC y RM con contraste Gd–EOB–DTPA así como al registro de imágenes de TC y MR para planificación de radioterapia intraoperatoria rectal. Adicionalmente, se ha desarrollado un algoritmo de apoyo de segmentación 3D basado en Level-Sets para la incorporación de la información de órgano en el registro. El algoritmo de segmentación se ha diseñado específicamente para la estimación volumétrica de hígado sano funcional y ha demostrado un buen funcionamiento en un conjunto de imágenes de TC abdominales. Los resultados muestran una mejora estadísticamente significativa de OF-MI comparada con la Información Mutua clásica en las medidas de calidad de los registros; tanto con datos simulados (p<0.001) como con datos reales en registro hepático de TC y RM con contraste Gd– EOB–DTPA y en registro para planificación de radioterapia rectal usando OF-MI multi-órgano (p<0.05). Adicionalmente, OF-MI presenta resultados más estables con menor dispersión que la Información Mutua y un comportamiento más robusto con respecto a cambios en la relación señal-ruido y a la variación de parámetros. La métrica OF-MI propuesta en esta tesis presenta siempre igual o mayor precisión que la clásica Información Mutua y consecuentemente puede ser una muy buena alternativa en aplicaciones donde la robustez del método y la facilidad en la elección de parámetros sean particularmente importantes. Abstract Pre-operative planning has become an essential task in complex surgeries and therapies, especially for those affecting soft tissue. One example where soft tissue preoperative planning is of high interest is liver surgery. It involves the accurate detection and identification of individual liver lesions and vessels as well as the proper functional liver segmentation and volume estimation. This process is very important because it determines whether the patient is a suitable candidate for surgical therapy and the type of procedure. Soft tissue radiation therapy is a second example where planning is required for both conventional external and intraoperative radiotherapy. It involves the segmentation of the tumor target and vulnerable organs and the estimation of the planned dose. Functional liver segmentations and volume estimations for surgery planning are commonly estimated from computed tomography (CT) images. Similarly, in radiation therapy planning, targets to be irradiated and healthy and vulnerable tissues to be protected from irradiation are commonly delineated on CT scans. However, developments in magnetic resonance imaging (MRI) technology are progressively offering advantages. For instance, the hepatic metastasis detection rate has been found to be significantly higher in Gd–EOB–DTPAenhanced MRI than in CT. Therefore, recent studies highlight the importance of combining the information from CT and MRI to achieve the highest level of accuracy in radiotherapy and to facilitate accurate liver lesion description. In order to improve those two soft tissue pre operative planning scenarios, an accurate nonrigid image registration algorithm is clearly required. However, the vast majority of commercial systems only provide rigid registration. Voxel intensity measures have been shown to be robust measures of image similarity, and among them, Mutual Information (MI) is always the first candidate in multimodal registrations. However, one of the main drawbacks of Mutual Information is the absence of spatial information and the assumption that statistical relationships between images are the same over the whole domain of the image. The hypothesis of the present thesis is that incorporating spatial organ information into the registration process may improve the registration robustness and quality, taking advantage of the clinical segmentations availability. In this work, a multimodal nonrigid 3D registration framework using a new Organ- Focused Mutual Information metric (OF-MI) is proposed, validated and compared to the classical formulation of the Mutual Information (MI). It allows improving registration results in problematic areas by adding regional information into the similitude criterion taking advantage of actual segmentations availability in standard clinical protocols and allowing the statistical dependence between the two modalities differ among organs or regions. The proposed method is applied to CT and T1 weighted delayed Gd–EOB–DTPA-enhanced MRI registration as well as to register CT and MRI images in rectal intraoperative radiotherapy planning. Additionally, a 3D support segmentation algorithm based on Level-Sets has been developed for the incorporation of the organ information into the registration. The segmentation algorithm has been specifically designed for the healthy and functional liver volume estimation demonstrating good performance in a set of abdominal CT studies. Results show a statistical significant improvement of registration quality measures with OF-MI compared to MI with both simulated data (p<0.001) and real data in liver applications registering CT and Gd–EOB–DTPA-enhanced MRI and in registration for rectal radiotherapy planning using multi-organ OF-MI (p<0.05). Additionally, OF-MI presents more stable results with smaller dispersion than MI and a more robust behavior with respect to SNR changes and parameters variation. The proposed OF-MI always presents equal or better accuracy than the classical MI and consequently can be a very convenient alternative within applications where the robustness of the method and the facility to choose the parameters are particularly important.

Generación de modelos físicos personalizados y segmentación de la pared arterial y el trombo intraluminal en imágenes de aneurismas aórticos abdominales

Esta tesis doctoral está encuadrada dentro del marco general de la ingeniería biomédica aplicada al tratamiento de las enfermedades cardiovasculares, enfermedades que provocan alrededor de 1.9 millones (40%) de muertes al año en la Unión Europea. En este contexto surge el proyecto europeo SCATh-Smart Catheterization, cuyo objetivo principal es mejorar los procedimientos de cateterismo aórtico introduciendo nuevas tecnologías de planificación y navegación quirúrgica y minimizando el uso de fluoroscopía. En particular, esta tesis aborda el modelado y diagnóstico de aneurismas aórticos abdominales (AAA) y del trombo intraluminal (TIL), allí donde esté presente, así como la segmentación de estas estructuras en imágenes preoperatorias de RM. Los modelos físicos específicos del paciente, construidos a partir de imágenes médicas preoperatorias, tienen múltiples usos, que van desde la evaluación preoperatoria de estructuras anatómicas a la planificación quirúrgica para el guiado de catéteres. En el diagnóstico y tratamiento de AAA, los modelos físicos son útiles a la hora de evaluar diversas variables biomecánicas y fisiológicas de las estructuras vasculares. Existen múltiples técnicas que requieren de la generación de modelos físicos que representen la anatomía vascular. Una de las principales aplicaciones de los modelos físicos es el análisis de elementos finitos (FE). Las simulaciones de FE para AAA pueden ser específicas para el paciente y permiten modelar estados de estrés complejos, incluyendo los efectos provocados por el TIL. La aplicación de métodos numéricos de análisis tiene como requisito previo la generación de una malla computacional que representa la geometría de interés mediante un conjunto de elementos poliédricos, siendo los hexaédricos los que presentan mejores resultados. En las estructuras vasculares, generar mallas hexaédricas es un proceso especialmente exigente debido a la compleja anatomía 3D ramificada. La mayoría de los AAA se encuentran situados en la bifurcación de la arteria aorta en las arterias iliacas y es necesario modelar de manera fiel dicha bifurcación. En el caso de que la sangre se estanque en el aneurisma provocando un TIL, éste forma una estructura adyacente a la pared aórtica. De este modo, el contorno externo del TIL es el mismo que el contorno interno de la pared, por lo que las mallas resultantes deben reflejar esta particularidad, lo que se denomina como "mallas conformadas". El fin último de este trabajo es modelar las estructuras vasculares de modo que proporcionen nuevas herramientas para un mejor diagnóstico clínico, facilitando medidas de riesgo de rotura de la arteria, presión sistólica o diastólica, etc. Por tanto, el primer objetivo de esta tesis es diseñar un método novedoso y robusto para generar mallas hexaédricas tanto de la pared aórtica como del trombo. Para la identificación de estas estructuras se utilizan imágenes de resonancia magnética (RM). Deben mantenerse sus propiedades de adyacencia utilizando elementos de alta calidad, prestando especial atención al modelado de la bifurcación y a que sean adecuadas para el análisis de FE. El método tiene en cuenta la evolución de la línea central del vaso en el espacio tridimensional y genera la malla directamente a partir de las imágenes segmentadas, sin necesidad de reconstruir superficies triangulares. Con el fin de reducir la intervención del usuario en el proceso de generación de las mallas, es también objetivo de esta tesis desarrollar un método de segmentación semiautomática de las distintas estructuras de interés. Las principales contribuciones de esta tesis doctoral son: 1. El diseño, implementación y evaluación de un algoritmo de generación de mallas hexaédricas conformadas de la pared y el TIL a partir de los contornos segmentados en imágenes de RM. Se ha llevado a cabo una evaluación de calidad que determine su aplicabilidad a métodos de FE. Los resultados demuestran que el algoritmo desarrollado genera mallas conformadas de alta calidad incluso en la región de la bifurcación, que son adecuadas para su uso en métodos de análisis de FE. 2. El diseño, implementación y evaluación de un método de segmentación automático de las estructuras de interés. La luz arterial se segmenta de manera semiautomática utilizando un software disponible a partir de imágenes de RM con contraste. Los resultados de este proceso sirven de inicialización para la segmentación automática de las caras interna y externa de la pared aórtica utilizando métodos basado en modelos de textura y forma a partir de imágenes de RM sin contraste. Los resultados demuestran que el algoritmo desarrollado proporciona segmentaciones fieles de las distintas estructuras de interés. En conclusión, el trabajo realizado en esta tesis doctoral corrobora las hipótesis de investigación postuladas, y pretende servir como aportación para futuros avances en la generación de modelos físicos de geometrías biológicas. ABSTRACT The frame of this PhD Thesis is the biomedical engineering applied to the treatment of cardiovascular diseases, which cause around 1.9 million deaths per year in the European Union and suppose about 40% of deaths per year. In this context appears the European project SCATh-Smart Catheterization. The main objective of this project is creating a platform which improves the navigation of catheters in aortic catheterization minimizing the use of fluoroscopy. In the framework of this project, the specific field of this PhD Thesis is the diagnosis and modeling of abdominal aortic aneurysm (AAAs) and the intraluminal thrombus (ILT) whenever it is present. Patient-specific physical models built from preoperative imaging are becoming increasingly important in the area of minimally invasive surgery. These models can be employed for different purposes, such as the preoperatory evaluation of anatomic structures or the surgical planning for catheter guidance. In the specific case of AAA diagnosis and treatment, physical models are especially useful for evaluating pressures over vascular structures. There are multiple techniques that require the generation of physical models which represent the target anatomy. Finite element (FE) analysis is one the principal applications for physical models. FE simulations for AAA may be patient-specific and allow modeling biomechanical and physiological variables including those produced by ILT, and also the segmentation of those anatomical structures in preoperative MR images. Applying numeric methods requires the generation of a proper computational mesh. These meshes represent the patient anatomy using a set of polyhedral elements, with hexahedral elements providing better results. In the specific case of vascular structures, generating hexahedral meshes is a challenging task due to the complex 3D branching anatomy. Each patient’s aneurysm is unique, characterized by its location and shape, and must be accurately represented for subsequent analyses to be meaningful. Most AAAs are located in the region where the aorta bifurcates into the iliac arteries and it is necessary to model this bifurcation precisely and reliably. If blood stagnates in the aneurysm and forms an ILT, it exists as a conforming structure with the aortic wall, i.e. the ILT’s outer contour is the same as the wall’s inner contour. Therefore, resulting meshes must also be conforming. The main objective of this PhD Thesis is designing a novel and robust method for generating conforming hexahedral meshes for the aortic wall and the thrombus. These meshes are built using largely high-quality elements, especially at the bifurcation, that are suitable for FE analysis of tissue stresses. The method accounts for the evolution of the vessel’s centerline which may develop outside a single plane, and generates the mesh directly from segmented images without the requirement to reconstruct triangular surfaces. In order to reduce the user intervention in the mesh generation process is also a goal of this PhD. Thesis to develop a semiautomatic segmentation method for the structures of interest. The segmentation is performed from magnetic resonance image (MRI) sequences that have tuned to provide high contrast for the arterial tissue against the surrounding soft tissue, so that we determine the required information reliably. The main contributions of this PhD Thesis are: 1. The design, implementation and evaluation of an algorithm for generating hexahedral conforming meshes of the arterial wall and the ILT from the segmented contours. A quality inspection has been applied to the meshes in order to determine their suitability for FE methods. Results show that the developed algorithm generates high quality conforming hexahedral meshes even at the bifurcation region. Thus, these meshes are suitable for FE analysis. 2. The design, implementation and evaluation of a semiautomatic segmentation method for the structures of interest. The lumen is segmented in a semiautomatic way from contrast filled MRI using an available software. The results obtained from this process are used to initialize the automatic segmentation of the internal and external faces of the aortic wall. These segmentations are performed by methods based on texture and shape models from MRI with no contrast. The results show that the algorithm provides faithful segmentations of the structures of interest requiring minimal user intervention. In conclusion, the work undertaken in this PhD. Thesis verifies the investigation hypotheses. It intends to serve as basis for future physical model generation of proper biological anatomies used by numerical methods.