Search

Background: Extraction of mandibular third molars (M3Ms) is a routine procedure in oral and maxillofacial surgery, often associated with postoperative symptoms like pain, facial swelling, and trismus. This study aimed to introduce a standardized and automated protocol for swelling analysis following M3M surgery, presenting results regarding clinical conditions immediately and one-week after surgery., Methods: In a prospective study, 35 patients were enrolled (mean age: 24.4 ± 5.8 years) for removal of 54 M3Ms. Facial swelling was evaluated through 3D facial scans before surgery (T0), at three days (T1), and seven days (T2) post-surgery. The open-source software 3DSlicer facilitated automated analysis, including data anonymization, orientation, surface registration, qualitative comparisons, linear measurements, and volumetric quantification. Pairwise superimposition of facial models enabled qualitative, vectorial, and quantitative assessments, comparing initial conditions with swelling development at T1 and T2. Additionally, changes between T1 and T2 were also evaluated. Secondary outcomes encompassed clinical evaluations of pain, trismus (maximum mouth opening), and surgery time. Statistical analysis involved the paired Student t-test to assess longitudinal changes and analysis of variance to evaluate outcome variables concerning difficulty scores. Linear regression models correlated primary outcome variables with secondary study variables (α < 0.05)., Results: Longitudinal analysis demonstrated significant but variable facial swelling, pain, and trismus at T1, followed by improvement at T2 (p < 0.001). Linear and volumetric differences correlated positively with surgery time (p < 0.05). A direct proportionality between linear and volume differences was observed, higher values at T1 correlated with higher values at T2 (p < 0.05)., Conclusions: An innovative digital workflow for precise quantification of postoperative facial changes was implemented, incorporating volumetric measurements that surpass linear assessments. Clinical conditions demonstrated a direct correlation with surgery time, deteriorating immediately and improving one-week after surgery., Competing Interests: Declarations Ethics approval and consent to participate The experimental protocol was approved by the ethical committee of Magna Graecia University of Catanzaro, Italy (n. 465/2020). Informed consent was diligently acquired from all subjects and/or their legal guardians. Consent for publication Not applicable. Competing interests The authors declare no competing interests., (© 2024. The Author(s).)

Background: Longitudinal Displacement (LD) is the relative motion of the intima-media upon adventitia of the arterial wall during the cardiac cycle, probably linked to atherosclerosis. It has a direction, physiologically first backward in its main components with respect to the arterial flow. Here, LD was investigated in various disease and in presence of a unilateral carotid stent., Methods: Carotid acquisitions were performed by ultrasound imaging on both body sides of 75 participants (150 Arteries). LD was measured in its percent quantity and direction., Results: Obesity (p = 0.001) and carotid plaques (p = 0.01) were independently associated to quantity decrease of LD in the whole population. In a subgroup analysis, it was respectively 143% in healthy (n = 48 carotids), 129% (n = 34) in presence of cardiovascular risk factors, 121% (n = 20) in MACE patients, 119% (n = 24) in the carotid contralateral to a stent, 110% (n = 24) in carotids with stents. Regarding the direction of LD, in a subgroup analysis an inverted movement was identified in aged (p = 0.001) and diseased (p = 0.001) participants who also showed less quantity of LD (p = 0.001), but independently with age only (p = 0.002) in the whole population., Conclusions: This observational study suggests that LD within carotid wall layers is lower additively with ageing, cardiovascular risk factors, cardiovascular diseases, and stent. Even if stent is surely beneficial, these data might shed some light on stent restenosis, emphasising the need for interventional studies.

Background: The global coronavirus disease 2019 (COVID-19) pandemic has posed substantial challenges for healthcare systems, notably the increased demand for chest computed tomography (CT) scans, which lack automated analysis. Our study addresses this by utilizing artificial intelligence-supported automated computer analysis to investigate lung involvement distribution and extent in COVID-19 patients. Additionally, we explore the association between lung involvement and intensive care unit (ICU) admission, while also comparing computer analysis performance with expert radiologists' assessments., Methods: A total of 81 patients from an open-source COVID database with confirmed COVID-19 infection were included in the study. Three patients were excluded. Lung involvement was assessed in 78 patients using CT scans, and the extent of infiltration and collapse was quantified across various lung lobes and regions. The associations between lung involvement and ICU admission were analysed. Additionally, the computer analysis of COVID-19 involvement was compared against a human rating provided by radiological experts., Results: The results showed a higher degree of infiltration and collapse in the lower lobes compared to the upper lobes (P<0.05). No significant difference was detected in the COVID-19-related involvement of the left and right lower lobes. The right middle lobe demonstrated lower involvement compared to the right lower lobes (P<0.05). When examining the regions, significantly more COVID-19 involvement was found when comparing the posterior vs . the anterior halves and the lower vs . the upper half of the lungs. Patients, who required ICU admission during their treatment exhibited significantly higher COVID-19 involvement in their lung parenchyma according to computer analysis, compared to patients who remained in general wards. Patients with more than 40% COVID-19 involvement were almost exclusively treated in intensive care. A high correlation was observed between computer detection of COVID-19 affections and the rating by radiological experts., Conclusions: The findings suggest that the extent of lung involvement, particularly in the lower lobes, dorsal lungs, and lower half of the lungs, may be associated with the need for ICU admission in patients with COVID-19. Computer analysis showed a high correlation with expert rating, highlighting its potential utility in clinical settings for assessing lung involvement. This information may help guide clinical decision-making and resource allocation during ongoing or future pandemics. Further studies with larger sample sizes are warranted to validate these findings., Competing Interests: Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://jtd.amegroups.com/article/view/10.21037/jtd-23-1150/coif). R.S.J.E. reports grants from the National Heart, Lung, and Blood Institute (NHLBI) for the present manuscript, with no time limit specified for this item; contract with Lung Biotechnology and Insmed to serve as Image Core for a study; sponsored research agreement with Boehringer Ingelheim; consulting fees from Leuko Labs and Mount Sinai; payment/honoraria Chiesi; patent pending in the area of lung cancer risk assessment using machine learning technology; he is co-founder and stockholder of Quantitative Imaging Solutions, a company specializing in imaging analytics in the lung cancer space. S.P. is a full-time employee of Isomics, Inc. and participated in the research as part of his normal work activities. The other authors have no conflicts of interest to declare., (2024 Journal of Thoracic Disease. All rights reserved.)

Purpose: Primary central nervous system lymphoma (PCNSL) is a rare, aggressive form of extranodal non-Hodgkin lymphoma. To predict the overall survival (OS) in advance is of utmost importance as it has the potential to aid clinical decision-making. Though radiomics-based machine learning (ML) has demonstrated the promising performance in PCNSL, it demands large amounts of manual feature extraction efforts from magnetic resonance images beforehand. deep learning (DL) overcomes this limitation., Methods: In this paper, we tailored the 3D ResNet to predict the OS of patients with PCNSL. To overcome the limitation of data sparsity, we introduced data augmentation and transfer learning, and we evaluated the results using r stratified k-fold cross-validation. To explain the results of our model, gradient-weighted class activation mapping was applied., Results: We obtained the best performance (the standard error) on post-contrast T1-weighted (T1Gd)-area under curve [Formula: see text], accuracy [Formula: see text], precision [Formula: see text], recall [Formula: see text] and F1-score [Formula: see text], while compared with ML-based models on clinical data and radiomics data, respectively, further confirming the stability of our model. Also, we observed that PCNSL is a whole-brain disease and in the cases where the OS is less than 1 year, it is more difficult to distinguish the tumor boundary from the normal part of the brain, which is consistent with the clinical outcome., Conclusions: All these findings indicate that T1Gd can improve prognosis predictions of patients with PCNSL. To the best of our knowledge, this is the first time to use DL to explain model patterns in OS classification of patients with PCNSL. Future work would involve collecting more data of patients with PCNSL, or additional retrospective studies on different patient populations with rare diseases, to further promote the clinical role of our model., (© 2023. The Author(s).)

A minimally-invasive manipulator characterized by hyper-redundant kinematics and embedded sensing modules is presented in this work. The bending angles (tilt and pan) of the robot tip are controlled through tendon-driven actuation; the transmission of the actuation forces to the tip is based on a Bowden-cable solution integrating some channels for optical fibers. The viability of the real-time measurement of the feedback control variables, through optoelectronic acquisition, is evaluated for automated bending of the flexible endoscope and trajectory tracking of the tip angles. Indeed, unlike conventional catheters and cannulae adopted in neurosurgery, the proposed robot can extend the actuation and control of snake-like kinematic chains with embedded sensing solutions, enabling real-time measurement, robust and accurate control of curvature, and tip bending of continuum robots for the manipulation of cannulae and microsurgical instruments in neurosurgical procedures. A prototype of the manipulator with a length of 43 mm and a diameter of 5.5 mm has been realized via 3D printing. Moreover, a multiple regression model has been estimated through a novel experimental setup to predict the tip angles from measured outputs of the optoelectronic modules. The sensing and control performance has also been evaluated during tasks involving tip rotations.

Clonogenic assays are routinely used to evaluate the response of cancer cells to external radiation fields, assess their radioresistance and radiosensitivity, estimate the performance of radiotherapy. However, classic clonogenic tests focus on the number of colonies forming on a substrate upon exposure to ionizing radiation, and disregard other important characteristics of cells such their ability to generate structures with a certain shape. The radioresistance and radiosensitivity of cancer cells may depend less on the number of cells in a colony and more on the way cells interact to form complex networks. In this study, we have examined whether the topology of 2D cancer-cell graphs is influenced by ionizing radiation. We subjected different cancer cell lines, i.e. H4 epithelial neuroglioma cells, H460 lung cancer cells, PC3 bone metastasis of grade IV of prostate cancer and T24 urinary bladder cancer cells, cultured on planar surfaces, to increasing photon radiation levels up to 6 Gy. Fluorescence images of samples were then processed to determine the topological parameters of the cell-graphs developing over time. We found that the larger the dose, the less uniform the distribution of cells on the substrate-evidenced by high values of small-world coefficient (cc), high values of clustering coefficient (cc), and small values of characteristic path length (cpl). For all considered cell lines, [Formula: see text] for doses higher or equal to 4 Gy, while the sensitivity to the dose varied for different cell lines: T24 cells seem more distinctly affected by the radiation, followed by the H4, H460 and PC3 cells. Results of the work reinforce the view that the characteristics of cancer cells and their response to radiotherapy can be determined by examining their collective behavior-encoded in a few topological parameters-as an alternative to classical clonogenic assays., (© 2023. The Author(s).)

Primary Central Nervous System Lymphoma (PCNSL) is an aggressive neoplasm with a poor prognosis. Although therapeutic progresses have significantly improved Overall Survival (OS), a number of patients do not respond to HD-MTX-based chemotherapy (15-25%) or experience relapse (25-50%) after an initial response. The reasons underlying this poor response to therapy are unknown. Thus, there is an urgent need to develop improved predictive models for PCNSL. In this study, we investigated whether radiomics features can improve outcome prediction in patients with PCNSL. A total of 80 patients diagnosed with PCNSL were enrolled. A patient sub-group, with complete Magnetic Resonance Imaging (MRI) series, were selected for the stratification analysis. Following radiomics feature extraction and selection, different Machine Learning (ML) models were tested for OS and Progression-free Survival (PFS) prediction. To assess the stability of the selected features, images from 23 patients scanned at three different time points were used to compute the Interclass Correlation Coefficient (ICC) and to evaluate the reproducibility of each feature for both original and normalized images. Features extracted from Z-score normalized images were significantly more stable than those extracted from non-normalized images with an improvement of about 38% on average ( p -value < 10-12). The area under the ROC curve (AUC) showed that radiomics-based prediction overcame prediction based on current clinical prognostic factors with an improvement of 23% for OS and 50% for PFS, respectively. These results indicate that radiomics features extracted from normalized MR images can improve prognosis stratification of PCNSL patients and pave the way for further study on its potential role to drive treatment choice.

Background: Time-resolved 4D cone beam-computed tomography (4D-CBCT) allows a daily assessment of patient anatomy and respiratory motion. However, 4D-CBCTs suffer from imaging artifacts that affect the CT number accuracy and prevent accurate proton dose calculations. Deep learning can be used to correct CT numbers and generate synthetic CTs (sCTs) that can enable CBCT-based proton dose calculations., Purpose: In this work, sparse view 4D-CBCTs were converted into 4D-sCT utilizing a deep convolutional neural network (DCNN). 4D-sCTs were evaluated in terms of image quality and dosimetric accuracy to determine if accurate proton dose calculations for adaptive proton therapy workflows of lung cancer patients are feasible., Methods: A dataset of 45 thoracic cancer patients was utilized to train and evaluate a DCNN to generate 4D-sCTs, based on sparse view 4D-CBCTs reconstructed from projections acquired with a 3D acquisition protocol. Mean absolute error (MAE) and mean error were used as metrics to evaluate the image quality of single phases and average 4D-sCTs against 4D-CTs acquired on the same day. The dosimetric accuracy was checked globally (gamma analysis) and locally for target volumes and organs-at-risk (OARs) (lung, heart, and esophagus). Furthermore, 4D-sCTs were also compared to 3D-sCTs. To evaluate CT number accuracy, proton radiography simulations in 4D-sCT and 4D-CTs were compared in terms of range errors. The clinical suitability of 4D-sCTs was demonstrated by performing a 4D dose reconstruction using patient specific treatment delivery log files and breathing signals., Results: 4D-sCTs resulted in average MAEs of 48.1 ± 6.5 HU (single phase) and 37.7 ± 6.2 HU (average). The global dosimetric evaluation showed gamma pass ratios of 92.3% ± 3.2% (single phase) and 94.4% ± 2.1% (average). The clinical target volume showed high agreement in D 98 between 4D-CT and 4D-sCT, with differences below 2.4% for all patients. Larger dose differences were observed in mean doses of OARs (up to 8.4%). The comparison with 3D-sCTs showed no substantial image quality and dosimetric differences for the 4D-sCT average. Individual 4D-sCT phases showed slightly lower dosimetric accuracy. The range error evaluation revealed that lung tissues cause range errors about three times higher than the other tissues., Conclusion: In this study, we have investigated the accuracy of deep learning-based 4D-sCTs for daily dose calculations in adaptive proton therapy. Despite image quality differences between 4D-sCTs and 3D-sCTs, comparable dosimetric accuracy was observed globally and locally. Further improvement of 3D and 4D lung sCTs could be achieved by increasing CT number accuracy in lung tissues., (© 2022 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

Radiation therapy (RT) is now considered to be a main component of cancer therapy, alongside surgery, chemotherapy and monoclonal antibody-based immunotherapy. In RT, cancer tissues are exposed to ionizing radiation causing the death of malignant cells and favoring cancer regression. However, the efficiency of RT may be hampered by cell-radioresistance (RR)-that is a feature of tumor cells of withstanding RT. To improve the RT performance, it is decisive developing methods that can help to quantify cell sensitivity to radiation. In acknowledgment of the fact that none of the existing methods to assess RR are based on cell graphs topology, in this work we have examined how 2D cell networks, within a single colony, from different human lung cancer lines (H460, A549 and Calu-1) behave in response to doses of ionizing radiation ranging from 0 to 8 Gy. We measured the structure of resulting cell-graphs using well-assessed networks-analysis metrics, such as the clustering coefficient (cc), the characteristic path length (cpl), and the small world coefficient (SW). Findings of the work illustrate that the clustering characteristics of cell-networks show a marked sensitivity to the dose and cell line. Higher-than-one values of SW coefficient, clue of a discontinuous and inhomogeneous cell spatial layout, are associated to elevated levels of radiation and to a lower radio-resistance of the treated cell line. Results of the work suggest that topology could be used as a quantitative parameter to assess the cell radio-resistance and measure the performance of cancer radiotherapy., (© 2022. The Author(s).)

Coronary Angiography (CA) is the standard of reference to diagnose coronary artery disease. Yet, only a portion of the information it conveys is usually used. Quantitative Coronary Angiography (QCA) reliably contributes to improving the measurable assessment of CA. In this work, we developed a new software, CoroFinder, able to automatically identify epicardial coronary arteries and to dynamically track the vessel profile in dye-free frames. The coronary tree is automatically segmented by Frangi's filter in the angiogram's frames where vessels are contrasted ("template frames"). Afterward, the image similarity among each template frame and the dye-free images is scored by cross-correlation. Finally, each dye-free image is associated with the most similar template frame, resulting in an estimation of vessel contour. CoroFinder allows locating the position of coronary arteries in absence of contrast dye. The developed algorithm is robust to diverse vessel curvatures, variation of vessel widths, and the presence of stenoses. This article describes the newly developed CoroFinder algorithm and the associated software and provides an overview of its potential application in research and for translation to the clinic.

Microarray is a powerful technology that enables the monitoring of expression levels for thousands of genes simultaneously, providing scientists with a full overview about DNA and RNA investigation. The process is made of three main phases: interaction with biological samples, data extraction, and data analysis. In particular, the data extraction phase strongly relies on image processing algorithms, since the expression levels are revealed by the interaction of light with fluorescent markers. More in detail, in order to extract quantitative information from probes image, three steps are required: (1) gridding, (2) segmentation, and (3) intensity quantification. Errors in one of these steps can deeply affect the process outcome. In this chapter each of the above mentioned steps will be analyzed and discussed. Software platforms dedicated to this purpose will be reported as well., (© 2022. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Purpose: Adaptive proton therapy (APT) of lung cancer patients requires frequent volumetric imaging of diagnostic quality. Cone-beam CT (CBCT) can provide these daily images, but x-ray scattering limits CBCT-image quality and hampers dose calculation accuracy. The purpose of this study was to generate CBCT-based synthetic CTs using a deep convolutional neural network (DCNN) and investigate image quality and clinical suitability for proton dose calculations in lung cancer patients., Methods: A dataset of 33 thoracic cancer patients, containing CBCTs, same-day repeat CTs (rCT), planning-CTs (pCTs), and clinical proton treatment plans, was used to train and evaluate a DCNN with and without a pCT-based correction method. Mean absolute error (MAE), mean error (ME), peak signal-to-noise ratio, and structural similarity were used to quantify image quality. The evaluation of clinical suitability was based on recalculation of clinical proton treatment plans. Gamma pass ratios, mean dose to target volumes and organs at risk, and normal tissue complication probabilities (NTCP) were calculated. Furthermore, proton radiography simulations were performed to assess the HU-accuracy of sCTs in terms of range errors., Results: On average, sCTs without correction resulted in a MAE of 34 ± 6 HU and ME of 4 ± 8 HU. The correction reduced the MAE to 31 ± 4HU (ME to 2 ± 4HU). Average 3%/3 mm gamma pass ratios increased from 93.7% to 96.8%, when the correction was applied. The patient specific correction reduced mean proton range errors from 1.5 to 1.1 mm. Relative mean target dose differences between sCTs and rCT were below ± 0.5% for all patients and both synthetic CTs (with/without correction). NTCP values showed high agreement between sCTs and rCT (<2%)., Conclusion: CBCT-based sCTs can enable accurate proton dose calculations for APT of lung cancer patients. The patient specific correction method increased the image quality and dosimetric accuracy but had only a limited influence on clinically relevant parameters., (© 2021 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

Recently,deep learning (DL)-based methods for the generation of synthetic computed tomography (sCT) have received significant research attention as an alternative to classical ones. We present here a systematic review of these methods by grouping them into three categories, according to their clinical applications: (i) to replace computed tomography in magnetic resonance (MR) based treatment planning, (ii) facilitate cone-beam computed tomography based image-guided adaptive radiotherapy, and (iii) derive attenuation maps for the correction of positron emission tomography. Appropriate database searching was performed on journal articles published between January 2014 and December 2020. The DL methods' key characteristics were extracted from each eligible study, and a comprehensive comparison among network architectures and metrics was reported. A detailed review of each category was given, highlighting essential contributions, identifying specific challenges, and summarizing the achievements. Lastly, the statistics of all the cited works from various aspects were analyzed, revealing the popularity and future trends and the potential of DL-based sCT generation. The current status of DL-based sCT generation was evaluated, assessing the clinical readiness of the presented methods., (© 2021 The Authors. Medical Physics published by Wiley Periodicals LLC on behalf of American Association of Physicists in Medicine.)

Breast cancer is the most frequent cancer in women worldwide and late diagnosis often adversely affects the prognosis of the disease. Radiotherapy is commonly used to treat breast cancer, reducing the risk of recurrence after surgery. However, the eradication of radioresistant cancer cells, including cancer stem cells, remains the main challenge of radiotherapy. Recently, lipid droplets (LDs) have been proposed as functional markers of cancer stem cells, also being involved in increased cell tumorigenicity. LD biogenesis is a multistep process requiring various enzymes, including Diacylglycerol acyltransferase 2 (DGAT2). In this context, we evaluated the effect of PF-06424439, a selective DGAT2 inhibitor, on MCF7 breast cancer cells exposed to X-rays. Our results demonstrated that 72 h of PF-06424439 treatment reduced LD content and inhibited cell migration, without affecting cell proliferation. Interestingly, PF-06424439 pre-treatment followed by radiation was able to enhance radiosensitivity of MCF7 cells. In addition, the combined treatment negatively interfered with lipid metabolism-related genes, as well as with EMT gene expression, and modulated the expression of typical markers associated with the CSC-like phenotype. These findings suggest that PF-06424439 pre-treatment coupled to X-ray exposure might potentiate breast cancer cell radiosensitivity and potentially improve the radiotherapy effectiveness.

The coronavirus disease 19 (COVID-19) pandemic is having a dramatic impact on society and healthcare systems. In this complex scenario, lung computerized tomography (CT) may play an important prognostic role. However, datasets released so far present limitations that hamper the development of tools for quantitative analysis. In this paper, we present an open-source lung CT dataset comprising information on 50 COVID-19-positive patients. The CT volumes are provided along with (i) an automatic threshold-based annotation obtained with a Gaussian mixture model (GMM) and (ii) a scoring provided by an expert radiologist. This score was found to significantly correlate with the presence of ground glass opacities and the consolidation found with GMM. The dataset is freely available in an ITK-based file format under the CC BY-NC 4.0 license. The code for GMM fitting is publicly available, as well. We believe that our dataset will provide a unique opportunity for researchers working in the field of medical image analysis, and hope that its release will lay the foundations for the successfully implementation of algorithms to support clinicians in facing the COVID-19 pandemic.

Cone-beam computed tomography (CBCT)- and magnetic resonance (MR)-images allow a daily observation of patient anatomy but are not directly suited for accurate proton dose calculations. This can be overcome by creating synthetic CTs (sCT) using deep convolutional neural networks. In this study, we compared sCTs based on CBCTs and MRs for head and neck (H&N) cancer patients in terms of image quality and proton dose calculation accuracy. A dataset of 27 H&N-patients, treated with proton therapy (PT), containing planning CTs (pCTs), repeat CTs, CBCTs and MRs were used to train two neural networks to convert either CBCTs or MRs into sCTs. Image quality was quantified by calculating mean absolute error (MAE), mean error (ME) and Dice similarity coefficient (DSC) for bones. The dose evaluation consisted of a systematic non-clinical analysis and a clinical recalculation of actually used proton treatment plans. Gamma analysis was performed for non-clinical and clinical treatment plans. For clinical treatment plans also dose to targets and organs at risk (OARs) and normal tissue complication probabilities (NTCP) were compared. CBCT-based sCTs resulted in higher image quality with an average MAE of 40 ± 4 HU and a DSC of 0.95, while for MR-based sCTs a MAE of 65 ± 4 HU and a DSC of 0.89 was observed. Also in clinical proton dose calculations, sCT CBCT achieved higher average gamma pass ratios (2%/2 mm criteria) than sCT MR (96.1% vs. 93.3%). Dose-volume histograms for selected OARs and NTCP-values showed a very small difference between sCT CBCT and sCT MR and a high agreement with the reference pCT. CBCT- and MR-based sCTs have the potential to enable accurate proton dose calculations valuable for daily adaptive PT. Significant image quality differences were observed but did not affect proton dose calculation accuracy in a similar manner. Especially the recalculation of clinical treatment plans showed high agreement with the pCT for both sCT CBCT and sCT MR.

Interaction between medical image platform and external environment is a desirable feature in several clinical, research, and educational scenarios. In this work, the integration between 3D Slicer package and Arduino board is introduced, enabling a simple and useful communication between the two software/hardware platforms. The open source extension, programmed in Python language, manages the connection process and offers a communication layer accessible from any point of the medical image suite infrastructure. Deep integration with 3D Slicer code environment is provided and a basic input-output mechanism accessible via GUI is also made available. To test the proposed extension, two exemplary use cases were implemented: (1) INPUT data to 3D Slicer, to navigate on basis of data detected by a distance sensor connected to the board, and (2) OUTPUT data from 3D Slicer, to control a servomotor on the basis of data computed through image process procedures. Both goals were achieved and quasi-real-time control was obtained without any lag or freeze, thus boosting the integration between 3D Slicer and Arduino. This integration can be easily obtained through the execution of few lines of Python code. In conclusion, SlicerArduino proved to be suitable for fast prototyping, basic input-output interaction, and educational purposes. The extension is not intended for mission-critical clinical tasks.

With the advent of Minimally Invasive Surgery (MIS), intra-operative imaging has become crucial for surgery and therapy guidance, allowing to partially compensate for the lack of information typical of MIS. This paper reviews the advancements in both classical (i.e. ultrasounds, X-ray, optical coherence tomography and magnetic resonance imaging) and more recent (i.e. multispectral, photoacoustic and Raman imaging) intra-operative imaging modalities. Each imaging modality was analyzed, focusing on benefits and disadvantages in terms of compatibility with the operating room, costs, acquisition time and image characteristics. Tables are included to summarize this information. New generation of hybrid surgical room and algorithms for real time/in room image processing were also investigated. Each imaging modality has its own (site- and procedure-specific) peculiarities in terms of spatial and temporal resolution, field of view and contrasted tissues. Besides the benefits that each technique offers for guidance, considerations about operators and patient risk, costs, and extra time required for surgical procedures have to be considered. The current trend is to equip surgical rooms with multimodal imaging systems, so as to integrate multiple information for real-time data extraction and computer-assisted processing. The future of surgery is to enhance surgeons eye to minimize intra- and after-surgery adverse events and provide surgeons with all possible support to objectify and optimize the care-delivery process.

Background: The targeting accuracy of proton therapy (PT) for moving soft-tissue tumours is expected to greatly improve by real-time magnetic resonance imaging (MRI) guidance. The integration of MRI and PT at the treatment isocenter would offer the opportunity of combining the unparalleled soft-tissue contrast and real-time imaging capabilities of MRI with the most conformal dose distribution and best dose steering capability provided by modern PT. However, hybrid systems for MR-integrated PT (MRiPT) have not been realized so far due to a number of hitherto open technological challenges. In recent years, various research groups have started addressing these challenges and exploring the technical feasibility and clinical potential of MRiPT. The aim of this contribution is to review the different aspects of MRiPT, to report on the status quo and to identify important future research topics., Methods: Four aspects currently under study and their future directions are discussed: modelling and experimental investigations of electromagnetic interactions between the MRI and PT systems, integration of MRiPT workflows in clinical facilities, proton dose calculation algorithms in magnetic fields, and MRI-only based proton treatment planning approaches., Conclusions: Although MRiPT is still in its infancy, significant progress on all four aspects has been made, showing promising results that justify further efforts for research and development to be undertaken. First non-clinical research solutions have recently been realized and are being thoroughly characterized. The prospect that first prototype MRiPT systems for clinical use will likely exist within the next 5 to 10 years seems realistic, but requires significant work to be performed by collaborative efforts of research groups and industrial partners.

In-room imaging is a prerequisite for adaptive proton therapy. The use of onboard cone-beam computed tomography (CBCT) imaging, which is routinely acquired for patient position verification, can enable daily dose reconstructions and plan adaptation decisions. Image quality deficiencies though, hamper dose calculation accuracy and make corrections of CBCTs a necessity. This study compared three methods to correct CBCTs and create synthetic CTs that are suitable for proton dose calculations. CBCTs, planning CTs and repeated CTs (rCT) from 33 H&N cancer patients were used to compare a deep convolutional neural network (DCNN), deformable image registration (DIR) and an analytical image-based correction method (AIC) for synthetic CT (sCT) generation. Image quality of sCTs was evaluated by comparison with a same-day rCT, using mean absolute error (MAE), mean error (ME), Dice similarity coefficient (DSC), structural non-uniformity (SNU) and signal/contrast-to-noise ratios (SNR/CNR) as metrics. Dosimetric accuracy was investigated in an intracranial setting by performing gamma analysis and calculating range shifts. Neural network-based sCTs resulted in the lowest MAE and ME (37/2 HU) and the highest DSC (0.96). While DIR and AIC generated images with a MAE of 44/77 HU, a ME of -8/1 HU and a DSC of 0.94/0.90. Gamma and range shift analysis showed almost no dosimetric difference between DCNN and DIR based sCTs. The lower image quality of AIC based sCTs affected dosimetric accuracy and resulted in lower pass ratios and higher range shifts. Patient-specific differences highlighted the advantages and disadvantages of each method. For the set of patients, the DCNN created synthetic CTs with the highest image quality. Accurate proton dose calculations were achieved by both DCNN and DIR based sCTs. The AIC method resulted in lower image quality and dose calculation accuracy was reduced compared to the other methods.

The intricate relationships between innate immunity and brain diseases raise increased interest across the wide spectrum of neurodegenerative and neuropsychiatric disorders. Barriers, such as the blood-brain barrier, and innate immunity cells such as microglia, astrocytes, macrophages, and mast cells are involved in triggering disease events in these groups, through the action of many different cytokines. Chronic inflammation can lead to dysfunctions in large-scale brain networks. Neurodegenerative diseases, such as Alzheimer's disease, Parkinson's disease, Huntington's disease, amyotrophic lateral sclerosis, and frontotemporal dementia, are associated with a substrate of dysregulated immune responses that impair the central nervous system balance. Recent evidence suggests that similar phenomena are involved in psychiatric diseases, such as depression, schizophrenia, autism spectrum disorders, and post-traumatic stress disorder. The present review summarizes and discusses the main evidence linking the innate immunological response in neurodegenerative and psychiatric diseases, thus providing insights into how the responses of innate immunity represent a common denominator between diseases belonging to the neurological and psychiatric sphere. Improved knowledge of such immunological aspects could provide the framework for the future development of new diagnostic and therapeutic approaches.

Purpose: The first aim of this work is to present a novel deep convolution neural network (DCNN) multiplane approach and compare it to single-plane prediction of synthetic computed tomography (sCT) by using the real computed tomography (CT) as ground truth. The second aim is to demonstrate the feasibility of magnetic resonance imaging (MRI)-based proton therapy planning for the brain by assessing the range shift error within the clinical acceptance threshold., Methods and Materials: The image database included 15 pairs of MRI/CT scans of the head. Three DCNNs were trained to estimate, for each voxel, the Hounsfield unit (HU) value from MRI intensities. Each DCNN gave an estimation in the axial, sagittal, and coronal plane, respectively. The median HU among the 3 values was selected to build the sCT. The sCT/CT agreement was evaluated by a mean absolute error (MAE) and mean error, computed within the head contour and on 6 different tissues. Dice similarity coefficients were calculated to assess the geometric overlap of bone and air cavities segmentations. A 3-beam proton therapy plan was simulated for each patient. Beam-by-beam range shift (RS) analysis was conducted to assess the proton-stopping power estimation. RS analysis was performed using clinically accepted thresholds of (1) 3.5% + 1 mm and (2) 2.5% + 1.5 mm of the total range., Results: DCNN multiplane statistically outperformed single-plane prediction of sCT (P < .025). MAE and mean error within the head were 54 ± 7 HU and -4 ± 17 HU (mean ± standard deviation), respectively. Soft tissues were very close to perfect agreement (11 ± 3 HU in terms of MAE). Segmentation of air and bone regions led to a Dice similarity coefficient of 0.92 ± 0.03 and 0.93 ± 0.02, respectively. Proton RS was always below clinical acceptance thresholds, with a relative RS error of 0.14% ± 1.11%., Conclusions: The multiplane DCNN approach significantly improved the sCT prediction compared with other DCNN methods presented in the literature. The method was demonstrated to be highly accurate for MRI-only proton planning purposes., (Copyright © 2019 Elsevier Inc. All rights reserved.)

The data of literature are discordant about the role of mast cells in different types of neoplasms. In this paper the authors propose the hypothesis that tumor-associated mast cells may switch to different polarization states, conditioning the immunogenic capacities of the different neoplasms. Anti-inflammatory polarized mast cells should express cytokines such as interleukin-10 (IL-10) and then mast cells number should be inversely related to the intensity of inflammatory infiltrate. On the contrary, when mast cells do not express anti-inflammatory cytokines their number should be directly related to the intensity of the inflammatory infiltrate. In this paper we briefly argue around feasible approaches, based on the retrospective studies of tumor tissue samples from neoplasms considered "immunologically hot" and neoplasms considered "immunologically cold", through immunohistochemistry and immunofluorescence techniques (confocal microscopy). The establishment of the actual existence of a polarization interchange of mast cells, could lead to a new vision in prognostic terms, useful to contrive new approaches in immunotherapy of tumors., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Purpose: In radiation therapy, a key step for a successful cancer treatment is image-based treatment planning. One objective of the planning phase is the fast and accurate segmentation of organs at risk and target structures from medical images. However, manual delineation of organs, which is still the gold standard in many clinical environments, is time-consuming and prone to inter-observer variations. Consequently, many automated segmentation methods have been developed., Methods: In this work, we train two hierarchical 3D neural networks to segment multiple organs at risk in the head and neck area. First, we train a coarse network on size-reduced medical images to locate the organs of interest. Second, a subsequent fine network on full-resolution images is trained for a final accurate segmentation. The proposed method is purely deep learning based; accordingly, no pre-registration or post-processing is required., Results: The approach has been applied on a publicly available computed tomography dataset, created for the MICCAI 2015 Auto-Segmentation challenge. In an extensive evaluation process, the best configurations for the trained networks have been determined. Compared to the existing methods, the presented approach shows state-of-the-art performance for the segmentation of seven different structures in the head and neck area., Conclusion: We conclude that 3D neural networks outperform the most existing model- and atlas-based methods for the segmentation of organs at risk in the head and neck area. The ease of use, high accuracy and the test time efficiency of the method make it promising for image-based treatment planning in clinical practice.

Recent studies have clarified many still unknown aspects related to innate immunity and the blood-brain barrier relationship. They have also confirmed the close links between effector immune system cells, such as granulocytes, macrophages, microglia, natural killer cells and mast cells, and barrier functionality. The latter, in turn, is able to influence not only the entry of the cells of the immune system into the nervous tissue, but also their own activation. Interestingly, these two components and their interactions play a role of great importance not only in infectious diseases, but in almost all the pathologies of the central nervous system. In this paper, we review the main aspects in the field of vascular diseases (cerebral ischemia), of primitive and secondary neoplasms of Central Nervous System CNS, of CNS infectious diseases, of most common neurodegenerative diseases, in epilepsy and in demyelinating diseases (multiple sclerosis). Neuroinflammation phenomena are constantly present in all diseases; in every different pathological state, a variety of innate immunity cells responds to specific stimuli, differentiating their action, which can influence the blood-brain barrier permeability. This, in turn, undergoes anatomical and functional modifications, allowing the stabilization or the progression of the pathological processes., Competing Interests: The authors declare no conflict of interest.

Purpose: To investigate advanced multimodal methods for pseudo-computed tomography (CT) generation from standard magnetic resonance imaging sequences and to validate the results by intensity-modulated radiation therapy (IMRT) and volumetric modulated arc therapy (VMAT) plans. We present 2 novel methods that employ key techniques to enhance pseudo-CTs and investigate the effect on image quality and applicability for IMRT and VMAT planning., Materials and Methods: The data set contains CT and magnetic resonance image scans from 15 patients who underwent cranial radiation therapy. For each patient, pseudo-CTs of the head were generated with a patch-based and a voxel-based algorithm. The accuracy of the pseudo-CTs in comparison to clinical CTs was evaluated by mean absolute error, bias, and the Dice coefficient (of bone). IMRT and VMAT plans were created for each patient. Dose distributions were calculated with both the pseudo-CT and the clinical CT scans and compared by gamma tests, dose-volume histograms, and isocenter doses., Results: The generated pseudo-CTs exhibited average mean absolute errors of 118.7 ± 10.4 HU for the voxel-based algorithm and 73.0 ± 6.4 HU for the patch-based algorithm. The dose calculations were in good agreement and showed gamma test (2 mm, 2%) pass rates for both beam setups (IMRT and VMAT) of over 99% for 14 patients and over 98% for 1 patient., Conclusions: We showed that the key techniques of our 2 novel algorithms advance the quality of pseudo-CT significantly and generate very competitive pseudo-CTs compared with previously published methods. This quality was confirmed by low dose error in comparison to the ground-truth CT. With the achieved level of accuracy, our patch-based algorithm especially is a candidate for clinical routine use in IMRT and VMAT planning., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Background: In radiotherapy, MR imaging is only used because it has significantly better soft tissue contrast than CT, but it lacks electron density information needed for dose calculation. This work assesses the feasibility of using pseudo-CT (pCT) generated from T1w/T2w MR for proton treatment planning, where proton range comparisons are performed between standard CT and pCT., Material and Methods: MR and CT data from 14 glioblastoma patients were used in this study. The pCT was generated by using conversion libraries obtained from tissue segmentation and anatomical regioning of the T1w/T2w MR. For each patient, a plan consisting of three 18 Gy beams was designed on the pCT, for a total of 42 analyzed beams. The plan was then transferred onto the CT that represented the ground truth. Range shift (RS) between pCT and CT was computed at R 80 over 10 slices. The acceptance threshold for RS was according to clinical guidelines of two institutions. A γ-index test was also performed on the total dose for each patient., Results: Mean absolute error and bias for the pCT were 124 ± 10 and -16 ± 26 Hounsfield Units (HU), respectively. The median and interquartile range of RS was 0.5 and 1.4 mm, with highest absolute value being 4.4 mm. Of the 42 beams, 40 showed RS less than the clinical range margin. The two beams with larger RS were both in the cranio-caudal direction and had segmentation errors due to the partial volume effect, leading to misassignment of the HU., Conclusions: This study showed the feasibility of using T1w and T2w MRI to generate a pCT for proton therapy treatment, thus avoiding the use of a planning CT and allowing better target definition and possibilities for online adaptive therapies. Further improvements of the methodology are still required to improve the conversion from MRI intensities to HUs.

Multi atlas based segmentation (MABS) uses a database of atlas images, and an atlas selection process is used to choose an atlas subset for registration and voting. In the current state of the art, atlases are chosen according to a similarity criterion between the target subject and each atlas in the database. In this paper, we propose a new concept for atlas selection that relies on selecting the best performing group of atlases rather than the group of highest scoring individual atlases. Experiments were performed using CT images of 50 patients, with contours of brainstem and parotid glands. The dataset was randomly split into two groups: 20 volumes were used as an atlas database and 30 served as target subjects for testing. Classic oracle selection, where atlases are chosen by the highest dice similarity coefficient (DSC) with the target, was performed. This was compared to oracle group selection, where all the combinations of atlas subgroups were considered and scored by computing DSC with the target subject. Subsequently, convolutional neural networks were designed to predict the best group of atlases. The results were also compared with the selection strategy based on normalized mutual information (NMI). Oracle group was proven to be significantly better than classic oracle selection (p  <  10 -5 ). Atlas group selection led to a median  ±  interquartile DSC of 0.740  ±  0.084, 0.718  ±  0.086 and 0.670  ±  0.097 for brainstem and left/right parotid glands respectively, outperforming NMI selection 0.676  ±  0.113, 0.632  ±  0.104 and 0.606  ±  0.118 (p  <  0.001) as well as classic oracle selection. The implemented methodology is a proof of principle that selecting the atlases by considering the performance of the entire group of atlases instead of each single atlas leads to higher segmentation accuracy, being even better then current oracle strategy. This finding opens a new discussion about the most appropriate atlas selection criterion for MABS.

The purpose of this work is to present and validate a novel approach for ultra-sound-based speckle tracking to measure the carotid artery longitudinal displacement, and to assess the apparent sliding between of Intima-Media Complex (IMC) and Adventitia (Ad) layers. This method utilizes feature detectors and descriptors to localize and track keypoints for local motion quantification. The procedure was tested and validated on an in silico dataset and on 18 heathy volunteers and 16 patients. Accuracy measured on in silico data gave a mean ± standard deviation of 23 ± 15 and 19 ± 18 μm for IMC and Ad respectively, and thus smaller than the pixel size (0.0925 mm). Robustness analysis was performed on in vivo images, obtaining a maximum variation coefficient, over 5 repeated measures, of 9.5 and 13.8% for IMC and Ad, respectively. The novel method capability for detecting the relative motion of IMC vs. Ad was compared with visual assessment performed by 2 physicians, leading to a correlation coefficient R of 0.7 in the worst case. (Healthy group scored by rater #1.) In conclusion, our results provide evidence that the novel method is able to accurately and reliably track carotid artery layer motion and that it overcomes limitations currently present in the literature, therefore providing an automatic tool for clinical evaluation of IMC vs. Ad relative displacement.

Purpose: Automated delineation of structures and organs is a key step in medical imaging. However, due to the large number and diversity of structures and the large variety of segmentation algorithms, a consensus is lacking as to which automated segmentation method works best for certain applications. Segmentation challenges are a good approach for unbiased evaluation and comparison of segmentation algorithms., Methods: In this work, we describe and present the results of the Head and Neck Auto-Segmentation Challenge 2015, a satellite event at the Medical Image Computing and Computer Assisted Interventions (MICCAI) 2015 conference. Six teams participated in a challenge to segment nine structures in the head and neck region of CT images: brainstem, mandible, chiasm, bilateral optic nerves, bilateral parotid glands, and bilateral submandibular glands., Results: This paper presents the quantitative results of this challenge using multiple established error metrics and a well-defined ranking system. The strengths and weaknesses of the different auto-segmentation approaches are analyzed and discussed., Conclusions: The Head and Neck Auto-Segmentation Challenge 2015 was a good opportunity to assess the current state-of-the-art in segmentation of organs at risk for radiotherapy treatment. Participating teams had the possibility to compare their approaches to other methods under unbiased and standardized circumstances. The results demonstrate a clear tendency toward more general purpose and fewer structure-specific segmentation algorithms., (© 2017 American Association of Physicists in Medicine.)

Objectives: Atlas-based automatic segmentation (ABAS) addresses the challenges of accuracy and reliability in manual segmentation. We aim to evaluate the contribution of specific-purpose in ABAS of breast cancer (BC) patients with respect to generic-purpose libraries., Materials and Methods: One generic-purpose and 9 specific-purpose libraries, stratified according to type of surgery and size of thorax circumference, were obtained from the computed tomography of 200 BC patients. Keywords about contralateral breast volume and presence of breast expander/prostheses were recorded. ABAS was validated on 47 independent patients, considering manual segmentation from scratch as reference. Five ABAS datasets were obtained, testing single-ABAS and multi-ABAS with simultaneous truth and performance level estimation (STAPLE). Center of mass distance (CMD), average Hausdorff distance (AHD) and Dice similarity coefficient (DSC) between corresponding ABAS and manual structures were evaluated and statistically significant differences between different surgeries, structures and ABAS strategies were investigated., Results: Statistically significant differences between patients who underwent different surgery were found, with superior results for conservative-surgery group, and between different structures were observed: ABAS of heart, lungs, kidneys and liver was satisfactory (median values: CMD<2 mm, DSC≥0.80, AHD<1.5 mm), whereas chest wall, breast and spinal cord obtained moderate performance (median values: 2 mm ≤ CMD<5 mm, 0.60 ≤ DSC<0.80, 1.5 mm ≤ AHD<4 mm) and esophagus, stomach, brachial plexus and supraclavicular nodes obtained poor performance (median CMD≥5 mm, DSC<0.60, AHD≥4 mm). The application of STAPLE algorithm generally yields higher performance and the use of keywords improves results for breast ABAS., Conclusion: The homogeneity in the selection of atlases based on multiple anatomical and clinical features and the use of specific-purpose libraries can improve ABAS performance with respect to generic-purpose libraries., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Purpose: Multiatlas based segmentation is largely used in many clinical and research applications. Due to its good performances, it has recently been included in some commercial platforms for radiotherapy planning and surgery guidance. Anyway, to date, a software with no restrictions about the anatomical district and image modality is still missing. In this paper we introduce plastimatch mabs, an open source software that can be used with any image modality for automatic segmentation., Methods: plastimatch mabs workflow consists of two main parts: (1) an offline phase, where optimal registration and voting parameters are tuned and (2) an online phase, where a new patient is labeled from scratch by using the same parameters as identified in the former phase. Several registration strategies, as well as different voting criteria can be selected. A flexible atlas selection scheme is also available. To prove the effectiveness of the proposed software across anatomical districts and image modalities, it was tested on two very different scenarios: head and neck (H&N) CT segmentation for radiotherapy application, and magnetic resonance image brain labeling for neuroscience investigation., Results: For the neurological study, minimum dice was equal to 0.76 (investigated structures: left and right caudate, putamen, thalamus, and hippocampus). For head and neck case, minimum dice was 0.42 for the most challenging structures (optic nerves and submandibular glands) and 0.62 for the other ones (mandible, brainstem, and parotid glands). Time required to obtain the labels was compatible with a real clinical workflow (35 and 120 min)., Conclusions: The proposed software fills a gap in the multiatlas based segmentation field, since all currently available tools (both for commercial and for research purposes) are restricted to a well specified application. Furthermore, it can be adopted as a platform for exploring MABS parameters and as a reference implementation for comparing against other segmentation algorithms.

Purpose: To improve the contouring of clinical target volume for the radiotherapy of neck Hodgkin/non-Hodgkin lymphoma by localizing the prechemotherapy gross target volume onto the simulation computed tomography using [18F]-fluorodeoxyglucose positron emission tomography/computed tomography., Material and Methods: The gross target volume delineated on prechemotherapy [18F]-fluorodeoxyglucose positron emission tomography/computed tomography images was warped onto simulation computed tomography using deformable image registration. Fifteen patients with neck Hodgkin/non-Hodgkin lymphoma were analyzed. Quality of image registration was measured by computing the Dice similarity coefficient on warped organs at risk. Five radiation oncologists visually scored the localization of automatic gross target volume, ranking it from 1 (wrong) to 5 (excellent). Deformable registration was compared to rigid registration by computing the overlap index between the automatic gross target volume and the planned clinical target volume and quantifying the V95 coverage., Results: The Dice similarity coefficient was 0.80 ± 0.07 (median ± quartiles). The physicians' survey had a median score equal to 4 (good). By comparing the rigid versus deformable registration, the overlap index increased from a factor of about 4 and the V95 (percentage of volume receiving the 95% of the prescribed dose) went from 0.84 ± 0.38 to 0.99 ± 0.10 (median ± quartiles)., Conclusion: This study demonstrates the impact of using deformable registration between prechemotherapy [18F]-fluorodeoxyglucose positron emission tomography/computed tomography and simulation computed tomography, in order to automatically localize the gross target volume for radiotherapy treatment of patients with Hodgkin/non-Hodgkin lymphoma., (© The Author(s) 2015.)

Purpose: To obtain a contrasted image of the tumor region during the setup for proton therapy in lung patients, by using proton radiography and x-ray computed tomography (CT) prior knowledge., Methods and Materials: Six lung cancer patients' CT scans were preprocessed by masking out the gross tumor volume (GTV), and digitally reconstructed radiographs along the planned beam's eye view (BEV) were generated, for a total of 27 projections. Proton radiographies (PR) were also computed for the same BEV through Monte Carlo simulations. The digitally reconstructed radiograph was subtracted from the corresponding proton image, resulting in a contrast-enhanced proton radiography (CEPR). Michelson contrast analysis was performed both on PR and CEPR. The tumor region was then automatically segmented on CEPR and compared to the ground truth (GT) provided by physicians in terms of Dice coefficient, accuracy, precision, sensitivity, and specificity., Results: Contrast on CEPR was, on average, 4 times better than on PR. For 10 lateral projections (±45° off of 90° or 270°), although it was not possible to distinguish the tumor region in the PR, CEPR offers excellent GTV visibility. The median ± quartile values of Dice, precision, and accuracy indexes were 0.86 ± 0.03, 0.86 ± 0.06, and 0.88 ± 0.02, respectively, thus confirming the reliability of the method in highlighting tumor boundaries. Sensitivity and specificity analysis demonstrated that there is no systematic over- or underestimation of the tumor region. Identification of the tumor boundaries using CEPR resulted in a more accurate and precise definition of GTV compared to that obtained from pretreatment CT., Conclusions: In most proton centers, the current clinical protocol is to align the patient using kV imaging with bony anatomy as a reference. We demonstrated that CEPR can significantly improve tumor visualization, allowing better patient set-up and permitting image guided proton therapy (IGPT)., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Purpose: Accurate delineation of organs at risk (OARs) is a precondition for intensity modulated radiation therapy. However, manual delineation of OARs is time consuming and prone to high interobserver variability. Because of image artifacts and low image contrast between different structures, however, the number of available approaches for autosegmentation of structures in the head-neck area is still rather low. In this project, a new approach for automated segmentation of head-neck CT images that combine the robustness of multiatlas-based segmentation with the flexibility of geodesic active contours and the prior knowledge provided by statistical appearance models is presented., Methods: The presented approach is using an atlas-based segmentation approach in combination with label fusion in order to initialize a segmentation pipeline that is based on using statistical appearance models and geodesic active contours. An anatomically correct approximation of the segmentation result provided by atlas-based segmentation acts as a starting point for an iterative refinement of this approximation. The final segmentation result is based on using model to image registration and geodesic active contours, which are mutually influencing each other., Results: 18 CT images in combination with manually segmented labels of parotid glands and brainstem were used in a leave-one-out cross validation scheme in order to evaluate the presented approach. For this purpose, 50 different statistical appearance models have been created and used for segmentation. Dice coefficient (DC), mean absolute distance and max. Hausdorff distance between the autosegmentation results and expert segmentations were calculated. An average Dice coefficient of DC = 0.81 (right parotid gland), DC = 0.84 (left parotid gland), and DC = 0.86 (brainstem) could be achieved., Conclusions: The presented framework provides accurate segmentation results for three important structures in the head neck area. Compared to a segmentation approach based on using multiple atlases in combination with label fusion, the proposed hybrid approach provided more accurate results within a clinically acceptable amount of time.

The aim of this work is to provide insights into multiple metrics clinical validation of deformable image registration and contour propagation methods in 4D lung radiotherapy planning. The following indices were analyzed and compared: Volume Difference (VD), Dice Similarity Coefficient (DSC), Positive Predictive Value (PPV) and Surface Distances (SD). The analysis was performed on three patient datasets, using as reference a ground-truth volume generated by means of Simultaneous Truth And Performance Level Estimation (STAPLE) algorithm from the outlines of five experts. Significant discrepancies in the quality assessment provided by the different metrics in all the examined cases were found. Metrics sensitivity was more evident in presence of image artifacts and particularly for tubular anatomical structures, such as esophagus or spinal cord. Volume Differences did not account for position and DSC exhibited criticalities due to its intrinsic symmetry (i.e. over- and under-estimation of the reference contours cannot be discriminated) and dependency on the total volume of the structure. PPV analysis showed more robust performance, as each voxel concurs to the classification of the propagation, but was not able to detect inclusion of propagated and ground-truth volumes. Mesh distances could interpret the actual shape of the structures, but might report higher mismatches in case of large local differences in the contour surfaces. According to our study, the combination of VD and SD for the validation of contour propagation algorithms in 4D could provide the necessary failure detection accuracy.

Purpose: To characterize the modulation transfer function (MTF) of proton/carbon radiography using Monte Carlo simulations. To assess the spatial resolution of proton/carbon radiographic imaging., Methods: A phantom was specifically modeled with inserts composed of two materials with three different densities of bone and lung. The basic geometry of the phantom consists of cube-shaped inserts placed in water. The thickness of the water, the thickness of the cubes, the depth of the cubes in the water, and the particle beam energy have all been varied and studied. There were two phantom thicknesses considered 20 and 28 cm. This represents an average patient thickness and a thicker sized patient. Radiographs were produced for proton beams at 230 and 330 MeV and for a carbon ion beam at 400 MeV per nucleon. The contrast-to-noise ratio (CNR) was evaluated at the interface of two materials on the radiographs, i.e., lung-water and bone-water. The variation in CNR at interface between lung-water and bone-water were study, where a sigmoidal fit was performed between the lower and the higher CNR values. The full width half-maximum (FWHM) value was then obtained from the sigmoidal fit. Ultimately, spatial resolution was defined by the 10% point of the modulation-transfer-function (MTF10%), in units of line-pairs per mm (lp/mm)., Results: For the 20 cm thick phantom, the FWHM values varied between 0.5 and 0.7 mm at the lung-water and bone-water interfaces, for the proton beam energies of 230 and 330 MeV and the 400 MeV/n carbon beam. For the 28 cm thick phantom, the FWHM values varied between 0.5 and 1.2 mm at the lung-water and bone-water interface for the same inserts and beam energies. For the 20 cm phantom the MTF10% for lung-water interface is 2.3, 2.4, and 2.8 lp/mm, respectively, for 230, 330, and 400 MeV/n beams. For the same 20 cm thick phantom but for the bone-water interface the MTF10% yielded 1.9, 2.3, and 2.7 lp/mm, respectively, for 230, 330, and 400 MeV/n beams. In the case of the thicker 28 cm phantom, the authors observed that at the lung-water interface the MTF10% is 1.6, 1.9, and 2.6 lp/mm, respectively, for 230, 330, and 400 MeV/n beams. While for the bone-water interface the MTF10% was 1.4, 1.9, and 2.9 lp/mm, respectively, for 230, 330, and 400 MeV/n beams., Conclusions: Carbon radiography (400 MeV/n) yielded best spatial resolution, with MTF10% = 2.7 and 2.8 lp/mm, respectively, at the lung-water and bone-water interfaces. The spatial resolution of the 330 MeV proton beam was better than the 230 MeV proton, because higher incident proton energy suffer smaller deflections within the patient and thus yields better proton radiographic images. The authors also observed that submillimeter resolution can be obtained with both proton and carbon beams.

The purpose of this work was to evaluate the intrapatient tumor position reproducibility in a deep inspiration breath-hold (DIBH) technique based on two infrared optical tracking systems, ExacTrac and ELITETM, in stereotactic treatment of lung and liver lesions. After a feasibility study, the technique was applied to 15 patients. Each patient, provided with a real-time visual feedback of external optical marker displacements, underwent a full DIBH, a free-breathing (FB), and three consecutive DIBH CT-scans centered on the lesion to evaluate the tumor position reproducibility. The mean reproducibility of tumor position during repeated DIBH was 0.5 ± 0.3 mm in laterolateral (LL), 1.0 ± 0.9 mm in anteroposterior (AP), and 1.4 ± 0.9 mm in craniocaudal (CC) direction for lung lesions, and 1.0 ± 0.6 mm in LL, 1.1 ± 0.5 mm in AP, and 1.2 ± 0.4 mm in CC direction for liver lesions. Intra- and interbreath-hold reproducibility during treatment, as determined by optical markers displacements, was below 1 mm and 3 mm, respectively, in all directions for all patients. Optically-guided DIBH technique provides a simple noninvasive method to minimize breathing motion for collaborative patients. For each patient, it is important to ensure that the tumor position is reproducible with respect to the external markers configuration.

The aim of this study was to assess the ability of metal artifact reduction (MAR) algorithm in restoring the CT image quality while correcting the tissue density information for the accurate estimation of the absorbed dose. A phantom filled with titanium (low-Z metal) and Cerrobend (high-Z metal) inserts was used for this purpose. The MAR algorithm was applied to phantom's CT dataset. Static intensity-modulated radiation therapy (IMRT) plans, including five beam angles, were designed and optimized on the uncorrected images to deliver 10 Gy on the simulated target. Monte Carlo dose calculation was computed on uncorrected, corrected, and ground truth image datasets. It was firstly verified that MAR methodology was able to correct HU errors due to the metal presence. In the worst situation (high-Z phantom), the image difference, uncorrected ground truth and corrected ground truth, went from -4.4 ± 118.8 HU to 0.4 ± 10.8 HU, respectively. Secondly, it was observed that the impact of dose errors estimation depends on the atomic number of the metal: low-Z inserts do not produce significant dose inaccuracies, while high-Z implants substantially influence the computation of the absorbed dose. In this latter case, dose errors in the PTV region were up to 23.56% (9.72% mean value) when comparing the uncorrected vs. the ground truth dataset. After MAR correction, errors dropped to 0.11% (0.10% mean value). In conclusion, it was assessed that the new MAR algorithm is able to restore image quality without distorting mass density information, thus producing a more accurate dose estimation.

The aim of this work was to develop and validate a 3D female breast deformation model for computer assisted breast surgery. Magnetic resonance (MR) image data of a patient undergoing breast biopsy, were acquired using two different protocols with the patient in prone position: (i) uncompressed breast and (ii) compressed breast, with lateral single breast compression, realized with a movable slab. The acquired images were then segmented using a semi-automatic procedure and from the extracted volumes of interest tetrahedral meshes representing skin, fat and mammary glands were generated. Tissue deformation was ruled by a mass-spring model: first, an iterative approximation algorithm was implemented to estimate the spring's rest length and stiffness, accounting for gravity force; then the resulting parameters were used to deform the uncompressed breast model in order to reach the real compressed one (ground truth). Results showed that gravity force applied to the mesh was properly compensated by the internal elastic forces, leading to a distance between the deformed mesh and the reference data of 0.036±0.092 mm (median±inter quartile range). The point to mesh residual distance between the deformed mesh and the ground truth was 1.224±2.202 mm (median±inter quartile range). Further investigation on a larger patient dataset is required for a more robust confirmation of model accuracy in predicting breast deformations., (Copyright © 2012 IPEM. Published by Elsevier Ltd. All rights reserved.)

Purpose: The purpose of this work was to develop and validate an efficient and automatic strategy to generate online virtual computed tomography (CT) scans for adaptive radiation therapy (ART) in head-and-neck (HN) cancer treatment., Method: We retrospectively analyzed 20 patients, treated with intensity modulated radiation therapy (IMRT), for an HN malignancy. Different anatomical structures were considered: mandible, parotid glands, and nodal gross tumor volume (nGTV). We generated 28 virtualCT scans by means of nonrigid registration of simulation computed tomography (CTsim) and cone beam CT images (CBCTs), acquired for patient setup. We validated our approach by considering the real replanning CT (CTrepl) as ground truth. We computed the Dice coefficient (DSC), center of mass (COM) distance, and root mean square error (RMSE) between correspondent points located on the automatically segmented structures on CBCT and virtualCT., Results: Residual deformation between CTrepl and CBCT was below one voxel. Median DSC was around 0.8 for mandible and parotid glands, but only 0.55 for nGTV, because of the fairly homogeneous surrounding soft tissues and of its small volume. Median COM distance and RMSE were comparable with image resolution. No significant correlation between RMSE and initial or final deformation was found., Conclusion: The analysis provides evidence that deformable image registration may contribute significantly in reducing the need of full CT-based replanning in HN radiation therapy by supporting swift and objective decision-making in clinical practice. Further work is needed to strengthen algorithm potential in nGTV localization., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Purpose: To compare infrared (IR) optical vs. stereoscopic X-ray technologies for patient setup in image-guided stereotactic radiotherapy., Methods and Materials: Retrospective data analysis of 233 fractions in 127 patients treated with hypofractionated stereotactic radiotherapy was performed. Patient setup at the linear accelerator was carried out by means of combined IR optical localization and stereoscopic X-ray image fusion in 6 degrees of freedom (6D). Data were analyzed to evaluate the geometric and dosimetric discrepancy between the two patient setup strategies., Results: Differences between IR optical localization and 6D X-ray image fusion parameters were on average within the expected localization accuracy, as limited by CT image resolution (3 mm). A disagreement between the two systems below 1 mm in all directions was measured in patients treated for cranial tumors. In extracranial sites, larger discrepancies and higher variability were observed as a function of the initial patient alignment. The compensation of IR-detected rotational errors resulted in a significantly improved agreement with 6D X-ray image fusion. On the basis of the bony anatomy registrations, the measured differences were found not to be sensitive to patient breathing. The related dosimetric analysis showed that IR-based patient setup caused limited variations in three cases, with 7% maximum dose reduction in the clinical target volume and no dose increase in organs at risk., Conclusions: In conclusion, patient setup driven by IR external surrogates localization in 6D featured comparable accuracy with respect to procedures based on stereoscopic X-ray imaging., (Copyright © 2012 Elsevier Inc. All rights reserved.)

Background: The purpose of this study is to investigate intra-fraction setup variability in hypo-fractionated cranial and body radiotherapy; this is achieved by means of integrated infrared optical localization and stereoscopic kV X-ray imaging., Method and Materials: We analyzed data coming from 87 patients treated with hypo-fractionated radiotherapy at cranial and extra-cranial sites. Patient setup was realized through the ExacTrac X-ray 6D system (BrainLAB, Germany), consisting of 2 infrared TV cameras for external fiducial localization and X-ray imaging in double projection for image registration. Before irradiation, patients were pre-aligned relying on optical marker localization. Patient position was refined through the automatic matching of X-ray images to digitally reconstructed radiographs, providing 6 corrective parameters that were automatically applied using a robotic couch. Infrared patient localization and X-ray imaging were performed at the end of treatment, thus providing independent measures of intra-fraction motion., Results: According to optical measurements, the size of intra-fraction motion was (median ± quartile) 0.3 ± 0.3 mm, 0.6 ± 0.6 mm, 0.7 ± 0.6 mm for cranial, abdominal and lung patients, respectively. X-ray image registration estimated larger intra-fraction motion, equal to 0.9 ± 0.8 mm, 1.3 ± 1.2 mm, 1.8 ± 2.2 mm, correspondingly., Conclusion: Optical tracking highlighted negligible intra-fraction motion at both cranial and extra-cranial sites. The larger motion detected by X-ray image registration showed significant inter-patient variability, in contrast to infrared optical tracking measurement. Infrared localization is put forward as the optimal strategy to monitor intra-fraction motion, featuring robustness, flexibility and less invasivity with respect to X-ray based techniques.

Purpose: Artifacts affect 4D CT images due to breathing irregularities or incorrect breathing phase identification. The purpose of this study is the reduction of artifacts in sorted 4D CT images. The assumption is that the use of multiple respiratory related signals may reduce uncertainties and increase robustness in breathing phase identification., Methods: Multiple respiratory related signals were provided by infrared 3D localization of a configuration of markers placed on the thoracoabdominal surface. Multidimensional K-means clustering was used for retrospective 4D CT image sorting, which was based on multiple marker variables, in order to identify clusters representing different breathing phases. The proposed technique was tested on computational simulations, phantom experimental acquisitions, and clinical data coming from two patients. Computational simulations provided a controlled and noise-free condition for testing the clustering technique on regular and irregular breathing signals, including baseline drift, time variant amplitude, time variant frequency, and end-expiration plateau. Specific attention was given to cluster initialization. Phantom experiments involved two moving phantoms fitted with multiple markers. Phantoms underwent 4D CT acquisition while performing controlled rigid motion patterns and featuring end-expiration plateau. Breathing cycle period and plateau duration were controlled by means of weights leaned upon the phantom during repeated 4D CT scans. The implemented sorting technique was applied to clinical 4D CT scans acquired on two patients and results were compared to conventional sorting methods., Results: For computational simulations and phantom studies, the performance of the multidimensional clustering technique was evaluated by measuring the repeatability in identifying the breathing phase among adjacent couch positions and the uniformity in sampling the breathing cycle. When breathing irregularities were present, the clustering technique consistently improved breathing phase identification with respect to conventional sorting methods based on monodimensional signals. In patient studies, a qualitative comparison was performed between corresponding breathing phases of 4D CT images obtained by conventional sorting methods and by the described clustering technique. Artifact reduction was clearly observable on both data set especially in the lower lung region., Conclusions: The implemented multiple point method demonstrated the ability to reduce artifacts in 4D CT imaging. Further optimization and development are needed to make the most of the availability of multiple respiratory related variables and to extend the method to 4D CT-PET hybrid scan.

The purpose of this study is to assess the temporal and reconstruction accuracy of a surface imaging system, the GateCT under ideal conditions, and compare the device with a commonly used respiratory surrogate: the Varian RPM. A clinical CT scanner, run in cine mode, was used with two optical devices, GateCT and RPM, to detect respiratory motion. A radiation detector, GM-10, triggers the X-ray on/off to GateCT system, while the RPM is directly synchronized with the CT scanner through an electronic connection. Two phantoms were imaged: the first phantom translated on a rigid plate along the anterior-posterior (AP) direction, and was used to assess the temporal synchronization of each optical system with the CT scanner. The second phantom, consisting of five spheres translating 3 cm peak-to-peak in the superior-inferior direction, was used to assess the quality of rebinned images created by GateCT and RPM. Calibration assessment showed a nearly perfect synchronization with the scanner for both the RPM and GateCT systems, thus demonstrating the good performance of the radiation detector. Results for the volume rebinning test showed discrepancies in volumes for the 3D reconstruction (compared to ground truth) of up to 36% for GateCT and up to 40% for RPM. No statistical difference was proven between the two systems in volume sorting. Errors are mainly due to phase detection inaccuracies and to the large motion of the phantom. This feasibility study assessed the consistency of two optical systems in synchronizing the respiratory signal with the image acquisition. A new patient protocol based on both RPM and GateCT will be soon started.