Objective grading of the anterior eye

Purpose. To convert objective image analysis of anterior ocular surfaces into recognisable clinical grades, in order to provide a more sensitive and reliable equivalent to current subjective grading methods; a prospective, randomized study correlating clinical grading with digital image assessment. Methods. The possible range of clinical presentations Of bulbar and palpebral hyperaemia, palpebral roughness and corneal staining were represented by 4 sets of 10 images. The images were displayed in random order and graded by 50 clinicians using both subjective CCLRU and Efron grading scales. Previously validated objective image analysis was performed 3 times oil each of the 40 images. Digital measures included edge-detection and relative-coloration components. Step-wise regression analysis determined correlations between the average subjective grade and the objective image analysis measures. Results. Average subjective grades Could be predicted by a combination of the objective image analysis components. These digital ``grades'' accounted for between 69%, (for Efron scale-graded palpebral redness) and 98% (for Efron scale-graded bulbar hyperaemia) of the subjective variance. Conclusions. The results indicate that clinicians may use a combination of vessel areas and overall hue in their judgment of clinical severity for certain conditions. Objective grading call take these aspects into account, and be used to predict an average ``objective grade'' to be used by a clinician in describing the anterior eye. These measures are more sensitive and reliable than subjective grading while still utilizing familiar terminology, and can be applied in research or practice to improve the detection, and monitoring of ocular surface changes.

The design of a real-time nowcasting system for localised weather

National meteorological offices are largely concerned with synoptic-scale forecasting where weather predictions are produced for a whole country for 24 hours ahead. In practice, many local organisations (such as emergency services, construction industries, forestry, farming, and sports) require only local short-term, bespoke, weather predictions and warnings. This thesis shows that the less-demanding requirements do not require exceptional computing power and can be met by a modern, desk-top system which monitors site-specific ground conditions (such as temperature, pressure, wind speed and direction, etc) augmented with above ground information from satellite images to produce `nowcasts'. The emphasis in this thesis has been towards the design of such a real-time system for nowcasting. Local site-specific conditions are monitored using a custom-built, stand alone, Motorola 6809 based sub-system. Above ground information is received from the METEOSAT 4 geo-stationary satellite using a sub-system based on a commercially available equipment. The information is ephemeral and must be captured in real-time. The real-time nowcasting system for localised weather handles the data as a transparent task using the limited capabilities of the PC system. Ground data produces a time series of measurements at a specific location which represents the past-to-present atmospheric conditions of the particular site from which much information can be extracted. The novel approach adopted in this thesis is one of constructing stochastic models based on the AutoRegressive Integrated Moving Average (ARIMA) technique. The satellite images contain features (such as cloud formations) which evolve dynamically and may be subject to movement, growth, distortion, bifurcation, superposition, or elimination between images. The process of extracting a weather feature, following its motion and predicting its future evolution involves algorithms for normalisation, partitioning, filtering, image enhancement, and correlation of multi-dimensional signals in different domains. To limit the processing requirements, the analysis in this thesis concentrates on an `area of interest'. By this rationale, only a small fraction of the total image needs to be processed, leading to a major saving in time. The thesis also proposes an extention to an existing manual cloud classification technique for its implementation in automatically classifying a cloud feature over the `area of interest' for nowcasting using the multi-dimensional signals.

The characterisation of clean metal surfaces

The most perfectly structured metal surface observed in practice is that of a field evaporated field-ion microscope specimen. This surface has been characterised by adopting various optical analogue techniques. Hence a relationship has been determined between the structure of a single plane on the surface of a field-ion emitter and the geometry of a binary zone plate. By relating the known focussing properties of such a zone plate to those obtained from the projected images of such planes in a field-ion micrograph, it is possible to extract new information regarding the local magnification of the image. Further to this, it has been shown that the entire system of planes comprising the field-ion imaging surface may be regarded as a moire pattern formed between over-lapping zone plates. The properties of such moire zone plates are first established in an analysis of the moire pattern formed between zone plates on a flat surface. When these ideas are applied to the field-ion image it becomes possible to deduce further information regarding the precise topography of the emitter. It has also become possible to simulate differently proJected field-ion images by overlapping suitably aberrated zone plates. Low-energy ion bombardment is an essential preliminary to much surface research as a means of producing chemically clean surfaces. Hence it is important to know the nature and distribution of the resultant lattice damage, and the extent to which it may be removed by annealing. The field-ion microscope has been used to investigate such damage because its characterisation lies on the atomic scale. The present study is concerned with the in situ sputtering of tungsten emitters using helium, neon, argon and xenon ions with energies in the range 100eV to 1keV, together with observations of the effect of annealing. The relevance of these results to surface cleaning schedules is discussed.

Application of thermal neutron radiography for the mass transport of moisture through freezing soil

This thesis reports on the development of a technique to evaluate hydraulic conductivities in a soil (Snowcal) subject to freezing conditions. The technique draws on three distinctly different disciplines, Nuclear Physics, Soil Physics and Remote Sensing to provide a non-destructive and reliable evaluation of hydraulic conductivity throughout a freezing test. Thermal neutron radiography is used to provide information on local water/ice contents at anytime throughout the test. The experimental test rig is designed so that the soil matrix can be radiated by a neutron beam, from a nuclear reactor, to obtain radiographs. The radiographs can then be interpreted, following a process of remote sensing image enhancement, to yield information on relative water/ice contents. Interpretation of the radiographs is accommodated using image analysis equipment capable of distinguishing between 256 shades of grey. Remote sensing image enhancing techniques are then employed to develop false colour images which show the movement of water and development of ice lenses in the soil. Instrumentation is incorporated in the soil in the form of psychrometer/thermocouples, to record water potential, electrical resistance probes to enable ice and water to be differentiated on the radiographs and thermocouples to record the temperature gradient. Water content determinations are made from the enhanced images and plotted against potential measurements to provide the moisture characteristic for the soil. With relevant mathematical theory pore water distributions are obtained and combined with water content data to give hydraulic conductivities. The values for hydraulic conductivity in the saturated soil and at the frozen fringe are compared with established values for silts and silty-sands. The values are in general agreement and, with refinement, this non-destructive technique could afford useful information on a whole range of soils. The technique is of value over other methods because ice lenses are actually seen forming in the soil, supporting the accepted theories of frost action. There are economic and experimental restraints to the work which are associated with the use of a nuclear facility, however, the technique is versatile and has been applied to the study of moisture transfer in porous building materials and could be further developed into other research areas.

Pupillary shapes in response to electrical stimulation of the sclera of peripheral cornea in cats and porcines

Purpose: We have reported that the changes in the pupillary shape in response to electrical stimulation of the branches of the ciliary nerves in cats. (Miyagawa et al. PLoS One, 2014). This study investigates the changes in the pupillary shapes in response to electrical stimulations of the sclera of peripheral cornea in cats and porcines. Methods: Two enucleated eyes of two cats and three enucleated porcine eyes were studied. Trains of biphasic pulses (current, 3 mA; duration, 2 ms/phase; frequency, 40 Hz) were applied using a tungsten electrode (0.3mm diameter). The stimulation was performed at every 45 degree over the entire circular region on the sclera near the cornea. The pupillary images were recorded before and 4 s (cat) and 10 s (pig) after the stimulation and the change in the pupil diameter (Δr) was quantified. The pupillary images were obtained with a custom-built compact wavefront aberrometer (Uday et al. J Cataract Refract Surg, 2013). Results: In a cat eye, the pupil was dilated by the electrical stimulation at six out of eight orientations (before stimulation pupil diameter r=10.10±0.49 mm, Δr=0.33±0.12 mm). The pupil dilated only toward the electrode (relative eccentricity of the pupil center to the pupil diameter change amount rdec=1.15±0.28). In the porcine eyes, the pupils were constricted by the electrical stimulations at the temporal and nasal orientations (r=10.04±0.57 mm, Δr=1.52±0.70 mm). The pupils contracted symmetrically (rdec=0.30±0.12). Conclusions: With electrical stimulation in the sclera of the peripheral cornea, asymmetric mydriasis in cat eyes and symmetrical miosis in porcine eyes were observed. Under the assumption that the electrical stimulation stimulated both muscles that contribute to the pupil control, our hypothesis proposed here is that the pupil dilator is stronger than the pupil sphincter in cat, and pupil sphincter is stronger than pupil dilator in porcine.

Transmission of image data over digital networks involving mobile terminals

There is a growing demand for data transmission over digital networks involving mobile terminals. An important class of data required for transmission over mobile terminals is image information such as street maps, floor plans and identikit images. This sort of transmission is of particular interest to the service industries such as the Police force, Fire brigade, medical services and other services. These services cannot be applied directly to mobile terminals because of the limited capacity of the mobile channels and the transmission errors caused by the multipath (Rayleigh) fading. In this research, transmission of line diagram images such as floor plans and street maps, over digital networks involving mobile terminals at transmission rates of 2400 bits/s and 4800 bits/s have been studied. A low bit-rate source encoding technique using geometric codes is found to be suitable to represent line diagram images. In geometric encoding, the amount of data required to represent or store the line diagram images is proportional to the image detail. Thus a simple line diagram image would require a small amount of data. To study the effect of transmission errors due to mobile channels on the transmitted images, error sources (error files), which represent mobile channels under different conditions, have been produced using channel modelling techniques. Satisfactory models of the mobile channel have been obtained when compared to the field test measurements. Subjective performance tests have been carried out to evaluate the quality and usefulness of the received line diagram images under various mobile channel conditions. The effect of mobile transmission errors on the quality of the received images has been determined. To improve the quality of the received images under various mobile channel conditions, forward error correcting codes (FEC) with interleaving and automatic repeat request (ARQ) schemes have been proposed. The performance of the error control codes have been evaluated under various mobile channel conditions. It has been shown that a FEC code with interleaving can be used effectively to improve the quality of the received images under normal and severe mobile channel conditions. Under normal channel conditions, similar results have been obtained when using ARQ schemes. However, under severe mobile channel conditions, the FEC code with interleaving shows better performance.

OCT retinal thickness ratios and frequency of follow-up in Wet AMD

Poster Introduction: In neovascular age-related macular degeneration (nAMD), optical coherence tomography (OCT) is an important tool to determine when intravitreal injections of ranibizumab should be administered. Current guidelines recommend that patients should be reviewed four weekly and OCT indications for further treatment include subretinal fluid and intraretinal fluid or cysts. Purpose: We have reviewed the OCT scans of subjects who have successfully responded to ranibizumab to look for factors that might predict which patients will not require injection and could have extended appointments. Method: This was a prospective study in which we observed for 6 consecutive months the OCT images of 28 subjects who had received intravitreal ranibizumab for nAMD and were judged to be clinically inactive at recruitment to the study. Ratios between full retinal thickness (FRT = neurosensory retina + outer reflective band) and outer reflective band (ORB) thickness at the fovea were calculated for each subject at the moment of entering the study and at each successive visit for 6 consecutive months. Results: Patients with lower FRT/ORB ratios were found to be less likely to require an additional injection of ranibizumab and no subject with a ratio of 1.75 or less needed further injections. Conclusion: This small pilot study suggests that on macular OCT, the FRT/ORB ratio, and in particular values of 1.75 or less, may prove to be a useful, practical tool when deciding the follow up period for subjects undergoing treatment with intravitreal ranibizumab for nAMD.

Consensus on recording of gas permeable contact lens fit

Purpose: To develop a new schematic scheme for efficiently recording the key parameters of gas permeable contact lens (GP) fits based on current consensus. Methods: Over 100 established GP fitters and educators met to discuss the parameters proposed in educational material for evaluating GP fit and concluded on the key parameters that should be recorded. The accuracy and variability of evaluating the fluorescein pattern of GP fit was determined by having 35 experienced contact lens practitioners from across the world, grading 5 images of a range of fits and the topographer simulation of the same fits, in random, order using the proposed scheme. The accuracy of the grading was compared to objective image analysis of the fluorescein intensity of the same images. Results: The key information to record to adequately describe the fit of an GP was agreed as: the manufacturer, brand and lens parameters; settling time; comfort on a 5 point scale; centration; movement on blink on a ±2 scale; and the Primary Fluorescein Pattern in the central, mid-peripheral and edge regions of the lens averaged along the horizontal and vertical lens axes, on a ±2 scale. On average 50-60% of practitioners selected the median grade when subjectively rating fluorescein intensity and this was correlated to objective quantification (r= 0.602, p< 0.001). Objective grading suggesting horizontal median fluorescein intensity was generally symmetrical, as was the vertical meridian, but this was not the case for subjective grading. Simulated fluorescein patterns were subjectively and objectively graded as being less intense than real photographs (p< 0.01). Conclusion: GP fit recording can be standardised and simplified to enhance GP practice. © 2013 British Contact Lens Association.

Advances in characterisation, calibration and data processing speed of optical coherence tomography systems

This thesis describes advances in the characterisation, calibration and data processing of optical coherence tomography (OCT) systems. Femtosecond (fs) laser inscription was used for producing OCT-phantoms. Transparent materials are generally inert to infra-red radiations, but with fs lasers material modification occurs via non-linear processes when the highly focused light source interacts with the materials. This modification is confined to the focal volume and is highly reproducible. In order to select the best inscription parameters, combination of different inscription parameters were tested, using three fs laser systems, with different operating properties, on a variety of materials. This facilitated the understanding of the key characteristics of the produced structures with the aim of producing viable OCT-phantoms. Finally, OCT-phantoms were successfully designed and fabricated in fused silica. The use of these phantoms to characterise many properties (resolution, distortion, sensitivity decay, scan linearity) of an OCT system was demonstrated. Quantitative methods were developed to support the characterisation of an OCT system collecting images from phantoms and also to improve the quality of the OCT images. Characterisation methods include the measurement of the spatially variant resolution (point spread function (PSF) and modulation transfer function (MTF)), sensitivity and distortion. Processing of OCT data is a computer intensive process. Standard central processing unit (CPU) based processing might take several minutes to a few hours to process acquired data, thus data processing is a significant bottleneck. An alternative choice is to use expensive hardware-based processing such as field programmable gate arrays (FPGAs). However, recently graphics processing unit (GPU) based data processing methods have been developed to minimize this data processing and rendering time. These processing techniques include standard-processing methods which includes a set of algorithms to process the raw data (interference) obtained by the detector and generate A-scans. The work presented here describes accelerated data processing and post processing techniques for OCT systems. The GPU based processing developed, during the PhD, was later implemented into a custom built Fourier domain optical coherence tomography (FD-OCT) system. This system currently processes and renders data in real time. Processing throughput of this system is currently limited by the camera capture rate. OCTphantoms have been heavily used for the qualitative characterization and adjustment/ fine tuning of the operating conditions of OCT system. Currently, investigations are under way to characterize OCT systems using our phantoms. The work presented in this thesis demonstrate several novel techniques of fabricating OCT-phantoms and accelerating OCT data processing using GPUs. In the process of developing phantoms and quantitative methods, a thorough understanding and practical knowledge of OCT and fs laser processing systems was developed. This understanding leads to several novel pieces of research that are not only relevant to OCT but have broader importance. For example, extensive understanding of the properties of fs inscribed structures will be useful in other photonic application such as making of phase mask, wave guides and microfluidic channels. Acceleration of data processing with GPUs is also useful in other fields.

Sucrose cryo-protection facilitates imaging of whole eye sections by MALDI mass spectrometry

Sucrose is used as a cryo-preservation agent on large mammalian eyes post formalin fixation and is shown to reduce freezing artefacts allowing the collection of 12-μm thick sections from these large aqueous samples. The suitability of this technique for use in MALDI imaging experiments is demonstrated by the acquisition of the first images of lipid distributions within whole sagittal porcine eye sections. © 2012 John Wiley & Sons, Ltd.

Visualisation and browsing of image databases

In this chapter we provide a comprehensive overview of the emerging field of visualising and browsing image databases. We start with a brief introduction to content-based image retrieval and the traditional query-by-example search paradigm that many retrieval systems employ. We specify the problems associated with this type of interface, such as users not being able to formulate a query due to not having a target image or concept in mind. The idea of browsing systems is then introduced as a means to combat these issues, harnessing the cognitive power of the human mind in order to speed up image retrieval.We detail common methods in which the often high-dimensional feature data extracted from images can be used to visualise image databases in an intuitive way. Systems using dimensionality reduction techniques, such as multi-dimensional scaling, are reviewed along with those that cluster images using either divisive or agglomerative techniques as well as graph-based visualisations. While visualisation of an image collection is useful for providing an overview of the contained images, it forms only part of an image database navigation system. We therefore also present various methods provided by these systems to allow for interactive browsing of these datasets. A further area we explore are user studies of systems and visualisations where we look at the different evaluations undertaken in order to test usability and compare systems, and highlight the key findings from these studies. We conclude the chapter with several recommendations for future work in this area. © 2011 Springer-Verlag Berlin Heidelberg.

On the relation between dichoptic masking and binocular rivalry

When our two eyes view incompatible images, the brain invokes suppressive processes to inhibit one image, and favor the other. Two phenomena are typically observed: dichoptic masking (reduced sensitivity to one image) for brief presentations, and binocular rivalry (alternation between the two images), over longer exposures. However, it is not clear if these two phenomena arise from a common suppressive process. We investigated this by measuring both threshold elevation in simultaneous dichoptic masking and mean percept durations in rivalry, whilst varying relative stimulus orientation. Masking and rivalry showed significant correlations, such that strong masking was associated with long dominance durations. A second experiment suggested that individual differences across both measures are also correlated. These findings are consistent with varying the magnitude of interocular suppression in computational models of both rivalry and masking, and imply the existence of a common suppressive process. Since dichoptic masking has been localised to the monocular neurons of V1, this is a plausible first stage of binocular rivalry.

A verifiable solution to the MEG inverse problem

Magnetoencephalography (MEG) is a non-invasive brain imaging technique with the potential for very high temporal and spatial resolution of neuronal activity. The main stumbling block for the technique has been that the estimation of a neuronal current distribution, based on sensor data outside the head, is an inverse problem with an infinity of possible solutions. Many inversion techniques exist, all using different a-priori assumptions in order to reduce the number of possible solutions. Although all techniques can be thoroughly tested in simulation, implicit in the simulations are the experimenter's own assumptions about realistic brain function. To date, the only way to test the validity of inversions based on real MEG data has been through direct surgical validation, or through comparison with invasive primate data. In this work, we constructed a null hypothesis that the reconstruction of neuronal activity contains no information on the distribution of the cortical grey matter. To test this, we repeatedly compared rotated sections of grey matter with a beamformer estimate of neuronal activity to generate a distribution of mutual information values. The significance of the comparison between the un-rotated anatomical information and the electrical estimate was subsequently assessed against this distribution. We found that there was significant (P < 0.05) anatomical information contained in the beamformer images across a number of frequency bands. Based on the limited data presented here, we can say that the assumptions behind the beamformer algorithm are not unreasonable for the visual-motor task investigated.