示例#1
0
        /// <summary>
        /// Function to resample all leads.
        /// </summary>
        /// <param name="samplesPerSecond">samples per second to resample towards</param>
        public void Resample(int samplesPerSecond)
        {
            foreach (Signal sig in this._Lead)
            {
                if ((this.RhythmSamplesPerSecond != 0) &&
                    (this.RhythmAVM != 0) &&
                    (sig.Rhythm != null))
                {
                    ECGTool.ResampleLead(sig.Rhythm, this.RhythmSamplesPerSecond, samplesPerSecond, out sig.Rhythm);

                    sig.RhythmStart = (int)(((long)sig.RhythmStart * (long)samplesPerSecond) / (long)this.RhythmSamplesPerSecond);
                    sig.RhythmEnd   = (int)(((long)sig.RhythmEnd * (long)samplesPerSecond) / (long)this.RhythmSamplesPerSecond);
                }

                if ((this.MedianSamplesPerSecond != 0) &&
                    (this.MedianAVM != 0) &&
                    (sig.Median != null))
                {
                    ECGTool.ResampleLead(sig.Median, this.MedianSamplesPerSecond, samplesPerSecond, out sig.Median);
                }
            }

            if (this.QRSZone != null)
            {
                foreach (QRSZone zone in this.QRSZone)
                {
                    zone.Start    = (int)(((long)zone.Start * (long)samplesPerSecond) / (long)this.MedianSamplesPerSecond);
                    zone.Fiducial = (int)(((long)zone.Fiducial * (long)samplesPerSecond) / (long)this.MedianSamplesPerSecond);
                    zone.End      = (int)(((long)zone.End * (long)samplesPerSecond) / (long)this.MedianSamplesPerSecond);
                }
            }

            if ((this.RhythmSamplesPerSecond != 0) &&
                (this.RhythmAVM != 0))
            {
                this.RhythmSamplesPerSecond = samplesPerSecond;
            }

            if ((this.MedianSamplesPerSecond != 0) &&
                (this.MedianAVM != 0))
            {
                this.MedianFiducialPoint = (ushort)(((long)this.MedianFiducialPoint * (long)samplesPerSecond) / (long)this.MedianSamplesPerSecond);

                this.MedianSamplesPerSecond = samplesPerSecond;
            }
        }
示例#2
0
        /// <summary>
        /// Function to encode data in this section.
        /// </summary>
        /// <param name="data">Rhythm data to encode</param>
        /// <param name="tables">Huffman table to use during enconding</param>
        /// <param name="leadDefinition">Lead Definitions to use for encoding</param>
        /// <param name="difference">difference to use durring decoding</param>
        /// <returns>0 on succes</returns>
        public int EncodeData(short[][] data, SCPSection2 tables, SCPSection3 leadDefinition, SCPSection4 qrsLocations, int medianFreq, byte difference)
        {
            int localFreq = getSamplesPerSecond();

            if ((data != null) &&
                (tables != null) &&
                (leadDefinition != null) &&
                (localFreq > 0))
            {
                if ((medianFreq <= 0) ||
                    (medianFreq == localFreq))
                {
                    medianFreq = 1;
                    localFreq  = 1;
                }

                if ((_Bimodal == 0x1) &&
                    (qrsLocations == null))
                {
                    return(2);
                }

                ushort nrleads = leadDefinition.getNrLeads();
                _Data       = new byte[nrleads][];
                _DataLength = new ushort[nrleads];
                for (int loper = 0; loper < nrleads; loper++)
                {
                    if (data[loper] == null)
                    {
                        return(4);
                    }

                    short[] temp = data[loper];

                    int time = (leadDefinition.getLeadLength(loper) * localFreq) / medianFreq;
                    if (localFreq != medianFreq)
                    {
                        int rate = (medianFreq / localFreq);
                        // Bimodal part might be buggy unable to test.
                        if ((_Bimodal == 0x1) &&
                            ((medianFreq % localFreq) == 0) &&
                            (rate > 0) &&
                            (rate < 5))
                        {
                            // Calculate nr of samples stored in section.
                            time = 0;

                            int nrzones = qrsLocations.getNrProtectedZones();
                            for (int zone = 0; zone < nrzones; zone++)
                            {
                                int begin = (qrsLocations.getProtectedStart(zone) >= leadDefinition.getLeadStart(loper) ? qrsLocations.getProtectedStart(zone) : leadDefinition.getLeadStart(loper));
                                int end   = (qrsLocations.getProtectedEnd(zone) <= leadDefinition.getLeadEnd(loper) ? qrsLocations.getProtectedEnd(zone) : leadDefinition.getLeadEnd(loper));

                                begin = (end > begin ? end - begin + 1 : 0);

                                time += begin + (rate - (begin % rate));
                            }

                            time += ((leadDefinition.getLeadLength(loper) - time) * localFreq) / medianFreq;

                            int leadLength = leadDefinition.getLeadLength(loper);

                            time += ((leadLength - time) * localFreq) / medianFreq;

                            // Restructure bimodal data to length set in Section3.
                            temp = new short[time];

                            int time2Offset = leadDefinition.getLeadStart(loper);
                            int time1       = 0;
                            int time2       = 0;

                            while ((time1 < temp.Length) &&
                                   (time2 <= leadLength) &&
                                   (time2 < data[loper].Length))
                            {
                                int zone = 0;
                                int end  = qrsLocations.getNrProtectedZones();
                                for (; zone < end; zone++)
                                {
                                    if ((qrsLocations.getProtectedLength(zone) > 0) &&
                                        (time2 + time2Offset >= qrsLocations.getProtectedStart(zone)) &&
                                        (time2 + time2Offset <= qrsLocations.getProtectedEnd(zone)))
                                    {
                                        break;
                                    }
                                }

                                if (zone < end)
                                {
                                    temp[time1] = data[loper][time2++];
                                }
                                else
                                {
                                    int Sum = 0;
                                    for (int sumLoper = 0; sumLoper < rate; sumLoper++)
                                    {
                                        Sum += data[loper][time2 + sumLoper];
                                    }
                                    temp[time1] = (short)(Sum / rate);
                                    time2      += rate;
                                }

                                time1++;
                            }
                        }
                        else
                        {
                            _Bimodal = 0;
                            ECGTool.ResampleLead(temp, medianFreq, localFreq, out temp);
                        }
                    }

                    _Difference  = difference;
                    _Data[loper] = tables.Encode(temp, time, 0, _Difference);
                    if (_Data[loper] == null)
                    {
                        _Data       = null;
                        _DataLength = null;
                        return(8);
                    }
                    _DataLength[loper] = (ushort)_Data[loper].Length;
                    if ((_DataLength[loper] & 0x1) == 0x1)
                    {
                        _DataLength[loper]++;
                    }
                }
                return(0);
            }
            return(1);
        }
示例#3
0
        /// <summary>
        /// Function to decode data in this section.
        /// </summary>
        /// <param name="tables">Huffman table to use during deconding</param>
        /// <param name="leadDefinition"></param>
        /// <returns>decoded leads</returns>
        public short[][] DecodeData(SCPSection2 tables, SCPSection3 leadDefinition, SCPSection4 qrsLocations, int medianFreq)
        {
            int localFreq = getSamplesPerSecond();

            if (Works() &&
                (tables != null) &&
                (tables.Works()) &&
                (leadDefinition != null) &&
                (leadDefinition.Works()) &&
                (leadDefinition.getNrLeads() == _Data.Length) &&
                (localFreq > 0))
            {
                if ((medianFreq <= 0) ||
                    (medianFreq == localFreq))
                {
                    medianFreq = 1;
                    localFreq  = 1;
                }

                if ((medianFreq < localFreq) &&
                    !leadDefinition.isMediansUsed() &&
                    (_Bimodal == 0x0))
                {
                    medianFreq = localFreq;
                }

                if (((medianFreq % localFreq) != 0) ||
                    ((medianFreq / localFreq) < 1) ||
                    ((medianFreq / localFreq) > 4))
                {
                    // make an exception for ECGs that don't use compression (like corpuls ECGs that are in violation of this rule).
                    if (_Bimodal == 0x0)
                    {
                        medianFreq = localFreq;
                    }
                    else
                    {
                        return(null);
                    }
                }

                if ((_Bimodal == 0x1) &&
                    (qrsLocations == null))
                {
                    return(null);
                }

                // Reset selected table.
                tables.ResetSelect();

                short[][] leads = new short[_Data.Length][];
                for (int loper = 0; loper < _Data.Length; loper++)
                {
                    int time = (leadDefinition.getLeadLength(loper) * localFreq) / medianFreq;

                    // Bimodal part might be buggy unable to test.
                    if ((localFreq != medianFreq) &&
                        (_Bimodal == 0x1))
                    {
                        int rate = medianFreq / localFreq;

                        // Calculate nr of samples stored in section.
                        time = 0;

                        int nrzones = qrsLocations.getNrProtectedZones();
                        for (int zone = 0; zone < nrzones; zone++)
                        {
                            int begin = (qrsLocations.getProtectedStart(zone) >= leadDefinition.getLeadStart(loper) ? qrsLocations.getProtectedStart(zone) : leadDefinition.getLeadStart(loper));
                            int end   = (qrsLocations.getProtectedEnd(zone) <= leadDefinition.getLeadEnd(loper) ? qrsLocations.getProtectedEnd(zone) : leadDefinition.getLeadEnd(loper));

                            begin = (end > begin ? end - begin + 1 : 0);

                            time += begin + (rate - (begin % rate));
                        }

                        time += ((leadDefinition.getLeadLength(loper) - time) * localFreq) / medianFreq;
                    }

                    leads[loper] = tables.Decode(_Data[loper], 0, _Data[loper].Length, time, _Difference);
                    // Check if lead was decoded
                    if (leads[loper] == null)
                    {
                        // return if lead decode failed.
                        return(null);
                    }

                    if (localFreq != medianFreq)
                    {
                        int rate = medianFreq / localFreq;
                        if (_Bimodal == 0x1)
                        {
                            int beginNonProtected = 0;
                            int endNonProtected   = qrsLocations.getProtectedStart(0);
                            // Restructure bimodal data to length set in Section3.
                            short[] temp = new short[leadDefinition.getLeadLength(loper)];

                            int time1Offset = leadDefinition.getLeadStart(loper);
                            int time1       = 0;
                            int time2       = 0;

                            int zone = 0;

                            while ((time1 < temp.Length) &&
                                   (time2 < leads[loper].Length))
                            {
                                if (((time1 + time1Offset) >= beginNonProtected) &&
                                    ((time1 + time1Offset) < endNonProtected))
                                {
                                    for (int begin = 0; begin < (rate >> 1); begin++)
                                    {
                                        temp[time1++] = leads[loper][time2];
                                    }

                                    if ((time2 + ((endNonProtected - (time1 + time1Offset)) / rate)) >= leads[loper].Length)
                                    {
                                        endNonProtected -= ((time2 + ((endNonProtected - (time1 + time1Offset)) / rate)) - leads[loper].Length) * rate;
                                    }

                                    endNonProtected -= rate + (rate >> 1);

                                    while ((time1 + time1Offset) < endNonProtected)
                                    {
                                        for (int i = 0; (i < rate) && (time1 < temp.Length); i++)
                                        {
                                            temp[time1++] = (short)(((leads[loper][time2 + 1] - leads[loper][time2]) / rate) * i + leads[loper][time2]);
                                        }
                                        time2++;
                                    }

                                    endNonProtected += rate + (rate >> 1);

                                    for (int end = 0; end < (rate >> 1); end++)
                                    {
                                        temp[time1++] = leads[loper][time2];
                                    }

                                    time2++;

                                    beginNonProtected = (zone == qrsLocations.getNrProtectedZones() ? temp.Length : qrsLocations.getProtectedEnd(zone));
                                }
                                else
                                {
                                    // This should never happen!!
                                    if (zone == qrsLocations.getNrProtectedZones())
                                    {
                                        break;
                                    }

                                    while (((time1 + time1Offset) < qrsLocations.getProtectedEnd(zone)) &&
                                           (time1 < temp.Length) &&
                                           (time2 < leads[loper].Length))
                                    {
                                        temp[time1++] = leads[loper][time2++];
                                    }
                                    zone++;
                                    endNonProtected = (zone == qrsLocations.getNrProtectedZones() ? temp.Length : qrsLocations.getProtectedStart(zone));
                                }
                            }

                            leads[loper] = temp;
                        }
                        else
                        {
                            ECGTool.ResampleLead(leads[loper], localFreq, medianFreq, out leads[loper]);
                        }
                    }
                }
                return(leads);
            }
            return(null);
        }