Statistical distribution of wave-surface elevation for second-order random directional ocean waves in finite water depth

Based on the second-order random wave solutions of water wave equations in finite water depth, a statistical distribution of the wave-surface elevation is derived by using the characteristic function expansion method. It is found that the distribution, after normalization of the wave-surface elevation, depends only on two parameters. One parameter describes the small mean bias of the surface produced by the second-order wave-wave interactions. Another one is approximately proportional to the skewness of the distribution. Both of these two parameters can be determined by the water depth and the wave-number spectrum of ocean waves. As an illustrative example, we consider a fully developed wind-generated sea and the parameters are calculated for various wind speeds and water depths by using Donelan and Pierson spectrum. It is also found that, for deep water, the dimensionless distribution reduces to the third-order Gram-Charlier series obtained by Longuet-Higgins [J. Fluid Mech. 17 (1963) 459]. The newly proposed distribution is compared with the data of Bitner [Appl. Ocean Res. 2 (1980) 63], Gaussian distribution and the fourth-order Gram-Charlier series, and found our distribution gives a more reasonable fit to the data. (C) 2002 Elsevier Science B.V. All rights reserved.

Synthesis, Characterization and Applications of Modified TiO2 Systems

This thesis is divided in to 9 chapters and deals with the modification of TiO2 for various applications include photocatalysis, thermal reaction, photovoltaics and non-linear optics. Chapter 1 involves a brief introduction of the topic of study. An introduction to the applications of modified titania systems in various fields are discussed concisely. Scope and objectives of the present work are also discussed in this chapter. Chapter 2 explains the strategy adopted for the synthesis of metal, nonmetal co-doped TiO2 systems. Hydrothermal technique was employed for the preparation of the co-doped TiO2 system, where Ti[OCH(CH3)2]4, urea and metal nitrates were used as the sources for TiO2, N and metals respectively. In all the co-doped systems, urea to Ti[OCH(CH3)2]4 was taken in a 1:1 molar ratio and varied the concentration of metals. Five different co-doped catalytic systems and for each catalysts, three versions were prepared by varying the concentration of metals. A brief explanation of physico-chemical techniques used for the characterization of the material was also presented in this chapter. This includes X-ray Diffraction (XRD), Raman Spectroscopy, FTIR analysis, Thermo Gravimetric Analysis, Energy Dispersive X-ray Analysis (EDX), Scanning Electron Microscopy(SEM), UV-Visible Diffuse Reflectance Spectroscopy (UV-Vis DRS), Transmission Electron Microscopy (TEM), BET Surface Area Measurements and X-ray Photoelectron Spectroscopy (XPS). Chapter 3 contains the results and discussion of characterization techniques used for analyzing the prepared systems. Characterization is an inevitable part of materials research. Determination of physico-chemical properties of the prepared materials using suitable characterization techniques is very crucial to find its exact field of application. It is clear from the XRD pattern that photocatalytically active anatase phase dominates in the calcined samples with peaks at 2θ values around 25.4°, 38°, 48.1°, 55.2° and 62.7° corresponding to (101), (004), (200), (211) and (204) crystal planes (JCPDS 21-1272) respectively. But in the case of Pr-N-Ti sample, a new peak was observed at 2θ = 30.8° corresponding to the (121) plane of the polymorph brookite. There are no visible peaks corresponding to dopants, which may be due to their low concentration or it is an indication of the better dispersion of impurities in the TiO2. Crystallite size of the sample was calculated from Scherrer equation byusing full width at half maximum (FWHM) of the (101) peak of the anatase phase. Crystallite size of all the co-doped TiO2 was found to be lower than that of bare TiO2 which indicates that the doping of metal ions having higher ionic radius into the lattice of TiO2 causes some lattice distortion which suppress the growth of TiO2 nanoparticles. The structural identity of the prepared system obtained from XRD pattern is further confirmed by Raman spectra measurements. Anatase has six Raman active modes. Band gap of the co-doped system was calculated using Kubelka-Munk equation and that was found to be lower than pure TiO2. Stability of the prepared systems was understood from thermo gravimetric analysis. FT-IR was performed to understand the functional groups as well as to study the surface changes occurred during modification. EDX was used to determine the impurities present in the system. The EDX spectra of all the co-doped samples show signals directly related to the dopants. Spectra of all the co-doped systems contain O and Ti as the main components with low concentrations of doped elements. Morphologies of the prepared systems were obtained from SEM and TEM analysis. Average particle size of the systems was drawn from histogram data. Electronic structures of the samples were identified perfectly from XPS measurements. Chapter 4 describes the photocatalytic degradation of herbicides Atrazine and Metolachlor using metal, non-metal co-doped titania systems. The percentage of degradation was analyzed by HPLC technique. Parameters such as effect of different catalysts, effect of time, effect of catalysts amount and reusability studies were discussed. Chapter 5 deals with the photo-oxidation of some anthracene derivatives by co-doped catalytic systems. These anthracene derivatives come underthe category of polycyclic aromatic hydrocarbons (PAH). Due to the presence of stable benzene rings, most of the PAH show strong inhibition towards biological degradation and the common methods employed for their removal. According to environmental protection agency, most of the PAH are highly toxic in nature. TiO2 photochemistry has been extensively investigated as a method for the catalytic conversion of such organic compounds, highlighting the potential of thereof in the green chemistry. There are actually two methods for the removal of pollutants from the ecosystem. Complete mineralization is the one way to remove pollutants. Conversion of toxic compounds to another compound having toxicity less than the initial starting compound is the second way. Here in this chapter, we are concentrating on the second aspect. The catalysts used were Gd(1wt%)-N-Ti, Pd(1wt%)-N-Ti and Ag(1wt%)-N-Ti. Here we were very successfully converted all the PAH to anthraquinone, a compound having diverse applications in industrial as well as medical fields. Substitution of 10th position of desired PAH by phenyl ring reduces the feasibility of photo reaction and produced 9-hydroxy 9-phenyl anthrone (9H9PA) as an intermediate species. The products were separated and purified by column chromatography using 70:30 hexane/DCM mixtures as the mobile phase and the resultant products were characterized thoroughly by 1H NMR, IR spectroscopy and GCMS analysis. Chapter 6 elucidates the heterogeneous Suzuki coupling reaction by Cu/Pd bimetallic supported on TiO2. Sol-Gel followed by impregnation method was adopted for the synthesis of Cu/Pd-TiO2. The prepared system was characterized by XRD, TG-DTG, SEM, EDX, BET Surface area and XPS. The product was separated and purified by column chromatography using hexane as the mobile phase. Maximum isolated yield of biphenyl of around72% was obtained in DMF using Cu(2wt%)-Pd(4wt%)-Ti as the catalyst. In this reaction, effective solvent, base and catalyst were found to be DMF, K2CO3 and Cu(2wt%)-Pd(4wt%)-Ti respectively. Chapter 7 gives an idea about the photovoltaic (PV) applications of TiO2 based thin films. Due to energy crisis, the whole world is looking for a new sustainable energy source. Harnessing solar energy is one of the most promising ways to tackle this issue. The present dominant photovoltaic (PV) technologies are based on inorganic materials. But the high material, low power conversion efficiency and manufacturing cost limits its popularization. A lot of research has been conducted towards the development of low-cost PV technologies, of which organic photovoltaic (OPV) devices are one of the promising. Here two TiO2 thin films having different thickness were prepared by spin coating technique. The prepared films were characterized by XRD, AFM and conductivity measurements. The thickness of the films was measured by Stylus Profiler. This chapter mainly concentrated on the fabrication of an inverted hetero junction solar cell using conducting polymer MEH-PPV as photo active layer. Here TiO2 was used as the electron transport layer. Thin films of MEH-PPV were also prepared using spin coating technique. Two fullerene derivatives such as PCBM and ICBA were introduced into the device in order to improve the power conversion efficiency. Effective charge transfer between the conducting polymer and ICBA were understood from fluorescence quenching studies. The fabricated Inverted hetero junction exhibited maximum power conversion efficiency of 0.22% with ICBA as the acceptor molecule. Chapter 8 narrates the third order order nonlinear optical properties of bare and noble metal modified TiO2 thin films. Thin films were fabricatedby spray pyrolysis technique. Sol-Gel derived Ti[OCH(CH3)2]4 in CH3CH2OH/CH3COOH was used as the precursor for TiO2. The precursors used for Au, Ag and Pd were the aqueous solutions of HAuCl4, AgNO3 and Pd(NO3)2 respectively. The prepared films were characterized by XRD, SEM and EDX. The nonlinear optical properties of the prepared materials were investigated by Z-Scan technique comprising of Nd-YAG laser (532 nm,7 ns and10 Hz). The non-linear coefficients were obtained by fitting the experimental Z-Scan plot with the theoretical plots. Nonlinear absorption is a phenomenon defined as a nonlinear change (increase or decrease) in absorption with increasing of intensity. This can be mainly divided into two types: saturable absorption (SA) and reverse saturable absorption (RSA). Depending on the pump intensity and on the absorption cross- section at the excitation wavelength, most molecules show non- linear absorption. With increasing intensity, if the excited states show saturation owing to their long lifetimes, the transmission will show SA characteristics. Here absorption decreases with increase of intensity. If, however, the excited state has strong absorption compared with that of the ground state, the transmission will show RSA characteristics. Here in our work most of the materials show SA behavior and some materials exhibited RSA behavior. Both these properties purely depend on the nature of the materials and alignment of energy states within them. Both these SA and RSA have got immense applications in electronic devices. The important results obtained from various studies are presented in chapter 9.

Blind equalization of nonminimum phase channels: higher order cumulant based algorithm

Higher order cumulant analysis is applied to the blind equalization of linear time-invariant (LTI) nonminimum-phase channels. The channel model is moving-average based. To identify the moving average parameters of channels, a higher-order cumulant fitting approach is adopted in which a novel relay algorithm is proposed to obtain the global solution. In addition, the technique incorporates model order determination. The transmitted data are considered as independently identically distributed random variables over some discrete finite set (e.g., set {±1, ±3}). A transformation scheme is suggested so that third-order cumulant analysis can be applied to this type of data. Simulation examples verify the feasibility and potential of the algorithm. Performance is compared with that of the noncumulant-based Sato scheme in terms of the steady state MSE and convergence rate.

Efeitos de variáveis abióticas na composição química de Gelidium Crinale (Gelidiaceae, Rhodophyta) em cultivo unialgal

Os efeitos individuais e interativos dos parâmetros ambientais físicos e químicos, como temperatura, intensidade luminosa, salinidade e concentração de fósforo inorgânico dissolvido na água do mar, na produção de proteínas, carboidratos, acúmulo de fósforo tecidual e taxa de absorção do fósforo inorgânico disponível no meio de cultura em Gelidium crinale (Turner) Lamouroux, foram investigados durante um período de sete dias de cultivo laboratorial, em condições controladas. A ação dos parâmetros abióticos foi analisada de três maneiras diferentes. A primeira avaliação integrou a ação de temperatura, intensidade luminosa e fósforo inorgânico dissolvido, mantendo-se fixa a salinidade em 25 ups, onde se constatou que em todos os componentes químicos algais ocorreram interações de terceira ordem. O incremento de 2,28 a 2,67 % nos teores de proteínas foram obtidos à temperatura de 25 °C e 12 μmol m-2 s-1 de intensidade luminosa, diminuindo com a elevação da intensidade luminosa para 40 μmol m-2 s-1. Para carboidratos, ocorreram interações significativas entre os três parâmetros, com um aumento de 6,85 % sendo registrado a 25 °C de temperatura, 24 μmol m-2 s-1 de intensidade luminosa e 10,0 μM de fósforo inorgânico. O aumento máximo na taxa de fósforo tecidual (0,56 %) ocorreu em talos cultivados nas menores temperatura e intensidade luminosa e na maior concentração de fósforo inorgânico dissolvido. Com relação à intensidade luminosa, foi observada uma correlação negativa entre proteínas e carboidratos. A segunda avaliação estabeleceu a ação independente e sinérgica de temperatura, salinidade e fósforo inorgânico disponível no meio de cultivo, fixando-se a intensidade luminosa em 24 μmol m-2s-1. A maior produção de proteínas ocorreu em cultivos onde a temperatura foi de 25 °C, com uma concentração de 5,0 e 10,0 μM de fósforo inorgânico dissolvido e salinidade entre 15 e 20 ups, cujos valores médios do incremento variaram entre 2,62 a 2,83 % peso seco de alga, resultando em uma interação de terceira ordem altamente significativa. Para carboidratos a elevação de 6,85 % em sua concentração está associada à maior temperatura (25 °C), maior salinidade (25 ups) e maior quantidade de fósforo inorgânico disponível no meio de cultivo (10,0 μM). Contudo, não foi observada uma interação de terceira ordem através da análise estatística. Para esta biomolécula observaram-se interações de segunda ordem altamente significativa (P < 0,005) entre temperatura e diferentes concentrações de fósforo inorgânico e entre temperatura e salinidade (P < 0,000). O acúmulo de fósforo nos talos da alga foi menor durante os cultivos em que a salinidade foi de 25 ups,nas temperaturas de 20 e 25 °C e concentração de fósforo disponível de 2,5 μM, com percentuais entre 0,08 a 0,11 % em peso de cinzas. O maior incremento ocorreu na menor temperatura, associada à baixa salinidade e alta concentração de fósforo inorgânico no meio. O coeficiente de correlação de Pearson revelou correlações positivas, altamente significativas (P < 0,001) entre teor de proteína, temperatura e disponibilidade de fósforo inorgânico no meio de cultivo. Para carboidratos, as correlações foram positivas com os três parâmetros abióticos. Para fósforo tecidual somente com o fósforo inorgânico disponível no cultivo foi que ocorreu uma relação positiva; com os outros dois parâmetros esta correlação foi negativa. Entre os componentes químicos encontrados nas algas, proteínas e carboidratos apresentaram uma relação positiva, porém fósforo tecidual apresentou uma correlação negativa com ambos, embora com proteínas esta relação não tenha sido significativa. A terceira avaliação estudou a ação individual e o sinergismo entre os parâmetros ambientais, temperatura, intensidade luminosa e salinidade, a uma concentração fixa de fósforo inorgânico disponível no meio de cultivo (10,0 μM), sobre a composição química, bem como na taxa de absorção de fósforo inorgânico disponível. Observou-se a ocorrência de interações de terceira ordem em todos as variáveis estudadas. O teor de proteínas apresentou um aumento de 3,72 % durante o período de cultivo, passando de 20,63 % antes do cultivo, para 24,35 % após o término do experimento, principalmente nas condições de 25 °C de temperatura, 12 μmol m-2s-1 de intensidade luminosa e 15 ups de salinidade. Para carboidratos, nas condições de baixa intensidade luminosa (12 μmol m-2s- 1), a uma temperatura de 20 °C e salinidades de 10 e 15 ups, foram registrados valores inferiores à amostra controle, caracterizando um consumo desta biomolécula por parte das algas. Nestas mesmas condições ambientais, foram registrados os maiores teores de fósforo tecidual, variando entre 0,86 a 1,09 % do peso das cinzas. As maiores taxas de absorção do fósforo do meio ocorreram na salinidade de 25 ups e 25 °C de temperatura, diminuindo da intensidade luminosa de 12 μmol m-2s-1 para 40 μmol m-2s-1. As maiores concentrações de fósforo inorgânico residual na água do meio de cultivo ocorreram nas salinidades de 10 e 15 ups, em todas as intensidade luminosas e temperaturas estudadas. Através do coeficiente de correlação de Pearson, observou-se que os teores de proteínas apresentaram uma forte correlação negativa com a intensidade luminosa e positiva com a temperatura e salinidade, embora com esta última não tenha sido significativa. Para carboidratos, as correlações com os parâmetros abióticos foram todas positivas. Correlações negativa e positiva, não significativas, foram observadas entre esta biomolécula e o teor de proteínas e a taxa de absorção de fósforo disponível no meio, respectivamente. Por outro lado, com fósforo tecidual, ocorreu uma correlação negativa, altamente significativa. Este estudo mostra o estado fisiológico de Gelidium crinale e contribui para o estabelecimento das melhores condições de cultivo para produção de proteína, carboidrato e fósforo tecidual e indicação do uso racional de nutrientes, fornecendo informações para a otimização de processos de maricultura, tanto em termos de cultivo bem sucedido de algas, quanto de redução no impacto sobre o ambiente.

Estimation of genetic parameters for partial egg production periods by means of random regression models

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Estimates of genetic parameters for scrotal circumference using random regression models in Nelore cattle

A total of 15,901 scrotal circumference (SC) records from 5300 Nelore bulls, ranging from 229 to 560 days of age, were used with the objective of estimating (co)variance functions for SC, using random regression models. Models included the fixed effects of contemporary group and age of dam at calving as covariable (linear and quadratic effects). To model the population mean trend, a third order Legendre polynomial on animal age was utilized. The direct additive genetic and animal permanent environmental random effects were modeled by Legendre polynomials on animal age, with orders of fit ranging from 1 to 5. Residual variances were modeled considering 1 (homogeneity of variance) or 4 age classes. Results obtained with the random regression models were compared to multi-trait analysis. (Co)variance estimates using multi-trait and random regression models were similar. The model considering a third- and fifth-order Legendre polynomials for additive genetic and animal permanent environmental effects, respectively, was the most adequate to model changes in variance of SC with age. Heritability estimates for SC ranged from 0.24 (229 days of age) to 0.47 (300 days of age), remained almost constant until 500 days of age (0.52), decreasing thereafter (0.44). In general, the genetic correlations between measures of scrotal circumference obtained from 229 to 560 days of age decreased with increasing distance between ages. For genetic evaluation scrotal circumference could be measured between 400 and 500 days of age. (C) 2010 Elsevier B.V. All rights reserved.

Genetic parameters for test-day yield of milk, fat and protein in buffaloes estimated by random regression models

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Comparison of random regression models to estimate genetic parameters for milk production in Guzerat (Bos indicus) cows

Random regression models have been widely used to estimate genetic parameters that influence milk production in Bos taurus breeds, and more recently in B. indicus breeds. With the aim of finding appropriate random regression model to analyze milk yield, different parametric functions were compared, applied to 20,524 test-day milk yield records of 2816 first-lactation Guzerat (B. indicus) cows in Brazilian herds. The records were analyzed by random regression models whose random effects were additive genetic, permanent environmental and residual, and whose fixed effects were contemporary group, the covariable cow age at calving (linear and quadratic effects), and the herd lactation curve. The additive genetic and permanent environmental effects were modeled by the Wilmink function, a modified Wilmink function (with the second term divided by 100), a function that combined third-order Legendre polynomials with the last term of the Wilmink function, and the Ali and Schaeffer function. The residual variances were modeled by means of 1, 4, 6, or 10 heterogeneous classes, with the exception of the last term of the Wilmink function, for which there were 1, from 0.20 to 0.33. Genetic correlations between adjacent records were high values (0.83-0.99), but they declined when the interval between the test-day records increased, and were negative between the first and last records. The model employing the Ali and Schaeffer function with six residual variance classes was the most suitable for fitting the data. © FUNPEC-RP.

Genetic parameters for body weight in buffaloes (Bubalus bubalis) in Colombia using random regression models

The objective of this research was to estimate (co) variance functions and genetic parameters for body weight in Colombian buffalo populations using random regression models with Legendre polynomials. Data consisted of 34,738 weight records from birth to 900 days of age from 7815 buffaloes. Fixed effects in the model were contemporary group and parity order of the mother. Random effects were direct and maternal additive genetic, as well as animal and maternal permanent environmental effects. A cubic orthogonal Legendre polynomial was used to model the mean curve of the population. Eleven models with first to sixth order polynomials were used to describe additive genetic direct and maternal effects, and animal and maternal permanent environmental effects. The residual was modeled considering five variance classes. The best model included fourth and sixth order polynomials for direct additive genetic and animal permanent environmental effects, respectively, and third-order polynomials for maternal genetic and maternal permanent environmental effects. The direct heritability increased from birth until 120 days of age (0.32 +/- 0.05), decreasing thereafter until one year of age (0.18 +/- 0.04) and increased again, reaching 0.39 +/- 0.09, at the end of the evaluated period. The highest maternal heritability estimates (0.11 +/- 0.05), were obtained for weights around weaning age (weaning age range is between 8 and 9.5 months). Maternal genetic and maternal permanent environmental variances increased from birth until about one year of age, decreasing at later ages. Direct genetic correlations ranged from moderate (0.60 +/- 0.060) to high (0.99 +/- 0.001), maternal genetic correlations showed a similar range (0.41 +/- 0.401 and 0.99 +/- 0.003), and all of them decreased as time between weighings increased. Direct genetic correlations suggested that selecting buffalos for heavier weights at any age would increase weights from birth through 900 days of age. However, higher heritabilities for direct genetic weights effects after 600 days of age suggested that selection for these effects would be more effective if done during this age period. A greater response to selection for maternal ability would be expected if selection used maternal genetic predictions for weights near weaning. (C) 2013 Elsevier B.V. All rights reserved.

Application of the Z-scan technique to determine the optical Kerr coefficient and two-photon absorption coefficient of magnetite nanoparticles colloidal suspension

We investigate the occurrence of the optical Kerr effect and two-photon absorption when an oil-based magnetic Fe3O4 nanoparticles colloidal suspension is illuminated with high intensity femtosecond laser pulses. The frequency of the pulses is controlled and the Z-scan technique is employed in our measurements of the nonlinear optical Kerr coefficient (n(2)) and two-photon absorption coefficient (beta). From these values it was possible to calculate the real and imaginary parts of the third-order susceptibility. We observed that increasing the pulse frequency, additional physical processes take place, increasing artificially the absolute values of n(2) and beta. The experimental conditions are discussed to assure the obtention of reliable values of these nonlinear optical parameters, which may be useful in all-optical switching and optical power limiting applications. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4723829]

Simulación de componentes no estacionarias de viento en orografías complejas

Realistic operation of helicopter flight simulators in complex topographies (such as urban environments) requires appropriate prediction of the incoming wind, and this prediction should be made in real time. Unfortunately, the wind topology around complex topographies shows time-dependent, fully nonlinear, turbulent patterns (i.e., wakes) whose simulation cannot be made using computationally inexpensive tools based on corrected potential approximations. Instead, the full Navier-Stokes plus some kind of turbulent modeling is necessary, which is quite computationally expensive. The complete unsteady flow depends on two parameters, namely the velocity and orientation of the free stream flow. The aim of this MSc thesis is to develop a methodology for the real time simulation of these complex flows. For simplicity, the flow around a single building (20 mx20 m cross section and 100 m height) is considered, with free stream velocity in the range 5-25 m/s. Because of the square cross section, the problem shows two reflection symmetries, which allows for restricting the orientations to the range 0° < a. < 45°. The methodology includes an offline preprocess and the online operation. The preprocess consists in three steps: An appropriate, unstructured mesh is selected in which the flow is sim¬ulated using OpenFOAM, and this is done for 33 combinations of 3 free stream intensities and 11 orientations. For each of these, the simulation proceeds for a sufficiently large time as to eliminate transients. This step is quite computationally expensive. Each flow field is post-processed using a combination of proper orthogonal decomposition, fast Fourier transform, and a convenient optimization tool, which identifies the relevant frequencies (namely, both the basic frequencies and their harmonics) and modes in the computational mesh. This combination includes several new ingredients to filter errors out and identify the relevant spatio-temporal patterns. Note that, in principle, the basic frequencies depend on both the intensity and the orientation of the free stream flow. The outcome of this step is a set of modes (vectors containing the three velocity components at all mesh points) for the various Fourier components, intensities, and orientations, which can be organized as a third order tensor. This step is fairly computationally inexpensive. The above mentioned tensor is treated using a combination of truncated high order singular value, decomposition and appropriate one-dimensional interpolation (as in Lorente, Velazquez, Vega, J. Aircraft, 45 (2008) 1779-1788). The outcome is a tensor representation of both the relevant fre¬quencies and the associated Fourier modes for a given pair of values of the free stream flow intensity and orientation. This step is fairly compu¬tationally inexpensive. The online, operation requires just reconstructing the time-dependent flow field from its Fourier representation, which is extremely computationally inex¬pensive. The whole method is quite robust.

Application of higher order statistics/spectra in biomedical signals : a review

For many decades correlation and power spectrum have been primary tools for digital signal processing applications in the biomedical area. The information contained in the power spectrum is essentially that of the autocorrelation sequence; which is sufficient for complete statistical descriptions of Gaussian signals of known means. However, there are practical situations where one needs to look beyond autocorrelation of a signal to extract information regarding deviation from Gaussianity and the presence of phase relations. Higher order spectra, also known as polyspectra, are spectral representations of higher order statistics, i.e. moments and cumulants of third order and beyond. HOS (higher order statistics or higher order spectra) can detect deviations from linearity, stationarity or Gaussianity in the signal. Most of the biomedical signals are non-linear, non-stationary and non-Gaussian in nature and therefore it can be more advantageous to analyze them with HOS compared to the use of second order correlations and power spectra. In this paper we have discussed the application of HOS for different bio-signals. HOS methods of analysis are explained using a typical heart rate variability (HRV) signal and applications to other signals are reviewed.

Ocular characteristics of anisometropia

Animal models of refractive error development have demonstrated that visual experience influences ocular growth. In a variety of species, axial anisometropia (i.e. a difference in the length of the two eyes) can be induced through unilateral occlusion, image degradation or optical manipulation. In humans, anisometropia may occur in isolation or in association with amblyopia, strabismus or unilateral pathology. Non-amblyopic myopic anisometropia represents an interesting anomaly of ocular growth, since the two eyes within one visual system have grown to different endpoints. These experiments have investigated a range of biometric, optical and mechanical properties of anisometropic eyes (with and without amblyopia) with the aim of improving our current understanding of asymmetric refractive error development. In the first experiment, the interocular symmetry in 34 non-amblyopic myopic anisometropes (31 Asian, 3 Caucasian) was examined during relaxed accommodation. A high degree of symmetry was observed between the fellow eyes for a range of optical, biometric and biomechanical measurements. When the magnitude of anisometropia exceeded 1.75 D, the more myopic eye was almost always the sighting dominant eye. Further analysis of the optical and biometric properties of the dominant and non-dominant eyes was conducted to determine any related factors but no significant interocular differences were observed with respect to best-corrected visual acuity, corneal or total ocular aberrations during relaxed accommodation. Given the high degree of symmetry observed between the fellow eyes during distance viewing in the first experiment and the strong association previously reported between near work and myopia development, the aim of the second experiment was to investigate the symmetry between the fellow eyes of the same 34 myopic anisometropes following a period of near work. Symmetrical changes in corneal and total ocular aberrations were observed following a short reading task (10 minutes, 2.5 D accommodation demand) which was attributed to the high degree of interocular symmetry for measures of anterior eye morphology, and corneal biomechanics. These changes were related to eyelid shape and position during downward gaze, but gave no clear indication of factors associated with near work that might cause asymmetric eye growth within an individual. Since the influence of near work on eye growth is likely to be most obvious during, rather than following near tasks, in the third experiment the interocular symmetry of the optical and biometric changes was examined during accommodation for 11 myopic anisometropes. The changes in anterior eye biometrics associated with accommodation were again similar between the eyes, resulting in symmetrical changes in the optical characteristics. However, the more myopic eyes exhibited slightly greater amounts of axial elongation during accommodation which may be related to the force exerted by the ciliary muscle. This small asymmetry in axial elongation we observed between the eyes may be due to interocular differences in posterior eye structure, given that the accommodative response was equal between eyes. Using ocular coherence tomography a reduced average choroidal thickness was observed in the more myopic eyes compared to the less myopic eyes of these subjects. The interocular difference in choroidal thickness was correlated with the magnitude of spherical equivalent and axial anisometropia. The symmetry in optics and biometrics between fellow eyes which have undergone significantly different visual development (i.e. anisometropic subjects with amblyopia) is also of interest with respect to refractive error development. In the final experiment the influence of altered visual experience upon corneal and ocular higher-order aberrations was investigated in 21 amblyopic subjects (8 refractive, 11 strabismic and 2 form deprivation). Significant differences in aberrations were observed between the fellow eyes, which varied according to the type of amblyopia. Refractive amblyopes displayed significantly higher levels of 4th order corneal aberrations (spherical aberration and secondary astigmatism) in the amblyopic eye compared to the fellow non-amblyopic eye. Strabismic amblyopes exhibited significantly higher levels of trefoil, a third order aberration, in the amblyopic eye for both corneal and total ocular aberrations. The results of this experiment suggest that asymmetric visual experience during development is associated with asymmetries in higher-order aberrations, proportional to the magnitude of anisometropia and dependent upon the amblyogenic factor. This suggests a direct link between the development of higher-order optical characteristics of the human eye and visual feedback. The results from these experiments have shown that a high degree of symmetry exists between the fellow eyes of non-amblyopic myopic anisometropes for a range of biomechanical, biometric and optical parameters for different levels of accommodation and following near work. While a single specific optical or biomechanical factor that is consistently associated with asymmetric refractive error development has not been identified, the findings from these studies suggest that further research into the association between ocular dominance, choroidal thickness and higher-order aberrations with anisometropia may improve our understanding of refractive error development.

Application of higher-order spectra for automated grading of diabetic maculopathy

Diabetic macular edema (DME) is one of the most common causes of visual loss among diabetes mellitus patients. Early detection and successive treatment may improve the visual acuity. DME is mainly graded into non-clinically significant macular edema (NCSME) and clinically significant macular edema according to the location of hard exudates in the macula region. DME can be identified by manual examination of fundus images. It is laborious and resource intensive. Hence, in this work, automated grading of DME is proposed using higher-order spectra (HOS) of Radon transform projections of the fundus images. We have used third-order cumulants and bispectrum magnitude, in this work, as features, and compared their performance. They can capture subtle changes in the fundus image. Spectral regression discriminant analysis (SRDA) reduces feature dimension, and minimum redundancy maximum relevance method is used to rank the significant SRDA components. Ranked features are fed to various supervised classifiers, viz. Naive Bayes, AdaBoost and support vector machine, to discriminate No DME, NCSME and clinically significant macular edema classes. The performance of our system is evaluated using the publicly available MESSIDOR dataset (300 images) and also verified with a local dataset (300 images). Our results show that HOS cumulants and bispectrum magnitude obtained an average accuracy of 95.56 and 94.39 % for MESSIDOR dataset and 95.93 and 93.33 % for local dataset, respectively.

Attitude estimation using moment invariants

The problem of estimating the three-dimensional rotational parameters of a rigid body from its monocular image data has been considered using the method of moment invariants. Second- and third-order moment invariants are used to construct the feature vector for the scale and orientation independent identification of the camera view axis direction in the body-fixed reference frame. The camera rotation angle about the view axis is derived from second-order central moments. The relative attitude of the rigid body is then expressed in terms of quaternion parameters to model the outputs of a video sensor in attitude control simulations. Experimental results and simulation outputs are presented using the mathematical model of a spacecraft.