Data was collected from a set of twenty participants with COPD over a period of approximately 8 weeks in a purely observational manner (i.e. the recorded data was not used to guide clinical treatment in any way, and data was not analyzed until after the end of the observation period). This study was approved by the St Vincent's University Medical Ethics Committee. Patients provided written informed consent for study participation. Subjects were enrolled at the point of discharge from St Vincent's University Hospital, Dublin following an in-hospital stay for an exacerbation of their underlying COPD condition. The inclusion criteria required that patients had experienced an ECOPD-related hospitalization within the previous 12 months prior to enrolment, were Gold Stage 2 to 4 (GOLD 2014), were capable of giving informed consent and, for logistic reasons, lived within a 50 mile radius of the hospital. Note that the subset of patients selected (i.e. patients with a COPD related hospitalization) may not be representative of the general class of Gold Stage 2 to 4 patients. The demographics of the patients who took part in the study are listed in table1.
Table 1.Demographic details of subjects.
Patients studied | N = 20 |
Male Sex, n (%) | 11 (55%) |
Age, years | 68 ± 10 |
Body mass index (m kg−2) | 27.2 ± 5.9 |
Spirometry | |
FEV1 (litres) | 1.27 ± 0.62 |
FEV1,predicted (%) | 51 ± 20 |
FVC (litres) | 2.79 ± 1.09 |
FEV1/FVC (%) | 48 ± 16 |
GOLD | |
I | 1 |
II | 8 |
III | 10 |
IV | 1 |
Comorbidities | |
Left ventricular failiure, n (%) | 7 (35%) |
Right ventricular failure, n (%) | 6 (30%) |
Prior diagnosis of SDB, n(%) | 5 (25%) |
Following recruitment, patients were issued with a SleepMinder device, which is a biomotion sensor designed to collect respiration rate information during a person's sleep (the SleepMinder is more fully described in section2.3 below). This device did not require any patient intervention, as it continuously logged respiration rate information on to an SD card on a nightly basis. The SleepMinder was set to record on an automated schedule from 3 to 2 PM the following day to guarantee that the sleeping period was recorded. Patients were followed for an approximate 2 month period using the SleepMinder recording in their home environment, and returned the device to the study coordinators at the end of the 2 month period. An automated algorithm was developed to only analyze the portion of the signal during the night time period when the person was present in bed. As part of the protocol, events associated with an exacerbation of COPD (ECOPD) were recorded by the study co-ordinator: ECOPD events were defined as either hospitalizations or ECOPD-related non-scheduled visits to the emergency room(ER) or primary care doctor. If a person was hospitalized, there is a loss of SleepMinder data for that subject during their period of admission; similarly, there will be periods of missing data if the subject is out of their home for other reasons.
The non-contact radio frequency biomotion sensor (SleepMinder, ResMed, Dublin) used in this study has been described elsewhere (de Chazal et al 2011, Zaffaroni et al 2013, Pallin et al 2014, O'Hare et al 2014), but for convenience we provide a brief description of the technology here. The SleepMinder system has been developed specifically to measure movement in sleeping human subjects. This sensor operates by transmitting two short pulses of radio-frequency energy at 5.8 GHz (each pulse is approximately 5 ns long). The first pulse acts as the main transmit pulse, and the second pulse is the mixer pulse. The first pulse reflects off nearby objects to create an echo pulse that is received back in the sensor. By multiplying ('mixing') the echo pulse with the mixer pulse inside the receiver, a signal proportional to any phase shift of the echo pulse is generated. Moving objects (such as a person breathing) generate a variable phase shift that can be detected by the electronics in the circuit. In addition, it employs a technique to overcome a well-known limitation in radio frequency sensing called the range correlation effect (Droitcour et al 2004), which refers to the fact that if a moving object is undergoing small movements centered about a location with a fixed phase relationship to the sensor which is in a sensitivity null, then the movements will be detected poorly if only a single channel is used. The technique, referred to as quadrature detection, effectively means there are two modes of transmission, which leads to two estimates of the movement, called I and Q channels. Figure1 shows the illustrative size and positioning of the non-contact sensor relative to the monitored subject. Figure2 shows an example of a typical signal acquired by this sensor over 60 s, in which the main movement is chest movement associated with respiratory effort, but with a larger body movement (associated with a change of bodily position) in the middle of the record, followed by a return to stable breathing.
As seen in figure2, movement signals received in this device contain a mixture of all movements within the field of the sensor. For example, if a person is moving both arms and head, the received signal will reflect components of all three movements. Fortunately, during sleep, most movement is related to respiratory effort, which is the key to its utility in respiration rate monitoring.
The sensor is also designed to be directional (only measures movement in front of the sensor) and range limited (it only responds to objects which are within 1.5 m of the sensor). The system operates at low power (total average emitted power is less than 1 mW), to comply with safety and regulatory guidelines for radio-frequency devices.
The sensor system outputs two analog baseband voltage signals which represent the I and Q signals of the sensor. Once the baseband signals are available, signal processing techniques are then required for robust respiration rate (RR) estimation. In the ideal case of a high SNR signal, with infrequent movement artefacts, then simple non-parametric spectral estimation techniques could be used (such as a periodogram (PG) estimator technique). In practice, there are some limitations to such a technique (a) the signal may have periods of poor SNR and high movement artifact, (b) a reliable estimate may require a long time window, which is longer than the time period over which RR is useful to measure—therefore a simple periodogram based approach may be incapable of capturing the fine dynamics of the respiratory rate. In practice, we would like to obtain reliable estimates of RR over a time period of 30 s or less, since it is quite reasonable that the RR will vary over such epochs. Given these limitations of a periodogram approach, we have investigated other alternate approaches to measure the RR.
The technique used in this paper for extracting respiration rate from the raw respiratory movement recorded by the SleepMinder is an adaptive notch filter (ANF) (Nehorai 1985, Ballal et al 2012). We adapt the modified approach suggested by Tan who chose an Infinite Impulse Response (IIR) filter architecture for use in the adaptive filter (Tan and Jiang 2009). IIR filters are preferred to FIR (finite impulse response) filters when a target performance is required at low computational cost. Contrary to the general case for multi-pole IIR adaptive filters, no multiple local minima occur with an adaptive (second order) IIR notch filter, only a single global minimum (Chicharo and Ng 1990, Tan and Jiang 2009). Also, no stability issues have been noticed during the experimental work. Frequency estimation and tracking using an adaptive second-order IIR notch filter is described in this section. Here we follow the approach of Tan and Jiang (2009) for the single-harmonic case. The second order transfer function of an IIR notch filter is given by
where is the filter notch frequency and r is the parameter that controls the bandwidth of the notch. From the transfer function, the filter output for an input x[n], where n denotes discrete time, can be written as
For the frequency estimation/tracking case, θ is not known. If x[n] is a narrowband signal (i.e. a breathing signal with a single strong breathing frequency so that signal is approximately sinusoidal) with a center frequency , the filter output will be minimum when . Therefore, an estimate of can be obtained from an N-sample segment of x[n] using
In the above equation, could be obtained by evaluating the expression for all sample values of This approach may yield good frequency estimates but it may not be computationally efficient for real-time frequency tracking, unless the density of is reduced. Reducing the density of leads to a poorer resolution. Therefore in practice, a more effective approach is to use an adaptive algorithm to adjust the parameter θ and provide estimates of the frequency without the resolution limitation. The well-known least mean squares (LMS) algorithm can be used to carry out this task. The LMS algorithm for updating the filter parameter θ, to track the input signal frequency (which coincides with minimizing ) can be derived as:
where μ is the adaptive step-size. For a good convergence rate, a suitable value of μ is required. In this work, the amplitude of the input signal varies with time. To achieve good performance, μ needs to be adapted to the input. Therefore, an implementation of the normalized LMS algorithm was used by varying the parameter μ as
where β is a constant parameter and is an estimate of the signal power at instant n, which is approximated as . The length L should be sufficiently small to capture input amplitude changes, and sufficiently large to provide a good local estimate of the signal power. To summarise, an estimate of the input signal frequency, starting from some value is updated on a sampling interval basis using (4). The output, y[n], is computed directly from (2), while the same equationcan be differentiated and evaluated to obtain an estimate of the derivative . Appropriate choice of the tracking parameter β will ensure that a normalized LMS (nLMS) algorithm will converge to a correct value of up to a small misadjustment parameter (recall also that there is a fundamental tradeoff between the speed of tracking and the root mean square error in the estimate of the respiratory frequency).
In practice, the idealized signal model above suffers from a number of practical challenges. The real signals do not consist of a single strong purely sinusoidal component corrupted by a small amount of uncorrelated noise. Rather, the signal received is a combination of a respiration signal and bodily movement signal (see figure2 as example), together with a small amount of intrinsic system noise (internal electronic circuit noise plus background mechanical noise in the environment). Moreover, the respiration signal is seldom sinusoidal in nature; rather it has multiple harmonic components. There will also be periods when the subject is absent (e.g. bathroom break) and in those situations, the algorithm may take on arbitrary values dependent on the random noise signal presented. Also, as noted above, the system outputs two quadrature signals, whose relative strength will depend on the relative distance for the subject and the particular reflection pattern they observe, so either the I or Q signal is often significantly higher in SNR than the other. So in practice, the value for may be attracted to a local minimum, in particular for the case of a breathing signal with strong harmonic components. Because of all these reasons and since IIR adaptive filters are slow in convergence, it is a good idea to try to initialize the algorithm as close as possible to the correct BR value, such that it converges quickly. Since the change in BR is slow, BR tracking can be achieved without problems following good initialization. We propose that, the initial value be derived from a high quality data window that is not dominated by noise or artefacts. The calculation of the initial value from such data window can be carried out using equation(3). This process is applied once and should not affect the computational simplicity of the tracking algorithm. To overcome these problems, we need to ensure that in parallel with the ANF algorithm, we must also (a) recognize whether a valid biomotion signal is present, (b) select the optimal channel between I and Q, and (c) ensure that the nLMS is operating close to the true respiration frequency as much as possible.
The first step, is to ensure that a signal with sufficient signal-to-noise ratio is presented to the ANF algorithm. The SNR may be low because there is low amplitude signal (i.e. the subject is at the range of detection of the sensor, or presents an unusually low reflected signal to the sensor), or conversely the noise may be high (due to lots of body movement, external mechanical movement, etc.). We have developed a simple signal selection algorithm as follows:
A data segment of length is chosen (in our case we will typically use TW = 15 s). The data segment is zero-meaned, and the following procedure is implemented on the data segment:
(a)
Find the sign (+/−) of each of the data samples.
(b)
Determine zero-crossings as the points where sign changes.
(c)
Calculate the time durations between consecutive zero-crossing pairs.
(d)
Ignore/remove durations smaller than a certain threshold (e.g. durations less than 0.2 s)
(e)
Combine pairs of durations to form new larger durations, . (where is the number of cycles)
(f)
Calculate the signal power associated with each new duration,
The above procedure is used to identify 'clean' respiration signals. Such signals are characterized by their relatively regular and small number of zero-crossings compared to the noise and movement signals. Clean breathing signals also tend to have approximately equal breathing rate power distribution between consecutive cycles. If the procedure above is executed on a clean respiration segment, the (corresponding to breathing cycles) produced in step 5 will have a very small variation over different i values, provided that the data window is sufficiently small. The same applies for obtained from step 6. The purpose of step 4 is to suppress the effect of any small disturbances around the zero-crossing points. From the above procedure, a number of parameters will be available that can be utilized as signal quality measures. From step 2, the number of zero-crossings () is directly obtained. The outcomes of step 5 and step 6 are used to calculate the coefficients of variation and of and respectively (the coefficient of variation is the standard deviation divided by the mean—the coefficient of variation will be close to zero for a very regular signal). Therefore, the smaller and the more likely the data is due to a clean respiration. In addition to these parameters, the average signal power is used as another parameter for quality evaluation. To detect the presence of a signal suitable for the ANF algorithm, a number of tests are performed. The first test is to detect the presence of a bio-signal in the data segment under test. For this purpose a certain power threshold, has to be exceeded. Data whose power does not exceed the threshold, is considered due to background noise—and absence of any bio-information is declared (i.e. no subject is in range). Since breathing signal power is quite variable, finding a strict value for is difficult. In practice, we use quite a low value of and an additional metric is applied for bio-signal/noise discrimination. The number of zero-crossings, is utilised for this purpose. It was mentioned previously that data consisting of background noise tends to cross the zero level more times than when bio-signals are present. A very small number of zero-crossings can also indicate absence of useful respiration data due to a large DC shift. Therefore, two thresholds, and are used to confirm the presence/absence of a bio-signal. The detection of a useful bio-signal is declared if:
and
If detection is declared, the data window needs to pass additional tests to qualify for being used for computation of the initial value of the ANF, . A data window is adequate for calculating if, in addition to satisfying (8), it satisfies
where , and are appropriately chosen thresholds. The operation involving in (8) is to ensure that the value is within the appropriate range determined by the practical BR range (e.g. for 15 segment data windows a range of would correspond to an allowed breathing rate range of 8–32 breaths min−1. The final operation in (8) is to ensure that the detected cycles occupies a sufficient proportion of the window. A data window that satisfies the criteria in (8) is classified adequate (for the initialization purpose). In the case where both I and Q signals are adequate, we select the optimal channel by choosing the one with the lower number of zero crossings, and in the event of a tie, select the channel with the lowest value of . Contrary to a periodogram approach, which processes a large data block each time a BR estimate is produced, the ANF tracking algorithm updates its output by processing only three consecutive samples of the input data (see equations(2) and (4)). Tracking is achieved by sliding the three-sample block one sample at a time to update the previous output This feature of the ANF approach enables it to capture fine changes in breathing rate.
Prior to analysis of the nocturnal respiration rates of the COPD subjects, we conducted a validation exercise on previously acquired clinical data to confirm that the ANF algorithm was producing accurate results. In a previously published study (Zaffaroni et al 2013), SleepMinder was used to collect biomotion data of patients who were referred to the St. Vincent's University Hospital Sleep Disorders Clinic, Dublin, Ireland. In addition to SleepMinder signals, full polysomnography (PSG) was performed using the Jaeger-Toennies 1000 E system. Movements of the ribcage and abdomen were monitored using respiratory inductance plethysmography sensors around these regions, as is standard in PSG. Ten full-night recordings (each of duration longer than 6 h) were selected from this data set. A 15 min sectionwas selected from each recorded night, based on the requirement that it is free of significant movement and disordered breathing to benchmark the ANF algorithm against a human scorer. The human scorer used a visualization tool to annotate the 15 min sectionswith end-points of inspiration and expiration (for both the chest and abdomen channels). An expert derived breathing rate was generated on a 1 s basis using data blocks of 15 s that were shifted every 1 s. Expert results were manually scored by an experienced respiratory scientist. Final expert results were formed by averaging the results obtained from each of the chest and abdomen data (see Shouldice et al 2010 for details). The data were analyzed using Matlab (Mathworks, Natick, MA). The system used the value T = 1 s, and a data segment length of 15 s. The ANF algorithm was applied separately to the I and Q channels. A combined estimate was formed by averaging the ANF-I and ANF-Q respiration rate, to form an ANF-C.