A NEW FORECAST METHOD OF THE MOVING OPTICAL OBJECT POSITION BASED ON THE PARALLEL HIERARCHICAL NETWORK

LEONID TIMCHENKO1*, NATALIYA KOKRIATSKAYA2, VIKTOR MELNIKOV3, MYCOLA PETROVSKYI4, GALINA KOSENKO5
1Department of Telecommunication technologies and automatics, State Economy and Technology University of Transport, Kyiv, Ukraine
2Department of Telecommunication technologies and automatics, State Economy and Technology University of Transport, Kyiv, Ukraine
3Department of Telecommunication technologies and automatics, State Economy and Technology University of Transport, Kyiv, Ukraine
4Department of Telecommunication technologies and automatics, State Economy and Technology University of Transport, Kyiv, Ukraine
5Department of Telecommunication technologies and automatics, State Economy and Technology University of Transport, Kyiv, Ukraine
* Corresponding Author : timchen@list.ru

Received : 08-12-2011     Accepted : 15-12-2011     Published : 19-12-2011
Volume : 2     Issue : 2       Pages : 54 - 59
J Signal Image Process 2.2 (2011):54-59

Cite - MLA : LEONID TIMCHENKO, et al "A NEW FORECAST METHOD OF THE MOVING OPTICAL OBJECT POSITION BASED ON THE PARALLEL HIERARCHICAL NETWORK." Journal of Signal and Image Processing 2.2 (2011):54-59.

Cite - APA : LEONID TIMCHENKO, NATALIYA KOKRIATSKAYA, VIKTOR MELNIKOV, MYCOLA PETROVSKYI , GALINA KOSENKO (2011). A NEW FORECAST METHOD OF THE MOVING OPTICAL OBJECT POSITION BASED ON THE PARALLEL HIERARCHICAL NETWORK. Journal of Signal and Image Processing, 2 (2), 54-59.

Cite - Chicago : LEONID TIMCHENKO, NATALIYA KOKRIATSKAYA, VIKTOR MELNIKOV, MYCOLA PETROVSKYI , and GALINA KOSENKO "A NEW FORECAST METHOD OF THE MOVING OPTICAL OBJECT POSITION BASED ON THE PARALLEL HIERARCHICAL NETWORK." Journal of Signal and Image Processing 2, no. 2 (2011):54-59.

Copyright : © 2011, LEONID TIMCHENKO, et al, Published by Bioinfo Publications. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution and reproduction in any medium, provided the original author and source are credited.

Abstract

This paper discusses principles necessary to develop a forecast method based on the parallel-hierarchical (PH) network and an idea of the hyperbolic smoothing of empiric data. Mathematic models for forecasting dynamic image series, as well as software for forecasting moving optical object positions on the example of energy centers of laser beam spot images are developed. The developed method is more efficient and less time consuming for real-time systems in realization of the forecasting operation of energy center positions of laser beam spots.

Keywords

parallel-hierarchical (PH) network, laser beam spot images, statistical observations, dynamic series, forecasting

Introduction

Rapidly growing requirements to modern computational media encourage development of new intelligent methods of information transfer and processing. Rigid requirements to real-time information processing systems force scientists to regularly create and upgrade data transfer systems. Today most internet channels cannot provide data exchange of the required quality between such systems, which, in turn, results in the congestion of those channels and formation of so-called digital bottlenecks. A possible solution of the problem of transfer of large volumes of information is to use a fiber-optic cable, but laying such cable is rather expensive, even on short distances. At the same time, this problem can be solved through application of the laser-based technologies [1-4] , one of the most promising models of information transfer for the near future. In this case, for instance, tons of full-length films and virtual worlds could be transferred to any part of the globe in a blink of an eye.
Most satellites transmit information, such as TV programs, by means of microwave radiation, while laser-based information transmission could be hundreds of times faster, which, in turn, will considerably increase the carrying capacity of the channel.
The laser-based information transmission requires that both a satellite and a receiving unit (RU) were located in a certain position. A position of the RU lens, whose diameter is only several centimeters, must be adjusted to one thousandth degree, otherwise the information transmission will not happen.
One of the main problems in the process of satellite monitoring by the receiving unit, which occurs at all stages of the system operation, is the forecast of the laser beam spot image location, namely of its geometrical characteristics, which are being distorted by the turbulence and air masses.
Significantly noise-distorted images, in their turn, can both considerably worsen results of the spot location forecast, and set a system into the state when it will not be able to respond adequately to variations in the monitored object location.
To solve a problem of the efficient forecast, frames of the laser beam spot (LBS) image sequence should be classified with a goal to filter a laser route from images considerably distorted by noise, thus forming a tunnel of reference images.
Development of computerized real-time forecasting systems gained new momentum with the emergence of high-efficiency automated systems of information collection, processing and storage. Such systems consist of a complex of information collection and processing tools which are able, based on a system of specially developed algorithms, to solve a problem of the classification and forecasting of the characteristics of respective objects, phenomena or real time processes, in most cases using geometric characteristics of object boundaries. Fragments of laser beam routes used in the optical communication, navigation, ranging, and military equipment, are studied in this article. Development of the PH network-based forecasting system allows solving the following problems:
- Automatic object control;
- Laser-based data transfer;
- Forecasting of the moving optical object behavior.

Principles of the ph network organization

An image analysis consists of the consecutive transformation of the matched image components and determination (filtration) of the unmatched in time image components during neural network elements transfer from the current energetic state with certain coordinates into a state with lower energy and other space coordinates. A condition of the image transfer components to the next level is a presence of the mutual coincidence intermediate results dynamics of time processing in the corresponding low level channels. An analysis result is formed of image components isolated in the space-time domain [5] .
Let us consider a mathematical model of the parallel decomposition of a set [4] , used in each branch of the PH network.



Where – dimension of this set, – elements of subsets consisting of identical elements, nk – a number of elements in the kth subset (i.e. the multiplicity of number ), – an element of the set selected at the jth step, .

Let us use an idea of the population coding [6] and construct a model of some end action performed by all current actions.
For the PH network, both a number of elements in the branch of each level, which are determined based on the model (1) and values of the element itself may be taken as averaged parameters. Obviously, at the neural network branch level, this end action corresponds to this network’s averaged parameters when the population coding is implemented [6] In this case the current image being recognized will be reflected by the current PH network and compared with the reference PH network with averaged parameters.
Having denoted an average value of the random element of the first level as , second level as , third level as , etc., the kth last level as , and an average number of elements of the first level as , second level as , third level as , etc. and the kth level as , one can form a PH network with averaged parameters.
A structure of the PH network with averaged parameters thus synthesized is shown at [Fig-1] .
A current image or a tested image prepared by the PH network with current parameters and a respective number of elements in branches of each level is being compared with the reference image prepared by the PH network with averaged parameters and an average number of elements in branches of the respective level



Using a preparation of elements of the PH network for each level [5] , it is possible to proceed from averaged parameters to their representation with binarized preparations (-1, 0, +1) on the basis of the three-level coding. Then arrays of difference of the element with an average brightness value of decomposition elements of the image (or of its fragment where the image element is located) are determined.
For a random averaged parameter, a three-level transition may be presented by three types of preparations: zero - , positive - , and negative - . In this case the PH network with the numerical count is transformed into the PH network with binarized counts [7] . Then a procedure of correlation comparison of binarized counts of the current and reference PH networks becomes substantially simpler.
In order to form reference images, training within a sample should be conducted. For this purpose, it is necessary to perform averaging by branch elements of each level, i.e. to form averaged elements , and then move to binarized preparations . Having completed manipulations described above, it is possible to form a PH network with reference parameters for current images. After that its correlation comparison with the PH network using current parameters can be performed. (Under the PH network with current parameters we understand a PH network with current values of its elements with a transfer to binarized preparations and a current number of elements in the branch of each level . It is extremely important that a correlation coefficient can be calculated not only for each two levels, but also for two PH networks, thus increasing the reliability of a recognition result.

Description of the ph network-based forecasting method

The statistical information analysis makes it possible to reveal cause-and-effect relations of the phenomena under study, to identify impacts and interactions of various factors, thus permitting to evaluate efficiency of managerial decisions made. By comparing generalizing statistical indicators of phenomena understudy, it is possible to determine quantitative attributes of their space dissemination and temporal progress, to identify characteristics of their relation and interconnectedness, and formulate scientific and practical conclusions.
To analyze statistics understudy, it should be systemized by constructing the chronological series, also known as dynamic series. Every dynamic series consists of periods, or time moments t, to which series levels belong, and statistical indices which characterize time levels. Extended laser route data, namely, energy center coordinates of laser beam spots, are used as statistical indicators [8,9] .
An analysis of empirical data smoothing methods conducted in [8] shows that irrespective of a smoothing method or a method of trend line determination, a distribution of “peaks” and “wells” remains the same. It can be concluded therefore that in ducat or fluctuations for most phenomena may be determined by any empirical data smoothing method. That is why the least time consuming method should be used. A hyperbolic smoothing method described below will be used for the second stage of the statistical research, i.e. for the statistical aggregation and grouping of the initial data.

The hyperbolic smoothing method of statistical data

One of the key methods of the dynamic series analysis and generalization is the identification of a main trend of the series. A dynamic series trend line indicates a time change in the phenomenon under study without certain deviations caused by various factors. Atrend of temporal development of the phenomenon can be found according to the methods of increasing intervals, moving average and analytical smoothing methods.
The most efficient is the analytical smoothing method, a complex method of the main trend detection [10] .
Let us consider the dynamic series levels as a time function (2):



A problem of smoothing comes down to finding such form of the function that ordinates of its points are the closest to values of the actual dynamic series.
Most common regularities describing phenomenon development trends are: the straight line, power function, parable of the second and third order, hyperbola, logistic function, the exponent, Fourier series, etc. The analysis shows that the hyperbolic smoothing method using is most effective, because it is the hyperbola that most accurately describes the dynamic series increase or decrease; this smoothing method is also less time consuming compared to other methods [11,12] .
A hyperbola equation may be written as follows:



Where are hyperbola equation parameters. To find those parameters with the least-squared method, combined normal equations of the following type should be used:



If the equality is true, the set (4) will produce:



Having found the logarithm of combined equations (5), we obtain:



The system (6) gives us hyperbola equation parameters of the type:



Development of the forecasting method

Let us discuss an application of the hyperbolic smoothing method of empiric data as described above on a practical example. This example characterizes a fluctuation dynamics of the production costs of a business. Business data are shown in [Table-1] [10] .
Having had the hyperbolic smoothing method (8) applied, and then having found the logarithm of the hyperbole equation, one can find hyperbole equation parameters by formulas (7) and (8):





Therefore, the hyperbole equation looks as follows:



The smoothing data are presented in [Table-2] . See [Table-2] Data from [Table-2] are illustrated by the graph [Fig-2] .
As [Fig-2] shows, the smoothed levels are very close to the empirical ones, therefore, this indicates that the hyperbola equation is appropriate for the trend representation.
Certain approximation of the smoothed series permits to forecast the trend for several steps. The qualitative indicators of the forecast can be found experimentally through determining an optimal discretization step and a number of forecasting steps, as well as the forecasting accuracy.
Consequently, input data – energy center coordinates of laser beam spot images – should be fed to the PH network. After that, a stage of PH network training with the developed combined method takes place.

The combined training method

An algorithm of the combined training method is based on the error correction and the PH network memory use in the image recognition:
1. Construction of the PH network structure of the selected image;
2. Application of a training method based on the memory usage. A number of PH network levels is used in this method as a classifier based on k closest neighbors.
3. Construction of the PH network structure of the reference image.
4. Application of a training method based on the error correction. A criterion used for verification in this method is a correlation coefficient at the zero level of the PH network.
5. Error correction for correlation coefficients is performed at the zero, first, middle, last but one, and third from the end PH network levels. The correction is conducted till network parameters come to an absolute equilibrium state (+1) or exceed it. If parameter values lead to (-1) state (or an absolutely non-equilibrium state), a conclusion is made that image couldn’t classified to current class.
The network parameters are arrays of X and Y coordinated at certain period of time, as well as hyperbole equation parameters a0, a1 at the same period. A structure chart of the developed PH network-based forecasting method is shown at [Fig-3] .

Experimental Research Results

To test the developed method and determine its efficiency, experiments were conducted on forecasting positions of energy centers of laser beam images based on known neural networks [13-15] . The forecasting was performed for one step.
Experimental results of the dynamic series forecasting are shown in [Table-3] .
Chart obtained for the dynamic series forecasting (for MP 5-5-4) is shown at [Fig-4] , [Fig-5] .
The conducted research indicate that average error is practically the same and amounts to 0,77%. The maximum error value of forecasting based on different neural networks varies in the range from 4,07- 4,24%.
Let us perform a forecast with the developed PH network-based software. [Fig-6] and [Fig-7] demonstrate charts obtained by forecasting the dynamic series.
Respective parameters are:
- Average forecast error: 0.552%;
- Maximum value of the forecast error: 1.23%.

Conclusion

Principles necessary to design a PH network-based forecasting method are discussed in this article. An analysis of the statistical information permits to reveal cause-and-effect relations of phenomena under study, to determine an effect and interaction of various factors, therefore permitting to evaluate efficiency of the managerial decisions made. For the second stage of the statistical research, i.e. aggregation and grouping of initial data, the hyperbolic smoothing method is used. It was found out that the hyperbolic smoothing method is the most efficient, because it is the hyperbole that most accurately describes increase and decrease of the dynamic series. In addition, this method is less time-consuming in comparison with other methods. If certain approximation of the smoothed series is conducted, a trend forecast can be conducted for several steps. Qualitative parameters of the forecast are determined experimentally through calculation of the optimal discretization step and a number of forecasting steps, as well as of the forecast accuracy.
To automatize the forecasting process and to measure accurately coordinates of energy centers of laser beam spot images, special software was developed. This software contains a window reflecting a picture of a dynamics of laser beam spots by X- and Y-axes, and a prognostic chart of laser beam movement by X- and Y-axes. Fixation of the gravity center occurs dynamically in real-time with fixation of the following parameters:
- Position of the beam energy centers by X and Y axes respectively in pixel with accuracy to three decimal digits;
- Deviation of the beam energy center by X and Y axes from average values in pixels accurate within third decimal digit;
- Frame number.

Acknowledgment

Authors thank to General Constructor of the CJSC KIA Systems (Moscow, Russia) Yu.F. Kutayev for problem formulation at the initial stage of the research and for video data of fragments of laser rout images made available for the scientific research. Authors express their sheer gratitude to the team of the Chair of the Telecommunication Technologies and Automatics of the State Economy and Technology University of Transport (Kyiv, Ukraine) for their valuable comments and recommendations expressed in the course of experimental research. Authors also thank Professor M.V. Makarenko, Vice Rector on Scientific Work, Associate Professor V.I. Danilevskiy, and Dean of the Department of Infrastructure and Rolling Stock of Railway Transport Professor A.D. Voznenko for provided facilities and equipment for scientific studies of this problem in the framework of the state-financed scientific research topics.

References

[1] Khosrofian J.M. and Garetz B.A. (1983) Appl. Opt. 22, 3406–3410.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[2] Basov N.G., Zemskov E.M., Kutaev Y.F., Kazansky V.M., Kirillov G.A., Kulikov S.M., Nosach O.Y., Sukharev S.A. (1998) Proc. GCL/HPL 98, SPIE Symposium, St-Petersburg (Russia), vol. 3574, 219-228.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[3] Timchenko L.I., Kutaev Y.F, Kozhemyako V.P. (2002) Proceedings of SPIE, vol. 4790, 465-479.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[4] Y.P. Gowramma, C.N. RaviKumar (2010) Journal of Signal and Image Processing, Vol. 1, Issue 1,.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[5] Timchenko L.I. (2000) Cybernetics and Systems Analysis, 36, No. 2, 251–267.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[6] Hinton J.E. (1992) Scientific American [Russian translation], Nos. 11–12, 103–110.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[7] Kozhemyako V.P., Kutaev Y.F., Timchenko L.I., Chepornyuk S.V., Hamdi R.R., Gertsiy A.A., Ivasyuk I.D. (1998) in: Proc. International ICSC/IFAC Symposium on NEURAL COMPUTATION–NC’98, Vienna, 287-291.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[8] Giovanni Mana, Enrico Massa, and Alessandro Rovera (2001) Appl. Opt., vol. 40, Issue 9, 1378-1385.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[9] Orlov D.A. and Neverova E.A. (2003) Measurement Techniques. Vol. 53, No 10, 1140-1146.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[10] Paolo Giudici (2003) Applied data mining: statistical methods for business and industry, John Wiley and Sons, 111-117.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[11] Brockwell J. and Davis R.A. (2002) Introduction to Time Series and Forecasting, second edition, Springer-Verlag, New York.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[12] Box George, Jenkins Gwilym (1994) Time series analysis: forecasting and control, rev. 3, Prentice Hall PTR Upper Saddle River, NJ, USA, 592.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[13] Haykin S. (1999) Neural Networks, 2nd ed., New Jersey: Prentice Hall, 1104.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[14] Robert J. Howlett, Jain L.C. (2001) Radial basis function networks 2: new advances in design, Springer, 360.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

[15] Dimitri P. Bertsekas (1991) MIT Press, 109-122.  
» CrossRef   » Google Scholar   » PubMed   » DOAJ   » CAS   » Scopus  

Images
Fig. 1- A structure of the PH network [4] with averaged parameters.
Fig. 2- Use of the hyperbolic smoothing for further forecasting
Fig. 3- The structure chart of the developed PH network-based forecasting method
Fig. 4- Forecasting with the MP 5-5-4 neural network
Fig. 5- Forecasting with the MP 5-5-4 neural network (enlarged fragment).
Fig. 6- Forecasting with the PH network
Fig. 7- Forecasting with the PH network (enlarged fragment
Table 1- Business production data
Table 2- Smoothing Data
Table 3- Experimental research results on dynamic series forecasting obtained with the known types of neural networks.