Quantitative analysis of multi-spectral fundus images

We have developed a new technique for extracting histological parameters from multi-spectral images of the ocular fundus. The new method uses a Monte Carlo simulation of the reflectance of the fundus to model how the spectral reflectance of the tissue varies with differing tissue histology. The model is parameterised by the concentrations of the five main absorbers found in the fundus: retinal haemoglobins, choroidal haemoglobins, choroidal melanin, RPE melanin and macular pigment. These parameters are shown to give rise to distinct variations in the tissue colouration. We use the results of the Monte Carlo simulations to construct an inverse model which maps tissue colouration onto the model parameters. This allows the concentration and distribution of the five main absorbers to be determined from suitable multi-spectral images. We propose the use of "image quotients" to allow this information to be extracted from uncalibrated image data. The filters used to acquire the images are selected to ensure a one-to-one mapping between model parameters and image quotients. To recover five model parameters uniquely, images must be acquired in six distinct spectral bands. Theoretical investigations suggest that retinal haemoglobins and macular pigment can be recovered with RMS errors of less than 10%. We present parametric maps showing the variation of these parameters across the posterior pole of the fundus. The results are in agreement with known tissue histology for normal healthy subjects. We also present an early result which suggests that, with further development, the technique could be used to successfully detect retinal haemorrhages.

Biochemical sensor based on a novel all-fibre cavity ring down spectroscopy technique incorporating a tilted fibre Bragg grating

A novel all-fibre cavity ring down spectroscopy technique is proposed where a tilt fibre Bragg grating (TFBG) or long-period grating (LPG) in the cavity provides sensitivity to surrounding medium. Such configuration with an LPG as the representative was theoretically analyzed. Two spectral bands were identified employable for sensing of surrounding refractive index for a weak LPG while only one band existed for a strong LPG. A TFBG, with enhanced sensitivity compared to usual LPGs, was used in a ring down cavity of 1 m constructed with 2 fibre Bragg gratings as the reflectors and the decay time changed from 220 to 450 ns when the TFBG was immersed into water from air.

The integration of multi-source data to improve the classification of remotely sensed images

The number of remote sensing platforms and sensors rises almost every year, yet much work on the interpretation of land cover is still carried out using either single images or images from the same source taken at different dates. Two questions could be asked of this proliferation of images: can the information contained in different scenes be used to improve the classification accuracy and, what is the best way to combine the different imagery? Two of these multiple image sources are MODIS on the Terra platform and ETM+ on board Landsat7, which are suitably complementary. Daily MODIS images with 36 spectral bands in 250-1000 m spatial resolution and seven spectral bands of ETM+ with 30m and 16 days spatial and temporal resolution respectively are available. In the UK, cloud cover may mean that only a few ETM+ scenes may be available for any particular year and these may not be at the time of year of most interest. The MODIS data may provide information on land cover over the growing season, such as harvest dates, that is not present in the ETM+ data. Therefore, the primary objective of this work is to develop a methodology for the integration of medium spatial resolution Landsat ETM+ image, with multi-temporal, multi-spectral, low-resolution MODIS \Terra images, with the aim of improving the classification of agricultural land. Additionally other data may also be incorporated such as field boundaries from existing maps. When classifying agricultural land cover of the type seen in the UK, where crops are largely sown in homogenous fields with clear and often mapped boundaries, the classification is greatly improved using the mapped polygons and utilising the classification of the polygon as a whole as an apriori probability in classifying each individual pixel using a Bayesian approach. When dealing with multiple images from different platforms and dates it is highly unlikely that the pixels will be exactly co-registered and these pixels will contain a mixture of different real world land covers. Similarly the different atmospheric conditions prevailing during the different days will mean that the same emission from the ground will give rise to different sensor reception. Therefore, a method is presented with a model of the instantaneous field of view and atmospheric effects to enable different remote sensed data sources to be integrated.

Analysis of magnetoencephaoographic data as a nonlinear dynamical system

This thesis presents the results from an investigation into the merits of analysing Magnetoencephalographic (MEG) data in the context of dynamical systems theory. MEG is the study of both the methods for the measurement of minute magnetic flux variations at the scalp, resulting from neuro-electric activity in the neocortex, as well as the techniques required to process and extract useful information from these measurements. As a result of its unique mode of action - by directly measuring neuronal activity via the resulting magnetic field fluctuations - MEG possesses a number of useful qualities which could potentially make it a powerful addition to any brain researcher's arsenal. Unfortunately, MEG research has so far failed to fulfil its early promise, being hindered in its progress by a variety of factors. Conventionally, the analysis of MEG has been dominated by the search for activity in certain spectral bands - the so-called alpha, delta, beta, etc that are commonly referred to in both academic and lay publications. Other efforts have centred upon generating optimal fits of "equivalent current dipoles" that best explain the observed field distribution. Many of these approaches carry the implicit assumption that the dynamics which result in the observed time series are linear. This is despite a variety of reasons which suggest that nonlinearity might be present in MEG recordings. By using methods that allow for nonlinear dynamics, the research described in this thesis avoids these restrictive linearity assumptions. A crucial concept underpinning this project is the belief that MEG recordings are mere observations of the evolution of the true underlying state, which is unobservable and is assumed to reflect some abstract brain cognitive state. Further, we maintain that it is unreasonable to expect these processes to be adequately described in the traditional way: as a linear sum of a large number of frequency generators. One of the main objectives of this thesis will be to prove that much more effective and powerful analysis of MEG can be achieved if one were to assume the presence of both linear and nonlinear characteristics from the outset. Our position is that the combined action of a relatively small number of these generators, coupled with external and dynamic noise sources, is more than sufficient to account for the complexity observed in the MEG recordings. Another problem that has plagued MEG researchers is the extremely low signal to noise ratios that are obtained. As the magnetic flux variations resulting from actual cortical processes can be extremely minute, the measuring devices used in MEG are, necessarily, extremely sensitive. The unfortunate side-effect of this is that even commonplace phenomena such as the earth's geomagnetic field can easily swamp signals of interest. This problem is commonly addressed by averaging over a large number of recordings. However, this has a number of notable drawbacks. In particular, it is difficult to synchronise high frequency activity which might be of interest, and often these signals will be cancelled out by the averaging process. Other problems that have been encountered are high costs and low portability of state-of-the- art multichannel machines. The result of this is that the use of MEG has, hitherto, been restricted to large institutions which are able to afford the high costs associated with the procurement and maintenance of these machines. In this project, we seek to address these issues by working almost exclusively with single channel, unaveraged MEG data. We demonstrate the applicability of a variety of methods originating from the fields of signal processing, dynamical systems, information theory and neural networks, to the analysis of MEG data. It is noteworthy that while modern signal processing tools such as independent component analysis, topographic maps and latent variable modelling have enjoyed extensive success in a variety of research areas from financial time series modelling to the analysis of sun spot activity, their use in MEG analysis has thus far been extremely limited. It is hoped that this work will help to remedy this oversight.

An investigation into the wavelength stability of polymer optical fibre Bragg gratings

The inscription of Bragg gratings has been demonstrated in PMMA-based polymer optical fibre. The water affinity of PMMA can introduce significant wavelength change in a polymer optical fibre Bragg grating (POFBG). In polymer optical fibre losses are much higher than with silica fibre. Very strong absorption bands related to higher harmonics of vibrations of the C-H bond dominate throughout the visible and near infrared. Molecular vibration in substances generates heat, which is referred to as the thermal effect of molecular vibration. This means that a large part of the absorption of optical energy in those spectral bands will convert into thermal energy, which eventually drives water content out of the polymer fibre and reduces the wavelength of POFBG. In this work we have investigated the wavelength stability of POFBGs in different circumstances. The experiment has shown that the characteristic wavelength of a POFBG starts decreasing after a light source is applied to it. This decrease continues until equilibrium inside the fibre is established, depending on the initial water content inside the fibre, the surrounding humidity, the optical power applied, and the fibre size. Our investigation has shown that POFBGs operating at around 850 nm show much smaller wavelength reduction than those operating at around 1550 nm in the same fibre; POFBGs with different diameters show different changes; POFBGs powered by a low level light source, or operating in a very dry environment are least affected by this thermal effect.

High-power versatile picosecond pulse generation from mode-locked quantum-dot laser diodes

This paper presents the current status of our research in mode-locked quantum-dot edge-emitting laser diodes, particularly highlighting the recent progress in spectral and temporal versatility of both monolithic and external-cavity laser configurations. Spectral versatility is demonstrated through broadband tunability and novel mode-locking regimes that involve distinct spectral bands, such as dual-wavelength mode-locking, and robust high-power wavelength bistability. Broad tunability of the pulse repetition rate is also demonstrated for an external-cavity mode-locked quantum-dot laser, revealing a nearly constant pulse peak power at different pulse repetition rates. High-energy and low-noise pulse generations are demonstrated for low-pulse repetition rates. These recent advances confirm the potential of quantum-dot lasers as versatile, compact, and low-cost sources of ultrashort pulses. © 2011 IEEE.

Random distributed feedback fibre lasers

The concept of random lasers exploiting multiple scattering of photons in an amplifying disordered medium in order to generate coherent light without a traditional laser resonator has attracted a great deal of attention in recent years. This research area lies at the interface of the fundamental theory of disordered systems and laser science. The idea was originally proposed in the context of astrophysics in the 1960s by V.S. Letokhov, who studied scattering with "negative absorption" of the interstellar molecular clouds. Research on random lasers has since developed into a mature experimental and theoretical field. A simple design of such lasers would be promising for potential applications. However, in traditional random lasers the properties of the output radiation are typically characterized by complex features in the spatial, spectral and time domains, making them less attractive than standard laser systems in terms of practical applications. Recently, an interesting and novel type of one-dimensional random laser that operates in a conventional telecommunication fibre without any pre-designed resonator mirrors-random distributed feedback fibre laser-was demonstrated. The positive feedback required for laser generation in random fibre lasers is provided by the Rayleigh scattering from the inhomogeneities of the refractive index that are naturally present in silica glass. In the proposed laser concept, the randomly backscattered light is amplified through the Raman effect, providing distributed gain over distances up to 100km. Although an effective reflection due to the Rayleigh scattering is extremely small (~0.1%), the lasing threshold may be exceeded when a sufficiently large distributed Raman gain is provided. Such a random distributed feedback fibre laser has a number of interesting and attractive features. The fibre waveguide geometry provides transverse confinement, and effectively one-dimensional random distributed feedback leads to the generation of a stationary near-Gaussian beam with a narrow spectrum. A random distributed feedback fibre laser has efficiency and performance that are comparable to and even exceed those of similar conventional fibre lasers. The key features of the generated radiation of random distributed feedback fibre lasers include: a stationary narrow-band continuous modeless spectrum that is free of mode competition, nonlinear power broadening, and an output beam with a Gaussian profile in the fundamental transverse mode (generated both in single mode and multi-mode fibres).This review presents the current status of research in the field of random fibre lasers and shows their potential and perspectives. We start with an introductory overview of conventional distributed feedback lasers and traditional random lasers to set the stage for discussion of random fibre lasers. We then present a theoretical analysis and experimental studies of various random fibre laser configurations, including widely tunable, multi-wavelength, narrow-band generation, and random fibre lasers operating in different spectral bands in the 1-1.6μm range. Then we discuss existing and future applications of random fibre lasers, including telecommunication and distributed long reach sensor systems. A theoretical description of random lasers is very challenging and is strongly linked with the theory of disordered systems and kinetic theory. We outline two key models governing the generation of random fibre lasers: the average power balance model and the nonlinear Schrödinger equation based model. Recently invented random distributed feedback fibre lasers represent a new and exciting field of research that brings together such diverse areas of science as laser physics, the theory of disordered systems, fibre optics and nonlinear science. Stable random generation in optical fibre opens up new possibilities for research on wave transport and localization in disordered media. We hope that this review will provide background information for research in various fields and will stimulate cross-disciplinary collaborations on random fibre lasers. © 2014 Elsevier B.V.

Multispectral fundus analysis

Purpose - To generate a reflectance model of the fundus that allows an accurate non-invasive quantification of blood and pigments. Methods - A Monte Carlo simulation was used to produce a mathematical model of light interaction with the fundus at different wavelengths. The model predictions were compared with fundus images from normal volunteers in several spectral bands (peaks at 507, 525, 552, 585, 596 and 611nm). Th e model was then used to calculate the concentration and distribution of the known absorbing components of the fundus. Results - The shape of the statistical distribution of the image data generally corresponded to that of the model data; the model however appears to overestimate the reflectance of the fundus in the longer wavelength region.As the absorption by xanthophyll has no significant eff ect on light transport above 534nm, its distribution in the fundus was quantified: the wavelengths where both shape and distribution of image and model data matched (<553nm) were used to train a neural network which was then applied to every point in the image data. The xanthophyll distribution thus found was in agreement with published literature data in normal subjects. Conclusion - We have developed a method for optimising multi-spectral imaging of the fundus and a computer image analysis capable of estimating information about the structure and properties of the fundus. Th e technique successfully calculates the distribution of xanthophyll in the fundus of healthy volunteers. Further improvement of the model is required to allow the deduction of other parameters from images; investigations in known pathology models are also necessary to establish if this method is of clinical use in detecting early chroido-retinopathies, hence providing a useful screening and diagnostic tool.

An investigation into the wavelength stability of polymer optical fibre Bragg gratings

The inscription of Bragg gratings has been demonstrated in PMMA-based polymer optical fibre. The water affinity of PMMA can introduce significant wavelength change in a polymer optical fibre Bragg grating (POFBG). In polymer optical fibre losses are much higher than with silica fibre. Very strong absorption bands related to higher harmonics of vibrations of the C-H bond dominate throughout the visible and near infrared. Molecular vibration in substances generates heat, which is referred to as the thermal effect of molecular vibration. This means that a large part of the absorption of optical energy in those spectral bands will convert into thermal energy, which eventually drives water content out of the polymer fibre and reduces the wavelength of POFBG. In this work we have investigated the wavelength stability of POFBGs in different circumstances. The experiment has shown that the characteristic wavelength of a POFBG starts decreasing after a light source is applied to it. This decrease continues until equilibrium inside the fibre is established, depending on the initial water content inside the fibre, the surrounding humidity, the optical power applied, and the fibre size. Our investigation has shown that POFBGs operating at around 850 nm show much smaller wavelength reduction than those operating at around 1550 nm in the same fibre; POFBGs with different diameters show different changes; POFBGs powered by a low level light source, or operating in a very dry environment are least affected by this thermal effect.

Raman gain and random distributed feedback generation in nitrogen doped silica core fiber

Recently, the concept of a random distributed feedback (DFB) lasing in optical fibers has been demonstrated [1], A number of different random DFB fiber lasers has been demonstrated so far including tunable, multiwalength, cascaded generation, generation in different spectral bands etc [2-7]. All systems are based on standard low-loss germanium doped silica core fibres having relatively low Rayleigh scattering coefficient. Thus, the typical length of random DFB fiber lasers is in the range from several kilometres to tens of kilometres to accumulate enough random feedback. Here we demonstrate for the first time to our knowledge the random DFB fiber laser based on a nitrogen doped silica core (N-doped) fiber. The fiber has several times higher Rayleigh scattering coefficient compared to standard telecommunication fibres. Thus, the generation is achieved in 500 meters long fiber only. © 2013 IEEE.

The spectral sensitivity of long period gratings fabricated in elliptical core D-shaped optical fibre

Long period gratings (LPGs) were written into a D-shaped optical fibre that has an elliptical core with a W-shaped refractive index profile and the first detailed investigation of such LPGs is presented. The LPGs’ attenuation bands were found to be sensitive to the polarisation of the interrogating light with a spectral separation of about 15 nm between the two orthogonal polarisation states. A finite element method was successfully used to model many of the behavioural features of the LPGs. In addition, two spectrally overlapping attenuation bands corresponding to orthogonal polarisation states were observed; modelling successfully reproduced this spectral feature. The spectral sensitivity of both orthogonal states was experimentally measured with respect to temperature and bending. These LPG devices produced blue and red wavelength shifts depending upon the orientation of the bend with measured maximum sensitivities of -3.56 and +6.51 nm m, suggesting that this type of fibre LPG may be useful as a shape/bend orientation sensor with reduced errors associated with polarisation dependence. The use of neighbouring bands to discriminate between temperature and bending was also demonstrated, leading to an overall curvature error of ±0.14 m-1 and an overall temperature error of ±0.3 °C with a maximum polarisation dependence error of ±8 × 10-2 m-1 for curvature and ±5 × 10-2 °C for temperature.

The intelligibility of noise-vocoded speech:spectral information available from across-channel comparison of amplitude envelopes

Noise-vocoded (NV) speech is often regarded as conveying phonetic information primarily through temporal-envelope cues rather than spectral cues. However, listeners may infer the formant frequencies in the vocal-tract output—a key source of phonetic detail—from across-band differences in amplitude when speech is processed through a small number of channels. The potential utility of this spectral information was assessed for NV speech created by filtering sentences into six frequency bands, and using the amplitude envelope of each band (=30 Hz) to modulate a matched noise-band carrier (N). Bands were paired, corresponding to F1 (˜N1 + N2), F2 (˜N3 + N4) and the higher formants (F3' ˜ N5 + N6), such that the frequency contour of each formant was implied by variations in relative amplitude between bands within the corresponding pair. Three-formant analogues (F0 = 150 Hz) of the NV stimuli were synthesized using frame-by-frame reconstruction of the frequency and amplitude of each formant. These analogues were less intelligible than the NV stimuli or analogues created using contours extracted from spectrograms of the original sentences, but more intelligible than when the frequency contours were replaced with constant (mean) values. Across-band comparisons of amplitude envelopes in NV speech can provide phonetically important information about the frequency contours of the underlying formants.

A near real-time water surface detection method based on HSV transformation of MODIS multi-Spectral time series data

In the face of global population growth and the uneven distribution of water supply, a better knowledge of the spatial and temporal distribution of surface water resources is critical. Remote sensing provides a synoptic view of ongoing processes, which addresses the intricate nature of water surfaces and allows an assessment of the pressures placed on aquatic ecosystems. However, the main challenge in identifying water surfaces from remotely sensed data is the high variability of spectral signatures, both in space and time. In the last 10 years only a few operational methods have been proposed to map or monitor surface water at continental or global scale, and each of them show limitations. The objective of this study is to develop and demonstrate the adequacy of a generic multi-temporal and multi-spectral image analysis method to detect water surfaces automatically, and to monitor them in near-real-time. The proposed approach, based on a transformation of the RGB color space into HSV, provides dynamic information at the continental scale. The validation of the algorithm showed very few omission errors and no commission errors. It demonstrates the ability of the proposed algorithm to perform as effectively as human interpretation of the images. The validation of the permanent water surface product with an independent dataset derived from high resolution imagery, showed an accuracy of 91.5% and few commission errors. Potential applications of the proposed method have been identified and discussed. The methodology that has been developed 27 is generic: it can be applied to sensors with similar bands with good reliability, and minimal effort. Moreover, this experiment at continental scale showed that the methodology is efficient for a large range of environmental conditions. Additional preliminary tests over other continents indicate that the proposed methodology could also be applied at the global scale without too many difficulties

The spectral sensitivity of long period gratings fabricated in elliptical core D-shaped optical fibre

Long period gratings (LPGs) were written into a D-shaped optical fibre that has an elliptical core with a W-shaped refractive index profile and the first detailed investigation of such LPGs is presented. The LPGs’ attenuation bands were found to be sensitive to the polarisation of the interrogating light with a spectral separation of about 15 nm between the two orthogonal polarisation states. A finite element method was successfully used to model many of the behavioural features of the LPGs. In addition, two spectrally overlapping attenuation bands corresponding to orthogonal polarisation states were observed; modelling successfully reproduced this spectral feature. The spectral sensitivity of both orthogonal states was experimentally measured with respect to temperature and bending. These LPG devices produced blue and red wavelength shifts depending upon the orientation of the bend with measured maximum sensitivities of -3.56 and +6.51 nm m, suggesting that this type of fibre LPG may be useful as a shape/bend orientation sensor with reduced errors associated with polarisation dependence. The use of neighbouring bands to discriminate between temperature and bending was also demonstrated, leading to an overall curvature error of ±0.14 m-1 and an overall temperature error of ±0.3 °C with a maximum polarisation dependence error of ±8 × 10-2 m-1 for curvature and ±5 × 10-2 °C for temperature.

Hypnotic induction is followed by state-like changes in the organization of EEG functional connectivity in the theta and beta frequency bands in high-hypnotically susceptible individuals

Altered state theories of hypnosis posit that a qualitatively distinct state of mental processing, which emerges in those with high hypnotic susceptibility following a hypnotic induction, enables the generation of anomalous experiences in response to specific hypnotic suggestions. If so then such a state should be observable as a discrete pattern of changes to functional connectivity (shared information) between brain regions following a hypnotic induction in high but not low hypnotically susceptible participants. Twenty-eight channel EEG was recorded from 12 high susceptible (highs) and 11 low susceptible (lows) participants with their eyes closed prior to and following a standard hypnotic induction. The EEG was used to provide a measure of functional connectivity using both coherence (COH) and the imaginary component of coherence (iCOH), which is insensitive to the effects of volume conduction. COH and iCOH were calculated between all electrode pairs for the frequency bands: delta (0.1-3.9 Hz), theta (4-7.9 Hz) alpha (8-12.9 Hz), beta1 (13-19.9 Hz), beta2 (20-29.9 Hz) and gamma (30-45 Hz). The results showed that there was an increase in theta iCOH from the pre-hypnosis to hypnosis condition in highs but not lows with a large proportion of significant links being focused on a central-parietal hub. There was also a decrease in beta1 iCOH from the pre-hypnosis to hypnosis condition with a focus on a fronto-central and an occipital hub that was greater in high compared to low susceptibles. There were no significant differences for COH or for spectral band amplitude in any frequency band. The results are interpreted as indicating that the hypnotic induction elicited a qualitative change in the organization of specific control systems within the brain for high as compared to low susceptible participants. This change in the functional organization of neural networks is a plausible indicator of the much theorized "hypnotic-state". © 2014 Jamieson and Burgess.