Monika Ewa Danielewska

Ultrasonic In Vivo Measurement of Ocular Surface Expansion

Our aim was to ascertain whether the ultrasonic measurement of longitudinal corneal apex displacements carried out in a proper headrest is a credible method of ocular pulse (OP) detection. To distinguish between longitudinal movements of the eye globe treated as a rigid body and ocular surface expansion caused by the variations of the eye-globe volume, two ultrasound distance sensors were applied to noninvasively measure displacements of cornea and sclera. The same sensors were used to examine the influence of the anterio-posterior movements of a fixed head on the registration of corneal apex pulsation. In both experiments, ECG signals were synchronically recorded. Time, spectral, and coherence analyses obtained for four healthy subjects showed that the ocular surface expansion due to pulsatile ocular blood flow (POBF) is the main component of longitudinal corneal displacement. Ocular surface pulsation is always affected by the head movement. However, there exist some unique properties of signals, which help to distinguish between head and eye movements. A rigid headrest and a bite bar are required to stabilize the head during OP measurement. Ultrasonic technique enables noninvasive and accurate in vivo measurement of corneal pulsation, which could be of interest for indirectly estimating intraocular pressure propagation and POBF component.

Phase dependencies between longitudinal corneal apex displacement and cardiovascular signals: is the ocular pulse influenced by the electrical activity of the heart?

The aim was to establish phase relationships between the principal harmonic, related to the heart rate, of synchronically registered longitudinal corneal apex displacement (LCAD), blood pulsation (BP) and electrical heart activity signals in a group of healthy subjects.

Glaucomatous and Age-Related Changes in Corneal Pulsation Shape. The Ocular Dicrotism

Purpose To ascertain whether the incidence of ocular dicrotic pulse (ODP) increases with age, it is more pronounced in glaucomatous than in normal eyes and whether it is related to cardiovascular activity. Methods 261 subjects aged 47 to 78 years were included in the study and classified into four groups: primary open angle glaucoma (POAG), primary angle-closure glaucoma (PACG), glaucoma suspects with glaucomatous optic disc appearance (GODA) and the controls (CG). Additionally, in each group, subjects with ODP were divided into two age subgroups around the median age. A non-contact ultrasonic method was used to measure corneal indentation pulse (CIP) synchronically with the acquisition of electrocardiography (ECG) and blood pulse signals. ODP was assessed from the acquired signals that were numerically processed in a custom written program. Results ODP incidence was about 78%, 66%, 66% and 84% for CG, GODA, POAG, and PACG group, respectively. With advancing age, the ODP incidence increased for all subjects (Δ = 12%), the highest being for the PACG and POAG groups (Δ = 30%). GODA group did not show an age-related increase in the incidence of ODP. Conclusions The ocular dicrotism, measured with non-contact ultrasonic method, was found to be a common phenomenon in elderly subjects. The increased ODP incidence in PACG and POAG group may correspond to either higher stiffness of glaucoma eyes, biochemical abnormalities in eye tissues, changes in ocular hemodynamics, may reflect the effect of medications or be a combination of all those factors. The results of GODA group suggest different mechanisms governing their ocular pulse that makes them less susceptible to generating ODP and having decreased predisposition to glaucoma.

Assessing Efficacy of Canaloplasty Using Continuous 24-Hour Monitoring of Ocular Dimensional Changes

PURPOSE We investigated whether 24-hour monitoring of corneoscleral limbus area (CSLA) with the Sensimed Triggerfish contact lens sensor (CLS) can be used clinically to assess midterm efficacy of canaloplasty and to assess the relationships of CSLA changes with the heart rhythm. METHODS Ten eyes of 10 patients, with POAG, which were qualified either to canaloplasty or canaloplasty and phacoemulsification, were included in this study. Eyes were washed out before the surgery and control visits were done at days 1, 7, and 3, 6, 12 months postoperatively, at which subjects were examined. We performed 24-hour monitoring of CSLA changes and Holter ECG at washout, and at 3- and 12-month postop visits. Raw CLS signals were processed to lead two parameters describing short-term 24-hour variability of CSLA (VAR, â). Heart activity parameters from CLS were compared to those acquired from Holter ECG. RESULTS Mean post washout IOP was 20.6 ± 4.7 and decreased to 14.2 ± 3.0 mm Hg 1 year after surgery (P < 0.01). A decreasing trend in VAR and â parameters were noted. Statistically significant differences were found between the washout and 3-month postop visit for VAR and â (P = 0.014 and P = 0.027, respectively) as well as between the washout and 12-month postop result for the parameter â (P = 0.031). No statistically significant differences were found between the 3- and 12-month postop results for both considered parameters. CONCLUSIONS Canaloplasty alone or combined with cataract surgery is a successful surgical method of lowering IOP in glaucoma patients. Canaloplasty decreases 24-hour CSLA fluctuation pattern measured with CLS.

Phase dependencies between longitudinal corneal apex displacement of human eye and cardiovascular system

Intraocular pressure (IOP) varies quasi-periodically due to blood pulsation in vessels inside the eye globe. This variations cause the eye deformations and displacements of the outer surface of the eye. The aim of this paper is to calculate the correlation between longitudinal corneal apex displacement and cardiovascular activity. Using ultrasound transducer at sampling frequency of 100Hz we have measured longitudinal corneal apex displacement (LCAD) of the left eye for 5 subjects. Synchronically we have registered ECG and blood pulsation signals at the same sampling frequency. Cross-correlation function was applied to investigate dependencies between these signals. To find time shift between LCAD and ECG or pulse, the time window of 3 seconds length have been chosen from all signals and had been shifting with the step of 0.01 seconds from 0 to 7s. For each shift the cross-correlation function and its extrema were calculated in the window area. We have obtained information about extrema position of cross-correlation function and its stability in time for particular subjects. The time shift between LCAD and ECG or pulse is individual feature of each subject. Such calculations may lead us to better understanding of pulse propagation in human eye and creation a non invasive method of eye hemodynamics and ocular diagnosis.

Assessing the Feasibility of the Use of Video Motion Magnification for Measuring Microdisplacements

There are different techniques for measuring microdisplacements. The purpose of this paper was to ascertain whether a method of video motion magnification (VMM) can be used for measuring such displacements. For this, standard video devices (a digital single-lens reflex camera and a webcam) were used to record subtle movements of an object, and the results of the VMM technique were contrasted with an air-coupled ultrasonic sensing method that could achieve submicrometer accuracy. The results of the VMM technique highly correlate with those achieved using the ultrasonic sensor, showing that the former can accurately measure displacements in the range from about 5 to 40

Wavelet Representation of the Corneal Pulse for Detecting Ocular Dicrotism

\mu \text{m}

Analysis of ocular hemodynamic using combined STdOCT and ultrasonic methods

from a distance of about 1 m. The temporal characteristics of the moving object were well preserved. The VMM technique is an alternative to other modalities for measuring microdisplacements and has the advantage of being noncontact, long-range, and relatively low-cost.

In-vivo corneal pulsation in relation to in-vivo intraocular pressure and corneal biomechanics assessed in-vitro. An animal pilot study

Purpose To develop a reliable and powerful method for detecting the ocular dicrotism from non-invasively acquired signals of corneal pulse without the knowledge of the underlying cardiopulmonary information present in signals of ocular blood pulse and the electrical heart activity. Methods Retrospective data from a study on glaucomatous and age-related changes in corneal pulsation [PLOS ONE 9(7),(2014):e102814] involving 261 subjects was used. Continuous wavelet representation of the signal derivative of the corneal pulse was considered with a complex Gaussian derivative function chosen as mother wavelet. Gray-level Co-occurrence Matrix has been applied to the image (heat-maps) of CWT to yield a set of parameters that can be used to devise the ocular dicrotic pulse detection schemes based on the Conditional Inference Tree and the Random Forest models. The detection scheme was first tested on synthetic signals resembling those of a dicrotic and a non-dicrotic ocular pulse before being used on all 261 real recordings. Results A detection scheme based on a single feature of the Continuous Wavelet Transform of the corneal pulse signal resulted in a low detection rate. Conglomeration of a set of features based on measures of texture (homogeneity, correlation, energy, and contrast) resulted in a high detection rate reaching 93%. Conclusion It is possible to reliably detect a dicrotic ocular pulse from the signals of corneal pulsation without the need of acquiring additional signals related to heart activity, which was the previous state-of-the-art. The proposed scheme can be applied to other non-stationary biomedical signals related to ocular dynamics.

Correlations between Longitudinal Corneal Apex Displacement, Head Movements and Pulsatile Blood Flow

The paper presents the way that colour can serve solving the problem of calibration points indexing in a camera geometrical calibration process. We propose a technique in which indexes of calibration points in a black-and-white chessboard are represented as sets of colour regions in the neighbourhood of calibration points. We provide some general rules for designing a colour calibration chessboard and provide a method of calibration image analysis. We show that this approach leads to obtaining better results than in the case of widely used methods employing information about already indexed points to compute indexes. We also report constraints concerning the technique. Nowadays we are witnessing an increasing need for camera geometrical calibration systems. They are vital for such applications as 3D modelling, 3D reconstruction, assembly control systems, etc. Wherever possible, calibration objects placed in the scene are used in a camera geometrical calibration process. This approach significantly increases accuracy of calibration results and makes the calibration data extraction process easier and universal. There are many geometrical camera calibration techniques for a known calibration scene [1]. A great number of them use as an input calibration points which are localised and indexed in the scene. In this paper we propose the technique of calibration points indexing which uses a colour chessboard. The presented technique was developed by solving problems we encountered during experiments with our earlier methods of camera calibration scene analysis [2]-[3]. In particular, the proposed technique increases the number of indexed points points in case of local lack of calibration points detection. At the beginning of the paper we present a way of designing a chessboard pattern. Then we describe a calibration point indexing method, and finally we show experimental results. A black-and-white chessboard is widely used in order to obtain sub-pixel accuracy of calibration points localisation [1]. Calibration points are defined as corners of chessboard squares. Assuming the availability of rough localisation of these points, the points can be indexed. Noting that differences in distances between neighbouring points in calibration scene images differ slightly, one of the local searching methods can be employed (e.g. [2]). Methods of this type search for a calibration point to be indexed, using a window of a certain size. The position of the window is determined by a vector representing the distance between two previously indexed points in the same row or column. However, experiments show that this approach has its disadvantages, as described below. * E-mail: A.Nowakowski@ire.pw.edu.pl Firstly, there is a danger of omitting some points during indexing in case of local lack of calibration points detection in a neighbourhood (e.g. caused by the presence of non-homogeneous light in the calibration scene). A particularly unfavourable situation is when the local lack of detection effects in the appearance of separated regions of detected calibration points. It is worth saying that such situations are likely to happen for calibration points situated near image borders. Such points are very important for the analysis of optical nonlinearities, and a lack of them can significantly influence the accuracy of distortion modelling. Secondly, such methods may give wrong results in the case of optical distortion with strong nonlinearities when getting information about the neighbouring index is not an easy task. Beside this, the methods are very sensitive to a single false localisation of a calibration point. Such a single false localisation can even result in false indexing of a big set of calibration points. To avoid the above-mentioned problems, we propose using a black-and-white chessboard which contains the coded index of a calibration point in the form of colour squares situated in the nearest neighbourhood of each point. The index of a certain calibration point is determined by colours of four nearest neighbouring squares (Fig.1). An order of squares in such foursome is important. Because the size of a colour square is determined only by the possibility of correct colour detection, the size of a colour square can be smaller than the size of a black or white square. The larger size of a black or white square is determined by the requirements of the exact localisation step which follows the indexing of calibration points [3]. In this step, edge information is extracted from a blackand-white chessboard. This edge information needs larger Artur Nowakowski, Wladyslaw Skarbek Institute of Radioelectronics, Warsaw University of Technology, Nowowiejska 15/19, 00-665 Warszawa, A.Nowakowski@ire.pw.edu.pl Received February 10, 2009; accepted March 27, 2009; published March 31, 2009 http://www.photonics.pl/PLP