Development of Techniques for the Automatic Extraction and Grade Detection of Glioma Tumors from Conventional Brain Magnetic Resonant Images

Cerebral glioma is the most prevalent primary brain tumor, which are classified broadly into low and high grades according to the degree of malignancy. High grade gliomas are highly malignant which possess a poor prognosis, and the patients survive less than eighteen months after diagnosis. Low grade gliomas are slow growing, least malignant and has better response to therapy. To date, histological grading is used as the standard technique for diagnosis, treatment planning and survival prediction. The main objective of this thesis is to propose novel methods for automatic extraction of low and high grade glioma and other brain tissues, grade detection techniques for glioma using conventional magnetic resonance imaging (MRI) modalities and 3D modelling of glioma from segmented tumor slices in order to assess the growth rate of tumors. Two new methods are developed for extracting tumor regions, of which the second method, named as Adaptive Gray level Algebraic set Segmentation Algorithm (AGASA) can also extract white matter and grey matter from T1 FLAIR an T2 weighted images. The methods were validated with manual Ground truth images, which showed promising results. The developed methods were compared with widely used Fuzzy c-means clustering technique and the robustness of the algorithm with respect to noise is also checked for different noise levels. Image texture can provide significant information on the (ab)normality of tissue, and this thesis expands this idea to tumour texture grading and detection. Based on the thresholds of discriminant first order and gray level cooccurrence matrix based second order statistical features three feature sets were formulated and a decision system was developed for grade detection of glioma from conventional T2 weighted MRI modality.The quantitative performance analysis using ROC curve showed 99.03% accuracy for distinguishing between advanced (aggressive) and early stage (non-aggressive) malignant glioma. The developed brain texture analysis techniques can improve the physician’s ability to detect and analyse pathologies leading to a more reliable diagnosis and treatment of disease. The segmented tumors were also used for volumetric modelling of tumors which can provide an idea of the growth rate of tumor; this can be used for assessing response to therapy and patient prognosis.

Quantile based stop-loss transform and its applications

Partial moments are extensively used in actuarial science for the analysis of risks. Since the first order partial moments provide the expected loss in a stop-loss treaty with infinite cover as a function of priority, it is referred as the stop-loss transform. In the present work, we discuss distributional and geometric properties of the first and second order partial moments defined in terms of quantile function. Relationships of the scaled stop-loss transform curve with the Lorenz, Gini, Bonferroni and Leinkuhler curves are developed

Automatic Detection and Classification of Glioma Tumors using Statistical Features

The characterization and grading of glioma tumors, via image derived features, for diagnosis, prognosis, and treatment response has been an active research area in medical image computing. This paper presents a novel method for automatic detection and classification of glioma from conventional T2 weighted MR images. Automatic detection of the tumor was established using newly developed method called Adaptive Gray level Algebraic set Segmentation Algorithm (AGASA).Statistical Features were extracted from the detected tumor texture using first order statistics and gray level co-occurrence matrix (GLCM) based second order statistical methods. Statistical significance of the features was determined by t-test and its corresponding p-value. A decision system was developed for the grade detection of glioma using these selected features and its p-value. The detection performance of the decision system was validated using the receiver operating characteristic (ROC) curve. The diagnosis and grading of glioma using this non-invasive method can contribute promising results in medical image computing

Preparation and Study of Optical Properties of CdS, CdS-TiO2 and CdS-Au Nanocomposites for Photonic Applications

Nanophotonics can be regarded as a fusion of nanotechnology and photonics and it is an emerging field providing researchers opportunities in fundamental science and new technologies. In recent times many new methodsand techniques have been developed to prepare materials at nanoscale dimensions. Most of these materials exhibit unique and interesting optical properties and behavior. Many of these have been found to be very useful to develop new devices and systems such as tracers in biological systems, optical limiters, light emitters and energy harvesters. This thesis presents a summary of the work done by the author in the field by choosing a few semiconductor systems to prepare nanomaterials and nanocomposites. Results of the study of linear and nonlinear optical properties of materials thus synthesized are also presented in the various chapters of this thesis. CdS is the material chosen here and the methods and the studies of the detailed investigation are presented in this thesis related to the optical properties of CdS nanoparticles and its composites. Preparation and characterization methods and experimental techniques adopted for the investigations were illustrated in chapter 2 of this thesis. Chapter 3 discusses the preparation of CdS, TiO2 and Au nanoparticles. We observed that the fluorescence behaviour of the CdS nanoparticles, prepared by precipitation technique, depends on excitation wavelength. It was found that the peak emission wavelength can be shifted by as much as 147nm by varyingthe excitation wavelengths and the reason for this phenomenon is the selective excitation of the surface states in the nanoparticles. This provided certain amount of tunability for the emission which results from surface states.TiO2 nanoparticle colloids were prepared by hydrothermal method. The optical absorption study showed a blue shift of absorption edge, indicating quantum confinement effect. The large spectral range investigated allows observing simultaneously direct and indirect band gap optical recombination. The emission studies carried out show four peaks, which are found to be generated from excitonic as well as surface state transitions. It was found that the emission wavelengths of these colloidal nanoparticles and annealed nanoparticles showed two category of surface state emission in addition to the excitonic emission. Au nanoparticles prepared by Turkevich method showed nanoparticles of size below 5nm using plasmonic absorption calculation. It was also found that there was almost no variation in size as the concentration of precursor was changed from 0.2mM to 0.4mM.We have observed SHG from CdS nanostructured thin film prepared onglass substrate by chemical bath deposition technique. The results point out that studied sample has in-plane isotropy. The relative values of tensor components of the second-order susceptibility were determined to be 1, zzz 0.14, xxz and 0.07. zxx These values suggest that the nanocrystals are oriented along the normal direction. However, the origin of such orientation remains unknown at present. Thus CdS is a promising nonlinear optical material for photonic applications, particularly for integrated photonic devices. CdS Au nanocomposite particles were prepared by mixing CdS nanoparticles with Au colloidal nanoparticles. Optical absorption study of these nanoparticles in PVA solution suggests that absorption tail was red shifted compared to CdS nanoparticles. TEM and EDS analysis suggested that the amount of Au nanoparticles present on CdS nanoparticles is very small. Fluorescence emission is unaffected indicating the presence of low level of Au nanoparticles. CdS:Au PVA and CdS PVA nanocomposite films were fabricated and optically characterized. The results showed a red-shift for CdS:Au PVA film for absorption tail compared to CdS PVA film. Nonlinear optical analysis showed a huge nonlinear optical absorption for CdS:Au PVA nanocomposite and CdS:PVA films. Also an enhancement in nonlinear optical absorption is found for CdS:Au PVA thin film compared to the CdS PVA thin film. This enhancement is due to the combined effect of plasmonic as well as excitonic contribution at high input intensity. Samples of CdS doped with TiO2 were also prepared and the linear optical absorption spectra of these nanocompositeparticles clearly indicated the influence of TiO2 nanoparticles. TEM and EDS studies have confirmed the presence of TiO2 on CdS nanoparticles. Fluorescence studies showed that there is an increase in emission peak around 532nm for CdS nanoparticles. Nonlinear optical analysis of CdS:TiO2 PVA nanocomposite films indicated a large nonlinear optical absorption compared to that of CdS:PVA nanocomposite film. The values of nonlinear optical absorption suggests that these nanocomposite particles can be employed for optical limiting applications. CdSe-CdS and CdSe-ZnS core-shell QDs with varying shell size were characterized using UV–VIS spectroscopy. Optical absorption and TEM analysis of these QDs suggested a particle size around 5 nm. It is clearly shown that the surface coating influences the optical properties of QDs in terms of their size. Fluorescence studies reveal the presence of trap states in CdSe-CdS and CdSe- ZnS QDs. Trap states showed an increase as a shell for CdS is introduced and increasing the shell size of CdS beyond a certain value leads to a decrease in the trap state emission. There is no sizeable nonlinear optical absorption observed. In the case of CdSe- ZnS QDs, the trap state emission gets enhanced with the increase in ZnS shell thickness. The enhancement of emission from trap states transition due to the increase in thickness of ZnS shell gives a clear indication of distortion occurring in the spherical symmetry of CdSe quantum dots. Consequently the nonlinear optical absorption of CdSe-ZnS QDs gets increased and the optical limiting threshold is decreased as the shell thickness is increased in respect of CdSe QDs. In comparison with CdSe-CdS QDs, CdSe-ZnS QDs possess much better optical properties and thereby CdSe-ZnS is a strong candidate for nonlinear as well as linear optical applications.

Orthogonal polynomials and recurrence equations, operator equations and factorization

This article surveys the classical orthogonal polynomial systems of the Hahn class, which are solutions of second-order differential, difference or q-difference equations. Orthogonal families satisfy three-term recurrence equations. Example applications of an algorithm to determine whether a three-term recurrence equation has solutions in the Hahn class - implemented in the computer algebra system Maple - are given. Modifications of these families, in particular associated orthogonal systems, satisfy fourth-order operator equations. A factorization of these equations leads to a solution basis.

Some New Classes of Orthogonal Polynomials and Special Functions

In dieser Dissertation präsentieren wir zunächst eine Verallgemeinerung der üblichen Sturm-Liouville-Probleme mit symmetrischen Lösungen und erklären eine umfassendere Klasse. Dann führen wir einige neue Klassen orthogonaler Polynome und spezieller Funktionen ein, welche sich aus dieser symmetrischen Verallgemeinerung ableiten lassen. Als eine spezielle Konsequenz dieser Verallgemeinerung führen wir ein Polynomsystem mit vier freien Parametern ein und zeigen, dass in diesem System fast alle klassischen symmetrischen orthogonalen Polynome wie die Legendrepolynome, die Chebyshevpolynome erster und zweiter Art, die Gegenbauerpolynome, die verallgemeinerten Gegenbauerpolynome, die Hermitepolynome, die verallgemeinerten Hermitepolynome und zwei weitere neue endliche Systeme orthogonaler Polynome enthalten sind. All diese Polynome können direkt durch das neu eingeführte System ausgedrückt werden. Ferner bestimmen wir alle Standardeigenschaften des neuen Systems, insbesondere eine explizite Darstellung, eine Differentialgleichung zweiter Ordnung, eine generische Orthogonalitätsbeziehung sowie eine generische Dreitermrekursion. Außerdem benutzen wir diese Erweiterung, um die assoziierten Legendrefunktionen, welche viele Anwendungen in Physik und Ingenieurwissenschaften haben, zu verallgemeinern, und wir zeigen, dass diese Verallgemeinerung Orthogonalitätseigenschaft und -intervall erhält. In einem weiteren Kapitel der Dissertation studieren wir detailliert die Standardeigenschaften endlicher orthogonaler Polynomsysteme, welche sich aus der üblichen Sturm-Liouville-Theorie ergeben und wir zeigen, dass sie orthogonal bezüglich der Fisherschen F-Verteilung, der inversen Gammaverteilung und der verallgemeinerten t-Verteilung sind. Im nächsten Abschnitt der Dissertation betrachten wir eine vierparametrige Verallgemeinerung der Studentschen t-Verteilung. Wir zeigen, dass diese Verteilung gegen die Normalverteilung konvergiert, wenn die Anzahl der Stichprobe gegen Unendlich strebt. Eine ähnliche Verallgemeinerung der Fisherschen F-Verteilung konvergiert gegen die chi-Quadrat-Verteilung. Ferner führen wir im letzten Abschnitt der Dissertation einige neue Folgen spezieller Funktionen ein, welche Anwendungen bei der Lösung in Kugelkoordinaten der klassischen Potentialgleichung, der Wärmeleitungsgleichung und der Wellengleichung haben. Schließlich erklären wir zwei neue Klassen rationaler orthogonaler hypergeometrischer Funktionen, und wir zeigen unter Benutzung der Fouriertransformation und der Parsevalschen Gleichung, dass es sich um endliche Orthogonalsysteme mit Gewichtsfunktionen vom Gammatyp handelt.

Approximate approximations for the Poisson and the Stokes equations

The method of approximate approximations is based on generating functions representing an approximate partition of the unity, only. In the present paper this method is used for the numerical solution of the Poisson equation and the Stokes system in R^n (n = 2, 3). The corresponding approximate volume potentials will be computed explicitly in these cases, containing a one-dimensional integral, only. Numerical simulations show the efficiency of the method and confirm the expected convergence of essentially second order, depending on the smoothness of the data.

Künstliche Randbedingungen und Algebraische Äquivalenzen in der Elastizitätstheorie

In dieser Arbeit werden zwei Aspekte bei Randwertproblemen der linearen Elastizitätstheorie untersucht: die Approximation von Lösungen auf unbeschränkten Gebieten und die Änderung von Symmetrieklassen unter speziellen Transformationen. Ausgangspunkt der Dissertation ist das von Specovius-Neugebauer und Nazarov in "Artificial boundary conditions for Petrovsky systems of second order in exterior domains and in other domains of conical type"(Math. Meth. Appl. Sci, 2004; 27) eingeführte Verfahren zur Untersuchung von Petrovsky-Systemen zweiter Ordnung in Außenraumgebieten und Gebieten mit konischen Ausgängen mit Hilfe der Methode der künstlichen Randbedingungen. Dabei werden für die Ermittlung von Lösungen der Randwertprobleme die unbeschränkten Gebiete durch das Abschneiden mit einer Kugel beschränkt, und es wird eine künstliche Randbedingung konstruiert, um die Lösung des Problems möglichst gut zu approximieren. Das Verfahren wird dahingehend verändert, dass das abschneidende Gebiet ein Polyeder ist, da es für die Lösung des Approximationsproblems mit üblichen Finite-Element-Diskretisierungen von Vorteil sei, wenn das zu triangulierende Gebiet einen polygonalen Rand besitzt. Zu Beginn der Arbeit werden die wichtigsten funktionalanalytischen Begriffe und Ergebnisse der Theorie elliptischer Differentialoperatoren vorgestellt. Danach folgt der Hauptteil der Arbeit, der sich in drei Bereiche untergliedert. Als erstes wird für abschneidende Polyedergebiete eine formale Konstruktion der künstlichen Randbedingungen angegeben. Danach folgt der Nachweis der Existenz und Eindeutigkeit der Lösung des approximativen Randwertproblems auf dem abgeschnittenen Gebiet und im Anschluss wird eine Abschätzung für den resultierenden Abschneidefehler geliefert. An die theoretischen Ausführungen schließt sich die Betrachtung von Anwendungsbereiche an. Hier werden ebene Rissprobleme und Polarisationsmatrizen dreidimensionaler Außenraumprobleme der Elastizitätstheorie erläutert. Der letzte Abschnitt behandelt den zweiten Aspekt der Arbeit, den Bereich der Algebraischen Äquivalenzen. Hier geht es um die Transformation von Symmetrieklassen, um die Kenntnis der Fundamentallösung der Elastizitätsprobleme für transversalisotrope Medien auch für Medien zu nutzen, die nicht von transversalisotroper Struktur sind. Eine allgemeine Darstellung aller Klassen konnte hier nicht geliefert werden. Als Beispiel für das Vorgehen wird eine Klasse von orthotropen Medien im dreidimensionalen Fall angegeben, die sich auf den Fall der Transversalisotropie reduzieren lässt.

Characterization theorem for classical orthogonal polynomials on non-uniform lattices: The functional approach

Using the functional approach, we state and prove a characterization theorem for classical orthogonal polynomials on non-uniform lattices (quadratic lattices of a discrete or a q-discrete variable) including the Askey-Wilson polynomials. This theorem proves the equivalence between seven characterization properties, namely the Pearson equation for the linear functional, the second-order divided-difference equation, the orthogonality of the derivatives, the Rodrigues formula, two types of structure relations,and the Riccati equation for the formal Stieltjes function.

Untersuchung der Adsorptionskinetik photoschaltbarer azobenzolfunktionalisierter selbstorganisierter Monolagen mit optischer Frequenzverdopplung

In dieser Arbeit wurde das Adsorptionsverhalten zweier azobenzolfunktionalisierter Liganden auf Goldoberflächen untersucht. Diese Liganden zeichnen sich dadurch aus, dass sie mithilfe von Licht bestimmter Wellenlängen zwischen zwei Isomerisierungszuständen – sowohl in Lösung als auch in den auf der Oberfläche resultierenden Monolage – hin und her geschaltet werden können. Somit ist es möglich, Oberflächen herzustellen, deren physikalische und chemische Eigenschaften zwischen zwei Zuständen variiert werden können. Die Messungen des Adsorptionsverhaltens wurden mittels optischer Frequenzverdopplung durchgeführt. Diese Messmethode ist höchst grenzflächensensitiv und ermöglicht es somit die Adsorption der Liganden in situ und ich Echtzeit zu verfolgen. Neben den Adsorptionsmessungen wurde auch die Phase des frequenzverdoppelten Signals über eine Interferenzmethode gemessen. Die Ergebnisse dieser Phasenmessungen ermöglichen es, eine Aussage über eine mögliche Nichtlinearität der untersuchten Moleküle zu treffen. An die in den Adsorptionsmessungen gewonnenen Messdaten wurden drei kinetische Standardmodelle angepasst. Beschreibt eines dieser Modelle den im Experiment bestimmten Adsorptionsverlauf, kann eine Aussage über die zugrunde liegenden Prozesse des Adsorptionsvorganges getroffen werden. Die Ergebnisse der Adsorptionsmessungen zeigen einen deutlichen Einfluss des Isomerisierungszustandes der Liganden auf den Verlauf der Adsorption. Liegen die Moleküle im geschalteten Zustand vor, so verläuft die Adsorption langsamer. Weiterhin konnte gezeigt werden, dass ebenso intermolekulare Wechselwirkungen über Wasserstoffbrückenbindungen einen verlangsamenden Einfluss auf die Adsorption der Liganden haben. In den durchgeführten Phasenmessungen zeigte sich darüber hinaus, dass Liganden, die über an die Azobenzolgruppe angebundene Amidgruppen verfügen, eine Nichtlinearität aufweisen. Diese Nichtlinearität ist zudem vom Isomerisierungszustand der Liganden abhängig. In den kinetischen Untersuchungen konnte darüber hinaus gezeigt werden, dass sich die Adsorption der Liganden bis auf eine Ausnahme durch die Langmuirkinetik 2. Ordnung beschreiben lässt. Somit handelt es ich bei der Adsorption der untersuchten Liganden um eine Reaktion, der eine Bindungsspaltung voran geht. Dieser Befund konnte durch Vergleich mit weiteren Experimenten bestätigt werden.

Time Discretization of the SST-generalized Navier-Stokes Equations: Positive and Negative Results

In the theory of the Navier-Stokes equations, the proofs of some basic known results, like for example the uniqueness of solutions to the stationary Navier-Stokes equations under smallness assumptions on the data or the stability of certain time discretization schemes, actually only use a small range of properties and are therefore valid in a more general context. This observation leads us to introduce the concept of SST spaces, a generalization of the functional setting for the Navier-Stokes equations. It allows us to prove (by means of counterexamples) that several uniqueness and stability conjectures that are still open in the case of the Navier-Stokes equations have a negative answer in the larger class of SST spaces, thereby showing that proof strategies used for a number of classical results are not sufficient to affirmatively answer these open questions. More precisely, in the larger class of SST spaces, non-uniqueness phenomena can be observed for the implicit Euler scheme, for two nonlinear versions of the Crank-Nicolson scheme, for the fractional step theta scheme, and for the SST-generalized stationary Navier-Stokes equations. As far as stability is concerned, a linear version of the Euler scheme, a nonlinear version of the Crank-Nicolson scheme, and the fractional step theta scheme turn out to be non-stable in the class of SST spaces. The positive results established in this thesis include the generalization of classical uniqueness and stability results to SST spaces, the uniqueness of solutions (under smallness assumptions) to two nonlinear versions of the Euler scheme, two nonlinear versions of the Crank-Nicolson scheme, and the fractional step theta scheme for general SST spaces, the second order convergence of a version of the Crank-Nicolson scheme, and a new proof of the first order convergence of the implicit Euler scheme for the Navier-Stokes equations. For each convergence result, we provide conditions on the data that guarantee the existence of nonstationary solutions satisfying the regularity assumptions needed for the corresponding convergence theorem. In the case of the Crank-Nicolson scheme, this involves a compatibility condition at the corner of the space-time cylinder, which can be satisfied via a suitable prescription of the initial acceleration.

Optical Antennas for Biosensing Applications

The aim of the thesis is to theoretically investigate optical/plasmonic antennas for biosensing applications. The full 3-D numerical electromagnetic simulations have been performed by using finite integration technique (FIT). The electromagnetic properties of surface plasmon polaritons (SPPs) and the localized surface plasmons (LSPs) based devices are studied for sensing purpose. The surface plasmon resonance (SPR) biosensors offer high refractive index sensitivity at a fixed wavelength but are not enough for the detection of low concentrations of molecules. It has been demonstrated that the sensitivity of SPR sensors can be increased by employing the transverse magneto-optic Kerr effect (TMOKE) in combination with SPPs. The sensor based on the phenomena of TMOKE and SPPs are known as magneto-optic SPR (MOSPR) sensors. The optimized MOSPR sensor is analyzed which provides 8 times higher sensitivity than the SPR sensor, which will be able to detect lower concentration of molecules. But, the range of the refractive index detection is limited, due to the rapid decay of the amplitude of the MOSPR-signal with the increase of the refractive indices. Whereas, LSPs based sensors can detect lower concentrations of molecules, but their sensitivity is small at a fixed wavelength. Therefore, another device configuration known as perfect plasmonic absorber (PPA) is investigated which is based on the phenomena of metal-insulator-metal (MIM) waveguide. The PPA consists of a periodic array of gold nanoparticles and a thick gold film separated by a dielectric spacer. The electromagnetic modes of the PPA system are analyzed for sensing purpose. The second order mode of the PPA at a fixed wavelength has been proposed for the first time for biosensing applications. The PPA based sensor combines the properties of the LSPR sensor and the SPR sensor, for example, it illustrates increment in sensitivity of the LSPR sensor comparable to the SPR and can detect lower concentration of molecules due to the presence of nanoparticles.

Support Vector Machines: Training and Applications

The Support Vector Machine (SVM) is a new and very promising classification technique developed by Vapnik and his group at AT&T Bell Labs. This new learning algorithm can be seen as an alternative training technique for Polynomial, Radial Basis Function and Multi-Layer Perceptron classifiers. An interesting property of this approach is that it is an approximate implementation of the Structural Risk Minimization (SRM) induction principle. The derivation of Support Vector Machines, its relationship with SRM, and its geometrical insight, are discussed in this paper. Training a SVM is equivalent to solve a quadratic programming problem with linear and box constraints in a number of variables equal to the number of data points. When the number of data points exceeds few thousands the problem is very challenging, because the quadratic form is completely dense, so the memory needed to store the problem grows with the square of the number of data points. Therefore, training problems arising in some real applications with large data sets are impossible to load into memory, and cannot be solved using standard non-linear constrained optimization algorithms. We present a decomposition algorithm that can be used to train SVM's over large data sets. The main idea behind the decomposition is the iterative solution of sub-problems and the evaluation of, and also establish the stopping criteria for the algorithm. We present previous approaches, as well as results and important details of our implementation of the algorithm using a second-order variant of the Reduced Gradient Method as the solver of the sub-problems. As an application of SVM's, we present preliminary results we obtained applying SVM to the problem of detecting frontal human faces in real images.

An Immersed Interface Method for the Incompressible Navier-Stokes Equations in Irregular Domains

We present an immersed interface method for the incompressible Navier Stokes equations capable of handling rigid immersed boundaries. The immersed boundary is represented by a set of Lagrangian control points. In order to guarantee that the no-slip condition on the boundary is satisfied, singular forces are applied on the fluid at the immersed boundary. The forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines. The strength of singular forces is determined by solving a small system of equations at each time step. The Navier-Stokes equations are discretized on a staggered Cartesian grid by a second order accurate projection method for pressure and velocity.