Electrooculogram Signals Analysis for Process Control Operator Based on Fuzzy c-means

Similar documents
n UL Listed and FM Approved for n Solenoid control n Quick pressure relief valve 73Q n Pressure sustaining & reducing valve 723

Load Calculation and Design of Roller Crowning of Truck Hub Bearing

SYMMETRY AND VARIABILITY OF VERTICAL GROUND REACTION FORCE AND CENTER OF PRESSURE IN ABLE-BODIED GAIT

ANALYSIS AND MODELING TIME HEADWAY DISTRIBUTIONS UNDER HEAVY TRAFFIC FLOW CONDITIONS IN THE URBAN HIGHWAYS: CASE OF ISFAHAN

SPEED OF SOUND MEASUREMENTS IN GAS-MIXTURES AT VARYING COMPOSITION USING AN ULTRASONIC GAS FLOW METER WITH SILICON BASED TRANSDUCERS

A Comparison of MOEA/D, NSGA II and SPEA2 Algorithms

THE EFFECTS OF COUPLED INDUCTORS IN PARALLEL INTERLEAVED BUCK CONVERTERS

P h o t o g r a p h i c L i g h t i n g ( 1 1 B )

THE LATENT DEMAND METHOD

ANALYSIS AND ASSESSMENT OF LITHUANIAN ROAD ACCIDENTS BY AHP METHOD

SPH4U Transmission of Waves in One and Two Dimensions LoRusso

Hazard Identificaiton of Railway Signaling System Using PHA and HAZOP Methods

The structure of the Fibonacci numbers in the modular ring Z 5

8.5. Solving Equations II. Goal Solve equations by balancing.

Influences of Teaching Styles and Motor Educability on Learning Outcomes of Volleyball

Chilled Mirror Dew Point Instrument

Real time lane departure warning system based on principal component analysis of grayscale distribution and risk evaluation model

Range St. Dev. n Mean. Total Mean % Competency. Range St. Dev. n Mean. Total Mean % Competency

number in a data set adds (or subtracts) that value to measures of center but does not affect measures of spread.

ASSESSMENT SCORING SYSTEM OF ROAD SAFETY INFRASTRUCTURE

Active Travel The Role of Self-Selection in Explaining the Effect of Built Environment on Active Travel

Catenary Analysis and Calculation Method of Track Rope of Cargo Cableway with Multiple Loads

Application of a Statistical Method to investigate patterns of beach evolution in the vicinity of a seawall

Available online at ScienceDirect. Procedia Engineering 113 (2015 )

DFC NIST DIGITAL MASS FLOW CONTROLLERS. DFC with optional LCD readout PROG RS485. Programmable Mass Flow Controller with Digital Signal Processing

Policy sensitivity analysis of Karachi commuters

Research Article. Relative analysis of Taekwondo back kick skills biomechanics based on 3D photograph parsing. Mingming Guo

A SECOND SOLUTION FOR THE RHIND PAPYRUS UNIT FRACTION DECOMPOSITIONS

A Comparative Investigation of Reheat In Gas Turbine Cycles

WIND TUNNEL EXPERIMENT ON THE EFFECT OF WIND ON SMOKE EXHAUST SYSTEMS FOR A HIGH RISE BUILDING

Controlling noise at work

ICC WORLD TWENTY ( WORLD CUP-2014 )- A CASE STUDY

The Analysis of Bullwhip Effect in Supply Chain Based on Strategic Alliance

The Prediction of Dynamic Strain in Leaf-Type Compressor Valves With Variable Mass and Stiffness

Traffic conflicts at roundabouts: risk analysis under car-following conditions

Seated valves (PN 16) VF 2-2-way valve, flange VF 3-3-way valve, flange

» WYOMING s RIDE 2013

University of California, Los Angeles Department of Statistics. Measures of central tendency and variation Data display

DEPARTMENT OF ELECTRICAL AND COMPUTER ENGINEERING, THE UNIVERSITY OF NEW MEXICO ECE-238L: Computer Logic Design Fall Notes - Chapter 6.

2) What s the Purpose of Your Project?

Basic Gas Spring Theory

Modelling Lane Changing Behaviour of Heavy Commercial Vehicles

Characterization of Refrigeration System Compressor Performance

Footwork is the foundation for a skilled basketball player, involving moves

The new name for... Mines Rescue Service

Obstacle Avoidance for Visually Impaired Using Auto-adaptive Thresholding on Kinect s Depth Image

XFM DIGITAL MASS FLOW METER. XFM with Profibus Interface. XFM without. Readout. XFM with. Readout. Option

Gait-Event-Based Synchronization Method for Gait Rehabilitation Robots via a Bio-inspired Adaptive Oscillator

M3P. Safety Data Sheet TABLE OF CONTENTS IDENTIFICATION OF THE SUBSTANCE/MIXTURE AND OF THE COMPANY/UNDERTAKING 2 SECTION 2 HAZARDS IDENTIFICATION 2

MINNESOTA DEER MANAGEMENT

Patrick Boston (Leeds University) and Mark Chapman (Edinburgh University)

securing your safety

Deng Shou Cheng 1,3, Wu Qing*,1, Chu Xiu Min 2 1 School of Logistics Engineering, Wuhan University of Technology,

GFC NIST MASS FLOW CONTROLLERS. Typical Stainless Steel GFC Mass Flow Controller. Design Features. General Description. Principles of Operation

This report presents an assessment of existing and future parking & traffic requirements, for the site based on the current development proposal.

Version IV: April a publication from

EQUIPEX NAOS WP5 : Deep oxygen floats in the North- Atlantic

1. Write down the ideal gas law and define all its variable and parameters. 2. Calculate the values and units of the ideal gas law constant R.

Travel Demand Management Implementation in Bandar Lampung

An Insect-based Detector of Bovine Manure Odors

GENETICS 101 GLOSSARY

operate regenerator top without boiling aq. amine solution.

-H- Note. Flow control valve. Key features

Climbing/Rappelling NATIONAL STANDARDS BOY SCOUTS OF AMERICA

PERFORMANCE TEAM EVALUATION IN 2008 BEIJING OLYMPIC GAMES

GFC NIST MASS FLOW CONTROLLERS. Typical Stainless Steel GFC Mass Flow Controller. Design Features. General Description. Principles of Operation

Analytical and numerical prediction of water sorptivity in cement-based materials

Characteristics of CNG Bubbles in Diesel Flow under the Influence of the Magnetic Field

A Data Envelopment Analysis Evaluation and Financial Resources Reallocation for Brazilian Olympic Sports

HYDRAULIC MOTORS MM APPLICATION CONTENTS GENERAL MOTORS

Capacity of Shared-Short Lanes at Unsignalised Intersections

Draft White-tailed Deer Management Policy for Ontario

THE RELATIONSHIP BETWEEN SELECTED KINEMATIC PARAMETERS AND LENGTH OF JUMPS OF THE SKI-FLYING COMPETITION

1 Bike MS: 2013 Proposal

Eagan YMCA Swim Lessons Schedule

"The twisting movement of any hoof should, for physiological reasons, not be hindered by Shoeing." (Lungwitz 1884)

Extensible Detection and Indexing of Highlight Events in Broadcasted Sports Video

Human-Robot Interaction: Group Behavior Level

draft final report NGSIM Arterial-Lane Selection Mode Federal Highway Administration Cambridge Systematics, Inc.

EMSBS/EMST. Drill For Machining Ultra-Deep Minute Holes FEATURES. For ultra-deep drilling of miniature holes. New chip stopper controls chip flow.

ELIGIBILITY / LEVELS / VENUES

ELIGIBILITY / LEVELS / VENUES

Sequential parimutuel games

DAMAGE ASSESSMENT OF FIBRE ROPES FOR OFFSHORE MOORING

Developing process safety indicators

Outline. Changing needs in Urban Traffic. Introduction The starting point Basic principles Analysis Facts Changes Context Solutions Key messages

Game Indicators Determining Sports Performance in the NBA

Andover YMCA Swim Lessons Schedule

» COLORADO s RIDE 2013

ELIGIBILITY / LEVELS / VENUES

Simulation Study of a Bus Signal Priority Strategy Based on GPS/AVL and Wireless Communications

ELIGIBILITY / LEVELS / VENUES

AN ESTIMATION OF INTER-ZONAL BUS DEMAND ON COLOMBO - KALUTARA TRANSPORT CORRIDOR

Lecture 13a: Chunks. Announcements. Announcements (III) Announcements (II) Project #3 Preview 4/18/18. Pipeline of NLP Tools

EFFICIENT ESTIMATION OF GAS LIQUID RATIOS FOR PLUNGER LIFT SYSTEMS IN PETROLEUM PRODUCTION OPERATIONS

Transit Development Plan December 2011

Daytime Road Marker Recognition Using Grayscale Histogram and Pixel Values

Intersleek Pro. Divers Manual. Our World is Water CONTENTS

Series 600 Accessories

Impacts of Typhoon on Wave Height at Bangkhuntien Shoreline

Transcription:

(IJACSA) Iteratioal Joural of Advaced Computer Sciece ad Applicatios, Vol. 6, No. 9, 5 Electrooculogram Sigals Aalysis for Process Cotrol Operator Based o Fuzzy c-meas Jiagwe Sog Raofe Wag* Guaghua Zhag Chaoxig Xiog Leya Zhag Cubag Su Abstract Biomedical sigals of huma ca reflect the body's task load, fatigue ad other psychological iformatio. Compared with other biomedical sigals, electrooculogram (EOG) has higher amplitude, less iterferece, ad is easy to detect. I this paper, the EOG sigals of operator s were aalyzed. Wavelet trasform was used to remove the highfrequecy artifacts. The fuzzy c-meas was adopted to detect the eye blik peak poits of EOG. After that, eye blik iterval () of operator was calculated. Four EOG features (the average of, variace of, stadard deviatio of ad variatio coefficiet of ) were extracted. Fially, the relatioship betwee EOG features ad operator s fatigue, effort, axiety ad task load were aalyzed. The experimetal results illustrate that EOG features had some relatio to the operator s fatigue, effort, axiety ad task load respectively. Keywords electrooculogram; fuzzy c-meas; operator fuctioal state; fatigue I. INTRODUCTION Fatigue is a atural physiological pheomeo, which is related to the operator s labor itesity, work eviromet, health, mood ad so o. It is a kid of self-adjustmet ad protectio fuctio of huma body. I recet years, with the rapid developmet of ecoomy ad the faster pace of life, may of the operators cotiue to work i the state of fatigue [,]. Cosequetly, may accidets occurred because of operator fatigue. Ad also the axiety ad effort level of the operators ca affect their work performace. The operator s performace degradatio may be the reaso of some serious disasters, particularly i safety-critical applicatios. The assessmet of operator metal workload would help predict the operator s periods of high operatioal risk. Early warig based o predictio of the operator s metal workload is oe of the effective methods to prevet such accidets. Therefore, the accurate assessmet of operator metal workload is the key to provide early warig successfully. Curret scietific research foud that, the huma body cotais may types of biological sigals which ca reflect i parts the operator s fatigue, axiety, metal workload, etc. [-5] After comparig the characteristics of electroecephalogram (EEG), electrocardiogram (ECG) ad electrooculogram (EOG), the authors aalyzed EOG sigals of operator s. Features of EOG sigals cotai eye blik rate, eye blik amplitude, closig time, eye blik variability, etc. [6] The researchers thought that EOG is associated with operator s metal workload because eye blik is associated with facial erve ad the awakeig is associated with cetral erve. Whe task load icreases, the visual demad is likely to icrease which may lead to eye blik rate decrease. Simultaeously, the eye activity may icrease. Whe the operator focuses o the tasks, his/her eyes blik rate may decrease ad the eye blik iterval may icrease [7,8]. I [9,], the authors aalyzed the driver s eye blik iterval. It is illustrated that whe the driver is fatigue, his/her eyelid ptosis arises, eye blik amplitude decreases. I other words eye blik iterval decreases while closig time icreases. This paper is to ivestigate the EOG sigals associated with the operator fuctioal state (OFS), like fatigue, axiety, effort ad task load. The experimetal eviromet of process cotrol was desiged, ad the EOG sigals of operators uder multi-level of task loads were sampled. Fuzzy c-meas (FCM) was adopted to detect the eye blik peak poits of EOG. Ad the EOG features to assess the huma fatigue, axiety, effort ad task load were foud through the study of EOG. Thereby, these idices will be used to predict the fuctioal state of the operator. Early warig system ca be costructed to prevet accidets which are caused by huma factors. 8 P a g e

(IJACSA) Iteratioal Joural of Advaced Computer Sciece ad Applicatios, Vol. 6, No. 9, 5 The rest of the paper is arraged as follows. Sectio described the data acquisitio experimet i details. I sectio, the Fuzzy c-meas algorithm is preseted. I sectio, the EOG sigal aalyzig process ad results are illustrated, ad the discussios are preseted. Sectio 5 makes the coclusios ad gives the future works. II. EXPERIMENT A. Data Acquisitio The experimets were carried out o Automatio-ehaced cabi air maagemet system (AUTOCAMS), developed by Hockey ad his colleagues [] ad modified by Lorez []. AUTOCAMS simulates a life support system of a spacecraft. It requires the subject to maage a semi-automatic system which is to maitai comfortable atmospheric coditio of the cabi. There are four sub-systems correspodig to four critical parameters of atmospheric coditio (oxyge cocetratio, carbo dioxide cocetratio, temperature, pressure). The subject was istructed to cotrol some of the sub-systems maually accordig to the experimetal desig. Subjects were traied o the maual process cotrol tasks for over 7 hours, followed by a further 5 hour experiece of cotrollig the system uder a rage of coditios. After that, the subjects are regarded as the adept at operatig AUTOCAMS. Three subjects (SA, SB, ad SC) participated i this experimet. The task of the subject was defied to cotrol some of the four subsystems so as to maitai their relevat variables withi target rages. There are 8 cotrol coditios (C~C8) i the experimet. The first (C) ad last (C8) coditios are five miute baselie testig coditios i which the subjects task was to moitor the operatio of AUTOCAMS. The amout of maually cotrolled subsystems i other cotrol coditios (C~C7) is,,,,,, respectively ad each coditio lasts for 5 miutes. For a experimetal sessio, the cyclic loadig method resulted i a loadig phase (C~C) ad a uloadig phase (C5~C7). Every subject was required to do the experimet for three times (e.g. SA, SA, SA) at the same time o differet days to avoid circadia effect. B. EOG Data The subjects vertical EOG sigals were recorded by Niho-Kode EEG physiological recorder from sites, which are placed above ad uder the subject s right eye ad refereced to liked earlobes. The EOG data were sampled at 5 Hz. III. METHODS A. Fuzzy c-meas Fuzzy c-meas is a data clusterig method i which each data poit belogs to a cluster, with a degree specified by a membership grade [-5]. FCM partitios a collectio of vectors X x, x,, x, where xk xk, xk,, x ks is the kth vector of X. FCM is based o miimizig a objective fuctio, which respects to fuzzy membership U, ad set of cluster cetroids, V. The FCM objective fuctio ca be described as, c m Jm ( U, V ) ik dik k i () c s. t. ; [,]; ik ik ik i k where, c is the umber of clusters. ik is the membership fuctio of vector x k to the ith cluster. d ik xk v i ca be either Euclidea distace or oe of its geeralizatios such as Mahalaobis distace. Euclidea distace is used i this study. The steps of the FCM method are explaied i brief: Firstly, the ceters of each cluster V(), fuzzy coefficiet m, umber of cluster c ad termiatio threshold are radomly selected. Secodly, the membership matrix with the followig equatio: ik ik jk j () t U is computed ( ) c m t d ( t) d ( t ). () Thirdly, the cth cluster cetroid is calculated with the followig equatio: m i ik k ik k k m v ( t ) ( t) x ( t ). () The process is stopped if ( t) ( t) V V. I this paper, the FCM algorithm was used to separate EOG sigals ito subsets, i.e. the higher amplitude subset ad the lower amplitude subset. The the eye blik peak of EOG ca be idetified. IV. RESULTS AND DISCUSSION The process of EOG aalysis is as follows: ) Deoisig. Wavelet Trasform was used to remove the high frequecy artifacts of EOG sigal. The parameter settigs of wavelet deoisig method are: wavelet base: DB, decompositio level:. Soft threshold deoisig method ad ubiased threshold estimatio method was adopted. ) Fuzzy c-meas. The EOG sigal was divided ito two clusters. The eye blik peak was foud i the higher amplitude cluster. Figure shows the eye blik peak poit of EOG sigals. It ca be see that the proposed method ca fid the eye blik peak correctly. The eye blik peaks for all the EOG sigals were detected by usig the same method. ) Feature extractio. The eye blik iterval was calculated. Figure shows the mea eye blik iterval of miute of SC. The teds to decrease (icrease) whe task load icrease (decrease). For the purpose of aalyzig the relatio of EOG to OFS, four features of, i.e. average of (A), variace of (V), stadard deviatio of (SD) ad variatio coefficiet of (VC), were extracted. The value of each feature i each coditio was 9 P a g e

average of (s) variatio coefficiet of (s) stadard deviatio of (s) variace of (s) (IJACSA) Iteratioal Joural of Advaced Computer Sciece ad Applicatios, Vol. 6, No. 9, 5 calculated to observe its chages i the whole experimet. Figure -6 show the chages of these four features i the whole experimet of SC. Figure 7 shows the umber of subsystems uder maual cotrol ad subjective measuremets i eight cotrol coditios of SC, respectively. - -..6.8...6.8 x Fig.. The results of fuzzy c-meas: lie EOG sigal;* eye blik peak poit 8 6 7 6 5 coditio Fig.. Variace of 8 7 6 5 coditio Fig. 5. Stadard deviatio of 6 8 time (mi) Fig.. per miute 5 coditio Fig.. Average of The average of, variace of ad stadard deviatio of decrease whe the operator s task load icreases from C to C. Ad these features show the same chages whe the operator s task load decrease from C5-C8. It is clear that the average of, variace of ad stadard deviatio of has a sigificat relatio to operator s task load..8.6...8 coditio Fig. 6. Variatio coefficiet of The variatio coefficiet of icreases from C-C5, which may be iflueced by the operator s icreasig effort ad axiety. I C ad C8, AUTOCAMS rus automatically. The operator was required to moitor the system. Thus there was t ay task load for the operator. I C-C7, the operator eeds to cotrol some of the subsystems. From all the four features, it could be see that there are sigificat differeces betwee o-load coditios (C ad C8) ad uder load coditios (C-C7). The differece of features betwee C ad C8 may be caused by the subject s fatigue. P a g e

subjective measures (%) umber of subsystems uder maual cotrol (IJACSA) Iteratioal Joural of Advaced Computer Sciece ad Applicatios, Vol. 6, No. 9, 5 Fig. 7. The umber of subsystems uder maual cotrol ad subjective measuremets These fids idicate that the chages of task load ad metal workload ca be reflected by features of EOG. The correlative aalysis is used to quatitatively aalyze the relatioship betwee EOG features ad operator task load, axiety, effort ad fatigue. Table shows the results of correlative aalysis betwee EOG features ad operator fatigue, axiety, effort ad task load of SC. Operator s fatigue ca be reflected by variatio coefficiet of. All these four features ca be idices of operator s axiety. Average of, variace of ad stadard deviatio of have sigificat relatios with operator s effort ad task. TABLE I. Operator Fuctioal State codito (a) 8 6 CORRELATION ANALYSIS BETWEEN EOG FEATURES AND OPERATOR FUNCTIONAL STATE Stadard Deviatio of Average of Variace of Fatigue -. -.6 -.56 -.68 Axiety -.7 -.77 -.8 -.76 Effort -.97 -.85 -.9 -.69 Task -.95 -.78 -.89 -.6 TABLE II. fatigue axiety effort codito (b) Variatio Coefficiet of SIGNIFICANT EOG FEATURES RELATED TO OPERATOR FUNCTIONAL STATE No. Fatigue Axiety Effort Task load SA VC VC VC SA V, SD VC V, SD SA A, SD A, SD A, SD SB V, SD VC SB A,V A SB V, SD A,V SC V,VC V,VC V,VC V, SD SC V,VC V, SD A, SD A, SD SC V,VC A,VC A,VC A,VC The same aalysis was doe to the EOG sigals for all the subjects. Table shows the top sigificat EOG features (P<.5) for each subject. Variace of is a good idex of operator fatigue. Stadard deviatio of ca reflect the operator s axiety, effort ad task load. It could be see that the saliet features for differet subjects are differet owig to idividual differece. For subject SA, stadard deviatio of is a saliet feature of operator axiety, effort ad task load. For subject SB, variace of ca sigificatly reflect the chages of operator fatigue. For subject SC, variace of ad variatio coefficiet of are sigificat features of fatigue ad axiety. Average of ca reflect the chages of effort ad task load. V. CONCLUSION This paper aalyzed the EOG sigal of operator. Fuzzy c- meas was used to detect the eye blik peak of EOG. Average of, variace of, stadard deviatio of ad variatio coefficiet of were extracted. The relatio of EOG features ad operators fatigue, axiety, effort ad task load was aalyzed. Saliet EOG features for each subject are successfully extracted. The experimet results show that it s ecessary to assess the operator fuctioal state usig persoal saliet features. The results of this paper could be used to idetify operator fuctioal sate i future works. Future works will focus o two areas. First, the authors will do more experimets ad further aalysis o EOG, i order to evaluate the idividual differeces. Secod, models that describe the relatio betwee EOG features ad operator fuctioal state will be costructed. ACKNOWLEDGMENTS This work is partially supported by the Shaghai Uiversity Youg Teachers' Traiig Scheme Fuds ZZGJD, the Natioal Nature Sciece Foudatio uder Grat No. 69 ad No. 68. REFERENCES [] G. Borghii, L. Astolfi, G. Vecchiato, D. Mattia, F. Babiloi, Measurig europhysiological sigals i aircraft pilots ad car drivers for the assessmet of metal workload, fatigue ad drowsiess, Neurosci. Biobehav. R. Kidligto, vol., pp. 58-75, July. [] R. Wag, J. Zhag, Y. Zhag, X. Wag, Assessmet of huma operator fuctioal state usig a ovel differetial evolutio optimizatio based adaptive fuzzy model, Biomed. Sigal Process. Cotrol, Kidligto, vol. 7, pp. 9-98, September. [] K. L. Saroj, C. Ashley, Driver fatigue : electroecephalography ad psychological assessmet, Psychophysiol. Malde MA, vol. 9, pp. -, May. [] I. Kather, S. C. Wriessegger, G. R. Muller-Putz, A. Kubler, S. Halder, Effects of metal workload ad fatigue o the P, alpha ad theta bad power durig operatio of a ERP (P) brai computer iterface, Biol. Psychol. Amsterdam, vol., pp. 8-9, October. [5] S. Rodriguez-Ruiz, P. M. Guerra, S. Moreo, M. Carme Feradez, J. Vila, Heart rate variability modulates eye-blik startle i wome with bulimic symptoms, J. Psychophysiol. Gottige, vol. 6, pp. -9, October. [6] H. A. Slagter, K. Georgopoulou, M. J. Frak, Spotaeous eye blik rate predicts learig from egative, but ot positive, outcomes, Neuropsychologia, Kidligto, vol. 7, pp. 6-, May 5. [7] M. J. Doughty, Further assessmet of geder- ad blik patterrelated differeces i the spotaeous eyeblik activity i primary gaze i P a g e

(IJACSA) Iteratioal Joural of Advaced Computer Sciece ad Applicatios, Vol. 6, No. 9, 5 youg adult humas, Optometry Visio Sci. Philadelphia, vol. 79, pp. 9-7, July. [8] P. L. Kaufma, A. ALM, Adler s physiology of the eye, th ed. St. Louis: Mosby,. [9] M. A. Recarte, L. M. Nues, Metal workload while drivig: effects o visual search, discrimiatio, ad decisio makig, J. Exp. Psychol. Appl. Washigto, vol. 9, o., pp. 9-7, Jue. [] M. Boyer, M. L. Cummigs, L. B. Spece, E. T. Solovey, Ivestigatig metal workload chages i a log duratio supervisory cotrol task, Iteract. Comput. Oxford, May 5. [] G. R. J. Hockey, D. G. Wastell, J. Sauer, Effects of sleep deprivatio ad user-iterface o complex performace: a multilevel aalysis of compesatory cotrol, Huma Factors. vol., o., pp. -5, Jue 998. [] B. Lorez, F. DiNocera, S. Rottger, R. Parasurama, Automated faultmaagemet i a simulated spaceflight micro-world, Aviat. Space Eviro. Med. Washigto, vol. 7, 886-897, September. [] J.C. Bezdek, Patter recogitio with fuzzy objective fuctio algorithms, New York: Pleum Press, 98. [] J. C. Bezdek, Patter Recogitio with Fuzzy Objective Fuctio Algorithms. Spriger Sciece & Busiess Media,. [5] N. R. Pal, K. Pal, J. M. Keller, J. C. Bezdek, A possibilistic fuzzy c- meas clusterig algorithm, IEEE Tras. Fuzzy Syst. vol., pp. 57-5, August 5. P a g e