DATA PROCESSING FOR SPO2 AND HR
Dependents: MAX30100_FirstTry MAX30100_V04
PROCESAMIENTO_SPO2_FC.cpp@0:264c403954d8, 2017-03-27 (annotated)
- Committer:
- cesarcazal
- Date:
- Mon Mar 27 09:02:23 2017 +0000
- Revision:
- 0:264c403954d8
- Child:
- 1:7e5b24227d0a
- Child:
- 2:2cd93c474a32
DATA PROCESSING FOR SPO2 AND HR
Who changed what in which revision?
User | Revision | Line number | New contents of line |
---|---|---|---|
cesarcazal | 0:264c403954d8 | 1 | /*int muestras_50hz[200]; |
cesarcazal | 0:264c403954d8 | 2 | int matriz_minimos[30][2], matriz_maximos[30][2],max_min_period[2]; |
cesarcazal | 0:264c403954d8 | 3 | int pico_pico,cresta, valle, cantidad; |
cesarcazal | 0:264c403954d8 | 4 | float frecuencia,periodo;*/ |
cesarcazal | 0:264c403954d8 | 5 | |
cesarcazal | 0:264c403954d8 | 6 | //int index, aux1, nmuestra; |
cesarcazal | 0:264c403954d8 | 7 | |
cesarcazal | 0:264c403954d8 | 8 | /*EJEMPLO DE UTILIZACION DEL SISTEMA DE PROCESAMIENTO. |
cesarcazal | 0:264c403954d8 | 9 | int main(int argc, char *argv[]) { |
cesarcazal | 0:264c403954d8 | 10 | promedio(muestras_100hz, muestras_50hz); |
cesarcazal | 0:264c403954d8 | 11 | periodo = valores(muestras_50hz,matriz_minimos,matriz_maximos,max_min_period); |
cesarcazal | 0:264c403954d8 | 12 | for(aux1=0;aux1<30;aux1++){ |
cesarcazal | 0:264c403954d8 | 13 | printf("%d,%d \n%d,%d \n",matriz_maximos[aux1][0],matriz_maximos[aux1][1],matriz_minimos[aux1][0],matriz_minimos[aux1][1]); |
cesarcazal | 0:264c403954d8 | 14 | //printf("%d , %d \n",matriz_minimos[aux1][0],matriz_minimos[aux1][1]); |
cesarcazal | 0:264c403954d8 | 15 | }*/ |
cesarcazal | 0:264c403954d8 | 16 | /* pico_pico=max_min_period[0]-max_min_period[1]; |
cesarcazal | 0:264c403954d8 | 17 | frecuencia=60/(periodo*0.02); |
cesarcazal | 0:264c403954d8 | 18 | //printf("el valor medio del maximo es %d\n", pico_pico); |
cesarcazal | 0:264c403954d8 | 19 | printf("el valor medio del pico_pico es %f\n", frecuencia); |
cesarcazal | 0:264c403954d8 | 20 | return 0; |
cesarcazal | 0:264c403954d8 | 21 | }*/ |
cesarcazal | 0:264c403954d8 | 22 | |
cesarcazal | 0:264c403954d8 | 23 | float valores(int muestras[200],int min_pos[30][2], int max_pos[30][2],int resultados[2]){//EMITE EL PERIODO, EL VALOR MAXIMO Y EL VALOR MINIMO COMO RESPUESTA. |
cesarcazal | 0:264c403954d8 | 24 | int a,b=0,posibilidad_anteriores=0,posibilidad_sucesores=0,indice=0; |
cesarcazal | 0:264c403954d8 | 25 | int posicion=0,Vpp,minimos=4194303, prom_max=0, prom_min=0, periodo_max=0, periodo_min=0,maximos=-4194303, magnitud=4194303, copy_min[30][2],copy_max[30][2]; |
cesarcazal | 0:264c403954d8 | 26 | float res_period,num_prom; |
cesarcazal | 0:264c403954d8 | 27 | for(a=6;a<194;a++){ |
cesarcazal | 0:264c403954d8 | 28 | if((muestras[a]<magnitud)&&(muestras[a]>15000)){//segun recomendacion de https://www.youtube.com/watch?v=36cELpNGJZY |
cesarcazal | 0:264c403954d8 | 29 | magnitud=muestras[a]; |
cesarcazal | 0:264c403954d8 | 30 | posicion=a; |
cesarcazal | 0:264c403954d8 | 31 | for (b=1;b<7;b++){ |
cesarcazal | 0:264c403954d8 | 32 | if((muestras[a]<=muestras[a-b])||(muestras[a-b]<=15000)){//segun recomendacion de https://www.youtube.com/watch?v=36cELpNGJZY |
cesarcazal | 0:264c403954d8 | 33 | posibilidad_anteriores++; |
cesarcazal | 0:264c403954d8 | 34 | } |
cesarcazal | 0:264c403954d8 | 35 | if((muestras[a]<=muestras[a+b])||(muestras[a-b]<=15000)){//segun recomendacion de https://www.youtube.com/watch?v=36cELpNGJZY |
cesarcazal | 0:264c403954d8 | 36 | posibilidad_sucesores++; |
cesarcazal | 0:264c403954d8 | 37 | } |
cesarcazal | 0:264c403954d8 | 38 | } |
cesarcazal | 0:264c403954d8 | 39 | if ((posibilidad_anteriores>5)&&(posibilidad_sucesores>5)){ |
cesarcazal | 0:264c403954d8 | 40 | min_pos[indice][0]=magnitud; |
cesarcazal | 0:264c403954d8 | 41 | min_pos[indice][1]=posicion; |
cesarcazal | 0:264c403954d8 | 42 | if(magnitud<minimos) |
cesarcazal | 0:264c403954d8 | 43 | minimos=magnitud; |
cesarcazal | 0:264c403954d8 | 44 | indice++; |
cesarcazal | 0:264c403954d8 | 45 | a=a+10; |
cesarcazal | 0:264c403954d8 | 46 | } |
cesarcazal | 0:264c403954d8 | 47 | magnitud=4194303; |
cesarcazal | 0:264c403954d8 | 48 | posibilidad_anteriores=0; |
cesarcazal | 0:264c403954d8 | 49 | posibilidad_sucesores=0; |
cesarcazal | 0:264c403954d8 | 50 | } |
cesarcazal | 0:264c403954d8 | 51 | } |
cesarcazal | 0:264c403954d8 | 52 | indice=0; |
cesarcazal | 0:264c403954d8 | 53 | magnitud=-4194303; |
cesarcazal | 0:264c403954d8 | 54 | posicion=0; |
cesarcazal | 0:264c403954d8 | 55 | for (a=6;a<194;a++){ |
cesarcazal | 0:264c403954d8 | 56 | if((muestras[a]>magnitud)&&(muestras[a]<125000)){//EL 125000 DEBE SER 45000 segun recomendacion de https://www.youtube.com/watch?v=36cELpNGJZY |
cesarcazal | 0:264c403954d8 | 57 | magnitud=muestras[a]; |
cesarcazal | 0:264c403954d8 | 58 | posicion=a; |
cesarcazal | 0:264c403954d8 | 59 | for (b=1;b<7;b++){ |
cesarcazal | 0:264c403954d8 | 60 | if((muestras[a]>=muestras[a-b])||(muestras[a-b]>=125000)){//EL 125000 DEBE SER 45000 segun recomendacion de https://www.youtube.com/watch?v=36cELpNGJZY |
cesarcazal | 0:264c403954d8 | 61 | posibilidad_anteriores++; |
cesarcazal | 0:264c403954d8 | 62 | } |
cesarcazal | 0:264c403954d8 | 63 | if((muestras[a]>=muestras[a+b])||(muestras[a-b]>=125000)){//EL 125000 DEBE SER 45000 segun recomendacion de https://www.youtube.com/watch?v=36cELpNGJZY |
cesarcazal | 0:264c403954d8 | 64 | posibilidad_sucesores++; |
cesarcazal | 0:264c403954d8 | 65 | } |
cesarcazal | 0:264c403954d8 | 66 | } |
cesarcazal | 0:264c403954d8 | 67 | if ((posibilidad_anteriores>5)&&(posibilidad_sucesores>5)){ |
cesarcazal | 0:264c403954d8 | 68 | max_pos[indice][0]=magnitud; |
cesarcazal | 0:264c403954d8 | 69 | max_pos[indice][1]=posicion; |
cesarcazal | 0:264c403954d8 | 70 | if(magnitud>maximos) |
cesarcazal | 0:264c403954d8 | 71 | maximos=magnitud; |
cesarcazal | 0:264c403954d8 | 72 | indice++; |
cesarcazal | 0:264c403954d8 | 73 | a=a+10; |
cesarcazal | 0:264c403954d8 | 74 | } |
cesarcazal | 0:264c403954d8 | 75 | posibilidad_anteriores=0; |
cesarcazal | 0:264c403954d8 | 76 | posibilidad_sucesores=0; |
cesarcazal | 0:264c403954d8 | 77 | magnitud=-4194303; |
cesarcazal | 0:264c403954d8 | 78 | } |
cesarcazal | 0:264c403954d8 | 79 | } |
cesarcazal | 0:264c403954d8 | 80 | a=0; |
cesarcazal | 0:264c403954d8 | 81 | b=0; |
cesarcazal | 0:264c403954d8 | 82 | indice=0; |
cesarcazal | 0:264c403954d8 | 83 | Vpp=maximos-minimos; |
cesarcazal | 0:264c403954d8 | 84 | if (max_pos[0][1]<min_pos[0][1]){ |
cesarcazal | 0:264c403954d8 | 85 | if(((max_pos[0][0]-min_pos[0][0])*2)<Vpp){ |
cesarcazal | 0:264c403954d8 | 86 | max_pos[0][1]=401; |
cesarcazal | 0:264c403954d8 | 87 | max_pos[0][0]=maximos; |
cesarcazal | 0:264c403954d8 | 88 | min_pos[0][1]=401; |
cesarcazal | 0:264c403954d8 | 89 | min_pos[0][0]=minimos; |
cesarcazal | 0:264c403954d8 | 90 | |
cesarcazal | 0:264c403954d8 | 91 | } |
cesarcazal | 0:264c403954d8 | 92 | a=1; |
cesarcazal | 0:264c403954d8 | 93 | } |
cesarcazal | 0:264c403954d8 | 94 | while ((max_pos[a][1])!=0){ |
cesarcazal | 0:264c403954d8 | 95 | if((((max_pos[a][0]-min_pos[b][0])*2)<Vpp)&&((max_pos[a][0]-min_pos[b][0])<(max_pos[a][0]-min_pos[b+1][0]))){ |
cesarcazal | 0:264c403954d8 | 96 | max_pos[a][1]=401; |
cesarcazal | 0:264c403954d8 | 97 | min_pos[b][1]=401; |
cesarcazal | 0:264c403954d8 | 98 | max_pos[a][0]=maximos; |
cesarcazal | 0:264c403954d8 | 99 | min_pos[b][0]=minimos; |
cesarcazal | 0:264c403954d8 | 100 | } |
cesarcazal | 0:264c403954d8 | 101 | else{ |
cesarcazal | 0:264c403954d8 | 102 | if(((max_pos[a][0]-min_pos[b+1][0])*2)<Vpp){ |
cesarcazal | 0:264c403954d8 | 103 | max_pos[a][1]=401; |
cesarcazal | 0:264c403954d8 | 104 | min_pos[b+1][1]=401; |
cesarcazal | 0:264c403954d8 | 105 | max_pos[a][0]=maximos; |
cesarcazal | 0:264c403954d8 | 106 | min_pos[b+1][0]=minimos; |
cesarcazal | 0:264c403954d8 | 107 | } |
cesarcazal | 0:264c403954d8 | 108 | } |
cesarcazal | 0:264c403954d8 | 109 | b++; |
cesarcazal | 0:264c403954d8 | 110 | a++; |
cesarcazal | 0:264c403954d8 | 111 | } |
cesarcazal | 0:264c403954d8 | 112 | b=0; |
cesarcazal | 0:264c403954d8 | 113 | for (a=0;a<30;a++){ |
cesarcazal | 0:264c403954d8 | 114 | if ((max_pos[a][1]>0)&&(max_pos[a][1]<400)){ |
cesarcazal | 0:264c403954d8 | 115 | copy_max[indice][0]=max_pos[a][0]; |
cesarcazal | 0:264c403954d8 | 116 | copy_max[indice][1]=max_pos[a][1]; |
cesarcazal | 0:264c403954d8 | 117 | prom_max=prom_max+max_pos[a][0]; |
cesarcazal | 0:264c403954d8 | 118 | if (b>0){ |
cesarcazal | 0:264c403954d8 | 119 | periodo_max=periodo_max+max_pos[a][1]-copy_max[indice-1][1]; |
cesarcazal | 0:264c403954d8 | 120 | } |
cesarcazal | 0:264c403954d8 | 121 | indice++; |
cesarcazal | 0:264c403954d8 | 122 | } |
cesarcazal | 0:264c403954d8 | 123 | if ((min_pos[a][1]>0)&&(min_pos[a][1]<400)){ |
cesarcazal | 0:264c403954d8 | 124 | copy_min[b][0]=min_pos[a][0]; |
cesarcazal | 0:264c403954d8 | 125 | copy_min[b][1]=min_pos[a][1]; |
cesarcazal | 0:264c403954d8 | 126 | prom_min=prom_min+min_pos[a][0]; |
cesarcazal | 0:264c403954d8 | 127 | if (b>0){ |
cesarcazal | 0:264c403954d8 | 128 | periodo_min=periodo_min+min_pos[a][1]-copy_min[b-1][1]; |
cesarcazal | 0:264c403954d8 | 129 | } |
cesarcazal | 0:264c403954d8 | 130 | b++; |
cesarcazal | 0:264c403954d8 | 131 | } |
cesarcazal | 0:264c403954d8 | 132 | if((max_pos[a+1][1]==0)&&(min_pos[a+1][1]==0)){ |
cesarcazal | 0:264c403954d8 | 133 | posicion=a+1; |
cesarcazal | 0:264c403954d8 | 134 | a=29; |
cesarcazal | 0:264c403954d8 | 135 | } |
cesarcazal | 0:264c403954d8 | 136 | } |
cesarcazal | 0:264c403954d8 | 137 | resultados[0]=prom_max/indice;//VALOR MAXIMO |
cesarcazal | 0:264c403954d8 | 138 | resultados[1]=prom_min/b;//VALOR MINIMO |
cesarcazal | 0:264c403954d8 | 139 | res_period=(periodo_max+periodo_min); |
cesarcazal | 0:264c403954d8 | 140 | num_prom=(b+indice-2); |
cesarcazal | 0:264c403954d8 | 141 | res_period=res_period/num_prom; |
cesarcazal | 0:264c403954d8 | 142 | //printf("minimo %d \n", period); |
cesarcazal | 0:264c403954d8 | 143 | if (indice>b) |
cesarcazal | 0:264c403954d8 | 144 | b=indice; |
cesarcazal | 0:264c403954d8 | 145 | for(a=0;a<30;a++){ |
cesarcazal | 0:264c403954d8 | 146 | if (a<b){ |
cesarcazal | 0:264c403954d8 | 147 | max_pos[a][0]=copy_max[a][0]; |
cesarcazal | 0:264c403954d8 | 148 | max_pos[a][1]=copy_max[a][1]; |
cesarcazal | 0:264c403954d8 | 149 | min_pos[a][0]=copy_min[a][0]; |
cesarcazal | 0:264c403954d8 | 150 | min_pos[a][1]=copy_min[a][1]; |
cesarcazal | 0:264c403954d8 | 151 | } |
cesarcazal | 0:264c403954d8 | 152 | else{ |
cesarcazal | 0:264c403954d8 | 153 | max_pos[a][0]=0; |
cesarcazal | 0:264c403954d8 | 154 | max_pos[a][1]=0; |
cesarcazal | 0:264c403954d8 | 155 | min_pos[a][0]=0; |
cesarcazal | 0:264c403954d8 | 156 | min_pos[a][1]=0; |
cesarcazal | 0:264c403954d8 | 157 | } |
cesarcazal | 0:264c403954d8 | 158 | } |
cesarcazal | 0:264c403954d8 | 159 | return(res_period); |
cesarcazal | 0:264c403954d8 | 160 | } |
cesarcazal | 0:264c403954d8 | 161 | |
cesarcazal | 0:264c403954d8 | 162 | void promedio(int *muestras400, int *muestras200){ |
cesarcazal | 0:264c403954d8 | 163 | int aux; |
cesarcazal | 0:264c403954d8 | 164 | for (aux=0; aux<400; aux=aux+2){ |
cesarcazal | 0:264c403954d8 | 165 | *(muestras200+aux/2)=(*(muestras400+aux)+*(muestras400+aux+1))/2; |
cesarcazal | 0:264c403954d8 | 166 | } |
cesarcazal | 0:264c403954d8 | 167 | } |