A three-dimensional representation method for noisy point clouds based on growing self-organizing maps accelerated on GPUs

The research described in this thesis was motivated by the need of a robust model capable of representing 3D data obtained with 3D sensors, which are inherently noisy. In addition, time constraints have to be considered as these sensors are capable of providing a 3D data stream in real time. This thesis proposed the use of Self-Organizing Maps (SOMs) as a 3D representation model. In particular, we proposed the use of the Growing Neural Gas (GNG) network, which has been successfully used for clustering, pattern recognition and topology representation of multi-dimensional data. Until now, Self-Organizing Maps have been primarily computed offline and their application in 3D data has mainly focused on free noise models, without considering time constraints. It is proposed a hardware implementation leveraging the computing power of modern GPUs, which takes advantage of a new paradigm coined as General-Purpose Computing on Graphics Processing Units (GPGPU). The proposed methods were applied to different problem and applications in the area of computer vision such as the recognition and localization of objects, visual surveillance or 3D reconstruction.

Protein disulfide isomerase as a target for besnoitiosis therapy : molecular characterization and studies of its role in infection and host immune response

Tese de Doutoramento em Ciências Veterinárias na especialidade de Sanidade Animal

Damage Identification of Structures through Vibration-Based Structural Monitoring Systems

The thesis has been carried out within the “SHAPE Project - Predicting Strength Changes in Bridges from Frequency Data Safety, Hazard, and Poly-harmonic Evaluation” (ERA-NET Plus Infravation Call 2014) which dealt with the structural assessment of existing bridges and laboratory structural reproductions through the use of vibration-based monitoring systems, for detecting changes in their natural frequencies and correlating them with the occurrence of damage. The main purpose of this PhD dissertation has been the detection of the variation of the main natural frequencies as a consequence of a previous-established damage configuration provided on a structure. Firstly, the effect of local damage on the modal feature has been discussed mainly concerning a steel frame and a composite steel-concrete bridge. Concerning the variation of the fundamental frequency of the small bridge, the increasing severity of two local damages has been investigated. Moreover, the comparison with a 3D FE model is even presented establishing a link between the dynamic properties and the damage features. Then, moving towards a diffused damage pattern, four concrete beams and a small concrete deck were loaded achieving the yielding of the steel reinforcement. The stiffness deterioration in terms of frequency shifts has been reconsidered by collecting a large set of dynamic experiments on simply supported R.C. beams discussed in the literature. The comparison of the load-frequency curves suggested a significant agreement among all the experiments. Thus, in the framework of damage mechanics, the “breathing cracks” phenomenon has been discussed leading to an analytical formula able to explain the frequency decay observed experimentally. Lastly, some dynamic investigations of two existing bridges and the corresponding FE Models are presented in Chapter 4. Moreover, concerning the bridge in Bologna, two prototypes of a network of accelerometers were installed and the data of a few months of monitoring have been discussed.

Urban air mobility: modelling the future transport system

The advances in the aviation field, particularly the development of electric flying vehicles, as UAV and eVTOL, paved the way for setting Urban Air Mobility (UAM) services. UAM would provide services for passengers, goods and emergencies and could offer faster trips than ground ones. It is expected that early UAM operations will be performed at Very Low-Level airspace as 0-500 m Above Ground Level. The purpose of this research is to both explore the main features of UAM and test an aerial network model, which could be integrated in a multimodal transport system where ground and aerial mobility services are provided. Analyses on UAM transport system involved two sub-systems: the transport demand sub-system, i.e., the mobility requirements, and the transport supply sub-system, i.e., the service and facilities enabling mobility. At first, the UAM demand levels and features for an Airport Shuttle service have been explored through a suitable survey, by combining Revealed and Stated Preference methodologies, and by calibrating some discrete mode choice models. Then, the focus has been on the transport supply model for UAM services, by focusing on both the ground access points (vertiports) and the aerial network model. A suitable three-dimensional urban aerial network (3D-UAN) model that could support fast aerial connections between O/D pairs has been proposed. Some tests have been implemented to verify the feasibility of the proposed model. Some flying vehicles supporting an Airport Shuttle service have been simulated on the aerial network, which has been specified in terms of both topological features and link transport costs. The preliminary results have showed that the proposed 3D-UAN model could be suitable for supporting UAM services. As for transport engineering, the UAM system framework proposed in this thesis paves the way for further research on air-ground multimodality in urban areas.

Analysis of the long-term subsidence in the urban area of Bologna

Land subsidence in urban areas represents a widespread geological hazard and a pressing challenge for modern society. This research focuses on the subsidence process affecting the city of Bologna (Italy). Since the 1960s, Bologna has experienced ground deformation due to aquifers overexploitation that peaked during the 1970s with rates of 10 cm/year. Despite a general reduction in these rates over the subsequent decades, thanks to groundwater regulations policies, recent data underscore a substantial subsidence resurgence. To reconstruct the subsurface stratigraphic architecture of Bologna’s urban area and generate a 3D geological model, a multidisciplinary approach centred on a stratigraphic analysis relying on the lithofacies criterion was adopted. The convergence of the analyses within this framework resulted in partitioning the study area into three geological domains exhibiting unique morphological features and depositional stacking patterns. Subsequently, since long-term data are crucial for a comprehensive understanding of ongoing subsidence, a methodology was developed to generate cumulative ground displacement time series and maps by integrating ground-based and remotely-sensed measurements. While the reconstructed long-term subsidence field consistently aligns with the primary geological variations summarised in the 3D model, the generated cumulative displacement curves systematically match pluriannual trends observed in groundwater level and pumping monitoring data. Lastly, to evaluate the expression of the observed relationships from a geotechnical perspective, a series of one-dimensional subsidence calculations were conducted considering the mechanical properties of the investigated deposits and piezometric data. These analyses provided valuable insight into the overall mechanical behaviour of the existing soils, as well as the post-pumping groundwater level and pore pressure distributions, consistent with field data. The methodological approach employed enables a comprehensive analysis of land subsidence in urban areas, allowing the exploration of individual factors governing the deformation process and their interactions, even within complex stratigraphic and hydrogeological environments such as Bologna’s urban area.

Stress distribution in the cervical region of an upper central incisor in a 3D finite element model

The aim of this study was to evaluate the stress distribution in the cervical region of a sound upper central incisor in two clinical situations, standard and maximum masticatory forces, by means of a 3D model with the highest possible level of fidelity to the anatomic dimensions. Two models with 331,887 linear tetrahedral elements that represent a sound upper central incisor with periodontal ligament, cortical and trabecular bones were loaded at 45º in relation to the tooth's long axis. All structures were considered to be homogeneous and isotropic, with the exception of the enamel (anisotropic). A standard masticatory force (100 N) was simulated on one of the models, while on the other one a maximum masticatory force was simulated (235.9 N). The software used were: PATRAN for pre- and post-processing and Nastran for processing. In the cementoenamel junction area, tensile forces reached 14.7 MPa in the 100 N model, and 40.2 MPa in the 235.9 N model, exceeding the enamel's tensile strength (16.7 MPa). The fact that the stress concentration in the amelodentinal junction exceeded the enamel's tensile strength under simulated conditions of maximum masticatory force suggests the possibility of the occurrence of non-carious cervical lesions such as abfractions.

Development of a new 3D OpenFOAM® solver to model the cooling stage in profile extrusion

[Excerpt] The advantages resulting from the use of numerical modelling tools to support the design of processing equipment are almost consensual. The design of calibration systems in profile extrusion is not an exception . H owever , the complex geome tries and heat exchange phenomena involved in this process require the use of numerical solvers able to model the heat exchange in more than one domain ( calibrator and polymer), the compatibilization of the heat transfer at the profile - calibrator interface and with the ability to deal with complex geometries. The combination of all these features is usually hard to find in commercial software. Moreover , the dimension of the meshes required to ob tain accurate results, result in computational times prohibitive for industrial application. (...)

Context-aware 3D model generation for biomedical applications

Magdeburg, Univ., Fak. für Informatik, Diss., 2014

Ammonium accumulation and cell death in a rat 3D brain cell model of glutaric aciduria type I.

Glutaric aciduria type I (glutaryl-CoA dehydrogenase deficiency) is an inborn error of metabolism that usually manifests in infancy by an acute encephalopathic crisis and often results in permanent motor handicap. Biochemical hallmarks of this disease are elevated levels of glutarate and 3-hydroxyglutarate in blood and urine. The neuropathology of this disease is still poorly understood, as low lysine diet and carnitine supplementation do not always prevent brain damage, even in early-treated patients. We used a 3D in vitro model of rat organotypic brain cell cultures in aggregates to mimic glutaric aciduria type I by repeated administration of 1 mM glutarate or 3-hydroxyglutarate at two time points representing different developmental stages. Both metabolites were deleterious for the developing brain cells, with 3-hydroxyglutarate being the most toxic metabolite in our model. Astrocytes were the cells most strongly affected by metabolite exposure. In culture medium, we observed an up to 11-fold increase of ammonium in the culture medium with a concomitant decrease of glutamine. We further observed an increase in lactate and a concomitant decrease in glucose. Exposure to 3-hydroxyglutarate led to a significantly increased cell death rate. Thus, we propose a three step model for brain damage in glutaric aciduria type I: (i) 3-OHGA causes the death of astrocytes, (ii) deficiency of the astrocytic enzyme glutamine synthetase leads to intracerebral ammonium accumulation, and (iii) high ammonium triggers secondary death of other brain cells. These unexpected findings need to be further investigated and verified in vivo. They suggest that intracerebral ammonium accumulation might be an important target for the development of more effective treatment strategies to prevent brain damage in patients with glutaric aciduria type I.

A MODEL OF CREATINE DEFICIENCY SYNDROMES IN 3D BRAIN CELL CULTURES BY KNOCKDOWN OF GAMT AND SLC6A8 GENES

La créatine joue un rôle essentiel dans le métabolisme cellulaire par sa conversion, par la creatine kinase, en phosphocreatine permettant la régénération de l'ATP. La synthèse de créatine, chez les mammifères, s'effectue par une réaction en deux étapes impliquant Γ arginine: glycine amidinotransférase (AGAT) et la guanidinoacétate méthyltransférase (GAMT). L'entrée de créatine dans les cellules s'effectue par son transporteur, SLC6A8. Les déficiences en créatine, dues au déficit en GAMT, AGAT ou SLC6A8, sont fréquentes et caractérisées par une absence ou une forte baisse de créatine dans le système nerveux central. Alors qu'il est connu que AGAT, GAMT et SLC6A8 sont exprimés par le cerveau, les conséquences des déficiences en créatine sur les cellules nerveuses sont peu comprises. Le but de ce travail était de développer de nouveaux modèles expérimentaux des déficiences en Cr dans des cultures 3D de cellules nerveuses de rat en agrégats au moyen de l'interférence à l'ARN appliquée aux gènes GAMT et SLC6A8. Des séquences interférentes (shRNAs) pour les gènes GAMT et SLC6A8 ont été transduites par des vecteurs viraux AAV (virus adéno-associés), dans les cellules nerveuses en agrégats. Nous avons ainsi démontré une baisse de l'expression de GAMT au niveau protéique (mesuré par western blot), et ARN messager (mesuré par qPCR) ainsi qu'une variation caractérisitique de créatine et guanidinoacétate (mesuré par spectrométrie de masse). Après avoir validé nos modèles, nous avons montré que les knockdown de GAMT ou SLC6A8 affectent le développement des astrocytes et des neurones ou des oligodendrocytes et des astrocytes, respectivement, ainsi qu'une augmentation de la mort cellulaire et des modifications dans le pattern d'activation des voies de signalisation impliquant caspase 3 et p38 MAPK, ayant un rôle dans le processus d'apoptose. - Creatine plays essential roles in energy metabolism by the interconversion, by creatine kinase, to its phosphorylated analogue, phosphocreatine, allowing the regeneration of ATP. Creatine is synthesized in mammals by a two step mechanism involving arginine:glycine amidinotransferase (AGAT) and guanidinoacetate methyltransferase (GAMT). Creatine is taken up by cells by a specific transporter, SLC6A8. Creatine deficiency syndromes, due to defects in GAMT, AGAT and SLC6A8, are among the most frequent inborn errors of metabolism, and are characterized by an absence or a severe decrease of creatine in central nervous system, which is the main tissue affected. While it is known that AGAT, GAMT and SLC6A8 are expressed in CNS, many questions remain on the specific effects of AGAT, GAMT and SLC6A8 deficiencies on brain cells. Our aim was to develop new experimental models of creatine deficiencies by knockdown of GAMT and SLC6A8 genes by RNAi in 3D organotypic rat brain cell cultures in aggregates. Specific shRNAs for the GAMT and SLC6A8 genes were transduced in brain cell aggregates by adeno-associated viruses (AAV). The AAV-transduced shRNAs were able to efficiently knockdown the expression of our genes of interest, as shown by a strong decrease of protein by western blotting, a decrease of mRNA by qPCR or characteristic variations of creatine and guanidinoacetate by tandem mass spectrometry. After having validated our experimental models, we have also shown that GAMT and SLC6A8 knockdown affected the development of astrocytes and neurons or oligodendrocytes and astrocytes, respectively. We also observed an increase of cell death and variations in activation pattern of caspase 3 and p38 MAPK pathways, involved in apoptosis, in our experimental model.

Development of model observers applied to 3D breast tomosynthesis microcalcifications and masses

The development of model observers for mimicking human detection strategies has followed from symmetric signals in simple noise to increasingly complex backgrounds. In this study we implement different model observers for the complex task of detecting a signal in a 3D image stack. The backgrounds come from real breast tomosynthesis acquisitions and the signals were simulated and reconstructed within the volume. Two different tasks relevant to the early detection of breast cancer were considered: detecting an 8 mm mass and detecting a cluster of microcalcifications. The model observers were calculated using a channelized Hotelling observer (CHO) with dense difference-of-Gaussian channels, and a modified (Partial prewhitening [PPW]) observer which was adapted to realistic signals which are not circularly symmetric. The sustained temporal sensitivity function was used to filter the images before applying the spatial templates. For a frame rate of five frames per second, the only CHO that we calculated performed worse than the humans in a 4-AFC experiment. The other observers were variations of PPW and outperformed human observers in every single case. This initial frame rate was a rather low speed and the temporal filtering did not affect the results compared to a data set with no human temporal effects taken into account. We subsequently investigated two higher speeds at 5, 15 and 30 frames per second. We observed that for large masses, the two types of model observers investigated outperformed the human observers and would be suitable with the appropriate addition of internal noise. However, for microcalcifications both only the PPW observer consistently outperformed the humans. The study demonstrated the possibility of using a model observer which takes into account the temporal effects of scrolling through an image stack while being able to effectively detect a range of mass sizes and distributions.

Model-based Segmentation and Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

Computed Tomography (CT) represents the standard imaging modality for tumor volume delineation for radiotherapy treatment planning of retinoblastoma despite some inherent limitations. CT scan is very useful in providing information on physical density for dose calculation and morphological volumetric information but presents a low sensitivity in assessing the tumor viability. On the other hand, 3D ultrasound (US) allows a highly accurate definition of the tumor volume thanks to its high spatial resolution but it is not currently integrated in the treatment planning but used only for diagnosis and follow-up. Our ultimate goal is an automatic segmentation of gross tumor volume (GTV) in the 3D US, the segmentation of the organs at risk (OAR) in the CT and the registration of both modalities. In this paper, we present some preliminary results in this direction. We present 3D active contour-based segmentation of the eye ball and the lens in CT images; the presented approach incorporates the prior knowledge of the anatomy by using a 3D geometrical eye model. The automated segmentation results are validated by comparing with manual segmentations. Then, we present two approaches for the fusion of 3D CT and US images: (i) landmark-based transformation, and (ii) object-based transformation that makes use of eye ball contour information on CT and US images.

Model-based Segmentation and Image Fusion of 3D Computed Tomography and 3D Ultrasound of the Eye for Radiotherapy Planning

For radiotherapy treatment planning of retinoblastoma inchildhood, Computed Tomography (CT) represents thestandard method for tumor volume delineation, despitesome inherent limitations. CT scan is very useful inproviding information on physical density for dosecalculation and morphological volumetric information butpresents a low sensitivity in assessing the tumorviability. On the other hand, 3D ultrasound (US) allows ahigh accurate definition of the tumor volume thanks toits high spatial resolution but it is not currentlyintegrated in the treatment planning but used only fordiagnosis and follow-up. Our ultimate goal is anautomatic segmentation of gross tumor volume (GTV) in the3D US, the segmentation of the organs at risk (OAR) inthe CT and the registration of both. In this paper, wepresent some preliminary results in this direction. Wepresent 3D active contour-based segmentation of the eyeball and the lens in CT images; the presented approachincorporates the prior knowledge of the anatomy by usinga 3D geometrical eye model. The automated segmentationresults are validated by comparing with manualsegmentations. Then, for the fusion of 3D CT and USimages, we present two approaches: (i) landmark-basedtransformation, and (ii) object-based transformation thatmakes use of eye ball contour information on CT and USimages.

Navigator-gated free-breathing 3D balanced FFE projection renal MRA: Comparison with contrast-enhanced breath-hold 3D MRA in a swine model.

A cardiac-triggered, free-breathing, 3D balanced FFE projection renal MR angiography (MRA) technique with a 2D pencil beam aortic labeling pulse for selective aortic spin tagging was developed. For respiratory motion artifact suppression during free breathing, a prospective real-time navigator was implemented for renal MRA. Images obtained with the new approach were compared with standard contrast-enhanced (CE) 3D breath-hold MRA in seven swine. Signal properties and vessel visualization were analyzed. With the presented technique, high-resolution, high-contrast renal projection MRA with superior vessel length visualization (including a greater visible number of distal branches of the renal arteries) compared to standard breath-hold CE-MRA was obtained. The present results warrant clinical studies in patients with renal artery disease.

Automatic Segmentation of the Eye in 3D Magnetic Resonance Imaging: A Novel Statistical Shape Model for Treatment Planning of Retinoblastoma.

PURPOSE: Proper delineation of ocular anatomy in 3-dimensional (3D) imaging is a big challenge, particularly when developing treatment plans for ocular diseases. Magnetic resonance imaging (MRI) is presently used in clinical practice for diagnosis confirmation and treatment planning for treatment of retinoblastoma in infants, where it serves as a source of information, complementary to the fundus or ultrasonographic imaging. Here we present a framework to fully automatically segment the eye anatomy for MRI based on 3D active shape models (ASM), and we validate the results and present a proof of concept to automatically segment pathological eyes. METHODS AND MATERIALS: Manual and automatic segmentation were performed in 24 images of healthy children's eyes (3.29 ± 2.15 years of age). Imaging was performed using a 3-T MRI scanner. The ASM consists of the lens, the vitreous humor, the sclera, and the cornea. The model was fitted by first automatically detecting the position of the eye center, the lens, and the optic nerve, and then aligning the model and fitting it to the patient. We validated our segmentation method by using a leave-one-out cross-validation. The segmentation results were evaluated by measuring the overlap, using the Dice similarity coefficient (DSC) and the mean distance error. RESULTS: We obtained a DSC of 94.90 ± 2.12% for the sclera and the cornea, 94.72 ± 1.89% for the vitreous humor, and 85.16 ± 4.91% for the lens. The mean distance error was 0.26 ± 0.09 mm. The entire process took 14 seconds on average per eye. CONCLUSION: We provide a reliable and accurate tool that enables clinicians to automatically segment the sclera, the cornea, the vitreous humor, and the lens, using MRI. We additionally present a proof of concept for fully automatically segmenting eye pathology. This tool reduces the time needed for eye shape delineation and thus can help clinicians when planning eye treatment and confirming the extent of the tumor.