SPINAL CORD SEGMENTATION AS OAR IN PLANNING CT FOR RADIOTHERAPY USING HISTOGRAM MATCHING, TEMPLATE MATCHING, AND U-NET

Radiotherapy is one of the major option used in cancer management. The treatment involves several steps, one of which is the construction of a computed tomography (CT) model of the patient so that the target tissues and organs at risk (OARs) surrounding that target can be evaluated. With the CT, the responsible physician delimits the OARs slice by slice, as the spinal cord that comprises almost all the tomography becomes more tiring to be segmented and thus susceptible to errors. Thus, this paper presents a method of spinal cord segmentation in planning CT for radiotherapy using template matching, histogram matching and a fully convolutional neural network. The result achieved an accuracy of 99.38\%, specificity of 99.12\%, sensitivity of 93.83\%, and dice index of 81.33\%, without any segmentation refinemen

Unified wavelet and gaussian filtering for segmentation of CT images; application in segmentation of bone in pelvic CT images

Background The analysis of pelvic CT scans is a crucial step for detecting and assessing the severity of Traumatic Pelvic Injuries. Automating the processing of pelvic CT scans could impact decision accuracy, decrease the time for decision making, and reduce health care cost. This paper discusses a method to automate the segmentation of bone from pelvic CT images. Accurate segmentation of bone is very important for developing an automated assisted-decision support system for Traumatic Pelvic Injury diagnosis and treatment. Methods The automated method for pelvic CT bone segmentation is a hierarchical approach that combines filtering and histogram equalization, for image enhancement, wavelet analysis and automated seeded region growing. Initial results of segmentation are used to identify the region where bone is present and to target histogram equalization towards the specific area. Speckle Reducing Anisotropic Didffusion (SRAD) filter is applied to accentuate the desired features in the region. Automated seeded region growing is performed to refine the initial bone segmentation results. Results The proposed method automatically processes pelvic CT images and produces accurate segmentation. Bone connectivity is achieved and the contours and sizes of bones are true to the actual contour and size displayed in the original image. Results are promising and show great potential for fracture detection and assessing hemorrhage presence and severity. Conclusion Preliminary experimental results of the automated method show accurate bone segmentation. The novelty of the method lies in the unique hierarchical combination of image enhancement and segmentation methods that aims at maximizing the advantages of the combined algorithms. The proposed method has the following advantages: it produces accurate bone segmentation with maintaining bone contour and size true to the original image and is suitable for automated bone segmentation from pelvic CT images

Assessment of the effects of different sample perfusion procedures on phase-contrast tomographic images of mouse spinal cord

Synchrotron X-ray Phase Contrast micro-Tomography (SXrPC\u3bcT) is a powerful tool in the investigation of biological tissues, including the central nervous system (CNS), and it allows to simultaneously detect the vascular and neuronal network avoiding contrast agents or destructive sample preparations. However, specific sample preparation procedures aimed to optimize the achievable contrast- and signal-to-noise ratio (CNR and SNR, respectively) are required. Here we report and discuss the effects of perfusion with two different fixative agents (ethanol and paraformaldehyde) and with a widely used contrast medium (MICROFIL\uae) on mouse spinal cord. As a main result, we found that ethanol enhances contrast at the grey/white matter interface and increases the contrast in correspondence of vascular features and fibres, thus providing an adequate spatial resolution to visualise the vascular network at the microscale. On the other hand, ethanol is known to induce tissue dehydration, likely reducing cell dimensions below the spatial resolution limit imposed by the experimental technique. Nonetheless, neurons remain well visible using either perfused paraformaldehyde or MICROFIL\uae compound, as these latter media do not affect tissues with dehydration effects. Paraformaldehyde appears as the best compromise: it is not a contrast agent, like MICROFIL\uae, but it is less invasive than ethanol and permits to visualise well both cells and blood vessels. However, a quantitative estimation of the relative grey matter volume of each sample has led us to conclude that no significant alterations in the grey matter extension compared to the white matter occur as a consequence of the perfusion procedures tested in this study

Segmentation automatique de la moelle épinière sur des images de résonance magnétique par propagation de modèles déformables

RÉSUMÉ Les lésions de la moelle épinière, induites par des traumas (e.g. accident de la route) ou par des maladies neurodégénératives, touchent plus 85 000 personnes au Canada avec environ 4250 nouveaux cas chaque année1. Elles ont de plus un impact majeur sur la vie quotidienne des personnes atteintes, en provoquant des pertes de sensibilité et de contrôle moteur dont la gravité dépend de la taille et de l’emplacement des lésions. Bien qu’il existe des approches thérapeutiques permettant d’améliorer la réhabilitation fonctionnelle des patients, toutes ces approches se heurtent à une inconnue majeure : l’étendue des dégâts causés par les lésions. Un diagnostic précoce et précis des maladies neurodégénératives touchant la moelle épinière permettrait d’améliorer grandement l’efficacité de leurs traitements. Depuis de nombreuses années, l’IRM a prouvé son potentiel dans le diagnostic et le pronostic des lésions de la moelle épinière (Cadotte, 2011; Cohen-Adad et al., 2011). Ce domaine manque cependant encore d’outils complètement automatisés permettant l’extraction et la comparaison de métriques cliniques reliées à la structure de la moelle (aire de section transverse, volume, etc.). La segmentation de la moelle épinière sur des images IRM anatomiques peut fournir des mesures d’aires et de volumes de la moelle (Losseff et al., 1996) et peut quantifier son atrophie en cas de maladies neurodégénératives telles que la sclérose en plaques (Chen et al., 2013) et la sclérose latérale amyotrophique (Cohen-Adad et al., 2011). Ce projet de maîtrise vise à développer une méthode de segmentation complètement automatique de la moelle épinière, fonctionnant sur plusieurs types d’images IRM (pondérées en T1 et en T2) et sur n’importe quel champ de vue (cervical ou thoracique), et permettant d’extraire et de comparer des mesures précises de la moelle épinière. La revue de la littérature a permis de mettre en évidence le manque de méthode de segmentation automatique de la moelle épinière fonctionnant sur n’importe quel type de contraste et de champ de vue. Elle a toutefois fait ressortir une série de propriétés intéressantes, dans les méthodes semi-automatiques existantes, pouvant être combinées pour former une méthode complètement automatisée.----------ABSTRACT Spinal cord lesions affects more than 85,000 people in Canada with about 4,250 new cases every year. Lesions can be caused by traumatic injuries or by neurodegenerative diseases such as multiple sclerosis. They have an important impact on a patient’s daily life, inducing loss of sensibility or motor control in the human body. The extent of damages caused by a lesion varies with the number of damaged spinal cord tracks, and depends on the size and the position of the lesion within the spinal cord. Although therapeutic approaches for patient functional rehabilitation exist, they all face an unknown variable: the extent of spinal cord lesions. A precise and early diagnosis of neurodegenerative diseases would improve their treatment efficiency. For a number of years, MRI has demonstrated its potential in the diagnosis and prognosis of spinal cord lesions (Cadotte, 2011; Cohen-Adad et al., 2010). However, this research field still lacks of fully automatized tools for the extraction and comparison of clinical metrics related to the spinal cord structure (e.g. cross-sectional area, volumes). Spinal cord segmentation on anatomical MR images can provide accurate area and volume measurements (Losseff et al., 1996) and could quantify spinal cord atrophy caused by neurodegenerative diseases such as multiple sclerosis (Chen et al., 2013) or amyotrophic lateral sclerosis (Cohen-Adad et al., 2011). The objective of this Master’s project is to develop a fully automatic spinal cord segmentation method, working on multiple MR contrasts and any field of view, able to extract and compare accurate spinal cord measurements. The literature review pointed out the lack of such a method but highlighted several interesting features in existing methods, that can be combined to develop a new automatic segmentation algorithm. The method developed in this project is based on the multi-resolution propagation of a deformable model. First, the spinal cord position and orientation is detected in the image using an elliptical Hough transform on multiple adjacent axial slices. A low-resolution tubular mesh is then build around the detection point and direction and deformed on spinal cord edges by minimizing an energy equation. An iterative process, composed by the duplication, translation, orientation and deformation of the mesh, propagates the surface along the spinal cord. Finally, a refinement and a global deformation of the surface provide accurate segmentation of the spinal cord. Measurements can be directly extracted from the segmentation surface. The spinal canal can also be segmented with our method by simply inversing the gradient in the image an

Automatic Segmentation of Anatomical Structures from CT Scans of Thorax for RTP

Modern radiotherapy techniques are vulnerable to delineation inaccuracies owing to the steep dose gradient around the target. In this aspect, accurate contouring comprises an indispensable part of optimal radiation treatment planning (RTP). We suggest a fully automated method to segment the lungs, trachea/main bronchi, and spinal canal accurately from computed tomography (CT) scans of patients with lung cancer to use for RTP. For this purpose, we developed a new algorithm for inclusion of excluded pathological areas into the segmented lungs and a modified version of the fuzzy segmentation by morphological reconstruction for spinal canal segmentation and implemented some image processing algorithms along with them. To assess the accuracy, we performed two comparisons between the automatically obtained results and the results obtained manually by an expert. The average volume overlap ratio values range between 94.30 ± 3.93% and 99.11 ± 0.26% on the two different datasets. We obtained the average symmetric surface distance values between the ranges of 0.28 ± 0.21 mm and 0.89 ± 0.32 mm by using the same datasets. Our method provides favorable results in the segmentation of CT scans of patients with lung cancer and can avoid heavy computational load and might offer expedited segmentation that can be used in RTP

Fusion and Analysis of Multidimensional Medical Image Data

Analýza medicínských obrazů je předmětem základního výzkumu již řadu let. Za tu dobu bylo v této oblasti publikováno mnoho výzkumných prací zabývajících se dílčími částmi jako je rekonstrukce obrazů, restaurace, segmentace, klasifikace, registrace (lícování) a fúze. Kromě obecného úvodu, pojednává tato disertační práce o dvou medicínsky orientovaných tématech, jež byla formulována ve spolupráci s Philips Netherland BV, divizí Philips Healthcare. První téma je zaměřeno na oblast zpracování obrazů subtrakční angiografie dolních končetin člověka získaných pomocí výpočetní X-Ray tomografie (CT). Subtrakční angiografie je obvykle využívaná při podezření na periferní cévní onemocnění (PAOD) nebo při akutním poškození dolních končetin jako jsou fraktury apod. Současné komerční metody nejsou dostatečně spolehlivé už v předzpracování, jako je například odstranění pacientského stolu, pokrývky, dlahy, apod. Spolehlivost a přesnost identifikace cév v subtrahovaných datech vedoucích v blízkosti kostí je v důsledku Partial Volume artefaktu rovněž nízká. Automatické odstranění kalcifikací nebo detekce malých cév doplňujících nezbytnou informaci o náhradním zásobení dolních končetin krví v případě přerušení hlavních zásobujících cév v současné době rovněž nesplňují kritéria pro plně automatické zpracování. Proto hlavním cílem týkající se tohoto tématu bylo vyvinout automatický systém, který by mohl současné nedostatky v CTSA vyšetření odstranit. Druhé téma je orientováno na identifikaci patologických změn na páteři člověka v CT obrazech se zaměřením na osteolytické a osteoblastické léze u jednotlivých obratlů. Tyto změny obvykle nastávají v důsledků postižení metastazujícím procesem rakovinového onemocnění. Pro detekci patologických změn je pak potřeba identifikace a segmentace jednotlivých obratlů. Přesnost analýzy jednotlivých lézí však závisí rovněž na správné identifikaci těla a zadních segmentů u jednotlivých obratlů a na segmentaci trabekulárního centra obratlů, tj. odstranění kortikální kosti. Během léčby mohou být pacienti skenováni vícekrát, obvykle s několika-mesíčním odstupem. Hodnocení případného vývoje již detekovaných patologických změn pak logicky vychází ze správné detekce patologií v jednotlivých obratlech korespondujících si v jednotlivých akvizicích. Jelikož jsou příslušné obratle v jednotlivých akvizicích obvykle na různé pozici, jejich fúze, vedoucí k analýze časového vývoje detekovaných patologií, je komplikovaná. Požadovaným výsledkem v tomto tématu je vytvoření komplexního systému pro detekci patologických změn v páteři, především osteoblastických a osteolytických lézí. Takový systém tedy musí umožnovat jak segmentaci jednotlivých obratlů, jejich automatické rozdělení na hlavní části a odstranění kortikální kosti, tak také detekci patologických změn a jejich hodnocení. Ačkoliv je tato disertační práce v obou výše zmíněných tématech primárně zaměřena na experimentální část zpracování medicínských obrazů, zabývá se všemi nezbytnými kroky, jako je předzpracování, registrace, dodatečné zpracování a hodnocení výsledků, vedoucími k možné aplikovatelnosti obou systému v klinické praxi. Jelikož oba systémy byly řešeny v rámci týmové spolupráce jako celek, u obou témat jsou pro některé konkrétní kroky uvedeny odkazy na doktorskou práci Miloše Malínského.Analysis of medical images has been subject of basic research for many years. Many research papers have been published in the field related to image analysis and focused on partial aspects such as reconstruction, restoration, segmentation and classification, registration (spatial alignment) and fusion. Besides the introduction of related general concepts used in medical image processing, this thesis deals with two specific medical problems formulated in cooperation with Philips Netherland BV, Philips Healthcare division. The first topic is focused on subtraction angiography in patients’ lower legs utilizing image data from X-Ray computed tomography (CT). CT subtraction angiography (CTSA) is typically used for indication of the Peripheral Artery Occlusive Disease (PAOD) and for examination of acute injuries of lower legs such as acute fractures, etc. Current methods in clinical praxis are not sufficient regarding the pre-processing such as masking of patient desk, cover, splint, etc. The subtraction of blood vessels adjacent to neighboring bones in lower legs is of low accuracy due to the Partial Volume artifact. Masking of calcifications and detection of tiny blood vessels complementing necessary information about the alternative blood supply in lower legs in case of obstruction in main arteries is also not reliable for fully automated process presently. Therefore, the main aim regarding this topic was to develop an automated framework that could overcome current shortcomings in CTSA examination. The second topic is oriented on the identification and evaluation of pathologic changes in human spine, focusing on osteolytic and osteoblastic lesions in individual vertebrae in CT images. Such changes occur typically as a consequence of metastasizing process of cancerous disease. For the detection of pathologic changes, an identification and segmentation of individual vertebrae is necessary. Moreover, the analysis of individual lesions in vertebrae depends also on correct identification of vertebral body and posterior segments of each vertebra, and on segmentation of their trabecular centers. Patients are typically examined more than once during their therapy. Then, the evaluation of possible tumorous progression is based on accurate detection of pathologies in individual vertebrae in the base-line and corresponding follow-up images. Since the corresponding vertebrae are in mutually different positions in the follow-up images, their fusion leading to the analysis of the lesion progression is complicated. The main aim regarding this topic is to develop a complex framework for detection of pathologic lesions on spine, with the main focus on osteoblastic and osteolystic lesions. Such system has to provide not only reliable segmentation of individual vertebrae and detection of their main regions but also the masking of their cortical bone, detection of their pathologic changes and their evaluation. Although this dissertation thesis is primarily oriented at the experimental part of medical image processing considering both the above mentioned topics, it deals with all necessary processing steps, i.e. preprocessing, image registration, post-processing and evaluation of results, leading to the future use of both frameworks in clinical practice. Since both frameworks were developed in a team, there are some chapters referring to the dissertation thesis of Milos Malinsky.