SYSTEM AND METHOD FOR OCT DEPTH CALIBRATION
The invention generally relates to methods for calibrating an OCT image by comparing a detected location of a feature within the image to an expected location of the feature, calculating a calibration value, and transforming the location of pixels within the image to provide a calibrated image.
Latest VOLCANO CORPORATION Patents:
- Device and Method for Determining the Likelihood of a Patient Having a Clinical Event or a Clinically Silent Event Based on Ascertained Physiological Parameters
- Vibrating guidewire torquer and methods of use
- Systems for indicating parameters in an imaging data set and methods of use
- Intravascular ultrasound imaging apparatus, interface, architecture, and method of manufacturing
- Retrieval snare device and method
This application claims the benefit of, and priority to, U.S. Provisional Patent Application No. 61/778,761, filed Mar. 13, 2013, which is incorporated by reference.
FIELD OF THE INVENTIONThe invention relates to methods for calibrating optical coherence tomography images.
BACKGROUNDIntravascular optical coherence tomography (OCT) systems operate via a rotating catheter that is introduced into a patient's blood vessels to image the vessel walls to detect atherosclerotic plaque and other life-threatening conditions. The catheter takes a picture by shooting beams of light as it rotates and translates within the vessel and detecting light that is reflected back. The light is detected in the form of a series of lines arranged in a helix. A computer processor can digitize the information and arrange it for presentation on a screen.
One problem in OCT is calibration. The lines are treated as if they originate at the center of the fiber optical catheter, but this is based on assumptions about the path length of light through the catheter. The catheter is, in-fact, subject to stresses that can distort the path length during imaging. Thus, if the catheter stretches by a millimeter during imaging, the radius to a point on each line is artificially shrunken by a millimeter. If a diameter is measured on that image, the diameter will be wrong by two millimeters.
Unfortunately, the fiber optical catheter is within a patient during imaging and the imaging proceeds rapidly, so it is not realistic for an attendant to stand by and measure stretching or distortion of the catheter. Since one measure of the severity of atherosclerosis is cross-sectional area of space within a blood vessel, path length changes in OCT imaging hardware interfere with accurately evaluating the severity of a patient's condition.
SUMMARYThe invention provides systems and methods for calibrating OCT images by detecting a feature within the image, where a position of the feature within the actual patient is known. For example, a surface of an imaging catheter can be detected automatically within the OCT image. If the diameter of the catheter is known a priori, the detected location can be compared to the expected location to determine a path length distortion. The image can be transformed to correct for the path length distortion, thereby repositioning all of the pixels into locations that accurately depict the bodily tissue. This provides a calibrated image in which apparent dimensions in the image display (e.g., as seen on a computer monitor) can be used to faithfully determine actual dimensions of medically significant features within a patient's body. Thus OCT imaging can be used to evaluate the severity of atherosclerosis by, for example, revealing the actual space through a blood vessel through which blood may flow.
In certain aspects, the invention provides a method of calibrating an imaging system by obtaining an image that includes a target and a reference item, detecting a location of the reference item within the image, determining a y-value of the location, comparing the determined y-value to a stored reference y-value, calculating a calibration value based on the comparison, and providing a calibrated image by shifting pixels in a y-direction according to the calibration value. The imaging system may be an optical coherence tomography system. The reference item may be an image of a catheter sheath. The detecting step may include a morphological image processing operation. The method may include digitally transforming image data to provide the calibrated image. Detecting the location may include peak searching, correlation, thresholding, or performing a pattern recognition algorithm. The method may involve obtaining an a priori range over which a path length displacement may occur. Determining the y-value may include averaging across the image. Calculating the calibration value may include taking a difference between the determined y-value and the stored reference y-value.
In some embodiments, a feature is also removed (cropped) from the image. For example, the detected reference item may be removed from the image. Removing may be by cropping out the feature across its extent in the y direction at all x positions. Cropping or removing a feature may include calculating an average value of pixels proximal to the cropped region and inserting the average values to replace the cropped pixels.
In certain embodiments, any or all of the foregoing steps are performed using a computer system comprising a tangible, non-transitory memory coupled to a processor. In related aspects, the invention provides a system for providing a calibrated image, in which the system includes a memory coupled to a processor, the system being operable to perform any of the foregoing steps.
The invention provides systems and methods for calibrating an imaging system. Systems and methods of the invention have application in imaging systems that require calibration to provide a scale. Exemplary systems include imaging and sensing systems based on principles of rotational imaging. In some embodiments, systems and applications contemplated for use with the invention include optical coherence tomography (OCT).
In OCT systems, a light source is used to provide a beam of coherent light. The light source can include an optical gain medium (e.g., laser or optical amplifier) to produce coherent light by stimulated emission. In some embodiments, the gain medium is provided by a semiconductor optical amplifier. A light source may further include other components, such as a tunable filter that allows a user to select a wavelength of light to be amplified. Wavelengths commonly used in medical applications include near-infrared light, for example between about 800 nm and about 1700 nm.
Generally, there are two types of OCT systems, common beam path systems and differential beam path systems, that differ from each other based upon the optical layout of the systems. A common beam path system sends all produced light through a single optical fiber to generate a reference signal and a sample signal whereas a differential beam path system splits the produced light such that a portion of the light is directed to the sample and the other portion is directed to a reference surface. Common beam path systems are further described for example in U.S. Pat. No. 7,999,938; U.S. Pat. No. 7,995,210; and U.S. Pat. No. 7,787,127 the contents of each of which are incorporated by reference herein in their entirety.
In a differential beam path system, the coherent light from the light source is input into an interferometer and split into a reference path and a sample path. The sample path is directed to the target and used to image the target. Reflections from the sample path are joined with the reference path and the combination of the reference-path light and the sample-path light produces interference patterns in the resulting light. The light, and thus the patterns, are converted to electric signals, which are then analyzed to produce depth-resolved images of the target tissue on a micron scale. Exemplary differential beam path interferometers are Mach-Zehnder interferometers and Michelson interferometers. Differential beam path interferometers are further described for example in U.S. Pat. No. 7,783,337; U.S. Pat. No. 6,134,003; and U.S. Pat. No. 6,421,164, the contents of each of which are incorporated by reference herein in its entirety.
Commercially available OCT systems are employed in diverse applications, including art conservation and diagnostic medicine, notably in ophthalmology where OCT can be used to obtain detailed images from within the retina. The detailed images of the retina allow one to identify diseases and trauma of the eye. Other applications of imaging systems of the invention include, for example, dermatology (e.g., to image subsurface structural and blood flow formations), dentistry (to image teeth and gum line), gastroenterology (e.g., to image the gastrointestinal tract to detect polyps and inflammation), and cancer diagnostics (for example, to discriminate between malignant and normal tissue).
In certain embodiments, systems and methods of the invention image within a lumen of tissue. Various lumen of biological structures may be imaged including, for example, blood vessels, including, but not limited, to vasculature of the lymphatic and nervous systems, various structures of the gastrointestinal tract including lumen of the small intestine, large intestine, stomach, esophagus, colon, pancreatic duct, bile duct, hepatic duct, lumen of the reproductive tract including the vas deferens, vagina, uterus and fallopian tubes, structures of the urinary tract including urinary collecting ducts, renal tubules, ureter, and bladder, and structures of the head and neck and pulmonary system including sinuses, parotid, trachea, bronchi, and lungs. Systems and methods of the invention have particular applicability in imaging veins and arteries such as, for example, the arteries of the heart. Since an OCT system can be calibrated to provide scale information, intravascular OCT imaging of the coronary arteries can reveal plaque build-up over time, change in dimensions of features, and progress of thrombotic elements. The accumulation of plaque within the artery wall over decades is the setup for vulnerable plaque which, in turn, leads to heart attack and stenosis (narrowing) of the artery. OCT images, if scaled or calibrated, are useful in determining both plaque volume within the wall of the artery and/or the degree of stenosis of the artery lumen. Intravascular OCT can also be used to assess the effects of treatments of stenosis such as with hydraulic angioplasty expansion of the artery, with or without stents, and the results of medical therapy over time.
As shown in
Light source 827, as discussed above, may use a laser or an optical amplifier as a source of coherent light. Coherent light is transmitted to interferometer 831.
An image is captured by introducing imaging catheter 826 into a target within a patient, such as a lumen of a blood vessel. This can be accomplished by using standard interventional techniques and tools such as a guide wire, guide catheter, or angiography system. Suitable imaging catheters and their use are discussed in U.S. Pat. No. 8,116,605 and U.S. Pat. No. 7,711,413, the contents of which are incorporated by reference in their entirety for all purposes.
Looking back at
After combining light from the sample, and reference paths, the combined light from splitter 919 is split into orthogonal polarization states, resulting in RF-band polarization-diverse temporal interference fringe signals. The interference fringe signals are converted to photocurrents using PIN photodiodes 929a, 929b, . . . on the OCB 851 as shown in
Data is collected from A-scans A11, A12, . . . , AN, as shown in
The invention provides methods and systems for detecting and using OCT image features, either intentionally generated or artifactual, for automatically adjusting the depth range in polar (“radar-like”) OCT images. Circular or cylindrical OCT scanning devices sample physical space in an inherently polar coordinate system (e.g. radius and angle rather than length and width). However, digital representations of images (i.e. arrays of pixels representing numeric values) are inherently rectangular as shown in
Polar OCT images must be converted from their rectangular representation (from
For example: y=0 (the top row of the rectangular image) maps to radius=0 (the center of the polar image) and y=ymax (the bottom row of the rectangular image) maps to radius=ymax (the perimeter of the polar image). Likewise, x=0 (the left column in the rectangular image) maps to angle=0° and x=xmax/2 maps to approximately 180° and x=xmax maps to an angle of approximately 359°. (It is noted that the assignment of x or y to rows or columns is arbitrary and non-limiting. Herein, “x” and “y” could be reversed in any example, meaning only that a rectangular image would be stored in the computer memory as a transpose relative to the other.)
For accurate quantitative dimensional measurement in polar images, pixels mapping to radius=0 must represent the actual physical space at the center of the axis of rotation of the imaging probe, otherwise the polar image will be artificially warped (expanded or contracted) in the radial direction. However, in an arbitrary OCT image, the pixels at y=0 do not necessarily satisfy this requirement and must be shifted in the y dimension until this is satisfied before mapping to a polar representation. Differential displacements (either controlled or uncontrolled) in the path length of the sample vs. reference arms of the interferometer will shift the pixels in the y dimension.
Displacements can occur when using cylindrical (typically actually helical) scanning fiber-optic OCT catheters. For example, when the catheter is pushed or pulled longitudinally, the fiber-optic cable can be compressed or stretched and thus a path length displacement is incurred.
The method disclosed herein is of automatic recognition of the uncontrolled displacement effect based on searching for image features which should be stationary (but are not due to uncontrollable displacement), and successive calibration of OCT image data so that polar representations can then be used for accurate dimensional measurements.
Additionally, a method is provided for subsequent removal of said features in image prior to display.
Image features used by the methods and systems disclosed herein are preferably generated within the catheter itself (i.e., not within the imaged subject or surroundings) and should appear somewhat stable in depth and consistent in intensity throughout the 360° rotation of the catheter. These include but are not limited to back reflections at interfaces between optical components (aka “ghost-lines” or “echo artifacts”, these occur along the optical axis of rotating parts and thus appear as uniform circles in the polar image when no differential path length displacement occurs over the course of one catheter rotation), or reflections from the boundaries of or from within the stationary (non-rotating) catheter sheath (if it is circular in cross-sectional profile and also mechanically concentric with the rotating portion).
In some embodiments, steps in the automatic recognition and calibration method could include first averaging an OCT image frame along the x—(i.e. angular) dimension. This selectively enhances the feature(s) which are rotationally stable in the y dimension (i.e. radius) instead of other image features generated by subject or surroundings. Efficacy of the method is improved if the image feature(s) used have high intensity relative to the surrounding pixels and if subject/environment features (noise) do not have strong circumferential symmetry.
Next, features can be detected within the image. Features can be found by any suitable method or algorithm such as, for example, peak searching, correlation, thresholding, or other pattern recognition algorithms known in the art. The efficacy of this method can be improved if the range over which uncontrolled path length displacements can occur is known a priori, thus limiting the required search space.
Then, the y-value(s) of feature(s) found in the previous step are compared to a pre-calibrated y-value which represents the actual physical location(s) of that feature(s) relative to the rotational axis, or to the location of a known “conjugate image” or “aliased image” of that feature(s) when using spectral-domain OCT.
Finally, the image is calibrated by shifting the OCT image pixels in the y dimension by the difference between searched feature(s) and pre-calibrated feature(s). Multiple features can be used to improved efficacy of the algorithm. After shifting the rectangular image in the y dimension, map to polar image coordinates. Radii measured to the center of the calibrated polar image will represent actual radii measured to the rotational axis in physical space.
Often image features due to the catheter are unwanted for effective and distraction-free display of the subject/environment features. For example, the catheter image features could overlap the subject/environment features.
Steps to remove (or make less noticeable) the image features could include cropping out the image feature(s) extent in the y, or radial, direction and in all columns, or angles. Removal (or diminishment) can further include calculating the average value of the pixels immediately inside and outside (above and below) of the cropped region for all columns, or angles, and inserting this averaged row, or circumference, in the cropped location.
As shown in
In some embodiments, calibration involves determining the position of ring 211 (or lines 217) in display 237 so that the system can calculate a calibration value based on a known position of calibration mark 215 (see
The system can additionally use a processor to perform an image processing operation to detect a feature such as sheath 211.
In some embodiments, the system detects a feature by starting with a defined area 227 (e.g., a pre-stored range of search values, or an arbitrary area, or the entire image).
The system searches for a feature (e.g., a sheath) within area 227 by performing a processing operation on the corresponding data. The processing operation can be any suitable search algorithm known in the art.
In some embodiments, a morphological image processing operation is used. Morphological image processing includes operations such as erosion, dilation, opening, and closing, as well as combination thereof. In some embodiments, these operations involve converting the image data to binary data giving each pixel a binary value. With pixels within area 227 converted to binary, each pixel of catheter sheath 211 will be black, and the background pixels will predominantly be white. In erosion, every pixel that is touching background is changed into a background pixel. In dilation, every background pixel that is adjacent to the non-background object pixels is changed into an object pixel. Opening is an erosion followed by a dilation, and closing is a dilation followed by an erosion. Morphological image processing is discussed in Smith, The Scientist and Engineer's Guide to Digital Signal Processing, 1997, California Technical Publishing, San Diego, Calif., pp. 436-442.
If sheath 211 is not found within area 227, area 227 can be increased and the increased area can be searched. This strategy can exploit the statistical properties of signal-to-noise ratio (SNR) by which the ability to detect an object is proportional to the square root of its area. See Smith, Ibid., pp. 432-436.
With continued reference to
While described above as detecting a reference item (e.g., catheter sheath 211) by receiving user input followed by using a processor to detect a location of the sheath, the steps can be performed in other orders. For example, the system can apply morphological processing operations to an entire image and detect every element, or every element that satisfies a certain quality criterion. In some embodiments, the system can receive user input that indicates a point within an image and the user can then choose the pre-detected element that is closest to that point within the image. Similarly, the steps can be performed simultaneously.
Using the methodologies herein, systems of the invention can detect an element within an image of an imaging system, such as an OCT system, with great precision, based on a search algorithm. Based on this detection, an actual location of a catheter sheath is determined and thus a precise calibration value for the catheter sheath and thus for the image (e.g., within a B-scan) is known. Where an expected z-coordinate Zc for the catheter sheath is known, based on information provided extrinsically, the a calibration value can be determined. For example, in
In some embodiments, the system calculates or uses the mean, median, or root-mean-squared distance of the sheath from the calibration mark to compute the calibration value. This may be advantageous in the event of interfering speckle noise, rough or acylindrical sheaths, non-uniform catheter rotation (NURD), angular displacement of a transducer within the sheath, off-center positioning of the transducer within the sheath, or a combination thereof. In certain embodiments, only a subset of the detected points are used, for example, for efficiency or performance optimization.
It will be appreciated that the foregoing description is applicable in live mode or review mode. If the imaging system is operating in live mode, capturing an image of tissue, the calibration can be put into effect either by changing the length of reference path 915 so that z-offset is zero or by transforming the dataset or on-screen image. The length of reference path 915 can be changed through the operation of the motor in the VDL. The distance Zc-Zs is converted into millimeters and the a command is sent to move the VDL to a new position.
If the dataset is to be transformed, either in live mode or while the system is operating in review mode, the system is digitally shifted, stretched, or a combination thereof.
In another aspect, the invention provides a method for calibrating an imaging system based on receipt of user input that indicates a “motion”, such as a click-and-drag operation on a computer screen.
While discussed above using a surface of a catheter sheath as a reference item which is used as a basis for calibration, other reference items are suitable. For example, any item that can be depicted such that its expected location and actual location can be compared in a display of an imaging system may be used. In some embodiments, a fiducial marker or calibration bar is introduced into the imaging target having a known dimension (e.g., 1 nm, 1 mm, 1 cm). The system operates to display a scale or a grid based on an expected appearance of the known dimension. The user then gives input indicating a point in the display near the reference item and the system also detects a location of the reference item in an area around the indicated point. Based on the expected and actual locations or dimensions of the reference item, a calibration value is calculated and a calibrated image is provided. User input, displays, and methods of receiving user input and performing calculations may be provided by one or more computers.
In certain embodiments, display 237 is rendered within a computer operating system environment, such as Windows, Mac OS, or Linux or within a display or GUI of a specialized system. Display 237 can include any standard controls associated with a display (e.g., within a windowing environment) including minimize and close buttons, scroll bars, menus, and window resizing controls. Elements of display 237 can be provided by an operating system, windows environment, application programming interface (API), web browser, program, or combination thereof (for example, in some embodiments a computer includes an operating system in which an independent program such as a web browser runs and the independent program supplies one or more of an API to render elements of a GUI). Display 237 can further include any controls or information related to viewing images (e.g., zoom, color controls, brightness/contrast) or handling files comprising three-dimensional image data (e.g., open, save, close, select, cut, delete, etc.). Further, display 237 can include controls (e.g., buttons, sliders, tabs, switches) related to operating a three dimensional image capture system (e.g., go, stop, pause, power up, power down).
In certain embodiments, display 237 includes controls related to three dimensional imaging systems that are operable with different imaging modalities. For example, display 237 may include start, stop, zoom, save, etc., buttons, and be rendered by a computer program that interoperates with OCT or IVUS modalities. Thus display 237 can display an image derived from a three-dimensional data set with or without regard to the imaging mode of the system.
A computer generally includes a processor for executing instructions and one or more memory devices for storing instructions, data, or both. Processors suitable for the execution of methods and operations described herein include, by way of example, both general and special purpose microprocessors (e.g., an Intel chip, an AMD chip, an FPGA). Generally, a processor will receive instructions or data from read-only memory, random access memory, or both. Generally, a computer will also include, or be operatively coupled, one or more mass storage devices for storing data that represent target such as bodily tissue. Any suitable computer-readable storage device may be used such as, for example, solid-state, magnetic, magneto-optical disks, or optical disks. Information carriers suitable for embodying computer program instructions and data include all forms of non-volatile memory, particularly tangible, non-transitory memory including by way of example semiconductor memory devices, (e.g., EPROM, EEPROM, NAND-based flash memory, solid state drive (SSD), and other flash memory devices); magnetic disks, (e.g., internal hard disks or removable disks); magneto-optical disks; and optical disks (e.g., CD and DVD disks).
INCORPORATION BY REFERENCEReferences and citations to other documents, such as patents, patent applications, patent publications, journals, books, papers, web contents, have been made throughout this disclosure. All such documents are hereby incorporated herein by reference in their entirety for all purposes.
EQUIVALENTSVarious modifications of the invention and many further embodiments thereof, in addition to those shown and described herein, will become apparent to those skilled in the art from the full contents of this document, including references to the scientific and patent literature cited herein. The subject matter herein contains important information, exemplification and guidance that can be adapted to the practice of this invention in its various embodiments and equivalents thereof.
Claims
1. A method of calibrating an imaging system, the method comprising:
- obtaining an image that includes a target and a reference item;
- detecting a location of the reference item within the image;
- determining a y-value of the location;
- comparing the determined y-value to a stored reference y-value;
- calculating a calibration value based on the comparison; and
- providing a calibrated image by shifting pixels in a y-direction according to the calibration value.
2. The method of claim 1, wherein the imaging system is an optical coherence tomography system.
3. The method of claim 1, wherein the reference item comprises an image of a catheter sheath.
4. The method of claim 1, wherein the detecting step comprises a morphological image processing operation.
5. The method of claim 1, further comprising digitally transforming image data to provide the calibrated image.
6. The method of claim 1, wherein detecting the location comprises one selected from the list consisting of peak searching, correlation, thresholding, and performing a pattern recognition algorithm.
7. The method of claim 1, further comprising obtaining an a priori range over which a path length displacement may occur.
8. The method of claim 1, wherein determining the y-value comprises averaging across the image.
9. The method of claim 1, wherein calculating the calibration value comprises taking a difference between the determined y-value and the stored reference y-value.
10. The method of claim 9, wherein the pixels are shifted by the difference.
11. The method of claim 1, further comprising removing a feature from the image.
12. The method of claim 1, further comprising removing the reference item from the image.
13. The method of claim 11, wherein removing comprises cropping out the feature across its extent in the y direction at all x positions.
14. The method of claim 13, wherein removing comprises calculating an average value of pixels proximal to the cropped region and inserting the average values to replace the cropped pixels.
15. The method claim 1, wherein the method is performed using a computer system comprising a tangible, non-transitory memory coupled to a processor.
16. An imaging system, the system comprising:
- a processor coupled to a non-transitory memory, wherein the system is operable to: obtain an image that includes a target and a reference item; detect a location of the reference item within the image; determine a y-value of the location; compare the determined y-value to a stored reference y-value; calculate a calibration value based on the comparison; and provide a calibrated image by shifting pixels in a y-direction per the calibration value.
Type: Application
Filed: Mar 12, 2014
Publication Date: Sep 18, 2014
Patent Grant number: 9301687
Applicant: VOLCANO CORPORATION (San Diego, CA)
Inventor: Nathaniel J. Kemp (Concord, MA)
Application Number: 14/206,485
International Classification: G06T 7/00 (20060101);