Wireless fetal monitoring system
09717412 · 2017-08-01
Assignee
Inventors
- Masoud Roham (San Diego, CA, US)
- Enrique Saldivar (Santee, CA, US)
- Srinivas Raghavan (San Diego, CA, US)
- Mehran Mehregany (San Diego, CA, US)
- Mitul Shah (San Diego, CA, US)
Cpc classification
A61B5/02438
HUMAN NECESSITIES
A61B5/0245
HUMAN NECESSITIES
A61B8/0866
HUMAN NECESSITIES
A61B5/0022
HUMAN NECESSITIES
International classification
A61B8/00
HUMAN NECESSITIES
A61B5/00
HUMAN NECESSITIES
Abstract
A wireless fetal and maternal monitoring system includes a fetal sensor unit adapted to receive signals indicative of a fetal heartbeat, the sensor optionally utilizing a Doppler ultrasound sensor. A short-range transmission unit sends the signals indicative of fetal heartbeat to a gateway unit, either directly or via an auxiliary communications unit, in which case the electrical coupling between the short-range transmission unit and the auxiliary communications unit is via a wired connection. The system includes a contraction actuator actuatable upon a maternal uterine contraction, which optionally is a EMG sensor. A gateway device provides for data visualization and data securitization. The gateway device provides for remote transmission of information through a data communication network. A server adapted to receive the information from the gateway device serves to store and process the data, and an interface system to permits remote patient monitoring.
Claims
1. A wireless fetal and maternal monitoring system comprising: a fetal sensor unit adapted to provide signals indicative of a fetal heartbeat, a contraction sensor adapted to provide signals indicative of a maternal uterine contraction upon a maternal uterine contraction, a digitization and control unit, the digitization and control unit adapted to receive the signals indicative of a fetal heart rate and the signals indicative of a maternal uterine contraction, the digitization and control unit including a heart rate calculation system to calculate the fetal heart rate, a fetal heart rate register for storing the calculated fetal heart rate, a microcontroller configured to refresh per fetal heartbeat the fetal heart rate register, a maternal uterine contraction register to store the signals indicative of maternal uterine contraction, the digitization and control unit providing data fusion by transmitting at least the calculated fetal heart rate from the fetal heart rate register and data based on the signals indicative of the maternal uterine contraction from the maternal uterine contraction register; a short-range transmission unit adapted to receive from the digitization and control unit the fused data signals and to retransmit the signals, a microphone comprising a maternal heartbeat sensor, and a processor including a comparator and a flag memory to eliminate as an erroneous measure a first sensed beat period by comparison to a second sensed beat period and to set a flag in the flag memory.
2. The wireless fetal and maternal monitoring system of claim 1 wherein the fetal sensor unit includes an ultrasound sensor.
3. The wireless fetal and maternal monitoring system of claim 2 wherein the ultrasound sensor is a Doppler ultrasound sensor.
4. The wireless fetal and maternal monitoring system of claim 1 wherein the fetal sensor unit is a sound sensor.
5. The wireless fetal and maternal monitoring system of claim 4 wherein the sound sensor includes a sensor array.
6. The wireless fetal and maternal monitoring system of claim 5 wherein the sensor array is a MEMS microphone array.
7. The wireless and maternal monitoring system of claim 5 wherein the sensor array is a directional sensor array.
8. The wireless fetal and maternal monitoring system of claim 1 wherein the fetal sensor unit includes an ECG sensor.
9. The wireless fetal and maternal monitoring system of claim 1 wherein the short range transmission unit utilizes the BLUETOOTH® standard.
10. The wireless fetal and maternal monitoring system of claim 1 wherein the short range transmission unit operates at a specific absorption rate (SAR) of less than or equal to 0.1 watts/kg.
11. The wireless fetal and maternal monitoring system of claim 1 wherein the short range transmission unit operates at a specific absorption rate (SAR) of less than or equal to 0.05 watts/kg.
12. The wireless fetal and maternal monitoring system of claim 1 wherein the short range transmission unit operates at a specific absorption rate (SAR) of less than or equal to 0.01 watts/kg.
13. The wireless fetal and maternal monitoring system of claim 1 further including an auxiliary communication receiver and transmitter adapted for wearing by the maternal user for receiving transmissions from the short-range transmission.
14. The wireless fetal and Maternal monitoring system of claim 13 wherein the transmissions between the short range transmission unit and the auxiliary communication unit are via a wired connection.
15. The wireless fetal and maternal monitoring system of claim 13 wherein the transmissions between the short range transmission unit and the auxiliary communication unit is wireless.
16. The wireless fetal and maternal monitoring system of claim 13 wherein the auxiliary communication unit is adapted to be worn on a neck.
17. The wireless fetal and maternal monitoring system of claim 1 wherein the contraction sensor is a tocodynamometer sensor.
18. The wireless fetal and maternal monitoring system of claim 1 wherein the contraction sensor includes a button.
19. The wireless fetal and maternal monitoring system of claim 1 wherein short-range transmission unit utilizes the Class 3 BLUETOOTH® standard.
20. The wireless fetal and maternal monitoring system of claim 1 further including a pack adapted to contain one or more of the wireless fetal and maternal monitoring system, a blood pressure device, a glucometer, a reagent dip stick, a cell phone gateway and a wireless hub.
21. The wireless fetal and maternal monitoring system of claim 1 further including a motion sensor.
22. The wireless fetal and maternal monitoring system of claim 1 further including a sensor unit adapted to provide maternal heart rate data.
23. The wireless fetal and maternal monitoring system of claim 1 wherein the data fusion includes data based on signals indicative of a second fetal heart rate.
24. The wireless fetal and maternal monitoring system of claims 1 further including a gateway device, the gateway device including a display to provide for data visualization.
25. The method for measuring fetal and maternal physiological parameters of claim 1 wherein the comparison is of a heartbeat period with a previously stored heartbeat period.
26. The method for measuring fetal and maternal physiological parameters of claim 1 wherein the signal is eliminated if the measured heartbeat period differs from the stored heartbeat period by more than +/−25%.
27. The method for measuring, fetal and maternal physiological parameters of claim 1 wherein the previously stored heart rate is a fetal heartrate.
28. The wireless fetal and maternal monitoring system of claim 1 including a pulse oximeter.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9)
(10)
(11)
(12)
(13)
(14)
(15)
(16)
(17)
(18)
(19)
(20)
(21)
(22)
(23)
(24)
DETAILED DESCRIPTION OF THE INVENTION
(25)
(26) The central unit includes a short range communication module. A gateway 20 is used for local data storage, visualization and to communicate with the mobile data network to transmit data to the server. The short range communication is employed for safety considerations so that radio frequency (RF) emission with high power (that is required for communication to the cellular network) gateway 20 is placed relatively far from the mother/baby 22. The short range wireless communication module implemented in the central unit 10 has low power RF emission thus it is very less likely to be harmful. The short-range transmission unit is a low power transmission unit, preferably having specific absorption rate (SAR) of less than or equal to 0.1 watts/kg, and more preferably less than 0.05 watts/kg, and most preferably less than or equal to 0.01 watts/kg. This level of SAR is implemented as known to those skilled in the art, such as through the use of BLUETOOTH® technology. Preferably class 3 BLUETOOTH® technology or otherwise the lowest radiation class is utilized. Optionally, radiofrequency shielding is utilized.
(27) One significant advantage of using a gateway in conjunction with the short range body sensor wireless link to the device against direct link from body worn sensor to mobile data network is reducing fetus and pregnant mother exposure to the RF radiation of wireless fetal monitor.
(28) Both wireless gateway and BLUETOOTH® module emit non-ionizing radiation at frequencies ranging in 1-2.5 GHz. The FCC limit on the Specific Absorption Rate (SAR), a measure of the rate of energy absorption by the body when exposed to an RF field (see, e.g., C. K. Choul, et al, “Radio Frequency Electromagnetic Exposure: Tutorial Review on Experimental Dosiinetry”, Bioelectro-magnetics, Vol. 17, Issue 3, pages 195-208 (1996)), for cellular telephones is 1.6 W/kg.
(29) The SAR rate of the gateway is comparable to typical smart phones, in the range of 0.5-1.5 W/kg (see, e.g., Electromagnetic Fields and Public Health: Mobile Phones“, World Health Organization, Fact Sheet No 193, May 2010) A BLUETOOTH® radio module configured in class II generates a SAR level of ˜0.01 W/kg. Therefore, by utilizing a gateway, placed relatively far from the pregnant woman the SAR level can be reduced by two orders-of-magnitude and well below FCC standards. Using the optional external BLUETOOTH® necklace, rather than the built-in module, further diminishes the undesired RF emission exposure to the fetus to an even less significant value.
(30) In order to eliminate any concern regarding absorption of radio frequency signals by the fetus, an auxiliary communication unit 16 is optionally utilized. In one variation, the auxiliary communication unit is in the form of a necklace, which locates the transmitter to the gateway 20 at a significant distance, such as at least two feet (though this distance will vary based on the height and physical structure of the mother) from the fetus. In this implementation, the communication from the central unit 10 to the auxiliary communication unit 16 may be wireless, but is preferably wired via connection 18. The wired, i.e., not wireless, communication from the central unit 10 minimizes radiation to the fetus.
(31)
(32) Separate belts are preferably used to hold the central unit and toco sensor so that during operation, position of sensors can be independently optimized. The central unit includes ultrasound transducers as well as control, processing and BLUETOOTH® communication circuitry. A toco pressure sensor, an optional audio feedback earphone and the optional external BLUETOOTH® necklace can be plugged in to the central unit.
(33)
(34)
(35)
(36) One particular implementation of the sensing hardware is described with reference to
(37) Due to motion artifacts and/or inappropriate positioning of transducers on a mothers abdomen, the heartbeat detector often misses one or more heartbeats. An algorithm for heartbeat to heart rate conversion, embedded on microcontroller, eliminates the erroneous measure via comparing input beat period with the previously stored value. In case that current reading is outside of ±25% of the stored value, the algorithm drops the new reading and raises a flag. If 6 consecutive readings are constantly out of that range the new reading is stored as updated measurement result.
(38) A low-cost disposable toco sensor (FeatherLiteToco, Ventrex) which consists of a pressure transducer configured in a Wheatstone bridge is used for contraction monitoring. An instrumentation amplifier with a gain of 100 amplifies the signal to the ADC input range. Further baseline subtraction and gain adjustment is implemented in the gateway software. The device makes an authenticated link with the gateway using a BLUETOOTH® module (RN-41, Roving Networks) configured in Serial Port Profile. The module's output RF power can be programmed for either class I, II or III. An optional external BLUETOOTH®, in a necklace form factor is designed so when it is plugged in to the unit, substitutes the internal BLUETOOTH®.
(39) Current consumption of the module is dominated by electronics driving ultrasound crystals and the BLUETOOTH® module, measured at 60 mA and 25 mA, respectively from the 3.3V regulated supply. The device is powered by two standard AAA batteries which results in approximately 8 hours of constant running time. By powering from separate up-converting voltage regulators, interference between the sensing interface electronics and the BLUETOOTH® module is minimized.
(40)
(41)
(42) The EMG signals occupies a different band in frequency and could be filtered out from other signals and used for tracking uterine contractions. The ultrasound transducers preferably are arranged as an array that enables the electronics process the signal to minimize the need to repositioning of the patch due to baby movements. This arrangement is described in co-pending U.S. Provisional Patent Application Ser. No. 61/327,975, entitled “Ultrasound Patch”, filed Apr. 26, 2010, incorporated herein by reference as if fully set forth herein.
(43)
(44) In a Doppler ultrasound, the measured shift in the frequency/phase of the received signal in comparison to the transmitted signal is of interest, even though it may be very small. This method is called continuous-wave (CW) Doppler, where the change in frequency and phase of the reflected ultrasound signal is measured. This technique is different from the traditional sonographic techniques and does not be used to create an image, but rather to measure the fetal heart rate, and optionally other parameters such as flow in blood vessels, veins, and arteries.
(45) Control circuitry is coupled to the transmission system and the receiver system. The control system may include analytical or analysis functions. A processor may be provided, either within the patch, or external to the patch, to perform analytical or analysis functions.
(46) In this patch embodiment, in addition to sensors for fetal heartbeat monitoring, dry electrodes are provided to record bio-potentials such as electromyogram (EMG).
(47) In one implementation of the patch, the wearable patch for use on a body is in the form of a planar pad. The preferred dimensions of the patch are 80 mm×25 mm and thickness 5 mm or less, and most preferably 60 mm×20 mm 3.5 mm or less. The patch should be light-weight, about 16 grams or preferably weighing 8 grams or less.
(48) The following detailed description has applicability to systems for multiple births, but also has general applicability for systems and methods for single births. The fetal monitoring device consists of two components, illustrated in
(49) Data is passed serially from the first (nearest toco) monitor in the daisy chain to the last. In all configurations, data is transmitted from the last FHR monitor in the chain to a nearby cellular gateway using a BLUETOOTH® communication module. In addition to the nominal (c) and twin (d) configurations shown in the figure, the FHR monitor may be used stand-alone (without toco sensor), or a 3rd FHR monitor may be connected at the end of the chain to be used as an off-body transmitter for a twin configuration.
(50) When fitted on the mother, the device would appear approximately as shown in FIG. 8. The architecture employed in the design of this fetal monitoring device could support any number of births, but it may be impractical to fit the monitor for more than twin births. By providing twin FHR monitors, the monitoring time of a mother can be cut in half.
(51) The device has been developed with usability in mind. The user must simply plug in components in order to activate power and data collection. The FHR monitor automatically detects the presence or lack of a connection, and its type. LED indicators illuminate to inform the user of the monitor status: green for a valid input connection, blue to signify that the BLUETOOTH® transmitter is operating, and flashing amber for the heartbeat.
(52) Additional features simplify the fitting procedure. When the toco belt is tightened, the green indicator flashes to signify that contraction threshold has been exceeded, and ceases to flash when the belt is loosened to produce strain below a slightly lower threshold. Also, the demodulated analog output from the Doppler signal processing employed by the FHR monitor is buffered and provided to a stereo audio jack so both mother and practitioner can listen to the sound of the heartbeat during fitting.
(53) To prevent data loss in the event that communication is lost during a monitoring session, the FHR monitor includes a back-up memory, by way of example a 4.5-hour backup memory, which can be implemented using a 1-Mbit serial EEPROM that is written and read using a SPI interface running at 1.5 Mbps. When the memory backup feature is enabled, each data packet that is transmitted to the serial daisy chain or to the BLUETOOTH® transceiver, is also written to the EEPROM.
(54) The backup memory is set up via the cellular gateway. During memory setup, the blue light flashes to indicate that data collection is suspended. A simple command language has devised in which an initial receipt of “M” by the BLUETOOTH® module switches operation from normal (acquisition) mode to memory mode, in which received characters are processed as commands to enable/disable memory backup, to set the memory address, and to upload data from the memory.
(55) As to possible circuit design, the FHR monitor is controlled using a low-cost, 8-bit microcontroller that includes all the analog-to-digital conversion, timing, and indicator drive required by the monitor, as illustrated in
(56) A piezoresistive Wheatstone bridge toco sensor is connected between the VsB and RET pins of the input connector, with its differential sensor output connected to the Vs+/Vs− pins. Alternatively, an FHR monitor may be connected to the input port, in which case the serial data output Tx1 connects to the serial data input Rx1, and the supply voltage VDD is used to bias the INA inputs.
(57) As illustrated in
(58) Upon start-up, the FHR monitor must determine what type of device is connected to its input port, i.e., a toco sensor, an FHR monitor, or a simple powering plug with no associated sensing device. This is accomplished through a combination of pull-up and pull-down resistors of the appropriate ratios (not shown), in addition to logic in the firmware of the embedded microprocessor. As was shown in
(59) There are only two types of output connections that must he detected, i.e., a twin FHR monitor or no connection. This is accomplished by providing a pull-down resistor on the serial data transmit line (Tx1). If no device is connected to the output connector, the logic level is pulled low. When an FHR monitor is plugged into the output connector, the pull-up resistor on its Rx input, having a 10× smaller value, results in a high logic level. Since output connections may be made or broken after start-up, this connection must be tested each time data is to be transmitted. If an FHR monitor is detected, the internal BLUETOOTH® module is disabled and data is sent to the serial daisy chain. If no connection is sensed, the data is sent to the BLUETOOTH® transmitter.
(60) For audio signal processing, the device preferably uses a precision 2.0-MHz sinusoid is derived from the 12-MHz master clock, and buffered to drive the transmitting ultrasonic transducer. The signal from the receiving transducer is first amplified using a tuned, JFET common-source amplifier, then demodulated using a chopping mixer. The baseband signal is then passed through a four-stage band-pass amplifier that passes the Doppler-shifted signal in the frequency range of 100-500 Hz. This audio signal is amplified using a PGA and input to the ADCs, and is also buffered to drive a stereo ear-piece. The total voltage gain may be varied from 64 dB to 106 dB.
(61) The differential input from the toco sensor is simply amplified by 46 dB using an instrumentation amplifier (INA), then input to its ADC and averaged over 120 samples (a half second) in the microprocessor. Additional baseline subtraction and gain adjustment is implemented in the gateway software, and as part of the fitting calibration procedure.
(62) For digital signal processing, the FHR is calculated using a robust algorithm that is based on autocorrelation, described in more detail, below. Given the requirement of a minimum FHR of 30 beats per minute (BPM), the autocorrelation window must be 2 seconds in duration. A preliminary examination of typical Doppler signals revealed that the 100-500 Hz signal (
(63) As to data format and daisy chain communication, the serial data chain could be extended indefinitely. The digital signals that originate with the first FHR monitor in the chain, i.e. the “primary”, are transmitted serially using RS-232 format. The toco sample would be dropped into the beginning of a data packet, and the value 0 could be used as a marker to indicate that the toco sensor is not present, as in stand-alone FHR monitoring. The primary FHR monitor would drop its FHR data into the next slot and marks all other slots in the packet as empty. Any additional FHR monitors in the chain would recognize that they are not the primary and would instead drop their FUR data into the first empty slot, then pass it up the chain. The final FHR monitor in the chain would transmit the data using its BLUETOOTH® module.
(64) While the concept could be extended indefinitely, it is limited by the chosen packet size. In the present implementation, illustrated in
(65) Including start and stop bits, a data packet consists of 40 bits, which is transmitted in 16.7 msec at 2400 bps. When the BLUETOOTH® module is enabled, data is transmitted wirelessly upon a query (“Q”) received from the module. When the wireless module is disabled, data is sent to the daisy chain Tx1 following each packed received from Rx1, or at regular update intervals (each half second in the current implementation). Since a twin FHR monitor could be connected/disconnected to/from Tx1 at any time after startup-up, the device is programmed to test the output connection before transmission of each data packet, which requires that the serial port circuits be temporarily disabled, then re-enabled prior to transmission.
(66) As illustrated in
(67) By ensuring synchronicity of timing between ADC samples and serial communication, interference from the communications circuitry can be minimized, as illustrated in
(68) Circuits may be fabricated on a printed circuit board (PCB) having dimension 115.5 mm by 95.0 mm, for ease of debug and test, then laid out for the final size and form factor, a double-side, oval PCB having dimensions 85.4 mm by 66.6 mm, of which 1480 mm2 are occupied by the rechargeable, lithium-polymer battery. Photographs of the assembled device components are provided in
(69) The test results were obtained using the assembled FHR monitor when possible, and from the PCB with increased form factor, when necessary. The schematic designs of the circuitry are equivalent in the two versions. Wireless sensor data was captured using the BLUETOOTH® transceiver of a laptop computer.
(70) Total measured current draw from the rechargeable, 4.2-V lithium-polymer battery is 112 mA, where 60 mA is drawn by the BLUETOOTH® module, and 13 mA is drawn by the transducer drive circuitry. The unit may therefore operate for almost 8 hrs before recharging the 850-mA battery.
(71) An overview of the signal processing algorithms is described. The digital signal processing approach for calculating fetal heart rate (FHR) from Doppler signal has several steps as illustrated in
(72) Preprocessing includes sampling, envelope detection and downsampling. Performance of the envelope detection depends on how accurately peaks on the Doppler signal are sampled. Therefore, the sampling frequency needs to be high enough to accurately sample peaks in the signal while maintaining the minimum requirement of satisfying the Nyquist criterion. The Doppler signal is sampled at fs=2400 sps to guarantee a precision of 92% in detecting peaks, given a nominal Doppler shift of fd=300 Hz.
(73) The sampled signal is passed through an envelope detection algorithm which detects the positive envelope of the signal. The envelope is then downsampled by a factor of 10, reducing the rate of input data to the autocorrelation algorithm to 240 sps, a sample rate adequate to track the nominal 20 Hz frequency of the envelope to a precision of 96%.
(74) Using autocorrelation, repetitive patterns are found from the Doppler ultrasound signals, and heart rate values are calculated according to the period of peaks in the autocorrelation results. Autocorrelation is a mathematical function that measures the similarity between different segments of a time series signal as a function of time-shift between the segments. Auto-correlation of a signal xt over a window of length W is given by
(75)
and is calculated for different values of time lag, τ. Window size is chosen in this work to be 480 samples to ensure that 2 seconds of Doppler data is considered in the autocorrelation calculation, permitting a minimum detectable heart rate of 30 bpm. While a normal fetal heart rate ranges between 110 and 160 bpm, abnormal rates can be as low as 30 bpm or as high as 240. Therefore, the window size used in autocorrelation algorithm needs to be long enough to accommodate at least one heart beat. Furthermore, the window is moved forward over the signal to find repeating patterns. The location of the repeating heart beats appear as peaks in the autocorrelation results which help in finding the duration and subsequently frequency of the heart rates. Thus, the window needs to be moved for a sufficiently long period of time in order to ensure that at least two repetitions of the slowest heart beat (30 bpm) appear in the autocorrelation data. Therefore, the autocorrelation is calculated for τ ranging from 1 to 480.
(76) Occurrence of repeating pattern in the original signal is manifested in the peaks of the autocorrelation results as shown in
(77) For preprocessing, the Doppler signal is sampled at 2400 sps and downsampled to 240 sps for input to the auto-correlation block. While particular design parameters are set forth herein, the particular design parameters may be set by those skilled in the art to achieve the functionality and operations of the inventions described herein.
(78) The choice of sampling frequency relies on two criteria that need to be met: 1) the sampling rate needs to be high enough to satisfy the Nyquist criterion, 2) it needs to be sufficiently high in order to precisely detect peaks of the Doppler signal, which will form the envelope of the signal in subsequent processing block. Studies have shown that in applications of Doppler ultra-sound for fetal heart rate monitoring, the Doppler-shifted signals in the range of 100 to 500 Hz are associated with the baby's heart movements. Therefore, any sampling frequency above 2×500 would satisfy the Nyquist criterion. In other words, fs≧1000.
(79) In order to explore the second criterion for sampling frequency, the peaks of the Doppler signals approximate a sinusoid of period 2T.sub.peak as shown in
(80)
(81) Thus, for a given value of precision, P, T.sub.sample can be calculated by
(82)
(83)
(84) As to the downsample rate, the preceding approach may be used. The input to the downsampling block is the envelope of the Doppler signal. Experimental data collected from real subjects shows that peaks on the envelope signal has a frequency range between 5 and 20 Hz. Choosing a sampling rate of 240 Hz for downsampled signal gives a downsampling rate of 10. The sampling rate of 240 is adequate to track the nominal 20 Hz frequency of the envelope to a precision of 96% as shown in
(85) Architecture of the autocorrelation block is illustrated in
(86) It is a semi-systolic array architecture with the main processing cells being Multiply-ACcumulate (MAC) units that hold partial autocorrelation results. The architecture is composed of 3 register arrays: envelope (top row), partial autocorrelation (middle row), and final autocorrelation (bottom row), each of which has a length of L associated with the maximum lag of τ. Each column of this architecture corresponds to the autocorrelation calculation for a specific τ. For example, the first column calculates autocorrelation for delay of τ=1, second column for delay of τ=2, etc. As shown in the figure, the maximum lag is L samples, which is considered to be L=480 for the experiments as discussed here.
(87) The envelope array stores the last L samples (2 seconds) from envelope and downsampling blocks. Each new downsampled data (x.sub.i) is multiplied by each sample in the envelope array and is added to a corresponding location in the partial autocorrelation array. The window size is controlled by the R.sub.i control signals which are activated sequentially (R.sub.1, R.sub.2, . . . R.sub.W, R.sub.1, . . . ). One element of the partial autocorrelation array is copied to its final autocorrelation location for each envelope sample, and partial autocorrelation cell is reset or ‘0’ in preparation for the next series of MAC operations. Given that the length of the autocorrelation array is L=480, each element of the partial correlation is copied/cleared every 2 seconds.
(88) A graphic of autocorrelation array (τ={(1), τ(2), . . . , τ(L)}) is shown in
(89) A weighted median approach is used to detect the center of each peak, in other words each peak is said to occur at the weighted median of all autocorrelation samples that exceed a certain threshold, as illustrated in the figure. Potentially, there might be more than one peak in the autocorrelation data. The time interval of the heart rate is thus calculated using
(90)
where τ(τ) refers to autocorrelation with lag τ, M=τ(1) denotes the autocorrelation value at τ=1 and ε is a parameter that specifies the threshold for peak detection, and N.sub.k denotes the peak number, i.e. the peak at 360 samples in the example is the 3rd peak. Peak number is identified by upward and downward threshold crossings, as indicated in the figure.
(91) As suggested by the above equation, the heart rate calculation requires scanning through the entire autocorrelation array of L elements. In the present implementation, one sample of the autocorrelation array (τ) is analyzed during the 2400-sps interrupts, so the entire τ array is scanned at a rate of 5 times per second, yielding an updated heart rate calculation every 0:2 seconds. Since 2 seconds are required to update the complete array, each heart rate calculation will be based on 10% “new” τ data and 90% “old” τ data, which provides a low-pass filter of sorts on the calculated heart rate.
(92) The algorithm for calculating heart rate from the autocorrelation data requires L iterations to complete a full scan of the autocorrelation array. At each iteration of the algorithm, the first autocorrelation value (τ(1)) is read and used to set the threshold since τ(1) will always have maximum correlation given that τ(0) is not calculated. A ‘peak’ is defined as a span of autocorrelation data that exceeds threshold thr=M(1=ε). Within each peak, a summation (S=Σ.sub.τr(τ)) and weighted summation (WS=τ×Σ.sub.τr(τ)) are calculated, as required for calculation of T.sub.peak, the weighted median of the peak.
(93) Several tests are then performed to test the validity of the peak. For example, a peak cannot occur less than ¼ of a second after the previous peak or from the beginning of the array, since the device is not sensitive to heart rates >240 bpm. The width of the peak should also exceed a minimum threshold, currently set to 3 samples—such false peaks can be the result of a noisy input. Finally, each valid peak is used to calculate the overall summation (S) and overall weighted summation (WS), which is to calculate T.sub.HB, and HeartRate (bpm)=60 s/m×240 sps/T.sub.HB.
(94) This formulation is equivalent to calculating the weighted median of each peak, T.sub.peak, and then calculating the average of T.sub.peak=N.sub.k, weighted by their strength S. The most memory consuming blocks in the signal processing pipeline include envelope detection, partial autocorrelation calculation, and final autocorrelation calculation as depicted in
(95) TABLE-US-00001 TABLE 1 Memory requirement of the algorithm Memory Array Length Unit Size (bits) Memory Usage (bytes) Envelope L = 480 7 420 Partial Autocorrelation L = 480 24 1440 Final Autocorrelation L = 480 24 1440 Total 3300
(96) MAC operations and τ data are as follows. As shown in Table 1, a total of 3300 bytes memory suffice to accommodate the entire autocorrelation calculation results. The amount of memory required for other computing blocks such as heart rate calculation algorithm described compared to the aforementioned storage is negligible.
(97) For time complexity, major operations which are needed for calculation of the final autocorrelation results are analyzed here. Table 2 shows the number of operations including multiplication, addition (24 bits) and register transfers (move) required to update each one of the
(98) TABLE-US-00002 TABLE 2 Number of instructions to process one envelope sample for the purpose of autocorrelation calculation Updated Array #Mult #Add24 #Mov Envelope 0 0 L Partial Autocorrelation L L 0 Final Autocorrelation 0 0 1 #Instructions L 6L L + 1 Total 8L + 1
arrays during autocorrelation calculations. The envelope array needs L number of register transfer operations in order to shift elements of envelope array upon receiving a new envelope data. Calculation of partial autocorrelation needs multiply-add operations as discussed before. Finally, only one element of final autocorrelation array is updated (transferring results from partial array) when a new envelope data is received. Given that envelope data are generated at a rate of 240 sps, the algorithm requires 921840=240×(8×480+1) instruction per second assuming a length of 480 for each one of the arrays.
(99) While described primarily herein for applications in fetal monitoring, as will be appreciated by those skilled in the art, the applications are much broader. The autocorrelation-based approach for estimating frequency of repeating patterns can be used for a variety of applications in addition to the Doppler ultrasound signal processing. In particular, this technique can be used to measure heart rates from ECG signals, gait parameters such as step rates from motion sensors, and respiration rate from photoplethysmograph (PPG) sensors.
(100)
(101)
(102) Various optional blood measurement systems are provided within the system or kit. Blood glucose strips are optionally included. If included, a blood draw tool such as a Lancet, holder, and sharps disposal unit are provided. If ultrasound is to be used, ultrasound gel is optionally provided. Preferably, blood pressure measuring apparatus is provided, including a blood pressure cuff (shown in the center of
(103)
(104)
(105) In designs described above, monitoring of electrical activity on the mother's belly could be used for detection of fetal heart beat. This technique eliminates the need for ultrasound transducer and it is less sensitive to the positioning of the device. It is a passive technique, meaning that, unlike ultrasound, the device does not emit any signal for heart beat detection, thus is suitable for continues monitoring.
(106) The other technique for monitoring of the fetal heart beat is using MEMS microphones or microphone arrays to detect the sound of fetal heart. See, e.g., R. R. Lahiji, M. Mehregany, “Microphone Arrays for Listening to Internal Organs of the Body”, U.S. Provisional Patent No. 61/258,082, filed November 2009, now published as U.S. Publication 2011-0137209, incorporated herein by reference as if fully set forth herein. This is a passive technique and is suitable for continues monitoring.
(107) Optionally, a manual entry is provided for recording contractions instead of, or in addition to, a toco or uterine EMG recording. If the mother herself senses the uterine contractions and presses a button or actuator on the gateway to record the contraction happening.
(108) All publications and patents cited in this specification are herein incorporated by reference as if each individual publication or patent application were specifically and individually indicated to be incorporated by reference. Although the foregoing invention has been described in some detail by way of illustration and example for purposes of clarity and understanding, it may be readily apparent to those of ordinary skill in the art in light of the teachings of this invention that certain changes and modifications may be made thereto without departing from the spirit or scope of the following claims.