IP Library Granted Patent US 12,217,434
Granted Patent B2
US 12,217,434 · App. 18/208,331 · Granted Feb 4, 2025

Automated right ventricle medical imaging and computation of clinical parameters

Inventor: Marina Yaacobi (Beer-Sheva, IL)
Assignee: KONINKLIJKE PHILIPS N.V.
G06T7/11G06T3/4007G06T2207/10132G06T2207/20081G06T2207/20084G06T2207/20104G06T2207/20112G06T2207/30048
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,217,434
App. No.
18/208,331
Granted
Feb 4, 2025
Kind
B2
Abstract

There is provided a method of processing 2D ultrasound images for computing clinical parameter(s) of a right ventricle (RV), comprising: selecting one 2D ultrasound image of 2D ultrasound images depicting the RV, interpolating an inner contour of an endocardial border of the RV for the selected 2D image, tracking the interpolated inner contour obtained for the one 2D ultrasound image over the 2D images over cardiac cycle(s), computing a RV area of the RV for each respective 2D image according to the tracked interpolated inner contour, identifying a first 2D image depicting an end-diastole (ED) state according to a maximal value of the RV area for the 2D images, and a second 2D US image depicting an end-systole (ES) state according to minimal value of the RV area for the 2D images, and computing clinical parameter(s) of the RV according to the identified first and second 2D images.

Claims (62)

1. A computer-implemented method of automatically processing two-dimensional (2D) ultrasound images for computing of at least one clinical parameter of a right ventricle (RV), comprising:

selecting one 2D ultrasound image of a plurality of 2D ultrasound images depicting at least a RV of a subject, sequentially captured over at least one cardiac cycle of the subject;

classifying the RV depicted in the selected one 2D ultrasound image into a predefined shape selected from a plurality of predefined shapes;

interpolating an inner contour of an endocardial border of the RV for the selected one 2D ultrasound image according to the classified predefined shape of the selected one 2D ultrasound image;

tracking the interpolated inner contour obtained for the one 2D ultrasound image over the plurality of 2D images over the at least one cardiac cycle;

computing, a RV area of the RV for each respective 2D ultrasound image of the plurality of 2D ultrasound images, according to the tracked interpolated inner contour;

identifying a first 2D ultrasound image depicting an end-diastole (ED) state according to a maximal value of the RV area for the plurality of 2D images, and a second 2D US image depicting an end-systole(ES) state according to minimal value of the RV area for the plurality of 2D images; and

computing at least one clinical parameter of the RV according to the identified first 2D ultrasound image depicting the ED state and the second 2D US image depicting the ES state.

2. The method of claim 1 , further comprising, wherein the at least one clinical parameter is computed according to the first 2D ultrasound image depicting the ED state and the second 2D US image depicting the ES state with interpolated inner contour.

3. The method of claim 2 , further comprising:

inputting the selected 2D image into a trained machine learning (ML) model trained on a training dataset including for each of a plurality of sample subjects, at least one 2D ultrasound image depicting a RV and a corresponding plurality of segmentation masks including: a RV cavity segmentation mask, a lateral tricuspid valve annulus segmentation mask, and a septal tricuspid valve annulus segmentation mask; and

obtaining, for the selected 2D image, as an outcome of the ML model, a plurality of segmentation masks including: a RV cavity segmentation mask, a lateral tricuspid valve annulus segmentation mask, and a septal tricuspid valve annulus segmentation mask,

extracting, an inner contour of an endocardial border of the RV for the selected one 2D ultrasound image based on the plurality of segmentation masks,

extracting a plurality of valve indications for the selected one 2D ultrasound image,

wherein the plurality of valve indications include a plurality of positions of a tricuspid valve extracted from a center region of the lateral tricuspid valve annulus segmentation mask, and extracted from a center region of the septal tricuspid valve annulus segmentation mask.

4. The method of claim 3 , further comprising automatically selecting a region of interest from the selected 2D image, the ROI including the RV cavity, the lateral tricuspid valve annulus, and the septal tricuspid valve annulus, wherein inputting the selected 2D image comprises inputting the ROI into the trained ML model, wherein the at least one 2D ultrasound image of the training dataset comprises the ROI image created from the at least one 2D ultrasound image of the training dataset.

5. The method of claim 3 , further comprising, for the selected 2D image of the plurality of 2D ultrasound images:

extracting a respective border of the RV depicted in the selected 2D image using the respective cavity mask; and

generating a binary image using the respective extracted RV border,

wherein the extracting of the inner contour for the selected 2D image is based on the binary image.

6. The method of claim 5 , further comprising, for each 2D image of the plurality of 2D ultrasound images:

segmenting around the respective extracted RV border on sampled points using local contrast enhancement, wherein the binary image is of the respective extracted RV border with segmentation.

7. The method of claim 1 , further comprising:

identifying a tricuspid valve of the inner contour;

identifying an apex of the RV on the inner contour,

dividing the inner contour into a lateral side and a septal side with respect to the apex and the tricuspid valve,

wherein interpolating the inner contour is done each of the lateral side and the septal side according to the classified shape.

8. The method of claim 7 , wherein the classified shape is selected from one of “triangle” and “square”, and the lateral side is interpolated using a first interpolation method when the shape is triangle, and the lateral side is interpolated using a second interpolation method when the shape is square.

9. The method of claim 1 , wherein the 2D ultrasound images exclude a correlation with electrocardiogram (ECG) data of the subject, and the ED state and the ES state are identified from the plurality of 2D ultrasound images without the correlation with ECG data.

10. The method of claim 1 , further comprising:

creating a model of a trajectory of a tricuspid annulus location on a free wall between each pair of ED and ES denoting systole, and ES and next ED 2D images denoting diastole;

comparing a predicted valve location to actual location from tracking and correct high deviation using the model and repeat interpolation to compute a deviation;

when the deviation is above a threshold, the tracked valve location is corrected according to the model, and a border point interpolation process is repeated taking into account a corrected location of a left annulus point to obtain a corrected inner contour;

calculate area of a chamber denoting the area enclosed in the tracked contours for each 2D ultrasound image based on the corrected inner contour;

wherein the at least one clinical parameter comprises the calculated area.

11. The computer implemented method of claim 1 , wherein tracking further comprises tracking a plurality of locations of a tricuspid valve obtained for the one 2D ultrasound image, and computing the RV area further comprises computing the RV area according to the tracked plurality of locations of the tricuspid valve.

12. The computer implemented method of claim 1 , wherein the at least one clinical parameter comprises at least one of: (i) a fractional area change (FAC) for the RV, (ii) a free wall strain (FWS) for the RV denoting the maximum strain of the free wall, wherein the FWS is computed for each respective cardiac cycle for the RV, RV basal diameter (RVD1), RV mid diameter (RVD2) and RV length (RVD3).

13. The computer implemented method of claim 1 , wherein the at least one clinical parameter comprises at least one of: (i) a tricuspid annular plane systolic excursion (TAPSE), and further comprising computing the TAPSE by:

computing the distance of lateral free wall tricuspid valve annulus from an apex in ED to obtain a first value,

computing the distance of lateral free wall tricuspid valve annulus from the apex in ES to obtain a second value, and computing a difference between the first value and the second value, and

(ii) a peak systolic velocity of a tricuspid annulus (S′), and further comprising computing the S′ by: computing a displacement curve by computing a respective distance of lateral free wall tricuspid valve annulus from an apex each of the plurality of 2D US images of a respective cardiac cycle and subtracting from each respective distance the relative displacement denoting a change in distance, wherein the change in distance is computed by deriving a distance curve of the respective 2D image with respect to frame time to obtain the velocity of movement in a corresponding location, wherein S′ denotes the maximum of the velocity value per respective cardiac cycle.

14. The method of claim 1 , wherein each of the at least one clinical parameters is computed for each cardiac cycle depicted in the plurality of 2D ultrasound images.

15. The method of claim 1 , further comprising treating the subject for a medical condition selected from the group consisting of: impaired heart function, cardiac injury due to COVID-19, pulmonary embolism, pulmonary hypertension (PH), congenital heart disease (CHD), coronary artery disease, left-sided heart failure (HF), right-sided heart failure, and valvular heart disease, using a treatment effective for the medical condition, selected from the group consisting of: fibrinolysis, thrombolytic, anticoagulant, pacemaker, beta blockers, nitrates, oxygen, morphine, ACE inhibitors, insulin, ATII receptor blockers (ARBs), diuretics, cardiac defibrillation, cardiac surgery, trans catheterization, stent placement, valve replacement.

16. The method of claim 1 , further comprising iteratively computing the at least one clinical parameter over a plurality of time intervals, and plotting, on a display, a trend line of the at least one clinical parameters for the plurality of time intervals.

17. The method of claim 1 , further comprising selecting a computational processes from a plurality of computational process according to the predefined shape selected from the plurality of predefined shapes.

18. A system for automatically processing two dimensional (2D) ultrasound images for computing of at least one clinical parameter of a right ventricle (RV), comprising:

at least one hardware processor executing a code for:

selecting one 2D ultrasound image of a plurality of 2D ultrasound images depicting at least a RV of a subject, sequentially captured over at least one cardiac cycle of the subject;

classifying the RV depicted in the selected one 2D ultrasound image into a predefined shape selected from a plurality of predefined shapes;

interpolating an inner contour of an endocardial border of the RV for the selected one 2D ultrasound image according to the classified predefined shape of the selected one 2D ultrasound image;

tracking the interpolated inner contour obtained for the one 2D ultrasound image over the plurality of 2D images over the at least one cardiac cycle;

computing, a RV area of the RV for each respective 2D ultrasound image of the plurality of 2D ultrasound images, according to the tracked interpolated inner contour;

identifying a first 2D ultrasound image depicting an end-diastole (ED) state according to a maximal value of the RV area for the plurality of 2D images, and a second 2D US image depicting an end-systole(ES) state according to minimal value of the RV area for the plurality of 2D images; and

computing at least one clinical parameter of the RV according to the identified first 2D ultrasound image depicting the ED state and the second 2D US image depicting the ES state.

19. A non-transitory computer readable medium for automatically processing two dimensional (2D) ultrasound images for computing of at least one clinical parameter indicative of a right ventricle (RV) comprising program instructions which, when executed by a processor, cause the processor to perform:

selecting one 2D ultrasound image of a plurality of 2D ultrasound images depicting at least a RV of a subject, sequentially captured over at least one cardiac cycle of the subject;

classifying the RV depicted in the selected one 2D ultrasound image into a predefined shape selected from a plurality of predefined shapes;

interpolating an inner contour of an endocardial border of the RV for the selected one 2D ultrasound image according to the classified predefined shape of the selected one 2D ultrasound image;

tracking the interpolated inner contour obtained for the one 2D ultrasound image over the plurality of 2D images over at least one cardiac cycle;

computing, a RV area of the RV for each respective 2D ultrasound image of the plurality of 2D ultrasound images, according to the tracked interpolated inner contour;

identifying a first 2D ultrasound image depicting an end-diastole (ED) state according to a maximal value of the RV area for the plurality of 2D images, and a second 2D US image depicting an end-systole(ES) state according to minimal value of the RV area for the plurality of 2D images; and

computing at least one clinical parameter of the RV according to the identified first 2D ultrasound image depicting the ED state and the second 2D US image depicting the ES state.

Assignments (2)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Mar 6, 2024
From: DIA IMAGING ANALYSIS, LTD.
To: KONINKLIJKE PHILIPS N.V.
Reel/Frame 066661/0819 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jun 19, 2023
From: YAACOBI, MARINA
To: DIA IMAGING ANALYSIS
Reel/Frame 063982/0941 →
Continuity (2)
Continuation 17145345 · Jan 10, 2021
Related Publication 20230326034A1 · Oct 12, 2023
References Cited (24)
US 20110021915A1 · Feng et al. · 2011 [cited by applicant]
US 20110262018A1 · Kumar et al. · 2011 [cited by applicant]
US 20150038846A1 · Abe et al. · 2015 [cited by applicant]
US 20150310612A1 · Beymer · 2015 [cited by examiner]
US 20180279997A1 · Abe · 2018 [cited by applicant]
US 20180357767A1 · Arakita · 2018 [cited by examiner]
US 20220012875A1 · Arnaout · 2022 [cited by applicant]
US 20220222825A1 · Yaacobi · 2022 [cited by applicant]
Notice of Allowance Dated Mar. 22, 2023 from the US Patent and Trademark Office Re. U.S. Appl. No. 17/145,345. (7 pages). [cited by applicant]
Office Action Dated Dec. 2, 2022 from the US Patent and Trademark Office Re. U.S. Appl. No. 17/145,345. (33 pages). [cited by applicant]
Addetia et al. “3-Dimensional Echocardiographic Analysis of the Tricuspid Annulus Provides New Insights Into Tricuspid Valve Geometry and Dynamics”, JACC Journals, 12(3): 401-412, Mar. 2019. [cited by applicant]
Beecy et al. “Development of Novel Machine Learning Model for Right Ventricular Quantification on Echocardiography—A Multimodality Validation Study”, Echocardiography, 37(5): 688-697, Published Online May 12, 2020. [cited by applicant]
Focardi et al. “Traditional and Innovative Echocardiographic Parameters for the Analysis of Right Ventricular Performance in Comparison With Cardiac Magnetic Resonance”, European Heart Journal—Cardiovascular Imaging, 16… [cited by applicant]
Han “MR-Based Synthetic CT Generation Using A Deep Convolutional Neural Network Method”, Medical Physics, 44(4): 1408-1419, Published Online Feb. 13, 2017. [cited by applicant]
Helbing et al. “On-Line Automated Border Detection for Echocardiographic Quantification of Right Ventricular Size and Function in Children”, Pediatric Cardiology, 18(4): 261-269, 1997. [cited by applicant]
Horton et al. “Assessment of the Right Ventricle by Echocardiography: A Primer for Cardiac Sonographers”, Journal of the American Society of Echocardiography, 22(7): 776-792, Jul. 2009. [cited by applicant]
Lang et al. “Recommendations for Cardiac Chahmber Quantification by Echocardiography in Adults: An Update From the American Society of Echocardiography and the European Association of Cardiovascular Imaging”, Journal of… [cited by applicant]
Li et al. “Prognostic Value of Right Ventricular Longitudinal Strain in Patients With COVID-19”, Journal of the American College of Cardiology, JACC: Cardovascular Imaging, 13(11): 2287-2299, Published Online Apr. 28, 2… [cited by applicant]
Motoji et al. “Efficacy of Right Ventricular Free-Wall Longitudinal Speckle-Tracking Strain for Predicting Long-Term Outcome in Patients With Pulmonary Hypertension”, Circulation Journal, 7(3): 756-763, 2013. [cited by applicant]
Qin et al. “Automatic Segmentation of Right Ventricular Ultrasound Images Using Sparse Matrix Transform and A Level Set”, Physics in Medicine and Biology, 58(21): 7609-7624, Published Online Oct. 10, 2013. [cited by applicant]
Quinones et al. “ACC/AHA Clinical Competence Statement on Echocardiography. A Report of the American College of Cardiology / American Heart Association / American College of Physicians—american Society of Internal Medic… [cited by applicant]
Ronneberger et al. “U-Net: Convolutional Networks for Biomedical Image Segmentation”, International Conference on Medical Image Computing and Computer-Assisted Intervention, MICCAI 2015, Medical Image Computing and Comp… [cited by applicant]
Vignon et al. “Quantitative Assessment of Regional Right Ventricular Function with Color Kinesis”, American Journal of Respiratory and Critical Care Medicine, 159(6): 1949-1959, 1999. [cited by applicant]
Yaacobi et al. “Simultaneous Left Atrium Volume Tracking From Echocardiographic Movies”, IEEE Convention of the Electrical and Electronic Engineers in Israel, 403-407, 2008. [cited by applicant]