Printer Friendly

Supervised expert system for wearable MEMS accelerometer-based fall detector.

1. Introduction

The problem of falls in the elderly has become a health care priority due to the related high social and economic costs [1]. In fact the European population aged 65 years or more, which may be in need of assistance is increasing. This trend asks care-holders institutions to employ more efficient and optimized methods in order to be able to grant the required service at lower costs. The consequences of falls in the elderly may lead to psychological trauma, physical injuries, hospitalization, and even death in the worst scenario [2-5]. The main reason that pushed for the development of the presented system is to allow noncompletely self-sufficient people (e.g., older people) to live safely in their own houses as long as possible. This is important not only for aspects of health regarding assisted people, but also for the consequent social advantages. The European community issued and funded various projects and consortia. The mission focuses on several purposes, all addressed to older people, varying from the assistance in case of need, to the prevention of dangerous or unhealthy situations. The purpose of the work described in this paper is to focus on people fall detection.

Many solutions have been proposed in the detection and prevention of falls, and some excellent review studies were presented [1, 6]. Basically, fall-detection solutions can be classified in three main classes: wearable devices, ambient devices, and camera-based devices. The first approach requires that the elderly holds some kind of devices (e.g., an assistive cane) or wears sensors like accelerometers and/or gyroscopes to detect the motion of the body. In particular, recent miniaturization and cost reduction of MEMS accelerometers and the availability of reliable wireless communication technologies enabled the realization of affordable wearable monitoring systems that can be worn by people performing their normal daily activities [7-10]. For these reasons, in the last few years, the use of portable devices in the health monitoring of chronic patients has increased considerably. However, these devices have some drawbacks: they are prone to be forgotten, worn in a wrong body position, or accidentally damaged. Regarding fall detection, with respect to vision or acoustic sensors, the accelerometer module has the advantage of not having to be set up and installed in all rooms of the "smart home," as it is required for instance for 3D video trackers or acoustic scene analyzers. On the other hand, the camera-based approach is a nonintrusive solution, since the sensor is a camera installed on a wall and is able to detect falls, switching off when the monitored person is no longer alone [11]. On the other hand, camera-based monitoring solutions suffer when large occluding objects (i.e., furniture) obstruct camera's viewing. Another significant drawback of any camera-based solution is the violation of monitored person's privacy.

In this paper a fall detector through a wearable triaxial MEMS accelerometer is presented. The proposed solution overcomes the limitation of well-known threshold-based approaches [12-17] for which the accelerometer-based device need, to be tuned in an appropriate way for each installation (i.e., the parameters setup may be different for different people). For these reasons, a machine-learning scheme [18] is used, showing high generalization capabilities in the fall detection discrimination process. The expert system uses robust features extracted taking into account important constraints and requirements of mobile solutions (workload). The extracted features are (quasi-) invariant both to specific characteristics of the mounting setup (device on chest, on waist, and on abdomen) and specific characteristics of the end users in terms of age, weight, height, and gender.

2. Materials and Methods

2.1. Triaxial Accelerometer Sensor System. The hardware used is a wearable device composed by commercial discrete circuits, according to the design proposed in [19]. The picture and the logical block diagram are shown in Figure 1. The system integrates an ST LIS3LV02DL tri-axial MEMS accelerometer with digital output, an FPGA for computing functionalities, and a ZigBee module for wireless communication up to 30 m. The power consumption is about 190 mW in streaming mode and 9 mW in idle.

The wearable device can operate in streaming (raw data are sent via ZigBee to an external computing platform for data analysis with a 10 Hz frequency) or in standalone modality, by running the threshold-based fall detection implementation on the on-board FPGA (the power consumption is limited since the ZigBee module is activated just when a fall event occurs). The LIS3LV02DL MEMS accelerometer is DC coupled, and it responds up to 0 Hz, with 16-bits resolution and a full scale in the range [+ or -] 2g. Data can be transmitted in hexadecimal format. The sensor measures both static and dynamic accelerations along the 3 axes and allows one to receive information on the 3D spatial relative position (compared to the Earth gravity vector) of the person who wears it.

In stationary conditions, assuming a particular axis, the component of the acceleration (amplitude A, rif. (1)) is defined according to the value of the sine of the angle a between the considered axis and the horizontal plane, which is perpendicular to the Earth gravity component (g):

A = g sin ([alpha]). (1)

In this way, if the accelerometer relative orientation is known, the resulting data can be used to determinate the angle of the user posture respect to the vertical direction.

2.2. Simulated Falls and ADL Tasks. In order to analyse the waveforms along each axis in the presence of falls and other kind of events (Activities of Daily Living, ADLs), a data collection has been defined in controlled (simulated) conditions by involving 11 healthy male and 2 healthy female volunteers. The simulated falls were performed by using a crash mat (height 2 cm) and knee/elbow pad protectors, meeting safety and ethical requirements. The range of actors age was 39.3 [+ or -] 12.3 years, weighing 73.7 [+ or -] 13.4 kg, and a height of 1.76 [+ or -] 0.1 m. 450 actions were simulated in which 250 were falls compliant to the specifications proposed by Noury et al. [20]. The following falls have been simulated for study:

(1) backward falls ending in the lying position,

(2) backward falls with recovery,

(3) forward falls ending in lying flat,

(4) forward falls with recovery,

(5) lateral falls.

During the data collection the wearable device was placed with an elastic band in a different positioning area on the upper part of the torso (on the chest, on the waist, and on the abdomen). Some ADLs were simulated in order to evaluate the ability of fall-detection algorithms to discriminate falls from ADLs. The simulated ADL tasks belong to the following categories:

(1) walking,

(2) sitting down on a chair and then standing up,

(3) lying down on a mat (height 33 cm) and then standing up,

(4) lying down on a mat (height 2 cm) and then standing up,

(5) kneeling on a mat (height 2 cm) and then standing up.

Each actor performed more than 15 simulated ADLs for a total of 205.

2.3. Preprocessing and System Calibration. The first four computational steps (Figure 2) deal with data acquisition, data preprocessing, noise filtering, and system calibration.

The acceleration data on three axes (Ax, Ay, and Az) is read out from the device worn by a user during the data collection. Data was stored into a portable computer and converted into gravitational units to represent acceleration data in the range [+ or -] 2g, in order to make it possible to extract the angle [alpha] (described previously) and for not having orders of magnitude too different in the features (during the training of the classifier used to detect the fall events). The samples coming from the device are filtered out by a low pass 8 order, 8 Hz cut-off FIR (Finite Impulse Response) filter to reduce the noise due to electronic components, environment, and human tremor.

In order to correctly handle preprocessed data, a calibration procedure was accomplished by recovering the initial conditions after the device mounting. During this step the correct placing of hardware is verified by checking if two acceleration axes are orthogonal to the Earth gravity g (Figure 3): the acceleration values measured on the two orthogonal components must be close to zero.

The calibration procedure is composed by the following steps:

(i) the user wears the device in a standing position for 10 seconds;

(ii) the calibration routine calculates the average of the acceleration on each axis over this period. These are the initial acceleration values [A.sub.x0], [A.sub.y0], and [A.sub.z0];

(iii) if [A.sub.x0], [A.sub.y0], and [A.sub.z0] will be close to those expected, they will be considered as references in the fall detector algorithm, otherwise a routine to compensate the sensed misplacement will be enabled.

With respect to Figure 3, the acceleration along x-component will be close to -1 and almost zero for the others according to (1). The expected values of acceleration on three axes [A.sub.x0], [A.sub.y0], and [A.sub.z0] are also reported in Table 1.

In practice, it is difficult to locate the device exactly in the right position, and the measured acceleration values will differ slightly from those expected. The calibration procedure is concluded if the following conditions are satisfied at the same time (and the values [A.sub.x0], [A.sub.y0], and [A.sub.z0] will be recorded and used as references in the feature extraction phase):


Since the values ([A.sup.ref.sub.x0], [A.sup.ref.sub.y0] and [A.sup.ref.sub.z0]) are (-1, 0, and 0);

[absolute value of [A.sub.x0] + l] < 0.3, [absolute value of [A.sub.y0]] < 0.3, [absolute value of [A.sub.z0]] < 0.3. (3)

Otherwise a routine to compensate the sensed misplacement is enabled, and the angles displacements of the sensor axes ([alpha][A.sub.x0], [alpha][A.sub.y0], and [alpha][A.sub.z0]) are calculated using the following trigonometric equation:


These values are stored and will be used for the correction of the misplacement during the feature extraction phase, described below.

2.4. Feature Extraction. Robust features are extracted in the time domain in a 5 sec sliding window by considering both quick and relevant acceleration changing along each axis (due to the fall) and by the change in position registered after the fall. The aim is to produce robust features taking into account all the information for the distinction of falls from other events. It is also important that such features have a low dependence on both the position of the sensor (whether it is placed on the waist, on the chest, or on the abdomen) and the human body characteristics of the user. Moreover, the computational cost must be limited for integration on embedded computing solutions.

The waveform in Figure 4 represents a typical acceleration signal of a forward fall on three axes and all phases are indicated according to the taxonomy proposed in [16]. For the features extraction process, both critical and postfall phases are of interest. In the former, the shock is measured due to the impact toward the ground, and a dynamic acceleration changing is registered. In the latter (the body is already lying on the floor) the static acceleration value records a great change due to the new position of the individual with respect to the calibration phase.

The 5 sec sliding window considered for features extraction is split in three parts as follow: 1 2

(1) from 0 sec to 0.5 sec: the maximum and minimum values of the acceleration are found, and the amplitude dynamics [DELTA]A is calculated. It is proportional to the shock strength;

(2) from 0.5 sec to 2 sec: the system is in a transitory regime (samples are discarded);

(3) from 2 sec to 5 sec: the static-averaged acceleration value [bar.A] is calculated to evaluate the new user position. Acceleration data is averaged to filter out tremors and the little movements of the end user.

During the system calibration phase it is verified that the device is worn correctly (see previous section), and the position changing is calculated as the difference between the value of the 3D-static acceleration after the fall (in the third part of the sliding window) and the one stored in the calibration phase. This difference, called Changing Position Offset (CPO), has a value in the range of 0 to 2, and it is proportional to the user displacement. In this way, for feature extraction, a study of posture was not made: only the relative varying posture analysis was considered causing a computational cost reduction and improving the robustness of the setup. The feature vector is made up of three parameters (one for each axis), coming from the modulation of the CPO and the dynamic acceleration peak, due to the impact of the fall.

If the routine to compensate the device misplacement is enabled, for the CPO calculation, the axes angles of the sensor in the initial condition need to be taken into account ([alpha][A.sub.x0], [alpha]x[A.sub.y0], and [alpha][A.sub.z0], stored during the calibration phase). Using (1), the CPO value for x-axis is obtained as

[CPO.sub.x] = [absolute value of sin ([alpha][[bar.A].sub.x] - [alpha][A.sub.x0])], (5)

where [alpha][[bar.A].sub.x] is the x-axis angle of the sensor in the third part of the sliding window; it is calculated considering the static averaged acceleration [[bar.A].sub.x] and the same trigonometric equation used in (4) as follows:


The same procedure must be done for the two other axes.

Pseudocode 1: Feature extraction pseudocode.

Let [A.sub.0] initial acceleration value measured
after the calibration phase

Let [A.sub.0] initial axis angle of the sensor measured
after the calibration phase

[DELTA]A [right arrow] shock dynamics due to the impact
calculated in the first part of the sliding window

A [right arrow] Static averaged acceleration in the third
part of the sliding window

If compensation device misplacement procedure is enabled

[alpha][bar.A] is calculated [right arrow] axis angle of
the sensor in the third part of the sliding window

CPO = [absolute value of sin ([alpha][bar.A] -
[alpha][A.sub.0])] [right arrow] Changing Position Offset


CPO = [absolute value of [bar.A] - [A.sub.0]] [right arrow]
Changing Position Offset

End if

Let [DELTA][A.sub.x], [DELTA][A.sub.y], [DELTA][A.sub.z], [CPO.sub.x],
[CPO.sub.y], [CPO.sub.z], [DELTA]A and CPO calculated for each axis

[X, Y, Z] = [[DELTA][A.sub.x], [CPO.sub.x], [DELTA][A.sub.y],
[CPO.sub.y], [DELTA][A.sub.z], [COP.sub.z] [right arrow] feature
vector input for SVM

Thanks to this routine, the change in the angle position of the user along the three axes is measured. In this way the information coming from the modulation of the CPO remains highly selective for fall detection recognition even if the device is placed incorrectly and the features remain nearly the same. In this way the system could work efficiently reducing the problems of device positioning at the expense of an increase in the computational costs.

The process of feature extraction, described previous, is summarized in the pseudocode shown in Pseudocode 1.

For the feature extraction, It makes sense to consider the acceleration signal on each axis singularly, because a fall event leads to a change in the value of the static acceleration in at least two of the three acceleration axes (due to the orientation change of sensing axes). This is evident comparing the positions A, B, and C in Figure 5 (postfall phase of the lateral, backward, and forward falls, resp., are depicted) with respect to the initial position (see Figure 3). On the other hand, when a sitting event (D in Figure 5) occurs, the change in static acceleration can be neglected with respect to that of a fall; it is possible to do the same considerations for kneeling and walking (E and F). Instead, the axes orientation changes for a lying event (position G), but the acceleration peak produced is slower and lower than a fall, and through the product between the value of the acceleration peak and CPO it is also possible to discriminate this ADL.

Table 2 shows the nominal values of acceleration along the three axes [A.sub.x], [A.sub.y], and [A.sub.z] and the value of CPO corresponding to the positions A-G of Figure 5.

When the orientation of acceleration does not change, the values of CPO are (close to) zero, and the elements of the features also become (close to) zero. In Figures 6 and 7 the waveforms of sitting and lying events are reported to make a comparison with that of the fall (shown in Figure 4), when the device is worn on the abdomen.

The difference among the features of the fall of Figure 4 and the sitting and lying events of Figures 6 and 7 is apparent in Figures 8, 9, and 10.

It was actually verified that all the consideration made thus so far remains valid whether the device is placed on the waist or on the chest. In Figures 11 and 12 the features of a sequence of falls and daily events, when the device is worn on the waist and on the chest, are shown. So it is also evident that the features obtained discriminate the falls from ADLs when the device is placed in other area of the torso.

It is important to highlight that the measurements of the actions as shown in Figures 11 and 12 have been carried out by simulating critical behaviours, that is, sitting with a strong impact on the chair; bending down, lying down, and standing up quickly; moderate backward fall.

2.5. Fall Detection Algorithm. Once features are extracted, the fall events are detected by a one-class support vector machine (OC-SVM). SVM is a robust classification tool (in the presence of outliers too) with a good generalization ability. Furthermore it is less computationally intensive than other algorithms like neural networks [21].

One class SVM divides all samples into an objective field and a nonobjective field and then nonlinearly maps those sample into a higher dimensional features space with some efficient operators, called kernel function. The target is to find a sphere that contains most of the normal data such that the corresponding radius R can be minimized as follows:


where a is the centre of the sphere and [v.sub.i] a positive sample set. The slack variables [[xi].sub.i] allow some data points to lie outside the sphere and the parameter C controls the trade-off between the volume of the sphere and the number of errors. The objective function is


The original data points are first mapped into a feature space, because some data are not spherically distributed in the input space. For this reason the kernel function k(., .) is used as follows:


The common types of kernel found in the literature are tested (in Table 3 the functions are shown). They were analyzed using the MATLAB tool STPRtool [22], changing the values of their parameters.

In the literature the Gaussian Radial Basis Function (GRBF) is often used in fall detectors [23-25]. For the features extracted, GRBF and polynomial kernel functions give the best results in terms of performance, even the polynomial kernel shows a lower computational cost (relative number of support vectors and relative execution time are considered). The comparison with other kernel functions will be shown in the following section.

The parameters of kernel were chosen to detect as many falls as possible, and a simple post-processing step has been implemented in order to reduce false negative events.

2.6. Postprocessing. In choosing the parameters of the tested kernel functions, the value of sensitivity was given priority. In this way some nonfall events are detected from SVM. To reduce this problem a simple filtering by voting was added: it detects an event such as a fall if the alarm coming from the SVM classifier remains high a given time t. With this approach the spikes or temporary anomalies are filtered too. The optimum time t for postprocessing is obtained from a parametric analysis. The results are shown in the next paragraph.

3. Results

This section will show a discussion about the measurements of the experimental data processed in MATLAB. To validate the implemented algorithm, the dataset previously described has been considered.

Since the output of the fall detection step is binary, there are four possible results.

(i) True Positive (TP): a fall happens, and the algorithm detects it.

(ii) False Positive (FP): a fall does not occur, and the algorithm reveals a fall.

(iii) True Negative (TN): a daily event is performed, and the algorithm does not detect it.

(iv) False Negative (FN): a fall occurs, but the algorithm does not detect it.

The performance of the system, with respect to both different kernel functions of SVM and their parameters, can be evaluated considering the following metrics:

sensitivity = TP/TP + FN, (10)

specificity = TN/TN + FP.

For completeness the results obtained with postprocessing elaboration have been added. The algorithm is tested when the device is placed on the waist, abdomen, or chest.

The one class SVM classifier has been trained by using about 40 falls and 50 ADLs belonging to a large dataset in which more than 250 falls and 200 daily events were performed. The remaining 210 falls and 150 ADLs have been used for testing. Table 4 shows the specificity and the sensitivity for the kernel functions already described. The complexity of the various kernel functions can be studied by considering the relative number of vectors required and also the relative execution time [26,27]. The kernel with the most significant values for the parameters is considered.

From Table 4 the polynomial kernel and GRBF kernel give better results. Their capacity to detect a fall is higher (more than 95% for sensitivity and specificity) than others. GRBF (with a = 3 and C = 2) and the polynomial (with P = 3 and C = 2.8) provide similar values of specificity and sensitivity, but the last one works faster, and its number of vector is slightly lower, thereby it was chosen for this work. Misclassifications are for falls presenting slow dynamics or falls with partial recovery. It is important to see the result of the linear function: its workload is very low (less support vectors number and very low execution time), but it gives many false positive, and its value of the specificity is low. The function of GRBF is the most common used in the literature for the detection of falls, due to its high ability to separate the classes, even in the presence of many outliers (at the expense of the convergence time). In this workit has been also possible to obtain a high accuracy with the polynomial function thanks to features extracted: they have a high degree of intrinsic separability.

The implemented OC-SVM shows improvements in the specificity and sensitivity with respect to a threshold-based approach, and this can be verified in comparison with both the frameworks presented in [17,28], choosing the algorithms where the same parameters of fall are used (the impact detection and posture monitoring). For the comparison of the fall detection algorithms were used the same hardware (shown in Figure 1), benchmark dataset and training/test sets described above. The values of thresholds were tuned according to the falls presented in the dataset used to train the OC-SVM scheme. The best results for threshold based algorithm were obtained with Bourke's algorithm and are compared with the proposed OC-SVM in Table 5: the latter seems more efficient, demonstrating its higher capacity to generalize, detecting a bigger number of falls with low impact magnitude than threshold-based algorithm. However, the computational cost of thresholdbased is lower.

Due to (a) the reduction in the number of the features, (b) the low computational cost of extracted features, and (c) the used kernel, the overall system workload implemented is compatible with an integration in embedded low-power solutions (DSP, FPGA, and microcontroller).

A reduction of false positive (and so an improving of specificity value), for the polynomial chosen, has been measured with filtering by voting with a temporal window of prefixed dimension (0.8 sec); see Table 6. Finally it was verified that the values of sensitivity and specificity remain almost the same even if the axes of the sensor are placed differently with respect to that shown in Figure 3 (the values of testing dataset have been changed simulating same misplacements of the device).

4. Conclusions

The proposed supervised scheme overcomes the limitation of well-known threshold-based approaches in which a heuristic choice of the parameters is accomplished. High performance in controlled conditions (events simulated) in terms of sensitivity and specificity was obtained using only the 20% of dataset for training. Performance metrics of different kernels in one class SVM are compared: best results are obtained with polynomial function and Gaussian Radial Basis Function. The polynomial kernel is used in order to limit the computational workload. A study of posture was not made, but only posture changing analysis was used, limiting the computational cost. Through the calibration step, the approach allows to generalize the detection of fall events leading invariance to physical characteristic of the end users. Future work will be devoted to validate the solution in real conditions, test the methodology with a large set of different MEMS accelerometers, and port the framework on embedded mobile solutions as FPGA or DSP.


[1] N. Noury, P. Rumeau, A. K. Bourke, G. OLaighin, and J. E. Lundy, "A proposal for the classification and evaluation of fall detectors" IRBM, vol. 29, no. 6, pp. 340-349, 2008.

[2] M. C. Chung, K. J. McKee, C. Austin et al., "Posttraumatic stress disorder in older people after a fall" International Journal of Geriatric Psychiatry, vol. 24, no. 9, pp. 955-964, 2009.

[3] S. Sadigh, A. Reimers, R. Andersson, and L. Laflamme, "Falls and fall-related injuries among the elderly: a survey of residential-care facilities in a Swedish municipality," Journal of Community Health, vol. 29, no. 2, pp. 129-140, 2004.

[4] A. Shumway-Cook, A. M. Ciol, J. Hoffman, J. B. Dudgeon, K. Yorkston, and L. Chan, "Falls in the medicare population: incidence, associated factors, and on health Care," Physical Therapy, vol. 89, no. 4, pp. 324-332, 2009.

[5] S. Elliott, J. Painter, and S. Hudson, "Living alone and fall risk factors in community-dwelling middle age and older adults," Journal of Community Health, vol. 34, no. 4, pp. 301-310, 2009.

[6] X. Yu, "Approaches and principles of fall detection for elderly and patient," in Proceedings of the 10th IEEE International Conference on E-Health Networking, Applications and Service (HEALTHCOM '08), pp. 42-47, Singapore, July 2008.

[7] D. M. Karantonis, M. R. Narayanan, M. Mathie, N. H. Lovell, and B. G. Celler, "Implementation of a real-time human movement classifier using a triaxial accelerometer for ambulatory monitoring," IEEE Transactions on Information Technology in Biomedicine, vol. 10, no. 1, pp. 156-167, 2006.

[8] J. Y. Hwang, J. M. Kang, Y. W. Jang, and H. C. Kim, "Development of novel algorithm and real-time monitoring ambulatory system using Bluetooth module for fall detection in the elderly," in Proceedings of the 26th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '04), pp. 2204-2207, San Francisco, Calif, USA, September 2004.

[9] H. J. Luinge and P. H. Veltink, "Inclination measurement of human movement using a 3-D accelerometer with autocalibration," IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 12, no. 1, pp. 112-121, 2004.

[10] G. Anania, A. Tognetti, N. Carbonaro et al., "Development of a novel algorithm for human fall detection using wearable sensors," in Proceedings of the IEEE Sensors (SENSORS '08), pp. 1336-1339, Lecce, Italy, October 2008.

[11] A. Leone, G. Diraco, and P. Siciliano, "Detecting falls with 3D range camera in ambient assisted living applications: a preliminary study," Medical Engineering and Physics, vol. 33, no. 6, pp. 770-781, 2011.

[12] J. Chen, K. Kwong, D. Chang, J. Luk, and R. Bajcsy, "Wearable sensors for reliable fall detection," in Proceedings of the 27th Annual International Conference of the Engineering in Medicine and Biology Society (IEEE-EMBS '05), pp. 3551-3554, Shanghai, China, September 2005.

[13] M. Kangas, A. Konttila, P. Lindgren, I. Winblad, and T. Jamsa, "Comparison of low-complexity fall detection algorithms for body attached accelerometers," Gait & Posture, vol. 28, no. 2, pp. 285-291, 2008.

[14] A. K. Bourke, J. V. O'Brien, and G. M. Lyons, "Evaluation of a threshold-based tri-axial accelerometer fall detection algorithm," Gait & Posture, vol. 26, no. 2, pp. 194-199, 2007.

[15] A. K. Bourke, P. Van de Ven, A. E. Chaya, G. M. Olaighin, and J. Nelson, "Testing of a long-term fall detection system incorporated into a custom vest for the elderly," in Proceedings of the 30th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC '08), pp. 2844-2847, Vancouver, Canada, August 2008.

[16] A. K. Bourke, P. van de Ven, M. Gamble et al., "Evaluation of waist-mounted tri-axial accelerometer based fall-detection algorithms during scripted and continuous unscripted activities," Journal of Biomechanics, vol. 43, no. 15, pp. 3051-3057, 2010.

[17] F. Bagala, C. Becker, A. Cappello et al., "Evaluation of accelerometer-based fall detection algorithms on real-world falls," PLoS One, vol. 7, Article ID e37062, 2012.

[18] S. -H. Liu and W. -C. Cheng, "Fall detection with the support vector machine during scripted and continuous unscripted activities," Sensor, vol. 12, no. 9, pp. 12301-12316, 2012.

[19] A. Leone, G. Diraco, C. Distante et al., "A multi-sensor approach for people fall detection in home environment," in Proceedings ofthe Workshop on Multi-Camera and Multi-Modal Sensor Fusion Algorithms and Applications (M2SFA2 '08), pp. 1-12, 2008.

[20] N. Noury, A. Fleury, P. Rumeau et al., "Fall detection--principles and methods," in Proceedings of the 29th Annual International Conference of IEEE-EMBS, Engineering in Medicine and Biology Society (EMBC '07), pp. 1663-1666, Lyon, France, August 2007.

[21] L. M. Manevitz and M. Yousef, "One-class SVMs for document classification," Journal of Machine Learning Research, vol. 2, no. 1, pp. 139-154, 2001.

[22] V. Franc and V. Hlavac, "Statistical Pattern Recognition Toolbox (STPRtool)".

[23] T. Zhang, J. Wang, L. Xu, and P. Liu, "Fall detection by wearable sensor and one-class SVM algorithm," Lecture Notes in Control and Information Sciences, vol. 345, pp. 858-863, 2006.

[24] J. Yin, Q. Yang, and J. J. Pan, "Sensor-based abnormal human-activity detection," IEEE Transactions on Knowledge and Data Engineering, vol. 20, no. 8, pp. 1082-1090, 2008.

[25] H. Cheng, H. Luo, and F. Zhao, "A fall detection algorithm based on pattern recognition and human posture analysis," in Proceedings of the IET International Conference on Communication Technology and Application (ICCTA '11), vol. 2011, pp. 853-857, Beijing, China.

[26] R. Sangeetha and B. Kalpana, "A comparative study and choice of an appropriate kernel for support vector machines," in Information and Communication Technologies, V. V. Das and R. Vijaykumar, Eds., vol. 101 of Communications in Computer and Information Science, pp. 549-553, Springer, Heidelberg, Germany, 2010.

[27] R. Sangeetha and B. Kalpana, "Performance evaluation of kernels in multiclass support vector machines," International Journal of Soft Computing and Engineering, vol. 1, no. 5, pp. 138-145, 2011.

[28] M. Grassi, A. Lombardi, G. Rescio et al., "A hardware-software framework for high-reliability people fall detection," in Proceedings of the IEEE Sensors (SENSORS '08), pp. 1328-1331, Lecce, Italy, October 2009.

Gabriele Rescio, Alessandro Leone, and Pietro Siciliano

Institute for Microelectronics and Microsystems, Italian National Research Council (CNR), Via Monteroni, c/o Campus Universita del Salento, Palazzina A3, 73100 Lecce, Italy

Correspondence should be addressed to Gabriele Rescio;

Received 8 February 2013; Revised 10 June 2013; Accepted 11 June 2013

Academic Editor: Andrea Cusano

TABLE 1: Acceleration values along the three axes
during the initial position of Figure 3.

           [A.sup.ref.   [A.sup.ref.   [A.sup.ref.
             sub.x0]       sub.y0]       sub.z0]

Initial        -1             0             0

TABLE 2: CPO coefficients by varying spatial position.

                    [A.sub.x]   [A.sub.y]   [A.sub.z]      CPO

Initial condition      -1           0           0
Position A              0          -1           0       (1, 1, 0)
Position B              0           0          -1       (1, 0, 1)
Position C              0           0           1       (1, 0, 1)
Position D             -1           0           0       (0, 0, 0)
Position E             -1           0           0       (0, 0, 0)
Position F             -1           0           0       (0, 0, 0)
Position G              0           0          -1       (1, 0, 1)

TABLE 3: Types of tested kernels.

Kernel       Function

Linear       K([v.sub.i], [v.sub.j]) =
               1 + [v.sup.T.sub.j][v.sub.j]
Polynomial   K([v.sub.i], [v.sub.j]) =
               [(1 + [v.sup.T.sub.j]
Gaussian     K([v.sub.i, [v.sub.j]) = exp
  Radial       (-[gamma][parallel][v.sub.i] -
  Basis        [v.sub.j][[parallel].sup.2])
Sigmoid      K([v.sub.i], [v.sub.j]) =
               tanh (k[v.sup.T.sub.j]
               [v.sub.j] - [delta])

TABLE 4: Performance by varying the kernel functions.

                     Sensitivity   Specificity   Relative   Relative
                         (%)           (%)       support    execution
                                                 vectors      time

Linear (C = 10)         98.6          52.3        0.51x       0.2x
Linear (C = 100)        97.2          46.4        0.44x       0.2x
Polynomial              97.5          81.25       0.94x       0.8x
  (C = 20, P = 2)
Polynomial              97.7          95.8          1x         1x
  (C = 2.8, P = 3)
GRBF (C = 2.6,          98.1          83.1         1.1x       1.3x
  [gamma] = 2)
GRBF (C = 2,            97.4          95.2         1.3x       1.5x
  [gamma] = 3)
Sigmoid (C = 1,         98.8           36          0.3x        4x
  k = 1,
  [delta] = 2)
Sigmoid (C = 100,       99.1           42          0.4x       4.1x
  k = 5,
  [delta] = 2)

TABLE 5: Comparison of the proposed OC-SVM method and
threshold-based algorithm.

                  Sensitivity (%)   Specificity (%)

OC-SVM                  977              94.8
Threshold based        89.2              85.7

TABLE 6: Performance after filtering by voting.

Polynomial         Sensitivity (%)   Specificity (%)
(C = 2.8, P = 3)

t = 0.4 s               97.7              94.8
t = 0.8 s               96.2              97.2
t = 1.2 s               86.4              98.1
COPYRIGHT 2013 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2013 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Rescio, Gabriele; Leone, Alessandro; Siciliano, Pietro
Publication:Journal of Sensors
Date:Jan 1, 2013
Previous Article:High-density fiber optical sensor and instrumentation for gas turbine operation condition monitoring.
Next Article:Force-sensor-based estimation of needle tip deflection in brachytherapy.

Terms of use | Privacy policy | Copyright © 2022 Farlex, Inc. | Feedback | For webmasters |