US20210085284A1 - System and method for electromechanical activation of arrhythmias - Google Patents
System and method for electromechanical activation of arrhythmias Download PDFInfo
- Publication number
- US20210085284A1 US20210085284A1 US16/871,355 US202016871355A US2021085284A1 US 20210085284 A1 US20210085284 A1 US 20210085284A1 US 202016871355 A US202016871355 A US 202016871355A US 2021085284 A1 US2021085284 A1 US 2021085284A1
- Authority
- US
- United States
- Prior art keywords
- heart
- wave
- electromechanical
- propagation
- map
- 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.)
- Abandoned
Links
- 238000000034 method Methods 0.000 title claims abstract description 85
- 206010003119 arrhythmia Diseases 0.000 title claims description 47
- 230000006793 arrhythmia Effects 0.000 title claims description 35
- 230000004913 activation Effects 0.000 title description 39
- 238000003384 imaging method Methods 0.000 claims description 43
- 238000002604 ultrasonography Methods 0.000 claims description 35
- 230000033764 rhythmic process Effects 0.000 claims description 26
- 230000008602 contraction Effects 0.000 claims description 25
- 238000013507 mapping Methods 0.000 claims description 23
- 230000002861 ventricular Effects 0.000 claims description 14
- 230000001746 atrial effect Effects 0.000 claims description 13
- 230000000694 effects Effects 0.000 claims description 11
- 238000006073 displacement reaction Methods 0.000 abstract description 83
- 230000033001 locomotion Effects 0.000 abstract description 33
- 238000004458 analytical method Methods 0.000 abstract description 14
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000001939 inductive effect Effects 0.000 abstract 1
- 230000000747 cardiac effect Effects 0.000 description 36
- 210000004165 myocardium Anatomy 0.000 description 27
- 210000002837 heart atrium Anatomy 0.000 description 25
- 210000001519 tissue Anatomy 0.000 description 25
- 230000006870 function Effects 0.000 description 19
- 206010003662 Atrial flutter Diseases 0.000 description 16
- 210000005246 left atrium Anatomy 0.000 description 16
- 206010003658 Atrial Fibrillation Diseases 0.000 description 13
- 230000001902 propagating effect Effects 0.000 description 13
- 210000005240 left ventricle Anatomy 0.000 description 11
- 239000000523 sample Substances 0.000 description 11
- 210000005245 right atrium Anatomy 0.000 description 10
- 230000002123 temporal effect Effects 0.000 description 9
- 230000001052 transient effect Effects 0.000 description 9
- 238000011282 treatment Methods 0.000 description 9
- 241000699666 Mus <mouse, genus> Species 0.000 description 8
- 208000009729 Ventricular Premature Complexes Diseases 0.000 description 8
- 238000001514 detection method Methods 0.000 description 8
- 208000028867 ischemia Diseases 0.000 description 8
- 238000005259 measurement Methods 0.000 description 8
- 238000012545 processing Methods 0.000 description 8
- 210000000709 aorta Anatomy 0.000 description 7
- 230000002600 fibrillogenic effect Effects 0.000 description 7
- 238000012285 ultrasound imaging Methods 0.000 description 7
- 206010003668 atrial tachycardia Diseases 0.000 description 6
- 230000002596 correlated effect Effects 0.000 description 6
- 230000000875 corresponding effect Effects 0.000 description 6
- 239000000835 fiber Substances 0.000 description 6
- 210000001367 artery Anatomy 0.000 description 5
- 238000001914 filtration Methods 0.000 description 5
- 230000006794 tachycardia Effects 0.000 description 5
- 206010003130 Arrhythmia supraventricular Diseases 0.000 description 4
- 241000699670 Mus sp. Species 0.000 description 4
- 208000001871 Tachycardia Diseases 0.000 description 4
- 238000002679 ablation Methods 0.000 description 4
- 206010061592 cardiac fibrillation Diseases 0.000 description 4
- 210000000038 chest Anatomy 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000002107 myocardial effect Effects 0.000 description 4
- 208000010125 myocardial infarction Diseases 0.000 description 4
- 208000031225 myocardial ischemia Diseases 0.000 description 4
- 230000000644 propagated effect Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 208000024891 symptom Diseases 0.000 description 4
- 241000282412 Homo Species 0.000 description 3
- 230000006399 behavior Effects 0.000 description 3
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013467 fragmentation Methods 0.000 description 3
- 238000006062 fragmentation reaction Methods 0.000 description 3
- 238000001727 in vivo Methods 0.000 description 3
- 238000012417 linear regression Methods 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 238000001356 surgical procedure Methods 0.000 description 3
- 206010047302 ventricular tachycardia Diseases 0.000 description 3
- 208000024172 Cardiovascular disease Diseases 0.000 description 2
- 206010061216 Infarction Diseases 0.000 description 2
- 206010028813 Nausea Diseases 0.000 description 2
- 208000002193 Pain Diseases 0.000 description 2
- 238000011298 ablation treatment Methods 0.000 description 2
- 230000036982 action potential Effects 0.000 description 2
- 210000000702 aorta abdominal Anatomy 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000009534 blood test Methods 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000004087 circulation Effects 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- CVSVTCORWBXHQV-UHFFFAOYSA-N creatine Chemical compound NC(=[NH2+])N(C)CC([O-])=O CVSVTCORWBXHQV-UHFFFAOYSA-N 0.000 description 2
- 230000002526 effect on cardiovascular system Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 208000019622 heart disease Diseases 0.000 description 2
- 230000004217 heart function Effects 0.000 description 2
- 230000001976 improved effect Effects 0.000 description 2
- 230000007574 infarction Effects 0.000 description 2
- 230000000302 ischemic effect Effects 0.000 description 2
- 230000010358 mechanical oscillation Effects 0.000 description 2
- 230000008693 nausea Effects 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000036407 pain Effects 0.000 description 2
- 210000003540 papillary muscle Anatomy 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 210000005241 right ventricle Anatomy 0.000 description 2
- 210000001013 sinoatrial node Anatomy 0.000 description 2
- 230000035900 sweating Effects 0.000 description 2
- 230000001960 triggered effect Effects 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 208000004998 Abdominal Pain Diseases 0.000 description 1
- 241000282465 Canis Species 0.000 description 1
- 208000031229 Cardiomyopathies Diseases 0.000 description 1
- 206010008479 Chest Pain Diseases 0.000 description 1
- 208000000059 Dyspnea Diseases 0.000 description 1
- 206010013975 Dyspnoeas Diseases 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 208000008454 Hyperhidrosis Diseases 0.000 description 1
- 208000002682 Hyperkalemia Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 102000036675 Myoglobin Human genes 0.000 description 1
- 108010062374 Myoglobin Proteins 0.000 description 1
- 102000013394 Troponin I Human genes 0.000 description 1
- 108010065729 Troponin I Proteins 0.000 description 1
- 206010047281 Ventricular arrhythmia Diseases 0.000 description 1
- 210000001015 abdomen Anatomy 0.000 description 1
- 210000003815 abdominal wall Anatomy 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000005856 abnormality Effects 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 238000011366 aggressive therapy Methods 0.000 description 1
- 239000000935 antidepressant agent Substances 0.000 description 1
- 229940005513 antidepressants Drugs 0.000 description 1
- 210000002376 aorta thoracic Anatomy 0.000 description 1
- 210000001765 aortic valve Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000036765 blood level Effects 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 238000009125 cardiac resynchronization therapy Methods 0.000 description 1
- 230000005779 cell damage Effects 0.000 description 1
- 208000037887 cell injury Diseases 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 210000004351 coronary vessel Anatomy 0.000 description 1
- 229960003624 creatine Drugs 0.000 description 1
- 239000006046 creatine Substances 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000012631 diagnostic technique Methods 0.000 description 1
- 230000003205 diastolic effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000002091 elastography Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000010247 heart contraction Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011503 in vivo imaging Methods 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 238000007914 intraventricular administration Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 210000004115 mitral valve Anatomy 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 230000008065 myocardial cell damage Effects 0.000 description 1
- 210000001087 myotubule Anatomy 0.000 description 1
- 230000001314 paroxysmal effect Effects 0.000 description 1
- 239000008177 pharmaceutical agent Substances 0.000 description 1
- 230000000144 pharmacologic effect Effects 0.000 description 1
- 239000002574 poison Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- HYTYHTSMCRDHIM-UHFFFAOYSA-M potassium;2-sulfanylacetate Chemical compound [K+].[O-]C(=O)CS HYTYHTSMCRDHIM-UHFFFAOYSA-M 0.000 description 1
- 238000007674 radiofrequency ablation Methods 0.000 description 1
- 230000007115 recruitment Effects 0.000 description 1
- 230000010410 reperfusion Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 210000002966 serum Anatomy 0.000 description 1
- 208000013220 shortness of breath Diseases 0.000 description 1
- 230000003319 supportive effect Effects 0.000 description 1
- 238000011477 surgical intervention Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 238000012731 temporal analysis Methods 0.000 description 1
- 230000008719 thickening Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
- YFDSDPIBEUFTMI-UHFFFAOYSA-N tribromoethanol Chemical compound OCC(Br)(Br)Br YFDSDPIBEUFTMI-UHFFFAOYSA-N 0.000 description 1
- 229950004616 tribromoethanol Drugs 0.000 description 1
- 210000000591 tricuspid valve Anatomy 0.000 description 1
- 230000035899 viability Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
- A61B8/0883—Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of the heart
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/485—Diagnostic techniques involving measuring strain or elastic properties
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
- A61B8/5223—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/895—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum
- G01S15/8956—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques characterised by the transmitted frequency spectrum using frequencies at or above 20 MHz
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8977—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using special techniques for image reconstruction, e.g. FFT, geometrical transformations, spatial deconvolution, time deconvolution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52036—Details of receivers using analysis of echo signal for target characterisation
- G01S7/52042—Details of receivers using analysis of echo signal for target characterisation determining elastic properties of the propagation medium or of the reflective target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52085—Details related to the ultrasound signal acquisition, e.g. scan sequences
- G01S7/52087—Details related to the ultrasound signal acquisition, e.g. scan sequences using synchronization techniques
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
- A61B8/461—Displaying means of special interest
- A61B8/463—Displaying means of special interest characterised by displaying multiple images or images and diagnostic data on one display
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8909—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
- G01S15/8915—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
- G01S15/8925—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being a two-dimensional transducer configuration, i.e. matrix or orthogonal linear arrays
Definitions
- This present disclosed subject matter relates to techniques for imaging and detecting the propagation of mechanical waves within a body structure of a patient.
- Certain medical conditions can be difficult to establish in their early stages when treatment is most effective.
- Patients suffering from myocardial ischemia can present to an emergency room or acute care facility with typical cardiac symptoms such as chest pain, described as tightness, pressure, or squeezing, but some patients can have other symptoms such as arm or chin pain, nausea, sweating, or abdominal pain.
- Certain techniques such as electrocardiogram often provide inconclusive findings regarding ischemia, and sometimes can even be unable to identify situations in which ischemia has progressed to cell damage and myocardial infarction (MI). Other techniques are available for diagnosing infarction relative to its predecessor, ischemia.
- CK-MB creatine kinase-MB
- Other serum markers include troponin I, and to a lesser extent, myoglobin.
- the blood levels of certain such compounds can take several hours to rise, so that diagnosis of MI can be delayed. Reliance on blood tests alone can result in a significant loss of time when early aggressive therapy is warranted.
- Certain less invasive diagnostic techniques have become available through the observation of mechanical properties of tissue via imaging techniques.
- Such evaluation of the function of the heart, cardiovascular tissue, or other body structures can be based on the mechanical interpretation of the movement of the these structures, such as, for example, the active contractions and passive relaxation of the myocardium.
- the evaluation of the heart function can be based on a single mechanical interpretation of myocardial deformation.
- the deformations of the myocardium can be quantified over a complete cardiac cycle in order to provide some information on the myocardial viability.
- Certain low frequency mechanical vibrations in the heart are known in human patients. Certain ultrasound techniques can be used to obtain pulsive mechanical vibrations around end-systole and end-diastole in the frequency range of 25 to 100 Hz.
- Atrial arrhythmias are a known and can cause of morbidity and mortality. Certain mechanical factors, such as chamber size and wall tension, can affect the onset and perpetuation of atrial arrhythmia. Certain echocardiographic measurements can also be used to characterize atrial arrhythmias. Yet, systems and techniques to analyze the 2-D spatio-temporal evolution of the local deformations of the atria during e.g., focal tachycardia, flutter, and fibrillation, would be beneficial.
- the present disclosure provides elasticity imaging techniques to evaluate mechanical wave propagation, and provide an estimation of electrical propagation in a noninvasive manner.
- the disclosed subject matter provides systems and methods for detecting wave propagation within the tissue of a patient in a series of image frames representing movement of such tissue of the body structure.
- Image data is acquired comprising a series of image frames corresponding to the movement of the tissue.
- the tissue can be the wave propagation in the myocardium.
- the movement of body tissue can be wave propagation in the arteries or the aorta.
- a correlation calculation can be performed on the image frames to generate a matrix with the location of correlation maxima representing the relative displacement between the first and second image frames, also referred to as a displacement map.
- a video can be generated comprising a series of displacement maps.
- the parameters of movement of the cardiac structure can be detected, such as velocity, attenuation, frequency, etc.
- the wave can be a shear wave, representative of the electrical wave propagation within the body structure.
- Image information of a heart of a subject can be obtained using an imaging device.
- a strain map of the heart can be generated from the image information.
- Occurrences of a first electromechanical event of the heart and a second electromechanical event can be determined from the strain map.
- a spatio-temporal map of atrial and ventricular mechanics of the heart can be generated by tracking the onset of the first and second event for each pixel of a heart wall of the subject identified from the image information.
- a representative mechanical cycle associated with a contraction of the heart can be identified using the spatio-temporal map.
- electromechanical activation mapping can characterize propagation patterns of electromechanical strains during focal and reentrant arrhythmias of the heart.
- regions in which the mechanical and electrical activities are decoupled can be identified by mapping the electromechanical activity of the heart.
- the spatio-temporal map can be generated by obtaining isochrones strongly correlated to electrical isochrones by tracking a propagation front of an end-diastole electromechanical activation of the heart.
- a type of cardiac arrhythmia present in the heart can be determined from the information collected in the ultrasound scan of the heart.
- an onset of ventricular contraction can be determined by identifying a first zero-crossing of an incremental strains occurring after an onset of a P-wave on an electrocardiogram.
- a high-resolution Fourier transform can be performed using a generalized Goertzel algorithm to interpolate strain signals in Fourier space for each individual pixel in an atria of the heart.
- a peak mechanical cycle length (MCL) map can be generated by selecting a MCL having a highest amplitude within the physiologically-relevant time range for each pixel of the ultrasound scan of the heart, such that the MCL map identifies, for each pixel of the atria, which cycle length is most present in a Fourier spectrum of cycle lengths.
- the cycle length best representing an atrial contraction of the heart can be determined.
- a phase corresponding to the determined cycle length can be determined to map a propagation of a mechanical oscillation of the heart at the determined cycle length.
- the ultrasound scan further can be performed by emitting a circular ultrasonic wave to instruct an ultrasound apparatus to perform a motion estimation sequence.
- a B-mode acquisition can be performed to capture heart anatomy of the heart.
- a plurality of beams can be generated to reconstruct frames from the motion estimation sequence using a delay-and-sum algorithm with a reconstructed sampling frequency.
- FIG. 1 is a diagram illustrating the system in accordance with the present disclosure.
- FIG. 2 is a diagram illustrating exemplary stages in a method in accordance with the present disclosure.
- FIG. 3 is a diagram illustrating a technique for measuring movement of structures within an image in accordance with the present disclosure.
- FIG. 4 is a chart representing the velocity of structures within an image in accordance with the present disclosure.
- FIG. 5 illustrates a technique of detecting wave propagation in accordance with a further embodiment of the present disclosure.
- FIG. 6 illustrates a technique of detecting wave propagation in accordance with yet another embodiment of the present disclosure.
- FIGS. 7( a )-11( b ) are images illustrating the propagation of a wave within a body structure in accordance with an exemplary embodiment of the present disclosure.
- FIGS. 12( a )-18( b ) are images illustrating the propagation of a wave within a body structure in accordance with another exemplary embodiment of the present disclosure.
- FIG. 19 is an ultrasound image of a mouse left ventricle in a parasternal long-axis view.
- FIG. 20( a ) is an axial displacement map overlaid to the grayscale B-mode image of the left ventricle during systole in accordance with the present disclosure.
- FIG. 20( b ) is an axial displacement map overlaid to the grayscale B-mode image of the left ventricle during diastole (relaxation phase) in accordance with the present disclosure.
- FIG. 20( c ) is an ECG indicating the time of the acquisition during the cardiac cycle of FIG. 20( a ) in accordance with the present disclosure.
- FIG. 20( d ) is an ECG indicating the time of the acquisition during the cardiac cycle of FIG. 20( b ) in accordance with the present disclosure.
- FIG. 21( a ) is a time plot illustrating the temporal variation of the axial displacements estimated on one central RF-line as line plotted on FIG. 20( b ) in accordance with the present disclosure.
- FIG. 21( b ) illustrates the frequency content of the displacement variation in the septum at the depth of 12.5 mm plotted as a function of time in accordance with the present disclosure.
- FIG. 21( c ) is a time plot illustrating the temporal variation of the axial displacements after bandpass filtering of the plot illustrated in FIG. 21( a ) showing the transient and high frequency components in accordance with the present disclosure.
- FIG. 21( d ) illustrates the ECG signal acquired simultaneously with the data illustrated in FIGS. 21( a )-( c ) in accordance with the present disclosure.
- FIGS. 22( a )-( d ) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the left ventricle around end-systole taken every 0.6 ms showing the propagation of a first mechanical wave front in the septum in accordance with the present disclosure.
- the arrows indicate the progression of the wave front in the septum.
- FIGS. 22( e )-( f ) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the left ventricle around end-systole taken every 0.6 ms showing the propagation of a second mechanical wave front in the septum in accordance with the present disclosure.
- the arrows indicate the progression of the wave front in the septum.
- FIGS. 22( g )-( l ) illustrate the ECG signal plotted below each respective image of FIGS. 22( a )-( f ) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure.
- FIG. 23 is a plot illustrating the distance of propagation as a function of the phase of the end-systolic wave at three frequencies in accordance with the present disclosure.
- FIGS. 24( a )-( f ) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the left ventricle around the beginning of systole taken every 2.8 ms, showing the propagation of a strong mechanical wave in the posterior wall in accordance with the present disclosure.
- the arrows indicate the progression of the wave front in the posterior wall.
- FIGS. 24( g )-( l ) illustrate the ECG signal plotted below each respective image of FIGS. 24( a )-( f ) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure.
- FIG. 25 is a plot illustrating the distance of propagation as a function of the phase of the wave at the frequency of 80 Hz during the beginning of systole transient motion in accordance with the present disclosure.
- FIGS. 26( a )-( e ) illustrate a sequence of axial displacement maps overlaid to the grayscale image (0.12 ms between successive frames) indicating an electromechanical wave propagating in the posterior wall of the mouse from the apex towards the base during pacing in the right atrium close to the sinoatrial node in accordance with the present disclosure.
- FIGS. 26( f )-( j ) illustrate the ECG signal plotted below each respective image of FIGS. 26( a )-( e ) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure.
- FIGS. 27( a )-( e ) illustrate a sequence of axial displacement maps overlaid to the grayscale image (0.07 ms between successive frames) indicating an electromechanical wave propagating in the posterior wall of the mouse from the base towards the apex during pacing in the right ventricle close to the base in accordance with the present disclosure.
- FIGS. 27( f )-( j ) illustrate the ECG signal plotted below each respective image of FIGS. 27( a )-( e ) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure.
- FIGS. 28( a )-( f ) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the aorta taken every 0.7 ms. Sequence of images showing the propagation of the pulse wave in the aorta. The arrows indicate the progression of the wave front in the aorta
- FIGS. 28( g )-( l ) illustrate the ECG signal plotted below each respective image of FIGS. 26( a )-( f ) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure.
- FIG. 29 illustrates the distance of propagation as a function of the phase of pulse wave at the frequency of 80 Hz.
- the slope of the curve gives the pulse wave velocity in accordance with the present disclosure.
- FIGS. 30( a )-( c ) illustrates examples of high temporal resolution strains during different types of arrhythmia in accordance with the present disclosure.
- FIGS. 31( a )-( d ) illustrates examples of propagating electromechanical activation in atria of normal subjects and a patient undergoing focal atrial tachycardia in accordance with the present disclosure.
- FIGS. 32( a )-( d ) illustrates analysis of reentrant arrhythmias using a single-frequency atrial flutter case using a peak cycle length map, cycle length histogram, a phase map, and electrogram in accordance with the present disclosure.
- FIG. 33( a ) illustrates peak cycle length maps in two atrial flutter patients in accordance with the present disclosure.
- FIG. 33( b ) illustrates a correlation between mechanical cycle length (MCL) and electrical cycle length in five atrial flutter patients in accordance with the present disclosure.
- FIG. 33( c ) illustrates a peak cycle length map during atrial fibrillation in accordance with the present disclosure.
- FIG. 33( d ) illustrates a histogram depicting spatial fragmentation of peak cycle length during atrial fibrillation in accordance with the present disclosure.
- the system and methods described herein can be useful for analyzing data obtained by an image generating device, such as an ultrasound transducer.
- the systems and methods can also be useful for measuring mechanical properties and estimating the electrical characteristics of a body tissue structure or organ, such as, for example, the myocardium or the aorta.
- the disclosed subject matter can be used in connection with imaging and characterizing the propagation of electromechanical waves in the heart.
- electrical waves propagate in the myocardium in order to induce its contraction. Contraction of the myocardial fibers results in a strong mechanical wave, which, since it results from the coupling of the electrical excitation and the mechanical properties of the myocardium, is referred to herein as an “electromechanical wave.”
- the speed of this wave is a function of the electrical and mechanical properties of the myocardium, and, according to the present disclosure, can be used to detect changes in these properties to diagnose heart diseases.
- System 100 can include an image detection device, such as ultrasound probe 102 , which is used to create images of the heart H or other organ or structure of the patient P.
- the image detection device does not induce discernible vibration in the body structure, and merely detects pre-existing motion.
- the signals detected by the probe 102 can be transferred to an ultrasound scanner 104 .
- the exemplary embodiments described herein are designed to work with conventional ultrasound scanners. For example, commercial portable scanners, such as Terason 2000, high frequency scanners, such as Visualsonics Vevo 770, and routinely used clinical scanners, such as GE System Five or GE Vivid Five or Seven, are useful for image acquisition in accordance with the exemplary embodiments.
- the raw data produced by the scanner 104 can be transferred to a computer 106 having a CPU 108 for processing the data.
- the computer and CPU would be Dell PC with a 2 GHz processor. It is understood that the computer and CPU can also be integrated with the ultrasound scanner 104 .
- Also useful in the system would be storage such as disk drives, for storing data on input files 110 and for writing output onto output files 112 .
- input files 110 can include information such as thresholds.
- Output files 112 can include the displacement maps, videos of myocardium displacements, or computed data, such as electromechanical wave properties. It is understood that a preprogrammed chip can be used to execute the algorithms described herein.
- an output device such as monitor 114
- an input device such as keyboard 116
- the methods described herein are particularly useful for imaging the propagation of electromechanical waves in the heart.
- a method for detecting the properties of the electromechanical wave are described herein and represented in FIG. 2 .
- raw imaging data of the body structure is acquired by image acquisition equipment such as the ultrasound probe 102 and scanner 104 .
- a set of N frames of raw ultrasound data of the heart is acquired during a cardiac cycle at high frame rate, e.g., higher than 100 fps, although frame rates of about 56 fps and 170 fps, etc., yield useful results ( 202 ).
- the selected frame rate should be commensurate with the speed of the propagation of the movement, such as the wave, being studied.
- the electrocardiogram (EKG) can also be recorded.
- the raw data can be digitized and stored in real-time in the scanner memory.
- the data can be transferred to a computer for processing ( 204 ).
- the transfer can occur using a protocol such as Ethernet TCP IP.
- the computer can be integrated with the scanner 104 .
- the raw data received from the image acquisition equipment is processed.
- the data processing computes an estimation of the displacement of particular objects in the images, such as the myocardium, between consecutive frames. Typically this processing occurs off-line; however, it is understood that this procedure can occur sequentially subsequent to receiving two consecutive frames.
- axial displacements in the direction of the transducer
- Lateral, or elevational, displacements perpendicular to the transducer
- N ⁇ 1 displacement 2D maps are computed through the correlation of two consecutive frames i and i+1 (1 ⁇ i ⁇ N ⁇ 1). Each frame is represented by a matrix of pixel values.
- the displacement maps provide an indication of the local axial movements between frames. Estimation of the axial displacements from the two consecutive frames is performed using a speckle tracking algorithm. In this algorithm, the time-shifts in the backscattered signals are determined between two consecutive frames through cross-correlation of small sliding windows over the entire ultrasound image. For each window, the signal of the frame i and the frame i+1 are cross-correlated. The maximum of the correlation coefficient gives an estimation of the time-shift between the two signals.
- This time-shift can be converted to a spatial displacement by assuming a constant speed of sound for the tissue.
- This technique can detect displacements on the order of 10 ⁇ m. Using small correlation windows of 7.5 mm, the resolution of the displacement maps is in the millimeter range.
- the cross-correlation algorithm suitable for estimating displacement between consecutive image frames is described in U.S. Provisional Patent Application No. 60/619,247, filed Oct. 15, 2004, which is incorporated by reference herein.
- a Matlab program Multiframe is used to compute the displacement maps for the complete sequence of frames obtained at 202 , above. Multiframe calls the Matlab routine FunCalculDispl to compute the displacements for the sequence of frames.
- Correlation2D.cpp is a C program that computes the displacement map between consecutive frames.
- Correlation2D.cpp uses small sliding windows to find the displacement which maximizes the correlation coefficient for each part of the image.
- auto-correlation calculations or coherence calculations can be performed.
- Two optional threshold procedures can be executed in the procedure 200 .
- a threshold can be applied on the energy of the signal, in order to remove the noise that is below a predetermined signal-to-noise ratio. Low energy ultrasound signals (e.g., noise in the cavity of the heart) can be removed from the displacement map according to this method.
- a threshold can be applied on the correlation coefficient to remove erroneous estimates in the displacements.
- the noise threshold and correlation-coefficient threshold can be implemented within the routine Correlation2D.cpp.
- the levels of the thresholds are determined experimentally and can be stored in an input data file 110 for processing on the CPU 108 .
- Procedures 206 , 208 and 210 are illustrated sequentially; however, it is understood that they can occur simultaneously or any other order to appropriately process the data. Moreover, one or more of these procedures can be omitted from the process described herein.
- a video of the sequence of N ⁇ 1 displacement maps can be assembled to create a video of the displacements of the body structure or tissue ( 212 ).
- a video of the myocardium displacements is created by this technique.
- the video of the displacement map of the myocardium will depict the propagation of the electromechanical wave.
- an observation and tracking of the wave propagation can be done manually, it can be difficult to discern the wave by visual observation and thus make accurate measurements.
- wave tracking can be performed by an algorithm, such as TrackPositionWave, a Matlab program which locates the position of the wave front by performing a zero-crossing calculation on consecutive displacement maps.
- the parameters of the electromechanical wave can be analyzed at 216 .
- the velocity of the electromechanical wave can be computed as a function of its position in the myocardium.
- the wall of the myocardium is approximated as circular with a radius R, and the origin of the spherical coordinate system was chosen at the center of the cavity.
- the wavefront of the electromechanical wave was then tracked by its angular coordinate ⁇ .
- the Matlab function Overlay can be used to compute the transformation of the raw image into polar coordinates.
- This routine can also display the displacement map superimposed on the ultrasound grayscale data.
- the speed of the electromechanical wave is shown on the FIG. 4 as a function of the angular position.
- the ultrasound imaging method described herein has the advantage of being completely non-invasive.
- the system described herein can be implemented in real-time on commercial scanners. It has been shown that the electrical conductivity is transversely isotropic with respect to fiber direction, with a longitudinal velocity of about 0.6 m/s and a transverse velocity of about 0.2 m/s (Roth, B. J. (2000), Influence of a perfusing bath on the foot of the cardiac action potential. Circulation Research 86, E19-E22; Spach, M. S., Heidlage, J. F., Dolber, P. C., and Barr, R. C. (1998), Extracellular discontinuities in cardiac muscle—Evidence for capillary effects on the action potential foot.
- FIG. 5 illustrates a transducer setup for 2D imaging of the longitudinal waves
- FIG. 6 illustrates a transducer setup for 2D imaging of transverse waves.
- the mechanical component of the electromechanical wave is related to the viscoelastic properties of the soft tissue.
- the elastic properties of the myocardium have been widely investigated.
- the stiffness of the myocardium has been shown to increase during ischemia and recovers after reperfusion.
- early detection of cardiovascular diseases such as ischemia and infarction can be strongly improved through non-invasive characterization of the local myocardial elasticity.
- V S ⁇ ⁇ ( 1 )
- a system can be implemented to provide early detection of ischemia through the measure of the velocity of the mechanical wave.
- the myocardium has also anisotropic mechanical properties and can be considered as a transverse isotropic medium.
- two shear waves of different velocities can propagate in the myocardium.
- Fast mechanical (shear) waves propagate in the direction of the fibers
- slow mechanical (shear) waves propagate in the direction perpendicular to the fibers.
- the measure of the two wave velocities can be achieved by using 3D Ultrasound imaging systems or multiple acquisitions of 2D images with a rotation of the transducer (see FIGS. 5-6 ).
- the wave velocities are related to two elastic constants, ⁇ // the shear modulus in direction of the fibers and ⁇ ⁇ the cross-fiber shear modulus.
- V // ⁇ / / ⁇ ( 2 )
- V ⁇ ⁇ ⁇ ⁇ ( 3 )
- the systems and methods described herein can potentially have different applications in the field of early detection of cardiovascular diseases and cardiac imaging.
- the measure of the electrical excitation propagation is of high interest in cardiology for early detection of heart diseases but also for pacing the heart when heartbeat is too slow or irregular.
- the purpose of an artificial pacemaker is to stimulate the heart when either the heart's natural pacemaker is not fast enough or if there are blocks in the electrical conduction system preventing the propagation of electrical impulses.
- the electrical propagation must be determined accurately.
- In vivo imaging of the electrical propagation in the heart can require implanting an electrode matrix (up to 500 electrodes) to measure extracellular potentials at the surface of the heart. This invasive and potentially precarious surgical procedure cannot be performed on human patients for diagnostic purposes.
- the present disclosure provides a means for determining electrical propagation in the myocardium of a subject, in the context of achieving an effective position of a pacemaker in the subject.
- Other methods known in the art involve optically-based techniques which also require invasive procedures, such as open-chest surgery.
- the present disclosure can be further used to create images and thereby detect myocardial ischemia in a subject either having symptoms (e.g., chest, arm, or chin pain, nausea, shortness of breath, and/or sweating) or a subject subjectively lacking such symptoms (e.g., “silent ischemia”), whereby a finding of increased electromechanical wave velocity (relative to control values) in a region of the myocardium of a subject is consistent with and supportive of a diagnosis of myocardial ischemia in that region.
- symptoms e.g., chest, arm, or chin pain, nausea, shortness of breath, and/or sweating
- a subject subjectively lacking such symptoms e.g., “silent ischemia”
- the present disclosure can also be used to diagnose, or assist in surgical intervention in, (i) conduction disturbances, such as re-entry phenomena, or associated with pharmaceutical agents, such as antidepressants or hyperkalemia, (ii) arrhythmias and dysrhythmias (e.g., surgical treatment of ventricular dysrhythmias, diagnosis of low-amplitude atrial fibrillation); and (iii) tissue abnormalities associated with cardiomyopathies or trauma, etc.
- conduction disturbances such as re-entry phenomena
- pharmaceutical agents such as antidepressants or hyperkalemia
- arrhythmias and dysrhythmias e.g., surgical treatment of ventricular dysrhythmias, diagnosis of low-amplitude atrial fibrillation
- tissue abnormalities associated with cardiomyopathies or trauma etc.
- the procedure described hereinabove was performed in an anesthetized open-chested dog.
- the transducer was placed on the anterior wall of the left ventricle of the heart, to obtain a short axis view.
- a sequence of three cardiac cycles was acquired during the experiment, with a frame rate of 56 fps.
- the 2D displacement maps were estimated using the cross-correlation method (window size: 5 mm, 90% overlap).
- the axial displacements were processed for the different sequences.
- two electromechanical waves were clearly detected, propagating in the posterior wall of the left ventricular, from the septum (left side of the images) to the lateral wall (right side). The propagation of the mechanical wave corresponds to the electrical activity shown on an associated EKG.
- the first electromechanical wave is found at the end-diastolic phase of the cardiac cycle (which corresponds to the beginning of the contraction).
- FIGS. 5( a ) through 9( a ) show five consecutive frames of the propagation of the wave.
- the location of the electromechanical wavefront is indicated by arrow W in each of FIGS. 7( a )-11( a ) .
- the displacements maps are overlaid to the grayscale ultrasound images (see, FIGS. 7( b )-11( b ) ). Blue displacements are in the direction of the transducer (top of the image), and red displacements are in the opposite direction.
- the contraction of the myocardium starts on the left side (septum) and propagates to the right side of the image.
- the blue region appears on the left side of the images (behind the wavefront), and the red region appears on the right side of the image (in front of the wavefront).
- the maximum displacements shown are 75 ⁇ m (dark blue and dark red), and the wave propagates within a few milliseconds. Therefore it is impossible to visually detect this electromechanical wave on the grayscale images.
- the wave speeds as measured using the techniques described above are represented in FIG. 4 .
- the wave velocity was found to be approximately 0.6 m/s in the posterior wall, which was corroborated by invasive electrophysiological measurements using a matrix of electrodes. Temporary regional ischemia was then induced by coronary artery ligation.
- the velocity of the electromechanical wave was found to increase up to approximately 1.7 m/s in the ischemic region. Although not entirely understood, this strong increase is believed to be due to an increase of the shear modulus in the ischemic region or a change in the conduction velocity, or both. (A second electromechanical wave has also been detected at the end-systole phase. However, due to its high propagation speed (related to the high contraction of the myocardium), the propagation was not caught with a sufficiently high frame rate. Some evidences of its propagation are detected in the human experiments, described herein.)
- the procedure 200 described hereinabove was performed on a young healthy patient.
- the transducer was placed on the patient's thorax in order to image the heart in the short axis view.
- a sequence of approximately four cardiac cycles was acquired at a very high frame rate of 170 fps using a Vingmed System Five for RF image acquisition.
- the axial displacements were processed for each frame.
- 2 electromechanical waves were clearly seen, propagating in the posterior wall of the left ventricular (not shown).
- FIGS. 12( a )-18( a ) which are consecutive displacement maps superimposed on the grayscale images ( FIGS.
- FIGS. 12( a )-18( a ) illustrate the propagation of the electromechanical wave at the end-systole phase. The speed was found to be 0.65 m/s in the posterior wall. The location of the electromechanical wavefront is indicated by arrow W in each of FIGS. 12( a )-18( a ) .
- mice were anesthetized with tribromoethanol.
- the hair was removed using potassium thioglycolate and the mouse was placed in the supine position on a heating stage (VisualSonics, Toronto ON, Canada) in order to keep the body temperature steady.
- ECG signal was obtained from the extremities.
- the ultrasound probe was placed on the chest or the abdominal wall using degassed ultrasound gel (Aquasonic 100, Parker Laboratories Inc., Fairfield N.J., USA) as a coupling medium.
- the high frequency ultrasound probe was composed of a single focused transducer working at 30 MHz, with a focal depth of 12.7 mm.
- the transducer was mechanically rotated and real-time 2D images could be acquired at a frame rate of up to 60 Hz.
- the field of view was 12 ⁇ 12 mm, the axial resolution was 50 microns, and the lateral resolution was 100 microns.
- a digitizer (2 channels, 200 MS/s, 14 bits, CS14200, Gage Applied Technologies, Lachine QC, Canada) mounted on a PC computer slot was connected to the analog RF-output of the ultrasound scanner. In addition, two TTL outputs were used to trigger the digitizer on the 2D frames. This setup allows the real-time acquisition of more than one thousand 2D RF-data, e.g., images.
- the ultrasound probe was placed on the chest in the parasternal position to obtain a longitudinal (long-axis) view of the left ventricle of the heart.
- the probe could also be positioned over the abdomen to obtain a longitudinal view of the abdominal aorta.
- a high frame rate acquisition mode (EKV) was provided on the scanner in the exemplary embodiment in order to allow detailed visualization of the heart contraction.
- the equipment can operate as quickly as 8000 frames per minute, although the user can see 1000 frames per minute due to dropped calls.
- the ultrasound acquisition of each RF-line was triggered on the mouse ECG.
- the transducer was slowly rotated and for each position of the transducer, ultrasound echo signals were recorded with a pulse repetition frequency (PRF) of 8000 pulses/s during several cardiac cycles.
- PRF pulse repetition frequency
- the ECG was simultaneously recorded and thus allowed for the synchronization of the RF-lines based on the R-wave peak, a reliable peak of the ECG during the cardiac cycle.
- the complete acquisition duration was approximately 5 min.
- RF-signals and ECG signals were digitized during the EKV acquisition and transferred to the computer in real-time.
- the data were then processed off-line, RF-lines were synchronized using the R-wave peak of the ECG signal, and a complete set of 2D ultrasound RF-data was reconstructed at 8000 fps for one complete cardiac cycle (approximately 150 ms).
- the motion of the tissue was estimated off-line using a well-known classical speckle tracking algorithm (Bonnefous, O. and P. Pesque. Time domain formulation of pulse-Doppler ultrasound and blood velocity estimation by cross correlation. Ultrason Imaging (1986) 8(2): 73-85).
- This technique was based on detecting the small local displacements of the tissue that occur between two consecutive frames. With the current method, only axial displacements (in the direction of the transducer) were computed. In this algorithm, the time-shifts in the backscatterered signals were determined between the two consecutive frames through cross-correlation of small sliding windows over the entire ultrasound image. This technique allowed the detection of very small displacements on the order of 1 ⁇ m or less (correlation windows of 150 ⁇ m, overlapping 90%). Finally, the movie of the axial displacements was processed at a frame rate up to 8000 frame/s for the entire cardiac cycle. It is understood the lateral displacement can be obtained using the same technique.
- the axial displacements were analyzed in the frequency domain as a function of the time.
- a sliding Blackman window (100 points, 25 ms) as is known in the art, was moved along the displacement variation at a fixed depth, in steps of 2 ms.
- the windowed signals were zero-padded to 8192 points and their FFT was calculated.
- the frequency content of the displacements was evaluated graphically by plotting these spectra as a function of time. Based on this frequency analysis, the transient and the slow motions of the tissues were separated using a digital filter.
- the phase velocity of the vibration was determined for an angular frequency co.
- the wave was assumed to propagate with a velocity c in a direction r that was arbitrarily determined on the image by the direction of the wall, and a set of measurement points was selected on this direction.
- the phase was measured as a function of the propagation distance r, using the Fourier Transform of the temporal displacements at the location r computed at the angular frequency co.
- the derivative of the phase of the wave with respect to distance was estimated using a linear regression fit on the set of measurements points, and the velocity of the wave at the frequency f was calculated:
- the theory of elastic wave propagation in soft biological tissue was considered in order to derive the Young's modulus of the tissue. Assuming that the medium is infinite and isotropic, the speed of shear waves propagation could be derived from general equations of the dynamic theory of elasticity. However, it is understood that the propagation of elastic waves in the myocardium can optionally take into account additional characteristics such as the active properties of the muscle fibers, the strong anisotropy of the tissue, and/or the geometry of the ventricles.
- c is the velocity of the wave
- E is the Young's modulus of the conduit wall
- h is the wall thickness
- p is the density of the fluid
- R the radius of the tube.
- FIG. 19 shows a B-mode image 1910 of a typical parasternal long-axis view obtained in a normal mouse.
- Image 1910 shows the main structures of the left ventricle: the intraventricular septum 1912 , the cavity of the left ventricle 1914 , the papillary muscle inside the cavity 1916 and the posterior wall 1918 which is visible due to strong reflections at the epicardium-lung interface.
- image 1910 is also shown in image 1910 .
- the right ventricle 1920 , aortic valve 1922 , aorta 1924 , mitral valve 1926 , and left atrium 1928 .
- the duration of the average cardiac cycle was 138 ms.
- Axial displacements were estimated for the complete set of data. In order to keep the displacements at appropriate magnitudes for the estimation (on the order of 1 ⁇ m and to reduce the amount of data, the number of frames was halved, which also reduced the frame rate to 4000 fps.
- FIGS. 20( a ) and 20( b ) show the color-coded axial displacements overlaid onto the grayscale B-mode image for two different phases of the cardiac cycle.
- the contraction of the myocardium is shown by positive displacements (red region) of the posterior wall 2018 and negative displacements of the septum 2012 (blue region) ( FIG. 20( a ) ).
- the directions of the displacements of the posterior wall 2018 and the septum 2012 are reversed during the relaxation ( FIG. 20( b ) ).
- FIG. 20( a ) The time of acquisition of FIG. 20( a ) is indicated at point t of FIG. 20( c ) .
- the time of acquisition of FIG. 20( b ) is indicated at point t of FIG. 20( d ) .
- FIG. 21( a ) A temporal analysis of the motion was performed for single RF lines of the image.
- the axial displacement along one central line of the image (indicated by the white, dotted vertical line 2050 on FIG. 20( b ) ) is shown as a function of time in FIG. 21( a ) with the corresponding ECG signal ( FIG. 21( d ) ).
- FIGS. 2( a )-( d ) are aligned on a temporal basis.
- the displacements of the septum, the papillary muscle and the posterior wall are shown in a M-mode format over two cardiac cycles.
- FIGS. 22( a )-( d ) show a sequence of axial displacements overlaid onto the grayscale B-mode images every 0.6 ms around end-systole. This sequence uncovers a strong mechanical wave W propagating in the longitudinal direction of the ventricle along the myocardium, from the base (right side of the images) to the apex (left side).
- a transverse wave propagates along the lateral direction (i.e., in-plane, perpendicular to the beam axis).
- a second wave W′ is shown in FIGS. 22( e )-( f ) .
- the mechanical wave i.e., generated by localized vibrations in the muscle ( FIG. 21( c ) ), was visible both in the posterior wall and the septum. Its amplitude was found to be eight times higher in the septum. Only the mechanical wave propagating in the septum is described herein. A set of 60 samples was selected in the septum along the propagation direction (lateral direction of the image), and the phase of the wave was computed at different frequencies. Three frequencies were selected for which the displacement amplitude was large enough to detect, e.g. with respect to noise level: (*) 82 Hz ( ⁇ ) 246 Hz ( ⁇ ) 410 Hz. The phase velocity of the wave was computed for these frequencies and a large dispersion was found.
- the distance of propagation was plotted in FIG. 23 as a function of the phase of the wave divided by the angular frequency.
- the phase velocity was found to be 1.20 m/s at 82 Hz, 3.02 m/s at 246 Hz and 4.21 m/s at 410 Hz.
- FIGS. 24( a )-( f ) show a sequence of axial displacements overlaid to the grayscale B-mode images every 2.8 ms around the beginning of systole. A strong vibration was found in the septum, but no wave propagation was visible in the image plane. Therefore, a mechanical wave can propagate in the perpendicular direction, but was not being observed with the equipment described herein.
- FIG. 24 shows a wave propagating in the posterior wall (see the white arrows W).
- the displacements were initiated at the apex (left side of the images) and then propagated towards the base (right side).
- the phase velocity was determined using the method previously described at the frequency of 80 Hz.
- the distance of propagation was plotted in FIG. 25 as a function of the phase of the wave divided by the angular frequency.
- the phase velocity of the wave was obtained using a linear regression fit and was estimated to be 0.44 m/s.
- mice were also scanned during right-atrial pacing (at 90 ms corresponding to a heart cycle at sinus rhythm of 100 ms period; FIGS. 26( a )-( e ) ) and right-ventricular pacing (also at 90 ms; FIGS. 27( a )-( e ) ).
- Pacing was achieved using catheterization through the right side of the heart, in which the catheter carried nine electrodes that could be separately activated for varying the pacing location.
- the catheter C was within the imaging field-of-view and allowed for imaging of the pacing wave during ventricular pacing ( FIGS. 26( c ) and 27( c ) ).
- the most pronounced wave propagating during atrial pacing was the contraction wave, or wave originating at the isovolumic contraction phase, that propagated along the longitudinal direction of the myocardium initiating radial thickening (or, positive (red) displacement) in its path.
- the contraction wave was very similar to the one during sinus rhythm ( FIGS. 24( a )-( f ) ), starting at the apex right at the QRS peak and then propagating along the posterior wall (generally from right to left in the figure.)
- Right-ventricular pacing FIGS.
- FIG. 27( a )-( e ) induced a reverse direction on the contraction wave that now started from the tip of the catheter (close to the base) with two waves propagating from base to apex, one along the septum and one along the posterior wall (generally from left to right in the figure) ( FIG. 27( a )-( e ) ). Since pacing occurred using the same mouse, same sonographic view and without affecting the function of the valves or the blood flow, the reverse direction of the propagation of the wave is concluded to be induced by the change in the origin of the electrical stimulus; thereby, confirming that the wave measured is electrically induced.
- FIGS. 28( a )-( f ) show a sequence of the axial displacements in color overlaid onto the grayscale B-mode image.
- a transverse wave W started propagating on the right side of the images (heart side) and then propagated towards the left side in less than 3 ms.
- This transverse wave was generated from the sudden pressure change of the blood bolus traveling through the vessel, known as the arterial pulsive wave (Nichols, W. and M. F. O'Rourke (1998). Vascular impedance. In McDonald's: blood flow in arteries: theoretical, experimental and clinical principles. E. Arnold. London).
- the phase velocity of the pulse wave was computed at the frequency of 200 Hz.
- the distance of propagation was plotted in FIG. 29 as a function of the phase of the wave divided by the angular frequency, the phase velocity was obtained using a linear regression fit and was found to be 3.08 m/s.
- Eq. 5 Moens-Korteweg equation
- systems and techniques are provided for electromechanical activation of arrhythmias, including non-transient electromechanical activation of paroxysmal and periodic arrhythmias in humans in vivo.
- Certain treatments of cardiac arrhythmias can be utilized in clinical practice but can lack a suitable noninvasive imaging modality to provide insight into the source or focus of an arrhythmia.
- Cardiac deformations can be imaged at high temporal and spatial resolution to elucidate electrical activation sequences in normal and paced human subjects non-invasively. In this manner, such imaging can be used to improve planning and monitoring of ablation-based arrhythmia treatments.
- aspects of the disclosed subject matter include techniques to quantitatively characterize focal and reentrant arrhythmias.
- spatio-temporal maps of a full-view of the atrial and ventricular mechanics can be obtained in a single heartbeat. Such maps can illustrate with suitable detail the electromechanical patterns of atrial flutter, fibrillation, and tachycardia.
- the electromechanical wave imaging techniques can be utilized to identify the location of the focal zone and the subsequent propagation of cardiac activation.
- exemplary electromechanical wave imaging techniques are described in International Application No.
- high frame rate ultrasound imaging of the heart can be used non-invasively and in real time to characterize lesser-known mechanical aspects of atrial and ventricular arrhythmias.
- Such techniques can also be used to assist treatment planning for intraoperative and longitudinal monitoring of arrhythmias.
- ultrasound imaging systems such as software-based systems can allow ultra-high frame rates, and thus ultrasound imaging can be used to allow unprecedented temporal resolution.
- ultrasound imaging systems can provide a five-fold improvement in the signal-to-noise ratio of cardiac motion and deformation mapping.
- frame rates up to 2000-5000 frames/s can be achieved by using defocussed transmissions, which can be suitable for depths utilized in transthoracic cardiac applications.
- ultrasound imaging techniques described herein can be used to map transient strains occurring in response to the electrical activation, (e.g., the electromechanical wave).
- the electrical activation e.g., the electromechanical wave
- such techniques can be used to map transmural activation sequences of normal and abnormal hearts and to locate pacing sites in patients undergoing cardiac resynchronization therapy.
- systems and techniques are provided to analyze and characterize the mechanical behavior of the atria.
- systems and techniques disclosed herein can be utilized to analyze and characterize the atria during certain types of cardiac arrhythmia, including and not limitation to, premature ventricular complex, focal tachycardia, atrial flutter, and atrial fibrillation.
- Electromechanical Wave Imaging can be suitable to characterize focal rhythms such as premature ventricular complex and focal tachycardia
- EWI can have difficulty describing and/or characterizing reentrant rhythms such as atrial flutter and fibrillation.
- systems and techniques described herein are provided to characterize electromechanical strains, including and without limitation, during reentrant rhythms based on Fourier analysis.
- Exemplary embodiments of the disclosed subject matter can include a single acquisition sequence of electromechanical activation mapping that can be used for standard EWI and/or for Fourier analysis of electromechanical strains.
- Electromechanical activation mapping can characterize electromechanical strain propagation patterns during both focal and reentrant arrhythmias. In this manner, systems and techniques described herein can determine that local deformations of the atria can be closely correlated with their electrical activation. As such, systems and techniques described herein can be used to determine characteristics of cardiac mechanics in arrhythmia, to plan ablation treatments, and to monitor the efficacy of such treatments non-invasively, longitudinally and in real-time.
- exemplary techniques for imaging electromechanical activation of arrhythmias are described.
- the systems and techniques described herein can be performed, for purpose of illustration and not limitation, on human subjects.
- the human subjects can undergo a diagnostic ultrasound scan, and as embodied herein, can occur a few minutes to a few hours prior to electroanatomic mapping and ablation.
- a normal human subject can be imaged as a control for purpose of comparison.
- strain maps can be generated, for example and without limitation using similar techniques as described herein for single-heartbeat electromechanical wave imaging (EWI).
- EWI electromechanical wave imaging
- a Verasonics system with a 2.5-MHz probe can be calibrated and customized to adhere to FDA standards, including measurements of mechanical index and of peak spatio-temporal-average intensity.
- the Verasonics system can be calibrated to have an acoustic power output that is similar to conventional clinical scanners.
- Such calibration can be performed by measuring the peak pressure and/or intensity (e.g., spatial-peak temporal average intensity, also referred to as Ispta) of the Verasonics system to ensure that its mechanical index (MI) is within FDA guidelines.
- Ispta spatial-peak temporal average intensity
- the ultrasound scan can include two sequences.
- a circular ultrasonic wave can be emitted with a virtual focus of 10.2 mm behind the probe at 2000 fps during 2 seconds.
- a standard B-mode acquisition can be performed during 1.5 seconds to depict the heart anatomy.
- Frames from the motion-estimation sequence can be reconstructed by generating a plurality of beams, for example and as embodied herein 128 beams, in post-processing using a delay-and-sum algorithm with a reconstructed sampling frequency of 20 MHz.
- the motion-estimation rate and the motion-sampling rate can be set to 1000 and 2000 fps, respectively.
- the window for motion-estimation can be 9.2 mm with an overlap of 95.8% (window shift of 0.3 mm), and the kernel strain estimation can be set to 4.9 mm.
- the techniques described herein for beamforming, motion-estimation, strain estimation, spatial moving-average of the strains (12 mm by 10 lines), and the automated contour tracking technique can be performed off-line on a graphics processing circuit (embodied herein as a Tesla graphics processing unit) and a Matlab parallel processing toolbox at a computing speed of 2.4 frames/s.
- FIG. 30A illustrates strains mapped in subjects having sinus rhythm.
- FIG. 30B illustrates strains mapped in subjects having atrial flutter.
- FIG. 30C illustrates strains mapped in subjects having atrial fibrillation.
- the strains in one location can present two representative events over time that correspond approximately to the beginning and the end of systole: end-systole, and end-diastole.
- end-systole By tracking the onset (e.g., the first zero-crossing) of these representative events for each pixel of the heart walls, isochrones maps can be generated. Isochrones correlated to electrical isochrones can be obtained, for example and without limitation, by tracking the propagation front of the end-diastole electromechanical activation. As illustrated for example in FIG.
- a similar location in the left atrium (LA) illustrates that strains can be periodic, which, in some cases, can be represented by a single frequency.
- FIG. 30C a plurality of frequencies can be observed in a patient with atrial fibrillation, and as such, analysis based on the Fourier transform might can be utilized.
- the strains can be chaotic and no period of zero strains, similar to the one in FIG. 30A , can be observed.
- the onset of contraction can be determined, for example in subjects who have focal rhythms such as sinus rhythm and focal tachycardia, as the first zero-crossing of the incremental strains occurring after the onset of the P-wave on the electrocardiogram (ECG), which can utilize the EWI techniques described herein.
- ECG electrocardiogram
- a high-resolution Fourier transform can be performed using a generalized Goertzel algorithm for interpolation in Fourier space on 1.5-second long incremental strains signals for each individual pixel in the atria.
- frequencies can be converted to cycle lengths, hereinafter referred to as mechanical cycle length (MCL).
- MCL mechanical cycle length
- peak MCL maps can be generated by selecting the MCL with the highest amplitude within the physiologically-relevant 100-330 ms range for each pixel.
- peak cycle lengths histograms can be constructed and compared to the electrical cycle length measured directly during the electroanatomic mapping and ablation.
- FIG. 31 illustrates an EWI cine-loop and isochrones during focal rhythms.
- FIG. 31A illustrates the atria of a normal subject, with propagation from the right atrium (RA) to the LA.
- the electromechanical activation regions 3110 can originate in the right atrium and propagate towards the left atrium as illustrated in the exemplary EWI cine-loop in FIG. 31A .
- FIG. 31B illustrates an EWI cine-loop depicting the atria of a subject undergoing a focal atrial tachycardia, which can have a focus located high in the left atrium (LA).
- LA left atrium
- FIGS. 31C and 31D each illustrate isochrones obtained from a patient with ventricular tachycardia.
- FIG. 31C illustrates an isochrone of the ventricular tachycardia patient during sinus rhythm.
- FIG. 31D illustrates an isochrone of the ventricular tachycardia patient during premature ventricular complex.
- EWI was performed during sinus rhythm and during pre-ventricular contraction.
- the EWI isochrones obtained during sinus rhythm depict propagation from the RA, into the LA and into the ventricles, as previously shown for purpose of illustration and comparison in normal patients.
- premature ventricular complex as illustrated for example in FIG.
- the region that was activated early in the ventricle during sinus rhythm triggered the entire electromechanical activation sequence.
- electromechanical activation can originate from the lateral wall, and can propagate toward the atria and into the atria. Early activation of the septum can indicate a potential recruitment of the Purkinje network.
- FIG. 32 illustrates the electromechanical behavior of a heart undergoing atrial flutter, including the analysis of reentrant arrhythmias using a single frequency flutter case.
- FIG. 32A illustrates an exemplary peak MCL map, and as depicted, a single MCL is representative.
- the peak cycle length map of FIG. 32A indicates, for each pixel of the atria, which cycle length is representative in the Fourier spectrum.
- FIG. 32B illustrates a histogram of the cycle length which can be used to determine, among all the pixels of the atria, which cycle length represents atrial contraction. With reference to FIG. 32B , as embodied herein, one peak cycle length of 294 ms can be identified.
- FIG. 32A illustrates an exemplary peak MCL map, and as depicted, a single MCL is representative.
- the peak cycle length map of FIG. 32A indicates, for each pixel of the atria, which cycle length is representative in the Fourier spectrum.
- FIG. 32C illustrates a phase map analyzing the phase of the MCL of FIG. 32B in Fourier space.
- a propagation pattern originates from the right atrium (RA) near the tricuspid valve towards the LA.
- the phase corresponding to the 294 ms cycle length can be retrieved from the Fourier spectrum and used to map the propagation of the mechanical oscillation at 294 ms. In this manner, the propagation direction can be determined.
- the electromechanical activation propagated from the RA to the LA.
- FIG. 32D illustrates the corresponding intracardiac electrograms obtained a few hours after the imaging procedure.
- the electrical cycle length was 283 ms.
- Atrial flutter cases can exhibit different patterns. Indeed, certain cases presented with two dominant frequencies can be separated between the left and right atria, whereas certain electrophysiological data can indicate that only one reentrant circuit was present.
- FIG. 33A illustrates two examples of such atrial flutter cases.
- FIG. 33A illustrates peak MCL maps of two exemplary atrial flutter patients. As shown for example in FIG. 33A , two representative frequencies can be identified in each patient, with the shorter cycle length located in the RA. A relationship between the MCL and the electrical cycle length can be obtained, for example and as embodied herein, by performing this analysis in multiple patients, and choosing the peak cycle length closest to the electrical cycle length. For purpose of illustration and confirmation of the disclosed subject matter, such an analysis of MCL and electrical cycle length was conducted for five exemplary patients.
- FIG. 33B is a graph illustrating the results of the five exemplary patients. With reference to FIG.
- FIG. 33B illustrates that electromechanical cycle length and/or MCL can be correlated with the electrical cycle length with a correlation of 0.96 and a slope of 1.1.
- FIGS. 33C and 33D together illustrate results from one exemplary patient undergoing atrial fibrillation.
- FIG. 33C illustrates a peak MCL map depicting multiple clustered dominant frequencies. For purpose of illustration and not limitation, the separation into these dominant frequencies is illustrated by the histogram shown in FIG. 33D . With reference to FIG. 33D , the peak cycle length map during atrial fibrillation depicts further spatial fragmentation of the peak cycle length.
- aspects of the present disclosed subject matter illustrate electromechanical activation mapping to identify the site of cardiac rhythm mechanisms during arrhythmia in humans and to characterize such cardiac rhythm mechanisms, which can lead to improved treatments and clinical management.
- Certain clinical practices utilize minimally invasive techniques to obtain precise maps of the activation of the atria and ventricles. Such techniques can be costly, time-consuming, and carry some degree of risk, and hence can be challenging to provide complete activation maps before and after treatment, as well as during catheter procedures.
- aspects of the present disclosed subject matter illustrate methods for electromechanical activation mapping during reentrant and focal arrhythmias.
- exemplary techniques are provided for imaging spatiotemporal mechanics of arrhythmias with high accuracy and spatial and temporal resolutions in a full field of view in humans. These exemplary techniques can provide for characterization of an electromechanical propagation pattern and/or representative mechanical cycle lengths, which can correspond with their electrophysiological equivalents.
- focal rhythms can behave similarly to paced rhythms.
- focal rhythms can have a single source of electromechanical activation located in the vicinity of the earliest electrical activation.
- EWI can be used to characterize the propagation of electromechanical activation, which can propagate from an atria's sinus node and from the bundle branch and which can terminate in the ventricles during ventricular pacing.
- Electromechanical activation propagation patterns similar to pacing can occur in a patient during premature ventricular complexes.
- the electromechanical activation sequence of the same patient during sinus rhythm can be similar to that of normal subjects.
- the electromechanical activation propagation pattern in a patient with atrial tachycardia, can indicate a source located near the roof of the LA, in accordance with electrical mapping.
- an exemplary application of non-invasive, ultrasound-based, electromechanical activation mapping is provided, which can be performed during or prior to invasive procedures.
- prior knowledge of an electromechanical source located in the LA can allow for clinical preparation.
- such prior knowledge can be used to determine whether transseptal access would be obtained during treatment and to perform risk-benefit analysis to determine the best course of treatment (e.g., pharmacological vs. ablation treatment).
- the electromechanical activation maps can be correlated with their electrical counterpart, at least in part of the atrial tissue, during atrial flutter.
- a single representative frequency can be identified, and the phase of that frequency can indicate a propagation direction from the cavotricuspid isthmus region to the RA and LA, which can occur during typical atrial flutters.
- other behaviors can be identified in the atria, including and without limitation, one part of the atria that contracts with the same frequency as the electrical activation and another region that does not contract with the same frequency.
- mapping the mechanics of the heart can identify regions of the heart in which the mechanical and electrical activities appear to be decoupled. Further spatial fragmentation of the periodicity of the mechanics of the atria can be observed during fibrillation. Such techniques can determine characteristics of the atrial mechanics during arrhythmia, including in the progression from flutter to fibrillation and vice versa.
- deformation of the atria caused by the onset of ventricular contraction and relaxation can affect certain aspects of the techniques described herein.
- a deformation can affect frequency analyses based on multiple activation cycles, which can be due at least in part to the relatively short acquisition time of these processes. Filtering and the development of longer acquisition sequences can inhibit or prevent such atrial deformation.
- Certain non-invasive electrical mapping techniques can be utilized to examine the epicardium, and can assume an immobilized heart function.
- a mechanical assessment of the atria can be utilized, for example, by electrophysiologists or interventional cardiologists to achieve the advantages described herein.
- echocardiograms can be performed on arrhythmia patients.
- Other non-invasive electrical mapping techniques can utilize on time-consuming and costly high resolution CT or MRI scans.
- the electromechanical activation mapping systems and techniques described herein can be obtained separately from, or in conjunction with, echocardiograms.
- mapping the electromechanical activity during arrhythmias non-invasively with real-time feedback can be used determine characteristics of atrial mechanics in the evolution and perpetuation of arrhythmias. Furthermore, and as embodied herein, such a mapping can be used to predict the origin site of arrhythmias and the mechanism and monitoring of intervention outcomes.
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Biophysics (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Veterinary Medicine (AREA)
- Radiology & Medical Imaging (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physiology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Cardiology (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
Description
- This application is a divisional of U.S. patent application Ser. No. 15/048,761, entitled “System and Method for Electromechanical Activation of Arrhythmias”, filed Feb. 19, 2016, now allowed, which is a continuation in part of U.S. patent application Ser. No. 14/449,820, entitled “Systems And Methods For Electromechanical Wave Imaging of Body Structures”, filed Aug. 1, 2014, which is a continuation of U.S. patent application Ser. No. 11/433,510, entitled “Systems And Methods For Electromechanical Wave Imaging of Body Structures”, filed May 12, 2006, which issued as U.S. Pat. No. 8,858,441 on Oct. 14, 2014, which claims priority to U.S. Provisional Application No. 60/680,081 entitled “Systems And Methods For Electromechanical Wave Imaging of Body Structures”, filed on May 12, 2005, each of which is incorporated herein by reference in its entirety and from each of which priority is claimed. This application also claims priority from U.S. Provisional Application No. 62/118,402, filed Feb. 19, 2015, which is incorporated by reference herein in its entirety.
- A portion of the disclosure of this patent document contains material which is subject to copyright protection. The copyright owner has no objection to the facsimile reproduction by anyone of any portion of the patent document, as it appears in any patent granted from the present application or in the Patent and Trademark Office file or records available to the public, but otherwise reserves all copyright rights whatsoever.
- This invention was made with government support from the National Institutes of Health under Grant Nos. R01EB006042, R21HL096094 and R01HL114358. The government has certain rights in the invention.
- This present disclosed subject matter relates to techniques for imaging and detecting the propagation of mechanical waves within a body structure of a patient.
- Certain medical conditions, such as diagnosis of myocardial ischemia, can be difficult to establish in their early stages when treatment is most effective. Patients suffering from myocardial ischemia can present to an emergency room or acute care facility with typical cardiac symptoms such as chest pain, described as tightness, pressure, or squeezing, but some patients can have other symptoms such as arm or chin pain, nausea, sweating, or abdominal pain. Certain techniques such as electrocardiogram often provide inconclusive findings regarding ischemia, and sometimes can even be unable to identify situations in which ischemia has progressed to cell damage and myocardial infarction (MI). Other techniques are available for diagnosing infarction relative to its predecessor, ischemia. For example, a blood test to measure the creatine kinase-MB (CK-MB) enzyme level is used for detection of myocardial cell damage. Other serum markers include troponin I, and to a lesser extent, myoglobin. However, the blood levels of certain such compounds can take several hours to rise, so that diagnosis of MI can be delayed. Reliance on blood tests alone can result in a significant loss of time when early aggressive therapy is warranted.
- Certain less invasive diagnostic techniques have become available through the observation of mechanical properties of tissue via imaging techniques. Such evaluation of the function of the heart, cardiovascular tissue, or other body structures can be based on the mechanical interpretation of the movement of the these structures, such as, for example, the active contractions and passive relaxation of the myocardium.
- Using certain imaging techniques, the evaluation of the heart function can be based on a single mechanical interpretation of myocardial deformation. By use of these techniques, the deformations of the myocardium can be quantified over a complete cardiac cycle in order to provide some information on the myocardial viability.
- Certain low frequency mechanical vibrations in the heart are known in human patients. Certain ultrasound techniques can be used to obtain pulsive mechanical vibrations around end-systole and end-diastole in the frequency range of 25 to 100 Hz.
- Additionally, atrial arrhythmias are a known and can cause of morbidity and mortality. Certain mechanical factors, such as chamber size and wall tension, can affect the onset and perpetuation of atrial arrhythmia. Certain echocardiographic measurements can also be used to characterize atrial arrhythmias. Yet, systems and techniques to analyze the 2-D spatio-temporal evolution of the local deformations of the atria during e.g., focal tachycardia, flutter, and fibrillation, would be beneficial.
- Accordingly, there is a need for a noninvasive imaging modality which provides insight into the source or focus of an arrhythmia.
- The present disclosure provides elasticity imaging techniques to evaluate mechanical wave propagation, and provide an estimation of electrical propagation in a noninvasive manner.
- In example embodiments, the disclosed subject matter provides systems and methods for detecting wave propagation within the tissue of a patient in a series of image frames representing movement of such tissue of the body structure. Image data is acquired comprising a series of image frames corresponding to the movement of the tissue. In an exemplary embodiment, the tissue can be the wave propagation in the myocardium. In another exemplary embodiment, the movement of body tissue can be wave propagation in the arteries or the aorta.
- A correlation calculation can be performed on the image frames to generate a matrix with the location of correlation maxima representing the relative displacement between the first and second image frames, also referred to as a displacement map. A video can be generated comprising a series of displacement maps. The parameters of movement of the cardiac structure can be detected, such as velocity, attenuation, frequency, etc. The wave can be a shear wave, representative of the electrical wave propagation within the body structure.
- According to another aspect of the present disclosure, systems and methods are provided for mapping electromechanical activity during an arrhythmia. Image information of a heart of a subject can be obtained using an imaging device. A strain map of the heart can be generated from the image information. Occurrences of a first electromechanical event of the heart and a second electromechanical event can be determined from the strain map. A spatio-temporal map of atrial and ventricular mechanics of the heart can be generated by tracking the onset of the first and second event for each pixel of a heart wall of the subject identified from the image information. A representative mechanical cycle associated with a contraction of the heart can be identified using the spatio-temporal map.
- For example, electromechanical activation mapping can characterize propagation patterns of electromechanical strains during focal and reentrant arrhythmias of the heart.
- Additionally, regions in which the mechanical and electrical activities are decoupled can be identified by mapping the electromechanical activity of the heart.
- Furthermore, the spatio-temporal map can be generated by obtaining isochrones strongly correlated to electrical isochrones by tracking a propagation front of an end-diastole electromechanical activation of the heart.
- In addition, a type of cardiac arrhythmia present in the heart can be determined from the information collected in the ultrasound scan of the heart. Upon identifying that the cardiac arrhythmia present in the heart includes focal rhythms, an onset of ventricular contraction can be determined by identifying a first zero-crossing of an incremental strains occurring after an onset of a P-wave on an electrocardiogram. Upon determining that the cardiac arrhythmia present in the heart is a type of reentrant arrhythmia, a high-resolution Fourier transform can be performed using a generalized Goertzel algorithm to interpolate strain signals in Fourier space for each individual pixel in an atria of the heart.
- For example, a peak mechanical cycle length (MCL) map can be generated by selecting a MCL having a highest amplitude within the physiologically-relevant time range for each pixel of the ultrasound scan of the heart, such that the MCL map identifies, for each pixel of the atria, which cycle length is most present in a Fourier spectrum of cycle lengths. The cycle length best representing an atrial contraction of the heart can be determined. A phase corresponding to the determined cycle length can be determined to map a propagation of a mechanical oscillation of the heart at the determined cycle length.
- Additionally, the ultrasound scan further can be performed by emitting a circular ultrasonic wave to instruct an ultrasound apparatus to perform a motion estimation sequence. A B-mode acquisition can be performed to capture heart anatomy of the heart. A plurality of beams can be generated to reconstruct frames from the motion estimation sequence using a delay-and-sum algorithm with a reconstructed sampling frequency.
- The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawings will be provided by the U.S. Patent and Trademark Office upon request and payment of the necessary fee.
-
FIG. 1 is a diagram illustrating the system in accordance with the present disclosure. -
FIG. 2 is a diagram illustrating exemplary stages in a method in accordance with the present disclosure. -
FIG. 3 is a diagram illustrating a technique for measuring movement of structures within an image in accordance with the present disclosure. -
FIG. 4 is a chart representing the velocity of structures within an image in accordance with the present disclosure. -
FIG. 5 illustrates a technique of detecting wave propagation in accordance with a further embodiment of the present disclosure. -
FIG. 6 illustrates a technique of detecting wave propagation in accordance with yet another embodiment of the present disclosure. -
FIGS. 7(a)-11(b) are images illustrating the propagation of a wave within a body structure in accordance with an exemplary embodiment of the present disclosure. -
FIGS. 12(a)-18(b) are images illustrating the propagation of a wave within a body structure in accordance with another exemplary embodiment of the present disclosure. -
FIG. 19 is an ultrasound image of a mouse left ventricle in a parasternal long-axis view. -
FIG. 20(a) is an axial displacement map overlaid to the grayscale B-mode image of the left ventricle during systole in accordance with the present disclosure. -
FIG. 20(b) is an axial displacement map overlaid to the grayscale B-mode image of the left ventricle during diastole (relaxation phase) in accordance with the present disclosure. -
FIG. 20(c) is an ECG indicating the time of the acquisition during the cardiac cycle ofFIG. 20(a) in accordance with the present disclosure. -
FIG. 20(d) is an ECG indicating the time of the acquisition during the cardiac cycle ofFIG. 20(b) in accordance with the present disclosure. -
FIG. 21(a) is a time plot illustrating the temporal variation of the axial displacements estimated on one central RF-line as line plotted onFIG. 20(b) in accordance with the present disclosure. -
FIG. 21(b) illustrates the frequency content of the displacement variation in the septum at the depth of 12.5 mm plotted as a function of time in accordance with the present disclosure. -
FIG. 21(c) is a time plot illustrating the temporal variation of the axial displacements after bandpass filtering of the plot illustrated inFIG. 21(a) showing the transient and high frequency components in accordance with the present disclosure. -
FIG. 21(d) illustrates the ECG signal acquired simultaneously with the data illustrated inFIGS. 21(a)-(c) in accordance with the present disclosure. -
FIGS. 22(a)-(d) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the left ventricle around end-systole taken every 0.6 ms showing the propagation of a first mechanical wave front in the septum in accordance with the present disclosure. The arrows indicate the progression of the wave front in the septum. -
FIGS. 22(e)-(f) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the left ventricle around end-systole taken every 0.6 ms showing the propagation of a second mechanical wave front in the septum in accordance with the present disclosure. The arrows indicate the progression of the wave front in the septum. -
FIGS. 22(g)-(l) illustrate the ECG signal plotted below each respective image ofFIGS. 22(a)-(f) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure. -
FIG. 23 is a plot illustrating the distance of propagation as a function of the phase of the end-systolic wave at three frequencies in accordance with the present disclosure. -
FIGS. 24(a)-(f) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the left ventricle around the beginning of systole taken every 2.8 ms, showing the propagation of a strong mechanical wave in the posterior wall in accordance with the present disclosure. The arrows indicate the progression of the wave front in the posterior wall. -
FIGS. 24(g)-(l) illustrate the ECG signal plotted below each respective image ofFIGS. 24(a)-(f) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure. -
FIG. 25 is a plot illustrating the distance of propagation as a function of the phase of the wave at the frequency of 80 Hz during the beginning of systole transient motion in accordance with the present disclosure. -
FIGS. 26(a)-(e) illustrate a sequence of axial displacement maps overlaid to the grayscale image (0.12 ms between successive frames) indicating an electromechanical wave propagating in the posterior wall of the mouse from the apex towards the base during pacing in the right atrium close to the sinoatrial node in accordance with the present disclosure. -
FIGS. 26(f)-(j) illustrate the ECG signal plotted below each respective image ofFIGS. 26(a)-(e) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure. -
FIGS. 27(a)-(e) illustrate a sequence of axial displacement maps overlaid to the grayscale image (0.07 ms between successive frames) indicating an electromechanical wave propagating in the posterior wall of the mouse from the base towards the apex during pacing in the right ventricle close to the base in accordance with the present disclosure. -
FIGS. 27(f)-(j) illustrate the ECG signal plotted below each respective image ofFIGS. 27(a)-(e) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure. -
FIGS. 28(a)-(f) illustrate a sequence of axial displacement maps overlaid to the grayscale B-mode image of the aorta taken every 0.7 ms. Sequence of images showing the propagation of the pulse wave in the aorta. The arrows indicate the progression of the wave front in the aorta -
FIGS. 28(g)-(l) illustrate the ECG signal plotted below each respective image ofFIGS. 26(a)-(f) indicating the time t of the acquisition during the cardiac cycle in accordance with the present disclosure. -
FIG. 29 illustrates the distance of propagation as a function of the phase of pulse wave at the frequency of 80 Hz. The slope of the curve gives the pulse wave velocity in accordance with the present disclosure. -
FIGS. 30(a)-(c) illustrates examples of high temporal resolution strains during different types of arrhythmia in accordance with the present disclosure. -
FIGS. 31(a)-(d) illustrates examples of propagating electromechanical activation in atria of normal subjects and a patient undergoing focal atrial tachycardia in accordance with the present disclosure. -
FIGS. 32(a)-(d) illustrates analysis of reentrant arrhythmias using a single-frequency atrial flutter case using a peak cycle length map, cycle length histogram, a phase map, and electrogram in accordance with the present disclosure. -
FIG. 33(a) illustrates peak cycle length maps in two atrial flutter patients in accordance with the present disclosure.FIG. 33(b) illustrates a correlation between mechanical cycle length (MCL) and electrical cycle length in five atrial flutter patients in accordance with the present disclosure.FIG. 33(c) illustrates a peak cycle length map during atrial fibrillation in accordance with the present disclosure.FIG. 33(d) illustrates a histogram depicting spatial fragmentation of peak cycle length during atrial fibrillation in accordance with the present disclosure. - Throughout the figures, the same reference numerals and characters, unless otherwise stated, are used to denote like features, elements, components or portions of the illustrated embodiments. Moreover, while the present disclosure will now be described in detail with reference to the figures, it is done so in connection with the illustrative embodiments.
- The system and methods described herein can be useful for analyzing data obtained by an image generating device, such as an ultrasound transducer. The systems and methods can also be useful for measuring mechanical properties and estimating the electrical characteristics of a body tissue structure or organ, such as, for example, the myocardium or the aorta.
- For example, the disclosed subject matter can be used in connection with imaging and characterizing the propagation of electromechanical waves in the heart. During the cardiac cycle, electrical waves propagate in the myocardium in order to induce its contraction. Contraction of the myocardial fibers results in a strong mechanical wave, which, since it results from the coupling of the electrical excitation and the mechanical properties of the myocardium, is referred to herein as an “electromechanical wave.” The speed of this wave is a function of the electrical and mechanical properties of the myocardium, and, according to the present disclosure, can be used to detect changes in these properties to diagnose heart diseases.
- An exemplary embodiment of the system is illustrated in
FIG. 1 and designatedsystem 100.System 100 can include an image detection device, such asultrasound probe 102, which is used to create images of the heart H or other organ or structure of the patient P. The image detection device does not induce discernible vibration in the body structure, and merely detects pre-existing motion. The signals detected by theprobe 102 can be transferred to anultrasound scanner 104. The exemplary embodiments described herein are designed to work with conventional ultrasound scanners. For example, commercial portable scanners, such as Terason 2000, high frequency scanners, such as Visualsonics Vevo 770, and routinely used clinical scanners, such as GE System Five or GE Vivid Five or Seven, are useful for image acquisition in accordance with the exemplary embodiments. The raw data produced by thescanner 104 can be transferred to acomputer 106 having aCPU 108 for processing the data. In the exemplary embodiment, the computer and CPU would be Dell PC with a 2 GHz processor. It is understood that the computer and CPU can also be integrated with theultrasound scanner 104. Also useful in the system would be storage such as disk drives, for storing data oninput files 110 and for writing output onto output files 112. As will be described herein, input files 110 can include information such as thresholds. Output files 112 can include the displacement maps, videos of myocardium displacements, or computed data, such as electromechanical wave properties. It is understood that a preprogrammed chip can be used to execute the algorithms described herein. Typically, an output device, such asmonitor 114, and an input device, such askeyboard 116, are also components of the system. - In accordance with an exemplary embodiment, the methods described herein are particularly useful for imaging the propagation of electromechanical waves in the heart. A method for detecting the properties of the electromechanical wave are described herein and represented in
FIG. 2 . In an early stage in the procedure, raw imaging data of the body structure is acquired by image acquisition equipment such as theultrasound probe 102 andscanner 104. In the exemplary embodiment, a set of N frames of raw ultrasound data of the heart is acquired during a cardiac cycle at high frame rate, e.g., higher than 100 fps, although frame rates of about 56 fps and 170 fps, etc., yield useful results (202). The selected frame rate should be commensurate with the speed of the propagation of the movement, such as the wave, being studied. The electrocardiogram (EKG) can also be recorded. The raw data can be digitized and stored in real-time in the scanner memory. - In a subsequent stage, the data can be transferred to a computer for processing (204). In an exemplary embodiment, the transfer can occur using a protocol such as Ethernet TCP IP. This is optional, as the computer can be integrated with the
scanner 104. - At 206, the raw data received from the image acquisition equipment is processed. In the exemplary embodiment, the data processing computes an estimation of the displacement of particular objects in the images, such as the myocardium, between consecutive frames. Typically this processing occurs off-line; however, it is understood that this procedure can occur sequentially subsequent to receiving two consecutive frames. According to the exemplary embodiment, axial displacements (in the direction of the transducer) are computed. Lateral, or elevational, displacements (perpendicular to the transducer) can also be computed using a similar technique, for example, as disclosed in Konofagou E. E. and Ophir, J., (1998), A New Elastographic Method for Estimation and Imaging of Lateral Strains, Corrected Axial Strains and Poison's Ratios in Tissues,” Ultrasound in Medicine and Biology 24(8), 1183-1199 (1998); Konofagou et al. (1998), Three-dimensional Motion estimation in Elastography, IEEE Proceedings of the Symposium of Ultrasonics, Ferroelectrics and Frequency Control in Sendai, Japan, 1745-1748. Korteweg, D. Uber die Fortpflanzungsgeschwindigkeit des Schalles in elastichen Rohren. Ann. Phys. Chem. (1879) 5: 525-37, the contents of which are incorporated herein.
- N−1
displacement 2D maps (also referred to as correlation matrices) are computed through the correlation of two consecutive frames i and i+1 (1<i<N−1). Each frame is represented by a matrix of pixel values. The displacement maps provide an indication of the local axial movements between frames. Estimation of the axial displacements from the two consecutive frames is performed using a speckle tracking algorithm. In this algorithm, the time-shifts in the backscattered signals are determined between two consecutive frames through cross-correlation of small sliding windows over the entire ultrasound image. For each window, the signal of the frame i and the frame i+1 are cross-correlated. The maximum of the correlation coefficient gives an estimation of the time-shift between the two signals. This time-shift can be converted to a spatial displacement by assuming a constant speed of sound for the tissue. This technique can detect displacements on the order of 10 μm. Using small correlation windows of 7.5 mm, the resolution of the displacement maps is in the millimeter range. The cross-correlation algorithm suitable for estimating displacement between consecutive image frames is described in U.S. Provisional Patent Application No. 60/619,247, filed Oct. 15, 2004, which is incorporated by reference herein. In the exemplary embodiment, a Matlab program Multiframe is used to compute the displacement maps for the complete sequence of frames obtained at 202, above. Multiframe calls the Matlab routine FunCalculDispl to compute the displacements for the sequence of frames. FunCalculDispl in turn calls the routine Correlation2D.cpp which is a C program that computes the displacement map between consecutive frames. As discussed above, Correlation2D.cpp uses small sliding windows to find the displacement which maximizes the correlation coefficient for each part of the image. In accordance with other embodiments of the present disclosure, auto-correlation calculations or coherence calculations, as are known in the art, can be performed. - Two optional threshold procedures can be executed in the
procedure 200. At 208, a threshold can be applied on the energy of the signal, in order to remove the noise that is below a predetermined signal-to-noise ratio. Low energy ultrasound signals (e.g., noise in the cavity of the heart) can be removed from the displacement map according to this method. At 210, a threshold can be applied on the correlation coefficient to remove erroneous estimates in the displacements. In the exemplary embodiment, the noise threshold and correlation-coefficient threshold can be implemented within the routine Correlation2D.cpp. The levels of the thresholds are determined experimentally and can be stored in an input data file 110 for processing on theCPU 108.Procedures - A video of the sequence of N−1 displacement maps can be assembled to create a video of the displacements of the body structure or tissue (212). In the exemplary embodiment, a video of the myocardium displacements is created by this technique.
- The video of the displacement map of the myocardium will depict the propagation of the electromechanical wave. Next, an observation and tracking of the wave propagation (214). Although such tracking can be done manually, it can be difficult to discern the wave by visual observation and thus make accurate measurements. Accordingly, wave tracking can be performed by an algorithm, such as TrackPositionWave, a Matlab program which locates the position of the wave front by performing a zero-crossing calculation on consecutive displacement maps.
- The parameters of the electromechanical wave, e.g., velocity, amplitude, attenuation, frequency, etc., can be analyzed at 216. For example, the velocity of the electromechanical wave can be computed as a function of its position in the myocardium. As illustrated in
FIG. 3 , the wall of the myocardium is approximated as circular with a radius R, and the origin of the spherical coordinate system was chosen at the center of the cavity. The wavefront of the electromechanical wave was then tracked by its angular coordinate θ. The Matlab function Overlay can be used to compute the transformation of the raw image into polar coordinates. This routine can also display the displacement map superimposed on the ultrasound grayscale data. As an example, the speed of the electromechanical wave is shown on theFIG. 4 as a function of the angular position. - The ultrasound imaging method described herein has the advantage of being completely non-invasive. In an exemplary embodiment, the system described herein can be implemented in real-time on commercial scanners. It has been shown that the electrical conductivity is transversely isotropic with respect to fiber direction, with a longitudinal velocity of about 0.6 m/s and a transverse velocity of about 0.2 m/s (Roth, B. J. (2000), Influence of a perfusing bath on the foot of the cardiac action potential. Circulation Research 86, E19-E22; Spach, M. S., Heidlage, J. F., Dolber, P. C., and Barr, R. C. (1998), Extracellular discontinuities in cardiac muscle—Evidence for capillary effects on the action potential foot. Circulation Research 83, 1144-1164). The electromechanical wave velocity noted herein was very close to the longitudinal velocity of the mechanical wave. The transverse velocity can be measured by using ultrasound imaging and displacement estimation using a 3D imaging probe or a rotational 2D imaging probe.
FIG. 5 illustrates a transducer setup for 2D imaging of the longitudinal waves, andFIG. 6 illustrates a transducer setup for 2D imaging of transverse waves. - The mechanical component of the electromechanical wave is related to the viscoelastic properties of the soft tissue. The elastic properties of the myocardium have been widely investigated. The stiffness of the myocardium has been shown to increase during ischemia and recovers after reperfusion. Thus, early detection of cardiovascular diseases such as ischemia and infarction can be strongly improved through non-invasive characterization of the local myocardial elasticity.
- Low frequency shear (mechanical) waves propagate in soft tissue at low velocity (0.5 to 50 m/s). For an isotropic and infinite medium, it has been shown that the velocity of the shear wave is related to the shear modulus p and the density p by:
-
- (Bercoff, J., Tanter, M., and Fink, M. (2004), Supersonic shear imaging: A new technique for soft tissue elasticity mapping. IEEE Transactions on Ultrasonics Ferroelectrics and Frequency Control 51, 396-409; Sarvazyan, A. P., O. V. Rudenko, S. D. Swanson, J. B. Fowlkes and S. Y. Emelianov, Shear wave elasticity imaging: A new ultrasonic technology of medical diagnostics. Ultrasound Med Biol (1998) 24(9): 1419-1435.) According to another exemplary embodiment of the disclosed subject matter, a system can be implemented to provide early detection of ischemia through the measure of the velocity of the mechanical wave.
- However, the myocardium has also anisotropic mechanical properties and can be considered as a transverse isotropic medium. As a consequence, two shear waves of different velocities can propagate in the myocardium. Fast mechanical (shear) waves propagate in the direction of the fibers, and slow mechanical (shear) waves propagate in the direction perpendicular to the fibers. The measure of the two wave velocities can be achieved by using 3D Ultrasound imaging systems or multiple acquisitions of 2D images with a rotation of the transducer (see
FIGS. 5-6 ). The wave velocities are related to two elastic constants, μ// the shear modulus in direction of the fibers and μ⊥ the cross-fiber shear modulus. -
- The systems and methods described herein can potentially have different applications in the field of early detection of cardiovascular diseases and cardiac imaging.
- For example, the measure of the electrical excitation propagation is of high interest in cardiology for early detection of heart diseases but also for pacing the heart when heartbeat is too slow or irregular. The purpose of an artificial pacemaker is to stimulate the heart when either the heart's natural pacemaker is not fast enough or if there are blocks in the electrical conduction system preventing the propagation of electrical impulses. Thus, in order to implant the artificial pacemaker at the correct location, the electrical propagation must be determined accurately. In vivo imaging of the electrical propagation in the heart can require implanting an electrode matrix (up to 500 electrodes) to measure extracellular potentials at the surface of the heart. This invasive and potentially precarious surgical procedure cannot be performed on human patients for diagnostic purposes. The present disclosure provides a means for determining electrical propagation in the myocardium of a subject, in the context of achieving an effective position of a pacemaker in the subject. Other methods known in the art involve optically-based techniques which also require invasive procedures, such as open-chest surgery.
- The present disclosure can be further used to create images and thereby detect myocardial ischemia in a subject either having symptoms (e.g., chest, arm, or chin pain, nausea, shortness of breath, and/or sweating) or a subject subjectively lacking such symptoms (e.g., “silent ischemia”), whereby a finding of increased electromechanical wave velocity (relative to control values) in a region of the myocardium of a subject is consistent with and supportive of a diagnosis of myocardial ischemia in that region. The present disclosure can also be used to diagnose, or assist in surgical intervention in, (i) conduction disturbances, such as re-entry phenomena, or associated with pharmaceutical agents, such as antidepressants or hyperkalemia, (ii) arrhythmias and dysrhythmias (e.g., surgical treatment of ventricular dysrhythmias, diagnosis of low-amplitude atrial fibrillation); and (iii) tissue abnormalities associated with cardiomyopathies or trauma, etc.
- The procedure described hereinabove was performed in an anesthetized open-chested dog. The transducer was placed on the anterior wall of the left ventricle of the heart, to obtain a short axis view. Approximately every two minutes, a sequence of three cardiac cycles was acquired during the experiment, with a frame rate of 56 fps. The 2D displacement maps were estimated using the cross-correlation method (window size: 5 mm, 90% overlap). The axial displacements were processed for the different sequences. On the displacement video, two electromechanical waves were clearly detected, propagating in the posterior wall of the left ventricular, from the septum (left side of the images) to the lateral wall (right side). The propagation of the mechanical wave corresponds to the electrical activity shown on an associated EKG.
- The first electromechanical wave is found at the end-diastolic phase of the cardiac cycle (which corresponds to the beginning of the contraction).
FIGS. 5(a) through 9(a) show five consecutive frames of the propagation of the wave. The location of the electromechanical wavefront is indicated by arrow W in each ofFIGS. 7(a)-11(a) . The displacements maps are overlaid to the grayscale ultrasound images (see,FIGS. 7(b)-11(b) ). Blue displacements are in the direction of the transducer (top of the image), and red displacements are in the opposite direction. As shown in these images, the contraction of the myocardium starts on the left side (septum) and propagates to the right side of the image. In the figures, the blue region appears on the left side of the images (behind the wavefront), and the red region appears on the right side of the image (in front of the wavefront). The maximum displacements shown are 75 μm (dark blue and dark red), and the wave propagates within a few milliseconds. Therefore it is impossible to visually detect this electromechanical wave on the grayscale images. The wave speeds as measured using the techniques described above are represented inFIG. 4 . The wave velocity was found to be approximately 0.6 m/s in the posterior wall, which was corroborated by invasive electrophysiological measurements using a matrix of electrodes. Temporary regional ischemia was then induced by coronary artery ligation. The velocity of the electromechanical wave was found to increase up to approximately 1.7 m/s in the ischemic region. Although not entirely understood, this strong increase is believed to be due to an increase of the shear modulus in the ischemic region or a change in the conduction velocity, or both. (A second electromechanical wave has also been detected at the end-systole phase. However, due to its high propagation speed (related to the high contraction of the myocardium), the propagation was not caught with a sufficiently high frame rate. Some evidences of its propagation are detected in the human experiments, described herein.) - The
procedure 200 described hereinabove was performed on a young healthy patient. The transducer was placed on the patient's thorax in order to image the heart in the short axis view. A sequence of approximately four cardiac cycles was acquired at a very high frame rate of 170 fps using a Vingmed System Five for RF image acquisition. In order to reach such a high frame rate, only a small part of the heart (the left ventricle) was imaged (80×40 mm). The axial displacements were processed for each frame. On the displacement video, 2 electromechanical waves were clearly seen, propagating in the posterior wall of the left ventricular (not shown).FIGS. 12(a)-18(a) , which are consecutive displacement maps superimposed on the grayscale images (FIGS. 12(b)-18(b) ), illustrate the propagation of the electromechanical wave at the end-systole phase. The speed was found to be 0.65 m/s in the posterior wall. The location of the electromechanical wavefront is indicated by arrow W in each ofFIGS. 12(a)-18(a) . - The procedure described hereinabove was performed on anesthetized mice. The mice were anesthetized with tribromoethanol. The hair was removed using potassium thioglycolate and the mouse was placed in the supine position on a heating stage (VisualSonics, Toronto ON, Canada) in order to keep the body temperature steady. ECG signal was obtained from the extremities. The ultrasound probe was placed on the chest or the abdominal wall using degassed ultrasound gel (
Aquasonic 100, Parker Laboratories Inc., Fairfield N.J., USA) as a coupling medium. - An ultrasound scanner specifically developed for imaging small animals (Vevo 770, Visualsonics, Toronto ON, Canada) was used in this exemplary embodiment. The high frequency ultrasound probe was composed of a single focused transducer working at 30 MHz, with a focal depth of 12.7 mm. The transducer was mechanically rotated and real-
time 2D images could be acquired at a frame rate of up to 60 Hz. The field of view was 12×12 mm, the axial resolution was 50 microns, and the lateral resolution was 100 microns. - A digitizer (2 channels, 200 MS/s, 14 bits, CS14200, Gage Applied Technologies, Lachine QC, Canada) mounted on a PC computer slot was connected to the analog RF-output of the ultrasound scanner. In addition, two TTL outputs were used to trigger the digitizer on the 2D frames. This setup allows the real-time acquisition of more than one thousand 2D RF-data, e.g., images.
- In the exemplary embodiment, the ultrasound probe was placed on the chest in the parasternal position to obtain a longitudinal (long-axis) view of the left ventricle of the heart. The probe could also be positioned over the abdomen to obtain a longitudinal view of the abdominal aorta.
- In addition to the real-time scanning mode, a high frame rate acquisition mode (EKV) was provided on the scanner in the exemplary embodiment in order to allow detailed visualization of the heart contraction. The equipment can operate as quickly as 8000 frames per minute, although the user can see 1000 frames per minute due to dropped calls. Using this technique, the ultrasound acquisition of each RF-line was triggered on the mouse ECG. The transducer was slowly rotated and for each position of the transducer, ultrasound echo signals were recorded with a pulse repetition frequency (PRF) of 8000 pulses/s during several cardiac cycles. The ECG was simultaneously recorded and thus allowed for the synchronization of the RF-lines based on the R-wave peak, a reliable peak of the ECG during the cardiac cycle. The complete acquisition duration was approximately 5 min.
- To compute the tissue motion, RF-signals and ECG signals were digitized during the EKV acquisition and transferred to the computer in real-time. The data were then processed off-line, RF-lines were synchronized using the R-wave peak of the ECG signal, and a complete set of 2D ultrasound RF-data was reconstructed at 8000 fps for one complete cardiac cycle (approximately 150 ms).
- The motion of the tissue was estimated off-line using a well-known classical speckle tracking algorithm (Bonnefous, O. and P. Pesque. Time domain formulation of pulse-Doppler ultrasound and blood velocity estimation by cross correlation. Ultrason Imaging (1986) 8(2): 73-85). This technique was based on detecting the small local displacements of the tissue that occur between two consecutive frames. With the current method, only axial displacements (in the direction of the transducer) were computed. In this algorithm, the time-shifts in the backscatterered signals were determined between the two consecutive frames through cross-correlation of small sliding windows over the entire ultrasound image. This technique allowed the detection of very small displacements on the order of 1 μm or less (correlation windows of 150 μm, overlapping 90%). Finally, the movie of the axial displacements was processed at a frame rate up to 8000 frame/s for the entire cardiac cycle. It is understood the lateral displacement can be obtained using the same technique.
- The axial displacements were analyzed in the frequency domain as a function of the time. A sliding Blackman window (100 points, 25 ms) as is known in the art, was moved along the displacement variation at a fixed depth, in steps of 2 ms. The windowed signals were zero-padded to 8192 points and their FFT was calculated. The frequency content of the displacements was evaluated graphically by plotting these spectra as a function of time. Based on this frequency analysis, the transient and the slow motions of the tissues were separated using a digital filter. The displacement estimates were temporally filtered using an FIR band-pass filter with cut-frequencies of f1=50 Hz and f2=500 Hz, which allows the removal of the low frequency components but also the high frequency noise.
- To analyze the propagation of the mechanical waves, the phase velocity of the vibration was determined for an angular frequency co. The wave was assumed to propagate with a velocity c in a direction r that was arbitrarily determined on the image by the direction of the wall, and a set of measurement points was selected on this direction. The wave number is k=ω/c, and the phase of the wave is φ(r)=kr along the direction of propagation. The phase was measured as a function of the propagation distance r, using the Fourier Transform of the temporal displacements at the location r computed at the angular frequency co. Finally, the derivative of the phase of the wave with respect to distance was estimated using a linear regression fit on the set of measurements points, and the velocity of the wave at the frequency f was calculated:
-
c(f)=2πf/(∂φ/∂r) (4) - The theory of elastic wave propagation in soft biological tissue was considered in order to derive the Young's modulus of the tissue. Assuming that the medium is infinite and isotropic, the speed of shear waves propagation could be derived from general equations of the dynamic theory of elasticity. However, it is understood that the propagation of elastic waves in the myocardium can optionally take into account additional characteristics such as the active properties of the muscle fibers, the strong anisotropy of the tissue, and/or the geometry of the ventricles.
- For the transverse wave on the artery wall, a simple model of the propagation of a pressure wave in a viscoelastic infinite thin conduit filled with an incompressible fluid is well described by the Moens-Korteweg equation:
-
- where c is the velocity of the wave, E is the Young's modulus of the conduit wall, h is the wall thickness, p is the density of the fluid and R the radius of the tube. According to this equation, the elasticity of the vessel wall can be derived from the measurement of the pulse wave velocity in the artery.
-
FIG. 19 shows a B-mode image 1910 of a typical parasternal long-axis view obtained in a normal mouse. Image 1910 shows the main structures of the left ventricle: theintraventricular septum 1912, the cavity of theleft ventricle 1914, the papillary muscle inside thecavity 1916 and theposterior wall 1918 which is visible due to strong reflections at the epicardium-lung interface. Also shown in image 1910 is theright ventricle 1920,aortic valve 1922,aorta 1924,mitral valve 1926, and leftatrium 1928. In this embodiment, the duration of the average cardiac cycle was 138 ms. Axial displacements were estimated for the complete set of data. In order to keep the displacements at appropriate magnitudes for the estimation (on the order of 1 μm and to reduce the amount of data, the number of frames was halved, which also reduced the frame rate to 4000 fps. -
FIGS. 20(a) and 20(b) show the color-coded axial displacements overlaid onto the grayscale B-mode image for two different phases of the cardiac cycle. During the systolic phase, the contraction of the myocardium is shown by positive displacements (red region) of theposterior wall 2018 and negative displacements of the septum 2012 (blue region) (FIG. 20(a) ). In the diastolic phase, the directions of the displacements of theposterior wall 2018 and the septum 2012 (and the colors associated with the direction of movement) are reversed during the relaxation (FIG. 20(b) ). It should be noted that even if a large part of the myocardium of the posterior wall is not visible, the motion of the epicardium undergoes similar motion. The time of acquisition ofFIG. 20(a) is indicated at point t ofFIG. 20(c) . The time of acquisition ofFIG. 20(b) is indicated at point t ofFIG. 20(d) . - A temporal analysis of the motion was performed for single RF lines of the image. The axial displacement along one central line of the image (indicated by the white, dotted
vertical line 2050 onFIG. 20(b) ) is shown as a function of time inFIG. 21(a) with the corresponding ECG signal (FIG. 21(d) ). (FIGS. 2(a)-(d) are aligned on a temporal basis.) On thisline 2050, the displacements of the septum, the papillary muscle and the posterior wall are shown in a M-mode format over two cardiac cycles. It shows the successive main phases of the cardiac cycle: the contraction of the myocardium (systole) indicated byarrow 2120 initiated at the R-wave peak of the ECG, followed by the relaxation phase (diastole) indicated byarrow 2130. The duration of the active contraction was approximately 50 ms, and that of the relaxation 35 ms. In addition to this slow and large motion, some rapid transient variations of a few ms were observed at the beginning and at the end of the systolic phase, in the septum and the posterior wall. - In order to separate the electromechanical wave from other mechanical waves generated by vibrations resulting from valve functions or blood flow, high-pass filtering was performed. The frequency content of the tissue displacements resulting from vibrations in the septum (at depth of 12.5 mm) was analyzed as a function of time and is shown in
FIG. 21(b) . During the contraction and the relaxation of the heart, the motion of tissue was found to be in the low frequency range of up to 60 Hz. However, during the transient motion at the end of systole, much larger frequency components were found that ranged between 50 Hz and 500 Hz. The same effect was found for the transient motion at the beginning of systole, but the frequency range was limited between 50 Hz and 250 Hz. Thus, it was possible to almost completely separate the transient part of the displacement by filtering out the low frequency component of the motion. After filtering the displacements using a FIR band-pass filter with cut-off frequencies f1=50 Hz and f2=500 Hz, the two vibrations were clearly visible and are shown onFIG. 21(c) asregions - In order to analyze spatially the vibration around end-systole, we considered the data between 52 ms and 70 ms after the peak of the R-wave.
FIGS. 22(a)-(d) show a sequence of axial displacements overlaid onto the grayscale B-mode images every 0.6 ms around end-systole. This sequence uncovers a strong mechanical wave W propagating in the longitudinal direction of the ventricle along the myocardium, from the base (right side of the images) to the apex (left side). In other words, as the tissue locally vibrates along the axial direction of the beam (i.e., along the beam axis), a transverse wave propagates along the lateral direction (i.e., in-plane, perpendicular to the beam axis). A second wave W′ is shown inFIGS. 22(e)-(f) . - The mechanical wave, i.e., generated by localized vibrations in the muscle (
FIG. 21(c) ), was visible both in the posterior wall and the septum. Its amplitude was found to be eight times higher in the septum. Only the mechanical wave propagating in the septum is described herein. A set of 60 samples was selected in the septum along the propagation direction (lateral direction of the image), and the phase of the wave was computed at different frequencies. Three frequencies were selected for which the displacement amplitude was large enough to detect, e.g. with respect to noise level: (*) 82 Hz (●) 246 Hz (⋄) 410 Hz. The phase velocity of the wave was computed for these frequencies and a large dispersion was found. The distance of propagation was plotted inFIG. 23 as a function of the phase of the wave divided by the angular frequency. The phase velocity was found to be 1.20 m/s at 82 Hz, 3.02 m/s at 246 Hz and 4.21 m/s at 410 Hz. - The same analysis was performed at the beginning of systole. The filtered data were processed between 0 ms and 20 ms from the peak of the R-wave.
FIGS. 24(a)-(f) show a sequence of axial displacements overlaid to the grayscale B-mode images every 2.8 ms around the beginning of systole. A strong vibration was found in the septum, but no wave propagation was visible in the image plane. Therefore, a mechanical wave can propagate in the perpendicular direction, but was not being observed with the equipment described herein. - However, the
FIG. 24 shows a wave propagating in the posterior wall (see the white arrows W). The displacements were initiated at the apex (left side of the images) and then propagated towards the base (right side). The phase velocity was determined using the method previously described at the frequency of 80 Hz. The distance of propagation was plotted in FIG. 25 as a function of the phase of the wave divided by the angular frequency. The phase velocity of the wave was obtained using a linear regression fit and was estimated to be 0.44 m/s. - In order to determine that the origin and direction of the wave W were electrically induced and driven, mice were also scanned during right-atrial pacing (at 90 ms corresponding to a heart cycle at sinus rhythm of 100 ms period;
FIGS. 26(a)-(e) ) and right-ventricular pacing (also at 90 ms;FIGS. 27(a)-(e) ). Pacing was achieved using catheterization through the right side of the heart, in which the catheter carried nine electrodes that could be separately activated for varying the pacing location. In some of the scans, the catheter C was within the imaging field-of-view and allowed for imaging of the pacing wave during ventricular pacing (FIGS. 26(c) and 27(c) ). - The most pronounced wave propagating during atrial pacing was the contraction wave, or wave originating at the isovolumic contraction phase, that propagated along the longitudinal direction of the myocardium initiating radial thickening (or, positive (red) displacement) in its path. At atrial pacing (
FIG. 26(a)-(e) ), the contraction wave was very similar to the one during sinus rhythm (FIGS. 24(a)-(f) ), starting at the apex right at the QRS peak and then propagating along the posterior wall (generally from right to left in the figure.) Right-ventricular pacing (FIGS. 27(a)-(e) ) induced a reverse direction on the contraction wave that now started from the tip of the catheter (close to the base) with two waves propagating from base to apex, one along the septum and one along the posterior wall (generally from left to right in the figure) (FIG. 27(a)-(e) ). Since pacing occurred using the same mouse, same sonographic view and without affecting the function of the valves or the blood flow, the reverse direction of the propagation of the wave is concluded to be induced by the change in the origin of the electrical stimulus; thereby, confirming that the wave measured is electrically induced. - A longitudinal view of the abdominal aorta of a mouse was imaged using the high frame rate technique. Axial displacements were calculated, and the movie of the motion was processed at 8000 fps for a complete cardiac cycle. During the cardiac cycle, the displacements of the artery wall were found to be very small except after the beginning of systole. Strong displacements of the wall started 10.3 ms after the R-wave peak of the ECG.
FIGS. 28(a)-(f) show a sequence of the axial displacements in color overlaid onto the grayscale B-mode image. A transverse wave W started propagating on the right side of the images (heart side) and then propagated towards the left side in less than 3 ms. This transverse wave was generated from the sudden pressure change of the blood bolus traveling through the vessel, known as the arterial pulsive wave (Nichols, W. and M. F. O'Rourke (1998). Vascular impedance. In McDonald's: blood flow in arteries: theoretical, experimental and clinical principles. E. Arnold. London). - The phase velocity of the pulse wave was computed at the frequency of 200 Hz. The distance of propagation was plotted in
FIG. 29 as a function of the phase of the wave divided by the angular frequency, the phase velocity was obtained using a linear regression fit and was found to be 3.08 m/s. The radius of the vessel R=0.47 mm and the wall thickness h=0.12 mm were approximately estimated from the B-mode images, and the blood density was assumed to be 1060 kg/m3 (Cutnell, J. and W. Kenneth (1998). Physics, Fourth Edition. New York). Using these parameters, the Young's modulus of the aorta wall E=78.8 kPa was derived from the Moens-Korteweg equation (Eq. 5), which is what has been typically reported for thoracic aorta moduli in biomechanics literature (Fung, Y. C. (1993). Biomechanics—Mechanical Properties of Living Tissues. New York). - According to another aspect of the disclosed subject matter, systems and techniques are provided for electromechanical activation of arrhythmias, including non-transient electromechanical activation of paroxysmal and periodic arrhythmias in humans in vivo.
- Certain treatments of cardiac arrhythmias, such as radio-frequency ablation, can be utilized in clinical practice but can lack a suitable noninvasive imaging modality to provide insight into the source or focus of an arrhythmia. Cardiac deformations can be imaged at high temporal and spatial resolution to elucidate electrical activation sequences in normal and paced human subjects non-invasively. In this manner, such imaging can be used to improve planning and monitoring of ablation-based arrhythmia treatments.
- Aspects of the disclosed subject matter include techniques to quantitatively characterize focal and reentrant arrhythmias. For purpose of illustration and not limitation, and as embodied herein, spatio-temporal maps of a full-view of the atrial and ventricular mechanics can be obtained in a single heartbeat. Such maps can illustrate with suitable detail the electromechanical patterns of atrial flutter, fibrillation, and tachycardia. For example and without limitation, during focal arrhythmias, such as premature ventricular complex and focal atrial tachycardia, the electromechanical wave imaging techniques can be utilized to identify the location of the focal zone and the subsequent propagation of cardiac activation. For purpose of illustration and not limitation, exemplary electromechanical wave imaging techniques are described in International Application No. PCT/US13/64377, filed Oct. 10, 2013, which is incorporated by reference herein in its entirety. Additionally or alternatively, and as embodied herein, during reentrant arrhythmias, such as atrial flutter and fibrillation, Fourier analysis of the strains can show correlated mechanical and electrical cycle lengths and propagation patterns.
- For purpose of illustration and application of the disclosed subject matter, high frame rate ultrasound imaging of the heart can be used non-invasively and in real time to characterize lesser-known mechanical aspects of atrial and ventricular arrhythmias. Such techniques can also be used to assist treatment planning for intraoperative and longitudinal monitoring of arrhythmias.
- Certain imaging systems, such as software-based systems can allow ultra-high frame rates, and thus ultrasound imaging can be used to allow unprecedented temporal resolution. For example, such ultrasound imaging systems can provide a five-fold improvement in the signal-to-noise ratio of cardiac motion and deformation mapping. For example and without limitation, as embodied herein, frame rates up to 2000-5000 frames/s can be achieved by using defocussed transmissions, which can be suitable for depths utilized in transthoracic cardiac applications. According to exemplary embodiments of the disclosed subject matter, ultrasound imaging techniques described herein can be used to map transient strains occurring in response to the electrical activation, (e.g., the electromechanical wave). For example and without limitation, and as embodied herein, such techniques can be used to map transmural activation sequences of normal and abnormal hearts and to locate pacing sites in patients undergoing cardiac resynchronization therapy.
- According aspects of the disclosed subject matter, systems and techniques are provided to analyze and characterize the mechanical behavior of the atria. For purpose of illustration and not limitation, and as embodied herein, systems and techniques disclosed herein can be utilized to analyze and characterize the atria during certain types of cardiac arrhythmia, including and not limitation to, premature ventricular complex, focal tachycardia, atrial flutter, and atrial fibrillation. While Electromechanical Wave Imaging (EWI) can be suitable to characterize focal rhythms such as premature ventricular complex and focal tachycardia, EWI can have difficulty describing and/or characterizing reentrant rhythms such as atrial flutter and fibrillation. Accordingly, systems and techniques described herein are provided to characterize electromechanical strains, including and without limitation, during reentrant rhythms based on Fourier analysis. Exemplary embodiments of the disclosed subject matter can include a single acquisition sequence of electromechanical activation mapping that can be used for standard EWI and/or for Fourier analysis of electromechanical strains. Electromechanical activation mapping can characterize electromechanical strain propagation patterns during both focal and reentrant arrhythmias. In this manner, systems and techniques described herein can determine that local deformations of the atria can be closely correlated with their electrical activation. As such, systems and techniques described herein can be used to determine characteristics of cardiac mechanics in arrhythmia, to plan ablation treatments, and to monitor the efficacy of such treatments non-invasively, longitudinally and in real-time.
- For purpose of illustration and confirmation of the disclosed subject matter, exemplary techniques for imaging electromechanical activation of arrhythmias are described. The systems and techniques described herein can be performed, for purpose of illustration and not limitation, on human subjects. The human subjects can undergo a diagnostic ultrasound scan, and as embodied herein, can occur a few minutes to a few hours prior to electroanatomic mapping and ablation. The cardiac arrhythmias of the patients can be confirmed during electroanatomic mapping and ablation to be, for example and without limitation, one or more of premature ventricular complex (n=1), atrial flutter (n=5), focal atrial tachycardia (n=1), and atrial fibrillation (n=1). Additionally, a normal human subject can be imaged as a control for purpose of comparison.
- Additionally, and as embodied herein, strain maps can be generated, for example and without limitation using similar techniques as described herein for single-heartbeat electromechanical wave imaging (EWI). For example and without limitation, as embodied herein, a Verasonics system with a 2.5-MHz probe can be calibrated and customized to adhere to FDA standards, including measurements of mechanical index and of peak spatio-temporal-average intensity. The Verasonics system can be calibrated to have an acoustic power output that is similar to conventional clinical scanners. Such calibration can be performed by measuring the peak pressure and/or intensity (e.g., spatial-peak temporal average intensity, also referred to as Ispta) of the Verasonics system to ensure that its mechanical index (MI) is within FDA guidelines. The ultrasound scan can include two sequences. As embodied herein, in a motion-estimation sequence, a circular ultrasonic wave can be emitted with a virtual focus of 10.2 mm behind the probe at 2000 fps during 2 seconds. Additionally or alternatively, as embodied herein, a standard B-mode acquisition can be performed during 1.5 seconds to depict the heart anatomy. Frames from the motion-estimation sequence can be reconstructed by generating a plurality of beams, for example and as embodied herein 128 beams, in post-processing using a delay-and-sum algorithm with a reconstructed sampling frequency of 20 MHz. As embodied herein, the motion-estimation rate and the motion-sampling rate can be set to 1000 and 2000 fps, respectively. The window for motion-estimation can be 9.2 mm with an overlap of 95.8% (window shift of 0.3 mm), and the kernel strain estimation can be set to 4.9 mm. For purpose of illustration and not limitation, and as embodied herein, the techniques described herein for beamforming, motion-estimation, strain estimation, spatial moving-average of the strains (12 mm by 10 lines), and the automated contour tracking technique can be performed off-line on a graphics processing circuit (embodied herein as a Tesla graphics processing unit) and a Matlab parallel processing toolbox at a computing speed of 2.4 frames/s.
- Furthermore, and as embodied herein, focal and reentrant arrhythmias can be analyzed differently for patients with different types of rhythms. For purpose of illustration and not limitation,
FIG. 30A illustrates strains mapped in subjects having sinus rhythm.FIG. 30B illustrates strains mapped in subjects having atrial flutter.FIG. 30C illustrates strains mapped in subjects having atrial fibrillation. - As illustrated for example in
FIG. 30A , for subjects with sinus rhythm, the strains in one location, (e.g., one pixel in the left atrium) can present two representative events over time that correspond approximately to the beginning and the end of systole: end-systole, and end-diastole. By tracking the onset (e.g., the first zero-crossing) of these representative events for each pixel of the heart walls, isochrones maps can be generated. Isochrones correlated to electrical isochrones can be obtained, for example and without limitation, by tracking the propagation front of the end-diastole electromechanical activation. As illustrated for example inFIG. 30B , for atrial flutter patients, a similar location in the left atrium (LA) illustrates that strains can be periodic, which, in some cases, can be represented by a single frequency. Alternatively, as illustrated for example inFIG. 30C , a plurality of frequencies can be observed in a patient with atrial fibrillation, and as such, analysis based on the Fourier transform might can be utilized. During atrial fibrillation, the strains can be chaotic and no period of zero strains, similar to the one inFIG. 30A , can be observed. - In addition, and as embodied herein, the onset of contraction can be determined, for example in subjects who have focal rhythms such as sinus rhythm and focal tachycardia, as the first zero-crossing of the incremental strains occurring after the onset of the P-wave on the electrocardiogram (ECG), which can utilize the EWI techniques described herein. Additionally or alternatively, as embodied herein, in atria with reentrant arrhythmia, during flutter and fibrillation, a high-resolution Fourier transform can be performed using a generalized Goertzel algorithm for interpolation in Fourier space on 1.5-second long incremental strains signals for each individual pixel in the atria. For purpose of illustration and comparison with conventional ECG measurements, and without limitation, frequencies can be converted to cycle lengths, hereinafter referred to as mechanical cycle length (MCL). As embodied herein, peak MCL maps can be generated by selecting the MCL with the highest amplitude within the physiologically-relevant 100-330 ms range for each pixel. Additionally or alternatively, and as embodied herein, peak cycle lengths histograms can be constructed and compared to the electrical cycle length measured directly during the electroanatomic mapping and ablation.
-
FIG. 31 illustrates an EWI cine-loop and isochrones during focal rhythms.FIG. 31A illustrates the atria of a normal subject, with propagation from the right atrium (RA) to the LA. Theelectromechanical activation regions 3110 can originate in the right atrium and propagate towards the left atrium as illustrated in the exemplary EWI cine-loop inFIG. 31A . -
FIG. 31B illustrates an EWI cine-loop depicting the atria of a subject undergoing a focal atrial tachycardia, which can have a focus located high in the left atrium (LA). With reference toFIG. 31B , as embodied herein, electrical mapping of this patient has not been completed in the LA. The EWI inFIG. 31B illustrateselectromechanical activation 3120 originating high in the LA and propagating into both atria, and further activation can be detected in the ventricles. -
FIGS. 31C and 31D each illustrate isochrones obtained from a patient with ventricular tachycardia.FIG. 31C illustrates an isochrone of the ventricular tachycardia patient during sinus rhythm.FIG. 31D illustrates an isochrone of the ventricular tachycardia patient during premature ventricular complex. EWI was performed during sinus rhythm and during pre-ventricular contraction. The EWI isochrones obtained during sinus rhythm, as shown for example inFIG. 31C , depict propagation from the RA, into the LA and into the ventricles, as previously shown for purpose of illustration and comparison in normal patients. When this patient underwent premature ventricular complex, as illustrated for example inFIG. 31D , the region that was activated early in the ventricle during sinus rhythm (e.g., from the ventricles to the atria) triggered the entire electromechanical activation sequence. During premature ventricular complex, electromechanical activation can originate from the lateral wall, and can propagate toward the atria and into the atria. Early activation of the septum can indicate a potential recruitment of the Purkinje network. -
FIG. 32 illustrates the electromechanical behavior of a heart undergoing atrial flutter, including the analysis of reentrant arrhythmias using a single frequency flutter case.FIG. 32A illustrates an exemplary peak MCL map, and as depicted, a single MCL is representative. The peak cycle length map ofFIG. 32A indicates, for each pixel of the atria, which cycle length is representative in the Fourier spectrum.FIG. 32B illustrates a histogram of the cycle length which can be used to determine, among all the pixels of the atria, which cycle length represents atrial contraction. With reference toFIG. 32B , as embodied herein, one peak cycle length of 294 ms can be identified.FIG. 32C illustrates a phase map analyzing the phase of the MCL ofFIG. 32B in Fourier space. As shown for example inFIG. 32C , as embodied herein, a propagation pattern originates from the right atrium (RA) near the tricuspid valve towards the LA. The phase corresponding to the 294 ms cycle length can be retrieved from the Fourier spectrum and used to map the propagation of the mechanical oscillation at 294 ms. In this manner, the propagation direction can be determined. With reference toFIG. 32C , as embodied herein, the electromechanical activation propagated from the RA to the LA.FIG. 32D illustrates the corresponding intracardiac electrograms obtained a few hours after the imaging procedure. With reference toFIG. 32D , as embodied herein, the electrical cycle length was 283 ms. - Additionally or alternatively, and as embodied herein, atrial flutter cases can exhibit different patterns. Indeed, certain cases presented with two dominant frequencies can be separated between the left and right atria, whereas certain electrophysiological data can indicate that only one reentrant circuit was present.
-
FIG. 33A illustrates two examples of such atrial flutter cases.FIG. 33A illustrates peak MCL maps of two exemplary atrial flutter patients. As shown for example inFIG. 33A , two representative frequencies can be identified in each patient, with the shorter cycle length located in the RA. A relationship between the MCL and the electrical cycle length can be obtained, for example and as embodied herein, by performing this analysis in multiple patients, and choosing the peak cycle length closest to the electrical cycle length. For purpose of illustration and confirmation of the disclosed subject matter, such an analysis of MCL and electrical cycle length was conducted for five exemplary patients.FIG. 33B is a graph illustrating the results of the five exemplary patients. With reference toFIG. 33B , and as embodied herein, at least one representative MCL was very close to the electrical cycle length.FIG. 33B illustrates that electromechanical cycle length and/or MCL can be correlated with the electrical cycle length with a correlation of 0.96 and a slope of 1.1. -
FIGS. 33C and 33D together illustrate results from one exemplary patient undergoing atrial fibrillation.FIG. 33C illustrates a peak MCL map depicting multiple clustered dominant frequencies. For purpose of illustration and not limitation, the separation into these dominant frequencies is illustrated by the histogram shown inFIG. 33D . With reference toFIG. 33D , the peak cycle length map during atrial fibrillation depicts further spatial fragmentation of the peak cycle length. - Aspects of the present disclosed subject matter illustrate electromechanical activation mapping to identify the site of cardiac rhythm mechanisms during arrhythmia in humans and to characterize such cardiac rhythm mechanisms, which can lead to improved treatments and clinical management. Certain clinical practices utilize minimally invasive techniques to obtain precise maps of the activation of the atria and ventricles. Such techniques can be costly, time-consuming, and carry some degree of risk, and hence can be challenging to provide complete activation maps before and after treatment, as well as during catheter procedures.
- Aspects of the present disclosed subject matter illustrate methods for electromechanical activation mapping during reentrant and focal arrhythmias. For purpose of illustration and not limitation, and as embodied herein, exemplary techniques are provided for imaging spatiotemporal mechanics of arrhythmias with high accuracy and spatial and temporal resolutions in a full field of view in humans. These exemplary techniques can provide for characterization of an electromechanical propagation pattern and/or representative mechanical cycle lengths, which can correspond with their electrophysiological equivalents.
- For purpose of illustration and not limitation, and as embodied herein, focal rhythms can behave similarly to paced rhythms. For example and without limitation, focal rhythms can have a single source of electromechanical activation located in the vicinity of the earliest electrical activation. As embodied herein, EWI can be used to characterize the propagation of electromechanical activation, which can propagate from an atria's sinus node and from the bundle branch and which can terminate in the ventricles during ventricular pacing. Electromechanical activation propagation patterns similar to pacing can occur in a patient during premature ventricular complexes. The electromechanical activation sequence of the same patient during sinus rhythm can be similar to that of normal subjects. For example, and as embodied herein, in a patient with atrial tachycardia, the electromechanical activation propagation pattern can indicate a source located near the roof of the LA, in accordance with electrical mapping. As such, an exemplary application of non-invasive, ultrasound-based, electromechanical activation mapping is provided, which can be performed during or prior to invasive procedures. For purpose of illustration and not limitation, prior knowledge of an electromechanical source located in the LA can allow for clinical preparation. For example, and as embodied herein, such prior knowledge can be used to determine whether transseptal access would be obtained during treatment and to perform risk-benefit analysis to determine the best course of treatment (e.g., pharmacological vs. ablation treatment).
- Additionally, and as embodied herein, the electromechanical activation maps can be correlated with their electrical counterpart, at least in part of the atrial tissue, during atrial flutter. For example, and as embodied herein, a single representative frequency can be identified, and the phase of that frequency can indicate a propagation direction from the cavotricuspid isthmus region to the RA and LA, which can occur during typical atrial flutters. Additionally or alternatively, and as embodied herein, other behaviors can be identified in the atria, including and without limitation, one part of the atria that contracts with the same frequency as the electrical activation and another region that does not contract with the same frequency. As such, mapping the mechanics of the heart can identify regions of the heart in which the mechanical and electrical activities appear to be decoupled. Further spatial fragmentation of the periodicity of the mechanics of the atria can be observed during fibrillation. Such techniques can determine characteristics of the atrial mechanics during arrhythmia, including in the progression from flutter to fibrillation and vice versa.
- For purpose of illustration and not limitation, deformation of the atria caused by the onset of ventricular contraction and relaxation can affect certain aspects of the techniques described herein. For example, and as embodied herein, such a deformation can affect frequency analyses based on multiple activation cycles, which can be due at least in part to the relatively short acquisition time of these processes. Filtering and the development of longer acquisition sequences can inhibit or prevent such atrial deformation.
- Certain non-invasive electrical mapping techniques can be utilized to examine the epicardium, and can assume an immobilized heart function. A mechanical assessment of the atria can be utilized, for example, by electrophysiologists or interventional cardiologists to achieve the advantages described herein. According to certain clinical practices, echocardiograms can be performed on arrhythmia patients. Other non-invasive electrical mapping techniques can utilize on time-consuming and costly high resolution CT or MRI scans. The electromechanical activation mapping systems and techniques described herein can be obtained separately from, or in conjunction with, echocardiograms.
- As embodied herein, mapping the electromechanical activity during arrhythmias non-invasively with real-time feedback can be used determine characteristics of atrial mechanics in the evolution and perpetuation of arrhythmias. Furthermore, and as embodied herein, such a mapping can be used to predict the origin site of arrhythmias and the mechanism and monitoring of intervention outcomes.
- It will be understood that the foregoing is only illustrative of the principles of the present disclosure, and that various modifications can be made by those skilled in the art without departing from the scope and spirit of the present disclosure.
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US16/871,355 US20210085284A1 (en) | 2005-05-12 | 2020-05-11 | System and method for electromechanical activation of arrhythmias |
Applications Claiming Priority (6)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US68008105P | 2005-05-12 | 2005-05-12 | |
US11/433,510 US8858441B2 (en) | 2005-05-12 | 2006-05-12 | System and method for electromechanical wave imaging of body structures |
US14/449,820 US20140343424A1 (en) | 2005-05-12 | 2014-08-01 | System and method for electromechanical wave imaging of body structures |
US201562118402P | 2015-02-19 | 2015-02-19 | |
US15/048,761 US10687785B2 (en) | 2005-05-12 | 2016-02-19 | System and method for electromechanical activation of arrhythmias |
US16/871,355 US20210085284A1 (en) | 2005-05-12 | 2020-05-11 | System and method for electromechanical activation of arrhythmias |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/048,761 Division US10687785B2 (en) | 2005-05-12 | 2016-02-19 | System and method for electromechanical activation of arrhythmias |
Publications (1)
Publication Number | Publication Date |
---|---|
US20210085284A1 true US20210085284A1 (en) | 2021-03-25 |
Family
ID=74880339
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/048,761 Active 2029-07-16 US10687785B2 (en) | 2005-05-12 | 2016-02-19 | System and method for electromechanical activation of arrhythmias |
US16/871,355 Abandoned US20210085284A1 (en) | 2005-05-12 | 2020-05-11 | System and method for electromechanical activation of arrhythmias |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/048,761 Active 2029-07-16 US10687785B2 (en) | 2005-05-12 | 2016-02-19 | System and method for electromechanical activation of arrhythmias |
Country Status (1)
Country | Link |
---|---|
US (2) | US10687785B2 (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10653388B2 (en) * | 2014-09-25 | 2020-05-19 | Koninklijke Philips N.V. | Device and method for automatic pneumothorax detection |
WO2018170440A1 (en) | 2017-03-17 | 2018-09-20 | The Trustees Of Columbia University In The City Of New York | Non-invasive systems and methods for rendering of cardiac electromechanical activation |
US11304682B2 (en) * | 2018-09-28 | 2022-04-19 | University Of South Carolina | Non-invasive estimation of the mechanical properties of the heart |
US20210244340A1 (en) * | 2020-02-11 | 2021-08-12 | Catheter Precision, Inc. | Method of providing ventricular arrhythmia localization and myocardium wall thickness within a 3d heart model |
US12016728B2 (en) * | 2021-06-24 | 2024-06-25 | Biosense Webster (Israel) Ltd. | Estimating strain on tissue using 4D ultrasound catheter |
US20240115239A1 (en) * | 2022-10-11 | 2024-04-11 | Biosense Webster (Israel) Ltd. | Systems and methods for tracking an anatomical focal point |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006124603A2 (en) * | 2005-05-12 | 2006-11-23 | The Trustees Of Columbia University In The City Of New York | System and method for electromechanical wave imaging of body structures |
US20120271140A1 (en) * | 2011-04-22 | 2012-10-25 | Topera, Inc. | Basket style cardiac mapping catheter having an atraumatic basket tip for detection of cardiac rhythm disorders |
Family Cites Families (189)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3598111A (en) | 1968-12-09 | 1971-08-10 | Health Technology Corp | Technique and apparatus for measuring and monitoring the mechanical impedance of body tissues and organ systems |
US4463608A (en) | 1979-05-07 | 1984-08-07 | Yokogawa Hokushin Electric Corp. | Ultrasound imaging system |
US4777599A (en) | 1985-02-26 | 1988-10-11 | Gillette Company | Viscoelastometry of skin using shear wave propagation |
DE3529195A1 (en) | 1985-08-14 | 1987-02-26 | Max Planck Gesellschaft | CONTRAST AGENTS FOR ULTRASONIC EXAMINATIONS AND METHOD FOR THE PRODUCTION THEREOF |
US4822679A (en) | 1985-08-26 | 1989-04-18 | Stemcor Corporation | Spray-applied ceramic fiber insulation |
AU597002B2 (en) | 1985-10-28 | 1990-05-24 | Board Of Regents, The University Of Texas System | Pulse centroid echo method and apparatus for enhanced sound velocity estimation in vivo |
US4858613A (en) | 1988-03-02 | 1989-08-22 | Laboratory Equipment, Corp. | Localization and therapy system for treatment of spatially oriented focal disease |
US4926675A (en) | 1988-03-18 | 1990-05-22 | Tennessee Valley Authority | Air detector for liquid-filled sensing lines |
US5038787A (en) | 1988-08-10 | 1991-08-13 | The Board Of Regents, The University Of Texas System | Method and apparatus for analyzing material properties using reflected ultrasound |
US5107837A (en) | 1989-11-17 | 1992-04-28 | Board Of Regents, University Of Texas | Method and apparatus for measurement and imaging of tissue compressibility or compliance |
US5457754A (en) | 1990-08-02 | 1995-10-10 | University Of Cincinnati | Method for automatic contour extraction of a cardiac image |
US6405072B1 (en) | 1991-01-28 | 2002-06-11 | Sherwood Services Ag | Apparatus and method for determining a location of an anatomical target with reference to a medical apparatus |
JP3109749B2 (en) | 1991-04-17 | 2000-11-20 | 株式会社東芝 | Ultrasound imaging device |
EP0766533A1 (en) | 1991-05-17 | 1997-04-09 | InnerDyne, Inc. | Method and device for thermal ablation |
JPH06511481A (en) | 1991-07-05 | 1994-12-22 | ユニバーシティ オブ ロチェスター | Ultra-fine non-agglomerated porous particles that incorporate air bubbles |
DE69431741T2 (en) | 1993-03-12 | 2003-09-11 | Toshiba Kawasaki Kk | Device for medical treatment with ultrasound |
US5435310A (en) | 1993-06-23 | 1995-07-25 | University Of Washington | Determining cardiac wall thickness and motion by imaging and three-dimensional modeling |
US5601084A (en) | 1993-06-23 | 1997-02-11 | University Of Washington | Determining cardiac wall thickness and motion by imaging and three-dimensional modeling |
US20020169394A1 (en) | 1993-11-15 | 2002-11-14 | Eppstein Jonathan A. | Integrated tissue poration, fluid harvesting and analysis device, and method therefor |
US6259943B1 (en) | 1995-02-16 | 2001-07-10 | Sherwood Services Ag | Frameless to frame-based registration system |
US5718241A (en) | 1995-06-07 | 1998-02-17 | Biosense, Inc. | Apparatus and method for treating cardiac arrhythmias with no discrete target |
US5662113A (en) | 1995-06-30 | 1997-09-02 | Siemens Medical Systems, Inc | Edge enhancement system for ultrasound images |
US6351659B1 (en) | 1995-09-28 | 2002-02-26 | Brainlab Med. Computersysteme Gmbh | Neuro-navigation system |
US5810731A (en) | 1995-11-13 | 1998-09-22 | Artann Laboratories | Method and apparatus for elasticity imaging using remotely induced shear wave |
US5606971A (en) | 1995-11-13 | 1997-03-04 | Artann Corporation, A Nj Corp. | Method and device for shear wave elasticity imaging |
WO1997029783A1 (en) | 1996-02-19 | 1997-08-21 | Nycomed Imaging A/S | Improvements in or relating to contrast agents |
JP4237256B2 (en) | 1996-02-29 | 2009-03-11 | シーメンス メディカル ソリューションズ ユーエスエイ インコーポレイテッド | Ultrasonic transducer |
JP3652791B2 (en) | 1996-06-24 | 2005-05-25 | 独立行政法人科学技術振興機構 | Ultrasonic diagnostic equipment |
US6026173A (en) | 1997-07-05 | 2000-02-15 | Svenson; Robert H. | Electromagnetic imaging and therapeutic (EMIT) systems |
US5752515A (en) | 1996-08-21 | 1998-05-19 | Brigham & Women's Hospital | Methods and apparatus for image-guided ultrasound delivery of compounds through the blood-brain barrier |
US8003705B2 (en) | 1996-09-23 | 2011-08-23 | Incept Llc | Biocompatible hydrogels made with small molecule precursors |
US5769790A (en) | 1996-10-25 | 1998-06-23 | General Electric Company | Focused ultrasound surgery system guided by ultrasound imaging |
US6090800A (en) | 1997-05-06 | 2000-07-18 | Imarx Pharmaceutical Corp. | Lipid soluble steroid prodrugs |
DE69838669T2 (en) | 1997-04-30 | 2008-10-30 | Point Biomedical Corp., San Carlos | MICROPARTICLES, SUITABLE AS A CONTRASTANT IN ULTRASOUND AND FOR ACTIVITY IN THE BLOOD CIRCUIT |
US6102864A (en) | 1997-05-07 | 2000-08-15 | General Electric Company | Three-dimensional ultrasound imaging of velocity and power data using average or median pixel projections |
US20020039594A1 (en) | 1997-05-13 | 2002-04-04 | Evan C. Unger | Solid porous matrices and methods of making and using the same |
US6287591B1 (en) | 1997-05-14 | 2001-09-11 | Inex Pharmaceuticals Corp. | Charged therapeutic agents encapsulated in lipid particles containing four lipid components |
US6095976A (en) | 1997-06-19 | 2000-08-01 | Medinol Ltd. | Method for enhancing an image derived from reflected ultrasound signals produced by an ultrasound transmitter and detector inserted in a bodily lumen |
US6106465A (en) | 1997-08-22 | 2000-08-22 | Acuson Corporation | Ultrasonic method and system for boundary detection of an object of interest in an ultrasound image |
US5928151A (en) | 1997-08-22 | 1999-07-27 | Acuson Corporation | Ultrasonic system and method for harmonic imaging in three dimensions |
DE19746060A1 (en) | 1997-10-17 | 1999-04-22 | Itt Mfg Enterprises Inc | Piston for hydraulic piston cylinder unit, especially for disc brake |
US6896659B2 (en) | 1998-02-06 | 2005-05-24 | Point Biomedical Corporation | Method for ultrasound triggered drug delivery using hollow microbubbles with controlled fragility |
US6200266B1 (en) | 1998-03-31 | 2001-03-13 | Case Western Reserve University | Method and apparatus for ultrasound imaging using acoustic impedance reconstruction |
JP4260920B2 (en) | 1998-05-13 | 2009-04-30 | 株式会社東芝 | Ultrasonic diagnostic equipment |
US6287765B1 (en) | 1998-05-20 | 2001-09-11 | Molecular Machines, Inc. | Methods for detecting and identifying single molecules |
US6270459B1 (en) | 1998-05-26 | 2001-08-07 | The Board Of Regents Of The University Of Texas System | Method for estimating and imaging of transverse displacements, transverse strains and strain ratios |
US6511426B1 (en) | 1998-06-02 | 2003-01-28 | Acuson Corporation | Medical diagnostic ultrasound system and method for versatile processing |
US6241675B1 (en) | 1998-06-09 | 2001-06-05 | Volumetrics Medical Imaging | Methods and systems for determining velocity of tissue using three dimensional ultrasound data |
US6425865B1 (en) | 1998-06-12 | 2002-07-30 | The University Of British Columbia | Robotically assisted medical ultrasound |
ATE253126T1 (en) | 1998-08-25 | 2003-11-15 | Univ Washington | RAPID QUANTITATIVE ANALYSIS OF PROTEINS OR PROTEIN FUNCTIONS IN COMPLEX MIXTURES |
US6425867B1 (en) | 1998-09-18 | 2002-07-30 | University Of Washington | Noise-free real time ultrasonic imaging of a treatment site undergoing high intensity focused ultrasound therapy |
US6246895B1 (en) | 1998-12-18 | 2001-06-12 | Sunnybrook Health Science Centre | Imaging of ultrasonic fields with MRI |
US6309355B1 (en) | 1998-12-22 | 2001-10-30 | The Regents Of The University Of Michigan | Method and assembly for performing ultrasound surgery using cavitation |
US6547730B1 (en) | 1998-12-31 | 2003-04-15 | U-Systems, Inc. | Ultrasound information processing system |
FR2791136B1 (en) | 1999-03-15 | 2001-06-08 | Mathias Fink | IMAGING METHOD AND DEVICE USING SHEAR WAVES |
AU768759B2 (en) | 1999-06-14 | 2004-01-08 | Exogen, Inc. | Method and kit for cavitation-induced tissue healing with low intensity ultrasound |
GB9916851D0 (en) | 1999-07-20 | 1999-09-22 | Univ Wales Bangor | Manipulation of particles in liquid media |
US6352507B1 (en) | 1999-08-23 | 2002-03-05 | G.E. Vingmed Ultrasound As | Method and apparatus for providing real-time calculation and display of tissue deformation in ultrasound imaging |
US6775400B1 (en) | 1999-10-29 | 2004-08-10 | Acuson Corporation | Medical diagnostic ultrasonic imaging method and apparatus for suppressing electronic noise |
US6447450B1 (en) | 1999-11-02 | 2002-09-10 | Ge Medical Systems Global Technology Company, Llc | ECG gated ultrasonic image compounding |
US6312382B1 (en) | 1999-11-15 | 2001-11-06 | Ronald Mucci | Method and apparatus for extracting cardiac information from acoustic information acquired with an ultrasound device |
US7374538B2 (en) | 2000-04-05 | 2008-05-20 | Duke University | Methods, systems, and computer program products for ultrasound measurements using receive mode parallel processing |
EP1313397A2 (en) | 2000-07-19 | 2003-05-28 | Nigel E. Sharrock | Non-invasive measurement of suprasystolic signals |
US6514221B2 (en) | 2000-07-27 | 2003-02-04 | Brigham And Women's Hospital, Inc. | Blood-brain barrier opening |
US6529770B1 (en) | 2000-11-17 | 2003-03-04 | Valentin Grimblatov | Method and apparatus for imaging cardiovascular surfaces through blood |
US6821274B2 (en) | 2001-03-07 | 2004-11-23 | Gendel Ltd. | Ultrasound therapy for selective cell ablation |
US6508768B1 (en) | 2000-11-22 | 2003-01-21 | University Of Kansas Medical Center | Ultrasonic elasticity imaging |
US7022077B2 (en) | 2000-11-28 | 2006-04-04 | Allez Physionix Ltd. | Systems and methods for making noninvasive assessments of cardiac tissue and parameters |
CA2428872C (en) | 2000-11-28 | 2013-01-08 | Allez Physionix Limited | Systems and methods for making non-invasive physiological assessments |
US6671541B2 (en) | 2000-12-01 | 2003-12-30 | Neomed Technologies, Inc. | Cardiovascular imaging and functional analysis system |
US6537221B2 (en) | 2000-12-07 | 2003-03-25 | Koninklijke Philips Electronics, N.V. | Strain rate analysis in ultrasonic diagnostic images |
US6491636B2 (en) | 2000-12-07 | 2002-12-10 | Koninklijke Philips Electronics N.V. | Automated border detection in ultrasonic diagnostic images |
WO2002064205A2 (en) | 2001-02-13 | 2002-08-22 | Quetzal Biomedical, Inc. | Multi-electrode apparatus and method for treatment of congestive heart failure |
US6689060B2 (en) | 2001-02-28 | 2004-02-10 | Siemens Medical Solutions Usa, Inc | System and method for re-orderable nonlinear echo processing |
US6488629B1 (en) | 2001-07-31 | 2002-12-03 | Ge Medical Systems Global Technology Company, Llc | Ultrasound image acquisition with synchronized reference image |
EP1983043A1 (en) | 2001-08-14 | 2008-10-22 | Washington University in St. Louis | Systems and methods for screening pharmaceutical chemicals |
US7192403B2 (en) | 2001-08-17 | 2007-03-20 | Russell Ted W | Methods, apparatus and articles-of-manufacture for noninvasive measurement and monitoring of peripheral blood flow, perfusion, cardiac output biophysic stress and cardiovascular condition |
US6537217B1 (en) | 2001-08-24 | 2003-03-25 | Ge Medical Systems Global Technology Company, Llc | Method and apparatus for improved spatial and temporal resolution in ultrasound imaging |
WO2003028801A2 (en) | 2001-10-04 | 2003-04-10 | Case Western Reserve University | Systems and methods for noninvasive electrocardiographic imaging (ecgi) using generalized minimum residual (gmres) |
FR2830936B1 (en) | 2001-10-16 | 2004-08-27 | Agronomique Inst Nat Rech | METHOD FOR MEASURING THE STATE OF TENSION OF A MATERIAL AND APPLICATIONS OF THIS METHOD |
US6855114B2 (en) | 2001-11-23 | 2005-02-15 | Karen Drukker | Automated method and system for the detection of abnormalities in sonographic images |
US6685641B2 (en) | 2002-02-01 | 2004-02-03 | Siemens Medical Solutions Usa, Inc. | Plane wave scanning reception and receiver |
US7166075B2 (en) | 2002-03-08 | 2007-01-23 | Wisconsin Alumni Research Foundation | Elastographic imaging of in vivo soft tissue |
CA2421468A1 (en) | 2002-03-14 | 2003-09-14 | Matsushita Electric Industrial Co., Ltd. | Image processing device and ultrasonic diagnostic device |
US6683454B2 (en) | 2002-03-28 | 2004-01-27 | Ge Medical Systems Global Technology Company, Llc | Shifting of artifacts by reordering of k-space |
US20030220556A1 (en) | 2002-05-20 | 2003-11-27 | Vespro Ltd. | Method, system and device for tissue characterization |
US7819806B2 (en) | 2002-06-07 | 2010-10-26 | Verathon Inc. | System and method to identify and measure organ wall boundaries |
US7043292B2 (en) | 2002-06-21 | 2006-05-09 | Tarjan Peter P | Single or multi-mode cardiac activity data collection, processing and display obtained in a non-invasive manner |
US6780152B2 (en) | 2002-06-26 | 2004-08-24 | Acuson Corporation | Method and apparatus for ultrasound imaging of the heart |
US7549985B2 (en) | 2002-06-26 | 2009-06-23 | The Regents Of The University Of Michigan | Method and system to create and acoustically manipulate a microbubble |
US20040049134A1 (en) | 2002-07-02 | 2004-03-11 | Tosaya Carol A. | System and methods for treatment of alzheimer's and other deposition-related disorders of the brain |
US7136518B2 (en) | 2003-04-18 | 2006-11-14 | Medispectra, Inc. | Methods and apparatus for displaying diagnostic data |
US7314446B2 (en) | 2002-07-22 | 2008-01-01 | Ep Medsystems, Inc. | Method and apparatus for time gating of medical images |
FR2843874B1 (en) | 2002-08-28 | 2004-11-05 | Centre Nat Rech Scient | NON-INVASIVE METHOD FOR OBTAINING A PREDETERMINED FIELD OF ACOUSTIC WAVES IN A SUBSTANTIALLY HOMOGENEOUS ENVIRONMENT MASKED BY A BONE BARRIER, IMAGING METHOD, AND DEVICE FOR CARRYING OUT SAID METHODS |
US8560063B2 (en) | 2002-09-10 | 2013-10-15 | Uab Research Foundation | Post-defibrillation pacing methods and devices |
US6749571B2 (en) | 2002-09-19 | 2004-06-15 | Wisconsin Alumni Research Foundation | Method and apparatus for cardiac elastography |
US7103400B2 (en) | 2002-11-08 | 2006-09-05 | Koninklijke Philips Electronics, N.V. | Artifact elimination in time-gated anatomical imaging |
US7697972B2 (en) | 2002-11-19 | 2010-04-13 | Medtronic Navigation, Inc. | Navigation system for cardiac therapies |
US7074187B2 (en) | 2002-12-13 | 2006-07-11 | Selzer Robert H | System and method for improving ultrasound image acquisition and replication for repeatable measurements of vascular structures |
JP3872424B2 (en) | 2002-12-20 | 2007-01-24 | アロカ株式会社 | Ultrasonic diagnostic equipment |
CA2552690C (en) | 2002-12-31 | 2014-12-09 | Ultra-Sonic Technologies, L.L.C. | Transdermal delivery using encapsulated agent activated by ultrasound and/or heat |
US6994673B2 (en) | 2003-01-16 | 2006-02-07 | Ge Ultrasound Israel, Ltd | Method and apparatus for quantitative myocardial assessment |
US7257244B2 (en) | 2003-02-24 | 2007-08-14 | Vanderbilt University | Elastography imaging modalities for characterizing properties of tissue |
US20040172081A1 (en) | 2003-02-28 | 2004-09-02 | Dai-Yuan Wang | Intracardiac pressure guided pacemaker |
US20040258760A1 (en) | 2003-03-20 | 2004-12-23 | Wheatley Margaret A. | Isolated nanocapsule populations and surfactant-stabilized microcapsules and nanocapsules for diagnostic imaging and drug delivery and methods for their production |
US7175599B2 (en) | 2003-04-17 | 2007-02-13 | Brigham And Women's Hospital, Inc. | Shear mode diagnostic ultrasound |
US7344509B2 (en) | 2003-04-17 | 2008-03-18 | Kullervo Hynynen | Shear mode therapeutic ultrasound |
US7601122B2 (en) | 2003-04-22 | 2009-10-13 | Wisconsin Alumni Research Foundation | Ultrasonic elastography with angular compounding |
US7052460B2 (en) | 2003-05-09 | 2006-05-30 | Visualsonics Inc. | System for producing an ultrasound image using line-based image reconstruction |
US20050277835A1 (en) | 2003-05-30 | 2005-12-15 | Angelsen Bjorn A | Ultrasound imaging by nonlinear low frequency manipulation of high frequency scattering and propagation properties |
CN100481096C (en) | 2003-06-25 | 2009-04-22 | 美国西门子医疗解决公司 | Automated regional myocardial assessment method for cardiac imaging |
US6984209B2 (en) | 2003-07-02 | 2006-01-10 | The Brigham And Women's Hospital, Inc. | Harmonic motion imaging |
US20050026262A1 (en) | 2003-07-30 | 2005-02-03 | Sonoenergy, Llc | Sonication-enhanced digestion process |
US7055378B2 (en) | 2003-08-11 | 2006-06-06 | Veeco Instruments, Inc. | System for wide frequency dynamic nanomechanical analysis |
US7358226B2 (en) | 2003-08-27 | 2008-04-15 | The Regents Of The University Of California | Ultrasonic concentration of drug delivery capsules |
US20050080469A1 (en) | 2003-09-04 | 2005-04-14 | Larson Eugene A. | Treatment of cardiac arrhythmia utilizing ultrasound |
US20050054930A1 (en) | 2003-09-09 | 2005-03-10 | The University Court Of The University Of Dundee | Sonoelastography using power Doppler |
ES2252649T3 (en) | 2003-09-30 | 2006-05-16 | Esaote S.P.A. | METHOD FOR ESTIMATING TISSUE SPEED VECTORS AND TISSUE DEFORMATION FROM DIAGNOSTIC IMAGE FORMATION DATA BY ULTRASOUNDS. |
US7421101B2 (en) | 2003-10-02 | 2008-09-02 | Siemens Medical Solutions Usa, Inc. | System and method for local deformable motion analysis |
US7753847B2 (en) | 2003-10-03 | 2010-07-13 | Mayo Foundation For Medical Education And Research | Ultrasound vibrometry |
US7896821B1 (en) | 2003-11-14 | 2011-03-01 | Perfusion Technology, LLC | Low intensity directed ultrasound (LODUS) mediated blood brain barrier disruption |
US20050175541A1 (en) | 2003-11-19 | 2005-08-11 | Lanza Gregory M. | Enhanced drug delivery |
US7318804B2 (en) | 2003-12-09 | 2008-01-15 | The Regents Of The University Of Michigan | Methods and systems for measuring mechanical property of a vascular wall and method and system for determining health of a vascular structure |
JPWO2005055831A1 (en) | 2003-12-10 | 2007-12-06 | 松下電器産業株式会社 | Ultrasonic diagnostic apparatus and ultrasonic diagnostic method |
ATE466596T1 (en) | 2004-01-20 | 2010-05-15 | Sunnybrook & Womens College | HIGH FREQUENCY ULTRASONIC IMAGING WITH CONTRAST AGENTS |
US7951083B2 (en) | 2004-02-05 | 2011-05-31 | Siemens Medical Solutions Usa, Inc. | Motion analysis improvements for medical diagnostic ultrasound |
WO2005098731A2 (en) | 2004-03-29 | 2005-10-20 | German Peter T | Systems and methods to determine elastic properties of materials |
US7809426B2 (en) | 2004-04-29 | 2010-10-05 | The Cleveland Clinic Foundation | Acquiring contrast-enhanced, T1 weighted, cine magnetic resonance images |
US7372984B2 (en) | 2004-05-05 | 2008-05-13 | California Institute Of Technology | Four-dimensional imaging of periodically moving objects via post-acquisition synchronization of nongated slice-sequences |
JP4465535B2 (en) | 2004-06-09 | 2010-05-19 | 株式会社日立メディコ | Elastic image display method and ultrasonic diagnostic apparatus |
JP2008509890A (en) | 2004-08-05 | 2008-04-03 | ベイラー・リサーチ・インスチチユート | Gene or drug delivery system |
US7699780B2 (en) | 2004-08-11 | 2010-04-20 | Insightec—Image-Guided Treatment Ltd. | Focused ultrasound system with adaptive anatomical aperture shaping |
US20060058651A1 (en) | 2004-08-13 | 2006-03-16 | Chiao Richard Y | Method and apparatus for extending an ultrasound image field of view |
US7678050B2 (en) | 2004-08-24 | 2010-03-16 | General Electric Company | Method and apparatus for detecting cardiac events |
US20060074315A1 (en) | 2004-10-04 | 2006-04-06 | Jianming Liang | Medical diagnostic ultrasound characterization of cardiac motion |
WO2006044996A2 (en) | 2004-10-15 | 2006-04-27 | The Trustees Of Columbia University In The City Of New York | System and method for automated boundary detection of body structures |
WO2006044997A2 (en) | 2004-10-15 | 2006-04-27 | The Trustees Of Columbia University In The City Of New York | System and method for localized measurement and imaging of viscosity of tissues |
US7223241B2 (en) | 2004-12-16 | 2007-05-29 | Aloka Co., Ltd. | Method and apparatus for elasticity imaging |
US7674229B2 (en) | 2005-03-07 | 2010-03-09 | The Brigham And Women's Hospital, Inc. | Adaptive ultrasound delivery system |
US20060241462A1 (en) | 2005-03-14 | 2006-10-26 | Yi-Hong Chou | Method of intracranial ultrasound imaging and related system |
US20070129652A1 (en) | 2005-11-15 | 2007-06-07 | Henry Nita | Methods and apparatus for intracranial ultrasound therapies |
US20110313328A1 (en) | 2005-06-24 | 2011-12-22 | Penumbra, Inc. | Methods and apparatus for dissolving blockages in intracranial catheters |
WO2007027584A2 (en) | 2005-08-30 | 2007-03-08 | University Of Virginia Patent Foundation | Deposit contrast agents and related methods thereof |
US7967763B2 (en) | 2005-09-07 | 2011-06-28 | Cabochon Aesthetics, Inc. | Method for treating subcutaneous tissues |
WO2007035721A2 (en) | 2005-09-19 | 2007-03-29 | The Trustees Of Columbia University In The City Of New York | Ultrasound method to open blood brain barrier |
US20070071683A1 (en) | 2005-09-27 | 2007-03-29 | The Regents Of The University Of California | Ultrasonic concentration of carrier particles |
US8257338B2 (en) | 2006-10-27 | 2012-09-04 | Artenga, Inc. | Medical microbubble generation |
JP5630958B2 (en) | 2005-11-02 | 2014-11-26 | ビジュアルソニックス インコーポレイテッド | High frequency array ultrasound system |
WO2007067987A2 (en) | 2005-12-09 | 2007-06-14 | The Trustees Of Columbia University In The City Ofnew York | Systems and methods for elastography imaging |
WO2007100731A2 (en) | 2006-02-24 | 2007-09-07 | Nanovibronix Inc. | System and method for surface acoustic wave treatment of medical devices |
JP4732925B2 (en) | 2006-03-09 | 2011-07-27 | イマグノーシス株式会社 | Medical image display method and program thereof |
ATE479166T1 (en) | 2006-06-23 | 2010-09-15 | Koninkl Philips Electronics Nv | METHOD, APPARATUS AND COMPUTER PROGRAM FOR THREE-DIMENSIONAL ULTRASOUND IMAGING |
ITBO20060593A1 (en) | 2006-08-04 | 2008-02-05 | Francesca Cavalieri | MICROBULTS MADE OF POLYVINYL ALCOHOL AND RELATIVE LOADING OF THE SAME WITH NITROGEN OXIDE |
WO2008027520A2 (en) | 2006-08-30 | 2008-03-06 | The Trustees Of Columbia University In The City Of New York | Systems and methods for composite elastography and wave imaging |
US20080081990A1 (en) * | 2006-09-28 | 2008-04-03 | The Research Foundation Of State University Of New York | Apparatus, system, kit and method for heart mapping |
US20100056924A1 (en) | 2006-11-20 | 2010-03-04 | Koninklijke Philips Electronics N.V. | Control and display of ultrasonic microbubble cavitation |
US8574157B2 (en) | 2007-02-14 | 2013-11-05 | General Electric Company | Method and apparatus for generating an ultrasound image of moving objects using deformable models |
US20080243214A1 (en) | 2007-03-26 | 2008-10-02 | Boston Scientific Scimed, Inc. | High resolution electrophysiology catheter |
WO2008131302A2 (en) | 2007-04-19 | 2008-10-30 | The Foundry, Inc. | Methods and apparatus for reducing sweat production |
CA2684752A1 (en) | 2007-04-20 | 2008-10-30 | Oregon Health & Science University | Ultrasound imaging with targeted microbubbles |
US20080311045A1 (en) | 2007-06-06 | 2008-12-18 | Biovaluation & Analysis, Inc. | Polymersomes for Use in Acoustically Mediated Intracellular Drug Delivery in vivo |
WO2008157422A1 (en) | 2007-06-13 | 2008-12-24 | Charles Thomas Hardy | Materials, methods, and systems for cavitation-mediated ultrasonic drug delivery |
US20080319355A1 (en) | 2007-06-20 | 2008-12-25 | Henry Nita | Ischemic stroke therapy |
WO2009043031A2 (en) | 2007-09-27 | 2009-04-02 | Children's Medical Center Corporation | Microbubbles and methods for oxygen delivery |
US20090247911A1 (en) | 2008-03-25 | 2009-10-01 | Petr Novak | Multiple-angle switched high intensity focused ultrasound |
US8208709B2 (en) | 2008-04-17 | 2012-06-26 | The Ohio State University Research Foundation | System and method for improved real-time cine imaging |
US8545405B2 (en) | 2008-04-23 | 2013-10-01 | Therataxis, Llc | Device, methods, and control for sonic guidance of molecules and other material utilizing time-reversal acoustics |
US9089278B2 (en) | 2008-07-10 | 2015-07-28 | Koninklijke Philips N.V. | Ultrasonic assessment of cardiac synchronicity and viability |
WO2010014977A1 (en) * | 2008-08-01 | 2010-02-04 | The Trustees Of Columbia University In The City Of New York | Systems and methods for matching and imaging tissue characteristics |
WO2010030819A1 (en) | 2008-09-10 | 2010-03-18 | The Trustees Of Columbia University In The City Of New York | Systems and methods for opening a tissue |
WO2010044385A1 (en) | 2008-10-14 | 2010-04-22 | 株式会社 日立メディコ | Ultrasonographic device and ultrasonographic display method |
US20100143241A1 (en) | 2008-10-22 | 2010-06-10 | Johnson G Allan | Method and apparatus for delivery of agents across the blood brain barrier |
FR2939512B1 (en) | 2008-12-04 | 2012-07-27 | Echosens | DEVICE AND METHOD FOR ELASTOGRAPHY |
US20100286527A1 (en) | 2009-05-08 | 2010-11-11 | Penrith Corporation | Ultrasound system with multi-head wireless probe |
US8628477B2 (en) | 2009-07-31 | 2014-01-14 | Nellcor Puritan Bennett Ireland | Systems and methods for non-invasive determination of blood pressure |
WO2011028690A1 (en) | 2009-09-01 | 2011-03-10 | The Trustees Of Columbia University In The City Of New York | Microbubble devices, methods and systems |
EP2480144B1 (en) | 2009-09-21 | 2024-03-06 | The Trustees of Columbia University in the City of New York | Systems for opening of a tissue barrier |
WO2011079177A1 (en) | 2009-12-22 | 2011-06-30 | The Trustees Of Columbia University In The City Of New York | A planning system for targeting tissue structures with ultrasound |
US8709451B2 (en) | 2010-01-20 | 2014-04-29 | University Of Utah Research Foundation | Stable nanoemulsions for ultrasound-mediated drug delivery and imaging |
US9256947B2 (en) | 2010-03-19 | 2016-02-09 | Koninklijke Philips N.V. | Automatic positioning of imaging plane in ultrasonic imaging |
US8717210B2 (en) | 2010-04-27 | 2014-05-06 | Technion Research & Development Foundation Ltd. | Multi-channel sampling of pulse streams at the rate of innovation |
US9119545B2 (en) * | 2010-05-27 | 2015-09-01 | St. Jude Medical, AB | Arrhythmia classification |
US9585631B2 (en) | 2010-06-01 | 2017-03-07 | The Trustees Of Columbia University In The City Of New York | Devices, methods, and systems for measuring elastic properties of biological tissues using acoustic force |
US8380305B2 (en) | 2010-07-01 | 2013-02-19 | DynaDx Corporation | System and method for predicting successful defibrillation for ventricular fibrillation cardiac arrest |
EP2600771A1 (en) | 2010-08-06 | 2013-06-12 | The Trustees of Columbia University in the City of New York | Medical imaging contrast devices, methods, and systems |
WO2012162664A1 (en) | 2011-05-26 | 2012-11-29 | The Trustees Of Columbia University In The City Of New York | Systems and methods for opening of a tissue barrier in primates |
JP5865050B2 (en) | 2011-12-15 | 2016-02-17 | キヤノン株式会社 | Subject information acquisition device |
CN102608212B (en) | 2012-02-29 | 2013-10-30 | 大连理工大学 | Method for measuring acoustic impedance and acoustic attenuation of thin layer based on sound pressure reflection coefficient power spectrum |
US9470728B2 (en) | 2012-05-09 | 2016-10-18 | Cardioinsight Technologies, Inc. | Channel integrity detection |
-
2016
- 2016-02-19 US US15/048,761 patent/US10687785B2/en active Active
-
2020
- 2020-05-11 US US16/871,355 patent/US20210085284A1/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006124603A2 (en) * | 2005-05-12 | 2006-11-23 | The Trustees Of Columbia University In The City Of New York | System and method for electromechanical wave imaging of body structures |
US20120271140A1 (en) * | 2011-04-22 | 2012-10-25 | Topera, Inc. | Basket style cardiac mapping catheter having an atraumatic basket tip for detection of cardiac rhythm disorders |
Also Published As
Publication number | Publication date |
---|---|
US10687785B2 (en) | 2020-06-23 |
US20160249880A1 (en) | 2016-09-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8858441B2 (en) | System and method for electromechanical wave imaging of body structures | |
US20210085284A1 (en) | System and method for electromechanical activation of arrhythmias | |
US8428687B2 (en) | Systems and methods for matching and imaging tissue characteristics | |
Hsu et al. | In vivo assessment of myocardial stiffness with acoustic radiation force impulse imaging | |
JP3875581B2 (en) | Ultrasound diagnostic system | |
Horton et al. | Assessment of the right ventricle by echocardiography: a primer for cardiac sonographers | |
US20200163646A1 (en) | Systems and methods for mechanical mapping of cardiac rhythm | |
US11678859B2 (en) | Systems and methods for rendering of cardiac electromechanical activation | |
Provost et al. | Electromechanical wave imaging for arrhythmias | |
Provost et al. | Electromechanical wave imaging of normal and ischemic hearts in vivo | |
US20190167233A1 (en) | Methods and systems for ultrasound elastography | |
Provost et al. | Assessing the atrial electromechanical coupling during atrial focal tachycardia, flutter, and fibrillation using electromechanical wave imaging in humans | |
Hollender et al. | Intracardiac acoustic radiation force impulse (ARFI) and shear wave imaging in pigs with focal infarctions | |
Konofagou et al. | Noninvasive electromechanical wave imaging and conduction-relevant velocity estimation in vivo | |
Pernot et al. | Electromechanical imaging of the myocardium at normal and pathological states | |
Kakkad et al. | Non-invasive measurement of dynamic myocardial stiffness using acoustic radiation force impulse imaging | |
Hsu et al. | Novel acoustic radiation force impulse imaging methods for visualization of rapidly moving tissue | |
Engel et al. | Cardiac atrial kick shear wave elastography with ultrafast diverging wave imaging: An in vivo pilot study | |
Konofagou et al. | 11B-1 Noninvasive Electromechanical Wave Imaging and Conduction Velocity Estimation In Vivo | |
Eyerly‐Webb et al. | Acoustic Radiation Force‐based Ultrasound Elastography for Cardiac Imaging Applications | |
Bradway et al. | 6B-6 Transthoracic cardiac acoustic radiation force impulse imaging: a feasibility study | |
Robert et al. | Spectral Analysis of Tissue Displacement for Cardiac Activation Mapping: Ex Vivo Working Heart and In Vivo Study | |
Luo et al. | 2I-6 Imaging the Electromechanical Wave Activation of the Left Ventricle in Vivo | |
Konofagou et al. | Evolving concepts in measuring ventricular strain in canine and human hearts: non-invasive imaging | |
Hsu | Acoustic radiation force impulse imaging of myocardial performance |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: NATIONAL INSTITUTES OF HEALTH (NIH), U.S. DEPT. OF HEALTH AND HUMAN SERVICES (DHHS), U.S. GOVERNMENT, MARYLAND Free format text: CONFIRMATORY LICENSE;ASSIGNOR:COLUMBIA UNIV NEW YORK MORNINGSIDE;REEL/FRAME:052819/0787 Effective date: 20200601 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: APPLICATION DISPATCHED FROM PREEXAM, NOT YET DOCKETED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |