IP Library Granted Patent US 9,600,891
Granted Patent B2
US 9,600,891 · App. 14/508,527 · Granted Mar 21, 2017

Web based fast query visualization of time-varying multi-variate vessel flow field by using uniform partition strategy

Inventors: Xiaoke Huang (Pittsburgh, PA); Jun Tao (Houghton, MI); Feng Qiu (Plainsboro, NJ); Daphne Yu (Yardley, PA)
Assignee: Siemens Healthcare GmbH
G06T7/0081A61B5/7425G06T19/20A61B5/026A61B5/055A61B6/504A61B8/06A61B2576/02G06F17/5009
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 9,600,891
App. No.
14/508,527
Filed
Oct 7, 2014
Granted
Mar 21, 2017
Kind
B2
Art Unit
2611
USPC
345/427
Abstract

A method for visualizing flow data from computation fluid dynamics (CFD) applications in 2-dimensions (2D) includes receiving a 3-dimensional (3D) image volume from a CFD simulation of fluids flowing through vessels in a patient that is a snapshot of a fluid flow in the vessels at a certain time, subdividing the 3D image volume into 3D data blocks, minimizing a sum over a matrix of energy interactions defined for each pair of data blocks in the 3D image volume, where the minimization preserves a local shape of the vessels, where minimizing the sum over the matrix of energy interactions is performed on a graphics processing unit (GPU), and using the minimized energy interaction matrix to display on a monitor a 2D sketch of the 3D image volume, where the 2D sketch is displayed in real-time with respect to the time scale of the CFD simulation.

Claims (106)

1. A method for visualizing flow data from computation fluid dynamics (CFD) applications in 2-dimensions (2D), comprising the steps of:

receiving a 3-dimensional (3D) image volume from a CFD simulation of fluids flowing through vessels in a patient that is a snapshot of a fluid flow in the vessels at a certain time;

subdividing the 3D image volume into 3D data blocks;

minimizing a sum over energy interactions defined for each pair of data blocks in the 3D image volume, wherein said minimization preserves a local shape of the vessels, wherein minimizing the sum over energy interactions is performed on a graphics processing unit (GPU), wherein minimizing the sum over the energy interactions comprises minimizing E=√{square root over (Σ ij E ij 2 )}

wherein

E

ij

=

{

W

n

×

t

ij

-

e

ij

,

e

ij

<

d

,

W

nn

×

t

ij

-

g

ij

,

g

ij

>

d

,

wherein w n and w nn are weights for neighboring and non-neighboring blocks, respectively, t ij =|v i −v j | is a target distance between blocks b i and b j in a 2D image wherein v i and v j are 2D positions of b i and b j , e ij is a Euclidean distance between b i and b j in the 3D volume, g ij is a geodesic distance between b i and b j in the vessels in the 3D image, and d is a threshold that separates neighboring blocks from non-neighboring blocks, wherein the minimization is performed by varying the target positions v i and v j ; and

using the minimized energy E to display on a monitor a 2D sketch of the 3D image volume, wherein the 2D sketch is displayed in real-time with respect to the time scale of the CFD simulation.

2. The method of claim 1 , further comprising loading the 3D data blocks into memory of a web application for display, wherein the web application is accessible over the World Wide Web.

3. The method of claim 2 , further comprising distinguishing non-empty 3D data blocks that contain vessels from empty 3D data blocks that do not contain vessels, and only loading the non-empty 3D data blocks.

4. The method of claim 3 , further comprising calculating two look-up-tables (LUTs) when loading the non-empty data blocks, wherein an energy interaction is defined for each pair of non-empty data blocks.

5. The method of claim 4 , wherein a first LUT is a mapping from a descriptor of the 3D image volume to an index that indicates a starting position in memory of the 3D image volume, and a second LUT indicates a memory position of each 3D data block of the image of the 3D image volume.

6. The method of claim 1 , further comprising receiving a selection of particles selected from the 2D sketch, and calculating and displaying statistical information about the selected particles.

7. The method of claim 6 , wherein the information includes one or more of velocity, vorticity, wall shear stress, inlet and pressure.

8. The method of claim 6 , wherein the selection is performed using a Cross filter.

9. The method of claim 6 , wherein the selection is performed using a Boolean filter.

10. The method of claim 1 , wherein the vessels are at least a part of a cardiovascular system of the patient.

11. The method of claim 1 , further comprising weighting each energy interaction E ij by a weight W n if a separation distance between block i and block j in 3D is less than the threshold distance d, and by a weight W nn if the separation distance between block i and block j in 3D is greater than or equal to the threshold distance d, wherein W n >W nn .

12. A non-transitory program storage device readable by a computer, tangibly embodying a program of instructions executed by the computer to perform the method steps for visualizing flow data from computation fluid dynamics (CFD) applications in 2-dimensions (2D), the method comprising the steps of:

receiving a 3-dimensional (3D) image volume from a CFD simulation of fluids flowing through vessels in a patient that is a snapshot of a fluid flow in the vessels at a certain time;

subdividing the 3D image volume into 3D data blocks;

minimizing a sum over energy interactions defined for each pair of data blocks in the 3D image volume, wherein said minimization preserves a local shape of the vessels, wherein minimizing the sum over energy interactions is performed on a graphics processing unit (GPU), wherein minimizing the sum over the energy interactions comprises minimizing E=√{square root over (Σ ij E ij 2 )}

wherein

E

ij

=

{

W

n

×

t

ij

-

e

ij

,

e

ij

<

d

,

W

nn

×

t

ij

-

g

ij

,

g

ij

>

d

,

wherein w n and w nn are weights for neighboring and non-neighboring blocks, respectively, t ij =|v i −v j | is a target distance between blocks b j and b j in a 2D image wherein v i and v j are 2D positions of b i and b j , e ij is a Euclidean distance between b i and b j in the 3D volume, g ij is a geodesic distance between b i and b j in the vessels in the 3D image, and d is a threshold that separates neighboring blocks from non-neighboring blocks, wherein the minimization is performed by varying the target positions v i and v j ; and

using the minimized energy E to display on a monitor a 2D sketch of the 3D image volume, wherein the 2D sketch is displayed in real-time with respect to the time scale of the CFD simulation.

13. The computer readable program storage device of claim 12 , the method further comprising loading the 3D data blocks into memory of a web application for display, wherein the web application is accessible over the World Wide Web.

14. The computer readable program storage device of claim 13 , the method further comprising distinguishing non-empty 3D data blocks that contain vessels from empty 3D data blocks that do not contain vessels, and only loading the non-empty 3D data blocks.

15. The computer readable program storage device of claim 14 , the method further comprising calculating two look-up-tables (LUTs) when loading the non-empty data blocks, wherein an energy interaction is defined for each pair of non-empty data blocks.

16. The computer readable program storage device of claim 15 , wherein a first LUT is a mapping from a descriptor of the 3D image volume to an index that indicates a starting position in memory of the 3D image volume, and a second LUT indicates a memory position of each 3D data block of the image of the 3D image volume.

17. The computer readable program storage device of claim 12 , the method further comprising receiving a selection of particles selected from the 2D sketch, and calculating and displaying statistical information about the selected particles.

18. The computer readable program storage device of claim 17 , wherein the information includes one or more of velocity, vorticity, wall shear stress, inlet and pressure.

19. The computer readable program storage device of claim 17 , wherein the selection is performed using a Cross filter.

20. The computer readable program storage device of claim 17 , wherein the selection is performed using a Boolean filter.

21. The computer readable program storage device of claim 12 , wherein the vessels are at least a part of a cardiovascular system of the patient.

22. The computer readable program storage device of claim 12 , the method further comprising weighting each energy interaction E ij by a weight W n if a separation distance between block i and block j in 3D is less than the threshold distance d, and by a weight W nn if the separation distance between block i and block j in 3D is greater than or equal to the threshold distance d, wherein W n >W nn .

Assignments (5)
CORRECTIVE ASSIGNMENT TO CORRECT THE ASSIGNEE PREVIOUSLY RECORDED AT REEL: 066088 FRAME: 0256. ASSIGNOR(S) HEREBY CONFIRMS THE ASSIGNMENT. Recorded Jan 17, 2024
From: SIEMENS HEALTHCARE GMBH
To: SIEMENS HEALTHINEERS AG
Reel/Frame 071178/0246 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Dec 20, 2023
From: SIEMENS HEALTHCARE GMBH
To: SIEMENS HEALTHINEERS AG
Reel/Frame 066088/0256 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jan 31, 2017
From: SIEMENS CORPORATION
To: SIEMENS AKTIENGESELLSCHAFT
Reel/Frame 041128/0591 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jan 31, 2017
From: SIEMENS AKTIENGESELLSCHAFT
To: SIEMENS HEALTHCARE GMBH
Reel/Frame 041128/0683 →
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jan 17, 2017
From: HUANG, XIAOKE; QIU, FENG; TAO, JUN; YU, DAPHNE
To: SIEMENS CORPORATION
Reel/Frame 040980/0925 →
Continuity (2)
Provisional Application 61888050 · Oct 8, 2013
Related Publication 20150097836A1 · Apr 9, 2015