US20010012328A1 - Method of imaging the blood flow as a function of time in an object to be examined - Google Patents

Method of imaging the blood flow as a function of time in an object to be examined Download PDF

Info

Publication number
US20010012328A1
US20010012328A1 US09/749,147 US74914700A US2001012328A1 US 20010012328 A1 US20010012328 A1 US 20010012328A1 US 74914700 A US74914700 A US 74914700A US 2001012328 A1 US2001012328 A1 US 2001012328A1
Authority
US
United States
Prior art keywords
ray
image data
images
data set
projection images
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.)
Granted
Application number
US09/749,147
Other versions
US6442235B2 (en
Inventor
Reiner Koppe
Erhard Klotz
Michael Kuhn
John Op De Beek
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Assigned to U.S. PHILIPS CORP. reassignment U.S. PHILIPS CORP. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: KUHN, MICHAEL HARALD, OP DE BEEK, JOHN, KLOTZ, ERHARD PAUL ARTUR, KOPPE, REINER HEINRICH
Publication of US20010012328A1 publication Critical patent/US20010012328A1/en
Assigned to KONINKLIJKE PHILIPS ELECTRONICS N.V. reassignment KONINKLIJKE PHILIPS ELECTRONICS N.V. ASSIGNMENT OF ASSIGNORS INTEREST (SEE DOCUMENT FOR DETAILS). Assignors: U.S. PHILIPS CORPORATION
Application granted granted Critical
Publication of US6442235B2 publication Critical patent/US6442235B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/48Diagnostic techniques
    • A61B6/481Diagnostic techniques involving the use of contrast agents
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/40Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis
    • A61B6/4007Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units
    • A61B6/4014Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment with arrangements for generating radiation specially adapted for radiation diagnosis characterised by using a plurality of source units arranged in multiple source-detector units
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/50Clinical applications
    • A61B6/504Clinical applications involving diagnosis of blood vessels, e.g. by angiography

Definitions

  • the invention relates to a method of imaging the blood flow as a function of time in an object to be examined, as well as to an X-ray device for carrying out this method.
  • EP 860 696 A2 discloses a method which enables the reproduction of distributed structures, for example a vascular system filled with a contrast medium, in a synthetic projection image which reproduces the parts of the vascular system that are present in a selectable sub-volume more clearly than the X-ray images that are acquired from different perspectives and wherefrom this projection image is derived. This method thus yields a time-averaged “frozen” image of the vascular system in which the blood flow as a function of time is not visible.
  • Three-dimensional rotation angiography produces a series of X-ray projection images of the object to be examined from different projection directions while a contrast medium is injected into the blood vessels of the object to be examined.
  • a three-dimensional image is derived from such X-ray projection images; this three-dimensional image reproduces the vascular system in space in a time-averaged manner.
  • two-dimensional angiography it is also known to form a two-dimensional “frozen” image of the vascular system.
  • Such a method should also take into account the fact that a contrast medium cannot be repeatedly injected into the same patient within a short period of time, so that the method should do, if at all possible, with a single contrast medium injection. Moreover, it should be possible to acquire the images in an as short as possible period of time, at an as small as possible expenditure and with an as high as possible resolution.
  • the invention is based on the recognition of the fact that an image data set, which may be a two-dimensional or a three-dimensional image data set and contains information concerning the course of the blood vessels in the object to be examined, can be encoded in time in such a manner that it also contains information concerning the blood flow as a function of time.
  • Such encoding in time is performed according to the invention in that the image data set is compared with a series of X-ray projection images; these X-ray projection images are formed successively in time and contain the information concerning the distribution of an injected contrast medium in the blood vessels at each time a different instant.
  • each X-ray projection image is individually compared with the image data set, that is, each image value of the image data set is compared with the image values of the individual X-ray projection images, it is quasi checked which parts of the vascular system contained in the image data set are filled with the contrast medium at the individual instants associated with the respective X-ray projection images.
  • the image data set thus encoded in time can be converted into one or more images which show the blood flow as a function of time.
  • the version disclosed in claim 2 is particularly suitable for two-dimensional rotation angiography. It already offers a two-dimensional X-ray image data set, for example a two-dimensional X-ray projection image which completely contains the vascular system filled with contrast medium and may be a previously formed or an instantaneous X-ray image data set.
  • the actual X-ray projection images derived from a fixed X-ray position during a contrast medium administration are subtracted from one another in conformity with this version, so that each difference image contains the information as regards the path followed by the contrast medium in the vascular system between the two instants at which the two subtracted X-ray projection images were acquired. Such difference images are then used for the time encoding of the X-ray image data set.
  • a further version of this method is disclosed in claim 3 .
  • This version represents a simple possibility for comparing the image data set with the difference images.
  • an X-ray image data sub-set is segmented from the X-ray image data set; this sub-set contains only the information concerning the course of the blood vessels.
  • an associated pixel sub-set in such a manner that the pixels of the individual difference images are compared with the pixels of the X-ray image data sub-set and that with each pixel sub-set there are associated those pixels of the X-ray image data sub-set for which the associated difference image includes corresponding pixels.
  • Each pixel sub-set thus contains the information concerning the distribution of the contrast medium at a given instant and one or more projection images which show the blood flow as a function of time can be derived from the pixel sub-sets.
  • the X-ray projection images are acquired from different positions.
  • the image data set then constitutes a three-dimensional X-ray image data set derived from such X-ray projection images.
  • a single series of X-ray projection images is acquired from different directions, said images also containing the time information.
  • pseudo-projection images are calculated from the X-ray image data set, utilizing the known imaging geometry of the X-ray device, so as to be compared with the actual X-ray projection images.
  • claim 5 discloses a particularly attractive further version which is also suitable for three-dimensional rotation angiography.
  • Two series of X-ray projection images are acquired therein, the first series being acquired from a fixed X-ray position whereas the second series is acquired from different X-ray positions.
  • This operation can be performed either successively in time, be it that two contrast medium injections are then required, or simultaneously by means of an X-ray device which includes two imaging units.
  • From the X-ray projection images of the first series there are derived difference images which, as described above, contain the time information whereas a three-dimensional X-ray image data set is acquired from the X-ray projection images of the first series by means of a known reconstruction algorithm.
  • This X-ray image data set is encoded in time by means of the difference images.
  • the claims 6 and 7 disclose preferred possibilities for the time encoding of the image data set and the comparison of the image data set with the segmented X-ray projection images.
  • the steps of the method as disclosed in claim 7 correspond to the method which is known from EP 880 109 A2 whose disclosure is explicitly referred to herein and is considered to be incorporated herein.
  • the latter publication describes a method of determining the spatial transformation between a three-dimensional object reproduced by a data set and the object itself.
  • a pseudo-projection image is then calculated for a part of the volume reproduced by the data set; this pseudo-projection image is compared with an X-ray projection image of the object itself.
  • the parameters on which the calculation of the pseudo-projection image is based are then varied until an optimum match is obtained.
  • This method can be advantageously used in an appropriate manner in conjunction with the present invention.
  • the claims 8 and 9 disclose advantageous further versions concerning the display of the blood flow as a function of time.
  • the time information can be converted into a color code so that the complete two-dimensional or three-dimensional data set can be reproduced as an image with the corresponding color code.
  • the encoded pixel sub-sets or voxel sub-sets can also be reproduced successively in time, for example as an endless loop, thus creating the impression of blood flowing through the blood vessels. It is also feasible to enable reproductions from arbitrary angles of observation and rotation of the image is also possible.
  • An X-ray device which is suitable for carrying out the method according to the invention is disclosed in claim 10 ; as is indicated in the embodiment of claim 11 , it may also include a second imaging unit.
  • FIG. 1 shows an X-ray device for carrying out the method according to the invention
  • FIG. 2 shows a flow chart illustrating a first version of the method according to the invention
  • FIG. 3 illustrates diagrammatically the encoding in time
  • FIG. 4 shows a flow chart of a second version of the method according to the invention.
  • FIG. 5 shows a flow chart of a third version of the method according to the invention.
  • FIG. 1 shows an X-ray device 1 which serves to form two-dimensional X-ray images of an object 3 to be examined, for example, a patient who is arranged on a table 4 .
  • the X-ray device 1 includes an X-ray source 12 and an X-ray detector 13 which are aligned relative to one another and mounted on a C-arm 10 which itself is journalized in a stand 11 which is only partly shown.
  • the C-arm 10 can be pivoted about a perpendicular axis on the one hand and be rotated around its center in the direction of the double arrow 20 on the other hand, for example through 180° C., by means of a motor drive (not shown).
  • a plurality of X-ray images can be acquired so as to image the object 3 to be examined from different, reproducible perspectives or X-ray positions r 1 , r 2 , r 3 of the imaging unit 12 , 13 .
  • a second X-ray source 12 ′ and a second X-ray detector 13 ′ which are mounted on a supporting device 11 ′ and are capable of forming X-ray images of the object 3 to be examined from a fixed X-ray position r 4 .
  • Each of the X-ray detectors 13 , 13 ′ may be formed by a respective X-ray image intensifier whereto there is connected a television chain whose output signals are digitized by an analog-to-digital converter 14 so as to be stored in a memory 15 .
  • r m (of which only the positions r 1 , r 2 , r 3 are explicitly shown in the drawing), can be processed by an image processing unit 16 so as to be reproduced on a monitor 18 either individually or as a series of images.
  • the X-ray projection images E 0 , E 1 , E 2 . . . , En, acquired by the imaging unit 12 ′, 13 ′ from the fixed X-ray position r′ at discrete instants t during a contrast medium bolus T, can also be processed by the image processing unit 16 so as to be displayed on the monitor 18 .
  • the individual components of the X-ray device are controlled by means of a control unit 17 .
  • an arithmetic unit 19 which receives the X-ray projection images D i and E j in order to derive therefrom images showing the variation in time of the blood flow in the relevant vessels of the object 3 to be examined. These images can be displayed on the monitor 18 again.
  • the invention will be illustrated hereinafter on the basis of the flow chart of a first version of the method according to the invention which is shown in FIG. 2.
  • the X-ray projection images Di constitute a three-dimensional X-ray image data set K.
  • the imaging unit 12 ′, 13 ′ acquires a second series of X-ray projection images E j from a fixed perspective.
  • correction is made for imaging errors which are due, for example, to the imperfection of the imaging device or to the mechanical deformation of the C-arm.
  • a three-dimensional reconstruction image R is formed from the X-ray projection images Di by means of a known reconstruction algorithm.
  • a reconstruction sub-image R′ is formed from said reconstruction image R, which reconstruction sub-image contains, mainly or exclusively, information concerning the course of the blood vessels in the region examined.
  • This reconstruction sub-image R′ contains q voxels V k which are characterized by their co-ordinates in space.
  • respective difference images F j are determined from each time two successively acquired X-ray projection images E j-1 and E j ; such difference images contain only the information concerning the path traveled by the contrast medium in the period of time elapsing between the formation of the two X-ray projection images E j-1 and E j .
  • each of the n difference images F j is compared with the reconstruction sub-image R′, so that the latter is encoded in time.
  • each of the q voxels V k of the reconstruction sub-image R′ is projected onto each difference image F j , that is, a pseudo-projection image is calculated from the reconstruction sub-image R′ so as to be compared with the individual difference images F j .
  • this pseudo-projection image it is necessary to know the geometry of the imaging unit 12 ′, 13 ′ in order to ensure that the co-ordinate system of the pseudo-projection image corresponds to the co-ordinate system of the X-ray projection images E j and the difference images F j .
  • the voxels of the reconstruction sub-image R′ are marked each time.
  • the marked voxels are combined so as to form a voxel sub-set L j in the step 107 .
  • n voxel sub-sets L j wherefrom one or more images B can be formed in the step 108 , said images representing the blood flow as a function of time.
  • the voxels of each voxel sub-set L j can be reproduced in a different color in an overall image.
  • the method is terminated in the step 109 .
  • the step 106 will be illustrated again with reference to FIG. 3.
  • the pixels V k of this reconstruction sub-image R′ are projected onto the difference images F j .
  • the voxel V 1 along the projection ray 21 is projected onto the pixel P 1 .
  • a voxel or a pixel represents a part of a vascular structure, as can be detected, for example, on the basis of the image value which lies above or below a given threshold value, it can be recognized that in the difference image F 1 a corresponding pixel P 1 exists for the voxel V 1 .
  • a reconstruction sub-image R′ can also be obtained in a manner other than that described with reference to FIG. 2, for example, by subtraction angiography.
  • two series of X-ray projection images are formed from different X-ray positions by means of the imaging unit 12 , 13 , that is, once with and once without administration of contrast medium.
  • the X-ray projection images of the two series that are formed from each time the same X-ray position are then subtracted from one another and the resultant projection images are reconstructed so as to form the desired reconstruction sub-image R′.
  • FIG. 4 shows the flow chart of a further version of the method according to the invention.
  • This flow chart deviates from that of the version shown in FIG. 2 first of all in that only a single series of X-ray projection images D i is acquired; no X-ray projection images E j are acquired from a fixed X-ray position, so that the second imaging unit 12 ′, 13 ′ (see FIG. 1) can also be dispensed with.
  • the steps 200 to 203 otherwise correspond without modification to the steps 100 to 104 with omission of the step 102 .
  • the X-ray projection images Di are segmented, that is, the pixels representing the vascular structures in the individual images D i are selected and reproduced exclusively in the segmented X-ray projection images D i ′, while other image elements such as, for example the background or other organs, are eliminated.
  • This operation can be performed, for example, by means of the image values which lie above or below a given limit value for the vessels filled with contrast medium whereas this is usually not the case for other image elements.
  • the successive segmented X-ray projection images Di are subsequently individually compared with the reconstruction sub-image R′ or the pixels P 1 of the segmented X-ray projection images D i ′ are compared with the voxels V k of the reconstruction sub-image R′.
  • This operation takes place in the same way as described with reference to FIG. 3 for the difference images F j .
  • the voxel sub-sets L j and the images B for reproducing the blood flow are formed in the steps 206 and 207 in the same way as already described above. The method is completed again in the step 208 .
  • the present version offers the advantages that only one imaging unit ( 12 , 13 ) is required and that only a single series of X-ray projection images Di need be acquired.
  • the segmentation of the images D i in the step 204 may be difficult in practice and could give rise to inaccuracies.
  • the information contents of the segmented X-ray projection images D i ′ differ from that of the difference images F j which contain the variation of the blood flow as a function of time in a time interval whereas the segmented X-ray projection images D i ′ show the entire region of the vascular system filled with the contrast medium.
  • a two-dimensional X-ray image data set H is acquired in the step 301 .
  • This may be a single X-ray projection image which reproduces the complete vascular system filled with a contrast medium, or an image composed of a plurality of individual X-ray projection images.
  • the data set H may also have been acquired already at an earlier instant.
  • a series of X-ray projection images E j is acquired (after contrast injection) from a fixed X-ray position by means of the imaging unit ( 12 , 13 ).
  • the two-dimensional image data set H is segmented in the step 303 , so that the segmented image data set H′ contains only the vessel structures.
  • step 304 respective difference images F j , having the already described information contents are determined from each time two temporally successively acquired X-ray projection images E j-1 , E j .
  • These difference images F j are then individually compared successively with the segmented image data set, that is, the pixels P 1 of each individual difference image F j are compared with the pixels P k of the segmented two-dimensional image data set H′.
  • the X-ray projection images E j and the two-dimensional image data set have been acquired by means of the same imaging unit and from the same X-ray position.
  • This comparison in the step 305 produces in an associated pixel sub-set M j for each difference image F j in the step 306 ; this pixel sub-set is comparable to the already described voxel sub-sets L j and contains the information as to which pixels of the segmented image data set H′ represent a part of the vascular system which has been filled with contrast medium in a given time interval or wherefrom the contrast medium has been drained during this time interval.
  • step 307 the entire segmented image H′ will have been encoded in time and one or more images can be displayed in order to visualize the blood flow in this image as a function of time (step 307 ), thus terminating the method (step 308 ).
  • the described X-ray device and the described versions of the method according to the invention have been given merely by way of example.
  • the X-ray device may also have a different construction.
  • Individual steps of the method for example the acquisition of X-ray projection images, may also be performed in a different manner in practice.

Abstract

The invention relates to a method of imaging the blood flow as a function of time in an object (3) to be examined, which method includes the following steps:
a) acquisition of a series of X-ray projection images (Di; Ej) during administration of a contrast medium to the blood vessels in the object (3) to be examined,
b) acquisition of an image data set (H; K) containing the course of the blood vessels in the object (3) to be examined,
c) segmentation of the regions of the blood vessels in the individual X-ray projection images (Di; Ej) that are filled with contrast medium,
d) encoding the image data set (H; K) in time by comparing the image data set (H; K) with the segmented X-ray projection images (Dj′; Fj), and
e) displaying one or more images (B) formed from the time-encoded image data set (H′; R′) and representing the blood flow as a function of time.
The invention also relates to a correspondingly constructed X-ray device.

Description

  • The invention relates to a method of imaging the blood flow as a function of time in an object to be examined, as well as to an X-ray device for carrying out this method. [0001]
  • EP 860 696 A2 discloses a method which enables the reproduction of distributed structures, for example a vascular system filled with a contrast medium, in a synthetic projection image which reproduces the parts of the vascular system that are present in a selectable sub-volume more clearly than the X-ray images that are acquired from different perspectives and wherefrom this projection image is derived. This method thus yields a time-averaged “frozen” image of the vascular system in which the blood flow as a function of time is not visible. [0002]
  • Three-dimensional rotation angiography produces a series of X-ray projection images of the object to be examined from different projection directions while a contrast medium is injected into the blood vessels of the object to be examined. Using a known reconstruction algorithm, for example the Feldkamp algorithm, a three-dimensional image is derived from such X-ray projection images; this three-dimensional image reproduces the vascular system in space in a time-averaged manner. In two-dimensional angiography it is also known to form a two-dimensional “frozen” image of the vascular system. [0003]
  • For various applications, such as the analysis of pathologies of the cerebral vessels, for example in the case of vascular anomalies (stenoses, arteriovenous deformations), however, it is important to know and reproduce the blood flow as a function of time. Therefore, it is an object of the invention to provide a method which enables the reproduction of the blood flow in an object to be examined as a function of time, and to provide an appropriate X-ray device for carrying out this method. [0004]
  • Such a method should also take into account the fact that a contrast medium cannot be repeatedly injected into the same patient within a short period of time, so that the method should do, if at all possible, with a single contrast medium injection. Moreover, it should be possible to acquire the images in an as short as possible period of time, at an as small as possible expenditure and with an as high as possible resolution. [0005]
  • This object is achieved by means of the method disclosed in [0006] claim 1 and the X-ray device disclosed in claim 10.
  • The invention is based on the recognition of the fact that an image data set, which may be a two-dimensional or a three-dimensional image data set and contains information concerning the course of the blood vessels in the object to be examined, can be encoded in time in such a manner that it also contains information concerning the blood flow as a function of time. Such encoding in time is performed according to the invention in that the image data set is compared with a series of X-ray projection images; these X-ray projection images are formed successively in time and contain the information concerning the distribution of an injected contrast medium in the blood vessels at each time a different instant. Because each X-ray projection image is individually compared with the image data set, that is, each image value of the image data set is compared with the image values of the individual X-ray projection images, it is quasi checked which parts of the vascular system contained in the image data set are filled with the contrast medium at the individual instants associated with the respective X-ray projection images. Using suitable reproduction methods, the image data set thus encoded in time can be converted into one or more images which show the blood flow as a function of time. [0007]
  • The version disclosed in [0008] claim 2 is particularly suitable for two-dimensional rotation angiography. It already offers a two-dimensional X-ray image data set, for example a two-dimensional X-ray projection image which completely contains the vascular system filled with contrast medium and may be a previously formed or an instantaneous X-ray image data set. The actual X-ray projection images derived from a fixed X-ray position during a contrast medium administration are subtracted from one another in conformity with this version, so that each difference image contains the information as regards the path followed by the contrast medium in the vascular system between the two instants at which the two subtracted X-ray projection images were acquired. Such difference images are then used for the time encoding of the X-ray image data set.
  • A further version of this method is disclosed in [0009] claim 3. This version represents a simple possibility for comparing the image data set with the difference images. In order to enhance the imaging precision, first an X-ray image data sub-set is segmented from the X-ray image data set; this sub-set contains only the information concerning the course of the blood vessels. Subsequently, for each difference image there is acquired an associated pixel sub-set in such a manner that the pixels of the individual difference images are compared with the pixels of the X-ray image data sub-set and that with each pixel sub-set there are associated those pixels of the X-ray image data sub-set for which the associated difference image includes corresponding pixels. Each pixel sub-set thus contains the information concerning the distribution of the contrast medium at a given instant and one or more projection images which show the blood flow as a function of time can be derived from the pixel sub-sets.
  • According to the preferred version disclosed in claim [0010] 4, the X-ray projection images are acquired from different positions. The image data set then constitutes a three-dimensional X-ray image data set derived from such X-ray projection images. Thus, according to this further version only a single series of X-ray projection images is acquired from different directions, said images also containing the time information. After segmentation of the blood vessels in said X-ray projection images they are compared with the X-ray image data set and encoded in time, for example, in that pseudo-projection images are calculated from the X-ray image data set, utilizing the known imaging geometry of the X-ray device, so as to be compared with the actual X-ray projection images.
  • claim [0011] 5 discloses a particularly attractive further version which is also suitable for three-dimensional rotation angiography. Two series of X-ray projection images are acquired therein, the first series being acquired from a fixed X-ray position whereas the second series is acquired from different X-ray positions. This operation can be performed either successively in time, be it that two contrast medium injections are then required, or simultaneously by means of an X-ray device which includes two imaging units. From the X-ray projection images of the first series there are derived difference images which, as described above, contain the time information whereas a three-dimensional X-ray image data set is acquired from the X-ray projection images of the first series by means of a known reconstruction algorithm. This X-ray image data set is encoded in time by means of the difference images. This further version, notably the formation of difference images as carriers of the time information, offers the advantage that the difference images always contain only the variation in time of the contrast medium flow (=the blood flow) in a given time interval whereas in the case of direct use of the X-ray projection images as carriers of the time information the overall distribution of the contrast medium in the vascular system is always contained therein and hence the information is significantly less exact. This is also due to the fact that the contrast medium propagates very quickly throughout the vascular system to be observed, so that the differences in the distribution of the contrast medium between two comparatively closely spaced instants are only comparatively small.
  • The claims [0012] 6 and 7 disclose preferred possibilities for the time encoding of the image data set and the comparison of the image data set with the segmented X-ray projection images. The steps of the method as disclosed in claim 7 correspond to the method which is known from EP 880 109 A2 whose disclosure is explicitly referred to herein and is considered to be incorporated herein. The latter publication describes a method of determining the spatial transformation between a three-dimensional object reproduced by a data set and the object itself. A pseudo-projection image is then calculated for a part of the volume reproduced by the data set; this pseudo-projection image is compared with an X-ray projection image of the object itself. The parameters on which the calculation of the pseudo-projection image is based are then varied until an optimum match is obtained. This method can be advantageously used in an appropriate manner in conjunction with the present invention.
  • The claims [0013] 8 and 9 disclose advantageous further versions concerning the display of the blood flow as a function of time. For example, the time information can be converted into a color code so that the complete two-dimensional or three-dimensional data set can be reproduced as an image with the corresponding color code. The encoded pixel sub-sets or voxel sub-sets can also be reproduced successively in time, for example as an endless loop, thus creating the impression of blood flowing through the blood vessels. It is also feasible to enable reproductions from arbitrary angles of observation and rotation of the image is also possible.
  • An X-ray device which is suitable for carrying out the method according to the invention is disclosed in [0014] claim 10; as is indicated in the embodiment of claim 11, it may also include a second imaging unit.
  • The invention will be described in detail hereinafter with reference to the drawings. Therein: [0015]
  • FIG. 1 shows an X-ray device for carrying out the method according to the invention, [0016]
  • FIG. 2 shows a flow chart illustrating a first version of the method according to the invention, [0017]
  • FIG. 3 illustrates diagrammatically the encoding in time, [0018]
  • FIG. 4 shows a flow chart of a second version of the method according to the invention, and [0019]
  • FIG. 5 shows a flow chart of a third version of the method according to the invention. [0020]
  • FIG. 1 shows an [0021] X-ray device 1 which serves to form two-dimensional X-ray images of an object 3 to be examined, for example, a patient who is arranged on a table 4. The X-ray device 1 includes an X-ray source 12 and an X-ray detector 13 which are aligned relative to one another and mounted on a C-arm 10 which itself is journalized in a stand 11 which is only partly shown. The C-arm 10 can be pivoted about a perpendicular axis on the one hand and be rotated around its center in the direction of the double arrow 20 on the other hand, for example through 180° C., by means of a motor drive (not shown). During this motion a plurality of X-ray images can be acquired so as to image the object 3 to be examined from different, reproducible perspectives or X-ray positions r1, r2, r3 of the imaging unit 12, 13. Also provided are a second X-ray source 12′ and a second X-ray detector 13′ which are mounted on a supporting device 11′ and are capable of forming X-ray images of the object 3 to be examined from a fixed X-ray position r4.
  • Each of the [0022] X-ray detectors 13, 13′ may be formed by a respective X-ray image intensifier whereto there is connected a television chain whose output signals are digitized by an analog-to-digital converter 14 so as to be stored in a memory 15. The X-ray projection images D1, . . . Di, . . . Dim, acquired by the imaging unit 12, 13 from different X-ray positions r1, . . . ri, . . . rm (of which only the positions r1, r2, r3 are explicitly shown in the drawing), can be processed by an image processing unit 16 so as to be reproduced on a monitor 18 either individually or as a series of images. The X-ray projection images E0, E1, E2 . . . , En, acquired by the imaging unit 12′, 13′ from the fixed X-ray position r′ at discrete instants t during a contrast medium bolus T, can also be processed by the image processing unit 16 so as to be displayed on the monitor 18. The individual components of the X-ray device are controlled by means of a control unit 17.
  • Also provided is an arithmetic unit [0023] 19 which receives the X-ray projection images Di and Ej in order to derive therefrom images showing the variation in time of the blood flow in the relevant vessels of the object 3 to be examined. These images can be displayed on the monitor 18 again.
  • The invention will be illustrated hereinafter on the basis of the flow chart of a first version of the method according to the invention which is shown in FIG. 2. After the initialization (step [0024] 100) and after a contrast medium injection, the C-arm 10 is step-wise rotated about its center and at the same time a series of m X-ray projection images Di (for example, m=100) is formed, which projection images represent the object 3 to be examined and the blood vessels that are present therein and are filled with a contrast medium, from different perspectives (step 101). The X-ray projection images Di constitute a three-dimensional X-ray image data set K.
  • During a [0025] step 102 which takes place at the same time or at a later instant (after a further contrast medium injection), the imaging unit 12′, 13′ acquires a second series of X-ray projection images Ej from a fixed perspective. During the subsequent step 103 on the one hand correction is made for imaging errors which are due, for example, to the imperfection of the imaging device or to the mechanical deformation of the C-arm. On the other hand a three-dimensional reconstruction image R is formed from the X-ray projection images Di by means of a known reconstruction algorithm. In the step 104 a reconstruction sub-image R′ is formed from said reconstruction image R, which reconstruction sub-image contains, mainly or exclusively, information concerning the course of the blood vessels in the region examined. This reconstruction sub-image R′ contains q voxels Vk which are characterized by their co-ordinates in space.
  • During a [0026] further step 105 respective difference images Fj are determined from each time two successively acquired X-ray projection images Ej-1 and Ej; such difference images contain only the information concerning the path traveled by the contrast medium in the period of time elapsing between the formation of the two X-ray projection images Ej-1 and Ej.
  • Subsequently, each of the n difference images F[0027] j is compared with the reconstruction sub-image R′, so that the latter is encoded in time. To this end, in the step 106 each of the q voxels Vk of the reconstruction sub-image R′ is projected onto each difference image Fj, that is, a pseudo-projection image is calculated from the reconstruction sub-image R′ so as to be compared with the individual difference images Fj. For the calculation of this pseudo-projection image it is necessary to know the geometry of the imaging unit 12′, 13′ in order to ensure that the co-ordinate system of the pseudo-projection image corresponds to the co-ordinate system of the X-ray projection images Ej and the difference images Fj.
  • Due to the comparison of the voxels V[0028] k of the reconstruction sub-image R′, or the pixels of the pseudo-projection image, with the pixels P1 of the individual difference images Fj, the voxels of the reconstruction sub-image R′ for which the corresponding difference image Fj comprises a corresponding pixel are marked each time. The marked voxels are combined so as to form a voxel sub-set Lj in the step 107.
  • In simplified form it may be stated that in the [0029] steps 106 and 107 it is checked to what segment of the vascular system the contrast medium has progressed in a time interval between the formation of two X-ray projection images Ej-1 and Ej, or from what segment of the vascular system the contrast medium has drained during this time interval, and that the voxels in the reconstruction sub-image R′ which represent the relevant vessel segment are marked accordingly and combined so as to form a voxel sub-set Lj. These steps are carried out n times, that is, once for every difference image Fj. This results in n voxel sub-sets Lj wherefrom one or more images B can be formed in the step 108, said images representing the blood flow as a function of time. For example, the voxels of each voxel sub-set Lj can be reproduced in a different color in an overall image. The method is terminated in the step 109.
  • Notably the [0030] step 106 will be illustrated again with reference to FIG. 3. After a reconstruction sub-image R′ has been determined from the X-ray projection images Di in several steps (which sub-image contains only the information concerning the vascular structure), the pixels Vk of this reconstruction sub-image R′ are projected onto the difference images Fj. Upon projection onto the first difference image F1, the voxel V1 along the projection ray 21 is projected onto the pixel P1. On the basis of the grey tone, intended to illustrate that a voxel or a pixel represents a part of a vascular structure, as can be detected, for example, on the basis of the image value which lies above or below a given threshold value, it can be recognized that in the difference image F1 a corresponding pixel P1 exists for the voxel V1. This means that the segment of the vascular system represented by the voxel V1 has been filled with contrast medium in the time interval corresponding to the difference image F1. Therefore, the voxel V1 is marked as symbolically indicated by the arrow 22.
  • The same holds for the voxel V[0031] 2 which is projected onto the pixel P2 and is also marked. The voxel V3, projected onto the pixel P3, however, is not marked because the pixel P3, or the part of the object to be examined which is represented by this pixel, has not been filled with contrast medium in the difference image F1.
  • In conformity with this method all voxels V[0032] k of the reconstruction sub-image R′, or at least all voxels reproducing the vascular system, are successively compared with the pixels of each individual difference image Fj so as to be marked. In the example shown the voxels V1 and V2, being the only voxels marked after the comparison with the difference image F1, form the voxel sub-set L1 which can be represented by a first color in the three-dimensional reconstruction image formed upon completion of the method. The voxel sub-set L1 thus contains all voxels into which contrast medium has flown in the time interval between the formation of the X-ray projection images E0 and E1.
  • A reconstruction sub-image R′ can also be obtained in a manner other than that described with reference to FIG. 2, for example, by subtraction angiography. To this end, two series of X-ray projection images are formed from different X-ray positions by means of the [0033] imaging unit 12, 13, that is, once with and once without administration of contrast medium. The X-ray projection images of the two series that are formed from each time the same X-ray position are then subtracted from one another and the resultant projection images are reconstructed so as to form the desired reconstruction sub-image R′.
  • FIG. 4 shows the flow chart of a further version of the method according to the invention. This flow chart deviates from that of the version shown in FIG. 2 first of all in that only a single series of X-ray projection images D[0034] i is acquired; no X-ray projection images Ej are acquired from a fixed X-ray position, so that the second imaging unit 12′, 13′ (see FIG. 1) can also be dispensed with. The steps 200 to 203 otherwise correspond without modification to the steps 100 to 104 with omission of the step 102.
  • During the [0035] subsequent step 204 the X-ray projection images Di are segmented, that is, the pixels representing the vascular structures in the individual images Di are selected and reproduced exclusively in the segmented X-ray projection images Di′, while other image elements such as, for example the background or other organs, are eliminated. This operation can be performed, for example, by means of the image values which lie above or below a given limit value for the vessels filled with contrast medium whereas this is usually not the case for other image elements.
  • The successive segmented X-ray projection images Di are subsequently individually compared with the reconstruction sub-image R′ or the pixels P[0036] 1 of the segmented X-ray projection images Di′ are compared with the voxels Vk of the reconstruction sub-image R′. This operation takes place in the same way as described with reference to FIG. 3 for the difference images Fj. The voxel sub-sets Lj and the images B for reproducing the blood flow are formed in the steps 206 and 207 in the same way as already described above. The method is completed again in the step 208.
  • In comparison with the version described with reference to FIG. 2 the present version offers the advantages that only one imaging unit ([0037] 12, 13) is required and that only a single series of X-ray projection images Di need be acquired. The segmentation of the images Di in the step 204, however, may be difficult in practice and could give rise to inaccuracies. Moreover, the information contents of the segmented X-ray projection images Di′ differ from that of the difference images Fj which contain the variation of the blood flow as a function of time in a time interval whereas the segmented X-ray projection images Di′ show the entire region of the vascular system filled with the contrast medium.
  • A third version of the method according to the invention will be described in detail with reference to FIG. 5. After the start in the [0038] step 300, a two-dimensional X-ray image data set H is acquired in the step 301. This may be a single X-ray projection image which reproduces the complete vascular system filled with a contrast medium, or an image composed of a plurality of individual X-ray projection images. The data set H may also have been acquired already at an earlier instant. During the step 302 a series of X-ray projection images Ej is acquired (after contrast injection) from a fixed X-ray position by means of the imaging unit (12, 13). At the same time, or subsequently, the two-dimensional image data set H is segmented in the step 303, so that the segmented image data set H′ contains only the vessel structures.
  • In the [0039] step 304 respective difference images Fj, having the already described information contents are determined from each time two temporally successively acquired X-ray projection images Ej-1, Ej. These difference images Fj are then individually compared successively with the segmented image data set, that is, the pixels P1 of each individual difference image Fj are compared with the pixels Pk of the segmented two-dimensional image data set H′. To this end it is necessary that the X-ray projection images Ej and the two-dimensional image data set have been acquired by means of the same imaging unit and from the same X-ray position. This comparison in the step 305 produces in an associated pixel sub-set Mj for each difference image Fj in the step 306; this pixel sub-set is comparable to the already described voxel sub-sets Lj and contains the information as to which pixels of the segmented image data set H′ represent a part of the vascular system which has been filled with contrast medium in a given time interval or wherefrom the contrast medium has been drained during this time interval.
  • After the [0040] steps 305 and 306 have been performed n times, the entire segmented image H′ will have been encoded in time and one or more images can be displayed in order to visualize the blood flow in this image as a function of time (step 307), thus terminating the method (step 308).
  • The described X-ray device and the described versions of the method according to the invention have been given merely by way of example. The X-ray device may also have a different construction. Individual steps of the method, for example the acquisition of X-ray projection images, may also be performed in a different manner in practice. [0041]

Claims (11)

1. A method of imaging the blood flow as a function of time in an object (3) to be examined, which method includes the following steps:
a) acquisition of a series of X-ray projection images (Di; Ej) during administration of a contrast medium to the blood vessels in the object (3) to be examined,
b) acquisition of an image data set (H; K) containing the course of the blood vessels in the object (3) to be examined,
c) segmentation of the regions of the blood vessels in the individual X-ray projection images (Di; Ej) that are filled with contrast medium,
d) encoding the image data set (H; K) in time by comparing the image data set (H; K) with the segmented X-ray projection images (Dj′; Fj), and
e) displaying one or more images (B) formed from the time-encoded image data set (H′; R′) and representing the blood flow as a function of time.
2. A method as claimed in
claim 1
, characterized in that the X-ray projection images (Ej) are acquired from a fixed X-ray position (r4), that the image data set is a two-dimensional X-ray image data set (H), that respective difference images (Fj) are formed from each time two temporally successive X-ray projection images (Ej) in order to segment the X-ray projection images (Ej), and that the encoding in time is performed by comparing the two-dimensional X-ray image data set (H) with the difference images (Fj).
3. A method as claimed in
claim 2
, characterized in that an X-ray image data sub-set (H′), containing the information concerning the course of the blood vessels, is segmented from the X-ray image data set (H), that the pixels (Pk) of the X-ray image data sub-set (H′) are projected onto the difference images (Fj), that associated pixel sub-sets (Mj) are acquired for the individual difference images (Fj), that those pixels (Pk) of the X-ray image data sub-set (H′) for which the associated difference image (Fj) comprises corresponding pixels (P1) are assigned to a pixel sub-set (Mj), and that one or more two-dimensional projection images (B) are formed from the pixel sub-sets (Pj).
4. A method as claimed in
claim 1
, characterized in that the X-ray projection images (Di) are acquired from different X-ray positions (r1, r2, r3), and that the image data set is a three-dimensional X-ray image data set (K) derived from said X-ray projection images (Di).
5. A method as claimed in
claim 1
, characterized in that the X-ray projection images (Ej) are acquired from a fixed X-ray position (r4), that a second series of X-ray projection images (Di) is acquired from different X-ray positions (r1, r2, r3), that the image data set is a three-dimensional X-ray image data set (K) derived from the second series of X-ray projection images (Di), and that respective difference images (Fj) are derived from each time two temporally successive X-ray projection images (Ej) in order to segment the X-ray projection images (Ej) of the first series.
6. A method as claimed in
claim 4
or
5
, characterized in that a three-dimensional reconstruction image (R) is derived from the X-ray image data set (K), that a reconstruction sub-image (R′) which contains essentially the blood vessels is derived from the reconstruction image (R), and that the encoding in time of the image data set (K) is performed by comparison of the reconstruction sub-image (R′) with the segmented X-ray projection images (Ej′; Fj).
7. A method as claimed in
claim 6
, characterized in that the voxels (Vk) of the reconstruction sub-image (R′) are projected onto the segmented X-ray projection images (Ej′; Fj) in order to compare the reconstruction sub-image (R′) with the segmented X-ray projection images (Ej′; Fj), that associated voxel sub-sets (Lj) are determined for the individual segmented X-ray projection images (Ej′; Fj), that those voxels (Vk) for which the associated X-ray projection image (Ej′; Fj) contains corresponding pixels (P1) are assigned to a voxel sub-set (Lj), and that one or more three-dimensional reconstruction images (B) are formed from the voxel sub-sets (Lj).
8. A method as claimed in
claim 3
or
7
, characterized in that the pixel sub-sets (Pj), or the voxel sub-sets (Lj), are reproduced in a different color in the projection image (B) or the reconstruction image (B).
9. A method as claimed in
claim 3
or
7
, characterized in that a respective reconstruction image (B) is formed from each pixel sub-set (Pj), or voxel sub-set (Lj), and that such reconstruction images (B) are successively displayed in rapid succession.
10. An X-ray device for carrying out the method claimed in
claim 1
, which device includes:
an imaging unit (1) with an X-ray source (12) and an X-ray detector (13) for the acquisition of a series of X-ray projection images (Di; Ej) during administration of a contrast medium to the blood vessels, and for the acquisition of an image data set (H; K) containing the course of blood vessels in the object (3) to be examined,
an arithmetic unit (19) for segmenting the regions of the blood vessels in the individual projection images (Di; Ej) that are filled with the contrast medium and for the encoding in time of the image data set (K) by comparison of the image data set (H; K) with the segmented X-ray projection images (Dj′; Fj), and
a display device for the display of one or more images (B) that are derived from the time encoded image data set (H′; R′) and represent the blood flow as a function of time.
11. An X-ray device as claimed in
claim 10
, characterized in that the imaging unit includes a first X-ray source (12) and a first X-ray detector (13) for the acquisition of a first series of two-dimensional X-ray images (Ej) from a fixed X-ray position (r4), and also includes a second X-ray source (12′) and a second X-ray detector (13′) for the acquisition of a second series of two-dimensional X-ray images (Ej) from different X-ray positions (r1, r2, r3).
US09/749,147 2000-01-05 2000-12-27 Method of imaging the blood flow as a function of time in an object to be examined Expired - Fee Related US6442235B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
DE10000185 2000-01-05
DE10000185A DE10000185A1 (en) 2000-01-05 2000-01-05 Method for displaying the time course of the blood flow in an examination object
DE10000185.8 2000-01-15

Publications (2)

Publication Number Publication Date
US20010012328A1 true US20010012328A1 (en) 2001-08-09
US6442235B2 US6442235B2 (en) 2002-08-27

Family

ID=7626779

Family Applications (1)

Application Number Title Priority Date Filing Date
US09/749,147 Expired - Fee Related US6442235B2 (en) 2000-01-05 2000-12-27 Method of imaging the blood flow as a function of time in an object to be examined

Country Status (5)

Country Link
US (1) US6442235B2 (en)
EP (1) EP1114615B1 (en)
JP (1) JP2001212125A (en)
CN (1) CN1212807C (en)
DE (2) DE10000185A1 (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070031018A1 (en) * 2005-08-03 2007-02-08 Siemens Aktiengesellschaft Operating method for an image-generating medical engineering assembly and articles associated herewith
US20070041625A1 (en) * 2005-08-18 2007-02-22 Siemens Aktiengesellschaft Image evaluation method for two-dimensional projection images and items corresponding thereto
US20080232548A1 (en) * 2007-03-20 2008-09-25 Kabushiki Kaisha Toshiba X-ray diagnostic system
US20090252287A1 (en) * 2008-04-02 2009-10-08 Siemens Aktiengesellschaft Operating method for a swiveling polyplanar imaging system for time-resolved imaging of an object being examined
WO2013038313A1 (en) * 2011-09-13 2013-03-21 Koninklijke Philips Electronics N.V. Vascular outlining with ostia visualization
US20130099127A1 (en) * 2011-10-24 2013-04-25 Siemens Aktiengesellschaft Method and Device for Detecting X-Ray Quanta
WO2014070191A1 (en) * 2012-11-02 2014-05-08 Analogic Corporation Volumetric and projection image generation

Families Citing this family (46)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4473358B2 (en) * 1999-01-21 2010-06-02 株式会社東芝 Diagnostic equipment
DE10100572A1 (en) * 2001-01-09 2002-07-11 Philips Corp Intellectual Pty Method for displaying blood flow in a vascular tree
EP1302163B1 (en) 2001-10-16 2006-07-05 Kabushiki Kaisha Toshiba Apparatus for calculating an index of local blood flows
US6512807B1 (en) * 2001-11-21 2003-01-28 Koninklijke Philips Electronics, N.V. Low signal correction for perfusion measurements
ATE376389T1 (en) * 2002-02-15 2007-11-15 Breakaway Imaging Llc GANTRY RING WITH REMOVABLE SEGMENT FOR MULTI-DIMENSIONAL X-RAY IMAGING
US7188998B2 (en) * 2002-03-13 2007-03-13 Breakaway Imaging, Llc Systems and methods for quasi-simultaneous multi-planar x-ray imaging
AU2003224711A1 (en) * 2002-03-19 2003-10-08 Breakaway Imaging, Llc Computer tomograph with a detector following the movement of a pivotable x-ray source
DE10224011A1 (en) * 2002-05-29 2003-12-24 Siemens Ag Computer-aided reconstruction method for a three-dimensional object
DE60315642T2 (en) * 2002-06-11 2008-06-05 Breakaway Imaging, LLC, Littleton OUTSTANDING GANTRY DEVICE FOR X-RAYING THROUGH X-RAYS
KR100639139B1 (en) * 2002-06-28 2006-10-30 후지쯔 가부시끼가이샤 Three-dimensional image comparing program, three-dimensional image comparing method, and three-dimensional image comparing device
US7338207B2 (en) * 2002-08-21 2008-03-04 Medtronic Navigation, Inc. Gantry positioning apparatus for X-ray imaging
AU2003262726A1 (en) * 2002-08-21 2004-03-11 Breakaway Imaging, Llc Apparatus and method for reconstruction of volumetric images in a divergent scanning computed tomography system
FR2847798B1 (en) * 2002-11-28 2006-02-10 Ge Med Sys Global Tech Co Llc METHOD FOR DETERMINING FUNCTIONAL PARAMETERS IN A FLUOROSCOPIC DEVICE
EP1617763A1 (en) * 2003-04-22 2006-01-25 Philips Intellectual Property & Standards GmbH Apparatus for angiographic x-ray photography
JP4308851B2 (en) * 2003-08-20 2009-08-05 コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ Method and apparatus for flow reconstruction
DE102004018499A1 (en) * 2004-04-14 2005-06-02 Siemens Ag 3D graphical image processing method, especially for medical applications, in which base data is processed to produce projection line data finally independent voxel data
CN101115989B (en) * 2005-02-04 2012-09-05 皇家飞利浦电子股份有限公司 System for the determination of vessel geometry and flow characteristics
EP1869643B1 (en) * 2005-04-07 2014-05-07 Koninklijke Philips N.V. Image processing device and method for blood flow imaging
DE102005042328A1 (en) * 2005-09-06 2007-03-08 Siemens Ag Method for determining the propagation behavior of a contrast agent bolus
DE102006003829A1 (en) * 2006-01-26 2007-08-16 Siemens Ag X-ray computed tomography and method of operating an X-ray CT scanner
DE102006025422B4 (en) 2006-05-31 2009-02-26 Siemens Ag Image evaluation method for two-dimensional projection images and objects corresponding thereto
DE102006025420B4 (en) 2006-05-31 2009-04-02 Siemens Ag Presentation method for two-dimensional projection images and objects corresponding thereto
CN101243980B (en) * 2006-12-04 2010-12-22 株式会社东芝 X-ray computed tomographic apparatus and medical image processing apparatus
DE102006058907A1 (en) * 2006-12-13 2008-06-26 Siemens Ag Angiography device for diagnosing vasculars of patient, filled with contrast agent, comprises X-ray source and associated detector with image processing unit, image display unit and control unit, which records images of cardiac cycle
US7680240B2 (en) * 2007-03-30 2010-03-16 General Electric Company Iterative reconstruction of tomographic image data method and system
DE102007024451B4 (en) * 2007-05-25 2010-07-08 Siemens Ag Determination method for spatially resolved three-dimensional occupancy distributions of a substance in a vascular system and facilities corresponding thereto
DE102007044406A1 (en) 2007-09-18 2009-03-19 Siemens Ag Registration method with circulation-dependent three-dimensional representation of a vascular tree
DE102007046281A1 (en) 2007-09-27 2009-04-09 Siemens Ag Method and arithmetic unit for measuring the flow rate of a contrast agent in a vessel of a patient
DE102007046514A1 (en) * 2007-09-28 2009-04-23 Siemens Ag Method for detecting and marking contrast medium in blood vessels of the lung using a CT examination and image evaluation unit of a CT system
US8615116B2 (en) * 2007-09-28 2013-12-24 The Johns Hopkins University Combined multi-detector CT angiography and CT myocardial perfusion imaging for the diagnosis of coronary artery disease
FR2924255A1 (en) * 2007-11-27 2009-05-29 Gen Electric METHOD FOR PROCESSING RADIOGRAPHIC CARDIAC IMAGES FOR OBTAINING A SUBTRACT AND RECALLED IMAGE
WO2009109887A1 (en) * 2008-03-06 2009-09-11 Philips Intellectual Property & Standards Gmbh Method for analyzing a tube system
DE102008023053A1 (en) 2008-05-09 2009-11-26 Siemens Aktiengesellschaft Evaluation method for two-dimensional fluoroscopic images of an examination object with time-coded representation of three-dimensional reconstructions
US8200466B2 (en) 2008-07-21 2012-06-12 The Board Of Trustees Of The Leland Stanford Junior University Method for tuning patient-specific cardiovascular simulations
DE102009004576B4 (en) 2009-01-14 2011-03-17 Siemens Aktiengesellschaft Catheter detection in volume data on probabilistic approach
US9405886B2 (en) 2009-03-17 2016-08-02 The Board Of Trustees Of The Leland Stanford Junior University Method for determining cardiovascular information
EP2408375B1 (en) 2009-03-20 2017-12-06 Orthoscan Incorporated Moveable imaging apparatus
US8155420B2 (en) * 2009-05-21 2012-04-10 Visiongate, Inc System and method for detecting poor quality in 3D reconstructions
DE102009031141B3 (en) * 2009-06-30 2010-12-23 Siemens Aktiengesellschaft Determination procedure for a color-coded evaluation image as well as corresponding objects
US8157742B2 (en) 2010-08-12 2012-04-17 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
US8315812B2 (en) 2010-08-12 2012-11-20 Heartflow, Inc. Method and system for patient-specific modeling of blood flow
WO2012082799A1 (en) 2010-12-13 2012-06-21 Orthoscan, Inc. Mobile fluoroscopic imaging system
US9107639B2 (en) * 2011-03-15 2015-08-18 Medicinsk Bildteknik Sverige Ab System for synchronously visualizing a representation of first and second input data
DE102011083685A1 (en) * 2011-09-29 2013-04-04 Siemens Aktiengesellschaft Imaging method for representation of flow velocities within objects in interventional angiographic examinations, involves subtracting several three-dimensional (3D) volume images of objects, for generation of 3D subtraction image
US8548778B1 (en) 2012-05-14 2013-10-01 Heartflow, Inc. Method and system for providing information from a patient-specific model of blood flow
US11069054B2 (en) 2015-12-30 2021-07-20 Visiongate, Inc. System and method for automated detection and monitoring of dysplasia and administration of immunotherapy and chemotherapy

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3018129C1 (en) * 1980-05-12 1981-10-01 Siemens AG, 1000 Berlin und 8000 München X-ray diagnostic facility for creating subtraction images
DE3632833A1 (en) * 1986-09-26 1988-03-31 Philips Patentverwaltung METHOD FOR PRODUCING A X-RAY LAYER IMAGE OF AN EXAMINATION AREA AND ARRANGEMENT FOR IMPLEMENTING THE METHOD
DE3710011A1 (en) * 1987-03-26 1988-10-06 Philips Patentverwaltung METHOD IN WHICH AN EXAMINATION AREA IS TRANSMITTED FROM SEVERAL RADIATION SOURCE POSITIONS
DE69325485T2 (en) * 1992-09-09 1999-10-28 Picker Int Inc Image forming method and apparatus
FR2708166A1 (en) * 1993-07-22 1995-01-27 Philips Laboratoire Electroniq A method of processing digitized images for the automatic detection of stenoses.
DE19620371A1 (en) * 1996-05-21 1997-12-04 Philips Patentverwaltung X-ray procedure
DE19705600A1 (en) * 1997-02-14 1998-08-20 Philips Patentverwaltung Spatial transformation determination method for X-ray imaging
DE19705599A1 (en) * 1997-02-14 1998-08-20 Philips Patentverwaltung X-ray imaging process with a series of exposures from different perspectives

Cited By (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7796796B2 (en) * 2005-08-03 2010-09-14 Siemens Aktiengesellschaft Operating method for an image-generating medical engineering assembly and articles associated herewith
US20070031018A1 (en) * 2005-08-03 2007-02-08 Siemens Aktiengesellschaft Operating method for an image-generating medical engineering assembly and articles associated herewith
US20070041625A1 (en) * 2005-08-18 2007-02-22 Siemens Aktiengesellschaft Image evaluation method for two-dimensional projection images and items corresponding thereto
US7729525B2 (en) * 2005-08-18 2010-06-01 Siemens Aktiengesellschat Image evaluation method for two-dimensional projection images and items corresponding thereto
US20080232548A1 (en) * 2007-03-20 2008-09-25 Kabushiki Kaisha Toshiba X-ray diagnostic system
US7668290B2 (en) * 2007-03-20 2010-02-23 Kabushiki Kaisha Toshiba X-ray diagnostic system
US20090252287A1 (en) * 2008-04-02 2009-10-08 Siemens Aktiengesellschaft Operating method for a swiveling polyplanar imaging system for time-resolved imaging of an object being examined
US7899151B2 (en) * 2008-04-02 2011-03-01 Siemens Aktiengesellschaft Operating method for a swiveling polyplanar imaging system for time-resolved imaging of an object being examined
WO2013038313A1 (en) * 2011-09-13 2013-03-21 Koninklijke Philips Electronics N.V. Vascular outlining with ostia visualization
US9256940B2 (en) 2011-09-13 2016-02-09 Koninklijke Philips N.V. Vascular outlining with ostia visualization
US20130099127A1 (en) * 2011-10-24 2013-04-25 Siemens Aktiengesellschaft Method and Device for Detecting X-Ray Quanta
US8981313B2 (en) * 2011-10-24 2015-03-17 Siemens Aktiengesellschaft Method and device for detecting x-ray quanta
WO2014070191A1 (en) * 2012-11-02 2014-05-08 Analogic Corporation Volumetric and projection image generation
US20150260875A1 (en) * 2012-11-02 2015-09-17 Analogic Corporation Volumetric and projection image generation
US9696452B2 (en) * 2012-11-02 2017-07-04 Analogic Corporation Volumetric and projection image generation

Also Published As

Publication number Publication date
CN1317292A (en) 2001-10-17
JP2001212125A (en) 2001-08-07
EP1114615A3 (en) 2003-05-14
DE10000185A1 (en) 2001-07-12
DE50109522D1 (en) 2006-05-24
EP1114615A2 (en) 2001-07-11
US6442235B2 (en) 2002-08-27
EP1114615B1 (en) 2006-04-19
CN1212807C (en) 2005-08-03

Similar Documents

Publication Publication Date Title
US6442235B2 (en) Method of imaging the blood flow as a function of time in an object to be examined
US6823204B2 (en) Method of imaging the blood flow in a vascular tree
JP4731476B2 (en) Reconstruction of moving object images from volumetric data
US8009885B2 (en) Image processing device and method for blood flow imaging
US6983182B2 (en) Time resolved computed tomography angiography
US7269246B2 (en) X-ray angiography apparatus
US7903856B2 (en) Method for post-processing a three-dimensional image data set of vessel structure
US7339587B2 (en) Method for medical imaging and image processing, computed tomography machine, workstation and computer program product
US8676297B2 (en) Method and device for separate three-dimensional presentation of arteries and veins in a part of body
EP1182613A2 (en) Diagnostic imaging
US20100201786A1 (en) Method and apparatus for reconstructing an image
US7852984B2 (en) Method and device for the separate three-dimensional representation of the arterial and venous vascular system using C-arm angiography systems
US8929633B2 (en) Diagnostic X-ray system and method
US5978439A (en) X-ray imaging method involving a series of images from different perspectives
US20030073892A1 (en) Imaging methods and apparatus particularly useful for two and three-dimensional angiography
US20100260404A1 (en) X-ray diagnosis apparatus and image reconstruction processing apparatus
US20090022271A1 (en) X-ray imaging apparatus
US6426994B1 (en) Image processing method
EP1658009B1 (en) Method and device for flow reconstruction
US6366635B1 (en) Method and apparatus for three-dimensional image-rendering of a spatial and tissue-based configuration through separating high contrast and injected contrast agents in multi-angular x-ray absorption measurement
JPH1119082A (en) Image reconstructing apparatus

Legal Events

Date Code Title Description
AS Assignment

Owner name: U.S. PHILIPS CORP., NEW YORK

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:KOPPE, REINER HEINRICH;KLOTZ, ERHARD PAUL ARTUR;KUHN, MICHAEL HARALD;AND OTHERS;REEL/FRAME:011637/0917;SIGNING DATES FROM 20010124 TO 20010129

AS Assignment

Owner name: KONINKLIJKE PHILIPS ELECTRONICS N.V., NETHERLANDS

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:U.S. PHILIPS CORPORATION;REEL/FRAME:013065/0041

Effective date: 20020701

FPAY Fee payment

Year of fee payment: 4

REMI Maintenance fee reminder mailed
LAPS Lapse for failure to pay maintenance fees
STCH Information on status: patent discontinuation

Free format text: PATENT EXPIRED DUE TO NONPAYMENT OF MAINTENANCE FEES UNDER 37 CFR 1.362

FP Lapsed due to failure to pay maintenance fee

Effective date: 20100827