Some applications of the present invention generally relate to medical imaging. Specifically, some applications of the present invention relate to medical imaging and analysis of such images when such images are acquired in the presence of a tool in the subject's body.
Vascular catheterizations, such as coronary catheterizations, are frequently-performed medical interventions. Such interventions are typically performed in order to diagnose the blood vessels for potential disease, and/or to treat diseased blood vessels. Typically, in order to enable observation of blood vessels, the catheterization is performed under extraluminal imaging. Additionally, for some procedures, an endoluminal data-acquisition device is used to perform endoluminal imaging and/or measurements. The extraluminal imaging and, where applicable, the endoluminal data are typically evaluated by the medical staff in combination with one another in the course of the intervention, as well as post procedurally.
In accordance with some applications of the present invention, a device is detected within a subject's body, and an indication of a device-specific parameter is generated on an output device, in response to detecting the device. Typically, the presence of a device within at least a portion of an image of a blood vessel, and the classification of the device as a given type of device are determined by a computer processor. Based upon the classification of the device as the given type of device, a parameter to be calculated is designated by the computer processor. In response thereto, the designated parameter is calculated automatically by the computer processor, and an output is generated on an output device (e.g., a display) in response to the calculated parameter. For example, in response to determining that a balloon is present inside the blood vessel, the minimum diameter of the blood vessel within a portion of the blood vessel in which the balloon is present may be designated as the parameter to be calculated. For some applications, the parameter that is designated to be calculated is the minimum diameter of the blood vessel within the portion of the blood vessel in which the balloon is present, at the occurrence of maximal inflation of the balloon within the blood vessel.
For some applications of the present invention, the location of a tip of a radiopaque portion of a wire within an image of a blood vessel is automatically determined by a computer processor. Typically, for each of the pixels within at least a portion of the image, a wireness measure (i.e., a measure of the extent to which the pixels has a wire-like characteristic) is determined by the computer processor, by measuring an extent to which the pixel, together with other pixels within the image, forms part of a long, thin set of pixels having a given characteristic, such as darkness, brightness, and/or a different characteristic (i.e., a wireness-indicating characteristic). In addition, for each of the pixels within at least the portion of the image the pixel intensity is determined by the computer processor. Subsequently, it is determined by the computer processor whether there is at least one pixel at which there is a change in the wireness measure by more than a threshold amount, relative to at least some of the pixels that belong to the set of pixels that have the wireness-indicating characteristic. In response to determining that the change in the wireness measure of at least one pixel within a given sampling region does exceed the threshold, the computer processor then determines whether at the at least one pixel there is a change, by more than a threshold amount, in the intensity of the pixel relative to the value of the intensity of at least some of the set of pixels that have the wireness-indicating characteristic. In response to determining that the change in intensity at the at least one pixel within the given sampling region does exceed the threshold, it is determined by the computer processor that the tip of the radiopaque portion of the wire is disposed within the given sampling region.
For some applications, within the designated sampling regions, the local directionality of the set of pixels that have the wireness-indicating characteristic is determined, and the changes in the wireness measure and/or intensity are measured along that direction. For some applications, it is first determined whether at at least one pixel within the sampling regions there is a change in intensity that exceeds the threshold, and, in response to determining that at at least one pixel within a given sampling region there is a change in intensity that exceeds the threshold, it is determined whether at the at least one pixel there is a change in the wireness measure that exceeds the threshold.
For some applications, in response to determining the location of the tip of the radiopaque portion of the wire within a given image, the image is aligned with a second image, by aligning the radiopaque portions of the guidewires in each of the images with one another. Alternatively or additionally, identification of the location of the tip of the radiopaque portion of the wire within a given image is used to facilitate the determination of the location of an endoluminal device within the lumen, for example, in accordance with techniques described in US 2012/0004537 to Tolkowsky, and/or WO 13/174472 to Steinberg, which are incorporated herein by reference. Further alternatively or additionally, identification of the location of the tip of the radiopaque portion of the wire within a given image is used to facilitate the determination of a transformation function for mapping between the image and a second image, for example, in accordance with techniques described in WO 13/174472 to Steinberg, which is incorporated herein by reference.
There is therefore provided, in accordance with some applications of the present invention, a method for use with an image of at least one blood vessel of a subject, including:
using at least one computer processor:
For some applications, determining the presence of the device within the portion of the blood vessel within the image includes receiving an input from a user that is indicative of the presence of the device within the portion of the blood vessel within the image.
For some applications, determining the classification of the device includes receiving an input from a user that is indicative of the device being the given type of device.
For some applications, determining the presence of the device within the portion of the blood vessel within the image includes automatically determining the presence of the device within the portion of the blood vessel within the image, by analyzing the image using the computer processor.
For some applications, determining the classification of the device includes automatically determining the classification of the device by analyzing the image using the computer processor.
For some applications, designating the parameter to be calculated includes, based upon the classification of the device as the given type of device, designating an event and designating a parameter to be calculated at the occurrence of the event.
For some applications, designating the parameter to be calculated includes, based upon the classification of the device as the given type of device, designating a dimension of the blood vessel to be calculated.
For some applications, designating the parameter to be calculated includes, based upon the classification of the device as the given type of device, designating a functional parameter of the blood vessel to be calculated.
For some applications, designating the parameter to be calculated includes, based upon the classification of the device as the given type of device, designating a dimension of the device to be calculated.
For some applications, calculating the parameter includes automatically calculating the parameter by analyzing the image using the processor.
For some applications, determining the classification of the device as the given type of device includes determining the classification of the device as a blood vessel filter, and designating the parameter to be calculated includes designating, as the parameter to be calculated, a maximum diameter of the blood vessel filter.
For some applications, determining the classification of the device as the given type of device includes determining the classification of the device as a hole-closure device, and designating the parameter to be calculated includes designating, as the parameter to be calculated, a maximum diameter of the hole-closure device.
For some applications:
determining the classification of the device as a given type of device includes determining that the device is a stent; and
designating the parameter to be calculated includes designating as the parameter to be calculated a minimum diameter of the blood vessel within a region of the blood vessel in which the stent is disposed.
For some applications:
determining the classification of the device as a given type of device includes determining that the device is a balloon; and
designating the parameter to be calculated includes designating as the parameter to be calculated a minimum diameter of the blood vessel within a region of the blood vessel in which the balloon is disposed.
For some applications, designating the parameter to be calculated includes designating as the parameter to be calculated the minimum diameter of the blood vessel within the region of the blood vessel in which the balloon is disposed, when maximal inflation of the balloon within the blood vessel occurs.
For some applications, determining the classification of the device includes receiving an input from a user that is indicative of the device being the given type of device, and determining the presence of the device within at least the portion of the blood vessel within the image includes automatically determining the presence of the device within at least the portion of the blood vessel within the image, by analyzing the image.
For some applications, determining the presence of the device within at least the portion of the blood vessel within the image includes determining the presence of the device within at least the portion of the blood vessel within the image, subsequent to receiving the input from the user that is indicative of the device being the given type of device.
There is further provided, in accordance with some applications of the present invention, a method including:
using at least one computer processor, detecting a device within a body of a subject; and
using the computer processor, generating an indication of a device-specific parameter on an output device, in response to detecting the device.
For some applications, detecting the device includes detecting a device within an image of the subject's body.
For some applications, detecting the device includes detecting the device by analyzing an image of at least a portion of the subject's body.
For some applications, generating the device-specific parameter includes generating an indication of a dimension of a portion of the subject's body.
For some applications, generating the device-specific parameter includes generating an indication of a functional parameter of a portion of the subject's body.
For some applications, generating the device-specific parameter includes generating an indication of a dimension of the device.
For some applications, generating the device-specific parameter includes calculating the parameter by analyzing an image of a portion of the subject's body, using the computer processor.
For some applications, detecting the device includes detecting the device and determining that the device is a blood vessel filter, and generating the device-specific parameter includes generating an indication of a maximum diameter of the blood vessel filter.
For some applications:
detecting the device includes determining that the device is a stent; and
generating the device-specific parameter includes generating an indication of a minimum diameter of a blood vessel within a region of the blood vessel in which the stent is disposed.
For some applications:
detecting the device includes determining that the device is a balloon; and
generating the device-specific parameter includes generating an indication of a minimum diameter of a blood vessel within a region of the blood vessel in which the balloon is disposed.
For some applications, generating the device-specific parameter includes generating an indication of a minimum diameter of the blood vessel within the region of the blood vessel in which the balloon is disposed, when maximal inflation of the balloon within the blood vessel occurs.
For some applications, detecting the device includes receiving an input from a user that is indicative of the device being a given type of device, and automatically determining the presence of the device within at least a portion of the subject's body.
For some applications, determining the presence of the device within the portion of the subject's body includes determining the presence of the device within the portion of the subject's body, subsequent to receiving the input from the user that is indicative of the device being the given type of device.
There is additionally provided, in accordance with some applications of the present invention, apparatus for use with an image of at least one blood vessel of a subject, including:
an output device; and
at least one computer processor configured to:
For some applications, the apparatus further includes a user interface, and the computer processor is configured to determine the presence of the device within the portion of the blood vessel within the image by receiving an input from a user, via the user interface, that is indicative of the presence of the device within the portion of the blood vessel within the image.
For some applications, the apparatus further includes a user interface, and the computer processor is configured to determine the classification of the device by receiving an input from a user, via the user interface, that is indicative of the device being the given type of device.
For some applications, the computer processor is configured to automatically determine the presence of the device within the portion of the blood vessel within the image, by analyzing the image.
For some applications, the computer processor is configured to automatically determine the classification of the device by analyzing the image.
For some applications, the computer processor is configured to designate the parameter to be calculated by, based upon the classification of the device as the given type of device, designating an event and designating a parameter to be calculated at the occurrence of the event.
For some applications, the computer processor is configured to designate the parameter to be calculated by, based upon the classification of the device as the given type of device, designating a dimension of the blood vessel to be calculated.
For some applications, the computer processor is configured to designate the parameter to be calculated by, based upon the classification of the device as the given type of device, designating a functional parameter of the blood vessel to be calculated.
For some applications, the computer processor is configured to designate the parameter to be calculated by, based upon the classification of the device as the given type of device, designating a dimension of the device to be calculated.
For some applications, the computer processor is configured to calculate the parameter by analyzing the image.
For some applications, the computer processor is configured to:
determine the classification of the device as the given type of device by determining that the device is a stent, and
designate the parameter to be calculated by designating, as the parameter to be calculated, a minimum diameter of the blood vessel within a region of the blood vessel in which the stent is disposed.
For some applications, the computer processor is configured to determine the classification of the device as the given type of device by determining the classification of the device as a hole-closure device, and the computer processor is configured to designate the parameter to be calculated by designating as the parameter to be calculated a maximum diameter of the hole-closure device.
For some applications, the computer processor is configured to determine the classification of the device as the given type of device by determining the classification of the device as a blood vessel filter, and the computer processor is configured to designate the parameter to be calculated by designating as the parameter to be calculated a maximum diameter of the blood vessel filter.
For some applications, the computer processor is configured to:
determine the classification of the device as the given type of device by determining that the device is a balloon, and
designate the parameter to be calculated by designating, as the parameter to be calculated, a minimum diameter of the blood vessel within a region of the blood vessel in which the balloon is disposed.
For some applications, the computer processor is configured to designate the parameter to be calculated by designating, as the parameter to be calculated, the minimum diameter of the blood vessel within the region of the blood vessel in which the balloon is disposed, when maximal inflation of the balloon within the blood vessel occurs.
For some applications:
the apparatus further includes a user interface, the computer processor is configured to determine the classification of the device by receiving an input from a user, via the user interface, that is indicative of the device being the given type of device, and
the computer processor is configured to automatically determine the presence of the device within the portion of the blood vessel within the image, by analyzing the image.
For some applications, the computer processor is configured to determine the presence of the device within at least the portion of the blood vessel within the image, subsequent to receiving the input from the user that is indicative of the device being the given type of device.
There is additionally provided, in accordance with some applications of the present invention, a computer software product, for use with an image of at least one blood vessel of a subject, the computer software product including a non-transitory computer-readable medium in which program instructions are stored, which instructions, when read by a computer cause the computer to perform the steps of determining a presence of a device within at least a portion of the blood vessel within the image; determining a classification of the device as a given type of device; based upon the classification of the device as the given type of device, designating a parameter to be calculated; automatically calculating the designated parameter; and generating an output in response to the calculated parameter.
There is further provided, in accordance with some applications of the present invention, a method for use with an image of a wire within a blood vessel of a subject, the wire including a radiopaque portion, the method including:
using at least one computer processor, automatically determining a location of a tip of the radiopaque portion of the wire within the image, by:
generating an output on an output device in response to the determined location of the tip of the radiopaque portion of the wire within the image.
For some applications, detecting that there is at least one pixel at which there is the change, by more than the threshold amount, in the wireness measure of the pixel includes detecting that there is at least one pixel disposed along a direction corresponding to a local direction of the length of the set of pixels within the region, at which there is the change, by more than the threshold amount, in the wireness measure of the pixel.
For some applications, detecting that there is at least one pixel at which there is the change, by more than the threshold amount, in the intensity of the pixel, includes detecting that there is at least one pixel disposed along a direction corresponding to a local direction of the length of the set of pixels within the region, at which there is the change, by more than the threshold amount, in the intensity of the pixel.
For some applications, the method further includes determining a location of the radiopaque portion of the wire within the image, based upon the determined location of the tip of the radiopaque portion of the wire, and
generating the output includes generating the output in response to the determined location of the radiopaque portion of the wire within the image.
For some applications, the apparatus further includes determining a location of a center of the radiopaque portion of the wire within the image, based upon the determined location of the tip of the radiopaque portion of the wire, and
generating the output includes generating the output in response to the determined location of the center of the radiopaque portion of the wire within the image.
For some applications, determining the wireness measure for each of the pixels within at least the portion of the image includes using machine-learning techniques.
For some applications, determining the wireness measure for each of the pixels within at least the portion of the image includes measuring, for each of the pixels, an extent to which the pixel, together with other pixels within the image, forms part of a continuous long, thin set of pixels having a given characteristic.
For some applications, determining the wireness measure for each of the pixels within at least the portion of the image includes analyzing eigenvalues of a multiscale second order local structure of at least the portion of the image.
For some applications, determining the wireness measure for each of the pixels within at least the portion of the image includes applying a filter that enhances curvilinear structures to at least the portion of the image.
For some applications, determining the wireness measure for each of the pixels within at least the portion of the image includes applying a filter that detects curvilinear structures to at least the portion of the image.
For some applications, determining the wireness measure for each of the pixels within at least the portion of the image includes applying a filter that segments curvilinear structures to at least the portion of the image.
For some applications, the method further includes aligning the image with a second image, based upon the determined location of the tip of the radiopaque portion of the wire within the image.
For some applications, generating the output includes, displaying the image and the second image in an image stream in which the image and the second image are aligned with one another.
For some applications, generating the output includes generating a composite image, based upon the alignment of the image and the second image.
For some applications, the method further includes determining a transformation function for mapping between the image and a second image at least partially in response to the determined location of the wire within the image, and
generating the output includes generating the output based upon the determined transformation function.
For some applications, the method further includes, based upon the determined transformation function, determining a location of an endoluminal device within the blood vessel, and generating the output includes generating the output in response to the determined location of the endoluminal device.
For some applications:
the endoluminal device includes an endoluminal data-acquisition device,
determining the location of the endoluminal device within the blood vessel includes determining a location within the blood vessel at which an endoluminal data point was acquired by the endoluminal data-acquisition device, and
generating the output includes generating the output based upon determining the location within the blood vessel at which the endoluminal data point was acquired by the endoluminal data-acquisition device.
For some applications, the method further includes, based upon the determined location of the wire within the image, determining a location of an endoluminal device within the blood vessel, and generating the output includes generating the output in response to the determined location of the endoluminal device.
For some applications:
the endoluminal device includes an endoluminal data-acquisition device,
determining the location of the endoluminal device within the blood vessel includes determining a location within the blood vessel at which an endoluminal data point was acquired by the endoluminal data-acquisition device, and
generating the output includes generating the output based upon determining the location within the blood vessel at which the endoluminal data point was acquired by the endoluminal data-acquisition device.
There is additionally provided, in accordance with some applications of the present invention, apparatus for use with an image of a wire within a blood vessel of a subject, the wire including a radiopaque portion, the apparatus including:
an output device; and
at least one computer processor configured to:
For some applications, the computer processor is configured to detect that there is at least one pixel at which there is the change, by more than the threshold amount, in the wireness measure of the pixel, by detecting that there is at least one pixel disposed along a direction corresponding to a local direction of the length of the set of pixels within the region, at which there is the change, by more than the threshold amount, in the wireness measure of the pixel.
For some applications, the computer processor is configured to detect that there is at least one pixel at which there is the change, by more than the threshold amount, in the intensity of the pixel, by detecting that there is at least one pixel disposed along a direction corresponding to a local direction of the length of the set of pixels within the region, at which there is the change, by more than the threshold amount, in the intensity of the pixel.
For some applications:
the computer processor is configured to determine a location of the radiopaque portion of the wire within the image, based upon the determined location of the tip of the radiopaque portion of the wire, and
the processor is configured to generate the output in response to the determined location of the radiopaque portion of the wire within the image.
For some applications:
the computer processor is configured to determine a location of a center of the radiopaque portion of the wire within the image, based upon the determined location of the tip of the radiopaque portion of the wire, and
the computer processor is configured to generate the output in response to the determined location of the center of the radiopaque portion of the wire within the image.
For some applications, the computer processor is configured to determine the wireness measure for each of the pixels within at least the portion of the image using machine-learning techniques.
For some applications, the computer processor is configured to determine the wireness measure for each of the pixels within at least the portion of the image by measuring, for each of the pixels, an extent to which the pixel, together with other pixels within the image, forms part of a continuous long, thin set of pixels having a given characteristic.
For some applications, the computer processor is configured to determine the wireness measure for each of the pixels within at least the portion of the image by analyzing eigenvalues of a multiscale second order local structure of at least the portion of the image.
For some applications, the computer processor is configured to determine the wireness measure for each of the pixels within at least the portion of the image by applying a filter that enhances curvilinear structures to at least the portion of the image.
For some applications, the computer processor is configured to determine the wireness measure for each of the pixels within at least the portion of the image by applying a filter that detects curvilinear structures to at least the portion of the image.
For some applications, the computer processor is configured to determine the wireness measure for each of the pixels within at least the portion of the image by applying a filter that segments curvilinear structures to at least the portion of the image.
For some applications, the computer processor is configured to align the image with a second image, based upon the determined location of the tip of the radiopaque portion of the wire within the image.
For some applications, the computer processor is configured to generate the output by displaying the image and the second image in an image stream in which the image and the second image are aligned with one another.
For some applications, the computer processor is configured to generate the output by generating a composite image, based upon the alignment of the image and the second image.
For some applications, the computer processor is configured to determine a transformation function for mapping between the image and a second image at least partially in response to the determined location of the wire within the image, and the computer processor is configured to generate the output based upon the determined transformation function.
For some applications, the computer processor is configured, based upon the determined transformation function, to determine a location of an endoluminal device within the blood vessel, and the computer processor is configured to generate the output in response to the determined location of the endoluminal device.
For some applications:
the endoluminal device includes an endoluminal data-acquisition device,
the computer processor is configured to determine the location of the endoluminal device within the blood vessel by determining a location within the blood vessel at which an endoluminal data point was acquired by the endoluminal data-acquisition device, and
the computer processor is configured to generate the output based upon determining the location within the blood vessel at which the endoluminal data point was acquired by the endoluminal data-acquisition device.
For some applications, the computer processor is configured, based upon the determined location of the wire within the image, to determine a location of an endoluminal device within the blood vessel, and the computer processor is configured to generate the output in response to the determined location of the endoluminal device.
For some applications:
the endoluminal device includes an endoluminal data-acquisition device,
the computer processor is configured to determine the location of the endoluminal device within the blood vessel by determining a location within the blood vessel at which an endoluminal data point was acquired by the endoluminal data-acquisition device, and
the computer processor is configured to generate the output based upon determining the location within the blood vessel at which the endoluminal data point was acquired by the endoluminal data-acquisition device.
There is further provided, in accordance with some applications of the present invention, a computer software product, for use with an image of a wire within a blood vessel of a subject, the wire including a radiopaque portion, the computer software product including a non-transitory computer-readable medium in which program instructions are stored, which instructions, when read by a computer cause the computer to perform the steps of automatically determining a location of a tip of a radiopaque portion of the wire within the image, by: for each pixel within at least a portion of the image determining a wireness measure of the pixel, by measuring an extent to which the pixel, together with other pixels within the image, forms part of a long, thin set of pixels having a given characteristic; for each pixel within at least the portion of the image, determining an intensity of the pixel; and determining that the tip of the radiopaque portion of the wire is located within a region within the portion of the image, by detecting that, within the region, there is at least one pixel at which there is a change, by more than respective threshold amounts, in both: the wireness measure of the pixel relative to the value of, the wireness measure of at least some of the set of pixels having the given characteristic, and the intensity of the pixel relative to the value of the intensity of at least some of the set of pixels having the given characteristic; and generating an output in response to the determined location of the tip of the radiopaque portion of the wire within the image.
The present invention will be more fully understood from the following detailed description of embodiments thereof, taken together with the drawings, in which:
Reference is now made to
It is noted that, for some applications, more than one processor is used. For some applications, more than one extraluminal imaging device is used with processor 20. For example, a first extraluminal imaging device may be used to acquire a first set of extraluminal images, and a second extraluminal imaging device may be used to acquire a second set of extraluminal images.
For some applications, endoluminal medical device 26 includes an endoluminal data-acquisition device that is configured to acquire data (e.g., functional data or images) from inside the subject's blood vessels. For some applications, the endoluminal data-acquisition device is an imaging probe. For some applications, the imaging probe is an IVUS probe, an EBUS probe, a different type of ultrasound probe, an OCT probe, an NIRS probe, an MR probe, a FACT probe, an OFDI probe, or any combination thereof. For some applications, the endoluminal data-acquisition device performs additional functions. For example, the endoluminal data-acquisition device may comprise a probe, such as the VIBE™ RX Vascular Imaging Balloon Catheter, marketed by Volcano Corporation (San Diego, USA), that includes both IVUS and coronary balloon functionalities. For some applications, the endoluminal data-acquisition device acquires data in a form other than images. For example, the data may include data related to pressure, flow, temperature, electrical activity, oxygenation, biochemical composition, or any combination thereof. For some applications, and typically when data are acquired with respect to a coronary vessel, the endoluminal data-acquisition device is a Fractional Flow Reserve (FFR) probe, and/or an instantaneous wave-free ratio (iFR) probe. For some applications, FFR and/or iFR measurements are determined by performing image-processing on extraluminal images, and the derived FFR and/or iFR measurements are co-registered with endoluminal images of the lumen, using techniques described herein. For some applications, FFR and/or iFR measurements are determined by performing image-processing on endoluminal images, and the derived FFR and/or iFR measurements are co-registered with extraluminal images of the lumen, using techniques described herein. For some applications, endoluminal images are co-registered with extraluminal images of the lumen, using techniques described herein, and FFR and/or iFR measurements are determined by performing image-processing on the co-registered images.
For some applications, endoluminal medical device 26 includes an endoluminal therapeutic device that is positioned and/or deployed at an anatomical feature that requires or potentially requires treatment, such as a partial or total occlusion, a native valve, an aneurism, a dissection, a malformation, a septal defect, a mass suspected of being malignant, a mass suspected of being inflammatory, etc. For example, the endoluminal therapeutic device may include a balloon (e.g., an angioplasty balloon), a stent, a valve, and/or a wire (e.g., a guide wire).
For some applications, apparatus and methods described herein are used with an endoluminal therapeutic device that is positioned and/or deployed at an implantation site of a previously-implanted device such as a stent, a graft or a replacement valve. The endoluminal data are determined at, and/or in the vicinity of, the implantation site. For example, the techniques described herein may be used during the placement of a new prosthetic aortic valve at the site of (e.g., inside) a previously implanted prosthetic aortic valve that is no longer functioning.
For some applications, apparatus and methods described herein are used with an endoluminal therapeutic device that is positioned and/or deployed at a defined location relative to a previously-implanted device such as a stent, a graft or a replacement valve. The endoluminal data are determined at, and in the vicinity of, the defined location. For example, the techniques described herein may be used during the placement of a coronary stent such that the new stent overlaps with or is adjacent to a previously-implanted stent, in order to treat a long lesion and/or a lesion that has diffused along a coronary artery.
For some applications, output device 40 is a display that is configured to display an extraluminal image 42 of a blood vessel (e.g., a fluoroscopic image), an endoluminal image of a blood vessel 44 (e.g., an IVUS image), and or a stack 46 of cross-sections of endoluminal images (e.g., a stack of IVUS images).
Reference is now made to
Typically, the quantitative vessel analysis is performed by determining the locations of vessel centerlines and/or edge lines, for example, using techniques such as those described in US 2012/0230565 to Steinberg, and/or US 2010/0222671 to Cohen, both of which are incorporated herein by reference. For some applications, a lesion is automatically detected in accordance with the following procedure.
Scan lines are generated perpendicular to the centerline of a segment of the vessel that is sampled. The image is resampled along the scan lines. Corresponding gray-level values are stored as columns of a rectangular matrix M, thereby resampling the segment of the vessel as a straightened vessel segment. For the straightened vessel segment, optimum upper and lower paths are determined (with respect to the middle row of M), which connect the first and last columns of M. The optimization criterion takes into account the changes in gray-level along columns of M, and the paths' slopes. The vessel edge lines are obtained via back projection of upper and lower optimal paths on the original image.
A shortest path algorithm (e.g., as described in an article by Dijkstra, entitled “A Note on Two Problems in Connexion with Graphs” (Numerische Mathematik 1, 269-271, 1959), which is incorporated herein by reference) is used in order to avoid irregularities, such as small gaps and loops, in the edge lines. For some applications, the centerline is corrected based upon the detected edge lines, and new scan lines are constructed. For each new scan line, vessel diameter is defined as a distance between the two points where the scan line intersects vessel boundaries.
For some applications, if the processor uses the algorithm described hereinabove for performing quantitative vessel analysis on the vessel with the balloon inside, the processor may identify one or both ends 58 of the balloon as being the location of the minimal lumen diameter, since the system may not differentiate between edge lines of the vessel and edge lines of the balloon. Therefore, for some applications, in order to avoid the processor identifying one or both ends 58 of the balloon as being the location of the minimal lumen diameter, the processor determines that the balloon is present within the blood vessel. The processor determines a parameter of the vessel responsively to determining the presence of the balloon within the vessel, for example, in accordance with the techniques described hereinbelow with reference to
Reference is now made to
As shown in
As shown in
Subsequent to designating the parameter to be calculated (step 65), the processor calculates the designated parameter (step 66) and generates an output in response thereto (step 68). For example, in the example of angioplasty balloon 55, shown in
As shown in
The processor designates as a region of the vessel in which the balloon is disposed, a region of the vessel that is between the radiopaque markers of the balloon, and/or a region of the vessel that is between the tapered pairs of edge lines at each longitudinal end of the balloon (step 74). The processor then determines the minimum lumen diameter within the region of the vessel in which the balloon is disposed (step 76). The minimum lumen diameter within the region of the vessel in which the balloon is disposed is the residual waist of the occlusion. The processor then generates an output that is indicative of the calculated residual waist (step 78).
For some applications, the detection and/or classification of the device (steps 63 and 64 of
For some applications, the processor designates the parameter to be calculated (step 65 of
Alternatively or additionally, in response to classifying the device as a stent, the processor may designate functional flow reserve (or another luminal flow-related index) at the location of the stent as the parameter to be calculated, in order to determine the effect of the stent on the functional flow reserve (or the other luminal flow-related index) of the vessel. For some applications, the processor designates the parameter to be calculated by designating a parameter of the device to be calculated. For example, in response to classifying the device as a stent, the processor may designate a maximum diameter of the stent, or a minimum diameter of the stent as the parameter to be calculated. For some applications, the stent is disposed around a balloon, and the processor determines the region of the blood vessel in which the stent is disposed by determining locations of the radiopaque markers of the balloon around which the stent is disposed.
For some applications, the processor designates an event and designates the parameter to be calculated by designating a parameter to be calculated at the occurrence of the event. For example, in the example described with reference to
Typically, the parameter is calculated (step 66 of
For some applications, in order to calculate the parameter, vessel and/or device edge lines are automatically identified, using techniques described herein. For example, scan lines may be generated perpendicular to the centerline of a segment of the vessel that is sampled. The image is resampled along the scan lines. Corresponding gray-level values are stored as columns of a rectangular matrix M, thereby resampling the segment of the vessel as a straightened vessel segment. For the straightened vessel segment, optimum upper and lower paths are determined (with respect to the middle row of M), which connect the first and last columns of M. The optimization criterion takes into account the changes in gray-level along columns of M, and the paths' slopes. The vessel edge lines are obtained via back projection of upper and lower optimal paths on the original image.
A shortest path algorithm (e.g., as described in an article by Dijkstra, entitled “A Note on Two Problems in Connexion with Graphs” (Numerische Mathematik 1, 269-271, 1959), which is incorporated herein by reference) is used in order to avoid irregularities, such as small gaps and loops, in the edge lines. For some applications, the centerline is corrected based upon the detected edge lines, and new scan lines are constructed. For each new scan line, vessel and/or device diameter is defined as a distance between the two points where the scan line intersects edge lines.
For some applications, the techniques described herein are performed with respect to other devices, and/or with respect to other portions of a subject's body to those described hereinabove.
For some applications, the techniques described herein are used in order to determine a parameter that relates to a hole-closure device. For example, the hole-closure device may be an atrial septal defect closure device, a left-atrial appendage closure device, and/or a hole-closure device that is used to close a surgically-created hole, such as in the apex of the subject's heart, and/or in a peripheral blood vessel, such as the femoral vein or the femoral artery. In response to determining the presence of a device within an image of a portion of the subject's body, and classifying the device as a hole-closure device, computer processor 28 may determine the maximum diameter of the hole-closure device, subsequent to the deployment of the hole-closure device. Alternatively or additionally, the techniques described herein may be used in order to determine a parameter that relates to an implantable valve (such as a prosthetic aortic valve, and/or a prosthetic mitral valve), e.g., the maximum diameter of the valve, subsequent to deployment of the valve. Further alternatively or additionally, the techniques described herein may be used in order to determine a parameter that relates to a blood vessel filter (e.g., a vena cava filter, such as the Crux® Vena Cava Filter, manufactured by Volcano Corporation (CA, USA)), e.g., the maximum diameter of the filter, subsequent to deployment of the filter within a blood vessel.
For some applications, the techniques described herein are used in order to determine a parameter that is related to a previously-implanted device (such as a stent, a graft or a replacement valve that was implanted prior to the present procedure being performed (e.g., at least one day prior to the present procedure being performed), in response to determining a presence of the previously-implanted device within an image of a portion of the subject's body, and in response to classifying the previously-implanted device as a given type of device.
Reference is now made to
Reference is now made to
For each pixel within at least a portion of the image, a wireness measure of the pixel is determined (step 90). (It is noted that as used herein, the term “pixel” should not be interpreted to be limited to the smallest controllable element of the picture represented on the screen. Rather, as used herein, the term “pixel” should be interpreted to mean a set of one or more such elements.) For some applications, wireness measures are determined for pixels within the entire image. Alternatively, wireness measures of pixels within only a portion of the image are determined. For example, wireness measures of pixels within a region of the image in which the radiopaque portion of the wire is expected to be disposed (e.g., based on locations of other features within the image) may be sampled. Alternatively or additionally, the processor may receive an input from the user indicating a region of the image in which the radiopaque portion of the wire is expected to be disposed.
The wireness measure is a measure of the extent to which each of the pixels has a wire-like characteristic. For some applications, the wireness measure of the pixel is determined using one or more algorithms described herein. Typically, for each of the selected pixels, the wireness measure is determined, by measuring an extent to which the pixel, together with other pixels within the image, forms part of a long, thin set of pixels having a given characteristic, such as darkness, brightness, and/or a different characteristic (i.e., a wireness-indicating characteristic). Typically, the wireness measure is indicative of the pixel, together with the other pixels within the set, corresponding to the wire.
For some applications, generally similar techniques to those described in US 2012/0230565 to Steinberg, and/or US 2010/0222671 to Cohen, both of which are incorporated herein by reference, for determining a vesselness measure of a pixel are used for determining the wireness measure of a pixel. For example, wireness may be determined by means of a Hessian filter, such as the filter described in the article by Frangi et al., entitled “Multiscale vessel enhancement filtering” (Medical Image Computing and Computer Assisted Intervention—MICCAI 1998—Lecture Notes in Computer Science, vol. 1496, Springer Verlag, Berlin, Germany, pp. 130-137), which is incorporated herein by reference, and/or by means of a filter that performs enhancement and/or detection and/or segmentation of curvilinear structures. For some applications, a filter is used that is similar to a Frangi filter, but that differs from a Frangi filter (a) in that wireness is a homogeneous function, and/or (b) in the multipliers employed for the normalization of scales.
For some applications, the wireness measure of a pixel is obtained by determining the extent to which the gradient of the pixel is orthogonal to the eigenvector of the Hessian matrix corresponding to the highest eigenvalue. For some applications, the determination is assisted by a voting function applied to pixels that are adjacent to those pixels that are eventually determined to constitute the wire itself.
For some applications, thresholding is applied to image pixels by means of hysteresis. For example, pixels the wireness values of which fall below the high threshold of the hysteresis, but yet above the low threshold of the hysteresis, are incorporated into the set of pixels if they are contiguous with pixels that fall at or above the high threshold of the hysteresis.
For some applications, the pixels which form the aforementioned set are determined by means of morphological operations. For example, such morphological operations may include the skeletonization of a thresholded vesselness image. For some applications, the threshold applied is adaptive according to the specific region in the image.
For some applications, machine-learning techniques are used to determine wireness measures of the pixels.
In the next step of the procedure, for each pixel within at least the portion of the image, an intensity of the pixel is determined (step 92). It is noted that, for some applications, steps 90 and 92 are performed in the reverse order to that shown in the flowchart in
Subsequent to the wireness measures and the intensities of the pixels within the portion of the image having been determined, computer processor 28 designates a first sampling region within the portion of the image (step 94). For some applications, the first sampling region is designated in accordance with one or more algorithms described herein. For example, the sampling region may include a single pixel, or a plurality of pixels. The first sampling region may be generated randomly, and/or in response to an input from a user. For some applications, the processor designates the first sampling region by designating a sampling region in which the tip of the guidewire is likely to be disposed. For example, the processor may designate the first sampling region by designating one or more regions that have high values of the wireness measure. Alternatively or additionally, the processor may designate the first sampling region in response to determining that a region of the image is likely to contain the tip of the wire, based upon machine-learning analysis of the image.
For some applications, step 94 is performed before steps 90 and 92, and steps 90 and 92 are only performed on the designated sampling region.
The processor determines whether, within the first sampling region, there is at least one pixel, at which there is a change in the wireness measure, by more than a threshold amount, relative to the value of the wireness measure of one or more pixels that are adjacent to the pixel and that have the wireness-indicating characteristic (step 96). For some applications, the processor performs step 96 using one or more algorithms described herein. For example, by way of illustration, the processor may determine that from one pixel to an adjacent pixel there is a decrease in the wireness measure that exceeds a threshold percentage decrease. Or, the processor may determine that at least one pixel has a wireness measure that is lower than the mean wireness measure of all of the pixels belonging to the set of pixels, by more than a threshold percentage, whereas one or more pixels that are adjacent to the pixel have wireness measure(s) that exceeds the threshold.
For some applications, within the designated region, the processor determines the local directionality of the set of pixels that have the wireness-indicating characteristic. The processor determines, along that direction, whether there is at least one pixel at which there is a change in the wireness measure by more than a threshold amount, relative to one or more pixels that are adjacent to the pixel and that have the wireness-indicating characteristic.
In response to determining that there is not a change in the wireness measure at at least one pixel that exceeds the threshold, within the first sampling region, the processor proceeds to the next sampling region (step 98), and step 96 is repeated at the second sampling region. The second sampling region is typically selected in a generally similar manner to the selection of the first sampling region, and/or based upon a spatial relationship to the first sampling region. In response to determining that the change in the wireness measure at at least one pixel within a sampling region does exceed the threshold, the processor then determines whether, at the at least one pixel there is a change, by more than a threshold amount, in the intensity of the pixel relative to the value of the intensity of at least some of the set of pixels having the wireness-indicating characteristic (step 100). For some applications, the processor performs step 100 using one or more algorithms described herein. For example, by way of illustration, the processor may determine that from one of the pixels to an adjacent pixel there is an increase in intensity that exceeds a threshold percentage increase. Or, the processor may determine that one of the pixels has an intensity that exceeds the mean intensity of all of the pixels belonging to the set of pixels, by more than a threshold percentage.
For some applications, within the designated region, the processor determines the local directionality of the set of pixels that have the wireness-indicating characteristic. The processor determines, along that direction, whether there is at least one pixel at which there is a change in intensity by more than a threshold amount, relative to at least some of the pixels that belong to the set of pixels that have the wireness-indicating characteristic.
In response to determining that there is not a change in the intensity at at least one pixel that exceeds the threshold, within the current sampling region, the processor proceeds to the next sampling region (step 98), and step 96 is repeated at the next sampling region. In response to determining that the change in intensity at at least one pixel within the current sampling region does exceed the threshold, it is determined that the tip of the radiopaque portion of the wire is disposed within the current sampling region (step 102). An output is generated in response to the determined location of the tip of the radiopaque portion of the wire within the image.
It is noted that, for some applications, steps 96 and 100 are performed in the reverse order, the processor first determining whether there is a change in the intensity at at least one pixel by more than a threshold amount, and, subsequently, determining whether there is a change in the wireness measure at the at least one pixel by more than the threshold amount.
Typically, an output is generated by computer processor 28 in response to the determined location of the tip of the radiopaque portion of the wire within the image. For some applications, the processor determines locations of both of the tips of the radiopaque portion of the wire, and thereby determines the location of the radiopaque portion of the wire, and/or the center of the radiopaque portion of the wire within the image. The output is generated in response to the determined location of the radiopaque portion of the wire, and/or in response to the determined location of the center of the radiopaque portion of the wire within the image.
For some applications, in response to determining the location of the tip of the radiopaque portion of the wire within a given image, the processor aligns the image with a second image, by aligning the radiopaque portions of the wires in each of the images with one another. In accordance with respective applications, the aligned images may be displayed in an image stream in which the images are aligned with one another, and/or a composite image may be generated based upon the alignment of the image and the second image. For example, the processor may average the image and the second image, subsequent to aligning the images with one another. In general, the identification of the location of the tip of the radiopaque portion of the wire within a given image may be used to perform any of the image stabilization and/or enhancement techniques that are described in any one of US 2008/0221440 to Iddan, US 2012/0230565 to Steinberg, and US 2010/0222671 to Cohen, all of which applications are incorporated herein by reference.
For some applications, identification of the location of the tip of the radiopaque portion of the wire within a given image may be used to facilitate the determination of the location of an endoluminal device within the lumen, for example, in accordance with techniques described in US 2012/0004537 to Tolkowsky, and/or WO 13/174472 to Steinberg, which are incorporated herein by reference. For example, the location within the blood vessel at which one or more endoluminal data points were acquired by the endoluminal data-acquisition device (e.g., an endoluminal imaging device, or an endoluminal data-acquisition device that is configured to acquire a plurality of functional endoluminal data points) may be determined. Based upon the determined locations within the blood vessel at which the endoluminal data points were acquired, the processor may generate an output, such as by generating an endoluminal imaging stack, and/or by generating an indication of the correspondence between an endoluminal data point and the location within the blood vessel at which the endoluminal data point was acquired.
For some applications, endoluminal data points are acquired by positioning an endoluminal data-acquisition device along a luminal segment of the blood vessel that includes a designated luminal site. Subsequently, while observing extraluminal images of the luminal segment, one or more locations along that segment are indicated by a user input device (e.g., user interface 30). In response to the indication of the one or more locations by the user input device, the corresponding, previously-acquired endoluminal images are displayed.
Typically, the designated luminal site includes a site being diagnosed, at which, subject to the outcome of the diagnosis, a therapeutic device will be positioned and deployed, e.g., the site of an anatomical feature, the implantation site of a previously-implanted device, and/or a site at a defined location with respect to the implantation site. For example, the designated luminal site may include a portion of the lumen that is narrow with respect to surrounding portions of the lumen, and/or the site of a lesion.
For some applications, endoluminal data points are acquired by positioning an endoluminal data-acquisition device at a designated luminal site. Subsequently, an endoluminal therapeutic device is positioned and deployed at the designated luminal site under extraluminal imaging, while concurrently viewing on-line the endoluminal data that were previously acquired by the endoluminal data-acquisition device at the current location of the therapeutic device. Typically, endoluminal data are acquired at respective endoluminal sites in the vicinity of the designated endoluminal site. Subsequently, when the endoluminal therapeutic device is placed inside the lumen, previously-acquired endoluminal data are displayed and updated, typically automatically and typically on-line, to correspond to the current location of the therapeutic device (or of a portion thereof), the location of the therapeutic device typically changing during the positioning of the therapeutic device.
For some applications, extraluminal imaging and the previously-acquired endoluminal data points are co-used such that it is as if the therapeutic device is being positioned and deployed under both real-time extraluminal imaging and real-time endoluminal data acquisition. This is because (a) the extraluminal imaging is performed in real-time, and (b), although the endoluminal data are not acquired in real-time, endoluminal data are displayed that correspond to the current location of the therapeutic device.
In accordance with some applications of the present invention, when the therapeutic device is disposed inside the lumen, the location of the device within the lumen is determined by performing image processing on the extraluminal image of the device inside the lumen.
For some applications, identification of the location of the tip of the radiopaque portion of the wire within a given image may be used to facilitate the determination of a transformation function for mapping between the image and a second image, for example, in accordance with techniques described in WO 13/174472 to Steinberg, which is incorporated herein by reference. For example, a transformation function may be determined for mapping a current fluoroscopic image to a previously acquired angiographic image, or vice versa. For some applications, a transformation function is determined for mapping a current fluoroscopic image to a previously acquired angiographic image, by comparing an arrangement of two or more features within the current fluoroscopic image to a shape of at least a portion of a pathway within the previously acquired angiographic image. For some applications, at least one of the features is the radiopaque portion of the guidewire in the current fluoroscopic image.
For some applications, based upon the determined transformation function, the processor determines the location of an endoluminal device within the lumen, for example, in accordance with techniques described in WO 13/174472 to Steinberg, which is incorporated herein by reference. For example, the location within the blood vessel at which one or more endoluminal data points were acquired by the endoluminal data-acquisition device (e.g., an endoluminal imaging device, or an endoluminal data-acquisition device that is configured to acquire a plurality of functional endoluminal data points) may be determined. Based upon the determined locations within the blood vessel at which the endoluminal data points were acquired, the processor may generate an output, such as by generating an endoluminal imaging stack, and/or by generating an indication of the correspondence between an endoluminal data point and the location within the blood vessel at which the endoluminal data point was acquired. Alternatively, based upon the determined location within the blood vessel at which the endoluminal data point was acquired, the processor may co-use endoluminal data points and extraluminal imaging using techniques described hereinabove, and/or as described in US 2012/0004537 to Tolkowsky, and/or WO 13/174472 to Steinberg, which are incorporated herein by reference.
For some applications, identification of the location of the tip of the radiopaque portion of the wire within a given image may be used to facilitate the classification of a portion of the image as being associated with the distal end of a guidewire, for example, in accordance with techniques described in WO 13/174472 to Steinberg, which is incorporated herein by reference. For example, classification of a portion of the image as being associated with the distal end of a guidewire may be used to facilitate the determination of a transformation function for mapping one image to another image, in accordance with techniques described in WO 13/174472 to Steinberg, which is incorporated herein by reference.
It is noted that although some techniques described herein are described primarily with respect to extraluminal fluoroscopic/angiographic images and endoluminal images, the scope of the present invention includes applying the techniques described herein to other forms of extraluminal and endoluminal images and/or data, mutatis mutandis. For example, the extraluminal images may include images generated by fluoroscopy, CT, MRI, ultrasound, PET, SPECT, other extraluminal imaging techniques, or any combination thereof. Endoluminal images may include images generated by intravascular ultrasound (IVUS) optical coherence tomography (OCT), near-infrared spectroscopy (NIRS), intravascular ultrasound (IVUS), endobronchial ultrasound (EBUS), magnetic resonance (MR), other endoluminal imaging techniques, or any combination thereof. Endoluminal data may include data related to pressure (e.g., fractional flow reserve), flow, temperature, electrical activity, or any combination thereof.
Although some techniques described herein are described primarily as being performed on a blood vessel, the scope of the present application includes performing similar techniques on a lumen in the vascular system, the respiratory tract, the digestive tract, the urinary tract, any other luminal structure within a patient's body, or any other suitable anatomical structure within a patient's body, mutatis mutandis. Examples of an anatomical structure to which the techniques described herein may be applied include a coronary vessel, a coronary lesion, a vessel, a vascular lesion, a lumen, a luminal lesion, and/or a valve.
Applications of the invention described herein can take the form of a computer program product accessible from a computer-usable or computer-readable medium providing program code for use by or in connection with a computer or any instruction execution system, such as computer processor 28. For the purposes of this description, a computer-usable or computer readable medium can be any apparatus that can comprise, store, communicate, propagate, or transport the program for use by or in connection with the instruction execution system, apparatus, or device. The medium can be an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system (or apparatus or device) or a propagation medium. Typically, the computer-usable or computer readable medium is a non-transitory computer-usable or computer readable medium.
Examples of a computer-readable medium include a semiconductor or solid state memory, magnetic tape, a removable computer diskette, a random access memory (RAM), a read-only memory (ROM), a rigid magnetic disk and an optical disk. Current examples of optical disks include compact disk-read only memory (CD-ROM), compact disk-read/write (CD-R/W) and DVD.
A data processing system suitable for storing and/or executing program code will include at least one processor (e.g., computer processor 28) coupled directly or indirectly to memory elements (e.g., memory 29) through a system bus. The memory elements can include local memory employed during actual execution of the program code, bulk storage, and cache memories which provide temporary storage of at least some program code in order to reduce the number of times code must be retrieved from bulk storage during execution. The system can read the inventive instructions on the program storage devices and follow these instructions to execute the methodology of the embodiments of the invention.
Network adapters may be coupled to the processor to enable the processor to become coupled to other processors or remote printers or storage devices through intervening private or public networks. Modems, cable modem and Ethernet cards are just a few of the currently available types of network adapters.
Computer program code for carrying out operations of the present invention may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C++ or the like and conventional procedural programming languages, such as the C programming language or similar programming languages.
It will be understood that each block of the flowcharts shown in
Computer processor 28 is typically a hardware device programmed with computer program instructions to produce a special purpose computer. For example, when programmed to perform the algorithms described with reference to
The scope of the present application includes combining the apparatus and methods described herein with apparatus and methods described in any one of the following applications, all of which are incorporated herein by reference:
It will be appreciated by persons skilled in the art that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and subcombinations of the various features described hereinabove, as well as variations and modifications thereof that are not in the prior art, which would occur to persons skilled in the art upon reading the foregoing description.
The present application claims priority from U.S. Provisional Patent Application 61/977,891 to Klaiman, filed Apr. 10, 2014, entitled “Image analysis in the presence of a medical device,” which is incorporated herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/IL15/05372 | 4/2/2015 | WO | 00 |
Number | Date | Country | |
---|---|---|---|
61977891 | Apr 2014 | US |