Classifications

• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
  • A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
  • A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
• • A—HUMAN NECESSITIES
  • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
  • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
  • A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
  • A61B3/0016—Operational features thereof
  • A61B3/0025—Operational features thereof characterised by electronic signal processing, e.g. eye models
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
  • G01B9/00—Measuring instruments characterised by the use of optical techniques
  • G01B9/02—Interferometers
  • G01B9/02083—Interferometers characterised by particular signal processing and presentation
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
  • G01B9/00—Measuring instruments characterised by the use of optical techniques
  • G01B9/02—Interferometers
  • G01B9/0209—Low-coherence interferometers
  • G01B9/02091—Tomographic interferometers, e.g. based on optical coherence
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T11/00—2D [Two Dimensional] image generation
  • G06T11/003—Reconstruction from projections, e.g. tomography
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T19/00—Manipulating 3D models or images for computer graphics
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T7/00—Image analysis
  • G06T7/20—Analysis of motion
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T7/00—Image analysis
  • G06T7/20—Analysis of motion
  • G06T7/215—Motion-based segmentation
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2207/00—Indexing scheme for image analysis or image enhancement
  • G06T2207/10—Image acquisition modality
  • G06T2207/10072—Tomographic images
  • G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2207/00—Indexing scheme for image analysis or image enhancement
  • G06T2207/30—Subject of image; Context of image processing
  • G06T2207/30004—Biomedical image processing
  • G06T2207/30041—Eye; Retina; Ophthalmic
• • G—PHYSICS
  • G06—COMPUTING OR CALCULATING; COUNTING
  • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
  • G06T2210/00—Indexing scheme for image generation or computer graphics
  • G06T2210/41—Medical

Definitions

• the present invention relates to medical imaging, and in particular to analysis and visualization methods for data acquired through optical coherence tomography.
• Optical coherence tomography is a noninvasive, noncontact imaging modality that uses coherence gating to obtain high-resolution cross-sectional images of tissue microstructure.
• FD-OCT Fourier domain OCT
• the interferometric signal between light from a reference and the back-scattered light from a sample point is recorded in the frequency domain rather than the time domain.
• a one-dimensional Fourier transform is taken to obtain an A-line spatial distribution of the object scattering potential.
• the spectral information discrimination in FD-OCT can be accomplished by using a dispersive spectrometer in the detection arm in the case of spectral-domain OCT (SD-OCT) or rapidly tuning a swept laser source in the case of swept-source OCT (SS-OCT).
• SD-OCT spectral-domain OCT
• SS-OCT swept-source OCT
• OCT Angiography to map the retinal vasculature or identify regions with flow in the tissue
• OCT Angiography provides correlation between microstructure and microvasculature of optic nerve head in human subjects
• J. Biomed. Opt. 17, 116018 (2012) Zhao et al, "Doppler standard deviation imaging for clinical monitoring of in vivo human skin blood flow”
• Optics Letters 25, 1358-1360 (2000) Fingler et al.
• Mobility and transverse flow visualization using phase variance contrast with spectral domain optical coherence tomography Optics Express. Vol. 15, No. 20.
• OCT Angiography provides a non-invasive technique to visualize and indirectly quantify the integrity of retinal circulation pathways.
• Anomalies in retinal circulation have a direct relation to ocular pathologies, especially within the macula, wherein compromised hemodynamics may not only be related to decreased visual acuity, but could also be a surrogate biomarker for ocular pathologies like retinal vein occlusion (RVO), diabetic retinopathy (DR), and intra retinal microvasculature abnormality (IRMA).
• RVO retinal vein occlusion
• DR diabetic retinopathy
• IRMA intra retinal microvasculature abnormality
• correlation between retinal vasculature and blood flow are attributes of interest in a number of ocular defects.
• DR and RVO are pathologies that could lead to early changes to the vascular structure and function, and may, in turn, be etiologic to numerous complications like macular edema, retinal ischemia and optic neuropathy. For these cases, quantification and
• ischemic regions in the retina can be mapped to evaluate the extent of damage and further management of the disease.
• vascular rich retina there is a small area in the macula, at the fovea, which is devoid of any capillaries. This is called the Foveal Avascular zone (FAZ), and abnormal changes in the size of this region are also indicative of pathologies like ischemic maculopathy and DR.
• FAZ Foveal Avascular zone
• the quantification of the FAZ and measuring changes in its size over time can be a clinically significant numerical score for disease presence and progression, especially for DR.
• Angiography No information of the depth structure of the vasculature is captured by this method.
• vascular images generated by examining the OCT intensity or phase signal are non-invasive, and provide comparable fidelity in capturing the existing vascular network with blood flow contrast along with its depth encoding.
• vascular and the typically avascular retina are important indicators of developing retinal pathologies.
• visualization of the vascular structure helps in boosting the diagnostic efficacy of this imaging technique, it can be further leveraged by augmenting the visualization with some salient quantifications and metrics derived from the identified vascular and avascular sections of the retina.
• the primary quantity of interest is the global or structure-specific retinal blood flow kinetics, which can be challenging to quantify because of low flow velocities relative to the temporal resolution of the technique, and the almost perpendicular orientation of the capillaries with respect to the probing beam.
• derived quantifiers from the angiography data which serve to aid in differentiating capillary networks in healthy and diseased eyes are also desirable.
• Avakian et al. demonstrated the use of fractal characterization of fluorescein angiography (FA) images of the human retina to distinguish between healthy and diseased retina (see for example Avakian, et al., "Fractal analysis of region-based vascular change in the normal and non-proliferative diabetic retina,” Curr. Eye Res. 24, 274-280, 2002). Schmoll et al.
• FA fluorescein angiography
• OCT angiography data is often displayed as 2D projections with the color encoded depth information (see Kim et al. "In vivo volumetric imaging of human retinal circulation with phase variance OCT,” Biomedical Optics Express, 2(6), 1504-1513 (2011)).
• 2D projections at least allow distinguishing capillary layers of different depths. They however lack the 3D impression and also don't provide easily accessible information of which larger retinal vessels feed and drain different capillary network regions.
• Retinal vessel connectivity measures are also known for fundus photography, they however only focus on a few major retinal vessels in 2D fundus images, rather than visualizing the supply of dense, complex parafoveal capillary networks (see for example Al-Diri et al.
• Ganesan et al. investigates the connectivity of vessels in mouse retinas from the largest vessels to the smallest capillaries in confocal microscopy images in order to develop a network model (see for example Ganesan et al. "Development of an Image-Based Network Model of Retinal Vasculature,” Annals of Biomedical Engineering 38(4) 1566-1585 (2010)). They however don't describe using this as a way to interactively visualize human angiography acquisitions.
• the method described herein is a non-invasive, computational technique to generate images of retinal vasculature (or blood flow) that are then used to either extract various
• the anatomical location of the vasculature is defined as locations where there is an appreciable motion contrast, which is typically due to flow of blood.
• OCT based methods to detect motion contrast such as Doppler OCT, speckle or intensity variance, and phase-resolved methods.
• the motion contrast is determined by obtaining at least two OCT measurements at approximately the same location, where the two measurements are separated in time by a pre-determined interval and by applying an algorithm to look at the changes in the complex OCT signal or its components such as intensity or phase.
• the accuracy of these measurements can be improved by minimizing motion related errors including but not limited to removing signal due to bulk motion of the sample in the axial direction.
• the vascular structure extracted is post- processed to remove outliers and smooth the vessel structure.
• This derived structure is then depth coded and displayed over the rendered retinal anatomy (magnitude image).
• the post-processed vascular structure within a specific depth range can be summed or integrated along the axial direction to generate a projection map that shows the vasculature as an en face view, devoid of any depth information.
• the proposed invention deviates from the known prior art by detailing a completely automated (no manual intervention) method to accurately determine the capillary devoid regions of retina, by examination of the statistical properties of the intensity content of the retinal image which preferentially contrasts vascular regions.
• avascular zone In addition to the avascular zone, other metrics, like the vasculature density, capillary density, vessel geometry, capillary diameter, inter-capillary distance, area bounded by capillary loops, etc. can be determined by standard mathematical models and tools. These metrics can be identified in the vicinity of the fovea, or in other areas of interest, such as the perifoveal or peripapillary regions, the papillomacular bundle, or within the optic nerve head. This technique could be further extended to automatically identify regions of retinal ischemia in pathologies such as branch retinal vein occlusion (BRVO) and central retinal vein occlusion (CRVO). Furthermore, these techniques could also assist in identifying intraretinal microvascular abnormalities (IRMA).
• BRVO branch retinal vein occlusion
• CRVO central retinal vein occlusion
• IRMA intraretinal microvascular abnormalities
• IRMA is typically a DR-related condition that results in areas of capillary dilatation and intraretinal formation of new capillary beds. Often the IRMA related new vessel formation occurs in retinal tissues to act as shunts through areas of nonperfusion. or ischemia. Change in the metrics defined above could be used as a criterion to monitor if there has been new growth of vasculature or change in the regions of nonperfusion or ischemia.
• a novel method for effectively visualizing OCT angiography acquisitions in a meaningful way and quantitatively characterizing vasculature networks is presented.
• the examiner could select a vessel within an OCT angiography acquisition and the program would show all connecting vessels down to the capillary network.
• the information about the connectivity of different retinal vessels may also be used to quantitatively evaluate OCT angiography acquisitions and compare them to a normative database.
• FIG. 1 shows a flow chart of various step involved with processing motion contrast OCT data according to the present invention.
• FIG. 2 illustrates a generalized ophthalmic OCT imaging system that could be used for collection of motion contrast data.
• FIG. 3 shows an en face vasculature image generated from OCT data using normalized vector difference variance.
• FIG. 4 shows an OCT image of the retina illustrating three different plexus and the layers that they include.
• FIG. 5a shows a 3D visualization of OCT angiography data
• FIG. 5b shows a
• FIG. 6 illustrates how a stereoscopic image pair could be generated from motion contrast data to enable a type of 3D visualization.
• FIG. 7 shows an intensity histogram of an en face vasculature image that can be used to identify the foveal avascular zone (FAZ).
• FAZ foveal avascular zone
• FIG. 8 shows the results of an isophote delineation of an en face vasculature image according to one aspect of the present invention.
• FIG. 9 illustrates how the FAZ can be isolated from the rest of the image data after the isophote delineation in FIG. 8.
• FIG. 10 shows a map of the density of vessels in radially distributed sectors around the fovea. DETAILED DESCRIPTION
• FIG. 1 Preferred and alternative embodiments for the processing of vasculature enhanced OCT data are illustrated in the schematic of FIG. 1.
• This figure illustrates three possible ways to generate and use vasculature data from complex OCT data 102 acquired and reconstructed in an OCT system 101.
• the text below refers to each possible combination as a "Process", A, B or C, depending on which process path has been taken to generate the input data for the embodiment under discussion.
• different portions of the complex data intensity only, phase only, or both intensity and phase
• undergo different processing steps layer segmentation, motion contrast, integration/summation for en face image generation
• OCT data can be collected with any type of OCT system employing a variety of scan patterns, for example, a spectral domain OCT system, or a swept source OCT system, employing laser sources of different wavelength like 840 nm or 1060 nm.
• a diagram of a generalized OCT system is shown in FIG. 2.
• Light from source 201 is routed, typically by optical fiber 205, to illuminate the sample 210, a typical sample being tissues in the human eye.
• the source 201 can be either a broadband light source with short temporal coherence length in the case of SD-OCT or a wavelength tunable laser source in the case of SS-OCT.
• the light is scanned, typically with a scanner 207 between the output of the fiber and the sample, so that the beam of light (dashed line 208) is scanned laterally (in x and y) over the area or volume to be imaged.
• Light scattered from the sample is collected, typically into the same fiber 205 used to route the light for sample illumination.
• Reference light derived from the same source 201 travels a separate path, in this case involving fiber 203 and retro- reflector 204 with an adjustable optical delay.
• a transmissive reference path can also be used and that the adjustable delay could be placed in the sample or reference arm of the interferometer.
• Collected sample light is combined with reference light, typically in a fiber coupler 202, to form light interference in a detector 220.
• the output from the detector is supplied to a processor 221.
• the results can be stored in the processor 221 or displayed on display 222.
• the processing and storing functions may be localized within the OCT instrument or functions may be performed on an external processing unit to which the collected data is transferred. This unit could be dedicated to data processing or perform other tasks which are quite general and not dedicated to the OCT device.
• the sample and reference arms in the interferometer could consist of bulk-optics, fiber-optics or hybrid bulk-optic systems and could have different architectures such as Michelson, Mach- Zehnder or common-path based designs as would be known by those skilled in the art.
• Light beam as used herein should be interpreted as any carefully directed light path.
• the reference arm needs to have a tunable optical delay to generate interference.
• Balanced detection systems are typically used in TD-OCT and SS-OCT systems, while spectrometers are used at the detection port for SD-OCT systems.
• the invention described herein could be applied to any type of OCT system capable of generating data for functional analysis.
• the interference causes the intensity of the interfered light to vary across the spectrum.
• the Fourier transform of the interference light reveals the profile of scattering intensities at different path lengths, and therefore scattering as a function of depth (z-direction) in the sample (see for example Leitgeb et al. "Ultrahigh resolution Fourier domain optical coherence tomography,” Optics Express 12(10):2156 (2004)).
• the Fourier transform results in complex data, and the absolute values of the complex data are tabulated to construct the intensity image.
• the complex OCT signal also encodes information related to the phase shifts arising from local sample motion, and can be used to deduce quantities related to physical motion of dominant scatterers in the sample with high sensitivity.
• A-scan The profile of scattering as a function of depth is called an axial scan (A-scan).
• a set of A-scans measured at neighboring locations in the sample produces a cross-sectional image (tomogram or B-scan) of the sample.
• B-scan cross-sectional image
• a collection of B-scans collected at different transverse locations on the sample makes up a data volume or cube.
• fast axis refers to the scan direction along a single B-scan whereas slow axis refers to the axis along which multiple B-scans are collected.
• Intensity based local searches or global approaches can be used on the magnitude images to extract prominent layers which are then used as boundaries for the summation of intensities.
• the boundaries extracted serve to include only that tissue extent in the summation which is known a-priori to have blood vessels.
• the result of this summing procedure is a flat view (projection) of the volume looking along (and into) the imaging axis, and the features in this projection (en face vasculature image) capture the vascular distribution (Process A in FIG. 1). Since the summation integrates out the depth along the axial direction, this view only captures the omnibus morphology of the vasculature, and not position in the depth direction.
• FIG. 3 shows the en face vasculature image generated by selective summation (FIG. 1, Process A) through the motion contrasted volume generated by performing normalized vector difference variance as described in US Patent Publication No. 2012/0277579, hereby incorporated by reference.
• the bright vasculature indicated by the arrows 301 and 302 stands out against the dark background.
• the absence of vessels in the central region indicates the foveal avascular zone (FAZ).
• FAZ foveal avascular zone
• FIG. 4 Anatomically the vasculature is distributed in the 3 major sections of the retina (FIG. 4): the superficial capillary plexus (SCP), intermediate capillary plexus (ICP) and the deep capillary plexus (DCP) (see for example Kim et al. "Noninvasive Imaging of the Foveal Avascular Zone with High-Speed, Phase- Variance Optical Coherence Tomography” Investigative Ophthalmology & Visual Science, 53 (1), 85 - 92 (2012) hereby incorporated by reference).
• SCP superficial capillary plexus
• ICP intermediate capillary plexus
• DCP deep capillary plexus
• Each of these plexus is made up of a finite number of retinal layers, which, in an alternate embodiment, can be used as bounding layers for selective summation to generate vasculature en face views showing the three specific types of vasculature networks.
• Selective summation is not required to generate an en face image and other examples of selective summations or other ways to represent a plurality of intensity values as a single representative value (e.g. integration, summing, minimum, maximum, median value, etc.) can be envisioned by those skilled in the art (see for example US Patent No. 7,301,644, US Patent Publication No. 2011/0034803 and US Patent Publication No. 2008/0100612 hereby incorporated by reference).
• the projection views are devoid of any depth information, because the hyperintense pixels, which signal the presence of vasculature are summed along the axial direction.
• the retinal vasculature is distributed in the 3 -dimensional space of the retinal tissue, with the distribution and characteristics of the vasculature varying also within each specific plexus (FIG. 4).
• the depth information in the view presented in FIG. 3 can be preserved by explicitly using the volumetric definition of the vasculature (after motion contrasting the phase data), and rendering the volume in three-dimensions (3D, process B in FIG. 1).
• this process which extends process B, would entail preserving the high intensity locations in the acquired volume using a thresholding or selection criteria, and performing some post-processing to enforce the anatomical connectivity of the vasculature.
• a schematic of the expected form of the 3D vasculature map is shown in FIG. 5a, where the black lines trace out the vessel path in 3D space.
• the corresponding 2D projection view is illustrated in FIG. 5b.
• visualization mode will provide a visual representation of the vascular architecture / distribution in the space of the imaged retina.
• the data included in the 3D representation could be limited to a particular plexus, using methods described above, to limit the data displayed to a particular location in the retina.
• the location in the retina can be defined by adjacency to a particular location, such as the centroid or a boundary identified by segmentation, or the location may be limited to be between two such boundaries. In this way the 3D vessel model of each plexus could be independently reviewed.
• the examiner could select a vessel within a volume using a data input device such as a mouse or touch screen interface and the program would then highlight the connecting vessels down to the capillary level.
• a data input device such as a mouse or touch screen interface
• Such visualization may improve the identification of blockages or leakages.
• Visualizing connected vasculature could be done, e.g. by only displaying the specific connected vasculature or by only highlighting it within the volume in order to contrast it from the other vasculature.
• Visualization of connected vasculature could also involve an image series or a movie, where the movie starts with only the initially selected vessel, which then grows until all the connected vessels are shown.
• the speed at which connecting vessels are added could be normalized by their vessel diameter in order to mimic the propagation speed of the blood within the network.
• the information created by the vessel segmentation algorithms may serve as additional quantitative parameters, which could be used for comparing data sets with a normative data base or for tracking changes in a particular patient over time.
• Such parameters may be total vessel length, number of bifurcations, capillary density of arteriole vs. venous capillaries, vessel diameter parent / daughter vessels, bifurcation angles, vascular tortuosity, capillary network volume vs. static tissue volume.
• the actual depths of the vascular locations in the retina can be used as disparity maps to generate a stereoscopic image pair of the vascular network that can be viewed by the clinician to get a better idea of the vascular distribution in space.
• a preferred scale is selected to map the range of possible distances (for example, from 0 to 2 mm), after which a left and right volume pair can be generated and fused into a stereoscopic pair as illustrated in FIG. 6. This can be rendered either through a special pair of 3D glasses, or via a 3D display technology.
• Extensions to this alternative embodiment also allows for augmenting this stereoscopic view by the retinal anatomy, depth encoded in the magnitude image to generate vascular maps in relation to the various layers of the retina.
• the vessels derived from the three plexus can be color coded preferentially to generate a more informative image fusion approach to visualization.
• the preferred embodiment employs the use of vasculature en face image to derive clinically significant information, the most important of which is the automated detection of the foveal avascular zone, or FAZ.
• the detection of the FAZ used the output of process A (FIG. 1), and can be accomplished by creating an intensity histogram of the vasculature en face image.
• the intensity histogram (FIG. 7) is a very popular graphing technique in image processing to generate a frequency distribution (represented along the ordinate of the graph) of the intensity values in the image which are represented along the abscissa.
• Many software systems contain standard and optimized libraries for the generation of the image histogram.
• the preferred embodiment starts the detection process by detecting the peaks of the histogram. As illustrated in FIG.
• histogram of the vascular en face image shows a characteristic bi-modal distribution, and, the transition between the FAZ and the rest of the map is located between these peaks A and B.
• the en- face image is interrogated for pixels which have this specific intensity value. Lines of similar intensity are called isophotes, and the trace of the isophote at intensity level i provides the preferred embodiment with a discriminating contour which maximally contains the FAZ on the inside, and the vascular retina on the outside (FIG. 8).
• the isophotes are generally delineated by a contour operation which entails thresholding the intensity levels on the image with a small neighborhood of the required isophote at level I as indicated by the arrow in FIG. 8.
• the results of the isophote delineation are shown in FIG. 8. Because of the noisy nature of the signal, some small isolated contours can also be seen, but the algorithm employed in the preferred embodiment ensures that the isophote that is convex and has the largest perimeter will always enclose the FAZ.
• the preferred embodiment reports multiple matches (small isolated contours in FIG. 8) as a list of a list of x- and y- coordinates of each contour. These coordinates can be used in a straightforward way to calculate the perimeter of each contour.
• the contour that best delineates the FAZ is the one with the highest perimeter.
• the region inside the FAZ contour can be isolated from the rest of the contours as illustrated in FIG. 9.
• morphological parameters of the resulting FAZ region can be calculated to quantify the FAZ on a case by case basis.
• the area is calculated as counting the number of "on" pixels, scaled by the area of each pixel, and is 5675 pixel units.
• the delineated FAZ has an eccentricity of 0.71 with a major axis of 104.7 pixel units (523.5 microns) and a minor axis of 73.1 pixel units (548.5 microns).
• the remainder of the vascular map (exterior to the delineated FAZ) can be used to calculate capillary density, especially as a function of distance from the center of the fovea, and generate a map of the density in radially distributed sectors around the fovea as illustrated in FIG. 10.
• radial sectors can be analyzed and the ratio of "on” to "off pixels can be used to generate a coarse sectorial vessel density map.
• the definition of the "on” and "off pixels can be established by a judicious threshold selection.
• an alternative embodiment (process C, FIG. 1) allows the use of prior knowledge of the fovea from the magnitude image using a fovea detection algorithm (see for example US Patent No. 8,079,711 hereby incorporated by reference) to create a localized histogram and hence confine the search for the vascular-avascular transition zone to a smaller area.
• a fovea detection algorithm see for example US Patent No. 8,079,711 hereby incorporated by reference
• the FAZ detection by the preferred embodiment on images generated by process A will most likely fail to qualify the maximum perimeter criteria.
• Such a case will be flagged by a null (zero) FAZ detection (i.e., a blank image in FIG. 9) after morphological cleaning steps.
• Other clinically significant quantifiable characteristics or metrics that can be derived from the en face vasculature image (from one or more plexus) or from a 3D OCT angiography data volume include total capillary volume (sum of all pixels above a given threshold intensity), as well as metrics that could be derived from the blood vessel patterns, such as tortuosity, regularity, segment length, total crossings, the number of bifurcations, vessel width parameters, ratio of small to large vessels, capillary density, capillary density ratio between arteriole and venous capillaries, capillary diameter, inter-capillary distance, area bounded by capillary loops, etc.
• a metric can be compared to a database of normal eyes or eyes with a known pathology to diagnose or to track progression of a particular disease or condition.
• Another way to use this information is to use the foveal avascular zone as a region over which to evaluate other features such as layer thickness or other parameters derived from the OCT image, or other registered images.
• the thickness of the photoreceptor layer within the FAZ should be specifically related to cones rather than rods.
• Registration to images that contain other information about the photoreceptors, such as adaptive optics images, might allow quantification of multiple metrics that affect the foveal region.
• the preferred embodiment described above refers to visualization of retinal capillary vessels and metrics associated with these. Similar methods could be applied to the choroidal vasculature, including the choriocapillaris, Sattler's layer and Haller's layer. Although various applications and embodiments that incorporate the teachings of the present invention have been shown and described in detail herein, those skilled in the art can readily devise other varied embodiments that still incorporate these teachings. Although the description of the present invention is discussed herein with respect to the sample being a human eye, the applications of this invention are not limited to eye and can be applied to any application using OCT.
• Leitgeb "Imaging of the parafoveal capillary network and its integrity analysis using fractal dimension,” Biomed. Opt. Express 2, 1159-1168 (2011).
• Leitgeb et al. "Real-time assessment of retinal blood flow with ultrafast acquisition by color Doppler FDOCT,” Optics Express, 11, 3116-3121 (2003).

Landscapes

• Engineering & Computer Science (AREA)
• Health & Medical Sciences (AREA)
• Physics & Mathematics (AREA)
• General Physics & Mathematics (AREA)
• Life Sciences & Earth Sciences (AREA)
• Theoretical Computer Science (AREA)
• General Health & Medical Sciences (AREA)
• Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
• Radiology & Medical Imaging (AREA)
• Medical Informatics (AREA)
• Veterinary Medicine (AREA)
• Ophthalmology & Optometry (AREA)
• Biomedical Technology (AREA)
• Heart & Thoracic Surgery (AREA)
• Molecular Biology (AREA)
• Surgery (AREA)
• Animal Behavior & Ethology (AREA)
• Biophysics (AREA)
• Public Health (AREA)
• Signal Processing (AREA)
• Multimedia (AREA)
• Computer Vision & Pattern Recognition (AREA)
• Computer Graphics (AREA)
• Computer Hardware Design (AREA)
• General Engineering & Computer Science (AREA)
• Software Systems (AREA)
• Eye Examination Apparatus (AREA)
• Apparatus For Radiation Diagnosis (AREA)
• Magnetic Resonance Imaging Apparatus (AREA)

Priority Applications (5)

Applications Claiming Priority (6)

Publications (1)

Family

ID=48470922

Family Applications (1)

Country Status (6)

Cited By (13)

Families Citing this family (80)

Citations (13)

Family Cites Families (60)

• 2013
  • 2013-02-28 US US13/781,375 patent/US9357916B2/en active Active
  • 2013-05-08 EP EP13724195.6A patent/EP2852317B1/en active Active
  • 2013-05-08 WO PCT/EP2013/059560 patent/WO2013167641A1/en not_active Ceased
  • 2013-05-08 JP JP2015510811A patent/JP6371762B2/ja active Active
  • 2013-05-08 CN CN201380023147.6A patent/CN104271031B/zh active Active
  • 2013-05-08 EP EP22185690.9A patent/EP4122377A1/en active Pending
  • 2013-05-08 IN IN8852DEN2014 patent/IN2014DN08852A/en unknown
• 2016
  • 2016-05-05 US US15/147,402 patent/US20160317029A1/en not_active Abandoned
• 2018
  • 2018-07-13 JP JP2018132927A patent/JP2018175888A/ja active Pending

Patent Citations (14)

Non-Patent Citations (27)

Also Published As

Similar Documents

Legal Events