Hi,
I am trying to send ECG datd to ThingSpeak, however, at one sample per 15seconds it is rather useless when I have 650k samples.
What I want to do is send a certain amount of samples before and after each R peak (I am usung find peaks to locate R-peaks after filtering the data). That way I can reduce the amount of samples sent and get more ECG signals inside the graph.
Can anyone advise?
PS - Patient 101 is data from MIT database.
Thank you in advance!
Gabriel N. Nicolau
My current code is below:
% clc
% clear all
% close all
%load ( 'patient100.csv');
dat=csvread('patient100.csv');
figure;
subplot(4,1,1)
plot(dat)
findpeaks(dat,'MinPeakHeight',0.5,'MinPeakDistance',200);
title('Original')
x=lowpass(dat,15,360);
y=highpass(x,10,360);
subplot(4,1,2)
plot(x)
findpeaks(x,'MinPeakHeight',0.5,'MinPeakDistance',200);
title('Low-Pass')
subplot(4,1,3)
plot(y)
findpeaks(y,'MinPeakHeight',0.3,'MinPeakDistance',200);
title('High-Pass')
y_diff=diff(y);
y_abs = y_diff.^2;
subplot(4,1,4)
plot(y_abs);
%for i=1: 1/length(dat)*100000: length(dat)*10000+ 1/length(dat)*100000;
%t=dat(dat(i)>dat(i+1));
%end
yy = y';
%%% DATA TO BE SENT TO THINGSPEAK CHANNEL
a = yy;
b = 0;
c = 650000;
x = [a];
chId = 622307;
writeKey = 'XXXYYYZZZ';
tOut = 99999999;
% tStamps = [datetime('now')-minutes(899):minutes(1):datetime('now')]';
for i = 1:c
pause(15)
b = x(i)
thingSpeakWrite(chId,[b],'WriteKey',writeKey)
end

1 commentaire

ekta yadav
ekta yadav le 20 Nov 2021
Hi I am working on the same thing but stuck on few things.
  1. Can you tell me how did you calculated number of samples of MIT database.
  2. How u sent numeric values of heart beat.
Thanks in advance.

Connectez-vous pour commenter.

 Réponse acceptée

GABRIEL NICOLAU
GABRIEL NICOLAU le 17 Nov 2018
Modifié(e) : GABRIEL NICOLAU le 18 Nov 2018

0 votes

I have resolved the issue.
Thank you.

2 commentaires

Ali Usman
Ali Usman le 15 Jan 2019
Can you please share the details that how you resolved the issue..?
GABRIEL NICOLAU
GABRIEL NICOLAU le 15 Jan 2019
yes. unless you pay to have an account you will not be able to send the data faster than at 15 sec intervals
i resolved the issue by calculating the current heart rate and sending the numerical value (ie 60BPM) every 15 seconds
this was done by analysing the R peaks at 15second intervals
gabriel.

Connectez-vous pour commenter.

Plus de réponses (0)

Communautés

Plus de réponses dans  ThingSpeak Community

Catégories

En savoir plus sur ThingSpeak dans Centre d'aide et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by