CN112494020B - Method for acquiring aortic pressure curve of interest and storage medium - Google Patents
Method for acquiring aortic pressure curve of interest and storage medium Download PDFInfo
- Publication number
- CN112494020B CN112494020B CN202011348579.2A CN202011348579A CN112494020B CN 112494020 B CN112494020 B CN 112494020B CN 202011348579 A CN202011348579 A CN 202011348579A CN 112494020 B CN112494020 B CN 112494020B
- Authority
- CN
- China
- Prior art keywords
- pressure
- aortic
- acquiring
- list
- curve
- 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
- 238000000034 method Methods 0.000 title claims abstract description 77
- 230000035488 systolic blood pressure Effects 0.000 claims description 21
- 230000035487 diastolic blood pressure Effects 0.000 claims description 20
- 238000004590 computer program Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000004321 preservation Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 5
- 239000003086 colorant Substances 0.000 claims description 2
- 210000000709 aorta Anatomy 0.000 claims 1
- 238000004458 analytical method Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 210000004351 coronary vessel Anatomy 0.000 description 7
- 230000006870 function Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 238000012545 processing Methods 0.000 description 4
- 208000031481 Pathologic Constriction Diseases 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000036262 stenosis Effects 0.000 description 3
- 208000037804 stenosis Diseases 0.000 description 3
- 201000000057 Coronary Stenosis Diseases 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000002586 coronary angiography Methods 0.000 description 2
- 208000029078 coronary artery disease Diseases 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 210000004165 myocardium Anatomy 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 206010002383 Angina Pectoris Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000036772 blood pressure Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 230000036770 blood supply Effects 0.000 description 1
- 150000001720 carbohydrates Chemical class 0.000 description 1
- 235000014633 carbohydrates Nutrition 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000000544 hyperemic effect Effects 0.000 description 1
- 208000028867 ischemia Diseases 0.000 description 1
- 150000002632 lipids Chemical group 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000002107 myocardial effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 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/021—Measuring pressure in heart or blood vessels
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
Abstract
The present application provides an acquisition of an aortic pressure curve of interestComprising: acquisition of aortic pressure P a Data; according to the aortic pressure P a Extracting a steady pressure waveform; generating an aortic pressure waveform image from the stationary pressure waveform; acquiring a time of interest; and selecting a corresponding aortic pressure curve from the aortic pressure waveform image according to the interesting moment. The method can obtain the interesting aortic pressure curve part from the continuous aortic pressure curve, and in this way, the method can avoid the influence of objective environment, and overcome the defect that the clinical analysis of the aortic pressure curve can only depend on observation and experience.
Description
Technical Field
The invention relates to the technical field of coronary artery medicine, in particular to a method for acquiring an interesting aortic pressure curve and a storage medium.
Background
The deposition of lipids and carbohydrates in human blood on the vessel wall will form plaque on the vessel wall, which in turn leads to stenosis of the vessel; especially, the stenosis of blood vessels around the coronary artery will lead to myocardial blood supply deficiency, induce coronary heart disease, angina pectoris and other diseases, and cause serious threat to human health. According to statistics, the number of patients with the existing coronary heart disease in China is about 1100 ten thousand, and the number of patients with cardiovascular interventional operation treatment is increased by more than 10% each year.
Although the conventional medical detection means such as Coronary Angiography (CAG) and Computed Tomography (CT) can show the severity of coronary stenosis of heart, the ischemia of the coronary artery cannot be accurately evaluated. In order to improve the accuracy of coronary blood vessel function evaluation, pijls in 1993 proposed a new index of calculating coronary blood vessel function by pressure measurement, namely fractional flow reserve (Fractional Flow Reserve, FFR), and FFR has become a gold standard for coronary stenosis function evaluation through long-term basic and clinical studies.
Fractional Flow Reserve (FFR) generally refers to the fractional flow reserve of the heart muscle, which is set as the ratio of the maximum blood flow that a diseased coronary can provide to the heart muscle to the maximum blood flow at which the coronary is completely normal, and studies have shown that the ratio of blood flow can be replaced with a pressure value at the maximum hyperemic state of the coronary. That is, the FFR value can be measured and then calculated by measuring the pressure at the distal end stenosis of the coronary artery and the proximal end pressure of the coronary artery by the pressure sensor in the maximum congestion state of the coronary artery.
In the current devices, the aortic pressure curve is obtained after acquisition by a pressure transducer and digital-to-analog conversion. In practical use, due to the influence of objective environment, it is impossible to directly obtain a smooth and stable aortic pressure curve through a sensor, and some unstable factors always exist, so that the pressure curve fluctuates, is distorted or is discontinuous. To obtain a pressure waveform that can be subjected to parameter calculations, the waveform must be processed by a series of means to obtain the aortic pressure curve of interest.
Disclosure of Invention
The invention provides a method and a storage medium for acquiring an interesting aortic pressure curve, which are used for solving the problem that a smooth and stable aortic pressure curve cannot be directly acquired through a sensor in an objective environment.
To achieve the above object, in a first aspect, the present application provides a method of acquiring an aortic pressure curve of interest, comprising:
acquisition of aortic pressure P a Data;
according to the aortic pressure P a Extracting a steady pressure waveform;
generating an aortic pressure waveform image from the stationary pressure waveform;
acquiring a time of interest;
and selecting a corresponding aortic pressure curve from the aortic pressure waveform image according to the interesting moment.
Optionally, the method for acquiring the aortic pressure curve of interest according to the above method for selecting the corresponding aortic pressure curve from the aortic pressure waveform image according to the time of interest includes:
acquiring pressure waveform segments within a number of heart cycles forward on the aortic pressure waveform image from the moment of interest;
acquiring an average value of aortic pressures on the pressure waveform segments in the plurality of heartbeat periods;
and removing the pressure value point with larger deviation according to the average value, and redrawing a smooth aortic pressure curve according to the Bezier curve.
Optionally, the above method for acquiring an aortic pressure curve of interest, wherein the pressure waveform segments in the plurality of heartbeat cycles are pressure waveform segments in 6-10 heartbeat cycles.
Optionally, the method for obtaining the aortic pressure curve of interest, removing the pressure value points with larger deviation according to the average value, and redrawing the smooth aortic pressure curve according to the bezier curve includes:
and removing the pressure value points with larger deviation according to the average value, and reserving the pressure curve with the most stable pressure in one or two heart cycles from the redrawn smooth aortic pressure curve as the aortic pressure curve of interest.
Optionally, a method of acquiring an aortic pressure curve of interest as described above, said method comprising the step of determining the aortic pressure P a The method for extracting the steady pressure waveform comprises the following steps:
setting a plurality of index lists, including: temporary storage list, preservation list, average value point index list, maximum value point index list and minimum value point index list;
setting a physiological parameter threshold, comprising: the heart cycle number C, the range from the diastolic pressure to the systolic pressure, the minimum pressure amplitude threshold, the fluctuation difference of the systolic pressure or the diastolic pressure, the heart rate range, the pressure number N threshold, the upper limit of the file number, the cutoff pressure threshold P cutoff, the cutoff heart rate threshold HR cutoff and the pressure transmission rate threshold V threshold;
calculating average pressure according to the index list and the physiological parameter threshold value and the aortic pressure data;
acquiring a pressure average value point according to the average pressure;
dividing the data period according to the pressure average value points, and obtaining pressure parameters in each heartbeat period;
and acquiring a stable pressure waveform according to the pressure parameter.
Optionally, a method for acquiring an aortic pressure curve of interest as described above, the method for calculating an average pressure from the aortic pressure data according to the index list and the physiological parameter threshold value includes:
obtaining a heart cycle time threshold according to the heart rate threshold;
acquiring a data point number range in a heartbeat period according to the heartbeat period number C and a pressure transmission rate threshold;
defining a time period according to the data point number range, counting the number of all data points after finishing one time period, and storing all acquired aortic pressures and acquired time in the temporary storage list in a one-to-one correspondence manner;
if the number of data in the temporary storage list is integer times of the N threshold and is larger than the minimum threshold value of the data point number range and smaller than the maximum threshold value of the data point number range, calculating the average value of the aortic pressure in all temporary storage lists, namely the average pressure.
Optionally, in the above method for acquiring an aortic pressure curve of interest, the method for calculating an average pressure according to the index list and the physiological parameter threshold value by using the aortic pressure data further includes:
if n aortic pressures P are consecutive in the temporary list a All smaller than the cutoff pressure threshold P, removing the n aortic pressures P from the temporary storage list a 。
Optionally, a method for acquiring an aortic pressure curve of interest as described above, the method for acquiring a pressure mean point according to the average pressure includes:
setting a mean point index list;
adding pressure average points into the average point index list;
generating a pressure average point curve according to the pressure average point;
filtering the pressure average value points on the pressure average value point curve;
judging whether the number of the pressure average points in the average point index list meets the heartbeat cycle requirement, and if not, adjusting the size of n until the heartbeat cycle requirement is met.
Optionally, the method for acquiring the aortic pressure curve of interest, the method for adding the pressure mean point into the mean point index list includes:
aortic pressure P in the temporary list a Respectively with the average pressureComparing;
if it isAnd->Will P K Adding the average value point index list into the average value point index list, wherein P K 、P K+1 Respectively representing the Kth and Kth+1st aortic pressures in the temporary storage list;
if it isAnd->Will P K And adding the average value point index list.
Optionally, the method for acquiring the aortic pressure curve of interest, and the method for filtering the pressure average point on the pressure average point curve includes:
according to D n ＝P n+1 -P n Wherein P is n 、P n+1 Respectively represent the nth and the (n+1) th aortic pressures in the mean point index list, and the D n Representing P n+1 、P n Is a pressure difference of (2);
according to HR n ＝30/(D n N) acquiring a heart rate value HR within the nth half-beat cycle n ；
If HR is n < HR cut, then delete the HR from the mean point index list n Corresponding mean pressure points;
if it isAnd at P n The posterior continuous m aortic pressures are all greater than +.>And->Then delete P from the mean point index list n And m pressure average points are continued.
Optionally, in the above method for acquiring an aortic pressure curve of interest, the determining whether the number of pressure mean points in the mean point index list meets the heartbeat cycle requirement, if not, repeating the claim 6, and adjusting the size of n until the heartbeat cycle requirement is met includes:
if the number N of the data in the average value point list is less than 2C+1, the heart cycle requirement is not met, and the size of N is adjusted until N is more than or equal to 2C+1.
Optionally, the method for acquiring the aortic pressure curve of interest, which divides the data period according to the pressure mean point, includes:
if it isExtracting P from the scratch list n Maximum pressure P in the heart cycle max And the maximum pressure P max Position S where it is located max The S is carried out max 、P max Adding the same into a diastolic blood pressure list in a one-to-one correspondence manner;
acquiring P in the diastolic blood pressure list max Obtaining a diastolic blood pressure P comfort;
if it isExtracting P from the scratch list n Minimum pressure P in the heart cycle min And the minimum pressure P min Position S where it is located min The S is carried out min 、P min Adding the data into a systolic pressure list in a one-to-one correspondence manner;
acquiring P in the systolic blood pressure list min Is used to obtain the systolic blood pressure pdiff.
Optionally, a method for acquiring an aortic pressure curve of interest according to the above method for acquiring a stationary pressure waveform according to the pressure parameter includes:
setting the offset as V threshold/a, and intercepting a starting point O and an ending point E;
and executing the mean value point index list, the diastolic pressure list and the systolic pressure list according to the offset, the starting point O and the ending point E, and outputting the steady pressure waveform, the diastolic pressure waveform and the systolic pressure waveform.
Optionally, a method for acquiring an aortic pressure curve of interest as described above, the method for generating an aortic pressure waveform image includes:
setting a display interval of the steady pressure waveform;
setting a data range of the stationary pressure waveform;
and generating an aortic pressure waveform image according to the display interval, the data range and the stable pressure waveform.
Optionally, the method for acquiring the aortic pressure curve of interest, the method for setting the display interval of the stationary pressure waveform includes:
setting the length W of the generated aortic pressure waveform image 1 And width H 1 ；
Setting the outer frame length W of the aortic pressure waveform image 2 And width H 2 ；
According to w=w 1 -2W 2 ，H＝H 1 -2H 2 The length W and width H of the display section of the stationary pressure waveform are obtained.
Optionally, the method for acquiring the aortic pressure curve of interest, the method for setting the data range of the stationary pressure waveform includes:
setting an X-axis duty cycle threshold of the stationary pressure waveform, and an interpolation number between 2 intersections of the stationary pressure waveform and the X-axis; setting sampling interval time; obtaining an X-axis data range of which the X-axis duty ratio is greater than or equal to the X-axis duty ratio threshold value of the steady pressure waveform by adjusting the interpolation number and the sampling interval time; and/or
Setting a Y-axis duty ratio threshold as r, and setting the proportion of the data range to the display space as a; the Y-axis data length is l= (pex-P-receive) + (pex-P-receive)/r; the Y-axis data minimum value Y min P receive- (P comfort-P receive)/ar; the Y-axis data maximum value Y max =pex+ (pex-pex)/ar; the Y-axis data range of the steady pressure is [ Y ] min ，y max ]。
Optionally, the method for acquiring the aortic pressure curve of interest, the method for generating the aortic pressure waveform image according to the display interval, the data range and the stationary pressure waveform includes:
newly building a picture space with length W and width H;
setting background colors or background pictures of the picture space, and adding frames;
setting the color of the (x, y) pixel points;
the X-axis display pixel range is [ W ] 2 ,W 1 -1-W 2 ]The Y-axis display pixel range is [ H ] 2 ,H 1 -1-H 2 ]；
And generating an aortic pressure waveform image according to the (x, y) pixel points.
In a second aspect, the present application provides a computer storage medium, which when executed by a processor, implements the above-described method of acquiring an aortic pressure curve of interest.
The beneficial effects brought by the scheme provided by the embodiment of the application at least comprise:
according to the method and the device, the aortic pressure waveform image required by the user is finally obtained through calculation and judgment of the original aortic pressure waveform. The invention can obtain the interesting aortic pressure curve, and overcomes the defect that the clinical analysis of the aortic pressure curve only depends on observation and experience.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the invention and do not constitute a limitation on the invention. In the drawings:
FIG. 1 is a flow chart of a method of acquiring an aortic pressure curve of interest of the present application;
FIG. 2 is a flow chart of the present application S200;
fig. 3 is a flowchart of the present application S220;
fig. 4 is a flowchart of the present application S230;
fig. 5 is a flowchart of the present application S232;
fig. 6 is a flowchart of the present application S234;
FIG. 7 is a flow chart of the present application S240;
FIG. 8 is a flow chart of the present application S250;
fig. 9 is a flowchart of the present application S300;
fig. 10 is a flowchart of the present application S310;
FIG. 11 is a flow chart of the present application S320;
FIG. 12 is a flow chart of the present application S330;
fig. 13 is a flowchart of the present application S500;
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be clearly and completely described below with reference to specific embodiments of the present invention and corresponding drawings. It will be apparent that the described embodiments are only some, but not all, embodiments of the invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Various embodiments of the invention are disclosed in the following drawings, in which details of the practice are set forth in the following description for the purpose of clarity. However, it should be understood that these practical details are not to be taken as limiting the invention. That is, in some embodiments of the invention, these practical details are unnecessary. Moreover, for the purpose of simplifying the drawings, some conventional structures and components are shown in the drawings in a simplified schematic manner.
As shown in fig. 1, to solve the above problem, the present application provides a method for acquiring an aortic pressure curve of interest, which includes:
the present application provides a method of acquiring an aortic pressure curve of interest, comprising:
s100, obtaining aortic pressure P a Data, comprising:
aortic inlet pressure is obtained by a disposable blood pressure sensor.
S200, as shown in FIG. 2, according to aortic pressure P a Extracting a plateau pressure waveform, comprising:
s210, setting a plurality of index lists and physiological parameter thresholds;
the plurality of index lists includes: temporary storage list, save list, mean value point index list, maximum value point index list and minimum value point index list.
The physiological parameter threshold includes: the heart cycle number C, the range from diastolic to systolic pressure, the minimum pressure amplitude threshold, the fluctuation difference of systolic or diastolic pressure, the heart rate range, the pressure number N threshold, the upper limit of the file number, the cutoff pressure threshold P cutoff, the cutoff heart rate threshold HR cutoff, and the pressure transmission rate threshold V threshold.
S220, as shown in fig. 3, calculating an average pressure from the index list, the physiological parameter threshold, and the aortic pressure data, including:
s221, obtaining a heartbeat cycle time threshold according to the heart rate threshold;
s222, acquiring the data point number range in a heartbeat cycle according to the heartbeat cycle number C and the pressure transmission rate threshold;
s223, defining a time period according to the number range of the data points, counting the number of all the data points after finishing one time period, and storing all the acquired aortic pressures and the acquired time in a temporary storage list in a one-to-one correspondence manner;
if n continuous aortic pressures P are in the temporary list a All of which are less than the cutoff pressure threshold P cut, removing the n aortic pressures P from the temporary list a The calculated amount is reduced, and the running speed is improved.
S224, if the number of data in the temporary storage list is integer times of the N threshold and is larger than the minimum threshold value of the number range of data points and smaller than the maximum threshold value of the number range of data points, calculating the average value of the aortic pressure in all temporary storage lists, namely the average pressure.
S230, as shown in FIG. 4, acquiring a pressure average point according to the average pressure includes:
s231, setting a mean point index list;
s232, as shown in FIG. 5, adding pressure average points into the average point index list comprises:
s2321, aortic pressure P in temporary storage list a Respectively with average pressureComparing;
s2322, ifAnd->Will P K Added to the mean point index list, where P K 、P K+1 Respectively representThe Kth and Kth+1st aortic pressures in the temporary storage list;
s2323, ifAnd->Will P K And adding the average value point index list.
S233, generating a pressure average point curve according to the pressure average point;
s234, as shown in fig. 6, filtering the pressure average point on the pressure average point curve, including:
s2341 according to D n ＝P n+1 -P n Wherein P is n 、P n+1 Respectively represent the nth and the (n+1) th aortic pressures in the mean point index list, D n Representing P n+1 、P n Is a pressure difference of (2);
s2342 according to HR n ＝30/(D n N) obtaining heart rate values HR in the nth half-beat cycle n ；
S2343 if HR n < HR cut, then delete HR from the mean point index list n Corresponding mean pressure points;
s2344, ifAnd at P n The posterior continuous m aortic pressures are all greater than +.>And-> Then delete P from the mean point index list n And m pressure average points are continued.
S235, judging whether the number of pressure mean points in the mean point index list meets the heartbeat cycle requirement, if not, repeating the method, and adjusting the size of n until the number meets the heartbeat cycle requirement, wherein the method comprises the following steps:
if the number N of the data in the average value point list is less than 2C+1, the heartbeat cycle requirement is not met, and the method is repeated;
and adjusting the size of N until N is equal to or greater than 2C+1.
S240, as shown in FIG. 7, the data period is divided according to the pressure mean point, and the pressure parameter in each heartbeat period is obtained, including:
s241, ifExtracting P from scratch list n Maximum pressure P in the heart cycle max And maximum pressure P max Position S where it is located max Will S max 、P max Adding the same into a diastolic blood pressure list in a one-to-one correspondence manner;
s242, obtaining P in diastolic blood pressure list max Obtaining a diastolic blood pressure P comfort;
if the P comfort is not in the range from the diastolic pressure to the systolic pressure, the temporary storage list is acquired again, wherein the temporary storage list comprises the following steps: intercepting the first half part of the temporary storage list into the preservation list, and if the number of the data is larger than the upper limit of the number of the preservation list, directly replacing the first half part of the data of the temporary storage list adopting the preservation list into the temporary storage list, and clearing the preservation list.
S243, ifExtracting P from scratch list n Minimum pressure P in the heart cycle min And a minimum pressure P min Position S where it is located min Will S min 、P min Adding the data into a systolic pressure list in a one-to-one correspondence manner;
s244, obtaining P in the systolic blood pressure list min Is used to obtain the systolic blood pressure pdiff.
S250, as shown in fig. 8, acquiring a stationary pressure waveform according to the pressure parameter, including:
s251, setting offset as V threshold/a, and intercepting a starting point O and an ending point E;
s252, executing the mean point index list, the diastolic pressure list and the systolic pressure list according to the offset, the starting point O and the ending point E, and outputting a steady pressure waveform, a diastolic pressure waveform and a systolic pressure waveform.
S300, as shown in fig. 9, generating an aortic pressure waveform image from the plateau pressure waveform, comprising:
s310, as shown in fig. 10, of setting a display section of a steady pressure waveform, including:
s311, setting the length W of the generated aortic pressure waveform image 1 And width H 1 ；
S312, setting the outer frame length W of the aortic pressure waveform image 2 And width H 2 ；
S313, according to w=w 1 -2W 2 ，H＝H 1 -2H 2 The length W and width H of the display section of the stationary pressure waveform are obtained.
S320, as shown in fig. 11, sets a data range of the stationary pressure waveform, including:
s321, setting an X-axis duty ratio threshold of the steady pressure waveform and interpolation numbers between 2 crossing points of the steady pressure waveform and the X axis; setting sampling interval time; obtaining the X-axis data range of which the X-axis duty ratio is greater than or equal to the X-axis duty ratio threshold value of the steady pressure waveform by adjusting the interpolation number and the sampling interval time;
s322, setting a Y-axis duty ratio threshold as r and setting the proportion of a data range to a display space as a; the Y-axis data length is l= (pex-P-receive) + (pex-P-receive)/r; minimum Y value of Y-axis data min P receive- (P comfort-P receive)/ar; y-axis data maximum Y max =pex+ (pex-pex)/ar; the Y-axis data range for plateau pressure is [ Y ] min ，y max ]。
S330, as shown in fig. 12, generating an aortic pressure waveform image according to the display section, the data range, and the stationary pressure waveform, including:
s331, newly building a picture space with length W and width H;
s332, setting the background color or background picture of the picture space and adding a frame;
s333, setting the color of the (x, y) pixel point;
s334, X-axis display pixel range is [ W ] 2 ,W 1 -1-W 2 ]The Y-axis display pixel range is [ H ] 2 ,H 1 -1-H 2 ]；
S335, generating an aortic pressure waveform image according to the (x, y) pixel points.
S400, acquiring the interested moment;
s500, as shown in FIG. 13, selecting a corresponding aortic pressure curve from the aortic pressure waveform image according to the moment of interest, including:
s510, starting from the interesting moment, acquiring pressure waveform segments in a plurality of heartbeat cycles forwards on an aortic pressure waveform image;
s520, obtaining an average value of aortic pressure on a pressure waveform segment in a plurality of heartbeat cycles, wherein the pressure waveform segment in 6-10 heartbeat cycles is generally selected for the purpose of comprehensive operation and reduced operation quantity;
s530, removing the pressure value point with larger deviation according to the average value, and redrawing a smooth aortic pressure curve according to the Bezier curve, wherein the method comprises the following steps:
and removing the pressure value points with larger deviation according to the average value, and re-drawing the pressure curve with one or two heart cycles with the most stable reserved pressure on the smooth aortic pressure curve as the aortic pressure curve of interest.
The present application provides a computer storage medium, which when executed by a processor implements the above-described method of acquiring an aortic pressure curve of interest.
Those skilled in the art will appreciate that the various aspects of the present invention may be implemented as a system, method, or computer program product. Accordingly, aspects of the invention may be embodied in the following forms, namely: an entirely hardware embodiment, an entirely software embodiment (including firmware, resident software, micro-code, etc.) or an embodiment combining hardware and software aspects may all generally be referred to herein as a "circuit," module "or" system. Furthermore, in some embodiments, aspects of the invention may also be implemented in the form of a computer program product in one or more computer-readable media having computer-readable program code embodied therein. Implementation of the methods and/or systems of embodiments of the present invention may involve performing or completing selected tasks manually, automatically, or a combination thereof.
For example, hardware for performing selected tasks according to embodiments of the invention could be implemented as a chip or circuit. As software, selected tasks according to embodiments of the invention could be implemented as a plurality of software instructions being executed by a computer using any suitable operating system. In an exemplary embodiment of the invention, one or more tasks according to exemplary embodiments of the method and/or system as herein, such as a computing platform for executing a plurality of instructions, are performed by a data processor. Optionally, the data processor comprises volatile storage for storing instructions and/or data and/or non-volatile storage for storing instructions and/or data, e.g. a magnetic hard disk and/or a removable medium. Optionally, a network connection is also provided. A display and/or a user input device such as a keyboard or mouse are optionally also provided.
Any combination of one or more computer readable may be utilized. The computer readable medium may be a computer readable signal medium or a computer readable storage medium. The computer readable storage medium can be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or a combination of any of the foregoing. More specific examples (a non-exhaustive list) of the computer-readable storage medium would include the following:
an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In this document, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.
The computer readable signal medium may include a propagated data signal with computer readable program code embodied therein, either in baseband or as part of a carrier wave. Such a propagated data signal may take any of a variety of forms, including, but not limited to, electro-magnetic, optical, or any suitable combination of the foregoing. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device.
Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to wireless, wireline, optical fiber cable, RF, etc., or any suitable combination of the foregoing.
For example, computer program code for carrying out operations for aspects of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, smalltalk, C ++ or the like and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of remote computers, the remote computer may be connected to the user computer through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computer (e.g., connected through the internet using an internet service provider).
It will be understood that each block of the flowchart illustrations and/or block diagrams, and combinations of blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, or other programmable data processing apparatus to produce a machine, such that the computer program instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable medium that can direct a computer, other programmable data processing apparatus, or other devices to function in a particular manner, such that the instructions stored in the computer-readable medium produce an article of manufacture including instructions which implement the function/act specified in the flowchart and/or block diagram block or blocks (article of manufacture).
The computer program instructions may also be loaded onto a computer (e.g., a coronary artery analysis system) or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer, other programmable apparatus or other devices to produce a computer implemented process such that the instructions which execute on the computer, other programmable device or other devices provide processes for implementing the functions/acts specified in the flowchart and/or block diagram block or blocks.
The foregoing embodiments of the present invention have been described in some detail by way of illustration of the principles of the invention, and it is to be understood that the invention is not limited to the specific embodiments of the invention but is intended to cover modifications, equivalents, alternatives and modifications within the spirit and principles of the invention.
Claims (16)
1. A method of acquiring an aortic pressure curve of interest, comprising:
acquisition of aortic pressure P a Data;
setting a plurality of index lists, including: temporary storage list, preservation list, average value point index list, maximum value point index list and minimum value point index list;
setting a physiological parameter threshold, comprising:
the heart rate cycle number C, the range from diastolic pressure to systolic pressure, the minimum pressure amplitude threshold value and the heart rate range are used for collecting heart rate data;
number of pressures N Threshold value The number of the corresponding pressures acquired in the time obtained through the heart rate and the heart beat period C is equal to or greater than 4;
the upper limit of the number of files is used for limiting the total amount of the collected pressure number;
wherein the cutoff pressure threshold P Cutting off Cut-off heart rate threshold HR Cutting off Pressure transmission rate threshold V Threshold value For collecting heart rate pressure waveforms, less than a cutoff pressure threshold P Cutting off Cut-off heart rate threshold HR Cutting off Pressure transmission rate threshold V Threshold value Is an ineffective pressure waveform;
according to the index list and the physiological parameter threshold, the aortic pressure P a The method for calculating the average pressure by the data comprises the following steps:
according to the number of heart beat cycles C and the pressure transmission rate threshold V Threshold value Acquiring a data point number range in the heartbeat period;
defining a time period according to the data point number range, counting the number of all data points after finishing one time period, and storing all acquired aortic pressures and acquired time in the temporary storage list in a one-to-one correspondence manner;
if the number of data in the temporary storage list is N Threshold value If the data point number is larger than the minimum threshold value of the data point number range and smaller than the maximum threshold value of the data point number range, calculating the average value of the aortic pressures in all temporary storage lists, namely the average pressure;
acquiring a pressure average value point according to the average pressure;
dividing the data period according to the pressure average value points, and obtaining pressure parameters in each heartbeat period;
acquiring a stable pressure waveform image according to the pressure parameter;
acquiring a time of interest;
and selecting a corresponding aortic pressure curve from the aortic pressure waveform image according to the interesting moment.
2. A method of acquiring an aortic pressure curve of interest according to claim 1, wherein selecting a corresponding aortic pressure curve from the aortic pressure waveform image according to the moment of interest comprises:
acquiring pressure waveform segments within a number of heart cycles forward on the aortic pressure waveform image from the moment of interest;
acquiring an average value of aortic pressures on the pressure waveform segments in the plurality of heartbeat periods;
and removing the pressure value point with the maximum heartbeat cycle deviation according to the average value, and redrawing a smooth aortic pressure curve according to the Bezier curve.
3. A method of deriving an aortic pressure curve of interest as claimed in claim 2 wherein the pressure waveform segments over the number of heart cycles are those over 6-10 heart cycles.
4. A method of deriving an aortic pressure curve of interest as claimed in claim 2 wherein said method of redrawing a smoothed aortic pressure curve from said bezier curve by removing pressure value points with maximum heart cycle deviation from said average value comprises:
and removing the pressure value point with the maximum deviation according to the average value, and reserving the pressure curve with the most stable pressure in one or two heart cycles from the redrawn smooth aortic pressure curve as the aortic pressure curve of interest.
5. The method of acquiring an aortic pressure curve of interest as set forth in claim 1, wherein the method of calculating an average pressure from the index list, the physiological parameter threshold, the aortic pressure data further comprises:
if n aortic pressures P are consecutive in the temporary list a Are all smaller than the cut-off pressure threshold value P Cutting off Removing the n aortic pressures P from the temporary list a 。
6. The method of acquiring an aortic pressure curve of interest as set forth in claim 5, wherein the method of acquiring a pressure mean point from the average pressure comprises:
setting a mean point index list;
adding pressure average points into the average point index list;
generating a pressure average point curve according to the pressure average point;
filtering the pressure average value points on the pressure average value point curve;
judging whether the number of the pressure average points in the average point index list meets the heartbeat cycle requirement, and if not, adjusting the size of n until the heartbeat cycle requirement is met.
7. The method of claim 6, wherein adding pressure mean points to the index list of mean points comprises:
aortic pressure P in the temporary list a Respectively comparing with the average pressure P;
if P K > P and P K+1 < P, then P K Adding the average value point index list into the average value point index list, wherein P K、 P K+1 Respectively representing the Kth and Kth+1st aortic pressures in the temporary storage list;
if P K < P and P K+1 > P, then P K And adding the average value point index list.
8. The method of acquiring an aortic pressure curve of interest as set forth in claim 7, wherein the method of filtering pressure mean points on the pressure mean point curve comprises:
according to D n ＝P n+1 -P n Wherein P is n 、P n+1 Respectively represent the nth and the (n+1) th aortic pressures in the mean point index list, and the D n Representing P n+1 、P n Is a pressure difference of (2);
according to HR n ＝30/(D n N) acquiring a heart rate value HR within the nth half-beat cycle n ；
If HR is n ＜HR Cutting off Deleting the HR from the mean point index list n Corresponding mean pressure points;
if P n > P, and at P n The pressure of the posterior m continuous aorta is greater than P, and P n+m+1 < P, then delete P from the mean point index list n And m pressure average points are continued.
9. The method of claim 8, wherein said determining whether the number of pressure mean points in the mean point index list meets the heart cycle requirement, and if not, repeating said claim 1, and adjusting the size of n until the heart cycle requirement is met comprises:
if the number N of data in the mean point list Are all If the heart rate is less than 2C+1, the heart rate is not satisfied, and the size of N is adjusted until N Are all ≥2C+1。
10. The method of acquiring an aortic pressure curve of interest as claimed in claim 9, wherein the method of acquiring pressure parameters in each heartbeat cycle by dividing data cycles according to the pressure mean points comprises:
if P n > P, extracting P from within the scratch list n In the heart cycleMaximum pressure P of (2) max And the maximum pressure P max Position S where it is located max The S is carried out max 、P max Adding the same into a diastolic blood pressure list in a one-to-one correspondence manner;
acquiring P in the diastolic blood pressure list max Average value of (2) to obtain diastolic pressure P Shu (Chinese character) ；
If P n < P, extracting P from the temporary list n Minimum pressure P in the heart cycle min And the minimum pressure P min Position S where it is located min The S is carried out min 、P min Adding the data into a systolic pressure list in a one-to-one correspondence manner;
acquiring P in the systolic blood pressure list min Average value of (2) to obtain the systolic blood pressure P Collecting and recovering 。
11. The method of acquiring an aortic pressure curve of interest as set forth in claim 10, wherein the method of acquiring a plateau pressure waveform from the pressure parameter comprises:
setting the offset to V Threshold value A, where a=4 or a=5, at V Threshold value A starting point and an ending point in the waveform range are offset, and a starting point O and an ending point E are intercepted;
and executing the mean value point index list, the diastolic pressure list and the systolic pressure list according to the offset, the starting point O and the ending point E, and outputting the steady pressure waveform, the diastolic pressure waveform and the systolic pressure waveform.
12. The method of acquiring an aortic pressure curve of interest as set forth in claim 11, wherein the method of generating an aortic pressure waveform image comprises:
setting a display interval of the steady pressure waveform;
setting a data range of the stationary pressure waveform;
and generating an aortic pressure waveform image according to the display interval, the data range and the stable pressure waveform.
13. The method of acquiring an aortic pressure curve of interest as set forth in claim 12, wherein the method of setting a display interval of the plateau pressure waveform comprises:
setting the length W of the generated aortic pressure waveform image 1 And width H 1 ；
Setting the outer frame length W of the aortic pressure waveform image 2 And width H 2 ；
According to w=w 1 -2W 2 ，H＝H 1 -2H 2 The length W and width H of the display section of the stationary pressure waveform are obtained.
14. The method of acquiring an aortic pressure curve of interest as set forth in claim 13, wherein the method of setting the data range of the plateau pressure waveform comprises:
setting an X-axis duty cycle threshold of the stationary pressure waveform, and an interpolation number between 2 intersections of the stationary pressure waveform and the X-axis; setting sampling interval time; obtaining an X-axis data range of which the X-axis duty ratio is greater than or equal to the X-axis duty ratio threshold value of the steady pressure waveform by adjusting the interpolation number and the sampling interval time; and/or
Setting a Y-axis duty ratio threshold as r and setting the proportion of the data range to the display space as a; the Y-axis data length is l= (P Shu (Chinese character) -P Collecting and recovering )+(P Shu (Chinese character) -P Collecting and recovering ) R; the Y-axis data minimum value Y min ＝P Collecting and recovering -(P Shu (Chinese character) -P Collecting and recovering ) Ar; the Y-axis data maximum value Y max ＝P Shu (Chinese character) +(P Shu (Chinese character) -P Collecting and recovering ) Ar; the Y-axis data range of the steady pressure is [ Y ] min ，y max ]。
15. The method of acquiring an aortic pressure curve of interest as set forth in claim 14, wherein the generating an aortic pressure waveform image from the display interval, the data range, and the plateau pressure waveform comprises:
newly building a picture space with length W and width H;
setting background colors or background pictures of the picture space, and adding frames;
setting the color of the (x, y) pixel points;
the X-axis display pixel range is [ W ] 2 ,W 1 -1-W 2 ]The Y-axis display pixel range is [ H ] 2 ,H 1 -1-H 2 ]；
And generating an aortic pressure waveform image according to the (x, y) pixel points.
16. A computer storage medium, characterized in that a computer program, when being executed by a processor, implements the method of acquiring an aortic pressure curve of interest as claimed in any one of claims 1 to 15.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011348579.2A CN112494020B (en) | 2020-11-26 | 2020-11-26 | Method for acquiring aortic pressure curve of interest and storage medium |
| PCT/CN2020/132275 WO2022110020A1 (en) | 2020-11-26 | 2020-11-27 | Method for obtaining aortic pressure curve of interest, and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011348579.2A CN112494020B (en) | 2020-11-26 | 2020-11-26 | Method for acquiring aortic pressure curve of interest and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112494020A CN112494020A (en) | 2021-03-16 |
| CN112494020B true CN112494020B (en) | 2024-01-23 |
Family
ID=74966271
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011348579.2A Active CN112494020B (en) | 2020-11-26 | 2020-11-26 | Method for acquiring aortic pressure curve of interest and storage medium |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN112494020B (en) |
| WO (1) | WO2022110020A1 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112494021A (en) * | 2020-11-26 | 2021-03-16 | 苏州润迈德医疗科技有限公司 | Method for generating aorta pressure waveform image and storage medium |
| CN112494019A (en) * | 2020-11-26 | 2021-03-16 | 苏州润迈德医疗科技有限公司 | Method for extracting flat stable pressure waveform based on aortic pressure and storage medium |
| CN115271256B (en) * | 2022-09-20 | 2022-12-16 | 华东交通大学 | Intelligent ordering method under multi-dimensional classification |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109475309A (en) * | 2016-07-13 | 2019-03-15 | 悦享趋势科技（北京）有限责任公司 | The method and apparatus for detecting physiological status |
| CN109770888A (en) * | 2019-03-19 | 2019-05-21 | 苏州润迈德医疗科技有限公司 | The method instantaneously without waveform ratio is calculated based on pressure sensor and contrastographic picture |
| CN109805949A (en) * | 2019-03-19 | 2019-05-28 | 苏州润迈德医疗科技有限公司 | The method for calculating blood flow reserve score based on pressure sensor and contrastographic picture |
| CN110226923A (en) * | 2018-03-05 | 2019-09-13 | 苏州润迈德医疗科技有限公司 | A method of blood flow reserve score is measured without vasodilator |
| CN110353645A (en) * | 2019-08-19 | 2019-10-22 | 苏州润迈德医疗科技有限公司 | Measurement is without waveform phase pressure, ratio approach, device, system and storage medium |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9314170B2 (en) * | 2010-05-07 | 2016-04-19 | Atcor Medical Pty Ltd | Brachial cuff |
| CN103284703B (en) * | 2013-06-18 | 2015-05-13 | 中国科学院合肥物质科学研究院 | Aortic pulse wave transfer time measuring method based on upper extremity artery information |
| CN104983412B (en) * | 2015-05-28 | 2017-08-11 | 中国科学院合肥物质科学研究院 | A kind of central aortic systolic phase average normalized blood flow waveform model and the method that aortic pulse wave passing time is obtained based on it |
| CN107992452B (en) * | 2017-12-12 | 2020-12-18 | 北京动亮健康科技有限公司 | Method, device, storage medium and equipment for calculating central hemodynamics index |
| CN110786840B (en) * | 2019-11-04 | 2021-06-08 | 苏州润迈德医疗科技有限公司 | Method, device and storage medium for obtaining blood vessel evaluation parameter based on physiological parameter |
-
2020
- 2020-11-26 CN CN202011348579.2A patent/CN112494020B/en active Active
- 2020-11-27 WO PCT/CN2020/132275 patent/WO2022110020A1/en active Application Filing
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109475309A (en) * | 2016-07-13 | 2019-03-15 | 悦享趋势科技（北京）有限责任公司 | The method and apparatus for detecting physiological status |
| CN110226923A (en) * | 2018-03-05 | 2019-09-13 | 苏州润迈德医疗科技有限公司 | A method of blood flow reserve score is measured without vasodilator |
| CN109770888A (en) * | 2019-03-19 | 2019-05-21 | 苏州润迈德医疗科技有限公司 | The method instantaneously without waveform ratio is calculated based on pressure sensor and contrastographic picture |
| CN109805949A (en) * | 2019-03-19 | 2019-05-28 | 苏州润迈德医疗科技有限公司 | The method for calculating blood flow reserve score based on pressure sensor and contrastographic picture |
| CN110353645A (en) * | 2019-08-19 | 2019-10-22 | 苏州润迈德医疗科技有限公司 | Measurement is without waveform phase pressure, ratio approach, device, system and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2022110020A1 (en) | 2022-06-02 |
| CN112494020A (en) | 2021-03-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112494020B (en) | Method for acquiring aortic pressure curve of interest and storage medium | |
| CN108186038B (en) | System for calculating coronary blood flow reserve fraction based on arteriography image | |
| CN108511075B (en) | Method and system for non-invasively acquiring fractional flow reserve | |
| JP6685319B2 (en) | Method and apparatus for quantitative flow analysis | |
| US20210236000A1 (en) | Method, device and system for acquiring blood vessel evaluation parameters based on angiographic image | |
| CN110786842B (en) | Method, device, system and storage medium for measuring diastolic blood flow velocity | |
| CN110786841B (en) | Method and device for adjusting maximum hyperemia state flow rate based on microcirculation resistance index | |
| CN112164020A (en) | Method, device, analysis system and storage medium for accurately extracting blood vessel center line | |
| CN112419276B (en) | Method for adjusting contour and center line of blood vessel and storage medium | |
| US20220261997A1 (en) | Method and apparatus for acquiring blood vessel evaluation parameter based on physiological parameter, and storage medium | |
| WO2022110095A1 (en) | Method for acquiring coronary artery vascular evaluation parameters and storage medium | |
| EP4005472A1 (en) | Method and apparatus for correcting blood flow velocity on the basis of interval time between angiographic images | |
| JP2022079550A (en) | Method of operating processor device | |
| CN112132882A (en) | Method and device for extracting blood vessel central line from coronary artery two-dimensional contrast image | |
| CN112472112B (en) | Method, system and storage medium for regulating vascular stenosis | |
| WO2022110019A1 (en) | Aortic pressure waveform image generation method and storage medium | |
| WO2020083390A1 (en) | Method, device and system for acquiring blood flow of large artery on heart surface, and storage medium | |
| CN110929604B (en) | Method, device, system and storage medium for screening flow velocity based on contrast image | |
| CN111815584B (en) | Method and system for acquiring heart gravity center based on CT sequence image | |
| WO2022109907A1 (en) | Method and storage medium for precise acquisition of stenotic lesion range | |
| WO2022110021A1 (en) | Method for extracting stable pressure waveform on the basis of aortic pressure, and storage medium | |
| CN113017667A (en) | Method, device and equipment for quantifying vascular stenosis and readable storage medium | |
| WO2022109909A1 (en) | Method and system for accurately acquiring stenotic lesion interval, and storage medium | |
| JP2015217113A (en) | Blood vessel analysis device, medical image diagnostic device, blood vessel analysis method, and blood vessel analysis program | |
| CN117809852A (en) | Method, device, equipment and medium for determining coronary blood flow reserve fraction |
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 |