Darien Figueroa / Mbed 2 deprecated repo3

repo time

Dependencies:   mbed MAX14720 MAX30205 USBDevice

HspGuiSourceV301/HSPGui/Presenter/RToRCalculator.cs@20:6d2af70c92ab, 2021-04-06 (annotated)

Committer:
darienf
Date:
Tue Apr 06 06:41:40 2021 +0000
Revision:
20:6d2af70c92ab
another repo

Who changed what in which revision?

UserRevisionLine numberNew contents of line
darienf 20:6d2af70c92ab 1 /*******************************************************************************
darienf 20:6d2af70c92ab 2 * Copyright (C) 2016 Maxim Integrated Products, Inc., All rights Reserved.
darienf 20:6d2af70c92ab 3 *
darienf 20:6d2af70c92ab 4 * This software is protected by copyright laws of the United States and
darienf 20:6d2af70c92ab 5 * of foreign countries. This material may also be protected by patent laws
darienf 20:6d2af70c92ab 6 * and technology transfer regulations of the United States and of foreign
darienf 20:6d2af70c92ab 7 * countries. This software is furnished under a license agreement and/or a
darienf 20:6d2af70c92ab 8 * nondisclosure agreement and may only be used or reproduced in accordance
darienf 20:6d2af70c92ab 9 * with the terms of those agreements. Dissemination of this information to
darienf 20:6d2af70c92ab 10 * any party or parties not specified in the license agreement and/or
darienf 20:6d2af70c92ab 11 * nondisclosure agreement is expressly prohibited.
darienf 20:6d2af70c92ab 12 *
darienf 20:6d2af70c92ab 13 * The above copyright notice and this permission notice shall be included
darienf 20:6d2af70c92ab 14 * in all copies or substantial portions of the Software.
darienf 20:6d2af70c92ab 15 *
darienf 20:6d2af70c92ab 16 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
darienf 20:6d2af70c92ab 17 * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
darienf 20:6d2af70c92ab 18 * MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
darienf 20:6d2af70c92ab 19 * IN NO EVENT SHALL MAXIM INTEGRATED BE LIABLE FOR ANY CLAIM, DAMAGES
darienf 20:6d2af70c92ab 20 * OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE,
darienf 20:6d2af70c92ab 21 * ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR
darienf 20:6d2af70c92ab 22 * OTHER DEALINGS IN THE SOFTWARE.
darienf 20:6d2af70c92ab 23 *
darienf 20:6d2af70c92ab 24 * Except as contained in this notice, the name of Maxim Integrated
darienf 20:6d2af70c92ab 25 * Products, Inc. shall not be used except as stated in the Maxim Integrated
darienf 20:6d2af70c92ab 26 * Products, Inc. Branding Policy.
darienf 20:6d2af70c92ab 27 *
darienf 20:6d2af70c92ab 28 * The mere transfer of this software does not imply any licenses
darienf 20:6d2af70c92ab 29 * of trade secrets, proprietary technology, copyrights, patents,
darienf 20:6d2af70c92ab 30 * trademarks, maskwork rights, or any other form of intellectual
darienf 20:6d2af70c92ab 31 * property whatsoever. Maxim Integrated Products, Inc. retains all
darienf 20:6d2af70c92ab 32 * ownership rights.
darienf 20:6d2af70c92ab 33 *******************************************************************************
darienf 20:6d2af70c92ab 34 */
darienf 20:6d2af70c92ab 35
darienf 20:6d2af70c92ab 36 using System;
darienf 20:6d2af70c92ab 37 using System.Collections.Generic;
darienf 20:6d2af70c92ab 38 using System.Linq;
darienf 20:6d2af70c92ab 39 using System.Text;
darienf 20:6d2af70c92ab 40
darienf 20:6d2af70c92ab 41 namespace HealthSensorPlatform.Presenter
darienf 20:6d2af70c92ab 42 {
darienf 20:6d2af70c92ab 43 /// <summary>
darienf 20:6d2af70c92ab 44 /// R To R Calculator for referring R to R code values to ECG samples
darienf 20:6d2af70c92ab 45 /// </summary>
darienf 20:6d2af70c92ab 46 public class RToRCalculator
darienf 20:6d2af70c92ab 47 {
darienf 20:6d2af70c92ab 48 const int RToROffset = 3;
darienf 20:6d2af70c92ab 49
darienf 20:6d2af70c92ab 50 double masterClockFrequency;
darienf 20:6d2af70c92ab 51 int rToRDifferentialDelay;
darienf 20:6d2af70c92ab 52
darienf 20:6d2af70c92ab 53 int ecgDelay;
darienf 20:6d2af70c92ab 54 int rToRDelay;
darienf 20:6d2af70c92ab 55
darienf 20:6d2af70c92ab 56 double sampleRateRToR;
darienf 20:6d2af70c92ab 57 double sampleRateEcg;
darienf 20:6d2af70c92ab 58
darienf 20:6d2af70c92ab 59 double rToRTimeMilliSecond;
darienf 20:6d2af70c92ab 60
darienf 20:6d2af70c92ab 61 /// <summary>
darienf 20:6d2af70c92ab 62 /// Error difference between rounded R To R in ECG samples and actual calculated value. Error used
darienf 20:6d2af70c92ab 63 /// in calculation of next R To R in ECG samples point.
darienf 20:6d2af70c92ab 64 /// </summary>
darienf 20:6d2af70c92ab 65 double ecgSampleError;
darienf 20:6d2af70c92ab 66
darienf 20:6d2af70c92ab 67 /// <summary>
darienf 20:6d2af70c92ab 68 /// FMSTR clock frequency. Array index by FMSTR[1:0]
darienf 20:6d2af70c92ab 69 /// </summary>
darienf 20:6d2af70c92ab 70 double[] clockFreq = { 32768, 32000, 32000, 32768 * 640 / 656 }; // 31968.7804878 = 32768 * 640 / 65
darienf 20:6d2af70c92ab 71 /// <summary>
darienf 20:6d2af70c92ab 72 /// ECG decimation delay. Array indexed by FMSTR, ECG_RATE, ECG_DLFF (0 for disabled, 1 for enabled)
darienf 20:6d2af70c92ab 73 /// </summary>
darienf 20:6d2af70c92ab 74 int[, ,] ecgDecimationDelay = new int[,,] { {{650, 1034}, {2922, 3690}, {3370, 4906}, {3370, 4906}},
darienf 20:6d2af70c92ab 75 {{650, 1034}, {2922, 3690}, {3370, 4906}, {3370, 4906}},
darienf 20:6d2af70c92ab 76 {{1242, 2202}, {1242, 2202}, {1242, 2202}, {1242, 2202}},
darienf 20:6d2af70c92ab 77 {{1242, 2202}, {1242, 2202}, {1242, 2202}, {1242, 2202}}
darienf 20:6d2af70c92ab 78 };
darienf 20:6d2af70c92ab 79
darienf 20:6d2af70c92ab 80 /// <summary>
darienf 20:6d2af70c92ab 81 /// ECG sample rate look up based on FMSTR and ECG_RATE
darienf 20:6d2af70c92ab 82 /// </summary>
darienf 20:6d2af70c92ab 83 double[][] ecgSampleRates = new double[][]
darienf 20:6d2af70c92ab 84 {
darienf 20:6d2af70c92ab 85 new double[] {512, 256, 128, 128},
darienf 20:6d2af70c92ab 86 new double[] {500, 250, 125, 125},
darienf 20:6d2af70c92ab 87 new double[] {200, 200, 200, 200},
darienf 20:6d2af70c92ab 88 new double[] {199.8, 199.8, 199.8, 199.8}
darienf 20:6d2af70c92ab 89 };
darienf 20:6d2af70c92ab 90
darienf 20:6d2af70c92ab 91 /// <summary>
darienf 20:6d2af70c92ab 92 /// ECG processing delay in FMSTR clocks
darienf 20:6d2af70c92ab 93 /// </summary>
darienf 20:6d2af70c92ab 94 public int EcgDecimationDelay
darienf 20:6d2af70c92ab 95 {
darienf 20:6d2af70c92ab 96 get
darienf 20:6d2af70c92ab 97 {
darienf 20:6d2af70c92ab 98 return ecgDelay;
darienf 20:6d2af70c92ab 99 }
darienf 20:6d2af70c92ab 100 }
darienf 20:6d2af70c92ab 101
darienf 20:6d2af70c92ab 102 /// <summary>
darienf 20:6d2af70c92ab 103 /// R To R path processing delay in FMSTR clocks
darienf 20:6d2af70c92ab 104 /// </summary>
darienf 20:6d2af70c92ab 105 public int RToRDelay
darienf 20:6d2af70c92ab 106 {
darienf 20:6d2af70c92ab 107 get
darienf 20:6d2af70c92ab 108 {
darienf 20:6d2af70c92ab 109 return rToRDelay;
darienf 20:6d2af70c92ab 110 }
darienf 20:6d2af70c92ab 111 }
darienf 20:6d2af70c92ab 112
darienf 20:6d2af70c92ab 113 /// <summary>
darienf 20:6d2af70c92ab 114 /// Last R To R detected time in millisecond from first ECG point
darienf 20:6d2af70c92ab 115 /// </summary>
darienf 20:6d2af70c92ab 116 public double RToRMilliSecond
darienf 20:6d2af70c92ab 117 {
darienf 20:6d2af70c92ab 118 get
darienf 20:6d2af70c92ab 119 {
darienf 20:6d2af70c92ab 120 return rToRTimeMilliSecond;
darienf 20:6d2af70c92ab 121 }
darienf 20:6d2af70c92ab 122 }
darienf 20:6d2af70c92ab 123
darienf 20:6d2af70c92ab 124 /// <summary>
darienf 20:6d2af70c92ab 125 /// Constructor for R to R calculator
darienf 20:6d2af70c92ab 126 /// </summary>
darienf 20:6d2af70c92ab 127 /// <param name="masterClockFrequencyField">Value of FMSTR[1:0]</param>
darienf 20:6d2af70c92ab 128 /// <param name="sampleRateEcgField">Value of ECG_RATE[1:0]</param>
darienf 20:6d2af70c92ab 129 /// <param name="ecgDigitalLowPass">Value of ECG_DLF[1:0]</param>
darienf 20:6d2af70c92ab 130 /// <param name="rTorWindowField">Value of RTOR_WNDW[3:0]</param>
darienf 20:6d2af70c92ab 131 public RToRCalculator(int masterClockFrequencyField, int sampleRateEcgField, int ecgDigitalLowPass, int rTorWindowField)
darienf 20:6d2af70c92ab 132 {
darienf 20:6d2af70c92ab 133 masterClockFrequency = clockFreq[masterClockFrequencyField];
darienf 20:6d2af70c92ab 134 sampleRateRToR = masterClockFrequency / 256;
darienf 20:6d2af70c92ab 135 sampleRateEcg = ecgSampleRates[masterClockFrequencyField][sampleRateEcgField];
darienf 20:6d2af70c92ab 136
darienf 20:6d2af70c92ab 137 ecgDelay = ecgDecimationDelay[masterClockFrequencyField, sampleRateEcgField, ecgDigitalLowPass > 0 ? 1 : 0];
darienf 20:6d2af70c92ab 138 rToRDelay = 5376 + 3370 + 256 * rTorWindowField;
darienf 20:6d2af70c92ab 139
darienf 20:6d2af70c92ab 140 rToRDifferentialDelay = rToRDelay - ecgDelay;
darienf 20:6d2af70c92ab 141
darienf 20:6d2af70c92ab 142 rToRTimeMilliSecond = 0;
darienf 20:6d2af70c92ab 143 }
darienf 20:6d2af70c92ab 144
darienf 20:6d2af70c92ab 145 /// <summary>
darienf 20:6d2af70c92ab 146 /// R To R code to beats per minute
darienf 20:6d2af70c92ab 147 /// </summary>
darienf 20:6d2af70c92ab 148 /// <param name="rToRCode">R To R Code</param>
darienf 20:6d2af70c92ab 149 /// <returns></returns>
darienf 20:6d2af70c92ab 150 public double BeatsPerMinute(int rToRCode)
darienf 20:6d2af70c92ab 151 {
darienf 20:6d2af70c92ab 152 return 60 * 1000 / Millisecond(rToRCode);
darienf 20:6d2af70c92ab 153 }
darienf 20:6d2af70c92ab 154
darienf 20:6d2af70c92ab 155 /// <summary>
darienf 20:6d2af70c92ab 156 /// R To R code to millisecond interval
darienf 20:6d2af70c92ab 157 /// </summary>
darienf 20:6d2af70c92ab 158 /// <param name="rToRCode">R To R Code</param>
darienf 20:6d2af70c92ab 159 /// <returns></returns>
darienf 20:6d2af70c92ab 160 public double Millisecond(double rToRCode)
darienf 20:6d2af70c92ab 161 {
darienf 20:6d2af70c92ab 162 return rToRCode * 512.0 * 1000 / (2 * masterClockFrequency);
darienf 20:6d2af70c92ab 163 }
darienf 20:6d2af70c92ab 164
darienf 20:6d2af70c92ab 165 public double Corrected(int rToR, bool first)
darienf 20:6d2af70c92ab 166 {
darienf 20:6d2af70c92ab 167 if (first)
darienf 20:6d2af70c92ab 168 rToRTimeMilliSecond = Millisecond(rToR - rToRDifferentialDelay / 256.0);
darienf 20:6d2af70c92ab 169 else
darienf 20:6d2af70c92ab 170 rToRTimeMilliSecond += Millisecond(rToR);
darienf 20:6d2af70c92ab 171
darienf 20:6d2af70c92ab 172 return rToRTimeMilliSecond;
darienf 20:6d2af70c92ab 173 }
darienf 20:6d2af70c92ab 174
darienf 20:6d2af70c92ab 175 /// <summary>
darienf 20:6d2af70c92ab 176 /// R To R code converted to ECG points
darienf 20:6d2af70c92ab 177 /// </summary>
darienf 20:6d2af70c92ab 178 /// <param name="rToR">R To R code</param>
darienf 20:6d2af70c92ab 179 /// <param name="first">first R To R sample include the processing delay between the ECG and R To R processing path</param>
darienf 20:6d2af70c92ab 180 /// <returns></returns>
darienf 20:6d2af70c92ab 181 public int EcgPoints(int rToR, bool first)
darienf 20:6d2af70c92ab 182 {
darienf 20:6d2af70c92ab 183 double rToRInEcgSamples;
darienf 20:6d2af70c92ab 184 int rToRInEcgSamplesInt;
darienf 20:6d2af70c92ab 185
darienf 20:6d2af70c92ab 186 if (first)
darienf 20:6d2af70c92ab 187 rToRInEcgSamples = (rToR - rToRDifferentialDelay / 256.0) * (sampleRateEcg / sampleRateRToR) + RToROffset;
darienf 20:6d2af70c92ab 188 else
darienf 20:6d2af70c92ab 189 rToRInEcgSamples = rToR * (sampleRateEcg / sampleRateRToR) - ecgSampleError;
darienf 20:6d2af70c92ab 190
darienf 20:6d2af70c92ab 191 rToRInEcgSamplesInt = (int)(rToRInEcgSamples + 0.5);
darienf 20:6d2af70c92ab 192 ecgSampleError = rToRInEcgSamplesInt - rToRInEcgSamples;
darienf 20:6d2af70c92ab 193
darienf 20:6d2af70c92ab 194 return rToRInEcgSamplesInt;
darienf 20:6d2af70c92ab 195 }
darienf 20:6d2af70c92ab 196 }
darienf 20:6d2af70c92ab 197 }