CN106821420B - The apparatus and method for carrying out the calibration of heart rate period are composed based on ultrasonic power - Google Patents
The apparatus and method for carrying out the calibration of heart rate period are composed based on ultrasonic power Download PDFInfo
- Publication number
- CN106821420B CN106821420B CN201710159512.6A CN201710159512A CN106821420B CN 106821420 B CN106821420 B CN 106821420B CN 201710159512 A CN201710159512 A CN 201710159512A CN 106821420 B CN106821420 B CN 106821420B
- Authority
- CN
- China
- Prior art keywords
- maximum
- minimum
- candidate
- doppler
- ultrasonic power
- 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.)
- Active
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/02—Measuring pulse or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Medical Informatics (AREA)
- Physics & Mathematics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cardiology (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
The present invention provides the apparatus and method for carrying out the calibration of heart rate period are composed based on ultrasonic power, it is related to clinical diagnosis field.This passes through the curved section for intercepting a cardiac cycle of Doppler's envelope curve based on the apparatus and method that ultrasonic power composes progress heart rate period calibration, extract the multiple candidate maximum and multiple candidate minimums of the curved section of the cardiac cycle of interception, according to calculating true minimum in the multiple candidate maximum and multiple candidate minimums, to eliminate the multiple candidate maximum and candidate minimum of the generation of the jitter as caused by noise or outside cause, improve the accuracy of period calculating, and it can be according to true maximum value and true minimum to the starting point of Doppler's envelope curve in a cycle, end point is demarcated, so that doctor is more intuitive to the observation of Doppler's envelope curve, conclusion to obtain is more efficient.
Description
Technical field
The present invention relates to clinical diagnosis fields, carry out heart rate period mark based on ultrasonic power spectrum in particular to one kind
Fixed apparatus and method.
Background technique
Currently, heart rate detection important in inhibiting in clinical diagnosis, to the detection of pregnant woman's perinatal period Fetal Heart Rate to raising
Delivery of fetus quality can reduce fetal mortality and abnormal rate.In ultrasonic doppler blood flow measurement, heart rate value and heart rate week
Diagnosis important in inhibiting of the calibration of phase for vascular diseases.Ultrasonic Doppler is detected due to lossless and highly sensitive
Feature has been widely used.By pressing on the organ-tissue (such as heart, arteria carotis) for being placed on human body that will pop one's head in
According to doppler principle, the frequency that the movement of reflecting interface will lead to echo changes compared with the frequency of source signal.Ultrasonic signal is met
To mobile organ-tissue, then Doppler frequency shift can be generated, and those metastable organ-tissues then will not.In addition, due to
It in ultrasonic doppler blood flow measurement, not only needs to calculate heart rate value, but also needs in Doppler's envelope curve (maximum frequency song
Line) in calibrate each cardiac cycle, facilitate doctor to check, accurately to understand the specific feelings of detected person using doctor is more wired
Condition.
In the prior art, when composing progress cycle detection to ultrasonic power, since the factors such as noise influence, it is easy to appear mistake
Difference, for example, the problem of period occur and halving or double.Simultaneously as signal width when being detected using prior art manner
The variation of value is too sensitive, and sign mutation leads to jitter once in a while as caused by noise or outside cause, it is easy to go out
The problem of mistake is now calculated, in addition there is presently no the method for calibration cardiac cycle, doctor is composed according to ultrasonic power to human body
The judgement of concrete condition is not intuitive enough.
Summary of the invention
In view of this, being designed to provide for the embodiment of the present invention is a kind of based on ultrasonic power spectrum progress heart rate period calibration
Apparatus and method.
In a first aspect, the embodiment of the invention provides a kind of device for composing progress heart rate period calibration based on ultrasonic power,
It is described based on ultrasonic power compose carry out the calibration of heart rate period device include:
Envelope recognition unit goes out Doppler's envelope curve of the ultrasonic power spectrum of input for identification；
Period computing unit, for calculating cardiac cycle according to Doppler's envelope curve；
Interception unit, the curved section of the cardiac cycle for intercepting Doppler's envelope curve；
Data extracting unit, multiple candidate maximum of the curved section of the cardiac cycle for extracting interception and multiple times
Select minimum；
Minimum computing unit, for true according to being calculated in the multiple candidate maximum and multiple candidate minimums
Minimum；
As a result output unit, for exporting the true minimum.
Second aspect, the embodiment of the invention also provides a kind of sides that progress heart rate period calibration is composed based on ultrasonic power
Method, the method for composing progress heart rate period calibration based on ultrasonic power include:
Identify Doppler's envelope curve of the ultrasonic power spectrum of input；
Cardiac cycle is calculated according to Doppler's envelope curve；
Intercept the curved section of a cardiac cycle of Doppler's envelope curve；
Extract the multiple candidate maximum and multiple candidate minimums of the curved section of the cardiac cycle of interception；
True minimum is calculated according to the multiple candidate maximum and multiple candidate minimums；
The true minimum is exported.
Compared with prior art, provided by the invention that device and the side for carrying out the calibration of heart rate period are composed based on ultrasonic power
Method, the curved section of the cardiac cycle by intercepting Doppler's envelope curve, extracts the song of the cardiac cycle of interception
The multiple candidate maximum and multiple candidate minimums of line segment, according to the multiple candidate maximum and multiple candidate minimums
It is middle to calculate true minimum, to eliminate multiple candidate poles of the generation of the jitter as caused by noise or outside cause
Big value and candidate minimum, improve the accuracy of period calculating, and can be according to true maximum value and true minimum to one
Starting point, the end point of Doppler's envelope curve in a period are demarcated, so that sight of the doctor to Doppler's envelope curve
It examines more intuitively, so that the conclusion obtained is more efficient.
To enable the above objects, features and advantages of the present invention to be clearer and more comprehensible, preferred embodiment is cited below particularly, and cooperate
Appended attached drawing, is described in detail below.
Detailed description of the invention
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.The present invention being usually described and illustrated herein in the accompanying drawings is implemented
The component of example can be arranged and be designed with a variety of different configurations.Therefore, below to the reality of the invention provided in the accompanying drawings
The detailed description for applying example is not intended to limit the range of claimed invention, but is merely representative of selected implementation of the invention
Example.Based on the embodiments of the present invention, obtained by those of ordinary skill in the art without making creative efforts
Every other embodiment, shall fall within the protection scope of the present invention.
Fig. 1 is the structural block diagram of server provided in an embodiment of the present invention；
Fig. 2 is the functional unit provided in an embodiment of the present invention composed based on ultrasonic power and carry out the device of heart rate period calibration
Schematic diagram；
Fig. 3 is the schematic diagram of gradient function provided in an embodiment of the present invention；
Fig. 4 is the schematic diagram of accumulative variance function provided in an embodiment of the present invention；
Fig. 5 is that the starting point of a cycle provided in an embodiment of the present invention to Doppler's envelope curve and end point carry out
The schematic diagram of calibration；
Fig. 6 is the flow chart provided in an embodiment of the present invention composed based on ultrasonic power and carry out the method for heart rate period calibration.
Icon: 100- composes the device for carrying out the calibration of heart rate period based on ultrasonic power；101- processor；102- memory；
103- storage control；104- Peripheral Interface；200- server；201- envelope recognition unit；202- period computing unit；203-
Interception unit；204- data extracting unit；205- minimum computing unit；206- result output unit.
Specific embodiment
Below in conjunction with attached drawing in the embodiment of the present invention, technical solution in the embodiment of the present invention carries out clear, complete
Ground description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Usually exist
The component of the embodiment of the present invention described and illustrated in attached drawing can be arranged and be designed with a variety of different configurations herein.Cause
This, is not intended to limit claimed invention to the detailed description of the embodiment of the present invention provided in the accompanying drawings below
Range, but it is merely representative of selected embodiment of the invention.Based on the embodiment of the present invention, those skilled in the art are not doing
Every other embodiment obtained under the premise of creative work out, shall fall within the protection scope of the present invention.
The apparatus and method for composing progress heart rate period calibration based on ultrasonic power provided by present pre-ferred embodiments can
Applied to server.The server may be, but not limited to, network server, database server, cloud server etc..
Fig. 1 shows a kind of structural block diagram of server 200 that can be applied in the embodiment of the present invention.The server-side packet
Include composed based on ultrasonic power carry out the device 100 of heart rate period calibration, processor 101, memory 102, storage control 103 and
Peripheral Interface 104.
The memory 102, storage control 103 and processor 101, each element are directly or indirectly electrical between each other
Connection, to realize the transmission or interaction of data.For example, these elements can pass through one or more communication bus or letter between each other
Number line, which is realized, to be electrically connected.It is described composed based on ultrasonic power carry out the calibration of heart rate period device 100 include at least one can be with
The form of software or firmware (firmware) is stored in the memory 102 or is solidificated in the operation system of the client 100
Software function module in system (operating system, OS).The processor 101 stores in memory 102 for executing
Executable module, for example, described compose the software function mould that includes of device 100 for carrying out the calibration of heart rate period based on ultrasonic power
Block or computer program.
Wherein, memory 102 may be, but not limited to, random access memory (Random Access Memory,
RAM), read-only memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only
Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM),
Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..
Wherein, memory 102 is for storing program, and the processor 101 executes described program after receiving and executing instruction, aforementioned
Method performed by the server-side that the stream process that any embodiment of the embodiment of the present invention discloses defines can be applied to processor 101
In, or realized by processor 101.
Processor 101 may be a kind of IC chip, the processing capacity with signal.Above-mentioned processor 101 can
To be general processor, including central processing unit (Central Processing Unit, abbreviation CPU), network processing unit
(Network Processor, abbreviation NP) etc.；Can also be digital signal processor (DSP), specific integrated circuit (ASIC),
Ready-made programmable gate array (FPGA) either other programmable logic device, discrete gate or transistor logic, discrete hard
Part component.It may be implemented or execute disclosed each method, step and the logic diagram in the embodiment of the present invention.General processor
It can be microprocessor or the processor 101 be also possible to any conventional processor 101 etc..
Various input/output devices are couple processor 101 and memory 102 by Peripheral Interface 104.In some implementations
In example, Peripheral Interface 104, processor 101 and storage control 103 can be realized in one single chip.In some other reality
In example, they can be realized by independent chip respectively.
It is appreciated that structure shown in FIG. 1 is only to illustrate, server 200 may also include than shown in Fig. 1 more or more
Few component, or with the configuration different from shown in Fig. 1.Each component shown in Fig. 1 can use hardware, software or its group
It closes and realizes.
Referring to Fig. 2, a kind of device that progress heart rate period calibration is composed based on ultrasonic power provided in an embodiment of the present invention,
It is described based on ultrasonic power compose carry out the calibration of heart rate period device include envelope recognition unit 201, period computing unit 202,
Interception unit 203, data extracting unit 204, minimum computing unit 205 and result output unit 206.
The envelope recognition unit 201 goes out Doppler's envelope curve of the ultrasonic power spectrum of input for identification.
In the present embodiment, the frequency spectrum that ultrasonic power spectrum is cardiac electrical cycle continuous signal x (t) is represented by discrete non-
Periodic sequence P (f), the sequence that square │ P (f) │ 2 of its amplitude spectrum is lined up, is just referred to as the " function of the periodic signal
Rate spectrum "；The amplitude of ultrasonic power spectrum is changed by low-frequency modulation signal, and the peak dot that ultrasonic power is composed is connected, so that it may
To obtain a curve corresponding with low-frequency modulation signal, this curve is exactly Doppler's envelope curve.Specifically, the packet
Network recognition unit 201 is used for according to percentage algorithm or geometric algorithm or improves geometric algorithm or thresholding algorithm or wavelet transformation calculation
Method or dynamic programming algorithm identify Doppler's envelope curve of the ultrasonic power spectrum of input.For example, Wavelet Transformation Algorithm
(wavelet transform, WT) is a kind of new signal transform analysis method, the think of to be localized using short time discrete Fourier transform
Think, be capable of providing " T/F " window with frequency shift, is the ideal work for carrying out signal time frequency analysis and processing
Tool.Of course, identify that the mode of Doppler's envelope curve of the ultrasonic power spectrum of input only can use the algorithm above,
It is merely illustrative herein.
The period computing unit 202 is used to calculate cardiac cycle according to Doppler's envelope curve.
In the present embodiment, calculating cardiac cycle includes but is not limited to following two mode:
As shown in figure 4, the first: the period computing unit 202 be used to calculate Doppler's envelope curve from phase
Function is closed, then calculates the accumulative variance function of auto-correlation function, first maximum point by accumulative variance function apart from zero point
Corresponding time point is as the period.
Specifically, according to formulaOrCalculate the auto-correlation letter of Doppler's envelope curve
Number R (n), wherein R (n) is auto-correlation function, and F (m) is Doppler's envelope curve, and N is the data length of Maximum frequency curve, μ
For the average value of F (m), calculated auto-correlation function needs to utilize accumulative variance function there may be certain error at this time
Or gradient function is modified auto-correlation function, when being modified using accumulative variance function to auto-correlation function, specifically
Mode be according to formulaCalculate accumulative variance function D (n)；
As shown in figure 3, auto-correlation function R (n) is defined as when being modified using gradient function to auto-correlation function
Function f (x, y), auto-correlation function R (n) are defined as function f (x, y) has the continuous partial derivative of single order in plane domain G, then right
In every bit P (x, y) ∈ G, a vector can be all madeThe vector is known as f (x, y) in the gradient of point P (x, y),
Gradient function is denoted as:
Wherein, gradient direction turns tangent of an angle to gradient for x-axisWith the side for obtaining maximum directional derivative
To it is consistent when, mould be directional derivative maximum value.
The local maximum that the present invention uses gradient rise method to find a function, the iterative process of algorithm is the mistake of one " upward slope "
Journey, each step selection maximum direction of change of slope are up walked, this direction is exactly function in this point gradient direction, finally
With iterations going on, gradient or constantly reduction, finally approach and zero.
The mode of calculating cycle can be with are as follows: on accumulative variance function D (n), to the forward lookup of X-axis the since zero point
One maximum point, the distance of the maximum point to zero point are regarded as periodic quantity, i.e. cardiac cycle T', then multiplied by super
Cardiac cycle T can be obtained in time interval between every spectral line of power sound spectrum, in addition, the mode of calculating cycle can be with
Are as follows: on gradient function, to first maximum point of the forward lookup of X-axis since zero point, the maximum point to zero point away from
From periodic quantity, i.e. cardiac cycle T' is regarded as, then multiplied by the time interval between every spectral line of ultrasonic power spectrum
Cardiac cycle T can be obtained.
Auto-correlation function is modified by the accumulative variance function or gradient function that calculate auto-correlation function, is obtained
Cardiac cycle T is more accurate effectively.In addition, heart rate refers under normal person's rest state by that can calculate heart rate cardiac cycle
The number of heartbeat per minute, is also quiet heart rate, generally 60~100 beats/min, can be because of age, gender or other physiologic factors
Individual difference is generated, since the unit in period is the second, the specific calculation of heart rate isWherein, HR is the heart
Rate.
The interception unit 203 is used to intercept the curved section of a cardiac cycle of Doppler's envelope curve.
Preferably, in order to obtain newest ultrasonic power modal data, interception unit 203 intercepts Doppler's envelope curve
A newest cardiac cycle curved section.
The data of a newest cardiac cycle are denoted as F'(f), wherein F'(f)=F (t+f), f ∈ [1, T], t are starting
Position offset.
The data extracting unit 204 be used to extract the curved section of the cardiac cycle of interception multiple candidate maximum and
Multiple candidate's minimums.
The wave crest of a cardiac cycle or trough is caused to there is recess due to influencing noise or external environment,
So as to cause the maxima and minima for the mistake extracted, it is therefore desirable to be rejected to the maxima and minima of mistake.
The minimum computing unit 205 is used to fall into a trap according to the multiple candidate maximum and multiple candidate minimums
Calculate true minimum.
Minimum computing unit 305 establishes subelement including reference threshold computation subunit, judgment sub-unit, candidate list
And extreme value demarcates subelement.
The reference threshold computation subunit is used to calculate according to formula TH1=maxVal- (maxVal-minVal)/K1
Maximum compares threshold value out, calculates minimum according to formula TH2=maxVal+ (maxVal-minVal)/K2 and compares threshold value,
Wherein, maxVal is current power maximum value, and minVal is current power minimum value, and K1, K2 are respectively presetting empirical value,
It is tested according to inventor, K1, the value of K2 is more than or equal to 3 preferably, for example, K1 is 4, k2 3, TH1 is that maximum compares threshold value,
TH2 is that minimum compares threshold value.
Judgment sub-unit, for judging it is very big whether N number of performance number of each candidate maximum two sides is respectively less than the candidate
Whether value and candidate's maximum, which are greater than the maximum, compares threshold value, and whether N number of performance number of each candidate's maximum two sides
It is all larger than candidate's minimum, and whether candidate's minimum is less than the minimum and compares threshold value.
For example, N is integer, and the value of N is bigger in the present embodiment, it is higher to the rejecting dynamics of the maximum value of mistake, but
It is that N value is bigger, the efficiency of screening is slower, therefore in the present embodiment, and the value of N is preferable 6~12.
Candidate list establishes subelement, if N number of performance number for one of candidate maximum two sides is respectively less than the time
When selecting maximum and candidate's maximum to be greater than the maximum to compare threshold value, it is candidate which is added to maximum
In list, if N number of performance number of one of candidate's minimum two sides is all larger than candidate's minimum and candidate's minimum is small
When the minimum compares threshold value, which is added in minimum candidate list；
Extreme value demarcates subelement, is demarcated as true pole for finding a maximum maximum from maximum candidate list
Big value, calculates separately the slope of line between each minimum and true maximum value in minimum candidate list, most by slope value
The corresponding minimum of small line is demarcated as true minimum.
If screening candidate maximum is F ' (P), judge whether F ' (P) meets formulaCondition, if wherein
One candidate maximum meets, then it is assumed that and the F ' (P) is true maximum value, similarly,
Judge whether F ' (P) meets formulaIf one of wait
Minimum is selected to meet, then it is assumed that the F ' (P) is true minimum.
The result output unit 206 is used to export the true minimum.
As shown in figure 5, the true minimum of output can be shown in the display screen of an intelligent terminal, in order to make doctor
More can intuitively observe period and the period starting point, end cycle point of ultrasonic power spectrum, can by the initial position in period and
The terminal position in period is demarcated.According to cardiac cycle T and wave trough position P', the first side of a newest cardiac cycle is obtained
Boundary location point T1 and the second boundary location point T2, wherein T1 and T2 is trough, and the distance between two adjacent troughs are
Cardiac cycle, to more accurately demarcate the size in period.
Referring to Fig. 6, the embodiment of the invention also provides a kind of sides for composing progress heart rate period calibration based on ultrasonic power
Method, it should be noted that the method for carrying out the calibration of heart rate period is composed provided by the present embodiment based on ultrasonic power, it is substantially former
Reason and the technical effect generated are identical with above-described embodiment, and to briefly describe, the present embodiment part does not refer to place, can refer to
Corresponding contents in the embodiment stated.It is described based on ultrasonic power compose carry out the calibration of heart rate period method include:
Step S601: Doppler's envelope curve of the ultrasonic power spectrum of input is identified.
It is to be appreciated that step S601 can be performed by envelope recognition unit 201.
In the present embodiment, step S601 be specifically as follows according to percentage algorithm or geometric algorithm or improve geometric algorithm or
Thresholding algorithm or Wavelet Transformation Algorithm or dynamic programming algorithm identify Doppler's envelope curve of the ultrasonic power spectrum of input.
Step S602: cardiac cycle is calculated according to Doppler's envelope curve.
It is to be appreciated that step S602 can be performed by period computing unit 202.
In the present embodiment, it is described according to Doppler's envelope curve calculate cardiac cycle include but is not limited to following two
Mode:
The first: calculating the auto-correlation function of Doppler's envelope curve, then calculates the accumulative variance of auto-correlation function
Function, first maximum point corresponding time point using accumulative variance function apart from zero point is as the period.
Second: the auto-correlation function of Doppler's envelope curve is calculated, then calculates the gradient function of auto-correlation function,
First maximum point corresponding time point using gradient function apart from zero point is as the period.
Step S603: the curved section of a cardiac cycle of Doppler's envelope curve is intercepted.
It is to be appreciated that step S603 can be performed by interception unit 203.
Step S604: multiple candidate maximum and the multiple candidates for extracting the curved section of the cardiac cycle of interception are minimum
Value.
It is to be appreciated that step S604 can be performed by data extracting unit 204.
Step S605: according to calculating true minimum in the multiple candidate maximum and multiple candidate minimums.
It is to be appreciated that step S605 can be performed by minimum computing unit 205.
The specific embodiment of step S605 can be with are as follows:
Maximum is calculated according to formula TH1=maxVal- (maxVal-minVal)/K1 and compares threshold value, according to formula
TH2=maxVal+ (maxVal-minVal)/K2 calculates minimum and compares threshold value, wherein maxVal is that current power is maximum
Value, minVal are current power minimum value, and K1, K2 are respectively presetting empirical value, and TH1 is that maximum compares threshold value, and TH2 is
Minimum compares threshold value.
Judge whether N number of performance number of each candidate maximum two sides is respectively less than candidate's maximum and candidate's maximum
Whether it is greater than the maximum and compares threshold value, it is minimum whether N number of performance number of each candidate's maximum two sides is all larger than the candidate
Value, and whether candidate's minimum is less than the minimum and compares threshold value.
If N number of performance number of one of candidate's maximum two sides is respectively less than candidate's maximum and candidate's maximum is big
When the maximum compares threshold value, which is added in maximum candidate list, if one of candidate pole
When N number of performance number of small value two sides is all larger than candidate's minimum and candidate's minimum is less than minimum comparison threshold value,
Candidate's minimum is added in minimum candidate list；
A maximum maximum is found from maximum candidate list and is demarcated as true maximum value, calculates separately minimum
In candidate list between each minimum and true maximum value line slope, by the corresponding minimum of the smallest line of slope value
It is demarcated as true minimum.
Step S606: the true minimum is exported.
It is to be appreciated that step S606 can be performed by result output unit 206.
In conclusion apparatus and method provided by the invention that carry out the calibration of heart rate period of being composed based on ultrasonic power are by cutting
The curved section for taking a cardiac cycle of Doppler's envelope curve, extracts the multiple of the curved section of the cardiac cycle of interception
It is calculated in candidate maximum and multiple candidate minimum foundations, the multiple candidate's maximum and multiple candidate minimums true
Minimum, to eliminate the multiple candidate maximum and time of the generation of the jitter as caused by noise or outside cause
Minimum is selected, improves the accuracy of period calculating, and can be according to true maximum value and true minimum in a cycle
Starting point, the end point of Doppler's envelope curve demarcate so that doctor is more straight to the observation of Doppler's envelope curve
It sees, so that the conclusion obtained is more efficient.
In several embodiments provided herein, it should be understood that disclosed device and method can also pass through
Other modes are realized.The apparatus embodiments described above are merely exemplary, for example, flow chart and block diagram in attached drawing
Show the device of multiple embodiments according to the present invention, the architectural framework in the cards of method and computer program product,
Function and operation.In this regard, each box in flowchart or block diagram can represent the one of a module, section or code
Part, a part of the module, section or code, which includes that one or more is for implementing the specified logical function, to be held
Row instruction.It should also be noted that function marked in the box can also be to be different from some implementations as replacement
The sequence marked in attached drawing occurs.For example, two continuous boxes can actually be basically executed in parallel, they are sometimes
It can execute in the opposite order, this depends on the function involved.It is also noted that every in block diagram and or flow chart
The combination of box in a box and block diagram and or flow chart can use the dedicated base for executing defined function or movement
It realizes, or can realize using a combination of dedicated hardware and computer instructions in the system of hardware.
In addition, each functional module in each embodiment of the present invention can integrate one independent portion of formation together
Point, it is also possible to modules individualism, an independent part can also be integrated to form with two or more modules.
It, can be with if the function is realized and when sold or used as an independent product in the form of software function module
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.It needs
Illustrate, herein, relational terms such as first and second and the like be used merely to by an entity or operation with
Another entity or operation distinguish, and without necessarily requiring or implying between these entities or operation, there are any this realities
The relationship or sequence on border.Moreover, the terms "include", "comprise" or its any other variant are intended to the packet of nonexcludability
Contain, so that the process, method, article or equipment for including a series of elements not only includes those elements, but also including
Other elements that are not explicitly listed, or further include for elements inherent to such a process, method, article, or device.
In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including the element
Process, method, article or equipment in there is also other identical elements.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any to repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.It should also be noted that similar label and letter exist
Similar terms are indicated in following attached drawing, therefore, once being defined in a certain Xiang Yi attached drawing, are then not required in subsequent attached drawing
It is further defined and explained.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
It should be noted that, in this document, relational terms such as first and second and the like are used merely to a reality
Body or operation are distinguished with another entity or operation, are deposited without necessarily requiring or implying between these entities or operation
In any actual relationship or order or sequence.Moreover, the terms "include", "comprise" or its any other variant are intended to
Non-exclusive inclusion, so that the process, method, article or equipment including a series of elements is not only wanted including those
Element, but also including other elements that are not explicitly listed, or further include for this process, method, article or equipment
Intrinsic element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that
There is also other identical elements in process, method, article or equipment including the element.
Claims (8)
1. a kind of compose the device for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described to be composed based on ultrasonic power
Carry out the calibration of heart rate period device include:
Envelope recognition unit goes out Doppler's envelope curve of the ultrasonic power spectrum of input for identification；
Period computing unit, for calculating cardiac cycle according to Doppler's envelope curve；
Interception unit, the curved section of the cardiac cycle for intercepting Doppler's envelope curve；
Data extracting unit, multiple candidate maximum of the curved section of the cardiac cycle for extracting interception and multiple candidate poles
Small value；
Minimum computing unit, for true minimum according to being calculated in the multiple candidate maximum and multiple candidate minimums
Value；The minimum computing unit includes:
Reference threshold computation subunit, for calculating maximum according to formula TH1=maxVal- (maxVal-minVal)/K1
Threshold value is compared, minimum is calculated according to formula TH2=maxVal+ (maxVal-minVal)/K2 and compares threshold value, wherein
MaxVal is current power maximum value, and minVal is current power minimum value, and K1, K2 are respectively presetting empirical value, and TH1 is
Maximum compares threshold value, and TH2 is that minimum compares threshold value；
Judgment sub-unit, judges whether N number of performance number of each candidate maximum two sides is respectively less than candidate's maximum and the time
It selects maximum whether to be greater than the maximum and compares threshold value, whether N number of performance number of each candidate's minimum two sides is all larger than this
Candidate minimum, and whether candidate's minimum is less than the minimum and compares threshold value；
Candidate list establishes unit, if N number of performance number for one of candidate maximum two sides is respectively less than the candidate greatly
When value and candidate's maximum are greater than maximum comparison threshold value, which is added to maximum candidate list
In, if N number of performance number of one of candidate's minimum two sides is all larger than candidate's minimum and candidate's minimum is less than institute
When stating minimum comparison threshold value, which is added in minimum candidate list；
Extreme value demarcates unit, is demarcated as true maximum value for finding a maximum maximum from maximum candidate list,
The slope for calculating separately line between each minimum and true maximum value in minimum candidate list, by the smallest company of slope value
The corresponding minimum of line is demarcated as true minimum；
As a result output unit, for exporting the true minimum.
2. according to claim 1 compose the device for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described
Period computing unit is used to calculate the auto-correlation function of Doppler's envelope curve, then calculates the accumulative variance of auto-correlation function
Function, first maximum point corresponding time point using accumulative variance function apart from zero point is as the period.
3. according to claim 1 compose the device for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described
Period computing unit is used to calculate the auto-correlation function of Doppler's envelope curve, then calculates the gradient letter of auto-correlation function
Number, first maximum point corresponding time point using gradient function apart from zero point is as the period.
4. according to claim 1 compose the device for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described
Envelope recognition unit is used for according to percentage algorithm or geometric algorithm or improves geometric algorithm or thresholding algorithm or wavelet transformation calculation
Method or dynamic programming algorithm identify Doppler's envelope curve of the ultrasonic power spectrum of input.
5. a kind of compose the method for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described to be composed based on ultrasonic power
Carry out the calibration of heart rate period method include:
Identify Doppler's envelope curve of the ultrasonic power spectrum of input；
Cardiac cycle is calculated according to Doppler's envelope curve；
Intercept the curved section of a cardiac cycle of Doppler's envelope curve；
Extract the multiple candidate maximum and multiple candidate minimums of the curved section of the cardiac cycle of interception；
True minimum is calculated according to the multiple candidate maximum and multiple candidate minimums；It is described according to the multiple time
The step of selecting maximum and multiple candidate minimums to calculate true minimum include:
Maximum is calculated according to formula TH1=maxVal- (maxVal-minVal)/K1 and compares threshold value, according to formula TH2=
MaxVal+ (maxVal-minVal)/K2 calculates minimum and compares threshold value, wherein and maxVal is current power maximum value,
MinVal is current power minimum value, and K1, K2 are respectively presetting empirical value, and TH1 is that maximum compares threshold value, and TH2 is pole
Small value compares threshold value；
Judge whether N number of performance number of each candidate maximum two sides be respectively less than candidate's maximum and candidate's maximum
Threshold value is compared greater than the maximum, whether N number of performance number of each candidate's minimum two sides is all larger than candidate's minimum, and
Whether candidate's minimum, which is less than the minimum, compares threshold value；
If N number of performance number of one of candidate's maximum two sides is respectively less than candidate's maximum and candidate's maximum is greater than institute
When stating maximum comparison threshold value, which is added in maximum candidate list, if one of candidate's minimum
When N number of performance number of two sides is all larger than candidate's minimum and candidate's minimum is less than minimum comparison threshold value, by this
Candidate minimum is added in minimum candidate list；
A maximum maximum is found from maximum candidate list and is demarcated as true maximum value, calculates separately minimum candidate
In list between each minimum and true maximum value line slope, the corresponding minimum of the smallest line of slope value is demarcated
For true minimum；
The true minimum is exported.
6. according to claim 5 compose the method for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described
Include: according to the step of Doppler's envelope curve calculating cardiac cycle
The auto-correlation function of Doppler's envelope curve is calculated, then calculates the accumulative variance function of auto-correlation function, will be added up
First maximum point corresponding time point of the variance function apart from zero point is as the period.
7. according to claim 5 compose the method for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described
The step of calculating cardiac cycle according to Doppler's envelope curve, comprising: calculate the auto-correlation of Doppler's envelope curve
Function, then the gradient function of auto-correlation function is calculated, first maximum point corresponding time by gradient function apart from zero point
Point is used as the period.
8. according to claim 5 compose the method for carrying out the calibration of heart rate period based on ultrasonic power, which is characterized in that described
Identify input ultrasonic power spectrum Doppler's envelope curve the step of include:
According to percentage algorithm or geometric algorithm or improve geometric algorithm or thresholding algorithm or Wavelet Transformation Algorithm or Dynamic Programming
Algorithm identifies Doppler's envelope curve of the ultrasonic power spectrum of input.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710159512.6A CN106821420B (en) | 2017-03-17 | 2017-03-17 | The apparatus and method for carrying out the calibration of heart rate period are composed based on ultrasonic power |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710159512.6A CN106821420B (en) | 2017-03-17 | 2017-03-17 | The apparatus and method for carrying out the calibration of heart rate period are composed based on ultrasonic power |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106821420A CN106821420A (en) | 2017-06-13 |
| CN106821420B true CN106821420B (en) | 2019-06-11 |
Family
ID=59145160
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710159512.6A Active CN106821420B (en) | 2017-03-17 | 2017-03-17 | The apparatus and method for carrying out the calibration of heart rate period are composed based on ultrasonic power |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106821420B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2023184951A1 (en) * | 2022-03-31 | 2023-10-05 | 武汉联影医疗科技有限公司 | Method and system for processing ultrasonic image |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1732854A (en) * | 2004-03-24 | 2006-02-15 | 株式会社东芝 | Ultrasonograph |
| CN103961141A (en) * | 2013-02-02 | 2014-08-06 | 中国人民解放军第四军医大学 | Imaged based backtracking intravascular ultrasound (IVUS) image heart-beating gating method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9724067B2 (en) * | 2009-10-27 | 2017-08-08 | Echosense Jersey Limited | Transthoracic pulmonary doppler ultrasound for evaluating the heart or lung via doppler shift power spectrum |
-
2017
- 2017-03-17 CN CN201710159512.6A patent/CN106821420B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1732854A (en) * | 2004-03-24 | 2006-02-15 | 株式会社东芝 | Ultrasonograph |
| CN103961141A (en) * | 2013-02-02 | 2014-08-06 | 中国人民解放军第四军医大学 | Imaged based backtracking intravascular ultrasound (IVUS) image heart-beating gating method |
Non-Patent Citations (1)
| Title |
|---|
| 超声多普勒子系统设计及计算机仿真;张琳;《四川大学硕士学位论文》;20070807;第21-22、30、34-36、51-52页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106821420A (en) | 2017-06-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10905345B2 (en) | Non-invasive method and system for characterizing cardiovascular systems | |
| US10980429B2 (en) | Method and system for cuffless blood pressure estimation using photoplethysmogram features and pulse transit time | |
| Gupta et al. | R-peak detection based chaos analysis of ECG signal | |
| US9968265B2 (en) | Method and system for characterizing cardiovascular systems from single channel data | |
| Khafaga et al. | Optimization of Electrocardiogram Classification Using Dipper Throated Algorithm and Differential Evolution. | |
| Udhayakumar et al. | Approximate entropy profile: a novel approach to comprehend irregularity of short-term HRV signal | |
| EP2578163A1 (en) | Malignant tissue tumor detection method and malignant tissue tumor detection device | |
| Rohila et al. | Phase entropy: A new complexity measure for heart rate variability | |
| EP3830837A1 (en) | Systems and methods for tissue assessment | |
| US20230309846A1 (en) | System and methods for micro impulse radar detection of physiological information | |
| Padmavathi et al. | Myocardial infarction detection using magnitude squared coherence and support vector machine | |
| CN106821420B (en) | The apparatus and method for carrying out the calibration of heart rate period are composed based on ultrasonic power | |
| CN116421163A (en) | Vital sign detection method and device | |
| CN109949897A (en) | Display methods, device and the terminal device of monitoring curve | |
| Keshtkar et al. | Discriminant analysis between myocardial infarction patients and healthy subjects using Wavelet Transformed signal averaged electrocardiogram and probabilistic neural network | |
| CN107890347B (en) | Filter out the control method and system of power frequency interference signals | |
| JP7468201B2 (en) | Apparatus, method and system for detecting biological activity | |
| Hayajneh et al. | Channel state information based device free wireless sensing for IoT devices employing TinyML | |
| WO2019096031A1 (en) | Systems and methods for multi-resolution discriminant analysis for ultrasound imaging | |
| Sikkandar et al. | Analysis of cardiac abnormalities using Hilbert-Huang transform | |
| CN112598033B (en) | Physiological signal processing method, device, equipment and storage medium | |
| Gudigar et al. | Directional-guided motion sensitive descriptor for automated detection of hypertension using ultrasound images | |
| US11064949B2 (en) | Method and apparatus to remove noise from electrocardiography (ECG) sensor signal | |
| US20230389850A1 (en) | Noninvasive infarct size determination | |
| Senyukova et al. | Universal multi-complexity measures for physiological state quantification in intelligent diagnostics and monitoring systems |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |