Outcome and Prognostic Factors in Adenosquamous Carcinoma of the Head and Neck - a Multicenter Rare Cancer Network Study

Background: Adenosquamous carcinoma (AC) of the head and neck is a distinct entity first described in 1968. Its natural history is more aggressive than squamous-cell carcinoma. The aim of this study was to assess the clinical profile, patterns of failure, and prognostic factors in patients with AC of the head and neck treated by radiation therapy (RT) with or without chemotherapy (CT).Materials and Methods: Data from 19 patients with stage I (n = 3), II (n = 1), III (n = 4), or IVa (n = 11) AC, treated between 1989 and 2009, were collected in a retrospective multicenter Rare Cancer Network study. Median age was 60 years (range, 48−73). Fifteen patients were male, and 4 female. Risk factors, including perineural invasion, lymphangitis, vascular invasion, positive margins were present in the majority (83%) of the patients. Tumour sites included oral cavity in 4, oropharynx in 4, hypopharynx in 2, larynx in 2, salivary glands in 2, nasal vestibule in 2, maxillary sinus in 2, and nasopharynx in 1 patient. Surgery (S) was performed in all but 5 patients. S alone was performed in only 1 patient, and definitive RT alone in 3 patients. Fifteen patients received combined modality treatment (S+RT in 11, RT+CT in 2, and all of the three modalities in 2 patients). Median RT dose to the primary and to the nodes was 66 Gy (range, 50−72) and 53 Gy (range, 44−66), respectively (1.8−2.0 Gy/fr., 5 fr./week). In 4 patients, the planning treatment volume included the primary tumour site only. Eight patients were treated with 2D RT, 7 with 3D conformal RT, and 2 with intensity-modulated RT.Results: After a median follow-up period of 39 months (range, 9−62), 9 patients developed distant metastases (lung, bone, mediastinum, and liver), 7 presented nodal recurrences, and only 4 had a local relapse at the primary site (all in-field recurrences). At last follow-up, 7 patients were alive without disease, 1 alive with disease, 9 died from progressive disease, and 2 died from intercurrent disease. The 3-year and median overall survival, disease-free survival (DFS), and locoregional control rates were 55% (95% confidence interval [CI]: 32−78%) and 39 months, 34% (95% CI: 12−56%) and 22 months, and 50% (95% CI: 22−78%) and 33 months, respectively. In multivariate analysis (Cox model), DFS was negatively influenced by the presence of extracapsular extension (p = 0.01) and advanced stage (IV versus I−III, p = 0.002).Conclusions: Overall prognosis of locoregionally advanced AC remains poor, and distant metastases and nodal relapse occur in almost half of the cases. However, local control is relatively better, and early stage AC patients had prolonged DFS when treated with combined-modality treatment.

Anàlisi de l'aplicació GeoMedia 3D

Treball de final de carrera sobre els sistemes d'informació geogràfica on s'analitza l'eina GeoMedia 3D i la seva aplicació als SIG 3D. L'anàlisi està fet desenvolupant un model 3D de la ciutat de Tarragona.

Anàlisi de l'aplicació GeoMedia 3D

Aquest treball es centra en la utilització i aplicació de les tecnologies dels sistemes d'informació geogràfica. En concret es realitza una pràctica sobre la ciutat de Lleida per poder veure certs edificis que compleixen una sèrie de premisses. En el projecte es fa una anàlisi de l'aplicació GeoMedia 3D.

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

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

Creación de una demo 3D mediante gráficos procedurales

Aquest projecte consisteix en el desenvolupament d’una demo 3D utilitzant exclusivament gràfics procedurals per tal d’avaluar la seva viabilitat en aplicacions més complexes com els videojocs. En aquesta aplicació es genera un terreny aleatori explorable amb vegetació i textures creades proceduralment.

Dose reconstruction in the context of tomotherapy

In vivo dosimetry is a way to verify the radiation dose delivered to the patient in measuring the dose generally during the first fraction of the treatment. It is the only dose delivery control based on a measurement performed during the treatment. In today's radiotherapy practice, the dose delivered to the patient is planned using 3D dose calculation algorithms and volumetric images representing the patient. Due to the high accuracy and precision necessary in radiation treatments, national and international organisations like ICRU and AAPM recommend the use of in vivo dosimetry. It is also mandatory in some countries like France. Various in vivo dosimetry methods have been developed during the past years. These methods are point-, line-, plane- or 3D dose controls. A 3D in vivo dosimetry provides the most information about the dose delivered to the patient, with respect to ID and 2D methods. However, to our knowledge, it is generally not routinely applied to patient treatments yet. The aim of this PhD thesis was to determine whether it is possible to reconstruct the 3D delivered dose using transmitted beam measurements in the context of narrow beams. An iterative dose reconstruction method has been described and implemented. The iterative algorithm includes a simple 3D dose calculation algorithm based on the convolution/superposition principle. The methodology was applied to narrow beams produced by a conventional 6 MV linac. The transmitted dose was measured using an array of ion chambers, as to simulate the linear nature of a tomotherapy detector. We showed that the iterative algorithm converges quickly and reconstructs the dose within a good agreement (at least 3% / 3 mm locally), which is inside the 5% recommended by the ICRU. Moreover it was demonstrated on phantom measurements that the proposed method allows us detecting some set-up errors and interfraction geometry modifications. We also have discussed the limitations of the 3D dose reconstruction for dose delivery error detection. Afterwards, stability tests of the tomotherapy MVCT built-in onboard detector was performed in order to evaluate if such a detector is suitable for 3D in-vivo dosimetry. The detector showed stability on short and long terms comparable to other imaging devices as the EPIDs, also used for in vivo dosimetry. Subsequently, a methodology for the dose reconstruction using the tomotherapy MVCT detector is proposed in the context of static irradiations. This manuscript is composed of two articles and a script providing further information related to this work. In the latter, the first chapter introduces the state-of-the-art of in vivo dosimetry and adaptive radiotherapy, and explains why we are interested in performing 3D dose reconstructions. In chapter 2 a dose calculation algorithm implemented for this work is reviewed with a detailed description of the physical parameters needed for calculating 3D absorbed dose distributions. The tomotherapy MVCT detector used for transit measurements and its characteristics are described in chapter 3. Chapter 4 contains a first article entitled '3D dose reconstruction for narrow beams using ion chamber array measurements', which describes the dose reconstruction method and presents tests of the methodology on phantoms irradiated with 6 MV narrow photon beams. Chapter 5 contains a second article 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations. A dose reconstruction process specific to the use of the tomotherapy MVCT detector is presented in chapter 6. A discussion and perspectives of the PhD thesis are presented in chapter 7, followed by a conclusion in chapter 8. The tomotherapy treatment device is described in appendix 1 and an overview of 3D conformai- and intensity modulated radiotherapy is presented in appendix 2. - La dosimétrie in vivo est une technique utilisée pour vérifier la dose délivrée au patient en faisant une mesure, généralement pendant la première séance du traitement. Il s'agit de la seule technique de contrôle de la dose délivrée basée sur une mesure réalisée durant l'irradiation du patient. La dose au patient est calculée au moyen d'algorithmes 3D utilisant des images volumétriques du patient. En raison de la haute précision nécessaire lors des traitements de radiothérapie, des organismes nationaux et internationaux tels que l'ICRU et l'AAPM recommandent l'utilisation de la dosimétrie in vivo, qui est devenue obligatoire dans certains pays dont la France. Diverses méthodes de dosimétrie in vivo existent. Elles peuvent être classées en dosimétrie ponctuelle, planaire ou tridimensionnelle. La dosimétrie 3D est celle qui fournit le plus d'information sur la dose délivrée. Cependant, à notre connaissance, elle n'est généralement pas appliquée dans la routine clinique. Le but de cette recherche était de déterminer s'il est possible de reconstruire la dose 3D délivrée en se basant sur des mesures de la dose transmise, dans le contexte des faisceaux étroits. Une méthode itérative de reconstruction de la dose a été décrite et implémentée. L'algorithme itératif contient un algorithme simple basé sur le principe de convolution/superposition pour le calcul de la dose. La dose transmise a été mesurée à l'aide d'une série de chambres à ionisations alignées afin de simuler la nature linéaire du détecteur de la tomothérapie. Nous avons montré que l'algorithme itératif converge rapidement et qu'il permet de reconstruire la dose délivrée avec une bonne précision (au moins 3 % localement / 3 mm). De plus, nous avons démontré que cette méthode permet de détecter certaines erreurs de positionnement du patient, ainsi que des modifications géométriques qui peuvent subvenir entre les séances de traitement. Nous avons discuté les limites de cette méthode pour la détection de certaines erreurs d'irradiation. Par la suite, des tests de stabilité du détecteur MVCT intégré à la tomothérapie ont été effectués, dans le but de déterminer si ce dernier peut être utilisé pour la dosimétrie in vivo. Ce détecteur a démontré une stabilité à court et à long terme comparable à d'autres détecteurs tels que les EPIDs également utilisés pour l'imagerie et la dosimétrie in vivo. Pour finir, une adaptation de la méthode de reconstruction de la dose a été proposée afin de pouvoir l'implémenter sur une installation de tomothérapie. Ce manuscrit est composé de deux articles et d'un script contenant des informations supplémentaires sur ce travail. Dans ce dernier, le premier chapitre introduit l'état de l'art de la dosimétrie in vivo et de la radiothérapie adaptative, et explique pourquoi nous nous intéressons à la reconstruction 3D de la dose délivrée. Dans le chapitre 2, l'algorithme 3D de calcul de dose implémenté pour ce travail est décrit, ainsi que les paramètres physiques principaux nécessaires pour le calcul de dose. Les caractéristiques du détecteur MVCT de la tomothérapie utilisé pour les mesures de transit sont décrites dans le chapitre 3. Le chapitre 4 contient un premier article intitulé '3D dose reconstruction for narrow beams using ion chamber array measurements', qui décrit la méthode de reconstruction et présente des tests de la méthodologie sur des fantômes irradiés avec des faisceaux étroits. Le chapitre 5 contient un second article intitulé 'Stability of the Helical TomoTherapy HiArt II detector for treatment beam irradiations'. Un procédé de reconstruction de la dose spécifique pour l'utilisation du détecteur MVCT de la tomothérapie est présenté au chapitre 6. Une discussion et les perspectives de la thèse de doctorat sont présentées au chapitre 7, suivies par une conclusion au chapitre 8. Le concept de la tomothérapie est exposé dans l'annexe 1. Pour finir, la radiothérapie «informationnelle 3D et la radiothérapie par modulation d'intensité sont présentées dans l'annexe 2.

Intracellular nanomanipulation by a photonic-force microscope with real-time acquisition of a 3D stiffness matrix.

A traditional photonic-force microscope (PFM) results in huge sets of data, which requires tedious numerical analysis. In this paper, we propose instead an analog signal processor to attain real-time capabilities while retaining the richness of the traditional PFM data. Our system is devoted to intracellular measurements and is fully interactive through the use of a haptic joystick. Using our specialized analog hardware along with a dedicated algorithm, we can extract the full 3D stiffness matrix of the optical trap in real time, including the off-diagonal cross-terms. Our system is also capable of simultaneously recording data for subsequent offline analysis. This allows us to check that a good correlation exists between the classical analysis of stiffness and our real-time measurements. We monitor the PFM beads using an optical microscope. The force-feedback mechanism of the haptic joystick helps us in interactively guiding the bead inside living cells and collecting information from its (possibly anisotropic) environment. The instantaneous stiffness measurements are also displayed in real time on a graphical user interface. The whole system has been built and is operational; here we present early results that confirm the consistency of the real-time measurements with offline computations.

Patient-specific computational hemodynamics of intracranial aneurysms from 3D rotational angiography and CT angiography: an in vivo reproducibility study

Patient-specific simulations of the hemodynamics in intracranial aneurysms can be constructed by using image-based vascular models and CFD techniques. This work evaluates the impact of the choice of imaging technique on these simulations

Calibration and Validation of SWAT for the Upper Maquoketa River Watershed, June 2005

A validation study has been performed using the Soil and Water Assessment Tool (SWAT) model with data collected for the Upper Maquoketa River Watershed (UMRW), which drains over 16,000 ha in northeast Iowa. This validation assessment builds on a previous study with nested modeling for the UMRW that required both the Agricultural Policy EXtender (APEX) model and SWAT. In the nested modeling approach, edge-offield flows and pollutant load estimates were generated for manure application fields with APEX and were then subsequently routed to the watershed outlet in SWAT, along with flows and pollutant loadings estimated for the rest of the watershed routed to the watershed outlet. In the current study, the entire UMRW cropland area was simulated in SWAT, which required translating the APEX subareas into SWAT hydrologic response units (HRUs). Calibration and validation of the SWAT output was performed by comparing predicted flow and NO3-N loadings with corresponding in-stream measurements at the watershed outlet from 1999 to 2001. Annual stream flows measured at the watershed outlet were greatly under-predicted when precipitation data collected within the watershed during the 1999-2001 period were used to drive SWAT. Selection of alternative climate data resulted in greatly improved average annual stream predictions, and also relatively strong r2 values of 0.73 and 0.72 for the predicted average monthly flows and NO3-N loads, respectively. The impact of alternative precipitation data shows that as average annual precipitation increases 19%, the relative change in average annual streamflow is about 55%. In summary, the results of this study show that SWAT can replicate measured trends for this watershed and that climate inputs are very important for validating SWAT and other water quality models.

Temporal diffeomorphic free-form deformation: application to motion and strain estimation from 3D echocardiography

This paper presents a new registration algorithm, called Temporal Di eomorphic Free Form Deformation (TDFFD), and its application to motion and strain quanti cation from a sequence of 3D ultrasound (US) images. The originality of our approach resides in enforcing time consistency by representing the 4D velocity eld as the sum of continuous spatiotemporal B-Spline kernels. The spatiotemporal displacement eld is then recovered through forward Eulerian integration of the non-stationary velocity eld. The strain tensor iscomputed locally using the spatial derivatives of the reconstructed displacement eld. The energy functional considered in this paper weighs two terms: the image similarity and a regularization term. The image similarity metric is the sum of squared di erences between the intensities of each frame and a reference one. Any frame in the sequence can be chosen as reference. The regularization term is based on theincompressibility of myocardial tissue. TDFFD was compared to pairwise 3D FFD and 3D+t FFD, bothon displacement and velocity elds, on a set of synthetic 3D US images with di erent noise levels. TDFFDshowed increased robustness to noise compared to these two state-of-the-art algorithms. TDFFD also proved to be more resistant to a reduced temporal resolution when decimating this synthetic sequence. Finally, this synthetic dataset was used to determine optimal settings of the TDFFD algorithm. Subsequently, TDFFDwas applied to a database of cardiac 3D US images of the left ventricle acquired from 9 healthy volunteers and 13 patients treated by Cardiac Resynchronization Therapy (CRT). On healthy cases, uniform strain patterns were observed over all myocardial segments, as physiologically expected. On all CRT patients, theimprovement in synchrony of regional longitudinal strain correlated with CRT clinical outcome as quanti ed by the reduction of end-systolic left ventricular volume at follow-up (6 and 12 months), showing the potential of the proposed algorithm for the assessment of CRT.

Dynamic estimation of three-dimensional cerebrovascular deformation from rotational angiography

Purpose: The objective of this study is to investigate the feasibility of detecting and quantifying 3D cerebrovascular wall motion from a single 3D rotational x-ray angiography (3DRA) acquisition within a clinically acceptable time and computing from the estimated motion field for the further biomechanical modeling of the cerebrovascular wall. Methods: The whole motion cycle of the cerebral vasculature is modeled using a 4D B-spline transformation, which is estimated from a 4D to 2D + t image registration framework. The registration is performed by optimizing a single similarity metric between the entire 2D + t measured projection sequence and the corresponding forward projections of the deformed volume at their exact time instants. The joint use of two acceleration strategies, together with their implementation on graphics processing units, is also proposed so as to reach computation times close to clinical requirements. For further characterizing vessel wall properties, an approximation of the wall thickness changes is obtained through a strain calculation. Results: Evaluation on in silico and in vitro pulsating phantom aneurysms demonstrated an accurate estimation of wall motion curves. In general, the error was below 10% of the maximum pulsation, even in the situation when substantial inhomogeneous intensity pattern was present. Experiments on in vivo data provided realistic aneurysm and vessel wall motion estimates, whereas in regions where motion was neither visible nor anatomically possible, no motion was detected. The use of the acceleration strategies enabled completing the estimation process for one entire cycle in 5-10 min without degrading the overall performance. The strain map extracted from our motion estimation provided a realistic deformation measure of the vessel wall. Conclusions: The authors' technique has demonstrated that it can provide accurate and robust 4D estimates of cerebrovascular wall motion within a clinically acceptable time, although it has to be applied to a larger patient population prior to possible wide application to routine endovascular procedures. In particular, for the first time, this feasibility study has shown that in vivo cerebrovascular motion can be obtained intraprocedurally from a 3DRA acquisition. Results have also shown the potential of performing strain analysis using this imaging modality, thus making possible for the future modeling of biomechanical properties of the vascular wall.

Magnetization transfer-based 3D visualization of foot peripheral nerves.

PURPOSE: To investigate magnetization transfer (MT) effects as a new source of contrast for imaging and tracking of peripheral foot nerves. MATERIALS AND METHODS: Two sets of 3D spoiled gradient-echo images acquired with and without a saturation pulse were used to generate MT ratio (MTR) maps of 260 μm in-plane resolution for eight volunteers at 3T. Scan parameters were adjusted to minimize signal loss due to T2 dephasing, and a dedicated coil was used to improve the inherently low signal-to-noise ratio of small voxels. Resulting MTR values in foot nerves were compared with those in surrounding muscle tissue. RESULTS: Average MTR values for muscle (45.5 ± 1.4%) and nerve (21.4 ± 3.1%) were significantly different (P < 0.0001). In general, the difference in MTR values was sufficiently large to allow for intensity-based segmentation and tracking of foot nerves in individual subjects. This procedure was termed MT-based 3D visualization. CONCLUSION: The MTR serves as a new source of contrast for imaging of peripheral foot nerves and provides a means for high spatial resolution tracking of these structures. The proposed methodology is directly applicable on standard clinical MR scanners and could be applied to systemic pathologies, such as diabetes.

Is spinal stenosis assessment dependent on slice orientation? A magnetic resonance imaging study.

INTRODUCTION: Lumbar spinal stenosis (LSS) treatment is based primarily on the clinical criteria providing that imaging confirms radiological stenosis. The radiological measurement more commonly used is the dural sac cross-sectional area (DSCA). It has been recently shown that grading stenosis based on the morphology of the dural sac as seen on axial T2 MRI images, better reflects severity of stenosis than DSCA and is of prognostic value. This radiological prospective study investigates the variability of surface measurements and morphological grading of stenosis for varying degrees of angulation of the T2 axial images relative to the disc space as observed in clinical practice. MATERIALS AND METHODS: Lumbar spine TSE T2 three-dimensional (3D) MRI sequences were obtained from 32 consecutive patients presenting with either suspected spinal stenosis or low back pain. Axial reconstructions using the OsiriX software at 0°, 10°, 20° and 30° relative to the disc space orientation were obtained for a total of 97 levels. For each level, DSCA was digitally measured and stenosis was graded according to the 4-point (A-D) morphological grading by two observers. RESULTS: A good interobserver agreement was found in grade evaluation of stenosis (k = 0.71). DSCA varied significantly as the slice orientation increased from 0° to +10°, +20° and +30° at each level examined (P < 0.0001) (-15 to +32% at 10°, -24 to +143% at 20° and -29 to +231% at 30° of slice orientation). Stenosis definition based on the surface measurements changed in 39 out of the 97 levels studied, whereas the morphology grade was modified only in two levels (P < 0.01). DISCUSSION: The need to obtain continuous slices using the classical 2D MRI acquisition technique entails often at least a 10° slice inclination relative to one of the studied discs. Even at this low angulation, we found a significantly statistical difference between surface changes and morphological grading change. In clinical practice, given the above findings, it might therefore not be necessary to align the axial cuts to each individual disc level which could be more time-consuming than obtaining a single series of axial cuts perpendicular to the middle of the lumbar spine or to the most stenotic level. In conclusion, morphological grading seems to offer an alternative means of assessing severity of spinal stenosis that is little affected by image acquisition technique.