Question

我从file.now读取ecg字节数组我想检测读取ecg字节的QRS。我怎样才能在java中实现这一点。

我从Lifegain defibrilator（一个ecg设备）获得字节数组。我从这些字节中绘制了关于android的ecg。现在我想检测QRS复合体（用于计算一个心跳波的时间和电压的术语）。 DATA = LeadData “：” - 284，-127，-122,17，-35,10,32,10,52,16,49,33,38,69,70,58,45,93,47,88， 58,90,149,5,82，-12，-4,40，-34,29，-29,5，-4，-17，-13，-29，-13，-4，-9，-9， -10，-20，-15，-22，-32，-25，-23，-2，-15，-7，-13，-19，-17，-28，-27，-27，-33 ，-20，-16，-13，-20，-10，-22，-20，-19，-28，-15，-19，-22，-21，-9，-3，-6， - 8，-6，-11，-8，-8，-5，-10，-5，-6，-9，-4，-6,3,20,3,14,7,11,10,5 ，11,5,10,2,10,13,14"

此致

沙

Answer 1

如果您拥有的数据是我认为您拥有的数据，则需要使用其中一种算法来检测您的QRS波群。

有很多算法用于检测QRS波群，其中最简单的算法之一是基于移动平均线的滤波系统及其在实时QRS检测中的应用由HC Chen和SW Chen提供（你可以在http://www.cinc.org/archives/2003/pdf/585.pdf上获得。）

阶段是：

高通过滤
低通过滤
描述制作阶段

从低通图片中你可以注意到，现在我们有了检测QRS波群所需的峰值。最后一个阶段是制定决策阶段。在本文中，您有一个实现此目的的公式。

我们需要知道QRS复合体何时开始，所以我们需要为此设置一个阈值。此实现中的公式为：

阈值= alpha * gamma * peak +（1 - alpha）* threshold

峰值是窗口中的局部最大值（我们通常用宽度为250的窗口搜索信号），阈值是初始值（第一个可以是第一个找到的峰值），alpha和gamma是随机创建的，α大于0且小于1，而γ为0.15或0.20。如果当前信号的值超过阈值，则找到QRS复合波。

以下是Java中用于低通，高通和决策的源代码：

// High pass filter
// y1[n] = 1/M * Sum[m=0, M-1] x[n-m]
// y2[n] = x[n - (M+1)/2]
public static float[] highPass(int[] sig0, int nsamp) {
    float[] highPass = new float[nsamp];
    int M = 5; // M is recommended to be 5 or 7 according to the paper
    float constant = (float) 1/M;

    for(int i=0; i<sig0.length; i++) {
        float y1 = 0;
        float y2 = 0;

        int y2_index = i-((M+1)/2);
        if(y2_index < 0) {
            y2_index = nsamp + y2_index;
        }
        y2 = sig0[y2_index];

        float y1_sum = 0;
        for(int j=i; j>i-M; j--) {
            int x_index = i - (i-j);
            if(x_index < 0) {
                x_index = nsamp + x_index;
            }
            y1_sum += sig0[x_index];
        }

        y1 = constant * y1_sum;
        highPass[i] = y2 - y1;

    }        

    return highPass;
}

// Low pass filter; na n-to mesto zapiši kvadrat 30ih števil v oknu
public static float[] lowPass(float[] sig0, int nsamp) {
    float[] lowPass = new float[nsamp];
    for(int i=0; i<sig0.length; i++) {
        float sum = 0;
        if(i+30 < sig0.length) {
            for(int j=i; j<i+30; j++) {
                float current = sig0[j] * sig0[j];
                sum += current;
            }
        }
        else if(i+30 >= sig0.length) {
            int over = i+30 - sig0.length;
            for(int j=i; j<sig0.length; j++) {
                float current = sig0[j] * sig0[j];
                sum += current;
            }
            for(int j=0; j<over; j++) {
                float current = sig0[j] * sig0[j];
                sum += current;
            }
        }

        lowPass[i] = sum;
    }

    return lowPass;

}

public static int[] QRS(float[] lowPass, int nsamp) {
    int[] QRS = new int[nsamp];

    double treshold = 0;

    for(int i=0; i<200; i++) {
        if(lowPass[i] > treshold) {
            treshold = lowPass[i];
        }
    }

    int frame = 250;

    for(int i=0; i<lowPass.length; i+=frame) {
        float max = 0;
        int index = 0;
        if(i + frame > lowPass.length) {
            index = lowPass.length;
        }
        else {
            index = i + frame;
        }
        for(int j=i; j<index; j++) {
            if(lowPass[j] > max) max = lowPass[j];
        }
        boolean added = false;
        for(int j=i; j<index; j++) {
            if(lowPass[j] > treshold && !added) {
                QRS[j] = 1;
                added = true;
            }
            else {
                QRS[j] = 0;
            }                
        }

        double gama = (Math.random() > 0.5) ? 0.15 : 0.20;
        double alpha = 0.01 + (Math.random() * ((0.1 - 0.01)));

        treshold = alpha * gama * max + (1 - alpha) * treshold;

    }

    return QRS;
}

Answer 2

请按照以下链接，我认为他们会帮助您。

http://www.cinc.org/archives/2008/pdf/0441.pdf

http://carg.site.uottawa.ca/doc/ELG6163GeoffreyGreen.pdf

http://www.eplimited.com/osea13.pdf

http://mirel.xmu.edu.cn/mirel/public/Teaching/QRSdetection.pdf

http://sourceforge.net/projects/ecgtoolkit-cs/files/ecgtoolkit-cs/ecgtoolkit-cs-2_2/

从ecg字节数组中检索java中的QRS

2 个答案: