CN115530791A - Intracardiac shunt electrical impedance image analysis method and system based on saline angiography - Google Patents

Intracardiac shunt electrical impedance image analysis method and system based on saline angiography Download PDF

Info

Publication number
CN115530791A
CN115530791A CN202211532638.0A CN202211532638A CN115530791A CN 115530791 A CN115530791 A CN 115530791A CN 202211532638 A CN202211532638 A CN 202211532638A CN 115530791 A CN115530791 A CN 115530791A
Authority
CN
China
Prior art keywords
saline
sequence diagram
image sequence
impedance
central point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211532638.0A
Other languages
Chinese (zh)
Other versions
CN115530791B (en
Inventor
何怀武
陈翔宇
隆云
池熠
翁利
杜斌
招展奇
袁思依
弥亮钰
刘旺林
薛月川
杨莹莹
刘圣均
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Peking Union Medical College Hospital Chinese Academy of Medical Sciences
Original Assignee
Peking Union Medical College Hospital Chinese Academy of Medical Sciences
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Peking Union Medical College Hospital Chinese Academy of Medical Sciences filed Critical Peking Union Medical College Hospital Chinese Academy of Medical Sciences
Priority to CN202211532638.0A priority Critical patent/CN115530791B/en
Publication of CN115530791A publication Critical patent/CN115530791A/en
Application granted granted Critical
Publication of CN115530791B publication Critical patent/CN115530791B/en
Priority to PCT/CN2023/099037 priority patent/WO2023237034A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/053Measuring electrical impedance or conductance of a portion of the body
    • A61B5/0536Impedance imaging, e.g. by tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Physics & Mathematics (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Signal Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The invention relates to a saline angiography-based intracardiac shunt electrical impedance image analysis method and a saline angiography-based intracardiac shunt electrical impedance image analysis system. The method comprises the following steps: acquiring a patient saline contrast impedance curve; obtaining an image sequence diagram based on a saline contrast impedance curve, wherein the image sequence diagram comprises a right cardiac phase image sequence diagram and/or a left cardiac phase image sequence diagram; extracting the central point of each image sequence diagram, and calculating the horizontal position relation between the central point and the central point of the first sequence diagram; and outputting the classification result of the intra-cardiac shunt of the patient based on the position relation. The method aims to realize the intra-cardiac shunt classification of the images based on the horizontal position relation between the central point and the central point of the first sequence diagram, and aims to discover the image analysis capability and the potential application value in the auxiliary diagnosis of the heart diseases.

Description

Intracardiac shunt electrical impedance image analysis method and system based on saline angiography
Technical Field
The present invention relates to the field of image analysis in clinical medicine, and more particularly, to a method, system, apparatus, computer-readable storage medium, and applications thereof for intra-cardiac shunt electrical impedance image analysis based on saline angiography.
Background
Some congenital diseases have defects in the interval between the left heart cavity and the right heart cavity, resulting in shunting between the left heart cavity and the right heart cavity, and the unequal volume of blood in systemic circulation and pulmonary circulation is called intracardiac shunting. Wherein, blood oxygen data obtained by invasive evaluation of right cardiac catheterization can determine the shunting position, shunting direction and shunting amount between left and right cardiac chambers. In addition, the oscillating saline contrast study is of great value for detecting intracardiac shunts.
At present, the electrical impedance tomography monitoring space resolution is low, and the electrical impedance tomography monitoring space resolution cannot be applied to ventricular radiography imaging all the time to assist diagnosis of severe heart diseases. The bedside ventricular imaging based on Electrical Impedance Technology (EIT) is of poor quality, with poor spatial resolution but high temporal resolution. The early research develops a saline angiography technology electrical impedance imaging technology method, saline is injected through a central venous catheter, apnea is carried out at the same time, the change of thoracic resistance signals is collected, a local resistance-time change curve of the saline angiography is established, and the pulmonary blood perfusion imaging is realized. On the basis, when the ventricular atrium has structural abnormalities such as deletion and the like, the abnormality of a saline angiography curve can occur, and on the basis, the method aims to provide a novel technical method which is beside a bed, non-invasive, non-radiative and more practical and is based on the saline angiography for intracardiac shunt electrical impedance image analysis, and the clinical application value of the method in heart disease diagnosis is expanded.
Disclosure of Invention
It is an object of the present application to provide a saline contrast based intracardiac shunt electrical impedance image analysis method, a saline contrast based intracardiac shunt electrical impedance image analysis system, a saline contrast based intracardiac shunt electrical impedance image analysis device, a computer readable storage medium and applications thereof, which are aimed at analyzing key features related to cardiopulmonary disorders by a saline enhanced electrical impedance ventricular contrast imaging method, and implementing intracardiac shunt classification of images based on the horizontal position relationship between the proposed center point and the center point of the first sequence diagram, so as to find out the potential application value thereof in assisted diagnosis of heart diseases, and to provide more sufficient support for selection of treatment decisions of patients.
According to a first aspect of the present application, an embodiment of the present application provides a saline contrast based intracardiac shunt electrical impedance image analysis method, comprising:
acquiring a patient saline contrast impedance curve;
obtaining an image sequence diagram based on the saline contrast impedance curve, wherein the image sequence diagram comprises a right cardiac phase image sequence diagram and/or a left cardiac phase image sequence diagram;
extracting the central point of each image sequence diagram, and calculating the horizontal position relation between the central point and the central point of the first sequence diagram;
and outputting the classification result of the intra-cardiac shunt of the patient based on the position relation.
Further, the obtaining of the image sequence diagram based on the saline contrast impedance curve is obtaining of the image sequence diagram based on the saline contrast impedance curve by image reconstruction, wherein the image reconstruction includes any one or more of the following methods: projection reconstruction, light and shade recovery shape, stereoscopic vision reconstruction and laser ranging reconstruction.
In one embodiment, the method for intracardiac shunt electrical impedance image analysis based on saline angiography further comprises:
acquiring a patient saline contrast impedance curve;
obtaining a right heart phase image sequence chart based on the saline contrast impedance curve;
extracting the central point of each image sequence chart in the right heart phase image sequence chart, and calculating the horizontal position relation between the central point and the central point of a first sequence chart in the right heart phase image sequence chart;
and outputting whether the patient has the classification result of right-to-left intracardiac shunting or not based on whether the position relation is on the left side or not.
In one embodiment, the method for intracardiac shunt electrical impedance image analysis based on saline angiography further comprises:
acquiring a patient saline contrast impedance curve;
obtaining a left cardiac phase image sequence diagram based on the saline angiography impedance curve;
extracting the central point of each image sequence diagram in the left cardiac phase image sequence diagram, and calculating the horizontal position relation between the central point and the central point of a first sequence diagram in the left cardiac phase image sequence diagram;
and outputting whether the patient has a left-to-right intracardiac shunt classification result or not based on whether the position relation is right or not.
Further, the right cardiac phase image sequence diagram is obtained by reconstructing a right ventricular angiography impedance curve of a time period from T0 to T1 in the saline angiography impedance curve; and the left cardiac phase image sequence chart is obtained by reconstructing a left ventricle blood flow contrast impedance curve of a T2-T3 time period in the saline contrast impedance curve.
In one embodiment, the center point of the sequence diagram is determined as follows:
Figure 259361DEST_PATH_IMAGE001
wherein i belongs to the heart region H, x i Is the abscissa, f, of pixel point i i The optimal impedance descent slope, coH (t), of the least-squares fit curve for the ith pixel point k ) Represents t k Time of dayPosition of the abscissa of the heart center, t k Is the kth time window.
Still further, the heart region H is a pixel point impedance decreasing slope f>20%×f max Pixel point of (f) max The slope of the pixel point with the maximum slope is calculated by the following formula:
Figure 550665DEST_PATH_IMAGE002
wherein,
Figure 286540DEST_PATH_IMAGE003
is composed of
Figure 406811DEST_PATH_IMAGE004
The relative impedance value of the time pixel i, e is the intercept.
According to a second aspect of the present application, an embodiment of the present application provides a saline contrast based intracardiac shunt electrical impedance image analysis system, comprising:
the acquisition module is used for acquiring a saline angiography impedance curve;
an image reconstruction module, configured to obtain an image sequence diagram based on the saline angiography impedance curve, where the image sequence diagram includes a right cardiac phase image sequence diagram and/or a left cardiac phase image sequence diagram;
the calculation module is used for extracting the central point of each image sequence diagram and calculating the horizontal position relation between the central point and the central point of the first sequence diagram;
and the output module is used for outputting the classification result of the intra-cardiac shunt of the patient based on the position relation.
According to a third aspect of the present application, an embodiment of the present application provides a saline contrast based intracardiac shunt electrical impedance image analysis apparatus, mainly comprising:
a memory and a processor;
the memory is to store program instructions;
the processor is configured to invoke program instructions that, when executed, perform the intracardiac shunt electrical impedance image analysis method based on saline contrast described above.
According to a fourth aspect of the present application, an embodiment of the present application provides a computer-readable storage medium having stored thereon a computer program for intra-cardiac shunt electrical impedance image analysis based on a saline contrast, which, when executed by a processor, implements the above-described intra-cardiac shunt electrical impedance image analysis based on a saline contrast.
The application of the device or the system in the intelligent prediction analysis of the intracardiac streaming situation; optionally, the intracardiac diversion condition includes two classification results, i.e., left-right intracardiac diversion and right-left intracardiac diversion.
The intracardiac left-to-right shunting results in shunting of blood from the left atrium to the right atrium.
The intracardiac right-to-left shunt is characterized by a decrease in pulmonary circulatory blood flow, a pressure difference between the enlarged right atrium and the functional right ventricle causes a light heart, moderately increases the decrease in pulmonary circulatory blood flow, and a pressure difference between the enlarged right atrium and the functional right ventricle causes a right-to-left intracardiac shunt.
The use of the apparatus or system described above for impedance radiography of a chamber of a strong brine ventricle at the bedside of a patient; optionally, the right cardiac phase and the left cardiac phase are determined by imaging of different sequences of impedance of the strong brine angiography, and the bedside rapid diagnosis of heart and lung diseases with serious symptoms such as intracardiac shunt is further performed based on image analysis and position relation judgment;
use of the above-described device or system for assisted diagnosis of a heart related disorder. The auxiliary diagnosis judges whether the intra-cardiac right-left shunt exists or not by monitoring whether the phenomenon that the saline aggregation center moves left exists in the right cardiac phase in the saline impedance radiography process or not; and further comprising the step of judging whether the intracardiac left-to-right shunt exists or not by monitoring whether the phenomenon that the saline water aggregation center moves to the right exists or not in the left cardiac phase impedance dilution curve. Specifically, the auxiliary diagnosis includes heart-related diseases such as atrial septal defect, ventricular septal defect, patent ductus arteriosus, and Ebbstein syndrome.
Among them, atrial septal defect, ventricular septal defect and patent ductus arteriosus are heart diseases shunted from left to right inside the heart. Taking the ventricular block as an example, left-to-right shunting may cause the blood that should be supplied to the whole body to be shunted from the ventricular septum to the right ventricle midway, reducing the blood flow supplied to the whole body. If the defect is large, the shunt flow is large, the pulmonary blood flow is increased, the chance of pulmonary infection is increased, pneumonia is caused, heart failure is induced, and right heart insufficiency, right heart failure and pulmonary hypertension are easy to occur.
The atrial septal defect refers to a defect that the development, absorption and fusion of the atrial septal are abnormal in the embryonic development process, so that a left patent is remained between the left atrium and the right atrium, and is one of the most common congenital heart diseases of adults.
The ventricular septal defect is a series of problems caused by the fact that when the ventricular septal hole appears, blood of the left ventricle and the right ventricle are fused with each other. The left and right ventricular septa of the heart are referred to as the ventricular septa, and a complete ventricular septum is capable of preventing left and right ventricular blood flow communication.
The Patent Ductus Arteriosus (PDA), an abnormal passage between the descending aorta and the pulmonary artery, is a common congenital heart disease, and is often found in premature infants.
The narrow form in the Ebbistatin syndrome, i.e., right-to-left intracardiac shunt. At present, diagnosis of the Ebbistatin syndrome, clinical symptoms and physical signs may give important hints, but accurate diagnosis depends on various auxiliary examinations of electrocardiogram, X-ray examination and echocardiogram, especially selective angiography.
The ventricular impedance imaging technology based on the saline angiography realizes the intracardiac shunt classification of the patient image through the horizontal position relation of the central point and the central point of the first sequence diagram, is a bedside, noninvasive, radiationless and more practical invention, realizes the automatic classification of intracardiac shunt through machine learning, has strong innovation, and has beneficial promotion effect on the analysis and research of cardiopulmonary disease image data.
The application has the advantages that:
1. the method obtains a right cardiac phase image sequence diagram and a left cardiac phase image sequence diagram respectively through saline angiography impedance curves of a patient in T0-T1 time period and T2-T3 time period, obtains a classification result of the intra-cardiac shunt of the patient based on the horizontal position relation between the center point of the right cardiac phase image sequence diagram and the center point of the first sequence diagram, and improves the precision and the depth of data analysis objectively;
2. the thoracic cavity impedance transformation monitoring system creatively utilizes a dynamic process that an imaging agent (strong brine) firstly passes through central circulation to continuously monitor the transformation of thoracic cavity impedance, obtains impedance series images of a heart region by analyzing the characteristics of impedance-dilution curves generated by strong brine radiography in different time phases, and has the advantages of no wound, no radiation and obvious aging;
3. the application creatively discloses a right cardiac phase image sequence diagram and a left cardiac phase image sequence diagram which are respectively obtained by a ventricular impedance imaging technology based on saline angiography, whether shunt results in the right left direction in the heart and/or the left right direction in the heart exist or not is obtained through the proposed position relation, and the application is more accurately applied to auxiliary analysis of occurrence and development of diseases related to cardiopulmonary disease image data in view of important research significance of the intra-cardiac shunt results on diagnosis and prevention and control of heart diseases.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic flow chart diagram of a saline contrast-based intracardiac shunt electrical impedance image analysis method according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a method for constructing a sequence diagram of a left and right cardiac phase image based on a saline contrast imaging according to an embodiment of the present invention;
FIG. 3 is a schematic right heart imaging impedance representation of an example atrial septal defect endocardial shunt patient provided by an embodiment of the present invention;
FIG. 4 is a schematic diagram of a sequence chart construction for different images based on a saline contrast curve according to an embodiment of the present invention;
FIG. 5 is a schematic diagram of a right cardiac phase image sequence chart analysis method based on a saline contrast medium according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of a left cardiac image sequence chart analysis method based on a saline contrast system according to an embodiment of the present invention;
fig. 7 is a schematic diagram of an intracardiac shunt electrical impedance image analysis device based on saline contrast provided by an embodiment of the invention.
Detailed Description
In order to make the technical solutions of the present invention better understood, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention.
In some of the flows described in the present specification and claims and in the above-described figures, a number of operations are included that occur in a particular order, but it should be clearly understood that these operations may be performed out of order or in parallel as they occur herein, with the order of the operations, e.g., S101, S102, etc., merely being used to distinguish between various operations, and the order of the operations itself does not represent any order of performance. Additionally, the flows may include more or fewer operations, and the operations may be performed sequentially or in parallel.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the application provides a saline angiography-based intracardiac shunt electrical impedance image analysis method, a saline angiography-based intracardiac shunt electrical impedance image analysis system, a saline angiography-based intracardiac shunt electrical impedance image analysis device and a computer-readable storage medium. The intra-cardiac shunt electrical impedance image analysis equipment based on the saline angiography can be equipment such as a terminal or a server. The terminal can be terminal equipment such as a smart phone, a tablet Computer, a notebook Computer, a Personal Computer (PC for short) and the like. The server may be an independent physical server, a server cluster or a distributed system formed by a plurality of physical servers, or a cloud server providing basic cloud computing services such as cloud service, a cloud database, cloud computing, cloud storage, network service, cloud communication, middleware service, domain name service, security service, content Delivery Network (CDN) and a big data and artificial intelligence platform. The terminal and the server may be directly or indirectly connected through wired or wireless communication, and the application is not limited herein.
Fig. 1 is a schematic flow chart of a saline contrast-based intracardiac shunt electrical impedance image analysis method, which specifically comprises the following steps:
s101: a patient saline contrast impedance curve is acquired.
In one embodiment, the patient's saline contrast impedance profile includes a series of electrical impedance changes based on the flow of saline through the superior vena cava → right atrium, right ventricle → pulmonary circulation → left atrium, left ventricle → aorta.
In a specific embodiment, the acquired patient saline contrast impedance curve is based on a patient inclusion study that follows the steps of:
(1) Breath holding test, the right heart and lung require at least 8 seconds; the influence of the cardiopulmonary transit time requires a minimum of more than 15 seconds for better construction of left-heart phase imaging, and the breath holding time is further prolonged if necessary;
(2) Injecting saline water to carry out ventricular blood flow radiography, and continuously acquiring the change of electrical impedance signals of the chest during breath holding;
further, the specific operation of step (1) is as follows: breath holding tests, requiring a minimum of more than 8 seconds (breath holding by expiration or inspiration for 8 seconds during mechanical ventilation by the ventilator; breath holding by the voluntary breathing patient for 8 seconds); after the breath holding test is passed, an EIT examination can be performed by a saline angiography. Wherein the concentration of the injection saline is 10%, and the injection amount is 10ml.
Still further, the specific operation of step (2) is as follows: connecting a test subject with a pulmonary electrical impedance monitoring instrument, preparing 10% NaCl 10ml, and confirming that the test subject establishes a central venous catheter; 1 second after breath holding, 10ml of 10% NaCl is synchronously and quickly injected from a central venous catheter to the body for pulmonary blood perfusion radiography; during which changes in the thoracic electrical impedance signal are continuously acquired.
S102: obtaining image sequence maps based on the saline contrast impedance curve, the image sequence maps including right cardiac phase image sequence maps and/or left cardiac phase image sequence maps.
In one embodiment, the right cardiac phase image sequence chart is reconstructed based on the right ventricle blood flow contrast impedance curve of the T0-T1 time segment in the saline contrast impedance curve; the left cardiac phase image sequence chart is reconstructed based on the left ventricular angiographic impedance curve of the T2-T3 time segment in the saline angiographic impedance curve.
In one embodiment, the acquisition of the right cardiac phase image sequence map is as follows: as shown in fig. 2, the overall resistance curve begins to decrease during breath holding as the starting point (T0) of saline entering the heart, and then the impedance of the lung region of a certain ROI begins to decrease as the time point (T1) of saline entering the pulmonary vessels reaching the lung region, therefore, the impedance curve of the T0-T1 time period mainly reflects the time period when saline enters the right heart, the resistance-time variation curve of the T0-T1 time period is applied to calculate the slope, and the right heart imaging composition is constructed by analyzing each time window in the T0-T1 time period. Specifically, the slope time window of the calculated resistance-time curve of the T0-T1 time period is 0.5 seconds, the step size is 0.5 seconds, and the time window division of the specific T0-T1 time period is from T0 to T0+0.5 seconds, from T0+0.5 to T0+1.0 seconds, and so on to T1-0.5 to T1, that is, by analyzing the right heart image constructed in each 0.5 second window in the T0-T1 time period. More specifically, the calculation of the optimal slope of the ith pixel point in the right heart phase image sequence diagram is obtained through an ith pixel point perfusion volume calculation formula:
Figure 544532DEST_PATH_IMAGE005
wherein, t n The nth time window of the time period T0-T1, a i Fitting for least squaresOptimum slope of the curve, Δ
Figure 6737DEST_PATH_IMAGE006
The relative impedance value of the pixel point i at a certain time, and b is the intercept.
In one embodiment, the left cardiac image sequence map is obtained as follows: as shown in fig. 2, when the saline flows back to the pixel point in the heart region, the position where the impedance-time curve starts to fall again is defined as time point T2, and the second valley time point is defined as T3, so that the impedance curve of the T2-T3 time period mainly reflects the time when the saline enters the gathered left heart, the resistance-time change curve of the T2-T3 time period is applied, the slope is calculated, and the left heart imaging is composed by analyzing the individual left heart images constructed in each time window of the T2-T3 time period. The heart region is determined by making a during T0-T1 max Defining the slope a of the pixel point for the slope of the pixel point with the maximum slope>20%×a max The pixel point of (a) is a heart region. During the period from T2 to T3, the impedance falling slope c of the ith pixel point in the left cardiac phase image sequence diagram is calculated by the following formula:
Figure 26645DEST_PATH_IMAGE007
wherein, t m The m-th time window, c, being the period of T2-T3 i Fitting the best slope of the curve, Δ, for least squares
Figure 216187DEST_PATH_IMAGE008
The relative impedance value of the pixel point i at a certain time in the time period from T2 to T3, and d is the intercept.
S103: and extracting the central point of each image sequence diagram, and calculating the horizontal position relation between the central point and the central point of the first sequence diagram.
In one embodiment, the horizontal position of the center point of the sequence diagram is determined as follows:
Figure 473993DEST_PATH_IMAGE009
wherein i belongs to the heart region H, x i Is the abscissa, f, of pixel point i i The optimal impedance descent slope, coH (t), of the least-squares fit curve for the ith pixel point k ) Represents t k The abscissa position of the heart center at time, t k Is the kth time window.
Still further, the heart region H is the pixel point impedance decreasing slope f>20%×f max Pixel point of (f) max Specifically, the impedance decreasing slope f is calculated by the following formula:
Figure 107100DEST_PATH_IMAGE010
wherein,
Figure 817567DEST_PATH_IMAGE003
is composed of
Figure 358270DEST_PATH_IMAGE004
The relative impedance value of the time pixel i, e is the intercept.
In one embodiment, FIG. 3 illustrates the determination of the positional relationship by calculating the horizontal positional offset of the center of the right heart imaging impedance.
In a more general embodiment, each image sequence is obtained from impedance images based on a series of saline contrast impedance curve constructions, and the constructed image sequence further includes different sequences of visualizations constructed based on the saline contrast time-impedance curve as shown in FIG. 4: saline focus right cardiac phase of right heart, pulmonary phase of bipulmonary perfusion, saline circulation through the lungs back to left cardiac phase of left cardiac series imaging phase.
S104: and outputting the classification result of the intra-cardiac shunt of the patient based on the position relation.
In one embodiment, the right heart imaging is used to determine whether there is right-to-left intracardiac shunting by determining whether there is a left shift in the right heart imaging impedance.
In one embodiment, the application of left heart visualization determines the presence of left-to-right intracardiac shunting by determining the presence of a saline focus shift to the right.
Fig. 3 shows that the right imaging impedance center shift is calculated according to the position relationship to determine whether there is right-left shunt, where the left side is left shift (+ 9.2%) of the right imaging impedance center of an atrial septal defect endocardial shunt patient, and the right side is not obvious left shift (-2.4%) of the normal control right imaging impedance center. Therefore, whether intracardiac right-to-left shunting exists is judged by monitoring whether the phenomenon that the saline concentration center moves left exists in the right heart phase in the saline impedance radiography process, and the method is bedside, non-invasive, non-radiative and more practical.
In a more complete embodiment, the intra-cardiac shunt electrical impedance image analysis method based on saline angiography shown in fig. 1 further comprises the analysis methods shown in fig. 5 and 6.
Fig. 5 is a schematic diagram of a right cardiac phase image sequence diagram analysis method based on saline contrast provided in an embodiment of the present invention, which includes the following main steps:
s501: a patient saline contrast impedance curve is acquired.
S502: right cardiac phase image sequence plots were obtained based on the saline contrast impedance curves.
More specifically, a resistance-time change curve of a time period from T0 to T1 is used for calculating slope time windows, and right heart images corresponding to the time windows are constructed to obtain a right heart phase image sequence diagram. Preferably, the calculation window of the slope of the resistance-time curve is 0.5 seconds, and the step size is 0.5 seconds, i.e., the time window of the time period T0-T1 is divided from T0 to T0+0.5 seconds, T0+0.5 to T0+1.0 seconds, and so on to T1-0.5 to T1.
S503: and extracting the central point of each image sequence chart in the right heart phase image sequence chart, and calculating the horizontal position relation between the central point and the central point of the first sequence chart.
Each image sequence in the right ventricle image sequence is corresponding to one right ventricle image constructed in each 0.5 second time window in the T0-T1 period by analyzing step S502.
Further, the first sequence diagram in step S503 refers to the right heart image constructed from the time window from T0 to T0+0.5 seconds, and the horizontal position of the center point of the first sequence diagram is calculated as follows:
Figure 725709DEST_PATH_IMAGE011
still further, the horizontal position of the center point of each image sequence in the right cardiac phase image sequence is calculated as follows:
Figure 529717DEST_PATH_IMAGE012
wherein x is i Is the abscissa of pixel i, which belongs to the heart region, coH
Figure 461901DEST_PATH_IMAGE013
Is t n Heart center abscissa at time t n The nth time window of the time period T0-T1.
Still further, the positional relationship is compared
Figure 275136DEST_PATH_IMAGE014
And
Figure 491223DEST_PATH_IMAGE015
is obtained relative to the offset position, i.e. when
Figure 466132DEST_PATH_IMAGE016
When, the position relation is the left side; when in use
Figure 151191DEST_PATH_IMAGE017
The positional relationship is not to the left.
S504: and outputting whether the patient has the classification result of right-to-left intracardiac shunting based on whether the position relation is on the left side.
When the position relation is on the left side, outputting the classification result of the patient with the intracardiac right-to-left shunt; when the position relation is not the left side, the classification result of the right-to-left intracardiac shunt of the patient is output.
Fig. 6 is a schematic diagram of a left cardiac phase image sequence chart analysis method based on saline contrast according to an embodiment of the present invention, which includes the following main steps:
s601: a patient saline contrast impedance curve is acquired.
S602: a left cardiac phase image sequence plot was obtained based on the saline contrast impedance curve.
Wherein, the left cardiac phase image sequence diagram is reconstructed based on the left ventricle blood flow contrast impedance curve of the T2-T3 time period in the saline contrast impedance curve.
S603: and extracting the central point of each image sequence diagram in the left cardiac phase image sequence diagram, and calculating the horizontal position relation between the central point and the central point of the first sequence diagram.
Each image sequence in the left cardiac image sequence corresponds to a single left cardiac image constructed in each time window in the T2-T3 period by analyzing step S602.
Further, the first sequence diagram in step S603 is the left-heart visualization constructed by the first time window, and the horizontal position of the center point of the first sequence diagram is the same as the horizontal position of the center point of the first sequence diagram
Figure 502538DEST_PATH_IMAGE018
Is calculated as follows:
Figure 307552DEST_PATH_IMAGE019
still further, the horizontal position of the center point of each image sequence in the left cardiac image sequence is calculated as follows:
Figure 718942DEST_PATH_IMAGE020
wherein x is i Is the abscissa of pixel i, which belongs to the heart region, coH
Figure 422456DEST_PATH_IMAGE021
Is t m Heart center abscissa at time t m The mth time window being the time period of T2-T3。
Still further, the positional relationship is compared
Figure 311914DEST_PATH_IMAGE022
And
Figure 253325DEST_PATH_IMAGE023
is obtained relative to the offset position, i.e. when
Figure 84884DEST_PATH_IMAGE024
When the position relation is not right; when in use
Figure 478956DEST_PATH_IMAGE025
The positional relationship is right.
S604: and outputting whether the patient has the classification result of left-to-right intracardiac shunt or not based on whether the position relation is right or not.
When the position relation is right, outputting the classification result of the patient with intracardiac left-to-right shunt; when the position relation is not right, outputting the classification result that the patient does not have intracardiac left-right shunt.
The method is feasible to be used for intracardiac shunt classification, and the same principle shows that whether intracardiac left-right shunt exists is judged by monitoring whether the phenomenon that the saline aggregation center moves to the right exists in the left cardiac phase impedance dilution curve or not, and whether intracardiac right-left shunt exists or not is judged by monitoring whether the phenomenon that the saline aggregation center moves to the left exists in the right cardiac phase in the saline impedance radiography process, so that more accurate information is provided for the prediction analysis of the intracardiac shunt of a patient, the specific situation and effect of the intracardiac shunt are better reflected, and the method is a bedside, noninvasive, radiationless and more practical method, so that the method is more favorable in the aspects of auxiliary diagnosis of the saline radiography data applied to heart diseases and auxiliary analysis of occurrence and development of the diseases.
The invention provides a saline angiography-based intracardiac shunt electrical impedance image analysis system, which comprises:
an acquisition module for acquiring a patient saline contrast impedance curve;
the image reconstruction module is used for obtaining an image sequence diagram based on a patient saline angiography impedance curve, and the image sequence diagram comprises a right cardiac phase image sequence diagram and/or a left cardiac phase image sequence diagram;
the calculation module is used for extracting the central point of each image sequence diagram and calculating the horizontal position relation between the central point and the central point of the first sequence diagram;
and the output module outputs the classification result of the intra-cardiac shunt of the patient based on the position relation.
Fig. 7 is an intracardiac shunt electrical impedance image analysis device based on saline contrast provided by an embodiment of the invention, which comprises:
a memory and a processor;
the apparatus may further include: an input device and an output device.
The memory, processor, input device and output device may be connected by a bus or other means, such as the bus connection shown in fig. 7;
wherein the memory is used for storing program instructions; the processor is configured to invoke program instructions that, when executed, are configured to perform the above-described saline contrast-based intracardiac shunt electrical impedance image analysis method.
The invention provides a computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, implements the above-mentioned intra-cardiac shunt electrical impedance image analysis method based on saline angiography.
It can be clearly understood by those skilled in the art that, for convenience and brevity of description, the specific working processes of the system, the apparatus and the module described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other manners. For example, the above-described apparatus embodiments are merely illustrative, and for example, a division of modules is merely a division of logical functions, and an actual implementation may have another division, for example, a plurality of modules or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or modules, and may be in an electrical, mechanical or other form.
Modules described as separate parts may or may not be physically separate, and parts displayed as modules may or may not be physical modules, may be located in one place, or may be distributed on a plurality of network modules. Specifically, some or all of the modules are selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional modules in the embodiments of the present invention may be integrated into one processing module, or each of the modules may exist alone physically, or two or more modules are integrated into one module. The integrated module can be realized in a hardware mode, and can also be realized in a software functional module mode.
Those skilled in the art will appreciate that all or part of the steps in the methods of the above embodiments may be implemented by associated hardware instructed by a program, which may be stored in a computer-readable storage medium, and the storage medium may include: read Only Memory (ROM), random Access Memory (RAM), magnetic or optical disks, and the like.
It will be understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiments may be implemented by hardware that is instructed to implement by a program, and the program may be stored in a computer-readable storage medium, where the above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.
While the invention has been described in detail with reference to specific embodiments thereof, it will be apparent to one skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention.

Claims (10)

1. An intracardiac shunt electrical impedance image analysis method based on saline angiography comprises the following steps:
acquiring a patient saline contrast impedance curve;
obtaining an image sequence diagram based on the saline contrast impedance curve, wherein the image sequence diagram comprises a right cardiac phase image sequence diagram and/or a left cardiac phase image sequence diagram;
extracting the central point of each image sequence diagram, and calculating the horizontal position relation between the central point and the central point of the first sequence diagram;
and outputting the classification result of the intra-cardiac shunt of the patient based on the position relation.
2. The method for intracardiac shunt electrical impedance image analysis based on saline contrast according to claim 1, characterized in that it comprises:
acquiring a patient saline contrast impedance curve;
obtaining a right heart phase image sequence chart based on the saline contrast impedance curve;
extracting the central point of each image sequence chart in the right heart phase image sequence chart, and calculating the horizontal position relation between the central point and the central point of a first sequence chart in the right heart phase image sequence chart;
and outputting whether the patient has the classification result of right-to-left intracardiac shunting or not based on whether the position relation is the left side or not.
3. The method for intracardiac shunt electrical impedance image analysis based on saline contrast according to claim 1, characterized in that it comprises:
acquiring a patient saline contrast impedance curve;
obtaining a left cardiac phase image sequence diagram based on the saline angiography impedance curve;
extracting the central point of each image sequence diagram in the left cardiac phase image sequence diagram, and calculating the horizontal position relation between the central point and the central point of a first sequence diagram in the left cardiac phase image sequence diagram;
and outputting whether the patient has a left-to-right intracardiac shunt classification result or not based on whether the position relation is right or not.
4. A saline contrast based intracardiac shunt electrical impedance image analysis method according to any of claims 1 to 3, characterized in that the determination of the central point of the sequence diagram is given by:
Figure 102047DEST_PATH_IMAGE001
wherein i belongs to the heart region H, x i Is the abscissa, f, of pixel point i i The optimal impedance descent slope, coH (t), of the least squares fit curve for the ith pixel point k ) Denotes t k The abscissa position of the heart center at time, t k Is the kth time window.
5. The intracardiac shunt electrical impedance image analysis method according to claim 4, wherein said heart region H is a pixel point impedance descent slope f>20%×f max Pixel point of (f) max The slope of the pixel point with the maximum slope is calculated by the following formula:
Figure 488029DEST_PATH_IMAGE002
wherein, delta
Figure 736608DEST_PATH_IMAGE003
Is composed of
Figure 584347DEST_PATH_IMAGE004
The relative impedance value of the time pixel point i; and e is the intercept.
6. The method for intracardiac shunt electrical impedance image analysis based on saline angiography of claim 1, wherein the right cardiac phase image sequence chart is obtained by reconstructing a right ventricular angiography impedance curve of a T0-T1 time period in the saline angiography impedance curve; and the left cardiac phase image sequence chart is obtained by reconstructing a left ventricle blood flow contrast impedance curve of a T2-T3 time period in the saline contrast impedance curve.
7. The method for intracardiac shunt electrical impedance image analysis based on saline contrast as claimed in claim 1, wherein said obtaining of image sequence plots based on said saline contrast impedance curves is obtaining of image sequence plots based on said saline contrast impedance curves using image reconstruction, said image reconstruction including any one or more of the following methods: projection reconstruction, light and shade recovery shape, stereoscopic vision reconstruction and laser ranging reconstruction.
8. A system for intra-cardiac shunt electrical impedance image analysis based on saline angiography, the system comprising:
an acquisition module for acquiring a patient saline contrast impedance curve;
an image reconstruction module for obtaining an image sequence diagram based on the patient saline angiography impedance curve, wherein the image sequence diagram comprises a right cardiac phase image sequence diagram and/or a left cardiac phase image sequence diagram;
the calculation module is used for extracting the central point of each image sequence diagram and calculating the horizontal position relation between the central point and the central point of the first sequence diagram;
and the output module outputs the intra-cardiac shunt classification result of the patient based on the position relation.
9. Intracardiac shunt electrical impedance image analysis device based on saline angiography, characterized in that it comprises: a memory and a processor;
the memory is to store program instructions;
the processor is configured to invoke program instructions, which when executed, perform the method of intra-cardiac shunt electrical impedance image analysis based on saline contrast as claimed in any one of claims 1 to 7.
10. A computer-readable storage medium, on which a computer program of a saline angiogram-based intracardiac shunt electrical impedance image analysis is stored, which, when being executed by a processor, carries out the saline angiogram-based intracardiac shunt electrical impedance image analysis method according to any one of claims 1-7.
CN202211532638.0A 2022-06-10 2022-12-02 Intracardiac shunt electrical impedance image analysis method and system based on saline angiography Active CN115530791B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211532638.0A CN115530791B (en) 2022-12-02 2022-12-02 Intracardiac shunt electrical impedance image analysis method and system based on saline angiography
PCT/CN2023/099037 WO2023237034A1 (en) 2022-06-10 2023-06-08 Image reconstruction and analysis method, system, and device based on saline radiography electrical impedance lung perfusion and cardiac imaging

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211532638.0A CN115530791B (en) 2022-12-02 2022-12-02 Intracardiac shunt electrical impedance image analysis method and system based on saline angiography

Publications (2)

Publication Number Publication Date
CN115530791A true CN115530791A (en) 2022-12-30
CN115530791B CN115530791B (en) 2023-03-10

Family

ID=84721595

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211532638.0A Active CN115530791B (en) 2022-06-10 2022-12-02 Intracardiac shunt electrical impedance image analysis method and system based on saline angiography

Country Status (1)

Country Link
CN (1) CN115530791B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023237034A1 (en) * 2022-06-10 2023-12-14 中国医学科学院北京协和医院 Image reconstruction and analysis method, system, and device based on saline radiography electrical impedance lung perfusion and cardiac imaging

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108550388A (en) * 2018-01-12 2018-09-18 深圳科亚医疗科技有限公司 The device and system of calculating vascular flow parameter based on angiography
US20190105104A1 (en) * 2001-12-04 2019-04-11 Atricure, Inc. Cardiac treatment devices and methods
CN111449657A (en) * 2020-04-15 2020-07-28 中国医学科学院北京协和医院 Bedside pulmonary ventilation-blood flow perfusion electrical impedance tomography method based on saline angiography
US11278261B1 (en) * 2017-06-01 2022-03-22 PFOmetrix, LLC Apparatus, system and method for the detection and quantification of conductance of right-to-left cardiac shunts
CN114519708A (en) * 2022-02-18 2022-05-20 济纶医工智能科技(南京)有限公司 CBIST imaging method and imaging system
CN114723845A (en) * 2022-06-10 2022-07-08 中国医学科学院北京协和医院 Saline contrast lung perfusion image reconstruction method, system and equipment with SPECT as standard

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190105104A1 (en) * 2001-12-04 2019-04-11 Atricure, Inc. Cardiac treatment devices and methods
US11278261B1 (en) * 2017-06-01 2022-03-22 PFOmetrix, LLC Apparatus, system and method for the detection and quantification of conductance of right-to-left cardiac shunts
CN108550388A (en) * 2018-01-12 2018-09-18 深圳科亚医疗科技有限公司 The device and system of calculating vascular flow parameter based on angiography
CN111449657A (en) * 2020-04-15 2020-07-28 中国医学科学院北京协和医院 Bedside pulmonary ventilation-blood flow perfusion electrical impedance tomography method based on saline angiography
CN114519708A (en) * 2022-02-18 2022-05-20 济纶医工智能科技(南京)有限公司 CBIST imaging method and imaging system
CN114723845A (en) * 2022-06-10 2022-07-08 中国医学科学院北京协和医院 Saline contrast lung perfusion image reconstruction method, system and equipment with SPECT as standard

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023237034A1 (en) * 2022-06-10 2023-12-14 中国医学科学院北京协和医院 Image reconstruction and analysis method, system, and device based on saline radiography electrical impedance lung perfusion and cardiac imaging

Also Published As

Publication number Publication date
CN115530791B (en) 2023-03-10

Similar Documents

Publication Publication Date Title
Secinaro et al. Cardiovascular magnetic resonance findings in repaired anomalous left coronary artery to pulmonary artery connection (ALCAPA)
Ferrari et al. Ultrafast three-dimensional contrast-enhanced magnetic resonance angiography and imaging in the diagnosis of partial anomalous pulmonary venous drainage
KR102402628B1 (en) Methods and systems for modeling the human heart and atria
ES2378208T3 (en) CARD�? ACOS DIAGNOSTICS THROUGH THE USE OF CARD� MAGNETIC RESONANCE ACA OF WALL MOVEMENT AND PERFUSION AND CARD�? ACO DIAGNOSTIC SYSTEMS.
CN111449657B (en) Image monitoring system and pulmonary embolism diagnosis system
CN101443812A (en) Method for determining and displaying at least one piece of information on a target volume
CN115530791B (en) Intracardiac shunt electrical impedance image analysis method and system based on saline angiography
CN115587971B (en) Organism reaction and hemodynamic monitoring method and system based on heart ultrasonic segment activity
WO2023237034A1 (en) Image reconstruction and analysis method, system, and device based on saline radiography electrical impedance lung perfusion and cardiac imaging
CN115530792B (en) Right heart failure image analysis method, system and equipment based on saline angiography
US20100135548A1 (en) Medical Imaging System
CN114723844B (en) Method, system and equipment for reconstructing pulsatile perfusion image corrected by saline contrast
Wajdan et al. Automatic detection of valve events by epicardial accelerometer allows estimation of the left ventricular pressure trace and pressure–displacement loop area
CN115035208A (en) Lung perfusion and region V/Q non-invasive imaging method, system and equipment
Binney The function of the heart is historically contingent
CN114129146A (en) Device, method and medium for analyzing impedance data of pulmonary perfusion
Donovan et al. Sinus venosus atrial septal defect as a cause of palpitations and dyspnea in an adult: a diagnostic imaging challenge
Krishnamurthy Pediatric cardiac MRI: anatomy and function
Crean Technical considerations for ACHD imaging
Muscogiuri et al. Multi-modality imaging approach in a challenging case of surgically corrected partial anomalous pulmonary venous return and atrial tachycardia treated with radiofrequency ablation
Hartmann et al. Validity of acoustic quantification colour kinesis for detection of left ventricular regional wall motion abnormalities: a transoesophageal echocardiographic study
Dual et al. Acute changes in preload and the QRS amplitude in advanced heart failure patients
JP2020171475A (en) Dynamic image analysis apparatus, dynamic image analysis method, and program
KR102607446B1 (en) Heart failure diagnosis method by measuring ventricular ejection fraction
Trimech et al. Scimitar syndrome with bicuspid aortic valve. a case report of cross-sectional non-invasive imaging allowing a complete anatomical and functional assessment

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20221230

Assignee: Dianqi biomedical technology (Suzhou) Co.,Ltd.

Assignor: PEKING UNION MEDICAL COLLEGE Hospital

Contract record no.: X2024980003899

Denomination of invention: A method and system for analyzing intracardiac shunt impedance images based on saline angiography

Granted publication date: 20230310

License type: Common License

Record date: 20240402