IP Library Granted Patent US 12,629,121
Granted Patent B2
US 12,629,121 · App. 18/136,006 · Granted May 19, 2026

Method of determining vessel fluid flow velocity

Inventors: Catalina Tobón Gómez (Haarlem, NL); Hua Ma (Rotterdam, NL); Johannes Petrus Janssen (Leiden, NL); Gianni Pedrizzetti (Florence, IT); Johan Hendrikus Christiaan Reiber (Rotterdam, NL)
Assignee: QFR Solutions B.V.
A61B6/504A61B6/481A61B6/486A61B6/5217G06T7/0016G06T2207/30104
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,629,121
App. No.
18/136,006
Granted
May 19, 2026
Kind
B2
Abstract

Image data of contrast dye in a vessel in a body is acquired for determining flow rate. Based on images acquired under and angle relative to one another, a three-dimensional model of the vessel is constructed and length of a vessel section is determined. A series of at least two images, apart in time, under a first angle is assessed for determining progress of a front of the dye bolus in the vessel in time. In the images, the vessel may be segmented and brightness or a derivative thereof over at least one of time and distance may be assessed to determine the front. Progress distance is mapped to the three-dimensional model, for example by mapping segments from the image to the model, to obtain a more accurate and natural distance of progress over time. Flow rate is determined by natural progress distance over progress time.

Claims (34)

1 . A method of determining vessel fluid flow velocity of a fluid in a vessel segment of a body of a human or another mammal, the method comprising:

obtaining first two-dimensional image data of the vessel segment under a first angle and obtaining second two-dimensional image data of the vessel segment under a second angle, to determine a natural length model of the vessel representing a length of the length of the vessel,

wherein the first two dimensional image data are obtained from a first series of two-dimensional images acquired consecutively in time and the second two-dimensional image data are obtained from a second series of two-dimensional images acquired consecutively in time;

obtaining a two-dimensional model of the vessel segment;

dividing the two-dimensional model in vessel sections;

selecting a set of images from the first series of two-dimensional images or from the second series of two-dimensional images, based on the series depicting the vessel with a longest length,

wherein the set of images includes image data of the vessel sections, each image corresponding to a different moment in time;

based on the images of the set, assigning, for each image of the set, at least one intensity value to each vessel section;

identifying, in a first image related to a first moment in time, a first vessel section, based on an intensity criterion;

identifying, in a second image related to a second moment in time, a second vessel section, based on the intensity criterion, the second moment in time being later than the first moment in time and the second vessel section being distal to the first vessel section;

obtaining a propagation length by relating the first vessel section and the second vessel section to the natural length; and

determining a vessel fluid flow velocity based on the propagation length, and a difference between the first moment in time and the second moment in time.

2 . The method according to claim 1 , wherein determining the natural length model further comprises:

constructing a three-dimensional model of the vessel based on the first two-dimensional image and the second two-dimensional image; and

determining a length of the centreline of the three-dimensional model of the vessel.

3 . The method according to claim 1 , wherein the set of images is acquired under the first angle or the second angle.

4 . The method according to claim 2 , wherein dividing the two-dimensional model in vessel sections comprises:

defining volume sections along the centreline of the of the three-dimensional model of the vessel; and

mapping the volume sections to the series of images as vessel sections.

5 . The method according to claim 1 , wherein dividing the two-dimensional model in vessel sections comprises:

defining the vessel sections in the two-dimensional model; and

associating the vessel sections to points in the natural length model.

6 . The method according to claim 5 , wherein determining the natural length model comprises:

obtaining first two-dimensional image data of the vessel segment under a first angle;

obtaining second two-dimensional image data of the vessel segment under a second angle;

constructing a three-dimensional model of the vessel based on the first two-dimensional image and the second two-dimensional image;

determining a length of the centreline of the three-dimensional model of the vessel; and

associating the vessel sections to points in the natural length model comprises associating the vessel sections with volume sections along the centreline of the three-dimensional model of the vessel.

7 . The method according to claim 1 , wherein the intensity criterion comprises a threshold value and the criterion is held to be met if the intensity of image data associated with a segment is at least one of more or less than the threshold value.

8 . The method according to claim 1 , wherein the intensity criterion comprises a time derivative threshold value of the intensity value over time.

9 . The method according to claim 1 , further comprising normalising, for each vessel section, the intensity value over the images of the set.

10 . The method according to claim 2 , wherein the first two-dimensional image data and the second two-dimensional image data are obtained using x-ray.

11 . A non-transitory computer-readable medium comprising instructions for causing a processor to carry out the method of claim 1 .

12 . A data processing device comprising a processor adapted to perform the method of claim 1 .

Assignments (2)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded May 19, 2025
From: MEDIS ASSOCIATED B.V.
To: QFR SOLUTIONS B.V.
Reel/Frame 071154/0586 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Nov 8, 2023
From: TOBÓN GÓMEZ, CATALINA; MA, HUA; JANSSEN, JOHANNES PETRUS; PEDRIZZETTI, GIANNI; REIBER, JOHAN HENDRIKUS CHRISTIAAN
To: MEDIS ASSOCIATED B.V.
Reel/Frame 065491/0721 →
Priority Claims (1)
NL 2026715 · Oct 20, 2020 · national
Continuity (1)
Related Publication 20250032079A1 · Jan 30, 2025
References Cited (27)
US 10687774B2 · Lautenschläger · 2020 [cited by applicant]
US 20040265393A1 · Unger et al. · 2004 [cited by applicant]
US 20090016587A1 · Strobel · 2009 [cited by examiner]
US 20110150274A1 · Patwardhan · 2011 [cited by examiner]
US 20130064435A1 · Taerum · 2013 [cited by examiner]
US 20150297373A1 · Schmitt · 2015 [cited by examiner]
US 20170140532A1 · Dascal · 2017 [cited by examiner]
US 20180099125A1 · Richer et al. · 2018 [cited by applicant]
US 20190029624A1 · Kunio · 2019 [cited by applicant]
US 20190365336A1 · Wagner · 2019 [cited by examiner]
US 20200175679A1 · Reiber · 2020 [cited by examiner]
US 20200219252A1 · Tsuyuki · 2020 [cited by examiner]
US 20200222018A1 · van Walsum et al. · 2020 [cited by applicant]
US 20220343494A1 · Dhatt et al. · 2022 [cited by applicant]
US 20220378396A1 · Yoshida · 2022 [cited by examiner]
CN 113887463A · 2022 [cited by applicant]
EP 3660858A1 · 2020 [cited by applicant]
JP 2002534204A · 2002 [cited by applicant]
JP 2016159116A · 2016 [cited by applicant]
JP 202295470A · 2022 [cited by applicant]
WO 2019210553A1 · 2019 [cited by applicant]
WO 2022086326A1 · 2022 [cited by applicant]
Mendonca, Ana Maria, and Aurelio Campilho. “Segmentation of retinal blood vessels by combining the detection of centerlines and morphological reconstruction.” IEEE transactions on medical imaging 25.9 (2006): 1200-1213.… [cited by examiner]
XP019043569 (Hrvoje Bogunovic et al): “Blood Flow and Velocity Estimation Based on Vessel Transit Time by Combining 2D and 3D X-Ray Angiography”, (Jan. 1, 2006)). [cited by applicant]
XP060008077 (Hentschke Clemens M et al): “Estimatingb lood flow velocity in angiographic image data”, (Mar. 3, 2011), (bladzijden 1-9). [cited by applicant]
International Search Report and Written Opinion—PCT/NL2021/050635—mailing date Apr. 28, 2022. [cited by applicant]
Gibson et al, Circulation, 93, 5 (1996) 879-888. [cited by applicant]