DE10125347A1 - Method for evaluation of a time sequence of discrete measurement values for e.g. electrocardiograph, in which the measurement values are divided into partial measurement series chains centered on a large or significant value - Google Patents
Method for evaluation of a time sequence of discrete measurement values for e.g. electrocardiograph, in which the measurement values are divided into partial measurement series chains centered on a large or significant valueInfo
- Publication number
- DE10125347A1 DE10125347A1 DE10125347A DE10125347A DE10125347A1 DE 10125347 A1 DE10125347 A1 DE 10125347A1 DE 10125347 A DE10125347 A DE 10125347A DE 10125347 A DE10125347 A DE 10125347A DE 10125347 A1 DE10125347 A1 DE 10125347A1
- Authority
- DE
- Germany
- Prior art keywords
- measured values
- measured value
- values
- attribute
- sequence
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005259 measurement Methods 0.000 title claims abstract description 21
- 238000011156 evaluation Methods 0.000 title claims abstract description 9
- 230000009466 transformation Effects 0.000 claims description 26
- 230000036772 blood pressure Effects 0.000 claims description 5
- 230000001360 synchronised effect Effects 0.000 claims 1
- 101100225106 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) EDC2 gene Proteins 0.000 description 14
- 238000002565 electrocardiography Methods 0.000 description 9
- 230000007774 longterm Effects 0.000 description 6
- 208000010125 myocardial infarction Diseases 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 4
- 208000010201 Exanthema Diseases 0.000 description 3
- 206010061216 Infarction Diseases 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 201000005884 exanthem Diseases 0.000 description 3
- 230000007574 infarction Effects 0.000 description 3
- 206010037844 rash Diseases 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000004083 survival effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 206010020772 Hypertension Diseases 0.000 description 1
- 108091008698 baroreceptors Proteins 0.000 description 1
- 230000035581 baroreflex Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 210000001715 carotid artery Anatomy 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007717 exclusion Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008288 physiological mechanism Effects 0.000 description 1
- 210000001774 pressoreceptor Anatomy 0.000 description 1
- 108020003175 receptors Proteins 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000003252 repetitive effect Effects 0.000 description 1
- 238000013517 stratification Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02405—Determining heart rate variability
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
- A61B5/318—Heart-related electrical modalities, e.g. electrocardiography [ECG]
- A61B5/346—Analysis of electrocardiograms
- A61B5/349—Detecting specific parameters of the electrocardiograph cycle
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D1/00—Measuring arrangements giving results other than momentary value of variable, of general application
- G01D1/14—Measuring arrangements giving results other than momentary value of variable, of general application giving a distribution function of a value, i.e. number of times the value comes within specified ranges of amplitude
Abstract
Description
Die Erfindung bezieht sich auf ein Verfahren zum Auswerten einer Folge von diskreten Meßwerten gemäß dem Oberbegriff des Patentanspruches 1.The invention relates to a method for Evaluation of a sequence of discrete measured values according to the Preamble of claim 1.
Die Untersuchung von Folgen diskreter Meßwerte ist in vielen Bereichen der Technik geläufig. So werden z. B. zur Funktionsprüfung von Maschinen mit rotierenden oder hin- und her bewegten Maschinenteilen Schwingungsmuster etwa mit Hilfe der Fourier-Analyse ausgewertet, um die in dem jeweiligen Schwingungsmuster enthaltenen Frequenzen zu bestimmen. Aus den ermittelten Frequenzen kann auf ansonsten nicht erkennbare Funktionsstörungen geschlossen werden, z. B. auf beginnende Schädigungen von Lagern etc.The investigation of sequences of discrete measured values is in familiar to many areas of technology. So z. B. for functional testing of machines with rotating or machine parts moved back and forth vibration patterns about evaluated using the Fourier analysis to get the contained in the respective vibration pattern To determine frequencies. From the determined frequencies can affect otherwise undetectable malfunctions be closed, e.g. B. on beginning damage of bearings etc.
Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zum Auswerten einer Folge von diskreten Meßwerten, so z. B. auf Schwingungsmustern, repetitiven Vorgängen oder sich zeitlich ändernden Vorgängen basierenden Meßwerten anzugeben, mit dem Parameter ermittelbar sind, die charakteristische Eigenschaften, auch verborgene, der Meßwertfolge definieren. The invention has for its object a method for evaluating a sequence of discrete measured values, e.g. B. on vibration patterns, repetitive processes or measured values based on temporally changing processes specify with which parameters can be determined, the characteristic properties, also hidden, of Define measurement sequence.
Diese Aufgabe ist gemäß der Erfindung durch die Merkmale des Anspruchs 1 gelöst.This object is according to the invention by the features of claim 1 solved.
Demnach hat das Verfahren zum Auswerten einer Folge von
diskreten Messwerten folgende Merkmale:
jedem Messwert wird ein Attribut zugeordnet, wobei
dieses Attribut für alle Messwerte die gleiche
Eigenschaft aufweist;
aus den derart mit Attributen belegten Messwerten werden
diejenigen ausgewählt, deren Attribut jeweils ein
bestimmtes gemeinsames Kriterium erfüllt, und jeweils
als bestimmter Messwert bzw. Ereignis bezeichnet;
aus der Folge der diskreten Messwerte wird eine Kette
bestimmter Länge ausgewählt;
aus dieser ausgewählten Kette werden Teilketten
gebildet, die jeweils ein ausgewähltes Ereignis als
Zentrum haben, wobei die übrigen Messwerte der
Teilketten in der ursprünglichen Folge zu beiden Seiten
des zentralen Ereignisses angeordnet werden;
die Messwerte der einzelnen Teilketten werden
entsprechend ihres Platzes in der jeweiligen Teilkette
einander zugeordnet, und zwar werden die Messwerte der
zentralen Ereignisse einander zugeordnet und die
Messwerte, die an dem gleichen Platz um das jeweilig
zentrale Ereignis stehen;
sämtliche zugeordneten Messwerte für den gleichen Platz
innerhalb der ausgewählten Teilketten werden gemittelt
und die Folge der so gemittelten Messwerte in der durch
die Kette der Messwerte vorgegebenen Reihenfolge
aufgezeichnet.Accordingly, the method for evaluating a sequence of discrete measured values has the following features:
an attribute is assigned to each measured value, this attribute having the same property for all measured values;
from the measured values thus occupied with attributes, those are selected whose attribute in each case fulfills a certain common criterion and are referred to as a certain measured value or event;
a chain of a certain length is selected from the sequence of discrete measured values;
partial chains are formed from this selected chain, each of which has a selected event as its center, the remaining measured values of the partial chains being arranged on both sides of the central event in the original sequence;
the measured values of the individual partial chains are assigned to one another in accordance with their position in the respective partial chain, specifically the measured values of the central events are assigned to one another and the measured values that are in the same place around the respective central event;
all assigned measurement values for the same location within the selected partial chains are averaged and the sequence of the measurement values averaged in this way is recorded in the order specified by the chain of measurement values.
Gemäß der Erfindung werden somit in einem ersten Schritt sämtliche Meßwerte des Signals mit einem gewissen Attribut belegt, das die Meßwerte in einer beliebigen Eigenschaft charakterisiert. Prinzipiell ist jede Art von Attribut denkbar, in einem einfachen Fall z. B. kann das Attribut aus der Differenz zum vorangegangen Meßwert bestehen. Das Attribut muss aber nicht zwangsläufig aus dem zu analysierenden Signal ableitbar sein. Es ist auch möglich, die Meßwerte des zu analysierenden Signals mit einem Attribut aus einem synchron aufgezeichneten zweiten Signal zu belegen, so z. B. die Meßwerte eines EKG-Signals mit Attributen aus einem zeitgleich aufgezeichneten Blutdrucksignal.According to the invention, in a first step all measured values of the signal with a certain Attribute proves that the measured values in any Characterized property. In principle, every type of attribute conceivable, in a simple case z. B. can the attribute from the difference to the previous measured value consist. However, the attribute does not necessarily have to be off be derivable from the signal to be analyzed. It is also possible to use the measured values of the signal to be analyzed an attribute from a synchronously recorded second signal to prove so. B. the measured values of a EKG signal with attributes from a simultaneous recorded blood pressure signal.
Aus den so mit Attributen belegten Meßwerten werden zumindest einige Meßwerte ausgewählt, deren Attribut ein gewisses Kriterium erfüllt. Wurden z. B. alle Meßwerte mit dem Attribut der Differenz zum vorangegangenen Meßwert belegt, so könnte das Auswahlkriterium sein, daß das Attribut negativ ist, d. h. der Meßwert größer als der vorangegangene ist. Die ausgewählten Meßwerte sollen im folgenden "Ereignisse" genannt werden. Für jedes Ereignis wird nun eine Kette von Meßwerten angelegt, deren zentrales Element das Ereignis selbst ist, umgeben von den Meßwerten, die im Ausgangssignal vor bzw. nach dem Ereignis angeordnet sind. Die auf diese Weise für jedes Ereignis gebildeten Ketten, deren Längen frei definierbar sind, werden "untereinander" in einer Tabelle angeordnet, wobei jede Reihe der Tabelle einer Kette von Meßwerten entspricht. Wichtig ist, daß die Meßwerte der gleichen Positionen bezüglich der Ereignisse jeweils untereinander stehen. Trifft das Kriterium auf einen relativ großen Teil der Meßwerte zu, so sind, in Abhängigkeit der gewählten Kettenlänge, Meßwerte des Ausgangssignals in der Tabelle mehrfach vertreten, d. h. die Tabelle enthält entsprechend der Methodik redundante Information. In einem nächsten Schritt werden die Mittelwerte einer jeden Säule der Tabelle gebildet, d. h. die Meßwerte der gleichen Position bezüglich des ausgewählten Meßwertes werden gemittelt.The measured values thus assigned attributes selected at least some measured values, their attribute certain criterion met. Have z. B. all measured values with the attribute of the difference to the previous one If the measured value is occupied, the selection criterion could be that the attribute is negative, d. H. the measured value is greater than the previous one is. The selected measured values should hereinafter referred to as "events". For each Event, a chain of measured values is now created, whose central element is the event itself of the measured values in the output signal before or after ordered by the event. That way for every event formed chains, their lengths free are definable, are "one below the other" in one Table arranged, each row of the table one Chain of measured values corresponds. It is important that the Measured values of the same positions with respect to the Events are one below the other. Is that true Criterion towards a relatively large part of the measured values, depending on the chosen chain length, Measured values of the output signal in the table several times represented, d. H. the table contains according to the Methodology redundant information. In a next one Step are the averages of each pillar of the Formed table, d. H. the measured values of the same Position with respect to the selected measured value averaged.
Auf diese Weise wird ein beliebig langes Ausgangssignal in eine relativ kurze Folge gemittelter Meßwerte transformiert, im folgenden "Schmidt-Bauer- Transformation (SBT)" genannt. Die SBT charakterisiert in Abhängigkeit des gewählten Attributs und des gewählten Kriteriums das beliebig lange Ausgangssignal, was einer Datenreduktion gleichkommt. Die SBT kann auf unterschiedliche Art ausgewertet werden, so z. B. in der Zeit- oder Frequenz-Domaine.In this way, an output signal of any length in a relatively short sequence of averaged measured values transformed, hereinafter "Schmidt-Bauer- Transformation (SBT) ". The SBT characterizes depending on the selected attribute and the selected criterion the output signal of any length, which is equivalent to data reduction. The SBT can be on different types are evaluated, e.g. B. in the Time or frequency domain.
Es hat sich herausgestellt, daß sich durch die mit der Erfindung vorgeschlagene Auswertung eine Vielzahl von Parametern ermitteln lässt, die aus der Folge der Meßwerte allein nicht ersichtlich ist.It has been found that the Evaluation proposed a variety of invention Lets determine parameters resulting from the sequence of Measured values alone cannot be seen.
So können z. B. aus einem Schwingungsmuster etwa einer Maschine mit drehenden bzw. hin- und her bewegbaren Maschinenteilen die dem Schwingungsmuster zugrunde liegenden Frequenzen, ähnlich wie bei einer Fourier- Analyse ermittelt werden, wobei dann auch Frequenzen ersichtlich werden, die z. B. "Störfrequenzen" sind, die auf Schädigungen der Maschine bzw. ihrer Teile hindeuten. Es hat sich gezeigt, daß mit dem Verfahren gemäß der Erfindung auch Frequenzen ermittelt werden, die bei der herkömmlichen schnellen Fourier- Transformation unentdeckt bleiben.So z. B. from a vibration pattern about one Machine with rotating or moving back and forth Machine parts that underlie the vibration pattern lying frequencies, similar to a Fourier Analysis can be determined, and then also frequencies be seen, the z. B. are "interference frequencies", the damage to the machine or its parts indicate. It has been shown that with the method frequencies are also determined according to the invention, which in the conventional fast Fourier Transformation remain undetected.
Besonders vorteilhaft kann das Verfahren gemäß der Erfindung in der Medizintechnik verwendet werden, z. B. zur Auswertung von Langzeit-EKGs, Langzeitblutdruckmessungen etc. Bezüglich der Auswertung von Langzeit-EKGs bei Infarktpatienten lässt sich das Verfahren folgendermaßen anwenden. Als Ausgangssignal dient die Folge der Schlag-zu-Schlag Intervalle des Langzeit-EKGs, welche die Grundlage der meisten Auswerteverfahren von Langzeit-EKGs ist. Jedem Meßwert wird zunächst ein Attribut zugeordnet. Als Attribut soll hier der Vergleich zum vorangegangenen Meßwert dienen, genauer der Quotient zwischen dem Meßwert selbst und dem vorangegangenen Meßwert.The method according to the Invention used in medical technology, e.g. B. for the evaluation of long-term ECGs, Long-term blood pressure measurements etc. Regarding the evaluation Long-term ECGs in infarct patients can do that Use the following procedure. As an output signal serves the sequence of beat-to-beat intervals of the Long-term ECGs, which are the basis of most Evaluation procedure of long-term ECGs is. Any reading an attribute is first assigned. As an attribute here the comparison to the previous measured value is used, more precisely the quotient between the measured value itself and the previous measurement.
Ist ein Meßwert z. B. 3% kürzer als der vorangegangene, so würde diesem Meßwert das Attribut "0.97" zugeordnet werden. Ist er um 3% länger, beträgt der Quotient 1,03. Das Kriterium für die Auswahl eines Meßwertes könnte ein Wert zwischen 1 und 1.05 sein. Nach diesem Kriterium würden alle Meßwerte als Ereignis ausgewählt, die zwar größer, aber nicht größer als 5% des vorangegangenen Meßwertes sind. In einem durchschnittlichen Langzeit-EKG entsprechen ungefähr 30% aller Meßwerte diesem Kriterium.Is a measured value z. B. 3% shorter than the previous one, the attribute "0.97" would be assigned to this measured value become. If it is 3% longer, the quotient is 1.03. The criterion for the selection of a measured value could be a Value between 1 and 1.05. According to this criterion all measured values would be selected as an event larger, but not larger than 5% of the previous one Measured value. In an average long-term ECG approximately 30% of all measured values correspond to this Criteria.
Die SBT entspricht dem durchschnittlichen Verlauf der Meßwerte vor, während und nach dem Ereignis (hier: Anstieg der Meßwerte im Ausgangssignal). Die SBT kann auf vielerlei Art quantifiziert werden: Für die Abschätzung des zukünftigen Risikos bei Patienten nach Herzinfarkt hat es sich als günstig, den zentralen Wert der SBT mit dem darauffolgenden Wert zu summieren und von dieser Summe die Summe der beiden vorangegangenen Meßwerte abzuziehen. Dieser sogenannte EDC2-Wert kann als Maß für den zentralen Anstieg angesehen werden und erlaubt eine Aussage hinsichtlich der Überlebenschance des Patienten. Je größer dieser EDC2-Wert ist, desto größer ist auch die Überlebenschance des Patienten. Der EDC2-Wert quantifiziert die SBT jedoch nur in unmittelbarer Umgebung um das Zentrum, es sind jedoch auch Parameter denkbar, die längerfristige Änderungen der SBT erfassen.The SBT corresponds to the average course of the Measured values before, during and after the event (here: Increase in the measured values in the output signal). The SBT can can be quantified in many ways: For the Estimation of future risk in patients Heart attack has turned out to be cheap, the central value sum the SBT with the subsequent value and from this sum the sum of the two previous ones Subtract measured values. This so-called EDC2 value can be regarded as a measure of the central increase and allows a statement regarding the chance of survival of the patient. The greater this EDC2 value, the more the patient's chance of survival is also greater. The However, the SBT only quantifies the EDC2 value in immediate surroundings around the center, but there are parameters also conceivable, the longer-term changes the SBT record.
Bei der gerade beschriebenen Auswertung wurde den Meßwerten des Ursprungssignals ein Attribut zugeordnet, das den Bezug zum unmittelbar vorangegangenen Meßwert herstellt, wodurch in der SBT vor allem die kürzeren Frequenzen des Ausgangssignal betont werden. Eine Betonung längerer Frequenzen kann dadurch erreicht werden, daß den Meßwerten Attribute zugeordnet werden, welche die Meßwerte über einen längeren Zeitraum vergleichen, so z. B. die Summe der Werte eines gewissen Zeitraumes nach dem Meßwert mit der Summe der Werten eines gewissen Zeitraumes vor dem Meßwert vergleichen.In the evaluation just described, the Assigned an attribute to measured values of the original signal, this is the reference to the immediately preceding measured value manufactures, which means that the shorter ones in the SBT Frequencies of the output signal are emphasized. A This can emphasize longer frequencies that attributes are assigned to the measured values, which are the measured values over a longer period of time compare, e.g. B. the sum of the values of a certain Period after the measured value with the sum of the values compare a certain period of time before the measured value.
Es ist darüber hinaus möglich, die Meßwerte einer Zeitreihe (z. B. die des EKGs) mit einem Attribut aus einer zeitgleich aufgezeichneten Zeitreihe (z. B. Blutdruckanstieg aus einer Blutdruckkurve) zu belegen. Wenn beide Zeitreihen verknüpft sind (in diesem Beispiel durch den Baroreflex: Blutdruckanstieg bewirkt Pulsverlangsamung, Blutdruckabfall bewirkt Pulsbeschleunigung), filtert die SBT deutliche Schwingungen aus dem EKG-Signal. Ist dies nicht der Fall, finden sich keine Schwingungen. Die SBT eignet sich also auch zum Nachweis bzw. Ausschluß von Koppelungen unterschiedlicher Signale.It is also possible to measure a Time series (e.g. that of the ECG) with an attribute a time series recorded at the same time (e.g. Increase in blood pressure from a blood pressure curve). If both time series are linked (in this example through the baroreflex: increased blood pressure Pulse slowdown, blood pressure drop causes Pulse acceleration), the SBT filters significantly Vibrations from the EKG signal. If not Case, there are no vibrations. The SBT is suitable thus also for the detection or exclusion of Couplings of different signals.
Weitere Ausgestaltungen der Erfindung gehen aus den Unteransprüchen hervor.Further refinements of the invention result from the Sub-claims emerge.
Die Erfindung ist in Ausführungsbeispielen anhand der Zeichnung näher erläutert. In dieser stellen dar:The invention is in exemplary embodiments based on the Drawing explained in more detail. In this represent:
Fig. 1 mit den Zeilen a bis f eine schematische Darstellung des Verfahrens zum Auswerten einer Folge von diskreten Meßwerten gemäß der Erfindung; Fig. 1 with the lines a to f is a schematic representation of the method for evaluating a sequence of discrete measured values according to the invention;
Fig. 2 ein sogenanntes Tachogramm, in dem die Zeitintervalle zwischen zwei aufeinander folgenden Herzschlägen über eine Dauer von 24 Stunden für einen Patienten aufgetragen sind; Fig. 2 is a so-called tachogram in which are plotted the time intervals between two successive heart beats over a period of 24 hours for a patient;
Fig. 3 das Frequenzspektrum des Tachogramms in Fig. 2, das über eine einfache Fourier-Transformation gebildet wurde; FIG. 3 shows the frequency spectrum of tachogram in Figure 2, which was formed through a simple Fourier transform.
Fig. 4 die Schmidt-Bauer-Transformation des Tachogramms in Fig. 2 gemäß der Erfindung; . Fig. 4 is the Schmidt-Bauer transform of tachogram in Figure 2 according to the invention;
Fig. 5 ein Ausschnitt aus Fig. 4, der die Berechnung des Parameters EDC2 veranschaulicht; FIG. 5 shows a detail from FIG. 4, which illustrates the calculation of the parameter EDC2;
Fig. 6 das Frequenzspektrum der Schmidt-Bauer- Transformation in Fig. 4 Fig. 6 shows the frequency spectrum of the Schmidt-Bauer transformation in Fig. 4
Fig. 7 die Schmidt-Bauer-Transformation ähnlich Fig. 4 eines risikogefährdeten Infarktpatienten; Fig. 7 is the Schmidt-Bauer transformation similar to Figure 4 of a risk-prone infarction patients.
Fig. 8 das Frequenzspektrum der Schmidt-Bauer- Transformation gemäß Fig. 7; Fig. 8, the frequency spectrum of the Schmidt-Bauer transformation according to Fig. 7;
Fig. 9 die Verteilungen des EDC2-Parameters für Infarktpatienten, die einen längeren Zeitraum überlebt haben, und für Infarktpatienten, die innerhalb eines gewissen Zeitraumes verstorben sind; Figure 9 shows the distributions of EDC2 parameter for heart attack patients who have survived a long time, and for heart attack patients who died within a certain period of time.
Fig. 10 eine Darstellung der Amplitude einer Schwingung über der Zeit, in diesem Falle einer Drehbeschleunigung eines Lastwagenrades mit Unwucht; FIG. 10 is a representation of the amplitude of vibration over time, in this case, a rotational acceleration of a wagon wheel with imbalance;
Fig. 11 eine einfache Fourier-Transformation der Schwingung gemäß Fig. 10; FIG. 11 shows a simple Fourier transformation of the oscillation according to FIG. 10;
Fig. 12 die Schmidt-Bauer-Transformation der Schwingung gemäß der Fig. 10 aufgrund einer Auswertung mit einem Verfahren gemäß der Erfindung; und FIG. 12 is the Schmidt-Bauer transformation of the vibration of Fig 10 as a result of analysis by a method according to the invention. and
Fig. 13 eine Fourier-Transformation der in Fig. 12 dargestellten Schmidt-Bauer-Transformation. FIG. 13 shows a Fourier transformation of the Schmidt-Bauer transformation shown in FIG. 12.
In Fig. 1 mit den Zeilen a), b), c), d), e) und f) ist das generelle Prinzip der Erfindung dargestellt.The general principle of the invention is shown in FIG. 1 with lines a), b), c), d), e) and f).
In Zeile a) ist die Folge der zu analysierenden Meßwerte aufgetragen.Line a) is the sequence of the measured values to be analyzed applied.
Gemäß Zeile b) wird jedem Meßwert ein gewisses frei definierbares Attribut zugeordnet, z. B. die Differenz zum vorangegangenen Meßwert.According to line b) a certain amount is released for each measured value assignable attribute assigned, e.g. B. the difference to the previous measured value.
Gemäß Zeile c) wird ein Kriterium definiert, das auf die den Meßwerten zugeordneten Attribute angewendet wird. In diesem Beispiel ist das Kriterium dann erfüllt (K = 1), wenn das Attribut positiv ist (A < 0). Meßwerte, deren Attribute das definierte Kriterium erfüllen, werden als Ereignisse definiert.According to line c), a criterion is defined that applies to the attributes assigned to the measured values are used. In In this example the criterion is then fulfilled (K = 1), if the attribute is positive (A <0). Readings whose Attributes that meet the defined criteria are marked as Events defined.
Gemäß Zeile d) werden die Ereignisse mitsamt jeweils vorangegangenen und folgenden Meßwerten als Ereignis- Ketten in eine Tabelle eingetragen, in der Art, daß Meßwerte der gleichen Position zu dem Ereignis untereinander stehen.According to line d), the events together with each previous and following measured values as event Chains entered in a table in such a way that Measured values of the same position for the event stand with each other.
Gemäß Zeile e) wird die Schmidt-Bauer-Transformation durch Mittelung der Ereignis-Ketten bestimmt. Hierdurch werden die Ereignisse selbst sowie die Meßwerte der jeweils gleichen Position bezüglich des Ereignisses gemittelt.According to line e), the Schmidt-Bauer transformation determined by averaging the event chains. hereby the events themselves as well as the measured values of the same position with respect to the event averaged.
In Fig. 2 ist ein Tachogramm, d. h. die Folge der Zeitintervalle zwischen jeweils zwei aufeinanderfolgenden Herzschlägen (= RR-Intervalle oder RRI), in diesem Falle ungefähr 80000 Zeitintervalle, gegen ihre Nummer im EKG aufgetragen. Dieses Tachogramm stammt von einem Patienten nach einem Myokardinfarkt, der diesen über zwei Jahre überlebt hat.In Fig. 2 is a tachogram, the sequence of time intervals that is, between each two sequential heartbeats (= RR intervals or RRI), in this case approximately 80,000 time intervals, plotted against their number in the ECG. This tachogram is from a patient after a myocardial infarction who has survived it for over two years.
In Fig. 3 ist das Frequenzspektrum des in Fig. 2 gezeigten Tachogramms dargestellt, das mit Hilfe einer einfachen Fourier-Transformation als Funktion der Amplitude in Abhängigkeit der Frequenz ermittelt wurde. FIG. 3 shows the frequency spectrum of the tachogram shown in FIG. 2, which was determined using a simple Fourier transformation as a function of the amplitude as a function of the frequency.
In Fig. 4 ist die Schmidt-Bauer-Transformation des Tachogramms in Fig. 2 gemäß der Erfindung dargestellt. Als Attribut wurde der Quotient zwischen dem Meßwertes selbst und dem vorangegangenen Meßwert, als Kriterium ein Wert zwischen 1 und 1.5 verwendet. Die Kettenlänge der Schmidt-Bauer-Transformation wurde hier auf 200 festgesetzt, d. h. 100 Werte vor und 100 Werte nach dem zentralen Ereignis. Die Werte der Transformation schwanken um den Wert 1000 ms, wobei um das Zentrum bei dem Wert Null eine erhebliche Dynamik festzustellen ist. Es lassen sich schnelle Modulationen der Schlag-zu- Schlag Intervalle von langsameren abgrenzen.In FIG. 4, the Schmidt-Bauer transform of tachogram in Fig. 2 is shown according to the invention. The quotient between the measured value itself and the previous measured value was used as an attribute, and a value between 1 and 1.5 was used as the criterion. The chain length of the Schmidt-Bauer transformation was set here to 200, ie 100 values before and 100 values after the central event. The values of the transformation fluctuate around the value 1000 ms, whereby considerable dynamics can be determined around the center at the value zero. Fast modulations of the beat-to-beat intervals can be distinguished from slower ones.
Fig. 5 zeigt einen Ausschnitt um das Zentrum der
Schmidt-Bauer-Transformation gemäß Fig. 4 und
veranschaulicht die Berechnung des Parameters EDC2, der
zur Risikostratifizierung verwendet wird. Eingetragen
sind die Werte RRI(-2), RRI(-1) sowie RRI(+1) um das
zentrale Ereignis RRI(0). Aus diesen Werten kann ein
charakteristischer Wert EDC2 errechnet werden, nämlich
FIG. 5 shows a section around the center of the Schmidt-Bauer transformation according to FIG. 4 and illustrates the calculation of the parameter EDC2, which is used for risk stratification. The values RRI (-2), RRI (-1) and RRI (+1) are entered around the central event RRI (0). A characteristic value EDC2 can be calculated from these values, namely
EDC2 = (RRI(0) + RRI(+1)) - (RRI(-1) + RRI(-2))EDC2 = (RRI (0) + RRI (+1)) - (RRI (-1) + RRI (-2))
Mit den Werten der dargestellten Kurven errechnet sich hierbei ein Wert von EDC2 = 32 ms.The values of the curves shown are calculated here a value of EDC2 = 32 ms.
In Fig. 6 ist eine einfache Fourier-Transformation der Schmidt-Bauer-Transformation gemäß Fig. 4 dargestellt. Man erkennt deutliche Ausschläge bei ungefähr 0.003 Hz, bei 0.03 Hz und bei 0.3 Hz, welche in der einfachen Fourier-Transformation des Ausgangssignals nicht erkennbar sind. Diese drei Frequenzen entsprechen bekanntermaßen verschiedenen physiologischen Mechanismen der Herzfrequenzregulation. So entspricht z. B. der Ausschlag bei 0.3 Hz der Atemmodulation des Herzschlages. Weiterhin erkennbar ist ein zusätzlicher Ausschlag zwischen 0.003 Hz und 0.03 Hz, welcher bisher nicht bekannt ist. FIG. 6 shows a simple Fourier transformation of the Schmidt-Bauer transformation according to FIG. 4. You can see clear swings at around 0.003 Hz, at 0.03 Hz and at 0.3 Hz, which are not recognizable in the simple Fourier transform of the output signal. These three frequencies are known to correspond to different physiological mechanisms of heart rate regulation. So corresponds to B. the rash at 0.3 Hz of the breath modulation of the heartbeat. There is also an additional deflection between 0.003 Hz and 0.03 Hz, which is not yet known.
In Fig. 7 ist eine Schmidt-Bauer-Transformation eines Tachogramms eines Patienten, der innerhalb von zwei Jahren nach Herzinfarkt plötzlich verstorben ist, dargestellt. Man erkennt das die Dynamik um das zentrale Ereignis deutlich eingeschränkt ist. Der EDC2-Wert fällt mit 6 ms deutlich geringer aus.In Fig. 7 is a Schmidt-Bauer Transformation of a speedometer gram of a patient who died within two years after a heart attack suddenly shown. You can see that the dynamic around the central event is clearly limited. The EDC2 value is significantly lower at 6 ms.
Fig. 8 zeigt die Fourier-Transformation der Schmidt- Bauer-Transformation gemäß Fig. 7 (zum Vergleich ist die Fourier-Transformation gemäß Fig. 6 als gepunktete Linie dargestellt). Der in Fig. 6 erkennbare Ausschlag bei 0.3 Hz ist nicht vorhanden, dafür jedoch ein Ausschlag bei ungefähr 0.75 Hz. Dieser entspricht der Modulation des Herzschlages durch die Druckrezeptoren der Halsschlagader, die sogenannten Barorezeptoren. FIG. 8 shows the Fourier transformation of the Schmidt-Bauer transformation according to FIG. 7 (for comparison, the Fourier transformation according to FIG. 6 is shown as a dotted line). The rash at 0.3 Hz shown in FIG. 6 is not present, but there is a rash at about 0.75 Hz. This corresponds to the modulation of the heartbeat by the pressure receptors of the carotid artery, the so-called baroreceptors.
Fig. 9 zeigt die Unterschiede der Verteilung der an über 1000 Patienten berechneten EDC2-Werte für Patienten, die ein Überwachungsintervall von zwei Jahren überlebt haben (Gruppe 0), und für Patienten, die während dieses Zeitraumes verstorben sind (Gruppe 1). Man erkennt, daß die EDC2 Wert der Patienten, die versterben, deutlich niedriger sind als die EDC2-Werte der Patienten, die diesen Zeitraum überleben. Die signifikante Vorhersagekraft der EDC2-Werte bisher an über 4000 Infarktpatienten getestet und bewiesen worden. FIG. 9 shows the differences in the distribution of EDC2 values calculated for over 1000 patients for patients who survived a monitoring interval of two years (group 0) and for patients who died during this period (group 1). It can be seen that the EDC2 values of the patients who died are significantly lower than the EDC2 values of the patients who survived this period. The significant predictive power of EDC2 values has been tested and proven on over 4000 infarction patients.
Die Fig. 10 zeigt das Zeitverhalten der Amplitude einer Schwingung, in diesem Falle der Schwingung im Zusammenhang mit der Drehbeschleunigung eines Lastwagenrades mit einer Unwucht. Aufgetragen ist in dieser Figur die Amplitude in Abhängigkeit der Zeit. Fig. 10 shows the timing of the amplitude of vibration in this case, the vibration in connection with the rotational acceleration of a wagon wheel with an unbalance. The amplitude as a function of time is plotted in this figure.
In Fig. 11 ist die Fourier-Transformation der Originalschwingung entsprechend Fig. 8 aufgetragen.The Fourier transformation of the original vibration corresponding to FIG. 8 is plotted in FIG. 11.
In Fig. 12 ist die Schmidt-Bauer-Transformation der Schwingung gemäß Fig. 10 dargestellt. Als Attribut wurde die Differenz zum vorangegangenen Meßwert, als Kriterium eine Wert über null verwendet. FIG. 12 shows the Schmidt-Bauer transformation of the vibration according to FIG. 10. The difference to the previous measured value was used as an attribute and a value above zero was used as the criterion.
In Fig. 13 ist die Fourier-Transformation der Schmidt- Bauer-Transformation gemäß Fig. 10 dargestellt. Man sieht, daß diese Darstellung im wesentlichen mit dem Diagramm gemäß Fig. 11 übereinstimmt. FIG. 13 shows the Fourier transformation of the Schmidt-Bauer transformation according to FIG. 10. It can be seen that this representation essentially corresponds to the diagram according to FIG. 11.
Ein entscheidender Vorteil ist, daß die Transformation ein relativ kurzes Signal darstellt, das jedoch im Vergleich zu langen Ausgangssignal keine Information verliert.A key advantage is that the transformation is a relatively short signal, but in No information compared to long output signal loses.
Claims (16)
jedem Messwert wird ein Attribut zugeordnet, wobei dieses Attribut für alle Messwerte die gleiche Eigenschaft aufweist;
aus den derart mit Attributen belegten Messwerten werden diejenigen ausgewählt, deren Attribut jeweils ein bestimmtes gemeinsames Kriterium erfüllt, und jeweils als bestimmter Messwert bzw. Ereignis bezeichnet;
aus der Folge der diskreten Messwerte wird eine Kette bestimmter Länge ausgewählt;
aus dieser ausgewählten Kette werden Teilketten gebildet, die jeweils ein ausgewähltes Ereignis als Zentrum haben, wobei die übrigen Messwerte der Teilketten in der ursprünglichen Folge zu beiden Seiten des zentralen Ereignisses angeordnet werden;
die Messwerte der einzelnen Teilketten werden entsprechend ihres Platzes in der jeweiligen Teilkette einander zugeordnet, und zwar werden die Messwerte der zentralen Ereignisse einander zugeordnet und die Messwerte, die an dem gleichen Platz um das jeweilig zentrale Ereignis stehen;
sämtliche zugeordneten Messwerte für den gleichen Platz innerhalb der ausgewählten Teilketten werden gemittelt und die Folge der so gemittelten Messwerte in der durch die Kette der Messwerte vorgegebenen Reihenfolge aufgezeichnet.1. Method for evaluating a sequence of discrete measured values, with the following features:
an attribute is assigned to each measured value, this attribute having the same property for all measured values;
from the measured values thus occupied with attributes, those are selected whose attribute in each case fulfills a certain common criterion and are referred to as a certain measured value or event;
a chain of a certain length is selected from the sequence of discrete measured values;
partial chains are formed from this selected chain, each of which has a selected event as its center, the remaining measured values of the partial chains being arranged on both sides of the central event in the original sequence;
the measured values of the individual partial chains are assigned to one another in accordance with their position in the respective partial chain, specifically the measured values of the central events are assigned to one another and the measured values that are in the same place around the respective central event;
all assigned measurement values for the same location within the selected partial chains are averaged and the sequence of the measurement values averaged in this way is recorded in the order specified by the chain of measurement values.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE10125347A DE10125347B4 (en) | 2001-05-23 | 2001-05-23 | Method for evaluating a sequence of discrete measured values |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE10125347A DE10125347B4 (en) | 2001-05-23 | 2001-05-23 | Method for evaluating a sequence of discrete measured values |
Publications (2)
Publication Number | Publication Date |
---|---|
DE10125347A1 true DE10125347A1 (en) | 2002-12-12 |
DE10125347B4 DE10125347B4 (en) | 2007-03-15 |
Family
ID=7685996
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE10125347A Expired - Fee Related DE10125347B4 (en) | 2001-05-23 | 2001-05-23 | Method for evaluating a sequence of discrete measured values |
Country Status (1)
Country | Link |
---|---|
DE (1) | DE10125347B4 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009117996A1 (en) * | 2008-03-27 | 2009-10-01 | Technische Universität München | Method and device for determining the bodily fitness of a person |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102013000704B4 (en) | 2013-01-17 | 2017-02-02 | FZI Forschungszentrum Informatik | Method and measuring arrangement for determining functional parameters of the heart rate response |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4006948A1 (en) * | 1990-03-06 | 1991-09-12 | Bayerische Motoren Werke Ag | Wear monitoring system for cyclically loaded component - compares static distribution difference signals at end of each loading cycle with stored reference signals |
DE4039648A1 (en) * | 1990-12-12 | 1992-07-16 | Rolf Wendler | Measurement value processing system for biological object - mathematically evaluates and compares with given measurement value structures |
DE4231314A1 (en) * | 1991-09-20 | 1993-04-01 | Nippon Selen Co Ltd | Spinner yarn break monitor - has detectors at each station with paired electrodes to give induced signals from static electricity at each yarn for data processing |
DE4134960A1 (en) * | 1991-10-23 | 1993-04-29 | Popp Fritz Albert Dr | METHOD FOR A HOLISTIC ANALYSIS OF THE HEALTH CONDITION |
-
2001
- 2001-05-23 DE DE10125347A patent/DE10125347B4/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE4006948A1 (en) * | 1990-03-06 | 1991-09-12 | Bayerische Motoren Werke Ag | Wear monitoring system for cyclically loaded component - compares static distribution difference signals at end of each loading cycle with stored reference signals |
DE4039648A1 (en) * | 1990-12-12 | 1992-07-16 | Rolf Wendler | Measurement value processing system for biological object - mathematically evaluates and compares with given measurement value structures |
DE4231314A1 (en) * | 1991-09-20 | 1993-04-01 | Nippon Selen Co Ltd | Spinner yarn break monitor - has detectors at each station with paired electrodes to give induced signals from static electricity at each yarn for data processing |
DE4134960A1 (en) * | 1991-10-23 | 1993-04-29 | Popp Fritz Albert Dr | METHOD FOR A HOLISTIC ANALYSIS OF THE HEALTH CONDITION |
Non-Patent Citations (1)
Title |
---|
DE-BUCH: PROFOS, P.: Handbuch der industriellen Meßtechnik, Essen, Vulkan-Verlag, Dr. W. Classen Nachf. GmbH & Co. KG, 1978, S. 168-189 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009117996A1 (en) * | 2008-03-27 | 2009-10-01 | Technische Universität München | Method and device for determining the bodily fitness of a person |
Also Published As
Publication number | Publication date |
---|---|
DE10125347B4 (en) | 2007-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP0355506B1 (en) | Arrangement for measuring local bioelectric currents in biological tissue | |
DE60210086T2 (en) | REAL TIME HIGH FREQUENCY QRS ELECTROCARDIOGRAPH | |
DE3843714C2 (en) | Method for determining electrocardiographic signals and arrangement for reproducing electrocardiographic signals | |
DE60124689T2 (en) | Monitoring a patient's condition under anesthesia or sedation | |
DE60114658T2 (en) | MONITORING DEVICE FOR MONITORING DIASTOLIC RELAXATION BY IMPEDANCE MEASUREMENT | |
EP1047987A2 (en) | Method and device for representing and monitoring functional parameters of a physiological system | |
EP1059875B1 (en) | Evaluation system for obtaining diagnostic information from the signals and data of medical sensor systems | |
CH632848A5 (en) | DEVICE FOR DETECTING SIGNALS, ESPECIALLY FETAL HEART SIGNALS. | |
DE2449606B2 (en) | Arrangement for the automatic processing of electrical heart action signals | |
WO2017050321A1 (en) | Method and device for quantifying a respiratory sinus arrhythmia and use of said type of method or said type of device | |
DE60204781T2 (en) | SYSTEM FOR THE ANALYSIS OF BIOLOGICAL SIGNAL DATA OF INTRA-CRANIAL PRESSURE AND BLOOD PRESSURE | |
DE3927709A1 (en) | HEART MONITORING DEVICE | |
WO1999027463A2 (en) | Method for determining at least one diagnostic piece of information from signal patterns of medical sensor systems | |
EP0150352B1 (en) | Method of determining the beginning and the end of a physiological signal displayed as a closed loop | |
DE4405827A1 (en) | ECG data compression in implanted monitor system | |
DE3511697C2 (en) | ||
DE10125347A1 (en) | Method for evaluation of a time sequence of discrete measurement values for e.g. electrocardiograph, in which the measurement values are divided into partial measurement series chains centered on a large or significant value | |
WO2004049190A2 (en) | Method for evaluating a sequence of discrete readings | |
EP1106205A2 (en) | Method of operation for an implanted cardiac device, in particular a cardiac pacemaker | |
EP1175865A2 (en) | Evaluation method for determining the temporal stationarity of measured physiological signals | |
DE4427991C2 (en) | Method and device for measuring and displaying changes in blood pressure | |
DE2744643B2 (en) | Processing of cardiac action signals stored in tape | |
EP1124186A2 (en) | Method for computing the heartrate variability for use in an ECG-monitor and an ECG-monitor comprising such a computational program | |
EP1192897B1 (en) | Risk monitoring | |
DE2516167A1 (en) | METHOD AND DEVICE FOR THE QUANTITATIVE DETERMINATION OF THE OCCURRENCE OF DELAYS IN A FOETAL HEARTBEAT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
OP8 | Request for examination as to paragraph 44 patent law | ||
8364 | No opposition during term of opposition | ||
R119 | Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee |