How can i visualize ecg contain 12 leads in 12 different figures

I have a .dat file and a .hea file which contains all 12 ecg leads.

i want to visualize each ecg lead in a different figure

I use the following code to store .dat data to a variable

fp= fopen('00001_lr.dat','rb');

ecg=fread(fp,'short');

frequency is 100hz and length is 10 secs and contain 12 leads (V1 to V6, I, II, III, aVR, aVL, and aVF )

thanks and regards.

[Ameer Hamza]

How the data is formatted? Can you share a data file?

[Naveen F]

Thank you very much Ameer Hamza for comment on my question,

I try to attach the .dat and .hea files to this comment but matlab says its unsupported file format so here i meantion the link where i get the data.

This is the link for the database:

https://www.physionet.org/content/ptb-xl/1.0.1/

Here i'll a screen shot of the ecg i take by visiting the physionet lightwave website

You can visualize these data by visiting this link

https://www.physionet.org/lightwave/?db=ptb-xl/1.0.1

[Naveen F]

This the text insie the .hea file

00001_lr 12 100 1000

00001_lr.dat 16 1000.0(0)/mV 16 0 -119 1508 0 I

00001_lr.dat 16 1000.0(0)/mV 16 0 -55 723 0 II

00001_lr.dat 16 1000.0(0)/mV 16 0 64 64758 0 III

00001_lr.dat 16 1000.0(0)/mV 16 0 86 64423 0 AVR

00001_lr.dat 16 1000.0(0)/mV 16 0 -91 1211 0 AVL

00001_lr.dat 16 1000.0(0)/mV 16 0 4 7 0 AVF

00001_lr.dat 16 1000.0(0)/mV 16 0 -69 63827 0 V1

00001_lr.dat 16 1000.0(0)/mV 16 0 -31 6999 0 V2

00001_lr.dat 16 1000.0(0)/mV 16 0 0 63759 0 V3

00001_lr.dat 16 1000.0(0)/mV 16 0 -26 61447 0 V4

00001_lr.dat 16 1000.0(0)/mV 16 0 -39 64979 0 V5

00001_lr.dat 16 1000.0(0)/mV 16 0 -79 832 0 V6

[FSK]

You need to download and configure the wfdb app for matlab:

https://archive.physionet.org/physiotools/matlab/wfdb-app-matlab/

and then run the command :

[signal, fs, tm] = rdsamp('test01_00s');

where test01_00s is the name of the .dat file you want to read. You can then plot and work with variable signal.

Hope it would help.

[Naveen F]

Hi FSK, Ameer Hamza,

I download and configure the wfdb app for matlab.

And I run the following code,

the sample code provided in the blog works like a charm

[sig, Fs, tm] = rdsamp('mitdb/100', 1);

plot(tm, sig);

---------------------------------------------------------------------

and then i run this code

with the extension

[sig, Fs, tm] = rdsamp('01000_lr.dat', 1);

plot(tm, sig);

without the extension

[sig, Fs, tm] = rdsamp('01000_lr', 1);

plot(tm, sig);

the database path from google cloud

[sig, Fs, tm] = rdsamp('ptb-xl/records100/00000/01000_lr');

plot(tm, sig);

after all these attempts,

I check it in the PhysioBank ATM database whether the ptb-xl is exsist. And i found its not in their database.

But they have publish ptb-xl dataset on April 24, 2020 in their https://www.physionet.org website

This is the link about the publication : https://www.physionet.org/content/ptb-xl/1.0.1/

However data can visualize in lightwave site : https://www.physionet.org/lightwave/?db=ptb-xl/1.0.1

so i my old code and whether the ecg contains all 12 leads using the following code,

fp= fopen('00001_lr.dat','rb');

ecg=fread(fp,'short');

It shows data only for a one lead since the lightwave site shows it contains 12 leads..

And here i attach a zip file which include both .dat file and .hea file i use.

I am so confused in this case please please help me.

thanks and best regards.

[FSK]

I have opened the field you attached using the follwoing command :

[signal, fs, tm] = rdsamp('00001_lr');

And it shows signal has 1000*12 values. So if I understand the question correctly, it does have 12 lead values.

Just make sure your matlab directory is also set to the location where your file .dat and .hea is located.

I am attaching the signal screenshot and all the signals retrieved for reference.

[Kevin T]

Hi, I have the same error when trying to read a local .dat file.

This is my matlab directory, my .dat file is located in F:\Open ECG file data.

Could you please help me solve this problem?

[Star Strider]

Kevin T —

I auggest that you post this as a new Question, and include a lot more detail, specifically uploading the ‘.dat’ file as an attachment.

Answer 1

[Star Strider]

This is the standard format for EKG record display:

EKGD = load('temp.mat');
EKG = EKGD.signal;

Leadc = {'I','II','III','aV_R','aV_L','aV_F','V_1','V_2','V_3','V_4','V_5','V_6'};
figure
kc = 1;
for k1 = 1:4
    for k2 = 1:4:9
        sbpt = (k1-1)+k2;
        subplot(3,4,sbpt)
        plot(EKG(1:300,kc))
        title(Leadc{kc})
        grid
        ylim([-1 1])
        kc = kc+1;
    end
end

producing:

Characteristically, only three cycles are displayed, except for the separate Lead II rhythm strip (not shown here, would be a separate plot).

Diagnosis: Abnormal EKG, with pronounced left-axis deviation,with an axis of about -45°. Impression: Consistent with left bundle-branch block.

[FSK]

Star! can you recommend some good source to understand the signal acquisition and analysis for EKG?

[Star Strider]

There are so many possibilities with respect to acquisition that I cannot recommend any. That has to depend on your particular application, computer infrastructure (including hardware and OS) and related resources. I recommend that you go online and see what best fits your needs, making certain that it can produce the sort of outputs that MATLAB requires (for example text files or specific interfaces). Choosing instrumentation is not a trivial problem, and will take some time to get the one that works best for you. MathWorks may also have resources with respect to instrumentation they have tested. I have no idea, because I have never looked for it.

With respect to analysis, I have written several Answers in thes area. There are likely other applications available online that might be able to do what you want, and that are rigorously tested and supported.

[Naveen F]

Thank you very much Star Strider, your code work perfectly .. FSK thank you too for helping me configure the wfdb. you guys are great. keep up the good work :)

[Star Strider]

Naveen Fernando —

My pleasure!

If my Answer helped you solve your problem, please Accept it!

.