Biomedical Signals. Signals and Images in Medicine Lecture 1 Dr Nabeel Anwar

Size: px

Start display at page:

Download "Biomedical Signals. Signals and Images in Medicine Lecture 1 Dr Nabeel Anwar"

Regina Dalton
5 years ago
Views:

1 Biomedical Signals Signals and Images in Medicine Lecture 1 Dr Nabeel Anwar

2 Books 1. The Scientist and Engineer's Guide to Digital Signal Processing By Steven W. Smith, Ph.D. This book is written in simple English with minimum mathematics. The examples are very easy to understand. Consult this book for understanding a topic. 2. Biomedical Signal Analysis: A Case-Study Approach by Rangaraj M. Rangayyan. This book has case study approach, most of the practical examples would be taken from this book.

3 What we know so far.. Living organisms are made up of many component systems: the human body includes several systems. Each system is made up of several subsystems that carry on many physiological processes. Most physiological processes are accompanied by signals of several types that reflect their nature and activities The first link in the biomedical measurement chain is typically a transducer or sensor, which measures signals (such as a heart valve sound, blood pressure, or X-ray absorption) Some such electrical signals occur spontaneously (e.g., the electroencephalogram, EEG, ECG, EMG); others can be observed upon stimulation (e.g., evoked potentials (EPs)).

4 Examples of Biomed Signals 1. The action potential: Action potential (AP): electrical signal from a single cell when stimulated by an electrical current (neural or external). Recording an action potential requires the isolation of a single cell, and microelectrodes with tips of the order of a few micrometers to stimulate the cell and record the response.

5 Examples of Biomed Signals The electro-neuro-gram (ENG) The ENG is an electrical signal observed as a stimulus and the associated nerve action potential propagate over the length of a nerve. Conduction velocity in a peripheral nerve measured by stimulating a motor nerve and measuring the related activity at two points at known distances along its course.

6 Examples of Biomed Signals EEG/ECoG and Evoked Potentials (EPs) The electroencephalogram (EEG) represents overall brain activity recorded from pairs of electrodes on the scalp. In special cases, brain activity may also be directly measured via electrodes on the cortical surface (the electrocorticogram, ECoG,) or via depth electrodes implanted in the brain. Another common class of neurophysiological signals used for clinical tests are auditory-, visual-, and somatosensory-evoked potentials (AEP, VEP, and SSEP, respectively). These signals represent the brain s response to a standard stimulus such as a tone burst, click, light flash, change of a visual pattern, or an electrical pulse delivered to a nerve.

7 Examples of Biomed Signals ECG (EKG) The activity of the heart is associated with a highly synchronized muscle contraction preceded by a wave of electrical activity. The ECG is usually characterized by several peaks, denoted alphabetically P-QRS-T. EMG The activity of the muscles. (for details use lecture slides from bioinstruments).

8 Examples of Biomed Signals The electro-gastro-gram (EGG) Electrical activity of the stomach: rhythmic waves of depolarization and repolarization of smooth muscle cells. Gastric myoelectric activity is composed of mainly two complementary rhythms; slow wave activity, responsible for muscle contraction timing, and electrical response activity, responsible for triggering peristaltic contractions.

9 Examples of Biomed Signals The phono-cardio-gram (PCG) PCG: vibration or sound signal related to the contractile activity of the cardiohemic system (heart and blood). Recording the PCG requires a transducer to convert the vibration or sound signal into an electronic signal: microphones, pressure transducers, or accelerometers.

10 Examples of Biomed Signals The carotid pulse (CP) Carotid pulse: pressure signal recorded over the carotid artery as it passes near the surface of the body at the neck. Pulse signal indicating the variations in arterial blood pressure and volume with each heart beat.

11 Examples of Biomed Signals The speech signal Speech produced by transmitting puffs of air from the lungs through the vocal tract as well as the nasal tract. Shape of vocal tract varied to produce different types of sound units or phonemes which form speech. Speech signal of the word safety uttered by a male speaker. Approximate time intervals of the various phonemes in the word are /S/: s; /E/: s; /F/: s; /T/: transient at 1.1 s; /I/: s.

12 Objectives of Biomedical Signal Analysis 1. Information gathering measurement of phenomena to interpret a system. 2. Diagnosis detection of malfunction, pathology, or abnormality. 3. Monitoring obtaining continuous or periodic information about a system. 4. Therapy and control modification of the behavior of a system based upon the outcome of the activities listed above to ensure a specific result. 5. Evaluation objective analysis to determine the ability to meet functional requirements, obtain proof of performance, perform quality control, or quantify the effect of treatment.

13 Suggested readings for Lecture 1. DSP guide, CH3 (35-49): main lecture DSP guide, CH2 (all) : for understanding few concepts of CH3 Lecture slides Homework: Chapter 2 of Rangaraj s book. Solve problem 2.2.1, 2.2.2, 2.3 using matlab You are provided with the.dat file readable in the matlab, plus the code to read it.

14 You need to have a basic background about mean, standard deviation and variance of a data. (read chapter 2 of DSPguide).

16 The Measurement chain

17 Analog and Digital Converting a continuously changing waveform (analog) into a series of discrete levels (digital)

18 Analog and Digital

19 Analog Signals Analog signals directly measurable quantities in terms of some other quantity Examples: Thermometer mercury height rises as temperature rises Car Speedometer Needle moves farther right as you accelerate Stereo Volume increases as you turn the knob.

20 Digital Signals Digital Signals have only two states. For digital computers, we refer to binary states, 0 and 1. 1 can be on, 0 can be off. Examples: Light switch can be either on or off Door to a room is either open or closed

21 Parts of an Analog Signal Period (T) Amplitude (peak-to-peak) Amplitude (peak) Frequency: 1 F Hz T

22 Amplitude Parts of a Digital Signal Amplitude: For digital signals, this will ALWAYS be 5 volts. Period: The time it takes for a periodic signal to repeat. (seconds) Frequency: A measure of the number of occurrences of the signal per second. (Hertz, Hz) Time High (t H ): The time, when the signal is at 5 v. Time Low (t L ): The time, when the signal is at 0 v. Duty Cycle: The ratio of t H to the total period (T). Rising Edge: A 0-to-1 transition of the signal. Falling Edge: A 1-to-0 transition of the signal. Time High (t H ) Period (T) Frequency: F 1 T Hz Time Low (t L ) t DutyCycle T Falling Edge Rising Edge H 100%

23 Periodic Signals They repeat themselves. A periodic signal can be decomposed into summation of harmonic sinusoids Summation of harmonic sinusoids results in a periodic signal We can use harmonic sinusoidal components to approximate any periodic signal.

24 Non-periodic Signals In order to synthesize a periodic signal, frequencies of the sinusoids must be harmonic, that is, integer multiples of a fundamental frequency; What happens if the sinusoids are not harmonic? A slight frequency change of a component will result in non-periodic signals

25 Non-periodic Signals Example: x h 2 2 ( t) 2cos(20 t ) cos(20 (3) t) cos(20 (5) t) 3 5 Frequencies are 10,30 and 50 respectively x ( t) 2cos(20 t ) 2 cos(20 3 8t) 2 cos( t 27 ) Frequencies are 10,28.8 and respectively

26 Periodic Signals

27 Non-periodic Signals

28 Non-periodic Signals Spectra are very similar!!

29 Analogue to Digital Microphones - take your voice varying pressure waves in the air and convert them into varying electrical signals. Strain Gages - determines the amount of strain (change in dimensions) when a stress is applied. Thermocouple temperature measuring device converts thermal energy to electric energy. What If we want to store/process the data with computers? We need to convert Analogue signal to Digital Signal.

30 A to D conversion

31 A to D conversion Analog-to-digital conversion (ADC) makes a continuous signal discrete, both in amplitude and time. We need to digitize the signal both in Amplitude Time

32 Analog Digital Conversion is a 2-Step Process: 1. Quantizing - breaking down analog value is a set of finite states 2. Encoding - assigning a digital word or number to each state and matching it to the input signal.

33 Step 1: Quantizing The discretization of the signal in the amplitude dimension is determined by the converter s input voltage range and the analog amplification of the signal input to it Output States Discrete Voltage Ranges (V) Example: You have 0-10V signals. Separate them into a set of discrete states with 1.25V increments. (How did we get 1.25V? See next slide )

34 Quantizing The number of possible states that the converter can output is: N=2 n where n is the number of bits in the AD converter Example: For a 3 bit A/D converter, N=2 3 =8. Analog quantization size: Q=(Vmax-Vmin)/N = (10V 0V)/8 = 1.25V

35 Quantizing Quantization noise At conversion, the amplitude of the analog signal is approximated by the discrete levels of the ADC. Depending on the type of converter, this approximation may behave numerically as a truncation or as a round-off of the continuous-valued signal to an integer. In both cases, one can consider the quantization as a source of noise in the measurement system, noise which is directly related to the resolution at the ADC.

36 Quantization Quantization errors can be reduced by increasing the number of bits Common for A-D converters to have 16 bit or better resolution However the accuracy of the reference voltage must be of the same precision

37 Encoding Here we assign the digital value (binary number) to each state for the computer to read. Output States Output Binary Equivalent

38 Quantization and encoding Output States Output Binary Equivalent Discrete Voltage Ranges (V)

39 Accuracy of A/D Conversion There are two ways to best improve accuracy of A/D conversion: increasing the resolution which improves the accuracy in measuring the amplitude of the analog signal. increasing the sampling rate which increases the maximum frequency that can be measured.

41 Resolution Resolution (number of discrete values the converter can produce) = Analog Quantization size (Q) (Q) = Vrange / 2^n, where Vrange is the range of analog voltages which can be represented In our previous example: Q = 1.25V, this is a high resolution. A lower resolution would be if we used a 2-bit converter, then the resolution would be 10/2^2 = 2.50V.

42 Digitizing noise Any one sample in the digitized signal can have a maximum error of ±½ LSB (Least Significant Bit). With standard deviation of 1/ sqrt(12) LSB (0.29 LSB). Example: Consider a signal of amplitude of 1.0 volt, and a random noise of 1.0 millivolt rms. Digitizing this signal to 8 bits results in 1.0 volt becoming digital number 255, and 1.0 millivolt becoming LSB. (The LSB will indicate, 1/255 = millivolt change in the signal. And 1.0 millivolt of noise is 1/ = LSB) The total noise of the digitized signal is Noise present inside the signal + quantization error = The range of LSB will determine the efficiency of the Digitizer. For a 1 volt signal 8 bit digitizer produce noise = 0.386*1/256 millivolt = 1.5 millivolt, it is 50% increase in the total noise.

43 Sampling The continuous signal is also discretized (sampled) in time. To obtain a reliable sampled representation of a continuous signals, the sample interval (Ts) or sample frequency (Fs = 1/Ts) must relate to the type of signal that is being recorded. Sampling at very fast rate and very slow rate, both cause problems.

44 Sampling: Ideally what we want is But sampling rate can damage our required output.

46 e.g. undersampling Sampling a 20-Hz sine wave at different rates Fs = 1/Ts.

47 This is called Aliasing It occurs when the input signal is changing much faster than the sample rate. In previous example, a 20 Hz sine wave being sampled at 24Hz would be reconstructed as a ~4 Hz (the aliased signal) sine wave. Alias frequency is the absolute value of the difference between the frequency of the input signal and the closest integer multiple of the sampling rate. Nyquist Rule: For a band limited signal, use a sampling frequency at least twice as high as the maximum frequency in the signal to avoid aliasing.

48 Example: Assume fs, the sampling frequency, is 100 Hz and that the input signal contains the following frequencies: 25 Hz, 70 Hz, 160 Hz, and 510 Hz. These frequencies are shown in the following figure. As shown in the following figure, frequencies below the Nyquist frequency (fs/2 = 50 Hz) are sampled correctly. Frequencies above the Nyquist frequency appear as aliases. For example, F1 (25 Hz) appears at the correct frequency, but F2 (70 Hz), F3 (160 Hz), and F4 (510 Hz) have aliases at 30 Hz, 40 Hz, and 10 Hz, respectively. Alias F2 = = 30 Hz Alias F3 = (2) = 40 Hz Alias F4 = (5) = 10 Hz

49 Aliasing Consider a band limited signal with maximum frequency content is 20 Hz. Nyquist Rule: The minimum sampling rate (in this example twice the maximum frequency, i.e. 40 Hz) is calledthe Nyquist limit. Thus, the sampling rate determines the highest frequency that can be represented by the sampled signal. This value (half the sample rate) is often indicated as the Nyquist frequency of the sampled signal.

50 Aliasing The Nyquist limit is a real bare minimum to capture the 20-Hz frequency component, and you can see in the figure (previous slides) that the wave morphology is already distorted at sample rates close to, but above, the Nyquist sampling frequency. Clearly the signal is seriously misrepresented below the Nyquist limit. This particular type of signal distortion is called aliasing.

51 In summary Aliasing occurs if Sampling below the Nyquist rate, Improper reconstruction, or Both

52 Antialiasing What can we do about aliasing? Sample more often Problem is huge dataset at the end Make the signal less wiggly Get rid of some high frequencies Will lose information But it s better than aliasing

53 Preventing aliasing Introduce lowpass filters: remove high frequencies leaving only safe, low frequencies choose lowest frequency in reconstruction (disambiguate)

54 Preventing aliasing To remove the effect of aliasing in digitized signals, the analog measurement chain must remove/attenuate all frequencies above the Nyquist frequency by using a filter (anti-aliasing filter).

55 Real life aliasing problems consider a continuous signal that is a pure sinusoid: x1 t = cos (0.4ωt) where ω = 2 pi f Now consider another signal x2 t = cos (2.4ωt) X2 has 6 times X1 s frequency. x2 t = cos (0.4ωt+2ωt), since Which gives: x2 t = cos (0.4ωt) This means that given a pair of sampled versions of the signals, one with the low frequency sinusoid and one with higher, you will have no way of distinguishing these signals from one another.

56 Real life aliasing problems Consider Two possible input sine waves are shown: one has a frequency of 110Hz, the other has a frequency of 1110Hz. Both are sampled at 1000Hz. With 1000Hz sampling rate we can exactly reproduce 110Hz signal. But what will happen with 1110Hz signal? By definition in previous slides the 1110Hz will be aliased into 110Hz signal as well. So we have a sampled signal consists of both 110Hz original signal and aliased signal of 1110Hz. This aliased signal will not only destroy the time domain signal but also its frequency components. You will never know which part is real and which part is aliased.

57 Setting the sampling rate So selecting a sampling rate is very critical. Intuitively, we are not capturing much information in the signal by only sampling at some random frequency Goal Capture a sufficient amount of information from the signal But: Don t sample too fast-- will have too much data

58 Setting the sampling rate Is Nyquist-Shannon rule is good enough? Lets try to find out with these three examples 1. I am going to sample at 8kHz, so I need to use a filter with a 4kHz cutoff. 2. I need to monitor the 60Hz power line, so I need to sample at 120Hz. 3. What is the spectrum of an EKG signal? I want to figure out the Nyquist rate.

59 1. I am going to sample at 8kHz, so I need to use a filter with a 4kHz cutoff. Let s take the human voice recording example, Human can hear upto 20k Hz sounds. However, the normal range of the usable frequency band of human speech is approximately 300 Hz to 4000 Hz. So the first step should be to use anti aliasing filter which can remove frequency components above 4000 Hz. Figure shows a first order filter applied to the signal. The solid curve shows the spectrum of the filtered speech, while the dashed curve shows the aliased speech. There is a big component of aliased speech, Which is bad.

60 1. I am going to sample at 8kHz, so I need to use a filter with a 4kHz cutoff. Can we increase the filter order to improve the results? The figure shows a sixth order filter, the situation improved but still aliasing is a problem. If you go to filter order of 3500 hundred you can eliminate the aliasing considerably. Question is whether that high order analogue filters are possible and feasible?

61 1. I am going to sample at 8kHz, so I need to use a filter with a 4kHz cutoff. Intelligible speech transmission is that you need to transmit frequencies up to 3kHz, and it gives a good quality sound. Now lets redraw the figure with sampling frequency of 8kHz with frequency of interest is under 3kHz. The figure shows a sixth order filter. The results are much better with almost negligible aliasing. If you want to stuck with 4kHz then you should consider increasing sampling frequency.

62 3. What is the spectrum of an ECG signal? I want to figure out the Nyquist rate. If the signal is unknown then one can find out the Nyquist rate by taking the spectra of the signal. In the case of ECG the spectra is mostly limited to 100Hz. So it seems OK if we apply a sharp low pass filter at 100Hz. Now what happened after reconstruction (blue) of the EKG, there is a delayed QRS peak with higher amplitude and some nonlinearity (ringing effect) at the dropping component.

63 3. What is the spectrum of an ECG signal? I want to figure out the Nyquist rate. If we change the type of filter, the ringing effect might be reduced, but still the amplitude difference remains. What, then, do we do? One obvious solution is that if you are working with an application that cannot tolerate these effects is to sample your signal faster. This will allow you to use a more open anti-aliasing filter.

64 3. What is the spectrum of an EKG signal? I want to figure out the Nyquist rate. Consider a rectangular pulse passes (blue) thorough a Anti-aliasing filter. The resultant (black) has two evident effects: a time delayed signal and ringing effect. The delay is constant then we can manage it by fitting it in overall delay of the system. But ringing effect is still present.

65 Finally you should be aware of the Nyquist rate when you re designing systems. First of all you must know your system and the limitations, e.g. decreasing sampling rate in the speech transfer over the phone might not be critical as compared to decreasing setting sampling rate for ECG. To really determine an appropriate sampling rate for a system, or to determine the necessary antialias and reconstruction filters for a system, you have to understand aliasing and filtering.

PHYS225 Lecture 22. Electronic Circuits

PHYS225 Lecture 22. Electronic Circuits PHYS225 Lecture 22 Electronic Circuits Last lecture Digital to Analog Conversion DAC Converts digital signal to an analog signal Computer control of everything! Various types/techniques for conversion