CN108564562B - Cardiac cycle acquisition method based on X-ray coronary angiography image - Google Patents
Cardiac cycle acquisition method based on X-ray coronary angiography image Download PDFInfo
- Publication number
- CN108564562B CN108564562B CN201810010148.1A CN201810010148A CN108564562B CN 108564562 B CN108564562 B CN 108564562B CN 201810010148 A CN201810010148 A CN 201810010148A CN 108564562 B CN108564562 B CN 108564562B
- Authority
- CN
- China
- Prior art keywords
- image
- guide catheter
- cardiac cycle
- ray
- coronary angiography
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20036—Morphological image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30048—Heart; Cardiac
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
The invention discloses a cardiac cycle acquisition method based on an X-ray coronary angiography image, which comprises the following steps: step S1: extracting a guide catheter moving image from the contrast image; step S2: carrying out image enhancement preprocessing on the guide catheter moving image; step S3: calculating an optimal segmentation threshold of the enhanced image by a maximum inter-class variance method, and performing threshold segmentation on the image; step S4: calculating the center of the guide catheter of each line in the segmentation image to finish the acquisition of the cardiac cycle; the method has the advantages that the corresponding cardiac cycle is obtained by analyzing the simple harmonic motion of the guide catheter, the method not only can reduce diagnosis and treatment equipment and simplify clinical operation, but also is more beneficial to coronary artery analysis based on X-ray radiography.
Description
Technical Field
The invention belongs to the technical field of medical image processing, and relates to a cardiac cycle acquisition method based on an X-ray coronary angiography image.
Background
The coronary artery is an artery supplying heart blood, is the first pair of branch arteries of a main artery which is emitted from the heart, is divided into a left branch and a right branch, is distributed on the outer surface of the heart and is divided into a plurality of small branches to form abundant capillary vessels for supplying myocardial blood. According to the investigation of the world health organization, cardiovascular diseases (such as hypertensive heart disease, coronary heart disease, arteriosclerosis and the like) are one of the main causes of death in developed countries, and the cardiovascular diseases are listed as one of the biggest diseases threatening human health.
X-ray coronary angiography is one of the common means for diagnosis and treatment of cardiovascular diseases, and has attracted great attention of researchers, and coronary artery analysis based on X-ray coronary angiography is the focus of research of vast researchers. The cardiac cycle is of great importance for coronary analysis, especially when calculating the blood flow rate of the coronary arteries, and in order to calculate it more accurately, the cardiac cycle must be acquired to determine the systolic phase and the diastolic phase of the heart.
Currently, cardiac cycles are typically acquired using an electrocardiographic gating device during X-ray coronary angiography. Since the cardiac gating device is not a standard configuration for the imaging system, additional equipment is required, which not only complicates the clinical procedure, but also may cause asynchrony between the electrocardiogram and the imaging images, which is not conducive to the analysis of the coronary arteries.
Disclosure of Invention
The invention aims at: the method can not only reduce diagnosis and treatment equipment and simplify clinical operation, but also is more beneficial to coronary artery analysis based on X-ray angiography.
The technical scheme of the invention is as follows: a cardiac cycle acquisition method based on an X-ray coronary angiography image comprises the following steps:
step S1: extracting a guide catheter moving image from the contrast image;
step S2: performing image enhancement preprocessing on the guide catheter moving image;
step S3: calculating the optimal segmentation threshold of the enhanced image by a maximum inter-class variance method, and performing threshold segmentation on the image;
step S4: and calculating the center of the guide catheter of each line in the segmentation image to finish the acquisition of the cardiac cycle.
As a preferred technical solution, the specific formula for extracting the guide catheter moving image in step S1 is as follows:
assuming that f (x, y, k) represents the k-th frame image, the guide catheter moving image I (x, y) is expressed as:
in the formula, width represents the width of an X-ray contrast image;
Ixz(x',y',k')、Iyz(x ', y ', k ') are respectively expressed as:
as a preferred technical solution, the specific method of the image enhancement preprocessing in step S2 is as follows:
step A: performing morphological closing operation on the guide catheter moving image by using a circular structural element;
and B: and carrying out image enhancement on the image after the closed operation processing, wherein an image enhancement calculation formula is as follows:
wherein cl (x, y) represents an image after the close operation; org (x, y) denotes a guide catheter moving image.
The invention has the advantages that:
1. the method for acquiring the cardiac cycle based on the X-ray coronary angiography image acquires the corresponding cardiac cycle by analyzing the simple harmonic motion of the guide catheter, can reduce diagnosis and treatment equipment, simplifies clinical operation, is more beneficial to the analysis of the coronary artery based on the X-ray angiography, and brings great convenience to the clinical operation.
Drawings
The invention is further described with reference to the following figures and examples:
FIG. 1 is a flow chart of the present invention;
FIG. 2 is a guide catheter motion image of the present invention;
fig. 3 is an image of a cardiac cycle of the present invention.
Detailed Description
Example (b): when X-ray coronary angiography is performed, a part of guide catheter enters a radiography image, the guide catheter is attached to the surface of a heart, the basic motion of the heart can be indirectly reflected by the simple harmonic motion of the guide catheter, and the analysis of the running condition of the guide catheter is facilitated due to the fact that the guide catheter belongs to rigid motion.
Referring to fig. 1, the method for acquiring a cardiac cycle based on an X-ray coronary angiography image of the present invention includes the following steps:
step S1: extracting a guide catheter moving image from the contrast image, as shown in fig. 2; assuming that f (x, y, k) represents the k-th frame image, the guide catheter moving image I (x, y) is represented as:
in the formula, width represents the width of an X-ray contrast image;
Ixz(x',y',k')、Iyz(x ', y ', k ') are respectively expressed as:
step S2: the image enhancement preprocessing is carried out on the guide catheter moving image, and the specific method comprises the following steps:
step A: performing morphological closing operation on the guide catheter moving image by using a circular structural element;
and B: and carrying out image enhancement on the image after the closed operation processing, wherein an image enhancement calculation formula is as follows:
in the formula, cl (x, y) represents an image after the closing operation; org (x, y) denotes a guide catheter moving image;
step S3: calculating the optimal segmentation threshold of the enhanced image by a maximum inter-class variance method, and performing threshold segmentation on the image;
step S4: the center of the guiding catheter of each line in the segmentation image is calculated, and the acquisition of the cardiac cycle is completed, which is shown in reference to fig. 3.
The foregoing embodiments are merely illustrative of the principles and utilities of the present invention and are not intended to limit the invention. Any person skilled in the art can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Accordingly, it is intended that all equivalent modifications or changes which can be made by those skilled in the art without departing from the spirit and technical spirit of the present invention be covered by the claims of the present invention.
Claims (3)
1. A cardiac cycle acquisition method based on an X-ray coronary angiography image is characterized by comprising the following steps:
step S1: extracting a guide catheter moving image from the contrast image;
step S2: image enhancement preprocessing is carried out on a guide catheter moving image, and the image enhancement preprocessing comprises the following steps: performing morphological closing operation on the guide catheter moving image by using a circular structural element;
step S3: calculating the optimal segmentation threshold of the enhanced image by a maximum inter-class variance method, and performing threshold segmentation on the image;
step S4: and calculating the center of the guide catheter of each line in the segmentation image to finish the acquisition of the cardiac cycle.
2. The method for acquiring a cardiac cycle based on an X-ray coronary angiography image according to claim 1, wherein the specific formula for extracting the motion image of the guide catheter in step S1 is as follows:
assuming that f (x, y, k) represents the k-th frame image, the guide catheter moving image I (x, y) is represented as:
in the formula, width represents the width of an X-ray contrast image;
Ixz(x',y',k')、Iyz(x ', y ', k ') are respectively expressed as:
3. the method for acquiring cardiac cycle based on X-ray coronary angiography image according to claim 1, wherein the method of image enhancement preprocessing in step S2 further comprises:
and B: and carrying out image enhancement on the image after the closed operation processing, wherein the image enhancement calculation formula is as follows:
in the formula, cl (x, y) represents an image after the closing operation; org (x, y) denotes a guide catheter moving image.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810010148.1A CN108564562B (en) | 2018-01-05 | 2018-01-05 | Cardiac cycle acquisition method based on X-ray coronary angiography image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810010148.1A CN108564562B (en) | 2018-01-05 | 2018-01-05 | Cardiac cycle acquisition method based on X-ray coronary angiography image |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108564562A CN108564562A (en) | 2018-09-21 |
CN108564562B true CN108564562B (en) | 2022-06-17 |
Family
ID=63529662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810010148.1A Active CN108564562B (en) | 2018-01-05 | 2018-01-05 | Cardiac cycle acquisition method based on X-ray coronary angiography image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108564562B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113768476B (en) * | 2021-09-13 | 2024-03-26 | 复旦大学附属中山医院 | Cardiac cycle phase detection method, apparatus and computer program product |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2848093A1 (en) * | 2002-12-06 | 2004-06-11 | Ge Med Sys Global Tech Co Llc | Detection of cardiac cycle uses analysis of successive x-ray images of coronary blood vessels to detect movement |
CN101283910A (en) * | 2008-06-05 | 2008-10-15 | 华北电力大学 | Method for obtaining the coronary artery vasomotion information |
CN101541245A (en) * | 2006-11-22 | 2009-09-23 | 皇家飞利浦电子股份有限公司 | Combining x-ray with intravascularly acquired data |
CN105741299A (en) * | 2016-02-02 | 2016-07-06 | 河北大学 | Coronary artery CT angiography image segmentation method |
-
2018
- 2018-01-05 CN CN201810010148.1A patent/CN108564562B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2848093A1 (en) * | 2002-12-06 | 2004-06-11 | Ge Med Sys Global Tech Co Llc | Detection of cardiac cycle uses analysis of successive x-ray images of coronary blood vessels to detect movement |
CN101541245A (en) * | 2006-11-22 | 2009-09-23 | 皇家飞利浦电子股份有限公司 | Combining x-ray with intravascularly acquired data |
CN101283910A (en) * | 2008-06-05 | 2008-10-15 | 华北电力大学 | Method for obtaining the coronary artery vasomotion information |
CN105741299A (en) * | 2016-02-02 | 2016-07-06 | 河北大学 | Coronary artery CT angiography image segmentation method |
Also Published As
Publication number | Publication date |
---|---|
CN108564562A (en) | 2018-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | A novel neural network model for blood pressure estimation using photoplethesmography without electrocardiogram | |
US11989652B2 (en) | Methods and systems to configure and use neural networks in characterizing physiological systems | |
CA2963866C (en) | Methods, systems and computer program products for determining hemodynamic status parameters using signals derived from multispectral blood flow and perfusion imaging | |
McDuff et al. | Remote detection of photoplethysmographic systolic and diastolic peaks using a digital camera | |
CN108364290B (en) | Method, medium, and system for analyzing a sequence of images of periodic physiological activity | |
US8488910B2 (en) | Image provision for registration | |
Revanur et al. | The first vision for vitals (v4v) challenge for non-contact video-based physiological estimation | |
CN109431531A (en) | Blood vessel segmentation method and device and computer installation based on Perfusion Imaging | |
US20230346232A1 (en) | Methods, Systems and Computer Program Products for Calculating MetaKG Signals for Regions Having Multiple Sets of Optical Characteristics | |
Tadi et al. | A new algorithm for segmentation of cardiac quiescent phases and cardiac time intervals using seismocardiography | |
Lan | Performer: A novel ppg-to-ecg reconstruction transformer for a digital biomarker of cardiovascular disease detection | |
CN113080991A (en) | Method, system, diagnosis device and storage medium for predicting and diagnosing heart failure based on CNN model and LSTM model | |
CN108564562B (en) | Cardiac cycle acquisition method based on X-ray coronary angiography image | |
Basak et al. | A novel deep learning technique for morphology preserved fetal ECG extraction from mother ECG using 1D-CycleGAN | |
WO2020198522A1 (en) | Electrocardiogram waveform reconstruction from photoplethysmogram | |
Bhattacharjee | Cuff-less blood pressure estimation from electrocardiogram and photoplethysmography based on VGG19-LSTM network | |
Hennersperger et al. | Vascular 3D+ T freehand ultrasound using correlation of doppler and pulse-oximetry data | |
Das et al. | Time-Frequency Learning Framework for rPPG Signal Estimation Using Scalogram Based Feature Map of Facial Video Data | |
O'Malley et al. | Image-based frame gating of IVUS pullbacks: A surrogate for ecg | |
Liang et al. | A deep learning method for continuous noninvasive blood pressure monitoring using photoplethysmography | |
Tapotee et al. | M2ECG: Wearable Mechanocardiograms to Electrocardiogram Estimation Using Deep Learning | |
Eranna et al. | Signal processing based diagnosis of cardiovascular anomalies | |
EP4260801A2 (en) | High resolution blood perfusion imaging using a camera and a pulse oximeter | |
Rahimi | Investigation of Trimodality Image Registration for Guiding Cardiovascular Diseases Treatment Procedure | |
US20230371827A1 (en) | Machine learning techniques for estimating carotid-femoral pulse wave velocity and/or vascular age from single-site arterial waveform measurements |
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 |