Edges and bars: where do people see features in 1-D images?

There have been two main approaches to feature detection in human and computer vision - based either on the luminance distribution and its spatial derivatives, or on the spatial distribution of local contrast energy. Thus, bars and edges might arise from peaks of luminance and luminance gradient respectively, or bars and edges might be found at peaks of local energy, where local phases are aligned across spatial frequency. This basic issue of definition is important because it guides more detailed models and interpretations of early vision. Which approach better describes the perceived positions of features in images? We used the class of 1-D images defined by Morrone and Burr in which the amplitude spectrum is that of a (partially blurred) square-wave and all Fourier components have a common phase. Observers used a cursor to mark where bars and edges were seen for different test phases (Experiment 1) or judged the spatial alignment of contours that had different phases (e.g. 0 degrees and 45 degrees ; Experiment 2). The feature positions defined by both tasks shifted systematically to the left or right according to the sign of the phase offset, increasing with the degree of blur. These shifts were well predicted by the location of luminance peaks (bars) and gradient peaks (edges), but not by energy peaks which (by design) predicted no shift at all. These results encourage models based on a Gaussian-derivative framework, but do not support the idea that human vision uses points of phase alignment to find local, first-order features. Nevertheless, we argue that both approaches are presently incomplete and a better understanding of early vision may combine insights from both. (C)2004 Elsevier Ltd. All rights reserved.

Depth propagation and surface construction in 3-D vision

In stereo vision, regions with ambiguous or unspecified disparity can acquire perceived depth from unambiguous regions. This has been called stereo capture, depth interpolation or surface completion. We studied some striking induced depth effects suggesting that depth interpolation and surface completion are distinct stages of visual processing. An inducing texture (2-D Gaussian noise) had sinusoidal modulation of disparity, creating a smooth horizontal corrugation. The central region of this surface was replaced by various test patterns whose perceived corrugation was measured. When the test image was horizontal 1-D noise, shown to one eye or to both eyes without disparity, it appeared corrugated in much the same way as the disparity-modulated (DM) flanking regions. But when the test image was 2-D noise, or vertical 1-D noise, little or no depth was induced. This suggests that horizontal orientation was a key factor. For a horizontal sine-wave luminance grating, strong depth was induced, but for a square-wave grating, depth was induced only when its edges were aligned with the peaks and troughs of the DM flanking surface. These and related results suggest that disparity (or local depth) propagates along horizontal 1-D features, and then a 3-D surface is constructed from the depth samples acquired. The shape of the constructed surface can be different from the inducer, and so surface construction appears to operate on the results of a more local depth propagation process.

Eigen-Adaptation and Distributed Representation of 2-D Phase, Energy, Scale and Orientation in Spatial Vision

Distributed representations (DR) of cortical channels are pervasive in models of spatio-temporal vision. A central idea that underpins current innovations of DR stems from the extension of 1-D phase into 2-D images. Neurophysiological evidence, however, provides tenuous support for a quadrature representation in the visual cortex, since even phase visual units are associated with broader orientation tuning than odd phase visual units (J.Neurophys.,88,455–463, 2002). We demonstrate that the application of the steering theorems to a 2-D definition of phase afforded by the Riesz Transform (IEEE Trans. Sig. Proc., 49, 3136–3144), to include a Scale Transform, allows one to smoothly interpolate across 2-D phase and pass from circularly symmetric to orientation tuned visual units, and from more narrowly tuned odd symmetric units to even ones. Steering across 2-D phase and scale can be orthogonalized via a linearizing transformation. Using the tiltafter effect as an example, we argue that effects of visual adaptation can be better explained by via an orthogonal rather than channel specific representation of visual units. This is because of the ability to explicitly account for isotropic and cross-orientation adaptation effect from the orthogonal representation from which both direct and indirect tilt after-effects can be explained.

Phakometric measurement of ocular surface radii of curvature, axial separations and alignment in relaxed and accommodated human eyes

Measurements (autokeratometry, A-scan ultrasonography and video ophthalmophakometry) of ocular surface radii, axial separations and alignment were made in the horizontal meridian of nine emmetropes (aged 20-38 years) with relaxed (cycloplegia) and active accommodation (mean ± 95% confidence interval: 3.7 ± 1.1 D). The anterior chamber depth (-1.5 ± 0.3 D) and both crystalline lens surfaces (front 3.1 ± 0.8 D; rear 2.1 ± 0.6 D) contributed to dioptric vergence changes that accompany accommodation. Accommodation did not alter ocular surface alignment. Ocular misalignment in relaxed eyes is mainly because of eye rotation (5.7 ± 1.6° temporally) with small amounts of lens tilt (0.2 ± 0.8° temporally) and decentration (0.1 ± 0.1 mm nasally) but these results must be viewed with caution as we did not account for corneal asymmetry. Comparison of calculated and empirically derived coefficients (upon which ocular surface alignment calculations depend) revealed that negligible inherent errors arose from neglect of ocular surface asphericity, lens gradient refractive index properties, surface astigmatism, effects of pupil size and centration, assumed eye rotation axis position and use of linear equations for analysing Purkinje image shifts. © 2004 The College of Optometrists.

Bars and edges: what makes a feature for human vision?

There have been two main approaches to feature detection in human and computer vision - luminance-based and energy-based. Bars and edges might arise from peaks of luminance and luminance gradient respectively, or bars and edges might be found at peaks of local energy, where local phases are aligned across spatial frequency. This basic issue of definition is important because it guides more detailed models and interpretations of early vision. Which approach better describes the perceived positions of elements in a 3-element contour-alignment task? We used the class of 1-D images defined by Morrone and Burr in which the amplitude spectrum is that of a (partially blurred) square wave and Fourier components in a given image have a common phase. Observers judged whether the centre element (eg ±458 phase) was to the left or right of the flanking pair (eg 0º phase). Lateral offset of the centre element was varied to find the point of subjective alignment from the fitted psychometric function. This point shifted systematically to the left or right according to the sign of the centre phase, increasing with the degree of blur. These shifts were well predicted by the location of luminance peaks and other derivative-based features, but not by energy peaks which (by design) predicted no shift at all. These results on contour alignment agree well with earlier ones from a more explicit feature-marking task, and strongly suggest that human vision does not use local energy peaks to locate basic first-order features. [Supported by the Wellcome Trust (ref: 056093)]

The primal sketch revisited: locating and representing edges in human vision via Gaussian-derivative filtering

Marr's work offered guidelines on how to investigate vision (the theory - algorithm - implementation distinction), as well as specific proposals on how vision is done. Many of the latter have inevitably been superseded, but the approach was inspirational and remains so. Marr saw the computational study of vision as tightly linked to psychophysics and neurophysiology, but the last twenty years have seen some weakening of that integration. Because feature detection is a key stage in early human vision, we have returned to basic questions about representation of edges at coarse and fine scales. We describe an explicit model in the spirit of the primal sketch, but tightly constrained by psychophysical data. Results from two tasks (location-marking and blur-matching) point strongly to the central role played by second-derivative operators, as proposed by Marr and Hildreth. Edge location and blur are evaluated by finding the location and scale of the Gaussian-derivative `template' that best matches the second-derivative profile (`signature') of the edge. The system is scale-invariant, and accurately predicts blur-matching data for a wide variety of 1-D and 2-D images. By finding the best-fitting scale, it implements a form of local scale selection and circumvents the knotty problem of integrating filter outputs across scales. [Supported by BBSRC and the Wellcome Trust]

Cascaded Bayesian processes: an account of bias in orientation perception

Following adaptation to an oriented (1-d) signal in central vision, the orientation of subsequently viewed test signals may appear repelled away from or attracted towards the adapting orientation. Small angular differences between the adaptor and test yield 'repulsive' shifts, while large angular differences yield 'attractive' shifts. In peripheral vision, however, both small and large angular differences yield repulsive shifts. To account for these tilt after-effects (TAEs), a cascaded model of orientation estimation that is optimized using hierarchical Bayesian methods is proposed. The model accounts for orientation bias through adaptation-induced losses in information that arise because of signal uncertainties and neural constraints placed upon the propagation of visual information. Repulsive (direct) TAEs arise at early stages of visual processing from adaptation of orientation-selective units with peak sensitivity at the orientation of the adaptor (theta). Attractive (indirect) TAEs result from adaptation of second-stage units with peak sensitivity at theta and theta+90 degrees , which arise from an efficient stage of linear compression that pools across the responses of the first-stage orientation-selective units. A spatial orientation vector is estimated from the transformed oriented unit responses. The change from attractive to repulsive TAEs in peripheral vision can be explained by the differing harmonic biases resulting from constraints on signal power (in central vision) versus signal uncertainties in orientation (in peripheral vision). The proposed model is consistent with recent work by computational neuroscientists in supposing that visual bias reflects the adjustment of a rational system in the light of uncertain signals and system constraints.

The separation of secondary oil-water dispersions in particulate beds

The literature relating to haze formation, methods of separation, coalescence mechanisms, and models by which droplets <100 μm are collected, coalesced and transferred, have been reviewed with particular reference to particulate bed coalescers. The separation of secondary oil-water dispersions was studied experimentally using packed beds of monosized glass ballotini particles. The variables investigated were superficial velocity, bed depth, particle size, and the phase ratio and drop size distribution of inlet secondary dispersion. A modified pump loop was used to generate secondary dispersions of toluene or Clairsol 350 in water with phase ratios between 0.5-6.0 v/v%.Inlet drop size distributions were determined using a Malvern Particle Size Analyser;effluent, coalesced droplets were sized by photography. Single phase flow pressure drop data were correlated by means of a Carman-Kozeny type equation. Correlations were obtained relating single and two phase pressure drops, as (ΔP2/μc)/ΔP1/μd) = kp Ua Lb dcc dpd Cine A flow equation was derived to correlate the two phase pressure drop data as, ΔP2/(ρcU2) = 8.64*107 [dc/D]-0.27 [L/D]0.71 [dp/D]-0.17 [NRe]1.5 [e1]-0.14 [Cin]0.26  In a comparison between functions to characterise the inlet drop size distributions a modification of the Weibull function provided the best fit of experimental data. The general mean drop diameter was correlated by: q_p q_p p_q /β      Γ ((q-3/β) +1) d qp = d fr  .α        Γ ((P-3/β +1 The measured and predicted mean inlet drop diameters agreed within ±15%. Secondary dispersion separation depends largely upon drop capture within a bed. A theoretical analysis of drop capture mechanisms in this work indicated that indirect interception and London-van der Waal's mechanisms predominate. Mathematical models of dispersed phase concentration m the bed were developed by considering drop motion to be analogous to molecular diffusion.The number of possible channels in a bed was predicted from a model in which the pores comprised randomly-interconnected passage-ways between adjacent packing elements and axial flow occured in cylinders on an equilateral triangular pitch. An expression was derived for length of service channels in a queuing system leading to the prediction of filter coefficients. The insight provided into the mechanisms of drop collection and travel, and the correlations of operating parameters, should assist design of industrial particulate bed coalescers.

Nucleoside and oligonucleotide approaches towards potential HIV chemotherapies

The use of oligonucleotides directed against the mRNA of HIV promises site-specific inhibition of viral replication. In this work, the effect of aralkyl substituents on oligonucleotide duplex stability was studied using model oligonucleotide sequences in an attempt to improve targeting of oligonucleotides to viral mRNA. Arakyl-substituted oligonucleotides were made by solid phase synthesis using either the appropriate aralkyl-substituted phosphoramidite or by post-synthetic substitution of a pentafluorophenoxy substituent by N-methylphenethylamine. The presence of phenethyl or benzoyl substituents invariably resulted in thermodynamic destabilisation of all duplexes studied. The methods which were developed for the synthesis of nucleoside intermediates for oligonucleotide applications were also used to prepare a series of nucleoside analogues derived from uridine, 2'-deoxyuridine and AZT. Crystal structures of six compounds were successfully determined. Anti-HIV activity was observed for most compounds in the series although none were without cytotoxicity. The most active compound of the series was the ribose nucleoside; 1-β-D-erythro-pentofuranosyl-4-pentafluorophenoxy-pyrimidine-2(1H)-one 95, derived directly from uridine. The same series of compounds also displayed very modest anti-cancer activity. To enable synthesis of prooligonucleotides and analogues for possible antisense applications, the properties of a new Silyl-Linked Controlled Pore Glass solid support were investigated. Synthesis of the sequences d(Tp)7T, d(Tps)7T and the base-sensitive d(Tp)3(CBzp)2(Tp)2T was achieved using the silyl-linked solid support in a fluoride-induced cleavage/deprotection strategy.

Cyclic and acyclic modifications of 5-aminoimidazole-4-carboxamide

The Introduction gives a brief resume' of the biologically important aspects of 5 -aminoimidazole -4 -carbozamide (1) and explores., in-depth, the synthetic routes to this imidazole. All documented reactions of 5 -aninoimidanole-4 -carboxamide are reviewed in detail, with particular emphasis on the preparation and subsequent coupling reactions of 5 –diazo-imidazole-4 -carboxamide (6). A series of thirteen novel amide 5-amino-2-arylazoimidazole-4-carboxamide derivatives (117-129) were prepared by the coupling of aryldiazonium salts with 5-aminoimidazole-4-carboxamide. Chemical modification of these azo-dyes resulted in the preparation of eight previously unknown acyl derivatives (136-143) Interaction of 5-amino-2-arylazoimidazole-4-carboxides with ethyl formate in sodium ethoxide effected pyrimidine ring closure to the novel 8-arylazohypoxanthines (144 and 145). Several reductive techniques were employed in an effort to obtain the elusive 2,5-diaminoimidazole-4-carboxamide (71),a candidate chemotherapeutic agent, from the arylazoiridazoles. No success can be reported although 5-amino-2-(3-aminoindazol-2-yl) imidazole-4-carboxamide (151) was isolated due to a partial reduction and intramolecular cyclisation of 5-amino72-(2-cyanaphenylazo)imidazole-4-carboxamide (122) .Further possible synthetic approaches to the diaminoimidazole are discussed in Chapter 4. An interesting degradation of a known unstable nitrohydrazone is described in Chapter 5.This resulted in formation of 1, 1-bis(pyrazol--3-ylazo)-1-nitroethane (164) instead of the expected cyclisation to a bicyclic tetrazine N-oxide. An improved preparation of 5-diazoinidazole-4-carboxamide has been achieved, and the diazo-azole formed cycloadducts with isocyanates to yield the hitherto unknown imidazo[5,1-d][1,2,3,5]tetrazin-7(6H)-ones. Eleven derivatives (167-177) of this new ring-system were prepared and characterised. Chemical and spectroscopic investigation showed this ring-system to be unstable under certain conditions, and a comparative study of stability within the group has been made. "Retro-cycloaddition" under protic and photolytic conditions was an unexpected property of 6-substituted imidazo[5,1-d][1,2,3,5]tetrazin--7(0)-ones.Selected examples of the imidazotetrazinone ring-system were tested for antitumour activity. The results of biological evaluation are given in Chapter 7, and have culminated in a Patent application by the collaborating body, May and Baker Ltd. One compound,3-carbamoyl-6-(2-chloro-ethyl)imidazo[5,1-d][1,2,3,5jtetrazin-7(6H)-one (175),shows striking anti-tumour activity in rodent test systems.

Antitumour imidazotetrazines: probes for the major groove of DNA

The antitumour imidazotetrazinones are believed to act as prodrugs for the triazene series of alkylating agents, showing a marked pteference for the alkylation of the middle guanine residue in a run of three or more contiguous guanines. However, the. exact nature of the interactions of imidazotetrazinones within the micro~environment of DNA are; as yet unknown. In order to examine such interactions a three pronged approach involving molecular modelling, synthetic chemistry and biological analysis has been undertaken during the course of this project. . Molecular modelling studies have shown that for the 8-carboxamido substituted imidazotetrazinones antitumour activity is dependent upon the. presence of a free NH group which can be involved in the formation of both intramolecular and intermolecular hydrogen bonds, and the presence of a non-bulky substituent with a small negative potential . volume. Modelling studies involving the docking of .mitozolomide into the major groove of DNA in the region of a triguanine sequence has shown that a number of hydrogen bonding interactions are feasible. A series of 8-substituted carboxamide derivatives of mitozolomide have been synthesised via the 8-acid chloride and 8-carboxylic acid derivatives including a number of peptide analogues. The peptide derivatives were based upon the key structural features of the helix-turn-helix motif of DNA-binding proteins with a view to developing agents that are capable of binding to DNA with greater selectivity. An examination of the importance of intramolecular hydrogen bonding in influencing the antitumour activity:of :the imidazotetrazinones has led to the synthesis of the novel pyrimido[4',5' :4,3]pyrazolo[5,1-d]-1,2,3,5-tetrazine ring system. In general, in vitro cytotoxicity assays showed that the new derivatives were less active against the TLX5 lymphoma cell line. than the parent compound mitozolomide despite an increased potential for hydrogen bonding interactions. Due to the high reactivity of the: tetrazinone ring system it is difficult to study the interactions between the imidazotetrazinones and DNA. Consequently a number of structural analogues that are stable under physiological conditions have been. prepared based upon the 1,2,3 triazin-4(3H)-one ring system fused with both benzene and pyrazole rings. Although the 3-methylbenzotriazinones failed to antagonise the cytotoxic activity of temozolomide encouraging results with a 3-methylpyrazolotriazinone may suggest the existence of an imidazotetrazinone receptor site within DNA. The potential of guanine rich sequences to promote the alkylating selectivity of imidazotetrazinones by acting as a catalyst for ring cleavage and thereby generation of the alkylating agent was examined. Experiments involving the monitoring: of the rate of breakdown of mitozolomide incubated in the presence of synthetic oIigonucleotides did not reveal any catalytic effect resulting from the DNA. However, it was noted that the breakdown of mitozolomide was dependent upon the type of buffer used in the incubations and this may indeed mask any catalysis by the oligonucleotides.

The Riesz transform and simultaneous representations of phase, energy and orientation in spatial vision

To represent the local orientation and energy of a 1-D image signal, many models of early visual processing employ bandpass quadrature filters, formed by combining the original signal with its Hilbert transform. However, representations capable of estimating an image signal's 2-D phase have been largely ignored. Here, we consider 2-D phase representations using a method based upon the Riesz transform. For spatial images there exist two Riesz transformed signals and one original signal from which orientation, phase and energy may be represented as a vector in 3-D signal space. We show that these image properties may be represented by a Singular Value Decomposition (SVD) of the higher-order derivatives of the original and the Riesz transformed signals. We further show that the expected responses of even and odd symmetric filters from the Riesz transform may be represented by a single signal autocorrelation function, which is beneficial in simplifying Bayesian computations for spatial orientation. Importantly, the Riesz transform allows one to weight linearly across orientation using both symmetric and asymmetric filters to account for some perceptual phase distortions observed in image signals - notably one's perception of edge structure within plaid patterns whose component gratings are either equal or unequal in contrast. Finally, exploiting the benefits that arise from the Riesz definition of local energy as a scalar quantity, we demonstrate the utility of Riesz signal representations in estimating the spatial orientation of second-order image signals. We conclude that the Riesz transform may be employed as a general tool for 2-D visual pattern recognition by its virtue of representing phase, orientation and energy as orthogonal signal quantities.

Molecular dynamics simulation of methane in sodium montmorillonite clay hydrates at elevated pressures and temperatures

Computer simulation has been used to study the structure and dynamics of methane in hydrated sodium montmorillonite clays under conditions encountered in sedimentary basins. Systems containing approximately one, two, three and four molecular layers of water have followed gradients of 150 bar km-1 and 30Kkm-1, to a maximum burial depth of 6 km (900 bar and 460 K). Methane is coordinated to approximately 19 oxygen atoms, of which typically 6 are provided by the clay surface. Only in the three- and four-layer hydrates is methane able to leave the clay surface. Diffusion depends strongly on the porosity (water content) and burial depth: self-diffusion coefficients are in the range 0.12 × 10-9m2s-1 for water and 0.04 × 10−9m2s−1 < D < 8.64 × 10−9m2s−1 for methane. Bearing in mind that porosity decreases with burial depth, it is estimated that maximum diffusion occurs at around 3 km. This is in good agreement with the known location of methane reservoirs in sedimentary basins.

Accuracy and optimal sampling in Monte Carlo solution of population balance equations

Implementation of a Monte Carlo simulation for the solution of population balance equations (PBEs) requires choice of initial sample number (N0), number of replicates (M), and number of bins for probability distribution reconstruction (n). It is found that Squared Hellinger Distance, H2, is a useful measurement of the accuracy of Monte Carlo (MC) simulation, and can be related directly to N0, M, and n. Asymptotic approximations of H2 are deduced and tested for both one-dimensional (1-D) and 2-D PBEs with coalescence. The central processing unit (CPU) cost, C, is found in a power-law relationship, C= aMNb0, with the CPU cost index, b, indicating the weighting of N0 in the total CPU cost. n must be chosen to balance accuracy and resolution. For fixed n, M × N0 determines the accuracy of MC prediction; if b > 1, then the optimal solution strategy uses multiple replications and small sample size. Conversely, if 0 < b < 1, one replicate and a large initial sample size is preferred. © 2015 American Institute of Chemical Engineers AIChE J, 61: 2394–2402, 2015

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

Purpose: We have reported that the changes in the accommodative response to electrical stimulation of the branches of the ciliary nerves in cats. (Miyagawa et al, PLoS One, 2014). We have also reported that no robust accommodative responses to the electrical stimulations of the sclera of peripheral cornea (SSPC) were observed in enucleated porcine eyes (Mihashi et al, VPOptics, 2014). In this study, accommodative responses to SSPC stimulation in cats and porcines were investigated. Methods: Two eyes of two cats under anesthesia and after they were sacrificed were studied. Three enucleated porcine eyes obtained from a local slaughterhouse were also studied. Trains of biphasic pulses (current, 3 mA; duration, 2 ms/phase; frequency, 40 Hz) were applied using a tungsten electrode (0.3mm diameter) from several orientations. Wavefront sensing with a compact wavefront aberrometer (Uday et al J Cataract Refract Surg, 2013) were performed before and 4 s (cat) and 10 s (pig) after the stimulations and wavefront aberrations including spherical errors were analyzed over a 4-mm pupil area. Results: In the first cat under anesthesia, at three out of seven stimulus positions, 0.2 D hyperopic accommodative responses were observed and in two orientations, myopic responses were observed. For the other cat, weak accommodative responses including astigmatic changes were observed. In the sacrificed condition of the second cat, 0.1 D myopic response was observed for one stimulus orientation and the smaller responses were observed at six out of eight stimulus positions. No accommodative responses were elicited for the enucleated porcine eyes. Conclusions: In the anesthetized cats, electrical stimulation of the SSPC induced accommodative responses; the responses were unstable and weaker than the responses by the ciliary nerve stimulations we observed in our previous study. Small accommodative responses were observed after one of two cats had been sacrificed, but no accommodative responses were detected in the enucleated porcine eyes. Further studies are needed to confirm difference in the accommodation functions in the two species.