US20100004526A1 - Abnormality finding in projection images - Google Patents
Abnormality finding in projection images Download PDFInfo
- Publication number
- US20100004526A1 US20100004526A1 US12/475,995 US47599509A US2010004526A1 US 20100004526 A1 US20100004526 A1 US 20100004526A1 US 47599509 A US47599509 A US 47599509A US 2010004526 A1 US2010004526 A1 US 2010004526A1
- Authority
- US
- United States
- Prior art keywords
- image
- lumen
- bolus
- mask
- 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.)
- Abandoned
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C07—ORGANIC CHEMISTRY
- C07D—HETEROCYCLIC COMPOUNDS
- C07D473/00—Heterocyclic compounds containing purine ring systems
- C07D473/02—Heterocyclic compounds containing purine ring systems with oxygen, sulphur, or nitrogen atoms directly attached in positions 2 and 6
- C07D473/04—Heterocyclic compounds containing purine ring systems with oxygen, sulphur, or nitrogen atoms directly attached in positions 2 and 6 two oxygen atoms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61P—SPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
- A61P25/00—Drugs for disorders of the nervous system
- A61P25/28—Drugs for disorders of the nervous system for treating neurodegenerative disorders of the central nervous system, e.g. nootropic agents, cognition enhancers, drugs for treating Alzheimer's disease or other forms of dementia
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61P—SPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
- A61P31/00—Antiinfectives, i.e. antibiotics, antiseptics, chemotherapeutics
- A61P31/12—Antivirals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61P—SPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
- A61P31/00—Antiinfectives, i.e. antibiotics, antiseptics, chemotherapeutics
- A61P31/12—Antivirals
- A61P31/14—Antivirals for RNA viruses
Definitions
- the present disclosure is directed to medical imaging systems. More specifically, the present disclosure is directed to systems and methods that alone or collectively facilitate real-time imaging.
- Coronary artery disease causes in excess of 1.5 million cases of myocardial infarction annually, and is the leading cause of death in the United States, resulting in more than 500,000 deaths per year.
- the accurate diagnosis and quantification of coronary artery disease is critical to subsequent treatment decisions.
- digital subtraction angiography remains the preferred procedure to help clinicians to make clinical-decisions in all kinds of vascular diseases.
- Digital subtraction angiography is a well-established modality for the visualization of blood vessels in the human body.
- a DSA procedure a first series of frames are taken which represents the body anatomy in static form. They are called mask frames. Then dye is injected in the body through a catheter to visualize the blood vessel. The dye-injected frames are called the bolus frames. A subtraction procedure is preformed between the averaged bolus frame and the averaged mask frame to get DSA image.
- the patient is positioned on the X-ray imaging system while an X-ray movie is acquired and a DSA image is generated for an interventional cardiologist or a radiologist.
- the DSA image mainly contains the dye-enhanced blood vessels and they appear dark. However, the artifacts due to patient motion and system noise frequently reduce the diagnostic value of the images.
- a pathological aspect is associated with a vascular area where there is a significant deviation from the diameter of the healthy vascular.
- a stenosis is associated with a significant narrowing of the vascular and is quantified by parameters such as the percentage of stenosis.
- Stenosis limits blood flow by raising the resistance to flow through the vessel.
- the contrast agent ensures that the outlines of the blood flow are revealed on the X-ray to indicate any narrowing of the blood vessel.
- the present invention includes a system and method (i.e., utility) for determining the severity of a stenosis in a blood vessel.
- a method for improving DSA image quality includes: (1) registration of the mask and bolus images prior to subtracting procedure to reduce the artifacts from misalignment; (2) enhancement of the registered DSA image; and (3) detecting the boundary of blood vessel and quantitatively measuring percentage stenosis, which may be done automatically. This last step is sometimes referred to as Quantitative coronary angiography (QCA). Aspects of the present invention allow such QCA to be performed by a computer with minimal user input.
- QCA Quantitative coronary angiography
- the utility allows for the semi-automated quantitative measurement (QCA) of a lumen (e.g., blood vessel or artery).
- QCA semi-automated quantitative measurement
- the utility involves a registration step of the mask and the bolus images prior to subtracting procedure and QCA measurement.
- registration is a motion compensation step can reduce the artifacts from misalignment and improve the accuracy of QCA measurement.
- movement of the heart is nearly constant and motion compensation significantly improves overall image quality for QCA purposes.
- the motion compensation is performed using an algorithm that involves an inverse-consistency constrain which implies that the correspondence provided by the registration in one direction matches closely with the correspondence in the opposite direction.
- This may entail a B-spline parameterization.
- Other image registration methods can also be applied in the proposed system to match the mask and bolus images.
- the registration is performed hierarchically using a multi-resolution strategy in both, spatial domain and in domain of basis functions.
- the registration can be performed at 1 ⁇ 4, 1 ⁇ 2 and full resolution using knot spacing of 8, 16 and 32.
- the multi-resolution strategy helps in improving the registration by matching global structures at lowest resolution and then matching local vessel structures as the resolution is refined. Due to the bi-directional approach, thin blood vessels edges are more prominent thereby avoiding the local minima both in iterative nature of thin vessel edge estimation and iterative nature of correction of bolus image.
- the subtraction process is applied on corrected bolus images.
- Image enhancement prior to QCA measurement also improves the overall QCA measurement.
- the registered DSA image is further enhanced by background diffusion to remove noise and nonlinear normalization for better visualization.
- This image enhancement increases the contrast between the blood vessels and the background. This results in much improved contrast and very crisp subtraction images, in which the regions of interest are easily identifiable. Enhancement of the registered DSA image can be performed using an anisotropic diffusion technique and a nonlinear normalization technique. Other similar image enhancement techniques can also be applied in the proposed utility to enhance DSA images before QCA measurement.
- the utility allows a user to input various types of lumen identification a center line method where a user identifies an approximate center of the lumen and an edge method where a user identifies initial edges of the lumen. Both methods can obtain percentage stenosis automatically with minimal user interaction.
- the QCA measurement may include one or more of the following three sub-processes: (1) initial edge detection; (2) edge refinement; and (3) lesion measurement. Initial edge detection is performed where a user inputs centerline information.
- the utility may refine the edges of the lumen.
- an active contour model algorithm may be used to refine the edges by deforming a contour to lock onto features of interest within in an image.
- Other edge detection or segmentation techniques can also be applied in the system for QCA measurement.
- stenosis or lesion For each stenosis or lesion, the following data are calculated: (a) minimum and maximum diameter; (b) lesion length (the beginning and end of the lesion are determined automatically); (c) Stenosis reference (“normal”) diameter; and (d) percent stenosis. Other measurements can also be included in the proposed system.
- the entire software architecture may be implemented on a GPU based framework, which makes visualization and computation faster by up to factor of 30.
- the entire utility may implemented substantially in real-time.
- the utility may be used with various imaging modalities and may be used with 2-D or 3-D images.
- FIG. 1 illustrates one embodiment of a X-ray DSA system.
- FIG. 2 illustrates a process flow diagram for performing an improved QCA measurement in accordance with various aspects of the present invention.
- FIG. 3 illustrates one non-limiting method for registering images.
- FIG. 4 illustrates one non-limiting method for enhancing an image.
- FIG. 5 illustrates a process flow diagram for performing a QCA measurement.
- FIGS. 6A-6C illustrate an image having a vein or artery with a lesion.
- FIG. 7 illustrates a process flow diagram of an edge detection, refinement and lesion measurement process.
- FIGS. 8A-8C illustrate the refinement of edge boundaries of the vein or artery of FIGS. 6A-6C .
- FIG. 9 illustrates a process flow sheet for an edge detection process for detecting edge surfaces of an object in an image.
- FIG. 10 illustrates a process flow sheet for an edge refinement process.
- FIG. 11 illustrates a process flow sheet for a lesion measurement process.
- angiography may be performed using a number of different medical imaging modalities, including biplane X-ray/DSA, magnetic resonance (MR), computed tomography (CT), ultrasound, and various combinations of these techniques.
- MR magnetic resonance
- CT computed tomography
- ultrasound various combinations of these techniques.
- the following description is presented for purposes of illustration and description. Furthermore, the description is not intended to limit the invention to the form disclosed herein.
- FIG. 1 shows one exemplary setup for a real-time imaging procedure for use during a contrast media/dye injection procedure.
- a patient is positioned on an X-ray imaging system 100 and an X-ray movie is acquired by a movie acquisition system 102 .
- An enhanced DSA image is generated by an enhancement system 104 for output to a display 106 and/or stenosis processing system that is accessible to an interventional radiologist.
- the projection images are acquired at different time instants and consist of a movie with a series of frames before, during and after the dye injection.
- the series of frames include mask images that are free of contrast-enhancing dye in their field of view 108 and bolus images that contain contrast-enhancing dye in their field of view 108 . That is, bolus frames are images that are acquired after injected dye has reached the field of view 108 .
- the movie acquisition system 102 is operative to detect the frames before and after dye injection automatically to make feasible a real-time acquisition system.
- One approach for identifying frames before and after dye injection is to find intensity differences between successive frames, such that a large intensity difference is detected between the first frame after dye has reached the field of view (FOV) and the frame acquired before it.
- FOV field of view
- the patient may undergo some motion during the image acquisition causing such an intensity difference between even successive mask images.
- the movie acquisition system 102 may align successive frames together, such that the motion artifacts are minimized.
- the first image acquired after the dye has reached the FOV will therefore cause a high intensity difference with the previous frame not containing the dye in FOV.
- the subtraction image or ‘DSA image’ obtained by subtracting a mask frame from a bolus frame (or vice versa) will contain a near-zero value everywhere if both images belong to background.
- the subtraction image or DSA image is obtained by computing a difference between pixel intensities of the mask image and the bolus image.
- the enhancement system 104 may then enhance the contrast of the subtraction image.
- Such enhancement may include rescaling the intensities of the pixels in the subtraction image and/or the removal of noise from the subtraction image.
- the acquisition and/or enhancement systems are computerized systems that may run application software and computer programs which can be used to control the system components, provide user interfaces, and/or provide features of the imaging system.
- the software may be originally provided on computer-readable media, such as compact disks (CDs), magnetic tape, or other mass storage medium. Alternatively, the software may be downloaded from electronic links such as a host or vendor website. The software is installed onto a hard drive and/or electronic memory of the system, and is accessed and controlled by the operating system. Software updates are also electronically available on mass storage media or downloadable from the host or vendor website.
- the software, as provided on the computer-readable media or downloaded from electronic links represents a computer program product usable with a programmable computer processor having computer-readable program code embodied therein.
- the software contains one or more programming modules, subroutines, computer links, and compilations of executable code, which perform the functions of the imaging system.
- the user interacts with the software via keyboard, mouse, voice recognition, and other user-interface devices (e.g., user I/O devices) connected to the computer system.
- FIG. 2 shows a method for improving DSA image quality and improving the accuracy of QCA measurements that includes three primary processes:
- the first process includes using an averaging system 202 to generate average mask frames 204 and average bolus frames 206 from an x-ray image or movie. These average images are motion compensated 208 .
- the averaged mask frame is registered with the averaged bolus frame to generate a registered mask frame 210 , which minimizes the motion artifacts before subtraction.
- the registered mask frame 210 is then subtracted 212 from the average bolus frame to generate a DSA image 214 .
- Various aspects of this process are set forth in co-pending U.S. application Ser. No.
- the second process uses the registered DSA image 214 as input.
- the result of this is the production of an enhanced DSA image 218 .
- the third process uses the enhanced DSA image 218 to do a QCA measurement 220 and generate a QCA report or value 222 .
- the motion compensation process is discussed herein relation to FIG. 3
- the image enhancement process is discussed in relation to FIG. 4
- the QCA measurement is discussed in relation to FIGS. 5-11 .
- FIG. 3 illustrates the process of registering the averaged mask frame with averaged bolus frame, to minimize the motion artifacts.
- this is done using an inverse consistent image registration using B-spline basis.
- image registration is used to find a point-wise correspondence between a pair or group of similar anatomical objects.
- Image registration is often posed as an optimization problem that minimizes an objective function representing the difference between two images 302 , 304 to be registered.
- the images are the average bolus and mask frames.
- the symmetric squared intensity difference is chosen as the driving function.
- regularization constraints are applied so that the deformation follows a model that matches closely with the deformation of real-world objects.
- the regularization is applied in the form of bending energy and inverse-consistency cost.
- Inverse-consistency implies that the correspondence provided by the registration in one direction matches closely with the correspondence in the opposite direction.
- Most image registration methods are uni-directional and therefore contain correspondence ambiguities originating from choice of direction of registration.
- the forward and reverse correspondences are computed simultaneously and bound together through an inverse consistency cost term.
- the inverse consistency cost term assigns higher cost to transformations deviating from being inverse-consistent. While inverse consistency minimizes the correspondence ambiguity, it also helps the transformation perform better by forcing it out of local minima.
- a cost function for performing image registration over the images is used:
- I 1 (x) and I 2 (x) represent the intensity of image at location x
- x represents the domain of the image.
- L is a differential operator and the second term in Eq. (1) represents an energy function.
- ⁇ , ⁇ and ⁇ are weights to adjust relative importance of the cost function.
- the first term represents the symmetric squared intensity cost function and represents the integration of squared intensity difference between deformed reference image and the target image in both directions.
- the second term represents the energy regularization cost term and penalizes high derivatives of u(x).
- the third and last term represents the inverse consistency cost function, which penalizes differences between transformation in one direction and inverse of transformation in opposite direction.
- the total cost 308 is computed 306 using Eq. 1 as a first step in registration.
- the optimization problem posed in Eq. (2) is solved by using a B-spline parameterization as is known in the art.
- B-splines are used due to ease of computation, good approximation properties and their local support. It is also easier to incorporate landmarks in the cost term if we use spatial basis function.
- the above optimization problem is solved by solving for B-spline coefficients c i 's, such that
- h ⁇ ( x ) x + ⁇ i ⁇ c i ⁇ ⁇ i ⁇ ⁇ ( x ) ( 2 )
- ⁇ i (x) represents the value of B-spline at location x, originating at index i.
- cubic B-splines are used.
- a gradient descent scheme is implemented based on the above parameterization.
- the total gradient cost is calculated and updated 310 with respect to the transformation parameters in every iteration.
- the transformation parameters are updated using the gradient descent update rule. Images are deformed into shape of one another using the updated correspondence and the cost function and gradient costs are calculated 314 until convergence 316 when the frames are registered 318 .
- the registration is performed hierarchically using a multi-resolution strategy in both, spatial domain and in domain of basis functions.
- the registration is performed at 1 ⁇ 4, 1 ⁇ 2 and full resolution using knot spacing of 8, 16 and 32.
- the multi-resolution strategy helps in improving the registration by matching global structures at lowest resolution and their matching local structures as the resolution is refined.
- these registered frames are subtracted to provide the DSA image for image enhancement. That is, after image registration, subtraction is performed between the registered mask image and bolus image, resulting in the motion compensated (registered) DSA image 214 .
- the system uses the registered DSA image as input, the system further improves image quality by the following two steps: 1) noise reduction through an anisotropic diffusion technique that can effectively remove the background noise without destroying the vessels and translucency of the contents inside the vessels; and 2) enhance the blood vessel image there by suppressing the background through nonlinear normalization. This procedure is shown in FIG. 4 .
- the DSA image 214 is provided for look-up-table (LUT) based diffusion 402 or anisotropic diffusion.
- LUT look-up-table
- anisotropic diffusion is described in relation to equations 3-6.
- PDE partial differential equation
- ⁇ ⁇ t ⁇ I ⁇ ( x , y , t ) c ⁇ ( x , y , t ) ⁇ ⁇ ⁇ ⁇ I + ⁇ c ⁇ ( x , y , t ) ⁇ ⁇ I ( 4 )
- ⁇ is the Laplacian operator.
- the diffusion coefficient c(x,y,t) is the key in the smoothing process and it should encourage homogenous-region smoothing and inhibit the smoothing across the boundaries. It is chosen as a function of the magnitude of the gradient of the brightness function, i.e.
- K is the diffusion constant which controls the edge magnitude threshold.
- CUDA compute unified device architecture
- GPU Graphics Processing Unit
- a nonlinear normalization method 406 is used. The principle is to force the background to suppress the non-vascular structures and noise, while gradually enhancing the foreground vascular structures.
- the nonlinear normalization 406 is currently implemented as a LUT filter that is based on the pre-defined parameters and is set forth in relation to equation 7.
- I in (x,y) be the input DSA image (after diffusion)
- I out (X,Y) is mathematically given as:
- I out ⁇ ( x , y ) ⁇ I t ⁇ ( I i ⁇ ⁇ n ⁇ ( x , y ) I t ) y 1 , I i ⁇ ⁇ n ⁇ ( x , y ) ⁇ [ 0 , I t ] ( I t + 1 ) + ( I t + 1 ) ⁇ ( I i ⁇ ⁇ n ⁇ ( x , y ) - ( I t + 1 ) I t + 1 ) y 2 , I i ⁇ ⁇ n ⁇ ( x , y ) ⁇ [ I t + 1 ⁇ , ⁇ 255 ] ( 7 )
- I i is a pre-defined threshold, y 1 >1.0 and y 2 ⁇ 1.0 for class-based contrast enhancement. So if the intensity range of dye lies in lower half of the image and the background lies in the higher half of the image, the blood vessels are enhanced in limits. The result of this non-linear normalization is the enhanced DSA image 218 .
- the overall process flow for QCA is as follows. Using a DSA image, which may be an enhanced image as set forth above, a user or operator may select 502 one frame from the DSA to use in stenosis estimation. This beings the QCA process 504 . Once the frame is selected a user selects a region of interest 506 . This may entail zooming in 508 on the region of interest to provide an enlarged view of the region of interest 510 . This is illustrated in FIG. 6A which shows a frame having a blood vessel 600 that includes a lesion region 602 or stenosis.
- the operator can then choose “measure from centerline” (center line method 512 ), or “measure from sides” (vessel edge method 514 ). If “measure from centerline” is chosen the operator needs to select 516 several points 604 along the centerline of vessel 600 as illustrated in FIG. 6B . If “measure from sides”, is chosen the operator needs to select 516 several points 604 along the both sides of vessel 600 as illustrated in FIG. 6C .
- the QCA process 518 is then performed by the system which subsequently generates a final report 222 that indicates a percentage of stenosis. The final stenosis results is computed and displayed automatically.
- the QCA measurement process is set forth in FIG. 7 .
- the QCA measurement includes the following three sub-processes: (1) initial edge detection 720 ; (2) edge refinement 740 ; and (3) lesion measurement 760 .
- initial edge detection in QCA is performed by calculating the gradients of in a direction perpendicular to the centerline across the lumen and finding their peak values. See FIG. 8A .
- the contrast between the dyed vessel and edge is usually significant and easily detectable. That is, the vessel is often dark or nearly black in an image whereas the vessel wall and surrounding area is grey or white. This process is performed along the centerline 610 until edge points 612 are determined along both edges of the vessel 600 through the stenosis region 602 . See FIG. 8B
- FIG. 9 shows the process of initial edge detection. Initially, after the centerline points 604 are selected and connected to form the centerline 610 , perpendicular line construction 902 is performed to form lines perpendicular 904 to the centerline. On each perpendicular line, the weighted sum of the first and second gradients is calculated 906 for the subject target line. The peak value is used to find the resulting initial boundary or initial edge points 612 . These initial edge points are detected by using Eq. (8).
- f edge ⁇ ( x , y ) arg ⁇ ⁇ max f ⁇ ( x , y ) ⁇ ( w 1 ⁇ ⁇ f ⁇ ( x , y ) + w 2 ⁇ ⁇ 2 ⁇ f ⁇ ( x , y ) ) ( 8 )
- initial points along the two boundary lines those points are used as the initial edge points and above initial edge detection step is omitted.
- FIG. 10 This process uses an active contour model algorithm, which will be known to one skilled in the art, that deforms a contour to lock onto features of interest within in an image. Usually the features are lines, edges, and/or object boundaries. Given an approximation of the boundary of an object in an image (e.g., initial point 612 ), an active contour model can be used to find the “actual” boundary.
- An active contour is an ordered collection of n points in the image plane:
- V ⁇ v 1 , v 2 . . . v n ⁇
- E int (v i ) is an energy function dependent on the shape of the contour
- E ext (v i ) is an energy function dependent on the image properties, such as the gradient, near point v i .
- ⁇ and ⁇ are constants providing the relative weighting of the energy terms.
- the internal energy function is intended to enforce a shape on the deformable contour and to maintain a constant distance between the points in the contour. Additional terms can be added to influence the motion of the contour.
- the internal energy function used herein is defined as follows:
- E con (v i ) is the continuity energy that enforces the shape of the contour
- E cur (v i ) is a curvature energy that causes the contour to grow or shrink.
- c and b provide the relative weighting of the energy terms.
- the external energy function attracts the deformable contour to interesting features, such as object boundaries, in an image.
- image gradient is used.
- the image gradient should be large at the object boundary ( ⁇ , ⁇ 0). Therefore, the following external energy function is investigated:
- This process iteratively updates the initial edge points 904 until the energy function is minimized or a maximum number of iterations are achieved
- the result are refined edge points 910 as illustrated in FIG. 8C . That is, the points are smoothed and can be fitted with a curve.
- the edges of the lumen are identifies a dimensions need to be calculated for lesion measurement 740 . See FIG. 7 . If dealing with a lesion, it is usually desirable to identify at least the minimum luminal diameter (MLD), the percent stenosis and perhaps the length of the lesion.
- MLD minimum luminal diameter
- One method to determine lumen dimensions is set forth in the process flow sheet of FIG. 11 .
- a centerline is constructed 1104 between the two edges that represents the initial centerline 1106 of the lumen.
- Diameters are then measured by constructing lines that are perpendicular the centerline. Then, it is determined where each diameter line intersects 1108 each of the two edges. This defines a set of corresponding points on the boundary 1110 .
- the distance is calculated 1112 between the two points of intersection for each line.
- the distance having the shortest length is the minimum luminal diameter.
- the MLD may be divided by the normal diameter of the vessel to provide a stenosis percentage 1114 .
- the following data are calculated: (a) Minimum and maximum diameter; (b) Lesion length (the beginning and end of the lesion are determined automatically); (c) Stenosis reference (“normal”) diameter; and (d) Percent stenosis. In any case, this information may be output to the operator.
- the proposed system provides a number of advantages. For instance, may be used in the quantitative measurement of a blood vessel and the system includes three sub-systems: motion compensation, image enhancement, and stenosis measurement. The whole system results in a more intelligent and more accurate system for improving QCA measurement.
- Another advantage is that the proposed system can involve a registration step to align the mask and bolus images prior to subtracting procedure and QCA measurement. This motion compensation step can reduce the artifacts from misalignment (due to patient movement) and improve the accuracy of QCA measurement.
- Another advantage is that the proposed system can involve an image enhancement step prior to QCA measurement; the resulted DSA image is further enhanced by background diffusion and nonlinear normalization for better visualization. This image enhancement increases the contrast between the blood vessels and the background. This results in much improved contrast and very crisp subtraction images, in which the regions of interest are easily identifiable.
- the entire software architecture in the proposed system may be implemented on a GPU based framework, which makes visualization and computation faster by up to factor of 30.
- the entire scheme may therefore be implemented as a real-time scheme.
Abstract
Description
- This application claims priority under 35 U.S.C. §119 to U.S. Provisional Application No. 61/057,725 having a filing date of May 30, 2008, the entire contents of which are incorporated by reference herein.
- The present disclosure is directed to medical imaging systems. More specifically, the present disclosure is directed to systems and methods that alone or collectively facilitate real-time imaging.
- Coronary artery disease causes in excess of 1.5 million cases of myocardial infarction annually, and is the leading cause of death in the United States, resulting in more than 500,000 deaths per year. The accurate diagnosis and quantification of coronary artery disease is critical to subsequent treatment decisions. Despite the emergence of new techniques for the visualization and analysis of vascular structures, digital subtraction angiography (DSA) remains the preferred procedure to help clinicians to make clinical-decisions in all kinds of vascular diseases. Digital subtraction angiography (DSA) is a well-established modality for the visualization of blood vessels in the human body.
- During a DSA procedure a first series of frames are taken which represents the body anatomy in static form. They are called mask frames. Then dye is injected in the body through a catheter to visualize the blood vessel. The dye-injected frames are called the bolus frames. A subtraction procedure is preformed between the averaged bolus frame and the averaged mask frame to get DSA image. The patient is positioned on the X-ray imaging system while an X-ray movie is acquired and a DSA image is generated for an interventional cardiologist or a radiologist. Ideally, the DSA image mainly contains the dye-enhanced blood vessels and they appear dark. However, the artifacts due to patient motion and system noise frequently reduce the diagnostic value of the images.
- In the produced DSA image, a pathological aspect is associated with a vascular area where there is a significant deviation from the diameter of the healthy vascular. In particular, a stenosis is associated with a significant narrowing of the vascular and is quantified by parameters such as the percentage of stenosis. Stenosis limits blood flow by raising the resistance to flow through the vessel. In X-ray angiography (DSA image), the contrast agent ensures that the outlines of the blood flow are revealed on the X-ray to indicate any narrowing of the blood vessel. Direct visual examination of cine film coronary angiograms and manual estimation of the degree of vascular stenosis were complicated and subject to a large inter- and intra-observer variability.
- The present invention includes a system and method (i.e., utility) for determining the severity of a stenosis in a blood vessel. In one aspect, a method for improving DSA image quality includes: (1) registration of the mask and bolus images prior to subtracting procedure to reduce the artifacts from misalignment; (2) enhancement of the registered DSA image; and (3) detecting the boundary of blood vessel and quantitatively measuring percentage stenosis, which may be done automatically. This last step is sometimes referred to as Quantitative coronary angiography (QCA). Aspects of the present invention allow such QCA to be performed by a computer with minimal user input.
- The utility allows for the semi-automated quantitative measurement (QCA) of a lumen (e.g., blood vessel or artery). The utility involves a registration step of the mask and the bolus images prior to subtracting procedure and QCA measurement. Generally, such registration is a motion compensation step can reduce the artifacts from misalignment and improve the accuracy of QCA measurement. As will be appreciate, in coronary applications movement of the heart is nearly constant and motion compensation significantly improves overall image quality for QCA purposes.
- In one arrangement, the motion compensation is performed using an algorithm that involves an inverse-consistency constrain which implies that the correspondence provided by the registration in one direction matches closely with the correspondence in the opposite direction. This may entail a B-spline parameterization. Other image registration methods can also be applied in the proposed system to match the mask and bolus images.
- In one arrangement, the registration is performed hierarchically using a multi-resolution strategy in both, spatial domain and in domain of basis functions. The registration can be performed at ¼, ½ and full resolution using knot spacing of 8, 16 and 32. In addition to being faster, the multi-resolution strategy helps in improving the registration by matching global structures at lowest resolution and then matching local vessel structures as the resolution is refined. Due to the bi-directional approach, thin blood vessels edges are more prominent thereby avoiding the local minima both in iterative nature of thin vessel edge estimation and iterative nature of correction of bolus image. The subtraction process is applied on corrected bolus images.
- Image enhancement prior to QCA measurement also improves the overall QCA measurement. In one arrangement the registered DSA image is further enhanced by background diffusion to remove noise and nonlinear normalization for better visualization. This image enhancement increases the contrast between the blood vessels and the background. This results in much improved contrast and very crisp subtraction images, in which the regions of interest are easily identifiable. Enhancement of the registered DSA image can be performed using an anisotropic diffusion technique and a nonlinear normalization technique. Other similar image enhancement techniques can also be applied in the proposed utility to enhance DSA images before QCA measurement.
- The utility allows a user to input various types of lumen identification a center line method where a user identifies an approximate center of the lumen and an edge method where a user identifies initial edges of the lumen. Both methods can obtain percentage stenosis automatically with minimal user interaction. With the initial points selected by the user, the QCA measurement may include one or more of the following three sub-processes: (1) initial edge detection; (2) edge refinement; and (3) lesion measurement. Initial edge detection is performed where a user inputs centerline information. In any case, the utility may refine the edges of the lumen. In one arrangement, an active contour model algorithm may be used to refine the edges by deforming a contour to lock onto features of interest within in an image. Other edge detection or segmentation techniques can also be applied in the system for QCA measurement.
- For each stenosis or lesion, the following data are calculated: (a) minimum and maximum diameter; (b) lesion length (the beginning and end of the lesion are determined automatically); (c) Stenosis reference (“normal”) diameter; and (d) percent stenosis. Other measurements can also be included in the proposed system.
- To improve the overall speed of the utility, the entire software architecture may be implemented on a GPU based framework, which makes visualization and computation faster by up to factor of 30. Thus the entire utility may implemented substantially in real-time. Further, the utility may be used with various imaging modalities and may be used with 2-D or 3-D images.
-
FIG. 1 illustrates one embodiment of a X-ray DSA system. -
FIG. 2 illustrates a process flow diagram for performing an improved QCA measurement in accordance with various aspects of the present invention. -
FIG. 3 illustrates one non-limiting method for registering images. -
FIG. 4 illustrates one non-limiting method for enhancing an image. -
FIG. 5 illustrates a process flow diagram for performing a QCA measurement. -
FIGS. 6A-6C illustrate an image having a vein or artery with a lesion. -
FIG. 7 illustrates a process flow diagram of an edge detection, refinement and lesion measurement process. -
FIGS. 8A-8C illustrate the refinement of edge boundaries of the vein or artery ofFIGS. 6A-6C . -
FIG. 9 illustrates a process flow sheet for an edge detection process for detecting edge surfaces of an object in an image. -
FIG. 10 illustrates a process flow sheet for an edge refinement process. -
FIG. 11 illustrates a process flow sheet for a lesion measurement process. - Reference will now be made to the accompanying drawings, which assist in illustrating the various pertinent features of the various novel aspects of the present disclosure. Although the present invention will now be described primarily in conjunction with angiography utilizing X-ray imaging, it should be expressly understood that aspects of the present invention may be applicable to other medical imaging applications. For instance, angiography may be performed using a number of different medical imaging modalities, including biplane X-ray/DSA, magnetic resonance (MR), computed tomography (CT), ultrasound, and various combinations of these techniques. In this regard, the following description is presented for purposes of illustration and description. Furthermore, the description is not intended to limit the invention to the form disclosed herein. Consequently, variations and modifications commensurate with the following teachings, and skill and knowledge of the relevant art, are within the scope of the present invention. The embodiments described herein are further intended to explain known modes of practicing the invention and to enable others skilled in the art to utilize the invention in such, or other embodiments and with various modifications required by the particular application(s) or use(s) of the present invention.
-
FIG. 1 shows one exemplary setup for a real-time imaging procedure for use during a contrast media/dye injection procedure. As shown, a patient is positioned on anX-ray imaging system 100 and an X-ray movie is acquired by amovie acquisition system 102. An enhanced DSA image, as will be more fully discussed herein, is generated by anenhancement system 104 for output to adisplay 106 and/or stenosis processing system that is accessible to an interventional radiologist. - The projection images (e.g., CT images) are acquired at different time instants and consist of a movie with a series of frames before, during and after the dye injection. The series of frames include mask images that are free of contrast-enhancing dye in their field of
view 108 and bolus images that contain contrast-enhancing dye in their field ofview 108. That is, bolus frames are images that are acquired after injected dye has reached the field ofview 108. Themovie acquisition system 102 is operative to detect the frames before and after dye injection automatically to make feasible a real-time acquisition system. One approach for identifying frames before and after dye injection is to find intensity differences between successive frames, such that a large intensity difference is detected between the first frame after dye has reached the field of view (FOV) and the frame acquired before it. However, the patient may undergo some motion during the image acquisition causing such an intensity difference between even successive mask images. To avoid this, themovie acquisition system 102 may align successive frames together, such that the motion artifacts are minimized. The first image acquired after the dye has reached the FOV will therefore cause a high intensity difference with the previous frame not containing the dye in FOV. The subtraction image or ‘DSA image’ obtained by subtracting a mask frame from a bolus frame (or vice versa) will contain a near-zero value everywhere if both images belong to background. - Generally, the subtraction image or DSA image is obtained by computing a difference between pixel intensities of the mask image and the bolus image. The
enhancement system 104 may then enhance the contrast of the subtraction image. Such enhancement may include rescaling the intensities of the pixels in the subtraction image and/or the removal of noise from the subtraction image. - The acquisition and/or enhancement systems are computerized systems that may run application software and computer programs which can be used to control the system components, provide user interfaces, and/or provide features of the imaging system. The software may be originally provided on computer-readable media, such as compact disks (CDs), magnetic tape, or other mass storage medium. Alternatively, the software may be downloaded from electronic links such as a host or vendor website. The software is installed onto a hard drive and/or electronic memory of the system, and is accessed and controlled by the operating system. Software updates are also electronically available on mass storage media or downloadable from the host or vendor website. The software, as provided on the computer-readable media or downloaded from electronic links, represents a computer program product usable with a programmable computer processor having computer-readable program code embodied therein. The software contains one or more programming modules, subroutines, computer links, and compilations of executable code, which perform the functions of the imaging system. The user interacts with the software via keyboard, mouse, voice recognition, and other user-interface devices (e.g., user I/O devices) connected to the computer system.
-
FIG. 2 shows a method for improving DSA image quality and improving the accuracy of QCA measurements that includes three primary processes: The first process includes using an averaging system 202 to generate average mask frames 204 and average bolus frames 206 from an x-ray image or movie. These average images are motion compensated 208. Specifically, the averaged mask frame is registered with the averaged bolus frame to generate a registered mask frame 210, which minimizes the motion artifacts before subtraction. The registered mask frame 210 is then subtracted 212 from the average bolus frame to generate aDSA image 214. Various aspects of this process are set forth in co-pending U.S. application Ser. No. 11/609,743 entitled “Medical Image Enhancement System” the entire contents of which are incorporated herein by reference. The second process uses the registeredDSA image 214 as input. This is animage enhancement process 216 that uses an anisotropic diffusion technique to reduce noise, followed by a nonlinear normalization method to enhance the image of the blood vessel. The result of this is the production of anenhanced DSA image 218. The third process uses theenhanced DSA image 218 to do aQCA measurement 220 and generate a QCA report orvalue 222. The motion compensation process is discussed herein relation toFIG. 3 , the image enhancement process is discussed in relation toFIG. 4 and the QCA measurement is discussed in relation toFIGS. 5-11 . -
FIG. 3 illustrates the process of registering the averaged mask frame with averaged bolus frame, to minimize the motion artifacts. In this embodiment, this is done using an inverse consistent image registration using B-spline basis. In medical imaging, image registration is used to find a point-wise correspondence between a pair or group of similar anatomical objects. Image registration is often posed as an optimization problem that minimizes an objective function representing the difference between twoimages -
- where, I1(x) and I2(x) represent the intensity of image at location x, x represents the domain of the image. hi,j(x)=x+ui,j(x) represents the transformation that maps image Ii to image Ij in the Eulerian frame of reference and u(x) represents the displacement field. L is a differential operator and the second term in Eq. (1) represents an energy function. σ, ρ and χ are weights to adjust relative importance of the cost function.
- In Equation (1), the first term represents the symmetric squared intensity cost function and represents the integration of squared intensity difference between deformed reference image and the target image in both directions. The second term represents the energy regularization cost term and penalizes high derivatives of u(x). L is represented as a Laplacian operator mathematically given as: L=∇2. The third and last term represents the inverse consistency cost function, which penalizes differences between transformation in one direction and inverse of transformation in opposite direction. The
total cost 308 is computed 306 using Eq. 1 as a first step in registration. - The optimization problem posed in Eq. (2) is solved by using a B-spline parameterization as is known in the art. B-splines are used due to ease of computation, good approximation properties and their local support. It is also easier to incorporate landmarks in the cost term if we use spatial basis function. The above optimization problem is solved by solving for B-spline coefficients ci's, such that
-
- where, βi(x) represents the value of B-spline at location x, originating at index i. In this registration method, cubic B-splines are used. A gradient descent scheme is implemented based on the above parameterization. The total gradient cost is calculated and updated 310 with respect to the transformation parameters in every iteration. The transformation parameters are updated using the gradient descent update rule. Images are deformed into shape of one another using the updated correspondence and the cost function and gradient costs are calculated 314 until
convergence 316 when the frames are registered 318. The registration is performed hierarchically using a multi-resolution strategy in both, spatial domain and in domain of basis functions. The registration is performed at ¼, ½ and full resolution using knot spacing of 8, 16 and 32. In addition to being faster, the multi-resolution strategy helps in improving the registration by matching global structures at lowest resolution and their matching local structures as the resolution is refined. - Referring again to
FIG. 2 , these registered frames are subtracted to provide the DSA image for image enhancement. That is, after image registration, subtraction is performed between the registered mask image and bolus image, resulting in the motion compensated (registered)DSA image 214. Using the registered DSA image as input, the system further improves image quality by the following two steps: 1) noise reduction through an anisotropic diffusion technique that can effectively remove the background noise without destroying the vessels and translucency of the contents inside the vessels; and 2) enhance the blood vessel image there by suppressing the background through nonlinear normalization. This procedure is shown inFIG. 4 . - Initially, the
DSA image 214 is provided for look-up-table (LUT) baseddiffusion 402 or anisotropic diffusion.Such diffusion 402 is described in relation to equations 3-6. Such an anisotropic diffusion and is based on partial differential equation (PDE) for noise smoothing. Given an image l(x,y,t) at time scale t, the diffusion process is expressed as: -
- where ∇ is the gradient operator, div is the divergence operator, and c(x,y,t) is the diffusion coefficient at location (x,y) at time t, With applying the divergence operator, Eq. (3) can be rewritten as
-
- where Δ is the Laplacian operator. The diffusion coefficient c(x,y,t) is the key in the smoothing process and it should encourage homogenous-region smoothing and inhibit the smoothing across the boundaries. It is chosen as a function of the magnitude of the gradient of the brightness function, i.e.
-
c(x,y,t)=g(∥∇I(x,y,t)∥) (5) - The suggested functions for g(·) are the following two
-
- where K is the diffusion constant which controls the edge magnitude threshold. Generally speaking, a larger K produces a smoother result in a homogenous region than a smaller one. Here we apply diffusion technique on the input DSA images to smooth background thereby reducing the structured noise. Also, the algorithm may be implemented using compute unified device architecture (CUDA), which is developed by nVidia to run general purpose computations on a Graphics Processing Unit (GPU). The result of the
diffusion process 402 is a diffusedimage 404. - To increase the contrast between the blood vessels of and background in the DSA image, a nonlinear normalization method 406 is used. The principle is to force the background to suppress the non-vascular structures and noise, while gradually enhancing the foreground vascular structures. The nonlinear normalization 406 is currently implemented as a LUT filter that is based on the pre-defined parameters and is set forth in relation to
equation 7. - Letting Iin(x,y) be the input DSA image (after diffusion), the nonlinear normalization, defined as Iout(X,Y) is mathematically given as:
-
- Here Ii is a pre-defined threshold, y1>1.0 and y2<1.0 for class-based contrast enhancement. So if the intensity range of dye lies in lower half of the image and the background lies in the higher half of the image, the blood vessels are enhanced in limits. The result of this non-linear normalization is the
enhanced DSA image 218. - After these processes, the DSA image quality has been improved greatly. Such improvement in the image facilitates QCA measurement as set forth herein. However, it will be appreciated that the QCA measurement set forth herein is considered novel in and of itself.
- Set forth herein are two semi-automatic algorithms to compute percentage stenosis with minimum user intervention. As shown in
FIG. 5 , the overall process flow for QCA is as follows. Using a DSA image, which may be an enhanced image as set forth above, a user or operator may select 502 one frame from the DSA to use in stenosis estimation. This beings theQCA process 504. Once the frame is selected a user selects a region ofinterest 506. This may entail zooming in 508 on the region of interest to provide an enlarged view of the region ofinterest 510. This is illustrated inFIG. 6A which shows a frame having ablood vessel 600 that includes alesion region 602 or stenosis. - The operator can then choose “measure from centerline” (center line method 512), or “measure from sides” (vessel edge method 514). If “measure from centerline” is chosen the operator needs to select 516
several points 604 along the centerline ofvessel 600 as illustrated inFIG. 6B . If “measure from sides”, is chosen the operator needs to select 516several points 604 along the both sides ofvessel 600 as illustrated inFIG. 6C . TheQCA process 518 is then performed by the system which subsequently generates afinal report 222 that indicates a percentage of stenosis. The final stenosis results is computed and displayed automatically. - The QCA measurement process is set forth in
FIG. 7 . With the initial points selected by the user, the QCA measurement includes the following three sub-processes: (1)initial edge detection 720; (2)edge refinement 740; and (3) lesion measurement 760. - When the user chooses
initial points 604 along thecenterline 610, initial edge detection in QCA is performed by calculating the gradients of in a direction perpendicular to the centerline across the lumen and finding their peak values. SeeFIG. 8A . As will be appreciated, the contrast between the dyed vessel and edge is usually significant and easily detectable. That is, the vessel is often dark or nearly black in an image whereas the vessel wall and surrounding area is grey or white. This process is performed along thecenterline 610 until edge points 612 are determined along both edges of thevessel 600 through thestenosis region 602. SeeFIG. 8B -
FIG. 9 shows the process of initial edge detection. Initially, after the centerline points 604 are selected and connected to form thecenterline 610,perpendicular line construction 902 is performed to form lines perpendicular 904 to the centerline. On each perpendicular line, the weighted sum of the first and second gradients is calculated 906 for the subject target line. The peak value is used to find the resulting initial boundary or initial edge points 612. These initial edge points are detected by using Eq. (8). -
- If the user chooses initial points along the two boundary lines, those points are used as the initial edge points and above initial edge detection step is omitted.
- After initial edge point detection, a group of edge points are identified along two boundaries of the vessel. See
FIG. 8B . These initial points are refined 740 (seeFIG. 7 ) to get a smooth boundary curve. This process is set forth inFIG. 10 . This process uses an active contour model algorithm, which will be known to one skilled in the art, that deforms a contour to lock onto features of interest within in an image. Usually the features are lines, edges, and/or object boundaries. Given an approximation of the boundary of an object in an image (e.g., initial point 612), an active contour model can be used to find the “actual” boundary. An active contour is an ordered collection of n points in the image plane: -
V={v1, v2 . . . vn} -
v i=(x i ,y i), i=1, . . . n (9) - The points in the contour iteratively approach the boundary of an object through the solution of an energy minimization problem. For each point in the neighborhood of vi, an
energy term 902 is computed: -
E i =αE int(v i)+βE ext(v i) (10) - where Eint(vi) is an energy function dependent on the shape of the contour and Eext(vi) is an energy function dependent on the image properties, such as the gradient, near point vi. α and β are constants providing the relative weighting of the energy terms.
- The internal energy function is intended to enforce a shape on the deformable contour and to maintain a constant distance between the points in the contour. Additional terms can be added to influence the motion of the contour. The internal energy function used herein is defined as follows:
-
αE int(v i)=bE con(v i)+cE cur(v i) (11) - where Econ(vi) is the continuity energy that enforces the shape of the contour and Ecur(vi) is a curvature energy that causes the contour to grow or shrink. c and b provide the relative weighting of the energy terms.
- The external energy function attracts the deformable contour to interesting features, such as object boundaries, in an image. Here image gradient is used. The image gradient should be large at the object boundary (β, <0). Therefore, the following external energy function is investigated:
-
βE int(v i)=βE grad(v i) (12) - In summary, energy function at each vi is minimized:
-
E i =bE con(v i)+cE cur(v i)+βE grad(v i) (13) - In (8), b>0,c>0 and β<0.
- This process iteratively updates the initial edge points 904 until the energy function is minimized or a maximum number of iterations are achieved The result are refined edge points 910 as illustrated in
FIG. 8C . That is, the points are smoothed and can be fitted with a curve. - After these first two
processes lesion measurement 740. SeeFIG. 7 . If dealing with a lesion, it is usually desirable to identify at least the minimum luminal diameter (MLD), the percent stenosis and perhaps the length of the lesion. - One method to determine lumen dimensions is set forth in the process flow sheet of
FIG. 11 . Using the detectedboundaries 1102 as set forth by the refined edge points, a centerline is constructed 1104 between the two edges that represents theinitial centerline 1106 of the lumen. Diameters are then measured by constructing lines that are perpendicular the centerline. Then, it is determined where each diameter line intersects 1108 each of the two edges. This defines a set of corresponding points on theboundary 1110. The distance is calculated 1112 between the two points of intersection for each line. The distance having the shortest length is the minimum luminal diameter. The MLD may be divided by the normal diameter of the vessel to provide astenosis percentage 1114. Generally, for each lesion, the following data are calculated: (a) Minimum and maximum diameter; (b) Lesion length (the beginning and end of the lesion are determined automatically); (c) Stenosis reference (“normal”) diameter; and (d) Percent stenosis. In any case, this information may be output to the operator. - The proposed system provides a number of advantages. For instance, may be used in the quantitative measurement of a blood vessel and the system includes three sub-systems: motion compensation, image enhancement, and stenosis measurement. The whole system results in a more intelligent and more accurate system for improving QCA measurement. Another advantage is that the proposed system can involve a registration step to align the mask and bolus images prior to subtracting procedure and QCA measurement. This motion compensation step can reduce the artifacts from misalignment (due to patient movement) and improve the accuracy of QCA measurement. Another advantage is that the proposed system can involve an image enhancement step prior to QCA measurement; the resulted DSA image is further enhanced by background diffusion and nonlinear normalization for better visualization. This image enhancement increases the contrast between the blood vessels and the background. This results in much improved contrast and very crisp subtraction images, in which the regions of interest are easily identifiable.
- It will be appreciated that the entire software architecture in the proposed system may be implemented on a GPU based framework, which makes visualization and computation faster by up to factor of 30. The entire scheme may therefore be implemented as a real-time scheme.
- The foregoing description of the present invention has been presented for purposes of illustration and description. Furthermore, the description is not intended to limit the invention to the form disclosed herein. Consequently, variations and modifications commensurate with the above teachings, and skill and knowledge of the relevant art, are within the scope of the present invention. The embodiments described hereinabove are further intended to explain best modes known of practicing the invention and to enable others skilled in the art to utilize the invention in such, or other embodiments and with various modifications required by the particular application(s) or use(s) of the present invention. It is intended that the appended claims be construed to include alternative embodiments to the extent permitted by the prior art.
Claims (19)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/475,995 US20100004526A1 (en) | 2008-06-04 | 2009-06-01 | Abnormality finding in projection images |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US5872508P | 2008-06-04 | 2008-06-04 | |
US12/475,995 US20100004526A1 (en) | 2008-06-04 | 2009-06-01 | Abnormality finding in projection images |
Publications (1)
Publication Number | Publication Date |
---|---|
US20100004526A1 true US20100004526A1 (en) | 2010-01-07 |
Family
ID=41397674
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/475,995 Abandoned US20100004526A1 (en) | 2008-06-04 | 2009-06-01 | Abnormality finding in projection images |
US12/995,004 Active 2030-01-09 US8552011B2 (en) | 2008-06-04 | 2009-06-04 | 6-(C2-6alkylseleno)purines and methods for treating neural measles viral infection therewith |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/995,004 Active 2030-01-09 US8552011B2 (en) | 2008-06-04 | 2009-06-04 | 6-(C2-6alkylseleno)purines and methods for treating neural measles viral infection therewith |
Country Status (5)
Country | Link |
---|---|
US (2) | US20100004526A1 (en) |
EP (1) | EP2310390A4 (en) |
CN (1) | CN102112477A (en) |
CA (1) | CA2726587C (en) |
WO (1) | WO2009146542A1 (en) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110224542A1 (en) * | 2010-03-12 | 2011-09-15 | Sushil Mittal | Method and System for Automatic Detection and Classification of Coronary Stenoses in Cardiac CT Volumes |
US20110286645A1 (en) * | 2009-01-29 | 2011-11-24 | Koninklijke Philips Electronics N.V. | Transmural perfusion gradient image analysis |
US20130066197A1 (en) * | 2011-09-13 | 2013-03-14 | Celine Pruvot | System and method for blood vessel stenosis visualization and navigation |
CN103135853A (en) * | 2011-11-21 | 2013-06-05 | 原相科技股份有限公司 | Optical input device, detection method and method for optical input device |
US8526700B2 (en) | 2010-10-06 | 2013-09-03 | Robert E. Isaacs | Imaging system and method for surgical and interventional medical procedures |
US20140044332A1 (en) * | 2012-08-10 | 2014-02-13 | National Taiwan University | Transformation method for diffusion spectrum imaging using large deformation diffeomorphic metric mapping |
KR101401734B1 (en) | 2010-09-17 | 2014-05-30 | 한국전자통신연구원 | Apparatus and method for detecting examination tissue using ultrasound signal |
US9767557B1 (en) | 2016-06-23 | 2017-09-19 | Siemens Healthcare Gmbh | Method and system for vascular disease detection using recurrent neural networks |
US9785246B2 (en) | 2010-10-06 | 2017-10-10 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US20180279983A1 (en) * | 2017-03-28 | 2018-10-04 | Canon Medical Systems Corporation | Medical image processing apparatus, medical image processing method, and x-ray diagnostic apparatus |
US10568587B2 (en) * | 2017-03-28 | 2020-02-25 | Canon Medical Systems Corporation | X-ray diagnostic apparatus, image processing apparatus, and image processing method |
CN111477303A (en) * | 2020-04-15 | 2020-07-31 | 吉林大学第一医院 | Aligning 3D real-time display system based on digital subtraction technology |
US10762637B2 (en) | 2017-10-27 | 2020-09-01 | Siemens Healthcare Gmbh | Vascular segmentation using fully convolutional and recurrent neural networks |
US10977792B2 (en) * | 2017-04-05 | 2021-04-13 | Siemens Healthcare Gmbh | Quantitative evaluation of time-varying data |
US11231787B2 (en) | 2010-10-06 | 2022-01-25 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
Citations (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4692864A (en) * | 1985-05-23 | 1987-09-08 | Elscint Ltd. | Method of determining stenosis of blood vessels |
US4804250A (en) * | 1987-06-12 | 1989-02-14 | The United States Of America As Represented By The Secretary Of The Army | Optical processor for an adaptive pattern classifier |
US5151856A (en) * | 1989-08-30 | 1992-09-29 | Technion R & D Found. Ltd. | Method of displaying coronary function |
US5908389A (en) * | 1996-09-27 | 1999-06-01 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging of harmonic frequencies with speckle reduction processing |
US6377832B1 (en) * | 1998-03-20 | 2002-04-23 | Georgia Tech Research Corporation | System and method for analyzing a medical image |
US6430432B1 (en) * | 1998-12-21 | 2002-08-06 | Koninklijke Philips Electronics, N.V. | Deriving dimensions of a detail of an object |
US20040223636A1 (en) * | 1999-11-19 | 2004-11-11 | Edic Peter Michael | Feature quantification from multidimensional image data |
US20040243365A1 (en) * | 2003-03-20 | 2004-12-02 | University Of Washington | Computation of wall thickness |
US6889072B2 (en) * | 1993-06-07 | 2005-05-03 | Martin R. Prince | Method and apparatus for administration of contrast agents for use in magnetic resonance arteriography |
US6990368B2 (en) * | 2002-04-04 | 2006-01-24 | Surgical Navigation Technologies, Inc. | Method and apparatus for virtual digital subtraction angiography |
US6999811B2 (en) * | 2001-07-25 | 2006-02-14 | Koninklijke Philips Electronics N.V. | Method and device for the registration of two 3D image data sets |
US7012603B2 (en) * | 2001-11-21 | 2006-03-14 | Viatronix Incorporated | Motion artifact detection and correction |
US20060056695A1 (en) * | 2004-09-10 | 2006-03-16 | Min Wu | Method for concealing data in curves of an image |
US20060072823A1 (en) * | 2004-10-04 | 2006-04-06 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
US20060106428A1 (en) * | 2004-11-18 | 2006-05-18 | Cardiac Pacemakers, Inc. | Cardiac rhythm management device with neural sensor |
US20060182349A1 (en) * | 2005-01-12 | 2006-08-17 | Valadez Gerardo H | System and method for quantifying motion artifacts in perfusion image sequences |
US20070260143A1 (en) * | 2003-06-23 | 2007-11-08 | Siemens Medical Solutions Usa, Inc. | Method and system for simultaneously displaying relationships of measurements of features associated with a medical image |
US20080051648A1 (en) * | 2006-08-25 | 2008-02-28 | Suri Jasjit S | Medical image enhancement system |
US20080095422A1 (en) * | 2006-10-18 | 2008-04-24 | Suri Jasjit S | Alignment method for registering medical images |
US7545967B1 (en) * | 2002-09-18 | 2009-06-09 | Cornell Research Foundation Inc. | System and method for generating composite subtraction images for magnetic resonance imaging |
US7826884B2 (en) * | 2006-07-31 | 2010-11-02 | Siemens Medical Solutions Usa, Inc. | Live fluoroscopic roadmapping including targeted automatic pixel shift for misregistration correction |
US7949170B2 (en) * | 2003-11-12 | 2011-05-24 | Hitachi Medical Corporation | Image processing method, image processing device, computer aided detection, and method for filtering along the time axis |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CZ126799A3 (en) | 1996-10-16 | 1999-07-14 | Icn Pharmaceuticals | Purine l-nucleosides and their analogs as well as pharmaceutical composition containing thereof |
-
2009
- 2009-06-01 US US12/475,995 patent/US20100004526A1/en not_active Abandoned
- 2009-06-04 CA CA2726587A patent/CA2726587C/en active Active
- 2009-06-04 US US12/995,004 patent/US8552011B2/en active Active
- 2009-06-04 WO PCT/CA2009/000770 patent/WO2009146542A1/en active Application Filing
- 2009-06-04 EP EP09757004A patent/EP2310390A4/en not_active Withdrawn
- 2009-06-04 CN CN2009801307160A patent/CN102112477A/en active Pending
Patent Citations (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4692864A (en) * | 1985-05-23 | 1987-09-08 | Elscint Ltd. | Method of determining stenosis of blood vessels |
US4804250A (en) * | 1987-06-12 | 1989-02-14 | The United States Of America As Represented By The Secretary Of The Army | Optical processor for an adaptive pattern classifier |
US5151856A (en) * | 1989-08-30 | 1992-09-29 | Technion R & D Found. Ltd. | Method of displaying coronary function |
US6889072B2 (en) * | 1993-06-07 | 2005-05-03 | Martin R. Prince | Method and apparatus for administration of contrast agents for use in magnetic resonance arteriography |
US5908389A (en) * | 1996-09-27 | 1999-06-01 | Atl Ultrasound, Inc. | Ultrasonic diagnostic imaging of harmonic frequencies with speckle reduction processing |
US6377832B1 (en) * | 1998-03-20 | 2002-04-23 | Georgia Tech Research Corporation | System and method for analyzing a medical image |
US6430432B1 (en) * | 1998-12-21 | 2002-08-06 | Koninklijke Philips Electronics, N.V. | Deriving dimensions of a detail of an object |
US20040223636A1 (en) * | 1999-11-19 | 2004-11-11 | Edic Peter Michael | Feature quantification from multidimensional image data |
US6999811B2 (en) * | 2001-07-25 | 2006-02-14 | Koninklijke Philips Electronics N.V. | Method and device for the registration of two 3D image data sets |
US7012603B2 (en) * | 2001-11-21 | 2006-03-14 | Viatronix Incorporated | Motion artifact detection and correction |
US6990368B2 (en) * | 2002-04-04 | 2006-01-24 | Surgical Navigation Technologies, Inc. | Method and apparatus for virtual digital subtraction angiography |
US7545967B1 (en) * | 2002-09-18 | 2009-06-09 | Cornell Research Foundation Inc. | System and method for generating composite subtraction images for magnetic resonance imaging |
US20040243365A1 (en) * | 2003-03-20 | 2004-12-02 | University Of Washington | Computation of wall thickness |
US7353117B2 (en) * | 2003-03-20 | 2008-04-01 | University Of Washington | Computation of wall thickness |
US20070260143A1 (en) * | 2003-06-23 | 2007-11-08 | Siemens Medical Solutions Usa, Inc. | Method and system for simultaneously displaying relationships of measurements of features associated with a medical image |
US7949170B2 (en) * | 2003-11-12 | 2011-05-24 | Hitachi Medical Corporation | Image processing method, image processing device, computer aided detection, and method for filtering along the time axis |
US20060056695A1 (en) * | 2004-09-10 | 2006-03-16 | Min Wu | Method for concealing data in curves of an image |
US20060072823A1 (en) * | 2004-10-04 | 2006-04-06 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
US20060106428A1 (en) * | 2004-11-18 | 2006-05-18 | Cardiac Pacemakers, Inc. | Cardiac rhythm management device with neural sensor |
US20060182349A1 (en) * | 2005-01-12 | 2006-08-17 | Valadez Gerardo H | System and method for quantifying motion artifacts in perfusion image sequences |
US7826884B2 (en) * | 2006-07-31 | 2010-11-02 | Siemens Medical Solutions Usa, Inc. | Live fluoroscopic roadmapping including targeted automatic pixel shift for misregistration correction |
US20080051648A1 (en) * | 2006-08-25 | 2008-02-28 | Suri Jasjit S | Medical image enhancement system |
US20080095422A1 (en) * | 2006-10-18 | 2008-04-24 | Suri Jasjit S | Alignment method for registering medical images |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110286645A1 (en) * | 2009-01-29 | 2011-11-24 | Koninklijke Philips Electronics N.V. | Transmural perfusion gradient image analysis |
US8755575B2 (en) * | 2009-01-29 | 2014-06-17 | Koninklijke Philips N.V. | Transmural perfusion gradient image analysis |
US8526699B2 (en) | 2010-03-12 | 2013-09-03 | Siemens Aktiengesellschaft | Method and system for automatic detection and classification of coronary stenoses in cardiac CT volumes |
US20110224542A1 (en) * | 2010-03-12 | 2011-09-15 | Sushil Mittal | Method and System for Automatic Detection and Classification of Coronary Stenoses in Cardiac CT Volumes |
KR101401734B1 (en) | 2010-09-17 | 2014-05-30 | 한국전자통신연구원 | Apparatus and method for detecting examination tissue using ultrasound signal |
US10444855B2 (en) | 2010-10-06 | 2019-10-15 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US10139920B2 (en) | 2010-10-06 | 2018-11-27 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US10684697B2 (en) | 2010-10-06 | 2020-06-16 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US8792704B2 (en) | 2010-10-06 | 2014-07-29 | Saferay Spine Llc | Imaging system and method for use in surgical and interventional medical procedures |
US11231787B2 (en) | 2010-10-06 | 2022-01-25 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US11941179B2 (en) | 2010-10-06 | 2024-03-26 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US9785246B2 (en) | 2010-10-06 | 2017-10-10 | Nuvasive, Inc. | Imaging system and method for use in surgical and interventional medical procedures |
US8526700B2 (en) | 2010-10-06 | 2013-09-03 | Robert E. Isaacs | Imaging system and method for surgical and interventional medical procedures |
CN103054598A (en) * | 2011-09-13 | 2013-04-24 | 通用电气公司 | System and method for blood vessel stenosis visualization and navigation |
US20130066197A1 (en) * | 2011-09-13 | 2013-03-14 | Celine Pruvot | System and method for blood vessel stenosis visualization and navigation |
US20130194236A1 (en) * | 2011-11-21 | 2013-08-01 | Pixart Imaging Inc. | Optical input device and input detection thereof |
CN103135853A (en) * | 2011-11-21 | 2013-06-05 | 原相科技股份有限公司 | Optical input device, detection method and method for optical input device |
US9047695B2 (en) * | 2012-08-10 | 2015-06-02 | National Taiwan University | Transformation method for diffusion spectrum imaging using large deformation diffeomorphic metric mapping |
US20140044332A1 (en) * | 2012-08-10 | 2014-02-13 | National Taiwan University | Transformation method for diffusion spectrum imaging using large deformation diffeomorphic metric mapping |
US9881372B2 (en) | 2016-06-23 | 2018-01-30 | Siemens Healthcare Gmbh | Method and system for vascular disease detection using recurrent neural networks |
US9767557B1 (en) | 2016-06-23 | 2017-09-19 | Siemens Healthcare Gmbh | Method and system for vascular disease detection using recurrent neural networks |
US20180279983A1 (en) * | 2017-03-28 | 2018-10-04 | Canon Medical Systems Corporation | Medical image processing apparatus, medical image processing method, and x-ray diagnostic apparatus |
US11793478B2 (en) * | 2017-03-28 | 2023-10-24 | Canon Medical Systems Corporation | Medical image processing apparatus, medical image processing method, and X-ray diagnostic apparatus |
US10568587B2 (en) * | 2017-03-28 | 2020-02-25 | Canon Medical Systems Corporation | X-ray diagnostic apparatus, image processing apparatus, and image processing method |
US10977792B2 (en) * | 2017-04-05 | 2021-04-13 | Siemens Healthcare Gmbh | Quantitative evaluation of time-varying data |
US10762637B2 (en) | 2017-10-27 | 2020-09-01 | Siemens Healthcare Gmbh | Vascular segmentation using fully convolutional and recurrent neural networks |
CN111477303A (en) * | 2020-04-15 | 2020-07-31 | 吉林大学第一医院 | Aligning 3D real-time display system based on digital subtraction technology |
Also Published As
Publication number | Publication date |
---|---|
WO2009146542A8 (en) | 2010-03-11 |
US20110086867A1 (en) | 2011-04-14 |
EP2310390A4 (en) | 2012-03-21 |
WO2009146542A1 (en) | 2009-12-10 |
CA2726587C (en) | 2017-09-19 |
US8552011B2 (en) | 2013-10-08 |
CN102112477A (en) | 2011-06-29 |
CA2726587A1 (en) | 2009-12-10 |
EP2310390A1 (en) | 2011-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20100004526A1 (en) | Abnormality finding in projection images | |
US20220265239A1 (en) | Systems and methods for assessing the severity of plaque and/or stenotic lesions using contrast distribution predictions and measurements | |
US20220284572A1 (en) | Systems, methods, and devices for medical image analysis, diagnosis, risk stratification, decision making and/or disease tracking | |
Jerman et al. | Enhancement of vascular structures in 3D and 2D angiographic images | |
US9754390B2 (en) | Reconstruction of time-varying data | |
Fung et al. | Cerebral blood flow with [15O] water PET studies using an image-derived input function and MR-defined carotid centerlines | |
Wan et al. | Automated coronary artery tree segmentation in X-ray angiography using improved Hessian based enhancement and statistical region merging | |
Frangi et al. | Quantitative analysis of vascular morphology from 3D MR angiograms: in vitro and in vivo results | |
Galassi et al. | 3D reconstruction of coronary arteries from 2D angiographic projections using non-uniform rational basis splines (NURBS) for accurate modelling of coronary stenoses | |
Athanasiou et al. | Three-dimensional reconstruction of coronary arteries and plaque morphology using CT angiography–comparison and registration with IVUS | |
US11127138B2 (en) | Automatic detection and quantification of the aorta from medical images | |
Wang et al. | Segmentation of lumen and outer wall of abdominal aortic aneurysms from 3D black-blood MRI with a registration based geodesic active contour model | |
Caldairou et al. | A surface patch-based segmentation method for hippocampal subfields | |
Gupta et al. | Cardiac MR perfusion image processing techniques: a survey | |
US8577104B2 (en) | Liver lesion segmentation | |
US9569839B2 (en) | Image processing apparatus, method and medical image device | |
Manniesing et al. | Robust CTA lumen segmentation of the atherosclerotic carotid artery bifurcation in a large patient population | |
Xu et al. | An improved algorithm for vessel centerline tracking in coronary angiograms | |
Lalys et al. | Generic thrombus segmentation from pre-and post-operative CTA | |
Park | Connectivity-based local adaptive thresholding for carotid artery segmentation using MRA images | |
Sakellarios et al. | Novel methodology for 3D reconstruction of carotid arteries and plaque characterization based upon magnetic resonance imaging carotid angiography data | |
He et al. | 3-D B-spline wavelet-based local standard deviation (BWLSD): its application to edge detection and vascular segmentation in magnetic resonance angiography | |
Hameeteman et al. | Carotid wall volume quantification from magnetic resonance images using deformable model fitting and learning-based correction of systematic errors | |
CN112819800A (en) | DSA image recognition method, device and storage medium | |
CN112784928A (en) | DSA image recognition method, device and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: EIGEN, INC., CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:WEI, LIYANG;SURI, JASJIT S.;REEL/FRAME:022774/0042;SIGNING DATES FROM 20090601 TO 20090603 |
|
AS | Assignment |
Owner name: KAZI MANAGEMENT VI, LLC, VIRGIN ISLANDS, U.S. Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:EIGEN, INC.;REEL/FRAME:024652/0493 Effective date: 20100630 |
|
AS | Assignment |
Owner name: KAZI, ZUBAIR, VIRGIN ISLANDS, U.S. Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:KAZI MANAGEMENT VI, LLC;REEL/FRAME:024929/0310 Effective date: 20100630 |
|
AS | Assignment |
Owner name: KAZI MANAGEMENT ST. CROIX, LLC, VIRGIN ISLANDS, U. Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:KAZI, ZUBAIR;REEL/FRAME:025013/0245 Effective date: 20100630 |
|
AS | Assignment |
Owner name: IGT, LLC, CALIFORNIA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:KAZI MANAGEMENT ST. CROIX, LLC;REEL/FRAME:025132/0199 Effective date: 20100630 |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |