Image Tracking for the High Similarity Drug Tablets Based on Light Intensity Reflective Energy and Artificial Neural Network

It is obvious that tablet image tracking exerts a notable influence on the efficiency and reliability of high-speed drug mass production, and, simultaneously, it also emerges as a big difficult problem and targeted focus during production monitoring in recent years, due to the high similarity shape and random position distribution of those objectives to be searched for. For the purpose of tracking tablets accurately in random distribution, through using surface fitting approach and transitional vector determination, the calibrated surface of light intensity reflective energy can be established, describing the shape topology and topography details of objective tablet. On this basis, the mathematical properties of these established surfaces have been proposed, and thereafter artificial neural network (ANN) has been employed for classifying those moving targeted tablets by recognizing their different surface properties; therefore, the instantaneous coordinate positions of those drug tablets on one image frame can then be determined. By repeating identical pattern recognition on the next image frame, the real-time movements of objective tablet templates were successfully tracked in sequence. This paper provides reliable references and new research ideas for the real-time objective tracking in the case of drug production practices.


Introduction
Nowadays, the automation level of drug mass production and quality inspection still remains undeveloped; most drug tablets in random position distribution were manually inspected for picking out those unqualified ones and then hand-encapsulated after amount counting or bottle packaging by ocular estimation. All drug tablets were highly similar or identical to each other in shape, size, surface topography, and color, which imposes colossal operation overload and time exhaustion on the working operators. In the interests of dealing with these difficulties, tablet tracking by computer video was greatly hampered by their identical size and similar shape, keeping a far distance from the traditional applications of image recognition, such as face recognition or vehicle video monitoring. Nevertheless, light intensity reflective energy and ANN recognition have remarkable superiority in locating those scattered tablet objectives, which undoubtedly confirms its reliable results and high efficient performance. Based on this theoretical hypothesis, the labor intensities will be greatly reduced, and simultaneously the high precision of tracking processes can also be ensured.
According to the recently published literatures, it can be learned that image tracking or video tracking were generally employed in such conditions as moving vehicle monitoring or people tracking. Wang and Wang [1] have proposed an approach based on computer version and image processing technique for automatic grain tracking; through an efficient box-tracking-based method, presented by Fukuda et al. [2], mammalian cells (fibroblast) on microcarrier have been clearly recognized. At the same time, Dill et al. [3] studied the tracking cells supported by digital image processing and obtained a rather satisfactory result for multiobjective positioning during their high-speed movement. Kuba et al. [4] also made their contributions on automatic particle detection by one-class SVM from microscope image. Similar performances have been obtained by González et al. and so forth [5] on tracking olive trees in high-resolution satellite images.
On the other hand, Guarino et al. [6] have optimized the traditional image analysis approach. But their optimization 2 Computational and Mathematical Methods in Medicine process needs a much more computation storage and time interval than that of other traditional ones, which restricts its direct usages in real-time drug mass production; Marçal [7] presented new methods based on mathematical morphology and were suitable for grains of circular shape; Kim et al. [8] have successfully proposed a real-time approach for tracking the number of passing people by using one single camera. Similar research founding proposed by Professor Luque-Baena et al. can also be learned from [9,10] as well. In the area of moving-objective image detection and its automatic classification, Wong et al. [11] used an example of outdoor people tracking for tourist-flow estimation in a constrained environment. According to the research result of Buczkowski et al. [12], the box-tracking method (BCM) shows a very high percentage of error in experimental practice, due to the identical shape and high-speed movement of those drug tablets to be studied. In order to cover this deficiency, Wang et al. [13] have proposed a new and efficient framework for pedestrian analysis and tracking, which consists of rule induction classifier and linear regression model. Agustin and Oh [14] preprocessed the detected foreground objects to eliminate pixel noise and small artifacts by performing opening morphology operation. Besides, Lien et al. [15] also paid high attention on a novel vehicle detection method without background modeling simultaneously. Zhan and Luo [16] made their contribution on system design of realtime vehicle recognition based on video for Windows (AVI) files; the development of a block-based real-time tracking system has been concluded by Park et al. [17] and Głowacz et al. [18] and Park et al. [19] have paid their attentions to the optical flow calculation and the area-based decision rule, respectively, which helps in acquiring reliable image tracking results. When discussing the cross-correlation in pattern recognition, Hong [20] used correlation coefficient of fuzzy numbers under arithmetic operations. Park et al. [21] have studied the correlation coefficient of interval-valued intuitionistic probability sets to multiple attribute group decisions. Ye [22] and Son [23], respectively, proposed probability decision-making method and fuzzy entropy determination for investigating the mutual-relationships among different factors. Other related research results can be learned from [24][25][26]. Furthermore, literature [27][28][29] also presented their latest progresses in cross-correlation mechanisms while such evaluation approaches as intuitionistic probability, intervalvalued probability sets, and gamma rank correlation coefficient were, respectively, described in [30][31][32].
Since these commonly used approaches were found being influenced by traditional limitations, the real-time tablet tracking cannot be ensured easily; more importantly, due to the fact that most existing image tracking technologies focus on human tracking or vehicle monitoring, which can be distinguished easily by the different shapes and the widely divergent sizes of those targeted objects in practice [33][34][35], the high similarity objective tracking method suitable for high-speed drug mass production still remains unstudied and undeveloped, that becomes a research margin of advanced video inspection in these years. This paper is structured as follows. Section 1 outlines the importance and necessity of objective tablet tracking.

CCD imaging head with light intensity sensor array The project light ray
The reflective light ray Figure 1: The light reflective effect caused by different surface topography areas, which composes an original property demonstration for calibrating the inflective-energy surface properties of drug tablets.
Section 2 describes the theoretical foundation of light intensity reflective energy surface. Section 3 presents mathematical properties for describing the reflective energy surfaces of drug tablets, Section 4 discusses the detailed tracking processes, with tracking discussions and performance comparisons mentioned in Section 5, and finally Section 6 concludes this paper as respected.

The Theoretical Foundation of Light Intensity Reflective Energy Surface
During the imaging process, the light was projected vertically thereafter a clear image can be obtained for showing the surface topography and geometrical shape of objective tablet. But, at the same time, the unsmooth surface topography on those drug tablets can also be identified due to their different light reflective effects. Figure 1 shows the light reflective effect according to different topography areas, which composes an original property demonstration for calibrating the inflectiveenergy of objective tablet in this experiment. Since different topographic areas on objective tablet present different reflective effects of light, an average-distributed reflective energy condition of light intensity can be obtained on the smooth part of tablet surface. When the lighting ray casts on the objective surface to be observed, it will be reflected to CCD unit in the form of roughly parallel light routes; thus, an image with uniform distribution of light intensity will be gotten and a low concentration of lightness energy can also be shown. On the contrary, the light illumination on the rugged areas of tablet surface brings about the diffuse reflection of lighting rays, which results in the ununiformed concentration of lightness energy. Based on this theoretical presumption, the topography properties of objective drug tablet can be clearly classified with the reflective light intensity measurement, the computation of illumination energy, and the characteristic modeling of energy distribution.

Surface Mathematical Properties
Since the reflective energy of light intensity demonstrated by one given image pixel can be regarded as the height value in the axis, thereafter the reflective energy distribution of the whole tablet image can be described in the form of free-form surface, from the perspective of spatial surface fitting domain. Based on this theoretical foundation the following physical Computational and Mathematical Methods in Medicine   3 properties have been proposed, for the sake of calibrating the established reflective energy surface of objective tablet in a quantitative way [36,37].

Property 1.
As it was well known that surface elasticity demonstrates the transformation resisting capability of one objective spatial surface shape under the force impact caused by external force loading, it becomes a typical index to calibrate the shape properties of spatial surfaces in a mathematical coordinate system, and simultaneously surface elasticity shows the concentricity level of characteristic distribution in the form of elasticity value as well. In this experiment, Property 1 demonstrates the elasticity distribution concentrating on one certain topography section, which makes it being conveniently used to describe topography properties from the perspective of elasticity variance distribution: Here, denotes an objective topography surface in the form of B-spline basis function; , V , , VV , and V are the partial derivatives of in the first order, second order, and hybrid state of , V axes, respectively; 1 , 2 , 1 , and 2 are the given coefficients, ( , V) denotes a given function of surface vector, and finally , denote the order amounts of surface vector in , V axes.

Property 2.
Energy of surface plays a prominent role as fairness function in the occasions of geometric modeling or surface microproperty analysis. Therefore, energy distribution can be computed and quantified to denote its belonged experimental conditions and geometric properties. This property defines the scattering level of surface energy in a mathematical sense, for a high value shows a more decentralized scattering of surface energy: Here, ( ), ( ), V (V ), VV (V ), and V ( V ) denote the first order, second order, and hybrid derivatives of objective surface ( , V) in , V axes. Property 3. Since surface construction highly depends on external loading effect, the fitted result will be deformed in an obvious scale. In the purpose of quantifying the influential effects caused by external loading, this property was newly proposed to calibrate the difference deviation and variation principle between the forced and original surfaces, as high amendment quantity demonstrates a relative critical deformation of resultant topography, with the definition shown as follows: Here, , ( ), , V (V), and , ( V) denote the B-spline boundary control surfaces in , V axes, respectively, with , denoting the transitional vector between two adjacent control vertexes impacted by external loading.

Property 5.
As the traditional Zernike moments have been extensively used, they receive much research attention in a number of fields: object recognition, terrain reconstruction, roughness segmentation, edge detection, and biomedical measurement. They have already become one of the most widely used families of orthogonal mathematical moments, owing to their extraordinary properties of being invariant or insensitive to any arbitrary rotation of the objective surface in the three-dimensional (3D) coordinate system. Based on the radial polynomial computation of one coordinate point ( , , ) mentioned before, an improved Zernike moment was presented by focusing on their magnitude values as follows: Property 6. ( , V) (0 < , V < 1) was supposed to be the objective surface labeled by control vertexes in direction, with the node vectors being ( 0 , 0 , 0 , 0 , . . . , −4 , −3 , , 0 = 0, −3 = 1, and control vertexes in V direction, with the node vectors being Two signless integral numbers 1 , 2 were employed for Here, 1 , 2 ( , V) denotes the approximated surface of ( , V) with 1 × 2 control vertexes, which can be obtained from the frequency decomposition of ( , V) in its 3D coordinate lattice.

Experiment for Drug Tablet Tracking
4.1. The Arrangement of Experimental Device. Stainless steel material was employed to manufacture the experimental platform, with the purpose of avoiding the high light reflective index and mirror effect of background material, causing these two physical phenomena to reduce the contrast ratio and blur the objective clearance in video image remarkably [38,39]. On the stainless steel pathways, the drug tablets were scattered into the monitoring field for high-speed video detection. Figure 2 illustrates the video inspection platform employed, including the video device for tablet tracing and the platform pathways for tablet transmission. For the convenience of tracing these tablet objectives with high accuracy and efficiency, we subdivided the whole imaging field into several blocks by the size of 100 mm × 100 mm and then used the DMIRM CCD high-speed imaging system manufactured by German LEICA to capture the moving processes of drug tablets at the rate of 24 frames per second (FPS). It is noteworthy that this video imaging system consists of three essential components: a precise high-speed CCD head equipped with light intensity sensor array, a PC equipped with video processing board, and a set of multidimensional triangulation processing software for ensuring high clearances of video image frames. Since the light source being used is a visible coaxial LED (220 V; 20 W; 700 nm) combined with diffraction grating that generates a fixed number of stripes, which enables the system to capture more than 450,000 pixels in one single image frame. The range measurements were computed by using triangulation operation between the binary digital images and the calibrated imaging plane locations. The high-speed CCD microscopic camera was able to capture image pixel data in 0.01 ms (enabling image acquisition from a continuously moving manipulator) and process data array in 10 ms and has a measurement deviation of 1/20000 in the field of view. This CCD imaging head has been attached to a DMIRM 3-axis directly driven configuration manipulator for self-adaptively video imaging; thus, it eliminates image distortion caused by inhomogeneous light and keeps a rapid reaction and stable tracing performance simultaneously. For the purpose of improving the detection precision or dynamic tracing of drug tablets to the precise scale (1∼5 mm), we predetermined the imaging depth as 0.3 m; therefore, the planar platform can be focused. On the other hand, the position overlap of drug tablets can be effectively avoided due to the continual highfrequency vibration of the experimental platform pathways, which reduces the probability of erroneous judgment to the maximum extent. In this experiment, the exposure time for digital imaging was supposed to be 0.5 s whereas the moving speed of drug tablet is about 3.5 m/s; therefore, one given tablet moves 1.75 m during the whole imaging process. In order to reduce the blur effect on the captured images, image sharpening was used to keep the high resolution of those studied objective tablets and thus accurate pattern recognition can be ensured.

The Computation of Transitional Vectors between Surface Control Points.
In order to establish surface model for describing the light intensity reflective energy of drug tablets, the obtained energy value at the specific coordinate position of one pixel can be demonstrated as the -axis value in the vertical direction; therefore, it was considered as the spatial coordinate of one control point. Under usual illumination, the external surface of drug tablet presents a constant situation of light intensity and uniform distribution of grey level. Based on this precondition, the three-dimensional (3D) transitional vectors can be computed.
During the capturing of digital images, the observer sight line keeps orthogonal intersecting with the imaging plane called . The light intensity of one surface point ( , ) can be denoted by ( , ), its project vector was determined as [ , , −1], and then the normal vector of its reflective energy surface can be described as [ , , −1]. In this paper, denotes the included angle that lies between the project vector of lighting and the normal vector of reflective energy surface. The reflective intensity of lighting energy in the direction of normal vector can be computed as [40,41] ( , ) = ( , ) cos .
Here, ( , ) denotes the reflective energy intensity of the image pixel ( , ), whereas ( , ) denotes the incident illumination intensity and is supposed as a constant number; Computational and Mathematical Methods in Medicine denotes the light reflectivity level of objective tablet surface (constant number).
Then, the cosine value of this included angle can be expressed as [42] cos = + + The energy intensity of one image pixel (or called as the reflective light intensity at the objective pixel) can be identified as Figure 3 shows the transitional vectors between surface control points and their spatial distribution. Through presenting a specific image pixel at ( , ), this research obtained the normal vector of energy surface: ( , ) = [ 1 , 1 , 1 ], together with other normal ones at its neighboring: On this basis, the tangent vectors , V , and V at the positions of ( + 1, ), ( , + 1), and ( + 1, + 1) can be determined, respectively: The secondary partial derivative vectors were identified as = ( 11 , 22 , 33 ) and VV = (V 11 , V 22 , V 33 ), which demonstrates the variance ratio of those tangent vectors in and V directions, respectively. Therefore, the finite difference method can be used to determine these two items: . . .
Based on the above-mentioned process, the normal vector ( , ) = [ 1 , 1 , 1 ] of one specific surface point can be determined the same as its partial derivative vectors at the 1st and 2nd order, respectively. Consider

The Establishment of Light Intensity
Reflective Energy Surface. Figure 4 shows some typical image frames for objective 10 th frame 20th frame 30th frame 40th frame 50th frame 60th frame 70th frame 80th frame 90th frame 100th frame 110th frame 120th frame tracking, the amount of drug tablets being seen in the monitoring area was strictly limited under 20 by predetermined setting, in the interest of keeping clear shape description and instantaneous tracking of those objective tablets. Through characteristic statistical analysis of drug tablets in geometric shapes, 12 tablet templates in typical geometric shapes and different laying angles were selected out for the following pattern recognition, with all possible shape morphologies and variant position arrangements during monitoring process being fully included, as Figure 5 demonstrates. Compared with other mathematical criterions, using tablet template facilitates the following determination of their geometrical characteristics in a clearer way. Since energy optimization fitting becomes a frequently used modeling method supported by mathematical programming and structure optimization, it regards the studied reflective energy surface providing the minimum physical distorting energy as its ultimate objective; therefore, various means of energy restriction and employed loadings were applied to achieve this goal in practice [43][44][45]. Energy optimization method regarding the laminose elastic deformation equation of elastic mechanics as its mathematical references: Here, denotes the objective surface represented in , V axes; , V , , and VV denote the first-order and secondorder derivative vectors of in and V axes, respectively.
V denotes a mixing derivative vector, and were given parameters, and denotes a given vector function to be approximated.
In practical modeling, and often denote the referential parameters of material characteristic; denotes the external employed loading used for energy surface control. Under normal condition, the influence results exerted by , , and  will be generally neglected and thereafter valuated as 1, 1, and 0, respectively. Consider By using the basis function of B-spline surface, the abovementioned function was transferred into Here, denotes a control point of surface ( , V); Substituting (16) into (14) and letting = 0, then 8 Computational and Mathematical Methods in Medicine Then, Thus, (17) will be simplified into Here, , , , denotes the integrating operation of oneknown basis function; thus, the mathematical model of energy optimization surface can be transferred into the quadratic function of control vertexes [46]. Based on this surface fitting method and the tablet templates shown in Figure 5, their corresponding reflective energy surfaces of light intensity can then be established with the resultant details being demonstrated in Figure 6.

The Application of ANN for Recognizing Surface Properties.
In this study, back propagation (BP) learning algorithm, which has a unique learning principle, generally called delta rule, is used. As we all know, the back propagation learning algorithm is a common method for training artificial neural networks. Firstly, forward propagation of input vectors was realized to generate the propagation output activations, and then ANN uses the pattern target to generate the deviation of all output and hidden neurons. Based on this precondition, multiplying the calculation between the output deviation and the input activation was operated for obtaining the gradient vectors of those weigh values, which results into the percentage subtraction of their gradient vectors. Repeat these steps again until the performance of ANN network reaches a satisfied level. The three layers of the network architecture include the input layer, middle layer (hidden layer), and output layer. Layers include several processing units known as neurons. They connected with each other by variable weights to be determined. In this network, the input layer receives information from external source and passes this information to the network for data processing. The middle layer receives data from the input layer and does all information analysis. The output layer receives the processed information from the middle layer and sends the results to an external receptor [47].
A neuron in ANN network produces its input by processing the net input through nonlinear activation (transfer) function. The sigmoidal activation function is the most utilized one, which updates the weight and derivative values of ANN according to the resilient back propagation algorithm; therefore, it is usually trained for updating ANN networks. In this training algorithm, the update value for ANN weight is increased whenever the derivative value of its performance function has the same sign for two successive calculate iterations; on the other hand it is decreased according to the derivative value with respect to that weight changes sign from the previous iteration. If the derivative value is zero, then the update value for ANN weight remains the same. If the objective weight value continues to change in the same direction for several training iterations, then the magnitude of the weight change will be increased accordingly. In this experiment, we supposed that the input layer has NI neurons, with the middle layer and output layer having NJ and NK neurons, respectively, whereas the input elements for the th neuron at the middle layer can then be described as Here, denotes the weight calibrating from the th neuron on the input layer to the th neuron on the middle layer, with demonstrating the output result of the th neuron on the input layer. Meanwhile, the resultant input of the th neuron on the output layer can also be identified as net = ∑ =1 ; = 1, 2, . . . , .
Here, denotes the weight calibrating from the th neuron on the middle layer to the th neuron on the output layer, with demonstrating the output result of the th neuron on the middle layer. Through substituting the input neurons of the middle and output layer, as denoted by net Here, and denote the output threshold value used for assessing the computation results of the th neuron on the middle layer and the th neuron on the output layer, respectively. In this experiment, the gradient-based learning rule was used to train this established network, with the minimum mean square error (MSE) between the computation results caused by the network output and tablet template being predetermined as the ultimate train objective. The template amount for training was supposed to be , and the input elements can then be denoted by 1 , 2 , . . . , , with the output elements being denoted by 1 , 2 , . . . , . Here, this experiment used 1 , 2 , . . . , to denote the data vector consisted by training samples. Thus, MSE of the th sample can be obtained as follows: Here, denotes the reference value of the th sample for the th output neutron, while denotes the practical output in the same condition. Therefore, the weight on the output layer can then be adjusted to Here, represents the learning rate, denotes the momentum factor, and denotes the deviation value between and , which shows the difference between the reference value and the practical output of the th sample  caused by the th output neutron [48]. Therefore, the weights of middle layer can then be adjusted into When focusing on the ANN Testing process, it should be noted that the idea which we have used for testing neural networks and generally applying test vectors is that we use a multicase test vector to determine the decision strength of this neural network. For the final evaluation, we captured 100 images describing the random movement cases of drug tablets. Such a large testing case amount was used to minimize the uncertainty of the objective classifier's performance estimation in this experiment. Then, according to the reflective energy surfaces of those drug tablets, we input their corresponding property vectors to the established ANN and check the obtained output vectors in sequence. Based on the testing operations, this neural network reaches the desired performance; thus, a high improvement in the video detection of drug tablets can be observed.
The weight values in ANN network markedly affect the classification of tablet templates, which keep a close correlation with the computation or tracking errors in the proposed experimental conditions; they were also easily impacted by the arrangement of neural neutrons and input/output vectors as well as weight values provide correlation among network layers, which ensures the feasibility of recognizing the geometrical shapes of tablet template and the property distribution of its reflective energy surface with the selfadaptive control of ANN. As Figure 7 shows, the detailed setup of ANN network is used, and it was observed that, for the purpose of realizing the improvements in network recognition and objective classification, the weights between input and middle layers should be kept in a relatively stable state, which provides a useful tool to markedly eliminate the external error-training interference caused by noisy sample vectors or error distribution gradient of input data. Simultaneously, the weights between middle and output layers should be kept in a radically distinction state, which ensures a clearer classification of tablet properties, especially when the selfadaptive learning rule was being employed. Based on these preconditions, ANN network was gradually optimized for improving its recognition performance [49,50].

Tracking Process Discussions and Performance Comparisons
With the trained neural network and the surface properties computed by (1)-(6), they were used as the input vector for recognizing the tablet templates: Then, the output vector can be predetermined as [ 1 , 2 , 3 , 4 ] by using forward propagation algorithm; the maximum possible reiteration times were determined as 15000, in order to prevent the case of training threshold does not be met; the node number of the hidden layer was 2, the learning rate coefficient was supposed to be 0.45, the momentum factor was supposed to be 0.065, the train step was supposed to be 0.15, and the interval illustration factor for recognition process was predetermined as 30. The error function was determined by (27), and it can also be employed to describe the ANN performance capability: This index examines if the prerequirement was satisfied or not: all < ? If it was confirmed, then the training processes of ANN can be terminated or it will be iterated for the next time until this prerequirement is satisfied.
When discussing the determination of , through establishing the weights between network layers, then  Here, denotes a distinguishing coefficient located in [0, 1]. Δ denotes the self-adaptive deviation value between the weight values calibrated from the th neuron on the middle layer to the th neuron on the output; denotes the deviation value between and . By using this threshold screening mechanism, the accuracy of pattern recognition can be ensured, and thereafter the overtraining phenomenon in ANN network could also be prevented. The greater this threshold value is, clearer the objective recognition result would be. Table 1 demonstrates the original and corrected weights of ANN. It is noteworthy that the original weights were tentatively selected with reference to (20)- (21), in the hope of reducing the computation deviation between the desired and practical output vectors to the minimum scale [51,52]. Table 2 demonstrates the energy surface properties of tablet templates on the 35th image frame. In the interest of clearly describing the detailed process of tablet template recognition and its position justification, a program flow schematic was provided, with the details shown in Figure 8. Table 3 Does the output vector match the predetermined vector result at the current x-axis position?
The next matching position in x-axis direction  Predetermined vector data comparison, the obtained results keep a close difference with those predetermined one, which undoubtedly confirms their effectiveness and accuracy. As one typical pattern recognition, Figure 9 shows the coordinate positions of surface areas best matching their corresponding tablet templates on the 35th image frame, with the black boxes showing the probable positions of the targeted tablets to be studied. In this figure, the most possible matching positions have been magnified and displayed in surrounding images for clearer description; the coordinate positions of template centers were highlighted by the focus marks. For better illustration of the tablet tracking result, 10 key tablets have been labeled with numbers, especially those close ones that were highly focused, as Figure 10 shows. We can observe that the practical operation with this algorithm can be smoothly identified, and those tablets in close position can also be distinguished clearly. Besides, as the interval distance between two adjacent dashed lines was predetermined as 30 mm, the instantaneous moving speeds of objective tablets can be calculated. In order to prove the accuracy and validity of this newly proposed method, several commonly used tracking methods were investigated with reference to the published literatures, including the maximum likelihood estimation, the Bayesian estimation, the prior probability-density statistics, the nonprior information estimation, the invariant prior probabilitydensity statistics, the Jeffreys statistics, the box-tracing-based method, and the sensor tracing method. Typical evaluation indexes were employed when an identical experimental condition has been considered: such as the computation time (the accurate time consumption in which the whole experiment uses the identical tracing platform; this research measures the computation time by means of the automated time counting instrument), the distributed Hash degree (the decentralized distributed complexity level that provides a lookup assessment index for algorithm evaluation similar to traditional Hash table, with a small value associated with the highefficiently structure of tracing algorithm, and vice versa. This research calculates them through the performance process of one given algorithm and then makes a definite assessment on its performance capabilities as respected), the tracking error (the tracking ratio between the amount of error-tracking objectives to total ones, by which the tracing accuracy can be quantitatively identified. It can be calibrated by counting the amount of wrong-located tablets and total ones, resp., and then the ratio between them can be used as the tracking error), and the computation storage (the internal memory  Template 2   Template 4  Template 5  Template 6   Template 7  Template 8  Template 12 Template 11 x-axis y-axis O Figure 9: The position distributions of surface areas to best match their corresponding tablet templates on the 35th image frame. capacity used for the whole experimental computations; we calibrate them from the real-time performance indicator of the video tracking system that is being involved during the whole process of tablet tracking) [44][45][46][47][48]. After determining these evaluation indexes, their average performance capabilities were compared to each other, with the results shown in Table 4 after data equalization and value regularization. Since this experiment is repeatable, we take a detailed notes about the practical performances for 10, 20, 30, 40, and 50 times, and thereafter a statistical evaluation on their mean index values can be made by employing those mentioned algorithms, respectively. It can be seen from Figure 11 that, with the increment of experiment times, most of these tracking algorithms show a growth tendency in the average computation times. But it is noteworthy that a relativestable performance capability can be ensured by using our newly proposed algorithm, which is kept at about 0.443 s throughout the whole process rather than changing radically. Similar performance evaluations can also be observed in the variation tendency schematics of Figure 12 with the average Hash degree (from 59.8% to 5.66%), the average tracking error (from 2.3% to 2.445%) in Figure 13, and the average computation storage (from 797.5 kb to 804.5 kb) in Figure 14.
From Table 4, it can be observed that the maximum likelihood estimation or the sensor tracing method has a relative superiority in the tracking error and distributed Hash degree; they can be widely used to investigate a stable statistical condition, by concerning the dynamic moving characteristics of identical tablets in huge cluster; the Bayesian estimation obtains a good performance result in the computation time; therefore, it will be more suitable to fast assess a simpler mutual relationship mechanism for the approximate shape properties of those drug tablets; the prior probabilitydensity statistics method or the Jeffreys statistics approach have excellent computation performances when the tracking error or distributed Hash degree were highly emphasized, which ensure the accurate demonstrations of data relationship analysis and thereafter make a series of remarkable progresses when compared with other traditional methods; the nonprior information estimation method and the invariant prior probability-density statistics approach have extraordinary capabilities in the distributed Hash degree, which show a more precise algorithm complexity of data distribution in multidimensional statistical domain. Finally a good performance evaluation result can be obtained by using this newly proposed algorithm, especially in the computation time, computation storage, or distributed Hash degree. Performance comparison proves its validation and efficiency when similar shape or identical size of objective tablets should be taken into account, and the working estimation on how well these alternative methods perform on average can be facilitated as well.

Conclusions
A new image tracking approach for high similarity drug tablets based on light intensity reflective energy and ANN recognition has been investigated. Considering the practical condition and precision requirement, the reflective energy  surfaces were established for modeling the high similarity geometric characteristics and surface topography of those targeted drug tablets, which results into the maximum deletion of external optical signal interferences or calibration error in practice, and simultaneously an accurate description of tablet objectives in the light intensity domain can also be ensured. Thereafter, ANN was used for recognizing those studied objective tablets with the computed surface properties, and then their instantaneous positions on one image frame can be identified clearly. Through repeating these steps on the sequential image frames, a series of tablet tracking results can be obtained. After being compared   with other tracking methods in performance indexes and determined results, it can be learned that this newly proposed method solves numerous difficult problems characterized by monitoring precision, complicated calculations or tracking errors, and other signal interferences caused by mathematical calibration or pattern recognition. Therefore, the quantitative analysis of the objective tracking can be simplified, and the searching precision or computation efficiency can also be greatly improved, accordingly. Experimental computation and result analysis confirmed the validation and accuracy of this new method; new research ideas for real-time objective recognition or tablet tracking can be provided as well.