System and method for measuring the instantaneous period of a quasi-periodic signal

• US 8,491,481 B2
• Filed: 01/29/2009
• Issued: 07/23/2013
• Est. Priority Date: 01/29/2009
• Status: Active Grant

• Current Assignee: General Electric Company
• Original Assignee: General Electric Company
• Inventors: Kabakov, Serguei
• Primary Examiner(s): Chen, Tse
• Assistant Examiner(s): Brutus, Joel F
• Application Number:
  US12/361,924
• Publication Number:
  US 20100191118A1
• Time in Patent Office: 1,636 Days
• Field of Search: 600/453, 600/509, 600/521, 600516-519, 600523-524
• US Class Current: 600/453
• CPC Class Codes: A61B 5/02411 of foetuses
• CPC Class Codes: A61B 5/344 Foetal cardiography
• CPC Class Codes: A61B 5/35 by template matching
• CPC Class Codes: A61B 5/4362 Assessing foetal parameters
• CPC Class Codes: A61B 8/02 Measuring pulse or heart rate
• CPC Class Codes: A61B 8/0866 involving foetal diagnosis; pre-natal or peri-natal diagnosis of the baby
• CPC Class Codes: A61B 8/5238 for combining image data of patient, e.g. merging several images from different acquisition modes into one image for image registration algorithms G06T7/30
• CPC Class Codes: G01R 23/15 Indicating that frequency of pulses is either above or below a predetermined value or within or outside a predetermined range of values, by making use of non-linear or digital elements (indicating that pulse width is above or below a certain limit)
• CPC Class Codes: G01S 15/52 Discriminating between fixed and moving objects or between objects moving at different speeds
• CPC Class Codes: G01S 15/88 Sonar systems specially adapted for specific applications seismic or acoustic prospecting or detecting G01V1/00

Abstract

A method of measuring the instantaneous period of a quasi-periodic signal includes determining an initial template, performing a cross-correlation between the initial template and a quasi-periodic signal and selecting portions of the signal that correspond to the peaks of the correlation signal. The selected portions are then averaged and another cross correlation is performed between the quasi-periodic signal and a template including the averaged selected portions. The instantaneous period is then measured from the correlation signal.

Claims

1. A method of measuring the instantaneous period of a quasi-periodic signal, the method comprising:

obtaining a quasi-periodic signal;

determining an initial template of a determined size and shift with a computer by;

performing a cross-correlation between the quasi-periodic signal and a template segment of the quasi-periodic signal, the template segment being of a pre-determined size and a pre-determined shift, to obtain a selection coefficient;

incrementing the pre-determined size and the pre-determined shift of the template segment to create an incremented template segment;

repeating the cross-correlation between the quasi-periodic signal and the incremented template segment through a plurality of increments to obtain selection coefficients for each of a plurality of incremented template segments; and

selecting the incremented template with a maximum selection coefficient;

performing a cross-correlation between the initial template and the quasi-periodic signal to obtain an initial correlation signal;

establishing a first correlation threshold;

selecting portions of the quasi-periodic signal that correspond to peaks of the initial correlation signal greater than the first correlation threshold;

counting a number of peaks in the initial correlation signal greater than the first correlation threshold;

averaging the selected portions to create a second template;

performing a cross-correlation between the second template and the quasi-periodic signal to obtain a second correlation signal;

establishing a second correlation threshold; and

counting the number of peaks in the second correlation signal greater than the second correlation threshold;

wherein if the number of peaks counted in the initial correlation signal is greater than or equal to the number of peaks counted in the second correlation signal, measuring the instantaneous period between the peaks in the second correlation signal and

wherein if the number of peaks counted in the initial correlation signal is less than the number of peaks counted in the second correlation signal;

selecting the portions of the quasi-periodic signal that correspond to the peaks of the second correlation signal that are greater than the second correlation threshold;

averaging the selected portions of the quasi-periodic signal to create a third template; and

performing a cross-correlation between the third template and the quasi-periodic signal to obtain a third correlation signal.

2. The method of claim 1 wherein the quasi-periodic signal comprises a first segment and a second segment and the step of determining an initial template is performed using the first segment of the quasi-periodic signal.

3. The method of claim 2 wherein the steps of performing a cross-correlation between the initial template and the quasi-periodic signal and performing a cross-correlation between the second template and the quasi-periodic signal are performed using the second segment of the quasi-periodic signal, and measuring the instantaneous period of the second segment.

4. The method of claim 1, wherein the quasi-periodic signal comprises a first segment, a second segment, and a third segment and the step of determining an initial template is performed using the first segment.

5. The method of claim 4 further comprising:

measuring the instantaneous period in the second segment by determining the initial template using the first segment, cross correlating the initial template with the second segment and using the cross-correlation between the initial template and the second segment to obtain a correlation signal of the second segment; and

measuring the instantaneous period in the third segment by determining a second initial template from the second segment, cross correlating the second initial template with the third segment, and using the cross-correlation between the second initial template and the third segment to obtain a correlation signal of the third segment.

6. The method of claim 4, further comprising:

performing a cross-correlation between the initial template and the first segment to obtain an initial correlation signal;

establishing a first correlation threshold;

selecting portions of the first segment that correspond to peaks of the initial correlation signal greater than the first correlation threshold;

averaging the selected portions to create a second template;

measuring the instantaneous period of the second segment using the second template;

performing a cross-correlation between the initial template and the second segment to obtain a second correlation signal;

establishing a second correlation threshold;

selecting portions of the second segment that correspond to peaks of the second correlation signal greater than the second correlation threshold;

averaging the selected portions to create a third template; and

measuring the instantaneous period of the third segment using the third template.

7. The method of claim 4 further comprising:

cross-correlating the initial template with the first segment to create an initial correlation signal of the first segment;

determining a first final template by;

counting the number of peaks in a previous correlation signal, wherein the previous correlation signal is initially the initial correlation signal of the first segment; and

counting a number of peaks in a current correlation signal, wherein the current correlation signal is initially the second correlation signal of the first segment;

wherein if the number of peaks counted in the previous correlation signal is greater than or equal to the number of peaks counted in the current correlation signal;

establishing the current template as the first final template; and

measuring the instantaneous period, between the peaks in the second correlation signal;

wherein if the number of peaks counted in the previous correlation signal is less than the number of peaks counted in the current correlation signal;

selecting the portions of the first segment that correspond to the peaks of the current correlation signal;

averaging the selected portions of the first segment to create a current template; and

cross-correlating the current template with the first segment to determine a new current correlation signal of the first segment;

repeating the above steps until the number of peaks counted in the current correlation signal is greater than or equal to the number of peaks counted in the previous correlation signal;

cross-correlating the first final template with the second segment to create an initial correlation signal of the second segment; and

measuring the instantaneous period in the second segment using the initial correlation signal of the second segment.

8. The method of claim 7 further comprising:

determining a second final template by;

counting the number of peaks in a previous correlation signal, wherein the previous correlation signal is initially the initial correlation signal of the second segment; and

counting a number of peaks in a current correlation signal, wherein the current correlation signal is initially the second correlation signal of the second segment;

wherein if the number of peaks counted in the previous correlation signal is greater than or equal to the number of peaks counted in the current correlation signal;

establishing the current template as the second final template; and

measuring the instantaneous period between the peaks in the current correlation signal;

wherein if the number of peaks counted in the previous correlation signal is less than the number of peaks counted in the current correlation signal;

selecting the portions of the second segment that correspond to the peaks of the current correlation signal;

averaging the selected portions of the first segment to create a current template; and

cross-correlating the current template with the second segment to determine a new current correlation signal of the second segment;

repeating the above steps until the number of peaks counted in the current correlation signal is greater than or equal to the number of peaks counted in the previous correlation signal;

cross-correlating the second final template with the third segment to create an initial correlation signal of the third segment; and

measuring the instantaneous period in the third segment using the initial correlation signal of the third segment.

9. A method of measuring an instantaneous period of a quasi-periodic signal, the method comprising:

obtaining a quasi-periodic signal;

determining an initial template of a determined size and shift with a computer by;

performing a cross-correlation between the quasi-periodic signal and a template segment of the quasi-periodic signal to obtain a selection coefficient, the template segment being of a pre-determined size and a pre-determined shift; and

incrementing the pre-determined size and a pre-determined shift of the template segment to create an incremented template segment;

repeating the cross-correlation between the quasi-periodic signal and the incremented template segment through a plurality of increments to obtain selection coefficients for each of a plurality of incremented template segments; and

selecting the incremented template with a maximum selection coefficient as the initial template;

performing a cross-correlation between an initial template and the quasi-periodic signal to obtain an initial correlation signal;

counting peaks of the initial correlation signal greater than a first correlation threshold;

selecting portions of the quasi-periodic signal that correspond to the peaks of the initial correlation signal greater than the first correlation threshold;

averaging the selected portions of the quasi-periodic signal to create a second template;

performing a cross-correlation between the second template and the quasi-periodic signal to obtain a second correlation signal; and

counting peaks of the second correlation signal greater than the correlation threshold;

wherein if the second correlation signal has more counted peaks than the initial correlation signal, using the second template as the initial template and iteratively repeating the method until the initial correlation signal has a same number or more counted peaks than the second correlation signal; and

wherein if the initial correlation signal has the same number or more counted peaks than the second correlation signal, measuring the instantaneous period from the first correlation signal.

10. The method of claim 9, wherein the quasi-periodic signal is an ultrasound signal and the instantaneous period is a fetal heart rate.