Cancer cell viscoelasticity measurement by quantitative phase and flow stress induction.

Cell viscoelastic properties are affected by the cell cycle, differentiation, and pathological processes such as malignant transformation. Therefore, evaluation of the mechanical properties of the cells proved to be an approach to obtaining information on the functional state of the cells. Most of the currently used methods for cell mechanophenotyping are limited by low robustness or the need for highly expert operation. In this paper, the system and method for viscoelasticity measurement using shear stress induction by fluid flow is described and tested. Quantitative phase imaging (QPI) is used for image acquisition because this technique enables one to quantify optical path length delays introduced by the sample, thus providing a label-free objective measure of morphology and dynamics. Viscosity and elasticity determination were refined using a new approach based on the linear system model and parametric deconvolution. The proposed method allows high-throughput measurements during live-cell experiments and even through a time lapse, whereby we demonstrated the possibility of simultaneous extraction of shear modulus, viscosity, cell morphology, and QPI-derived cell parameters such as circularity or cell mass. Additionally, the proposed method provides a simple approach to measure cell refractive index with the same setup, which is required for reliable cell height measurement with QPI, an essential parameter for viscoelasticity calculation. Reliability of the proposed viscoelasticity measurement system was tested in several experiments including cell types of different Young/shear modulus and treatment with cytochalasin D or docetaxel, and an agreement with atomic force microscopy was observed. The applicability of the proposed approach was also confirmed by a time-lapse experiment with cytochalasin D washout, whereby an increase of stiffness corresponded to actin repolymerization in time.

The Effect of Rhodamine-Derived Superparamagnetic Maghemite Nanoparticles on the Motility of Human Mesenchymal Stem Cells and Mouse Embryonic Fibroblast Cells.

Nanoparticles have become popular in life sciences in the last few years. They have been produced in many variants and have recently been used in both biological experiments and in clinical applications. Due to concerns over nanomaterial risks, there has been a dramatic increase in investigations focused on safety research. The aim of this paper is to present the advanced testing of rhodamine-derived superparamagnetic maghemite nanoparticles (SAMN-R), which are used for their nontoxicity, biocompatibility, biodegradability, and magnetic properties. Recent results were expanded upon from the basic cytotoxic tests to evaluate cell proliferation and migration potential. Two cell types were used for the cell proliferation and tracking study: mouse embryonic fibroblast cells (3T3) and human mesenchymal stem cells (hMSCs). Advanced microscopic methods allowed for the precise quantification of the function of both cell types. This study has demonstrated that a dose of nanoparticles lower than 20 µg·cm-2 per area of the dish does not negatively affect the cells' morphology, migration, cytoskeletal function, proliferation, potential for wound healing, and single-cell migration in comparison to standard CellTracker™ Green CMFDA (5-chloromethylfluorescein diacetate). A higher dose of nanoparticles could be a potential risk for cytoskeletal folding and detachment of the cells from the solid extracellular matrix.

Accurate micro-computed tomography imaging of pore spaces in collagen-based scaffold.

In this work we have used X-ray micro-computed tomography (µCT) as a method to observe the morphology of 3D porous pure collagen and collagen-composite scaffolds useful in tissue engineering. Two aspects of visualizations were taken into consideration: improvement of the scan and investigation of its sensitivity to the scan parameters. Due to the low material density some parts of collagen scaffolds are invisible in a µCT scan. Therefore, here we present different contrast agents, which increase the contrast of the scanned biopolymeric sample for µCT visualization. The increase of contrast of collagenous scaffolds was performed with ceramic hydroxyapatite microparticles (HAp), silver ions (Ag(+)) and silver nanoparticles (Ag-NPs). Since a relatively small change in imaging parameters (e.g. in 3D volume rendering, threshold value and µCT acquisition conditions) leads to a completely different visualized pattern, we have optimized these parameters to obtain the most realistic picture for visual and qualitative evaluation of the biopolymeric scaffold. Moreover, scaffold images were stereoscopically visualized in order to better see the 3D biopolymer composite scaffold morphology. However, the optimized visualization has some discontinuities in zoomed view, which can be problematic for further analysis of interconnected pores by commonly used numerical methods. Therefore, we applied the locally adaptive method to solve discontinuities issue. The combination of contrast agent and imaging techniques presented in this paper help us to better understand the structure and morphology of the biopolymeric scaffold that is crucial in the design of new biomaterials useful in tissue engineering.

Assessment of retinal vein pulsation through video-ophthalmoscopy and simultaneous biosignals acquisition.

The phenomenon of retinal vein pulsation is still not a deeply understood topic in retinal hemodynamics. In this paper, we present a novel hardware solution for recording retinal video sequences and physiological signals using synchronized acquisition, we apply the photoplethysmographic principle for the semi-automatic processing of retinal video sequences and we analyse the timing of the vein collapse within the cardiac cycle using of an electrocardiographic signal (ECG). We measured the left eyes of healthy subjects and determined the phases of vein collapse within the cardiac cycle using a principle of photoplethysmography and a semi-automatic image processing approach. We found that the time to vein collapse (Tvc) is between 60 ms and 220 ms after the R-wave of the ECG signal, which corresponds to 6% to 28% of the cardiac cycle. We found no correlation between Tvc and the duration of the cardiac cycle and only a weak correlation between Tvc and age (0.37, p = 0.20), and Tvc and systolic blood pressure (-0.33, p = 0.25). The Tvc values are comparable to those of previously published papers and can contribute to the studies that analyze vein pulsations.

In Vivo Repeatability and Multiscanner Reproducibility of MRI Radiomics Features in Patients With Monoclonal Plasma Cell Disorders: A Prospective Bi-institutional Study.

OBJECTIVES: Despite the extensive number of publications in the field of radiomics, radiomics algorithms barely enter large-scale clinical application. Supposedly, the low external generalizability of radiomics models is one of the main reasons, which hinders the translation from research to clinical application. The objectives of this study were to investigate reproducibility of radiomics features (RFs) in vivo under variation of patient positioning, magnetic resonance imaging (MRI) sequence, and MRI scanners, and to identify a subgroup of RFs that shows acceptable reproducibility across all different acquisition scenarios. MATERIALS AND METHODS: Between November 30, 2020 and February 16, 2021, 55 patients with monoclonal plasma cell disorders were included in this prospective, bi-institutional, single-vendor study. Participants underwent one reference scan at a 1.5 T MRI scanner and several retest scans: once after simple repositioning, once with a second MRI protocol, once at another 1.5 T scanner, and once at a 3 T scanner. Radiomics feature from the bone marrow of the left hip bone were extracted, both from original scans and after different image normalizations. Intraclass correlation coefficient (ICC) was used to assess RF repeatability and reproducibility. RESULTS: Fifty-five participants (mean age, 59 ± 7 years; 36 men) were enrolled. For T1-weighted images after muscle normalization, in the simple test-retest experiment, 110 (37%) of 295 RFs showed an ICC ≥0.8: 54 (61%) of 89 first-order features (FOFs), 35 (95%) of 37 volume and shape features, and 21 (12%) of 169 texture features (TFs). When the retest was performed with different technical settings, even after muscle normalization, the number of FOF/TF with an ICC ≥0.8 declined to 58/13 for the second protocol, 29/7 for the second 1.5 T scanner, and 49/7 for the 3 T scanner, respectively. Twenty-five (28%) of the 89 FOFs and 6 (4%) of the 169 TFs from muscle-normalized T1-weighted images showed an ICC ≥0.8 throughout all repeatability and reproducibility experiments. CONCLUSIONS: In vivo, only few RFs are reproducible with different MRI sequences or different MRI scanners, even after application of a simple image normalization. Radiomics features selected by a repeatability experiment only are not necessarily suited to build radiomics models for multicenter clinical application. This study isolated a subset of RFs, which are robust to variations in MRI acquisition observed in scanners from 1 vendor, and therefore are candidates to build reproducible radiomics models for monoclonal plasma cell disorders for multicentric applications, at least when centers are equipped with scanners from this vendor.

Deep Learning for Automatic Bone Marrow Apparent Diffusion Coefficient Measurements From Whole-Body Magnetic Resonance Imaging in Patients With Multiple Myeloma: A Retrospective Multicenter Study.

OBJECTIVES: Diffusion-weighted magnetic resonance imaging (MRI) is increasingly important in patients with multiple myeloma (MM). The objective of this study was to train and test an algorithm for automatic pelvic bone marrow analysis from whole-body apparent diffusion coefficient (ADC) maps in patients with MM, which automatically segments pelvic bones and subsequently extracts objective, representative ADC measurements from each bone. MATERIALS AND METHODS: In this retrospective multicentric study, 180 MRIs from 54 patients were annotated (semi)manually and used to train an nnU-Net for automatic, individual segmentation of the right hip bone, the left hip bone, and the sacral bone. The quality of the automatic segmentation was evaluated on 15 manually segmented whole-body MRIs from 3 centers using the dice score. In 3 independent test sets from 3 centers, which comprised a total of 312 whole-body MRIs, agreement between automatically extracted mean ADC values from the nnU-Net segmentation and manual ADC measurements from 2 independent radiologists was evaluated. Bland-Altman plots were constructed, and absolute bias, relative bias to mean, limits of agreement, and coefficients of variation were calculated. In 56 patients with newly diagnosed MM who had undergone bone marrow biopsy, ADC measurements were correlated with biopsy results using Spearman correlation. RESULTS: The ADC-nnU-Net achieved automatic segmentations with mean dice scores of 0.92, 0.93, and 0.85 for the right pelvis, the left pelvis, and the sacral bone, whereas the interrater experiment gave mean dice scores of 0.86, 0.86, and 0.77, respectively. The agreement between radiologists' manual ADC measurements and automatic ADC measurements was as follows: the bias between the first reader and the automatic approach was 49 × 10 -6 mm 2 /s, 7 × 10 -6 mm 2 /s, and -58 × 10 -6 mm 2 /s, and the bias between the second reader and the automatic approach was 12 × 10 -6 mm 2 /s, 2 × 10 -6 mm 2 /s, and -66 × 10 -6 mm 2 /s for the right pelvis, the left pelvis, and the sacral bone, respectively. The bias between reader 1 and reader 2 was 40 × 10 -6 mm 2 /s, 8 × 10 -6 mm 2 /s, and 7 × 10 -6 mm 2 /s, and the mean absolute difference between manual readers was 84 × 10 -6 mm 2 /s, 65 × 10 -6 mm 2 /s, and 75 × 10 -6 mm 2 /s. Automatically extracted ADC values significantly correlated with bone marrow plasma cell infiltration ( R = 0.36, P = 0.007). CONCLUSIONS: In this study, a nnU-Net was trained that can automatically segment pelvic bone marrow from whole-body ADC maps in multicentric data sets with a quality comparable to manual segmentations. This approach allows automatic, objective bone marrow ADC measurements, which agree well with manual ADC measurements and can help to overcome interrater variability or nonrepresentative measurements. Automatically extracted ADC values significantly correlate with bone marrow plasma cell infiltration and might be of value for automatic staging, risk stratification, or therapy response assessment.

Prediction of Bone Marrow Biopsy Results From MRI in Multiple Myeloma Patients Using Deep Learning and Radiomics.

OBJECTIVES: In multiple myeloma and its precursor stages, plasma cell infiltration (PCI) and cytogenetic aberrations are important for staging, risk stratification, and response assessment. However, invasive bone marrow (BM) biopsies cannot be performed frequently and multifocally to assess the spatially heterogenous tumor tissue. Therefore, the goal of this study was to establish an automated framework to predict local BM biopsy results from magnetic resonance imaging (MRI). MATERIALS AND METHODS: This retrospective multicentric study used data from center 1 for algorithm training and internal testing, and data from center 2 to 8 for external testing. An nnU-Net was trained for automated segmentation of pelvic BM from T1-weighted whole-body MRI. Radiomics features were extracted from these segmentations, and random forest models were trained to predict PCI and the presence or absence of cytogenetic aberrations. Pearson correlation coefficient and the area under the receiver operating characteristic were used to evaluate the prediction performance for PCI and cytogenetic aberrations, respectively. RESULTS: A total of 672 MRIs from 512 patients (median age, 61 years; interquartile range, 53-67 years; 307 men) from 8 centers and 370 corresponding BM biopsies were included. The predicted PCI from the best model was significantly correlated ( P ≤ 0.01) to the actual PCI from biopsy in all internal and external test sets (internal test set: r = 0.71 [0.51, 0.83]; center 2, high-quality test set: r = 0.45 [0.12, 0.69]; center 2, other test set: r = 0.30 [0.07, 0.49]; multicenter test set: r = 0.57 [0.30, 0.76]). The areas under the receiver operating characteristic of the prediction models for the different cytogenetic aberrations ranged from 0.57 to 0.76 for the internal test set, but no model generalized well to all 3 external test sets. CONCLUSIONS: The automated image analysis framework established in this study allows for noninvasive prediction of a surrogate parameter for PCI, which is significantly correlated to the actual PCI from BM biopsy.

Collaborative and individual learning of geography in immersive virtual reality: An effectiveness study.

Many university-taught courses moved to online form since the outbreak of the global pandemic of coronavirus disease (COVID-19). Distance learning has become broadly used as a result of the widely applied lockdowns, however, many students lack personal contact in the learning process. A classical web-based distance learning does not provide means for natural interpersonal interaction. The technology of immersive virtual reality (iVR) may mitigate this problem. Current research has been aimed mainly at specific instances of collaborative immersive virtual environment (CIVE) applications for learning. The fields utilizing iVR for knowledge construction and skills training with the use of spatial visualizations show promising results. The objective of this study was to assess the effectiveness of collaborative and individual use of iVR for learning geography, specifically training in hypsography. Furthermore, the study's goals were to determine whether collaborative learning would be more effective and to investigate the key elements in which collaborative and individual learning were expected to differ-motivation and use of cognitive resources. The CIVE application developed at Masaryk University was utilized to train 80 participants in inferring conclusions from cartographic visualizations. The collaborative and individual experimental group underwent a research procedure consisting of a pretest, training in iVR, posttest, and questionnaires. A statistical comparison between the geography pretest and posttest for the individual learning showed a significant increase in the score (p = 0.024, ES = 0.128) and speed (p = 0.027, ES = 0.123), while for the collaborative learning, there was a significant increase in the score (p<0.001, ES = 0.333) but not in speed (p = 1.000, ES = 0.000). Thus, iVR as a medium proved to be an effective tool for learning geography. However, comparing the collaborative and individual learning showed no significant difference in the learning gain (p = 0.303, ES = 0.115), speed gain (p = 0.098, ES = 0.185), or performance motivation (p = 0.368, ES = 0.101). Nevertheless, the collaborative learning group had significantly higher use of cognitive resources (p = 0.046, ES = 0.223) than the individual learning group. The results were discussed in relation to the cognitive load theories, and future research directions for iVR learning were proposed.

Combining Deep Learning and Radiomics for Automated, Objective, Comprehensive Bone Marrow Characterization From Whole-Body MRI: A Multicentric Feasibility Study.

OBJECTIVES: Disseminated bone marrow (BM) involvement is frequent in multiple myeloma (MM). Whole-body magnetic resonance imaging (wb-MRI) enables to evaluate the whole BM. Reading of such whole-body scans is time-consuming, and yet radiologists can transfer only a small fraction of the information of the imaging data set to the report. This limits the influence that imaging can have on clinical decision-making and in research toward precision oncology. The objective of this feasibility study was to implement a concept for automatic, comprehensive characterization of the BM from wb-MRI, by automatic BM segmentation and subsequent radiomics analysis of 30 different BM spaces (BMS). MATERIALS AND METHODS: This retrospective multicentric pilot study used a total of 106 wb-MRI from 102 patients with (smoldering) MM from 8 centers. Fifty wb-MRI from center 1 were used for training of segmentation algorithms (nnU-Nets) and radiomics algorithms. Fifty-six wb-MRI from 8 centers, acquired with a variety of different MRI scanners and protocols, were used for independent testing. Manual segmentations of 2700 BMS from 90 wb-MRI were performed for training and testing of the segmentation algorithms. For each BMS, 296 radiomics features were calculated individually. Dice score was used to assess similarity between automatic segmentations and manual reference segmentations. RESULTS: The "multilabel nnU-Net" segmentation algorithm, which performs segmentation of 30 BMS and labels them individually, reached mean dice scores of 0.88 ± 0.06/0.87 ± 0.06/0.83 ± 0.11 in independent test sets from center 1/center 2/center 3-8 (interrater variability between radiologists, 0.88 ± 0.01). The subset from the multicenter, multivendor test set (center 3-8) that was of high imaging quality was segmented with high precision (mean dice score, 0.87), comparable to the internal test data from center 1. The radiomic BM phenotype consisting of 8880 descriptive parameters per patient, which result from calculation of 296 radiomics features for each of the 30 BMS, was calculated for all patients. Exemplary cases demonstrated connections between typical BM patterns in MM and radiomic signatures of the respective BMS. In plausibility tests, predicted size and weight based on radiomics models of the radiomic BM phenotype significantly correlated with patients' actual size and weight ( P = 0.002 and P = 0.003, respectively). CONCLUSIONS: This pilot study demonstrates the feasibility of automatic, objective, comprehensive BM characterization from wb-MRI in multicentric data sets. This concept allows the extraction of high-dimensional phenotypes to capture the complexity of disseminated BM disorders from imaging. Further studies need to assess the clinical potential of this method for automatic staging, therapy response assessment, or prediction of biopsy results.

Self-supervised pretraining for transferable quantitative phase image cell segmentation.

In this paper, a novel U-Net-based method for robust adherent cell segmentation for quantitative phase microscopy image is designed and optimised. We designed and evaluated four specific post-processing pipelines. To increase the transferability to different cell types, non-deep learning transfer with adjustable parameters is used in the post-processing step. Additionally, we proposed a self-supervised pretraining technique using nonlabelled data, which is trained to reconstruct multiple image distortions and improved the segmentation performance from 0.67 to 0.70 of object-wise intersection over union. Moreover, we publish a new dataset of manually labelled images suitable for this task together with the unlabelled data for self-supervised pretraining.

Parametric Deconvolution for Cancer Cells Viscoelasticity Measurements from Quantitative Phase Images.

In this contribution, we focused on optimising a dynamic flow-based shear stress system to achieve a reliable platform for cell shear modulus (stiffness) and viscosity assessment using quantitative phase imaging. The estimation of cell viscoelastic properties is influenced by distortion of the shear stress waveform, which is caused by the properties of the flow system components (i.e., syringe, flow chamber and tubing). We observed that these components have a significant influence on the measured cell viscoelastic characteristics. To suppress this effect, we applied a correction method utilizing parametric deconvolution of the flow system's optimized impulse response. Achieved results were compared with the direct fitting of the Kelvin-Voigt viscoelastic model and the basic steady-state model. The results showed that our novel parametric deconvolution approach is more robust and provides a more reliable estimation of viscosity with respect to changes in the syringe's compliance compared to Kelvin-Voigt model.

Learning-based vertebra localization and labeling in 3D CT data of possibly incomplete and pathological spines.

BACKGROUND AND OBJECTIVE: We present a fully automatic system based on learning approaches, which aims to localization and identification (labeling) of vertebrae in 3D computed tomography (CT) scans of possibly incomplete spines in patients with bone metastases and vertebral compressions. METHODS: The framework combines a set of 3D algorithms for i) spine detection using a convolution neural network (CNN) ii) spinal cord tracking based on combination of a CNN and a novel growing sphere method with a population optimization, iii) intervertebral discs localization using a novel approach of spatially variant filtering of intensity profiles and iv) vertebra labeling using a CNN-based classification combined with global dynamic optimization. RESULTS: The proposed algorithm has been validated in testing databases, including also a publicly available dataset. The mean error of intervertebral discs localization is 4.4 mm, and for vertebra labeling, the average rate of correctly identified vertebrae is 87.1%, which can be considered a good result with respect to the large share of highly distorted spines and incomplete spine scans. CONCLUSIONS: The proposed framework, which combines several advanced methods including also three CNNs, works fully automatically even with incomplete spine scans and with distorted pathological cases. The achieved results allow including the presented algorithms as the first phase to the fully automated computer-aided diagnosis (CAD) system for automatic spine-bone lesion analysis in oncological patients.

The 3D hype: Evaluating the potential of real 3D visualization in geo-related applications.

The use of 3D visualization technologies has increased rapidly in many applied fields, including geovisualization, and has been researched from many different perspectives. However, the findings for the benefits of 3D visualization, especially in stereoscopic 3D forms, remain inconclusive and disputed. Stereoscopic "real" 3D visualization was proposed as encouraging the visual perception of shapes and volume of displayed content yet criticised as problematic and limited in a number of ways, particularly in visual discomfort and increased response time in tasks. In order to assess the potential of real 3D visualization for geo-applications, 91 participants were engaged in this study to work with digital terrain models in different 3D settings. The researchers examined the effectivity of stereoscopic real 3D visualization compared to monoscopic 3D (or pseudo 3D) visualization under static and interactive conditions and applied three tasks with experimental stimuli representing different geo-related phenomena, i.e. objects in the terrain, flat areas marked in the terrain and terrain elevation profiles. The authors explored the significant effects of real 3D visualization and interactivity factors in terms of response time and correctness. Researchers observed that the option to interact (t = -10.849, p < 0.001) with a virtual terrain and its depiction with real 3D visualization (t = 4.64, p < 0.001) extended the participants' response times. Counterintuitively, the data demonstrated that the static condition increased response correctness (z = 5.38, p < 0.001). Regarding detailed analysis of data, an interactivity factor was proposed as a potential substitute for real 3D visualization in 3D geographical tasks.

Digital restoration of fragmentary human skeletal remains: Testing the feasibility of virtual reality.

Experts in forensic anthropology and medicine have become gradually accustomed to examining components of the human body in the virtual workspace. While the computer-assisted approach offers numerous benefits, the interactions with digital three-dimensional biological objects are often problematic, particularly if conducted with mouse, keyboard and flat-panel screen. The study focusses on feasibility of a virtual reality (VR) system for virtual restoration of fragmentary skeletal remains. The VR system was confronted with three cases of fragmentary remains. The cases were reassembled manually by twenty participants using a HTC Vive headset combined with an in-house application A.R.T. The same task was performed using a CloudCompare software in conjunction with a desktop peripheral. The two systems were compared in terms of time efficiency, the geometric properties of the resulting restorations, and convenience of use. Restoration using the VR system took approximately half the time the desktop set-up did. The VR system also yielded a lower error rate when a severely fragmented skull was reassembled. Ultimately, although the efficiency of the reassembling was shown to be strongly dependent on the operator's experience, the use of the VR system balanced out the uneven levels of proficiency in computer graphics. The current generation of virtual reality headsets has a strong potential to facilitate and improve tasks relating to the virtual restoration of fragmented skeletal remains. A VR system offers an intuitive digital working environment which is less affected by an operator's computer skills and practical understanding of the technology than the desktop systems are.

Deep-learning-based fully automatic spine centerline detection in CT data.

In this contribution, we present a fully automatic approach, that is based on two convolution neural networks (CNN) together with a spine tracing algorithm utilizing a population optimization algorithm. Based on the evaluation of 130 CT scans including heavily distorted and complicated cases, it turned out that this new combination enables fast and robust detection with almost 90% of correctly determined spinal centerlines with computing time of fewer than 20 seconds.

Iterative machine learning based rotational alignment of brain 3D CT data.

The optimal rotational alignment of brain Computed Tomography (CT) images to a required standard position has a crucial importance for both automatic and manual diagnostic analysis. In this contribution, we present a novel two-step iterative approach for the automatic 3D rotational alignment of brain CT data. The angles of axial and coronal rotations are determined by an unsupervised by localisation of the Midsagittal Plane (MSP) method. This includes detection and pairing of medially symmetrical feature points. The sagittal rotation angle is subsequently estimated by regression convolutional neural network (CNN). The proposed methodology has been evaluated on a dataset of CT data manually aligned by radiologists. It has been shown that the algorithm achieved the low error of estimated rotations (≈1 degree) and in a significantly shorter time than the experts (≈2 minutes per case).

Deep convolutional neural network-based segmentation and classification of difficult to define metastatic spinal lesions in 3D CT data.

This paper aims to address the segmentation and classification of lytic and sclerotic metastatic lesions that are difficult to define by using spinal 3D Computed Tomography (CT) images obtained from highly pathologically affected cases. As the lesions are ill-defined and consequently it is difficult to find relevant image features that would enable detection and classification of lesions by classical methods of texture and shape analysis, the problem is solved by automatic feature extraction provided by a deep Convolutional Neural Network (CNN). Our main contributions are: (i) individual CNN architecture, and pre-processing steps that are dependent on a patient data and a scan protocol - it enables work with different types of CT scans; (ii) medial axis transform (MAT) post-processing for shape simplification of segmented lesion candidates with Random Forest (RF) based meta-analysis; and (iii) usability of the proposed method on whole-spine CTs (cervical, thoracic, lumbar), which is not treated in other published methods (they work with thoracolumbar segments of spine only). Our proposed method has been tested on our own dataset annotated by two mutually independent radiologists and has been compared to other published methods. This work is part of the ongoing complex project dealing with spine analysis and spine lesion longitudinal studies.