A Fully-Integrated IC with 0.85-uW/Channel Consumption for Epileptic iEEG Detection

A Fully-Integrated IC with 0.85-uW/Channel Consumption for Epileptic iEEG Detection
of 5
All materials on our website are shared by users. If you have any questions about copyright issues, please report us to resolve them. We are always happy to assist you.
  114 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—II: EXPRESS BRIEFS, VOL. 62, NO. 2, FEBRUARY 2015 A Fully Integrated IC With 0.85- µ W/ChannelConsumption for Epileptic iEEG Detection Mahsa Shoaran,  Student Member, IEEE  , Claudio Pollo, Kaspar Schindler, and Alexandre Schmid,  Member, IEEE   Abstract —Feature extraction from a multichannel compressedneural signal is introduced in this brief. Compressive sensing (CS)is an efficient method for reducing the transmission data rate of sparse biological signals and lowering the power consumption of resource-constrained sensor nodes. However, recovering the src-inal signal from compressed measurements is typically achievedby relatively complex and optimization-based algorithms, whichis hardly suitable for real-time applications. The previously pro-posed multichannel CS scheme enables the area-efficient imple-mentation of CS. In this brief, a low-power feature extractionmethod based on line length is directly applied in the compresseddomain. This approach exploits the spatial sparsity of the signalsrecorded by adjacent electrodes of a sensor array and detectsthe seizure onset for every sixteen channels of the array. Theproposed circuit architecture is implemented in a UMC 0.18- µ mCMOS technology. Extensive performance analysis and design op-timization enable a low-power and compact implementation. Theproposed feature extractor reaches a perfect sensitivity of 100%for 420 h of clinical data containing 23 seizures from four patients,with an average false alarm rate of 0.34  h − 1 for artifact-freechannels, consuming 0.85 µ W of power/channel at a compressionrate of 16.  Index Terms —Compressive sensing (CS), feature extraction, in-tracranial electroencephalography(iEEG), seizure onsetdetection. I. I NTRODUCTION R ECORDING with high spatiotemporal resolution enablesneuroscientists to investigate the fine-scale electrophysi-ological activities within the brain and provides them with apreviously unexplored insight into the mechanisms underlyingneurological disorders. In epileptic patients, the submillimeterscale of high-frequency oscillations involved in seizure genera-tion [1] motivates the wide-bandwidth recording of intracranialelectroencephalography (iEEG) signals using high-spatial-resolution electrodes. As opposed to standard widely spaced( ≈  1 cm) electrodes with large surface areas, the spacing anddiameter of dense intracranial electrodes can be as small as0.5 mm and 0.3 mm, respectively [2]. In addition, high- Manuscript received October 3, 2014; accepted November 28, 2014. Dateof publication January 5, 2015; date of current version February 7, 2015. Thiswork was supported by the Swiss National Science Foundation under Grant200021-130166. This brief was recommended by Associate Editor R. Rieger.M. Shoaran and A. Schmid are with the Microelectronic Systems Laboratory,Swiss Federal Institute of Technology (EPFL), 1015 Lausanne, Switzerland(e-mail: mahsa.shoaran@epfl.ch).C. Pollo is with the Department of Neurosurgery, Bern University Hospital,3010 Bern, Switzerland.K. Schindler is with the Department of Neurology, Inselspital, Bern Univer-sity Hospital and University of Bern, 3010 Bern, Switzerland.Color versions of one or more of the figures in this paper are available onlineat http://ieeexplore.ieee.org.Digital Object Identifier 10.1109/TCSII.2014.2387652 resolution electrophysiology is beneficial for patients withpharmacoresistant epilepsy, as 25% of epileptic patients. De-veloping a wireless implantable device capable of deliveringelectrical stimulation or delivering a drug upon detecting animminent seizure can make a great impact on the quality of lifeof these patients [3]. It is expected that the inclusion of signalsrecorded from dense electrodes allows the system to improvethe accuracy of seizure detection, as compared to traditionaldetection methods based on low-precision electrode arrays [4].In order to preserve the benefits of high-resolution recordingand satisfy the power requirements of the implantable device,appropriate data reduction methods are employed. Discretewavelet transform [5], spike detection [6], and compressivesensing (CS) [7], [8] are the most popular approaches proposedin literature. Among these methods, CS is a computationallysimple approach, which is successfully translated into circuitlevel [7], [8], with relatively low power consumption and highcompression rate.In order to alleviate the area cost of implementing CS in adense array with stringent area constraints on each channel,a multichannel compression method has been recently pro-posed [8]. However, the speed and computational complexityof the recovery algorithm is a major drawback of CS-basedapproaches.In this brief, spatial-domain CS in combination with featureextractionisproposedasanefficientmethodtolimittheamountof data acquired by a seizure detector device. As opposed tothe existing feature extraction methods, which are separatelyapplied to each channel, the proposed method (see Figs. 1(a)and (b)) is applied to the compressed signal of a group of channels. The proposed method is compared to the raw data-domain feature extraction, particularly focusing on the seizuredetection performance.II. C URRENT  S CHEMES AND  P ROPOSED  A RCHITECTURE The previously explored on-chip solution to seizure detec-tion consists of extracting the spectral energy of each elec-troencephalography (EEG) recording channel within severalfrequency bins [9]. The most challenging part in integratinga feature extraction block on-chip relates to minimizing thelarge required area and power consumption of the filter banks.Although useful techniques are introduced [10] in order tolimit the area of the filter banks, the area limitation becomesextremely stringent in a high-density system recording fromseveral channels and particularly at a high sample rate, whichnecessitates a higher number of filters. In this work, we propose 1549-7747 © 2015 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission.See http://www.ieee.org/publications_standards/publications/rights/index.html for more information.  SHOARAN  et al. : FULLY INTEGRATED IC WITH 0.85- µ W/CHANNEL CONSUMPTION FOR EPILEPTIC IEEG DETECTION 115 Fig. 1. (a) System-level view of the proposed method; signals recorded byevery 16 electrodes are processed by a microelectronic chip. (b) Block diagramof the system; offline recovery of the compressed data provides full access tothe srcinal signals. a method that includes the information from all channels of an electrode array and simultaneously improves the area andpower efficiency by applying a low-complexity feature ex-traction method to a group of channels (see Fig. 1) in thecompressed domain (CD).  A. MCS Technique The proposed feature extraction method builds upon themultichannel CS (MCS) system initially presented in [8].The CS transform is generally expressed by  y  = Φ x , wherethe  d -dimensional signal  x  (e.g., EEG signal, which is shownsparse in the Gabor basis [8]) is sampled as an  m -dimensionalmeasurement vector  y  with a compression ratio (CR) of   d/m ,i.e.,  m  d . In the MCS model [8], the random projection y  = Φ x  is applied to the spatially sampled values from  N  channels of the array. Assuming a CR of   N  , the individualsamples of   N   channels are multiplied by a random value of 0 or 1 and summed together as a single measurement.Employing MCS in hardware leads to a significant datareduction of the system and an area-efficient implementation.This approach circumvents the need to load each channelwith multiple random sequences and avoids placing additionalmemory on the chip.  B. Signal Recovery Cost  As opposed to the simple compression model applied inCS, the signal recovery is usually achieved by nonlinear andrelatively expensive optimization-based or iterative algorithms[11], generally transferring the following signal analysis (e.g.,seizure onset detection) into a base station. In spite of manyefforts aiming at improving the speed and accuracy of therecovery algorithms, the state of the art still lacks the capacityof real-time operation for many applications. However, in manysignal processing problems, only the specific features of thesignal are of interest, and the exact recovery is not necessary[11]. In order to leverage the benefits of data compression,developing algorithms that do not require the full data recoveryand directly perform feature extraction and processing in theCD is potentially interesting. C. Proposed CD Feature Extraction The complexity of making a reasonable decision based onthe single-channel detections in a dense recording array aswell as the area and power overhead of implementing thein-channel feature extraction circuits necessitate an improvedseizure detection strategy.In general, many software-based seizure detection ap-proaches employ sophisticated algorithms to increase the de-tection accuracy. However, the limited power and area budgetof an implantable device impose strict constraints on the choiceof appropriate seizure detection features. Extracting frequency-based features involves computationally intensive algorithmssuch as filtering and fast Fourier transform evaluation, which ishardly suitable for an on-chip implementation. Considering thehardware cost and detection accuracy, the coastline feature hasbeen selected and extracted from the multichannel iEEG signalprior to, and following, the compression. Coastline achievesthe best seizure detection performance among more than65 different time- and frequency-domain features [12]. Thisfeature is a measure of the line length between successive sam-ples and provides an appropriate characteristic of epileptiformiEEG since it increases at low-amplitude fast activities as wellas high-amplitude slow activities [13]. The coastline feature of block   i  is computed as C  i  =   1256  256 | x [ n ] − x [ n − 1] |  (1)where  x [ n ]  is the value of signal  x  at time  t n  =  n × T  S / H , with T  S / H  being the sample-and-hold (S/H) time of the analog-to-digital converter (ADC). Normalization is done by shifting thedigital feature output by 8 bits toward the LSB. The choice of the division factor (i.e., the summation window length) dependsonthesignalamplitude,thedetectionaccuracy,andtherequiredlatency of detection.III. M ICROELECTRONIC  A RCHITECTURE The circuit-level implementation of the proposed architec-ture is shown in Fig. 2. The neural signals recorded by every  N  electrodes of the array ( N   is equal to 16 in this design) are pro-cessed through amplification and filtering stages located insidethe channels. The differential outputs of the channels connectto the MCS block and randomly accumulate at the single-ended output of this stage. As shown in [8], the compressedsensing system is very robust against circuit nonidealities suchas noise. Thus, the compressed data stream is digitized using a  116 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—II: EXPRESS BRIEFS, VOL. 62, NO. 2, FEBRUARY 2015 Fig. 2. Schematic of the proposed multichannel feature extractor system, including the individual channels, an MCS block with 16 differential inputs and a singleoutput, an 8-bit successive approximation ADC, and a digital coastline extractor block. moderate-resolution successive approximation ADC. The serialoutput of the ADC can be used for a full recovery of the srcinaliEEG signal at the receiver side.The employed CS technique enables 16 times data reductionand power saving of the transmitter, as compared to the conven-tional noncompressed design. The digital output is delivered tothe 8-bit feature extractor block in parallel form. To improve thepower efficiency of the system, the bias current of the front-endamplifiers ( G 1  in Fig. 2) is doubled upon seizure detection atthe output comparator of the feature extraction block. This tech-nique avoids unnecessary power dissipation while normal (i.e.,seizure-free) iEEG is being recorded. Upon seizure detection,the current of the low-noise amplifier (LNA) increases, whichsimultaneously reduces the input-referred noise of the channeland enhances the achievable SNR at the output.  A. Signal Conditioning Blocks The in-channel signal conditioning circuit performs the sig-nal amplification and filtering in a three-stage topology. Thisarchitecture efficiently optimizesthedieareausage ofachannelat the given gain and bandwidth requirements. Compared to[8], stringent power requirements are considered in this design,and the power supply is reduced by 33%. In addition, the S/Hstage is embedded within the MCS block, alleviating the needto place a separate S/H and buffer circuit in each channel.To avoid aggressive increase of the input capacitors of theLNA, a topology similar to [14] is adapted in a differentialimplementation. The high-value feedback resistors are imple-mented by back-to-back MOS devices, which operate in thesubthreshold region. The total bias current of the operationaltransconductance amplifier(OTA)implementedin G 1  is0.7 µ Aand 1.4  µ A, in the two modes of operation (normal and seizurerecording). The current consumption of the OTA used in  G 2  is0.3  µ A, due to the relaxed noise requirements at this stage. Afolded-cascode OTA is implemented in both  G 1  and  G 2 . LargePMOS input pair limits the flicker noise at  G 1 , while smallertransistors are employed at the input pair of   G 2 ,A source-follower-based low-pass filter consuming a totalbias current of 13.8 nA is used to limit the high cutoff fre-quency. The total midband gain of the channel is 45.3 dB.The input signal is amplified in the frequency range of 30 Hz–1.7 kHz. This frequency range covers both ripplesand fast ripples, as two types of signals that are indicatorsof epilepsy. These signals are discernible using high-spatial-resolution electrodes. The total input-referred noise of thechannel integrated from 1 Hz to 100 kHz is 5.74  µV  rms  and7.07  µV  rms , at the high- and low-biasing conditions,respectively.  B. MCS Block  The differential outputs of the channels connect to themultiple-input single-ended summing stage (see Fig. 2), as thecore of the MCS block. This stage is controlled by a samplingsignal  (Φ S  )  and  N   random sampling signals  (Φ R 1 ,..., Φ RN  ) .The random sampling signals are generated using an on-chipmultiple-output pseudorandom binary sequence generator withminimal power and area overhead [8]. During sampling, thedifferential voltages at the outputs of the channels are sampledacross the two sampling capacitors  ( C  S  ) . During summation( Φ S   is off), the charge stored on the sampling capacitors of thechannels with a random value equal to one are transferred tothe capacitor in the feedback path  ( C  f  ) . Thus, the outputs of the channels are randomly accumulated at the output of theMCS block. A two-stage OTA, which provides a rail-to-railoutput swing (see Fig. 2), is used within the summing stage of the MCS block. A variable voltage gain is achieved through the2-bit controlling switches in series with the feedback capacitors(shown as a variable feedback capacitor in Fig. 2). In order toset the appropriate value of the gain in a real  in vivo  experiment,a patient-specific training phase is required.  SHOARAN  et al. : FULLY INTEGRATED IC WITH 0.85- µ W/CHANNEL CONSUMPTION FOR EPILEPTIC IEEG DETECTION 117 C. ADC and Feature Extractor  An 8-bit binary-weighted capacitive array with attenuationcapacitor is employed, which enables a compact implementa-tion of the successive approximation ADC. All the requiredsignals are generated using a single external clock of 64 kHz.The S/H signal of the ADC is generated at a rate equal to 1/16of the external clock, which enables 8-bit conversion of thesummed value.An 8-bit feature extractor circuit computes the coastlineparameter at the output of the ADC. The block diagram of the coastline extractor is shown in Fig. 2. The multibit outputof the ADC is passed through register stages, which capturetwo successive compressed values. The absolute difference of the two inputs is generated at the output of the subtractorand accumulated over a window length of 256. Extra bitsare employed in the arithmetic hardware to prevent overflowduring the accumulation. An 8-bit synchronous binary countercontrols the accumulation window and resets the register after256 consecutive additions. Within the coastline generationchain, the absolute differential value of the successive samplesis shifted toward the least significant bit by 8 bits, enablingthe division operation in (1). At the rising edge of the counterMSB output (with a period equal to  256 /f  CLK ), the accumu-lated value is compared with a trained threshold input using adigital comparator stage. The detection signal associated withthe coastline parameter is raised upon exceeding the thresholdvalue in a window length of 64 ms (see Fig. 4), at an ADCsampling rate of 4 kS/s. A low-power and low-power-delay-product 16-transistor adder cell (see Fig. 2) is employed asthe main building block of the subtractor, accumulator, andcomparator stages. The S-generation branch is excluded forimplementing the multibit comparison. Considering the relaxedspeed requirements inherent to the application, a ripple-carryadder is implemented to enable the multibit operation. Thetotal power consumption of the ADC and feature extractor is0.54  µ W. The feature extractor occupies an area of   300  µ m × 85  µ m on the chip.IV. S YSTEM AND  C IRCUIT -L EVEL  V ALIDATION The proposed seizure detection method was evaluated usingiEEG data from four patients with medically refractory epi-lepsy. The patients underwent presurgical evaluation at theInselspital in Bern, Switzerland (see Table I). The iEEG signalswere acquired using a 128-channel system and sampled at32 kHz (Neuralynx, Inc.).The proposed MCS method is applied to the subblocks of neighboring 16-channel iEEG signals. The coastline feature isextracted from16 channels ineach subblock and alsointheCD.A moving average is applied to the extracted features. In total,103 subblocks containing 23 seizures from four patients havebeen processed. The total length of the analyzed iEEG data is420 h and includes preictal, ictal, and postictal activities.The average sensitivity and false alarm rate (FAR) versusCR are shown in Fig. 3. Sensitivity is computed as the ratioof correct seizure predictions to the total number of registeredseizures. FAR is measured by the number of false positive TABLE IP ATIENT  I NFORMATION AND  C HARACTERISTICS OF  E PILEPSY Fig. 3. (Left) Average sensitivity and FAR and (right) latency of the coastlinefeature versus CR. detections in 1 h of recording. The whisker plot presented inFig. 3 shows the latency to the seizure onset versus CR. Latencyis defined as the amount of time after (or before if negative)electrographic seizure onset taken by the selected feature tocross the threshold and trigger a detection. The point CR  = 1 corresponds to the channel-based detection. For this purpose,themedianofthesortedvectorformedbythe16detectiontimesassociated with the channels of each subblock is selected. Thecalculated sensitivity and FAR are averaged over subblocks andpatients.As shown in Fig. 3, the sensitivity of the traditional approachbased on multiple channels’ detections is worse than a CDdetector, which exclusively handles a single compressed datastream. The features extracted in the CD represent the groupedactivity of several channels. The compressed output is a linearcombination of samples recorded from adjacent channels andis a reasonable  intermediate  representative of the subblock’stotal activity. Using a higher number of measurements (e.g.,CR  = 8 ), the system is able to achieve a more accurate repre-sentation of the signal, which results in a lower false detectionrate. The proposed system is easily adaptable to operate atCR  = 8 , by doubling the sampling frequency at the randomsummation stage and ADC. Considering the latency graph, thechannel-based coastline feature is mostly concentrated aroundthe median (zero), since its calculation method is similar to theelectrographic onset definition [13] used in our analysis.In order to assess the performance of the proposed design,multichannel recorded signals of the slices of a rat somatosen-sory cortex are applied to the implemented circuit. This signalincludes epileptiform burst activity and extracellularly detectedspikes (see Fig. 4). The seizure detection performance of the  118 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—II: EXPRESS BRIEFS, VOL. 62, NO. 2, FEBRUARY 2015 Fig. 4. (a) and (b) Measured seizure detection performance of the proposedcircuit; a window of length 256 ms at the beginning of the epileptic burst isperiodically applied into the circuit. Using an FPGA platform, the preloadedserial iEEG data are transferred into two off-chip eight-channel D/A converters,andtheanalogoutputsareappliedtotherecordingchannels.(c)Diemicrographof the fabricated chip and (d) the measured compressed output and an on-chipgenerated random sequence.TABLE IIP ERFORMANCE  C OMPARISON  W ITH  P UBLISHED  W ORK fabricated chip is shown in Fig. 4. A window length of 256 msat the beginning of the epileptic burst is periodically appliedto the circuit. Every 64 ms of the iEEG signal is mappedinto a seizure or nonseizure event, at the rising edge of thecounter MSB signal. A seizure is detected after 128 ms, whichcorrectly represents the seizure onset time falling within thesecond subwindow of the signal.The performance summary and a comparison with existingworks are shown in Table II. The total current consumptionof the proposed system is 17  µ A, which is drawn from a0.8-V power supply. The proposed solution incurs lower com-putational complexity in the feature extraction stage, due tothe reduction in the number of samples processed in the CD.This computational advantage could be further enhanced bythe fact that even fewer samples are required for detectionpurposes in the CD, than for signal reconstruction [11]. Inaddition, the spatial CS enables proper offline data recovery of the srcinal physiological signal at the receiver. This feature ishighly desired by physicians who prefer to have access to theentire iEEG data for a thorough clinical evaluation.V. C ONCLUSION In the future, dense neural recording interfaces will integratea high number of acquisition channels at high sampling rates,exacerbating the need to perform some type of data reductionat the sensors. This work represents the first fully integratedcircuit that addresses the multichannel CD feature extraction.A performance comparable to the raw data-domain detectoris retained up to CRs as high as 16, consuming 0.85  µ W of power within a small die area. The proposed system enables anefficient real-time seizure onset detection, which is integratedinside an implantable iEEG recording system for the treatmentof epilepsy.R EFERENCES[1] M. Stead  et al. , “Microseizures and the spatiotemporal scales of humanpartial epilepsy,”  Brain , vol. 133, no. 9, pp. 2789–2797, Sep. 2010.[2] J. Viventi  et al. , “Flexible, foldable, actively multiplexed, high-densityelectrode array for mapping brain activity in vivo,”  Nat. Neurosci. ,vol. 14, pp. 1599–1605, 2011.[3] W. M. Chen  et al. , “A fully integrated 8-channel closed-loop neural-prosthetic SoC for real-time epileptic seizure control,”  Proc. IEEE Int.Solid-State Circuits Conf. Dig. Tech. Papers , pp. 286–287, Feb. 2013.[4] S. S. Talathi  et al. , “Non-parametric early seizure detection in an animalmodel of temporal lobe epilepsy,”  J. Neural Eng. , vol. 5, no. 1, pp. 85–98,Mar. 2008.[5] A. M. Kamboh, M. Raetz, K. G. Oweiss, and A. J. Mason, “Area-powerefficient VLSI implementation of multichannel DWT for data compres-sioninimplantableneuroprosthetics,”  IEEETrans.Biomed.CircuitsSyst. ,vol. 1, no. 2, pp. 128–135, Jun. 2007.[6] A. Rodriguez-Perez, J. Ruiz-Amaya, M. Delgado-Restituto, andA. Rodriguez-Vazquez, “A low-power programmable neural spike de-tection channel with embedded calibration and data compression,”  IEEE Trans. Biomed. Circuits Syst. , vol. 6, no. 2, pp. 87–100, Apr. 2012.[7] F. Chen, A. P. Chandrakasan, and V. Stojanovic, “Design and analysis of a hardware-efficient compressed sensing architecture for data compres-sion in wireless sensors,”  IEEE J. Solid-State Circuits , vol. 47, no. 3,pp. 744–756, Mar. 2012.[8] M. Shoaran, M. Hosseini Kamal, C. Pollo, P. Vandergheynst, andA. Schmid, “Compact low-power cortical recording architecture for com-pressive multichannel data acquisition,”  IEEE Trans. Biomed. CircuitsSyst. , vol. 8, no. 6, pp. 857–870, Dec. 2014.[9] N. Verma  et al. , “A micro-power EEG acquisition SoC with integratedfeature extraction processor for a chronic seizure detection system,”  IEEE  J. Solid-State Circuits , vol. 45, no. 4, pp. 804–816, Apr. 2010.[10] J. Yoo  et al. , “An 8-channel scalable EEG acquisition SoC with patient-specific seizure classification and recording processor,”  IEEE J. Solid-State Circuits , vol. 48, no. 1, pp. 214–228, Jan. 2013.[11] M. A. Davenport, P. T. Boufounos, M. B. Wakin, and R. G. Baraniuk,“Signal processing with compressive measurements,”  IEEE J. Sel. TopicsSignal Process. , vol. 4, no. 2, pp. 445–460, Apr. 2010.[12] L. Logesparan, A. J. Casson, and E. Rodriguez-Villegas, “Optimal fea-tures for online seizure detection,”  Med. Biol. Eng. Comput. , vol. 50,no. 7, pp. 659–669, Jul. 2012.[13] K. Schindler, H. Leung, C. E. Elger, and K. Lehnertz, “Assessing seizuredynamics by analysing the correlation structure of multichannel intracra-nial EEG,”  Brain , vol. 130, pp. 65–77, Jan. 2007.[14] K. A. Ng and Y. P. Xu, “A compact, low input capacitance neuralrecording amplifier,”  IEEE Trans. Biomed. Circuits Syst. , vol. 7, no. 5,pp. 610–620, Oct. 2013.
Similar documents
View more...
We Need Your Support
Thank you for visiting our website and your interest in our free products and services. We are nonprofit website to share and download documents. To the running of this website, we need your help to support us.

Thanks to everyone for your continued support.

No, Thanks

We need your sign to support Project to invent "SMART AND CONTROLLABLE REFLECTIVE BALLOONS" to cover the Sun and Save Our Earth.

More details...

Sign Now!

We are very appreciated for your Prompt Action!