Mathematical Problems in Engineering

Volume 2008 (2008), Article ID 127827, 8 pages

http://dx.doi.org/10.1155/2008/127827

## Border Figure Detection Using a Phase Oscillator Network with Dynamical Coupling

^{1}Pósgraduação em Engenharia Elétrica, Escola de Engenharia, Universidade Presbiteriana Mackenzie, Rua da Consolação 896, CEP 01302-907, São Paulo, SP, Brazil^{2}Departamento de Engenharia de Telecomunicações e Controle, Escola Politécnica, Universidade de São Paulo, Av. Prof. Luciano Gualberto 380, Travessa 3, CEP 05508-900, São Paulo, SP, Brazil

Received 13 January 2008; Revised 17 April 2008; Accepted 30 May 2008

Academic Editor: Jerzy Warminski

Copyright © 2008 L. H. A. Monteiro et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Oscillator networks have been developed in order to perform specific tasks related to image processing. Here we analytically investigate the existence of synchronism in a pair of phase oscillators that are short-range dynamically coupled. Then, we use these analytical results to design a network able of detecting border of black-and-white figures. Each unit composing this network is a pair of such phase oscillators and is assigned to a pixel in the image. The couplings among the units forming the network are also dynamical. Border detection emerges from the network activity.

#### 1. Introduction

The synchronous firing of neurons seems to be important for accomplishing cognitive tasks such as attention (e.g., [1]), comprehension (e.g., [2]), coordination (e.g., [3]), perception (e.g., [4]), and sensory segmentation (e.g., [5]). Such experimental findings have inspired works about image processing (e.g., recognition, segmentation, symmetry detection) based on the synchronism of coupled oscillators (e.g., [6–11]). In addition, several neural systems, such as spinal cord (e.g., [12, 13]), hippocampus (e.g., [14, 15]), and visual cortex (e.g., [16, 17]), have been modelled by phase oscillators—one of the simplest oscillator models. In spite of its simplicity, this approach is suitable because more complex models for neurons like pulse-coupled and Hodgkin-Huxley-type models can be transformed into a phase oscillator model through coordinate changes [18]. Networks of phase oscillators can be electronically built using phase-locked loops (PLLs) (e.g., [19–22]).

Here we analytically study a phase oscillator model with dynamical coupling, and then use it in order to form a network capable of successfully detecting border of black-and-white figures. Each oscillator of our network corresponds to a first-order PLL (e.g., [8, 11, 17, 23]), which is equivalent to an overdamped pendulum.

The aim of any border detection process is to capture the main structural properties of an image, which can be useful, for instance, for image compression or diagnosis in echocardiography (e.g., [24]). Usually, this kind of image processing involves the use of partial derivatives (gradient, Laplacian; e.g., [25]), which are discretized in the space domain and in the time domain. Here we propose a scheme based on a network where each unit is formed by a pair of dynamically coupled PLLs corresponding to a pixel of the image. Border detection is obtained from the network activity. Notice that such a model can be transformed into a dedicated hardware for executing this image processing task. Also it naturally presents space discretization, because it is composed of a finite number of PLLs. This feature plus the fact that the time domain does need to be discretized can reduce the stability problems related to the algorithms employed for numerically calculating partial derivatives.

#### 2. Model of a Single Unit

Consider that the temporal activities of two first-order PLLs with dynamical coupling are described by the following equations:where and are the phase and the natural frequency of the oscillator , respectively; ( with ) is the connection strength from the oscillator to the oscillator ; () is the coefficient related to the Hebbian connection modification; and represents a natural (exponential) decay. The parameters , , and are positive numbers. Thus, the connection strengths are enhanced when the oscillators are in phase, and they are weakened when they are out of phase. This kind of (synaptic) modification between two oscillators (neurons) was suggested by Hebb [26].

By defining ; ; ; , system (2.1) can be rewritten asThe variables and represent the phase difference and the average connection strength, respectively. Notice that (). The formation of synchronized clusters in networks of phase oscillators described by system (2.2) was investigated by Seliger et al. [27]. In this section we study the asymptotic behaviors of this model.

Synchronism occurs when both nodes oscillate in a common frequency, which means that or . Thus, synchronism can happen if there are constants so that if then and for all time . Notice that a synchronous solution is an equilibrium point of the dynamical system (2.2). Such a solution corresponds to an intersection point of the nullclines () and ().

The stability of an equilibrium point can usually be determined by calculating the eigenvalues of the Jacobian matrix corresponding to system (2.2) linearized around this point. It is asymptotically stable when all eigenvalues have negative real parts (e.g., [28]). For system (2.2), the eigenvalues are the roots of the polynomialAccording to the Routh-Hurwitz criterion, both roots of the polynomial have negative real parts if and (e.g., [29]).

When , the equilibrium points are , , , , where and are asymptotically stable, and and are unstable.

When , there are also four equilibrium solutions if . The points withare asymptotically stable; the ones withare unstable. When , there are only two equilibrium points because two saddle-node bifurcations (e.g., [28]) occur. When , there is not synchronism; however, there is a limit cycle, which corresponds to a closed and isolated trajectory in the state space .

For an autonomous two-dimensional system, Poincaré-Bendixson theorem (e.g., [28]) ensures that there is an asymptotically stable limit cycle in a region of the state space if the vector field points inward everywhere on the boundary of this region, which must not contain any equilibrium point. For system (2.2), it is easy to verify that if , then there is such an attracting trajectory inside the region , where is the rectangle given by for . In this case, oscillates and grows as the time goes by, as shown in Figure 2. Figure 1 illustrates the case with synchronous (equilibrium) asymptotic solution.

In order to analytically characterize the limit cycle, suppose that the asymptotic behavior of in this case is given by . By inserting this approximated (linear) solution in the equation for and integrating it, the asymptotic solution of iswhere and . By substituting this expression in the equation for in (2.2), an equation for calculating is found. In fact, is the root of the polynomialObserve that when , then the angular velocity can be written as , with . The value of is estimated byConsequently, the oscillation amplitude of when is

For instance, for , , and , these approximated expressions give and , which are in good agreement with the numerical solutions of system (2.2) shown in Figure 2.

#### 3. Network for Detecting Figure Border

Our two-dimensional network for detecting figure border is composed of units consisting of a pair of phase oscillators forced by an external input . In this network, the coefficients and can vary with the time. Thus, system (2.2) is rewritten as

Assume that a black-and-white figure will be presented to this network, as illustrated in Figure 3(a). The variable represents the input for each unit according to the following rules: corresponds to a white part of the figure; corresponds to a black part. Notice that the input can be translated into a new natural frequency by defining . Thus, system (3.1) is reduced to system (2.2) if is a constant.

Each unit is coupled with its four closest neighbors (left, right, up, and down), which is usually known as four-neighborhood in image processing literature (e.g., [30]) or von Neumann neighborhood in cellular automaton literature (e.g., [31]). As a consequence of this coupling, the value of for each cell is given byThe parameter is a settling time; where is a positive number (here ); where if ; and if . The value of is obtained bywhere is the time step (here ) of the integration method (here fourth-order Runge-Kutta method) used for numerically solving the dynamical system. The variable corresponds to the average connection strength of the closest neighbor . If the neighbor does not exist, then .

An image is presented to the network at . For , all the units behave as if they were isolated, because the value of for each unit is independent of the neighbor activity. By taking for all units , , , , , and , the units with will tend to a stationary solution and the units with will tend to a limit cycle, as explained in the last section. For these parameter values, the permanent regime was already reached when , as shown in Figures 1 and 2.

At , the value of for each unit can be changed. Expressions (3.2) and (3.3) imply what follows. If all of its four neighbors are in a limit cycle (), then the value of is increased from to ; if at least one neighbor is in a stationary state or if the unit does not have a complete von Neumann neighborhood (), then the value of is kept equal to . The limit cycle is characterized by oscillating with frequency and amplitude , and linearly growing with slope , as presented in the last section.

For and , and tend to a stationary activity, which in our model corresponds to the white color. Thus, at , the unique units remaining in oscillatory activity (the black color) are the ones in the figure border. Figure 3(b) shows the contour of Figure 3(a), after the units have reached the (new) permanent regime.

The higher the value of is, the shorter the transitory phase will be. Hence, border detection can be made faster by increasing the value of this parameter.

#### 4. Conclusion

Networks governed by differential (e.g., [6–8, 10, 11]) or difference equations (e.g., [9, 32]) have been employed for image processing. Here we used a network with local and dynamical coupling for identifying contour of black-and-white images. The unit of such a network is the variant of the Kuramoto model [33] for two oscillators (neurons), which was proposed by Seliger et al. [27]. In our scheme, after applying the input corresponding to the figure, the value of for each unit is switched or not, depending on the neighborhood activity. This scheme could be implemented using first-order PLLs. This investigation could be extended for higher-order PLLs, which could be designed in order to shorten the transient.

#### Acknowledgment

L. Monteiro and J. Piqueira are partially supported by CNPq.

#### References

- P. N. Steinmetz, A. Roy, P. J. Fitzgerald, S. S. Hsiao, K. O. Johnson, and E. Niebur, “Attention modulates synchronized neuronal firing in primate somatosensory cortex,”
*Nature*, vol. 404, no. 6774, pp. 187–190, 2000. View at Publisher · View at Google Scholar - M. A. Just, V. L. Cherkassky, T. A. Keller, and N. J. Minshew, “Cortical activation and synchronization during sentence comprehension in high-functioning autism: evidence of underconnectivity,”
*Brain*, vol. 127, no. 8, pp. 1811–1821, 2004. View at Publisher · View at Google Scholar - S. N. Baker, J. M. Kilner, E. M. Pinches, and R. N. Lemon, “The role of synchrony and oscillations in the motor output,”
*Experimental Brain Research*, vol. 128, no. 1-2, pp. 109–117, 1999. View at Publisher · View at Google Scholar - E. Rodriguez, N. George, J.-P. Lachaux, J. Martinerie, B. Renault, and F. J. Varela, “Perception's shadow: long-distance synchronization of human brain activity,”
*Nature*, vol. 397, no. 6718, pp. 430–433, 1999. View at Publisher · View at Google Scholar - A. K. Engel, P. König, A. K. Kreiter, T. B. Schillen, and W. Singer, “Temporal coding in the visual cortex: new vistas on integration in the nervous system,”
*Trends in Neurosciences*, vol. 15, no. 6, pp. 218–226, 1992. View at Publisher · View at Google Scholar - Y. Hayashi, “Oscillatory neural network and learning of continuously transformed patterns,”
*Neural Networks*, vol. 7, no. 2, pp. 219–231, 1994. View at Publisher · View at Google Scholar - D. L. Wang, “Emergent synchrony in locally coupled neural oscillators,”
*IEEE Transactions on Neural Networks*, vol. 6, no. 4, pp. 941–948, 1995. View at Publisher · View at Google Scholar - F. C. Hoppensteadt and E. M. Izhikevich, “Pattern recognition via synchronization in phase-locked loop neural networks,”
*IEEE Transactions on Neural Networks*, vol. 11, no. 3, pp. 734–738, 2000. View at Publisher · View at Google Scholar - R. de Oliveira and L. H. A. Monteiro, “Symmetry detection using global-locally coupled maps,” in
*Proceedings of the 12th International Conference on Artificial Neural Networks (ICANN '02)*, vol. 2415 of*Lecture Notes in Computer Science*, pp. 75–80, Madrid, Spain, August 2002. View at Zentralblatt MATH - M. Ursino, G.-E. La Cara, and A. Sarti, “Binding and segmentation of multiple objects through neural oscillators inhibited by contour information,”
*Biological Cybernetics*, vol. 89, no. 1, pp. 56–70, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - B. K. Ghosh, A. D. Polpitiya, and W. Wang, “Bio-inspired networks of visual sensors, neurons, and oscillators,”
*Proceedings of the IEEE*, vol. 95, no. 1, pp. 188–214, 2007. View at Publisher · View at Google Scholar - A. H. Cohen, P. J. Holmes, and R. H. Rand, “The nature of the coupling between segmental oscillators of the lamprey spinal generator for locomotion: a mathematical model,”
*Journal of Mathematical Biology*, vol. 13, no. 3, pp. 345–369, 1981. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. Taylor and P. Holmes, “Simple models for excitable and oscillatory neural networks,”
*Journal of Mathematical Biology*, vol. 37, no. 5, pp. 419–446, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. M. Borisyuk and F. C. Hoppensteadt, “Memorizing and recalling spatial-temporal patterns in an oscillator model of the hippocampus,”
*BioSystems*, vol. 48, no. 1–3, pp. 3–10, 1998. View at Publisher · View at Google Scholar - F. K. Skinner, C. Wu, and L. Zhang, “Phase-coupled oscillator models can predict hippocampal inhibitory synaptic connections,”
*European Journal of Neuroscience*, vol. 13, no. 12, pp. 2183–2194, 2001. View at Publisher · View at Google Scholar - Y. B. Kazanovich and R. M. Borisyuk, “Dynamics of neural networks with a central element,”
*Neural Networks*, vol. 12, no. 3, pp. 441–454, 1999. View at Publisher · View at Google Scholar - L. H. A. Monteiro, N. C. F. Canto, J. G. Chaui-Berlinck, F. M. Orsatti, and J. R. C. Piqueira, “Global and partial synchronism in phase-locked loop networks,”
*IEEE Transactions on Neural Networks*, vol. 14, no. 6, pp. 1572–1575, 2003. View at Publisher · View at Google Scholar - E. M. Izhikevich, “Weakly pulse-coupled oscillators, FM interactions, synchronization, and oscillatory associative memory,”
*IEEE Transactions on Neural Networks*, vol. 10, no. 3, pp. 508–526, 1999. View at Publisher · View at Google Scholar - J. R. C. Piqueira, F. M. Orsatti, and L. H. A. Monteiro, “Computing with phase locked loops: choosing gains and delays,”
*IEEE Transactions on Neural Networks*, vol. 14, no. 1, pp. 243–247, 2003. View at Publisher · View at Google Scholar - L. H. A. Monteiro, D. N. F. Filho, and J. R. C. Piqueira, “Bifurcation analysis for third-order phase-locked loops,”
*IEEE Signal Processing Letters*, vol. 11, no. 5, pp. 494–496, 2004. View at Publisher · View at Google Scholar - J. R. C. Piqueira, M. Q. Oliveira, and L. H. A. Monteiro, “Synchronous state in a fully connected phase-locked loop network,”
*Mathematical Problems in Engineering*, vol. 2006, Article ID 52356, 12 pages, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - J. R. C. Piqueira and M. D. C. Freschi, “Models for master-slave clock distribution networks with third-order phase-locked loops,”
*Mathematical Problems in Engineering*, vol. 2007, Article ID 18609, 17 pages, 2007. View at Publisher · View at Google Scholar - L. H. A. Monteiro, R. V. dos Santos, and J. R. C. Piqueira, “Estimating the critical number of slave nodes in a single-chain PLL network,”
*IEEE Communications Letters*, vol. 7, no. 9, pp. 449–450, 2003. View at Publisher · View at Google Scholar - J. N. Kirkpatrick, R. M. Lang, S. E. Fedson, A. S. Anderson, J. Bednarz, and K. T. Spencer, “Automated border detection on contrast enhanced echocardiographic images,”
*International Journal of Cardiology*, vol. 103, no. 2, pp. 164–167, 2005. View at Publisher · View at Google Scholar - G. Aubert and P. Kornprobst,
*Mathematical Problems in Image Processing*, vol. 147 of*Applied Mathematical Sciences*, Springer, New York, NY, USA, 2nd edition, 2006. View at Zentralblatt MATH · View at MathSciNet - D. O. Hebb,
*The Organization of Behavior*, John Wiley & Sons, New York, NY, USA, 1949. - P. Seliger, S. C. Young, and L. S. Tsimring, “Plasticity and learning in a network of coupled phase oscillators,”
*Physical Review E*, vol. 65, no. 4, Article ID 041906, 7 pages, 2002. View at Publisher · View at Google Scholar · View at MathSciNet - J. Guckenheimer and P. Holmes,
*Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields*, vol. 42 of*Applied Mathematical Sciences*, Springer, New York, NY, USA, 1983. View at Zentralblatt MATH · View at MathSciNet - K. Ogata,
*Modern Control Engineering*, Prentice-Hall, Upper Saddle River, NJ, USA, 2001. - R. C. Gonzalez and R. E. Woods,
*Digital Image Processing*, Addison-Wesley, New York, NY, USA, 2007. - S. Wolfram, Ed.,
*Theory and Applications of Cellular Automata*, vol. 1 of*Advanced Series on Complex Systems*, World Scientific, Singapore, 1986. View at Zentralblatt MATH · View at MathSciNet - L. Zhao and E. E. N. Macau, “A network of dynamically coupled chaotic maps for scene segmentation,”
*IEEE Transactions on Neural Networks*, vol. 12, no. 6, pp. 1375–1385, 2001. View at Publisher · View at Google Scholar - Y. Kuramoto,
*Chemical Oscillations, Waves, and Turbulence*, vol. 19 of*Springer Series in Synergetics*, Springer, Berlin, Germany, 1984. View at Zentralblatt MATH · View at MathSciNet