Imaging method for the representation of the results of intravascular imaging and CFD results and medical system for execution of the method

An imaging method for representing results of intravascular imaging (IVB) and CFD results of plaques of a vascular vessel tree of a patient and an imaging medical system for execution of the method are proposed. A 3-D image data set of the area under examination is acquired with x-ray beams for the generation of 3-D x-ray images. A 3-D reconstruction image of the area under examination is generated from the 3-D x-ray images. IVB images are generated by intravascular imaging. Larger calcified plaques in the 3-D x-ray images and IVB images are identified for determining orientation points or landmarks. The spatial position of the plaques are correlated and/or registered based on the orientation points or landmarks. The IVB image data and 3-D image data are fused and reproduced.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims priority of German application No. 10 2010 019 421.2 filed May 5, 2010, which is incorporated by reference herein in its entirety.

FIELD OF THE INVENTION

The invention relates to an imaging method for the representation of results of intravascular imaging (IVB) and CFD results of plaques of a vascular vessel tree of a patient and an imaging medical system for executing the method.

BACKGROUND OF THE INVENTION

In the visualization of the vascular vessel tree, the blood flow, the plaque formation, the plaque structure and further CFD parameters the following standard methods are used:

1. Imaging with 3-D Image Modalities

The 3-D image modalities serve to enable reconstruction of the vascular vessel tree. To this end a contrast medium is injected, in order to strengthen the blood vessels to be observed. The segmentation of the vessels from a 3-D data set represents the vascular lumen, which is used as the input for CFD methods.

The intravascular imaging, for example IVUS (Intravascular Ultrasound) or OCT (Optical Coherence Tomography), enables a mapping and analysis of the vessel wall, for example for the representation of plaque. The virtual histology function of IVUS provides automatic means for the detection and classification of plaque, for example the fat content, the fibrousness and the calcification.

The simulation of the blood flow by means of the CFD method (Computational Fluid Dynamics) delivers a three-dimensional distribution of the flow parameters, such as for example WSS (Wall Shear Stress), along the surface of the vessel lumen.

3-D DSA rotation angiography is a standard method used for the assessment of vascular anatomy before and during interventions. In digital subtraction angiography (DSA), after the creation of mask images, images without contrast medium, and filling images, images with contrast medium, these are subtracted from each other, so that only the temporal changes arising as a result of the contrast medium are obtained, which represent the vessels.

In neuroradiology especially, three-dimensional digital subtraction angiography (3-D DSA) is a routine tool for the planning and execution of minimally invasive procedures. Modem neurological operating theaters integrate this 3-D capability too with a rotating C-arm, in order to enable pre- and intra-procedural 3-D imaging of the cerebral blood vessels.

A C-arm x-ray machine for digital subtraction angiography represented by way of example inFIG. 1has for example a C-arm2rotatably mounted on a stand in the form of a six-axis industrial or articulated-arm robot1, on the ends of which are attached an x-ray emission source, for example a x-ray emitter3with x-ray tube and collimator, and an x-ray image detector4as an image acquisition unit.

By means of the articulated-arm robot1known, for example, from U.S. Pat. No. 7,500,784 B2, which preferably has six axes of rotation and thus six degrees of freedom, the C-arm2can be spatially adjusted at will, for example by being rotated about a center of rotation between the x-ray emitter3and the x-ray detector4. The inventive x-ray system1to4is in particular rotatable about centers of rotation and axes of rotation in the C-arm plane of the x-ray image detector4, preferably about the center point of the x-ray image detector4and about the axes of rotation bisecting the center point of the x-ray image detector4.

The known articulated-arm robot1has a base frame, which for example is fixedly mounted on a base. Thereupon is fixed a carousel which can be rotated about a first axis of rotation. Arranged pivotably about a second axis of rotation on the carousel is a robot pinion, to which is attached a robot arm, which can be pivoted about a third axis of rotation. On the end of the robot aim is arranged a robot hand, which is rotatable about a fourth axis of rotation. The robot hand has a fixing element for the C-arm2, which can be pivoted about a fifth axis of rotation and rotated about a sixth axis of rotation running perpendicular thereto.

The realization of the x-ray diagnostic device is not reliant on the industrial robot. Conventional C-arm devices can also be employed here.

The x-ray image detector4can be a rectangular or square, flat semiconductor detector, which is preferably created from amorphous silicon (a-Si). Integrating and possibly counting CMOS detectors can, however also be employed.

A patient6to be examined as the object under examination is located in the beam path of the x-ray emitter3on a patient couch5, for scanning a heart, for example. Connected to the x-ray diagnostic apparatus is a system control unit7with an image system8, which receives and processes the image signals from the x-ray image detector4(operating elements are, for example, not shown). The x-ray images can then be examined on a monitor9.

Other systems, for example for neuroradiology, use two C-arms. These are so-called biplane systems, as described in greater detail with reference toFIG. 2.

These essentially have two so-called planes, wherein the first plane10can comprise the x-ray diagnostic device shown inFIG. 1with C-arm2, x-ray emitter3and x-ray image detector4. Via a ceiling bracket11a ceiling-hung C-arm2′ with an x-ray emitter3′ and an x-ray image detector4′ of a second plane12can be provided. A monitor traffic light13with a first display14for the first plane10and a second display15for the second plane12can likewise be arranged on the ceiling. A high-voltage generator16is provided alongside the system control unit7.

In order to render the plaque more readily visible, a separate IVUS (Intravascular Ultrasound) catheter can be introduced into the vascular vessel tree of a patient. Such an IVUS system is for example described in DE 198 27 460 A1, from which a method for intravascular ultrasound mapping is known, in which an ultrasound signal transmitter and detector is introduced into a body lumen, and within which this can be moved. The ultrasound signal transmitter and receiver transmits ultrasound signals and captures reflected ultrasound signals, which contain information about the body lumen. A processor coupled to the ultrasound signal transmitter and receiver derives a first image series and a second image series from the captured ultrasound signals, and compares the second image series with the first image series. The processor can also be programmed to monitor the first and second image for cardiovascular periodicity, image quality, temporal change and vessel movement. It can also assign the first image series and the second image series to each other.

For intravascular imaging (IVB), a generally known OCT-catheter (Optical Coherence Tomography) can however also be introduced into the vessel.

CFD method

In DE 10 2008 014 792 B3 a method for the simulation of a blood flow in a vessel section is described, wherein an image acquisition of a vessel area encompassing the vessel section is obtained, a 3-D vessel section model is determined from the image acquisition, a number of blood flow parameters are read in, taking account of the, or of each blood flow parameter, the blood flow in the vessel section model is simulated and a number of hemodynamic parameters are output. It is here provided for that the image acquisition with an implant introduced into the vessel section is obtained in such a way that image data of the implant is included, and that the 3-D vessel section model is determined taking account of the image data of the implant employed. Further, a corresponding device for the simulation of a blood flow in a vessel section is specified.

In today's medical world, the dynamic behavior of the blood flow in an aneurysm is frequently regarded as an important factor for the pathogenesis of the aneurysm, that is for its occurrence and development.

As is known from the article “Image-Based Computational Simulation of Flow Dynamics in a Giant Intracranial Aneurysm” by D. A. Steinman, J. S. Milner, C. J. Morley, S. P. Lowie and D. W. Holdsworth from the American Journal of Neuroradiology (2003), Number 24, pp 559 through 566, a number of so-called hemodynamic parameters are connected with a growth and a rupture of the aneurysm. A hemodynamic parameter is in particular understood to be a parameter relating to hemodynamics, that is the flow dynamics of the blood. As hemodynamic parameters the article cited includes a pressure, a shear stress affecting the vessel wall and a flow rate.

In order to draw conclusions about hemodynamic parameters of this kind, for example, the blood flow in a vessel section, which for example includes the aneurysm, is for example simulated.

To this end, in the aforementioned article “Image-Based Computational Simulation of Flow Dynamic in a Giant Intracranial Aneurysm” a 3-D vessel section model is determined from a 3-D image acquisition, which was obtained by means of a rotation angiography. The blood flow in the 3-D vessel section model is simulated by means of the CFD method. The simulation is here performed assuming rigid vessel walls and a constant blood viscosity. CFD is a method of numeric flow simulation. The model equations used in the numeric flow mechanics are mostly based on a Nervier-Stokes equation, on a Euler or potential equation.

SUMMARY OF THE INVENTION

The object of the invention is to represent the results of intravascular imaging and CFD results of plaques of a vascular vessel tree of a patient, such as for example the plaque formation, the plaque structure and CFD parameters like the wall shear stress, and thus to improve the workflow in a simple manner.

The object is achieved by a method and by a device according to the features specified in the independent claims. Advantageous embodiments are specified in the dependent claims.

According to the invention the object is achieved for a method by the following steps:acquisition of a 3-D image data set of the area under examination with x-ray beams for the generation of 3-D x-ray images,generation of a 3-D reconstruction image of the area under examination from the 3-D x-ray images,generation of IVB images by means of intravascular imaging (IVB),identification of larger calcified plaques in the 3-D x-ray images and IVB images for the determining of orientation points or landmarks,correlation and/or registration of the spatial position of the plaques on the basis of the orientation points or landmarks andmerging and reproduction of the IVB image data and 3-D image data.

It is thereby possible in a simple manner to obtain a readily comprehensible representation of results of intravascular imaging and CFD results of plaques of a vascular vessel tree of a patient, such as for example the plaque formation, the plaque structure and CFD parameters such as wall shear stress.

According to the invention a minimizing of the distance of the voxels in the 3-D space can be performed for correlation and/or registration.

It has proved to be advantageous if the function for minimization of the distance of the voxels in the 3-D space takes the following form:

∑i⁢∑j⁢(riIVB-rjX)2
where riIVBor rjXrespectively determine the coordinates of the ithor jthvoxels respectively in the IVB or 3-D x-ray data set.

In an advantageous manner the identification of larger calcified plaques in the 3-D x-ray images and NB images can take place by means of identification of the voxels which contain calcium.

According to the invention the reproduction of the NB image data and 3-D image data can take place by means of color-coded 3-D rendering of CFD results derived from the IVB image data.

It has proved to be advantageous if the reproduction of the IVB image data and 3-D image data takes place based on corresponding cross-sectional information from the intravascular imaging (NB) and/or 3-D x-rays in a multiplanar format.

According to the invention the reproduction can take place in selected positions along the vessel, which for example displays abnormal WSS behavior.

According to the invention the NB images by means of intravascular imaging can take place with intravascular ultrasound (NUS) or optical coherence tomography (OCT).

According to the invention the object is achieved for a device by means ofa 3-D x-ray machine for the generation of a 3-D image data set,an IVB-device for the generation of at least one NB image,image storage means for buffering of the 3-D image data set and of the NB image,a 3-D image processing stage for the generation of a 3-D reconstruction image from the 3-D image data set and identification of blood vessels in the 3-D reconstruction image,an NB image processing stage for the identification of blood vessels in the NB image,a correlation or registration device for the correlation or registration of the 3-D reconstruction image and of the NB image with the aid of the data from the identification,a fusion device for precise vessel superimposition anda 3-D reproduction device for reproduction of the merged images.

In an advantageous manner the device for generation of an NB image can be an IVUS catheter and the device for generation of an IVB image an OCT catheter.

DETAILED DESCRIPTION OF THE INVENTION

FIG. 3shows an inventive medical system for the representation of results of intravascular imaging (IVB) and CFD results of plaques of a vascular vessel tree of a patient in conjunction with x-ray angiography with a 3-D x-ray machine20and an NB-device21. Connected to the 3-D x-ray machine20, for example for rotational angiography such as DynaCT, is a 3-D image memory22, in which is stored an acquired 3-D image data set, from which after acquisition of the complete 3-D image data set a 3-D reconstruction image is calculated by a 3-D image processing stage23, which is then likewise stored in the 3-D image memory22. Subsequently, likewise by means of the 3-D image processing stage23, the spatial position for example of the larger calcified plaques of the blood vessels in the 3-D reconstruction image, that is the voxels which contain calcium is identified and this data assigned to the 3-D image data set in the 3-D image memory22.

The NB image of the NB device21, for example an NUS or an OCT catheter, is buffered in an NB image memory24. By means of an NB image processing stage25the spatial position of the larger calcified plaques of the blood vessels in the NB image, and thus the voxels which contain calcium, are likewise identified, and this data assigned to the current NB image data set in the IVB image memory24.

By means of a correlation or registration device26, the 3-D reconstruction image is then registered from the 3-D image memory22and the NB image from the NB image memory24with the aid of the data from the identification of the larger calcified plaques.

For this the voxels which contain calcium have their distance in the 3-D space minimized. The function for minimization can especially take the following form:

∑i⁢∑j⁢(riIVB-rjX)2
where riIVBor rjXrespectively determine the coordinates of the ithor jthvoxels respectively in the IVB or 3-D x-ray data set.

Subsequently, the thus correlated or registered images are superimposed in a vessel-precise manner by means of a fusion device27and the 3-D reconstruction image overlaid with the NB image represented on a 3-D reproduction device28.

After such a spatial registration with the aid of the calcified orientation points or landmarks, the information can be reproduced in a merged manner in various ways.

One special possibility is the customary color-coded 3-D rendering of the CFD results (WWS). In selected positions along the vessel, which for example displays an abnormal WSS, the corresponding cross-section information from IVUS and/or 3-D x-rays can be reproduced in a multiplanar format, in order to show the composition of the vessel walls.

Alternatively, the 3-D data sets of the vessel tree captured by means of CT angiography with DynaCT can be automatically registered with 2-D angiography and fluoroscopy.

The known registration between 2-D x-rays and IVUS can then be used in order to register CTA (and CFD) with IVUS.

The inventive method procedure is now explained in greater detail on the basis ofFIG. 4. Firstly, in step S1), a complete 3-D image data set is acquired with the 3-D x-ray machine20. In step S2) a 3-D reconstruction image of the area under examination is calculated from the 3-D image data set by means of the 3-D image processing stage23, from which according to step S3) the spatial position of the larger calcified plaques of the blood vessels, that is to say the calcium-containing voxels, are identified, so that orientation points or landmarks are obtained.

In step S4), at least one IVB image is acquired with the NB-device21, which can be for example an IVUS or an OCT catheter. According to step S5) the spatial position of the larger calcified plaques of the blood vessels is identified from the IVB image for the determining of orientation points or landmarks, so that it is possible to define the voxels containing calcium.

According to step S6) a correlation and/or registration of the spatial position of the plaques takes place based on the orientation points or landmarks in the 3-D reconstruction images and NB images by means of the correlation or registration device26with the aid of a minimization of the distance of the defined voxels in the 3-D space.

Finally, the 3-D reconstruction images and NB images correlated or registered on the basis of the orientation points or landmarks are merged on a precise vessel basis in step S7) by means of the fusion device27and the 3-D reconstruction image data merged with the NB-image data represented on the 3-D reproduction device28according to step S8).

The basic idea of the present invention is to obtain a correlation between the plaque formation, the plaque structure and CFD parameters, for example the wall shear stress.

An approximation procedure, in order for example to correlate IVUS images and CFD results, is as follows:

For 3-D imaging with x-ray beams, the calcified plaque on the wall can be seen and detected, in the event where manual annotation may be necessary.

In principle the proposed method identifies the larger calcified plaques from NUS and 3-D x-ray data and correlates or registers their spatial position. For this the voxels which contain calcium are identified and their distance in the 3-D space minimized. The function for minimization can especially take the following form:

∑i⁢∑j⁢(riIVUS-rjX)2
where riIVUSor rjXrespectively determines the coordinates of the or voxels respectively in the IVUS or 3-D x-ray data set. After such a spatial registration with the aid of the calcified orientation points or landmarks, various paths can be realized, in order to represent the information in a merged manner.

One special possibility is the customary color-coded 3-D rendering of the CFD results (WWS). In selected positions along the vessel, which for example shows an abnormal WSS, the corresponding cross-section information from IVUS and/or 3-D x-rays can be represented in a multiplanar format, in order to show the composition of the vessel walls.

An alternative method of achieving this registration is the capturing of 3-D data sets of the vessel tree by means of rotation angiography, for example DynaCT. The CTA (CT angiography) can then be merged with DynaCT by means of the known 3-D/3-D fusion algorithms. DynaCT is automatically registered with 2-D angiography and fluoroscopy. The known registration between 2-D x-ray images and IVUS images can then be used to register CTA (and CFD) with IVUS.