US20150211838A1 - Method and apparatus for ultrahigh sensitive optical microangiography - Google Patents
Method and apparatus for ultrahigh sensitive optical microangiography Download PDFInfo
- Publication number
- US20150211838A1 US20150211838A1 US14/658,558 US201514658558A US2015211838A1 US 20150211838 A1 US20150211838 A1 US 20150211838A1 US 201514658558 A US201514658558 A US 201514658558A US 2015211838 A1 US2015211838 A1 US 2015211838A1
- Authority
- US
- United States
- Prior art keywords
- scans
- fast
- slow
- scan
- sample
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Abandoned
Links
- 230000003287 optical effect Effects 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims description 38
- 238000003384 imaging method Methods 0.000 claims abstract description 85
- 230000017531 blood circulation Effects 0.000 claims abstract description 34
- 230000003595 spectral effect Effects 0.000 claims abstract description 21
- 239000000523 sample Substances 0.000 claims description 51
- 210000004204 blood vessel Anatomy 0.000 claims description 16
- 238000012014 optical coherence tomography Methods 0.000 claims description 10
- 238000001727 in vivo Methods 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 6
- 238000011503 in vivo imaging Methods 0.000 claims description 4
- 238000002583 angiography Methods 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 abstract description 10
- 210000004207 dermis Anatomy 0.000 description 21
- 210000001525 retina Anatomy 0.000 description 19
- 210000003583 retinal pigment epithelium Anatomy 0.000 description 18
- 241000282414 Homo sapiens Species 0.000 description 15
- 239000002245 particle Substances 0.000 description 15
- 210000003161 choroid Anatomy 0.000 description 14
- 230000033001 locomotion Effects 0.000 description 13
- 210000001519 tissue Anatomy 0.000 description 12
- 210000003491 skin Anatomy 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 8
- 230000005653 Brownian motion process Effects 0.000 description 7
- 238000005537 brownian motion Methods 0.000 description 7
- 238000001514 detection method Methods 0.000 description 7
- 230000010287 polarization Effects 0.000 description 7
- 230000000875 corresponding effect Effects 0.000 description 6
- 238000013459 approach Methods 0.000 description 5
- 210000000601 blood cell Anatomy 0.000 description 5
- 230000008822 capillary blood flow Effects 0.000 description 5
- 230000001575 pathological effect Effects 0.000 description 5
- 230000003068 static effect Effects 0.000 description 5
- 239000002131 composite material Substances 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 230000004089 microcirculation Effects 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 210000001367 artery Anatomy 0.000 description 3
- 230000004069 differentiation Effects 0.000 description 3
- 210000002615 epidermis Anatomy 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 206010033675 panniculitis Diseases 0.000 description 3
- 230000010412 perfusion Effects 0.000 description 3
- 238000004366 reverse phase liquid chromatography Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- 208000012661 Dyskinesia Diseases 0.000 description 2
- GWEVSGVZZGPLCZ-UHFFFAOYSA-N Titan oxide Chemical compound O=[Ti]=O GWEVSGVZZGPLCZ-UHFFFAOYSA-N 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000000877 morphologic effect Effects 0.000 description 2
- 230000017311 musculoskeletal movement, spinal reflex action Effects 0.000 description 2
- 210000003733 optic disk Anatomy 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 230000000284 resting effect Effects 0.000 description 2
- 230000002207 retinal effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000002792 vascular Effects 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 240000008100 Brassica rapa Species 0.000 description 1
- 206010012667 Diabetic glaucoma Diseases 0.000 description 1
- 206010012689 Diabetic retinopathy Diseases 0.000 description 1
- 208000003098 Ganglion Cysts Diseases 0.000 description 1
- 108010010803 Gelatin Proteins 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 206010067193 Naevus flammeus Diseases 0.000 description 1
- 208000006787 Port-Wine Stain Diseases 0.000 description 1
- 208000000453 Skin Neoplasms Diseases 0.000 description 1
- 208000005400 Synovial Cyst Diseases 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 206010064930 age-related macular degeneration Diseases 0.000 description 1
- 210000002565 arteriole Anatomy 0.000 description 1
- 238000010009 beating Methods 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 230000008081 blood perfusion Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000002316 cosmetic surgery Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 206010012601 diabetes mellitus Diseases 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000003743 erythrocyte Anatomy 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 208000002026 familial multiple nevi flammei Diseases 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- 239000008273 gelatin Substances 0.000 description 1
- 229920000159 gelatin Polymers 0.000 description 1
- 235000019322 gelatine Nutrition 0.000 description 1
- 235000011852 gelatine desserts Nutrition 0.000 description 1
- 230000004886 head movement Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 208000002780 macular degeneration Diseases 0.000 description 1
- 239000008267 milk Substances 0.000 description 1
- 210000004080 milk Anatomy 0.000 description 1
- 235000013336 milk Nutrition 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009448 modified atmosphere packaging Methods 0.000 description 1
- 235000019837 monoammonium phosphate Nutrition 0.000 description 1
- 210000004126 nerve fiber Anatomy 0.000 description 1
- 210000001328 optic nerve Anatomy 0.000 description 1
- 108091008695 photoreceptors Proteins 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000004286 retinal pathology Effects 0.000 description 1
- 210000001210 retinal vessel Anatomy 0.000 description 1
- 210000003786 sclera Anatomy 0.000 description 1
- 201000000849 skin cancer Diseases 0.000 description 1
- 210000000106 sweat gland Anatomy 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 210000000264 venule Anatomy 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- 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
-
- 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
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/026—Measuring blood flow
- A61B5/0261—Measuring blood flow using optical means, e.g. infrared light
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/44—Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
- A61B5/441—Skin evaluation, e.g. for skin disorder diagnosis
- A61B5/444—Evaluating skin marks, e.g. mole, nevi, tumour, scar
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
- A61B5/0062—Arrangements for scanning
- A61B5/0066—Optical coherence imaging
Definitions
- Embodiments herein relate to the field of imaging, and, more specifically, to a method and apparatus for ultrahigh sensitive optical microangiography.
- the assessment of blood flow in living tissue provides important information for diagnostics, treatment and/or management of pathological conditions.
- the assessment of cutaneous (skin) microcirculations may provide important information for pathological conditions in dermatology, such as skin cancer, port wine stain treatment, diabetes, and plastic surgery.
- assessment of the ocular perfusion within the retina and choroid of the human eye is important in the diagnosis, treatment, and management of a number of pathological conditions in ophthalmology, such as age-related macular degeneration, diabetic retinopathy, and glaucoma.
- 3D three dimensional
- FIG. 1 illustrates a functional block diagram of an imaging apparatus in accordance with various embodiments of the present invention
- FIG. 2 illustrates an example of a 3D data set, in accordance with various embodiments
- FIG. 3 illustrates an example of drive signals for an x-scanner and y-scanner, in accordance with various embodiments
- FIG. 4 illustrates another example of drive signals for an x-scanner and y-scanner, in accordance with various embodiments
- FIGS. 5A-C show images of capillaries within human dermis, taken with a ultrahigh sensitive optical microangiography (UHS-OMAG) system according to various embodiments: (A) a Fourier domain optical coherence tomography (FDOCT) structural image; (B) a UHS-OMAG flow image; and (C) a phase resolved optical Doppler tomography (PRODT) cross-sectional flow image from the UHS-OMAG image of FIG. 5B ;
- FDOCT Fourier domain optical coherence tomography
- PRODT phase resolved optical Doppler tomography
- FIGS. 5D-F show images of capillaries within human dermis, taken with a prior optical microangiography (OMAG) system, for comparison with the images in FIGS. 5A-C :
- D a Fourier domain optical coherence tomography (FDOCT) structural image
- E an OMAG flow image
- F a phase resolved optical Doppler tomography
- FIG. 6A is a B-scan structural image of a flow phantom, and Hg. 6 B is a corresponding UHS-OMAG flow image, in accordance with various embodiments.
- FIG. 7A is a velocity image obtained from a flow phantom assessed by UHS-OMAG, in accordance with various embodiments.
- FIG. 7B is a plot of the velocity data across the capillary tube at the position shown as the horizontal line in FIG. 7A , in accordance with various embodiments;
- FIG. 7C and FIG. 7D are the corresponding results to FIGS. 7A and 7B , respectively, obtained by PRODT imaging of the same phantom, in accordance with various embodiments;
- FIG. 8A is a schematic diagram of the blood vessel system of the human skin
- FIG. 8B is a photograph of a palm showing the scanning area
- FIG. 8C is a three dimensional (3D) rendered OMAG image of blood vessels together with the 3D structures, in accordance with various embodiments;
- FIG. 8D is a cross sectional view of the imaged blood vessels, in accordance with various embodiments.
- FIGS. 9A-D show UHS-OMAG detailed projection views of microcirculation network at different depths of skin obtained from: (A) 400-450 ⁇ m (closely representing papillary dermis); (B) 450-650 ⁇ m; (C) 650-780 ⁇ m (closely representing reticular dermis); and (D) 780-1100 ⁇ m (part of hypodermis), respectively, in accordance with various embodiments;
- FIG. 10 shows a schematic diagram of a UHS-OMAG system, in accordance with various embodiments.
- FIGS. 11A-B show in vivo UHS-OMAG imaging, in accordance with various embodiments, of the posterior segment of eye near the macular region towards the optic nerve head: (A) an OMAG B-scan of microstructures showing morphological features, and (B) the corresponding OMAG blood flow image;
- FIGS. 11C-D show projection maps of blood flow distribution within: (C) retina and (D) choroid, obtained from one 3D scan of an area of ⁇ 3 ⁇ 3 mm 2 near the macular region, in accordance with various embodiments;
- FIGS. 12A-F show depth-resolved images, in accordance with various embodiments, of patent blood vessels within the retina ( FIGS. 12A-C ) and choroid ( FIGS. 12D-F ) at the land-marked depths annotated in FIG. 11B :
- R 1 beyond 425 ⁇ m above RPE
- B R 2 —between 300 and 425 ⁇ m above RPE
- C R 3 —between 50 and 300 ⁇ m above RPE
- D C 1 —between 0 to 70 ⁇ m below RPE
- E C 2 —between 70 to 200 ⁇ m below RPE
- C 3 beyond 200 ⁇ m below RPE.
- Coupled may mean that two or more elements are in direct physical or electrical contact. However, “coupled” may also mean that two or more elements are not in direct contact with each other, but yet still cooperate or interact with each other.
- a phrase in the form “NB” or in the form “A and/or B” means (A), (B), or (A and B).
- a phrase in the form “at least one of A, B, and C” means (A), (B), (C), (A and B), (A and C), (B and C), or (A, B and C).
- a phrase in the form “(A)B” means (B) or (AB) that is, A is an optional element.
- the description may use the terms “embodiment” or “embodiments,” which may each refer to one or more of the same or different embodiments.
- the terms “comprising,” “including,” “having,” and the like, as used with respect to embodiments are synonymous, and are generally intended as “open” terms (e.g., the term “including” should be interpreted as “including but not limited to,” the term “having” should be interpreted as “having at least,” the term “includes” should be interpreted as “includes but is not limited to,” etc.).
- a computing device may be endowed with one or more components of the disclosed apparatuses and/or systems and may be employed to perform one or more methods as disclosed herein.
- Embodiments herein provide an UHS-OMAG system that delivers high sensitivity with a relatively low data acquisition time.
- OMAG is an imaging modality that is a variation on optical coherence tomography (OCT).
- OCT optical coherence tomography
- the imaging is based on the optical signals scattered by the moving particles.
- the light backscattered from a moving particle may carry a beating frequency that may be used to distinguish scattering signals by the moving elements from those by the static elements. Accordingly, OMAG can be used to image the flow of particles, such as blood flow.
- UHS-OMAG system may include a light source, a sample arm, a reference arm and a detection arm. Illustrated in FIG. 1 is an exemplary embodiment of an UHS-OMAG apparatus 100 suitable for ultrahigh sensitive 2-D and 3-D flow imaging.
- the illustrated UHS-OMAG apparatus 100 may include some features known in the art, features which may not be explained in great length herein except where helpful in the understanding of embodiments of the present invention.
- UHS-OMAG apparatus 100 may include a light source 10 .
- Light source 10 may comprise any light source suitable for the purpose including, but not limited to, a broadband light source or a tunable laser source.
- a suitable broadband light source 10 may include a superluminescent diode.
- light source 10 comprises a superluminescent diode with a central wavelength of 1310 nanometers (nm) and a full-width-at-half-maximum bandwidth of 65 nm.
- light source 10 may be a light source having one or more longer/shorter wavelengths, which may allow for deeper imaging.
- light source 10 may comprise a tunable laser source such as, for example, a swept laser source.
- UHS-OMAG apparatus 100 may include optics 11 to couple the light from the light source 10 into the system.
- the apparatus 100 may include a beam splitter 12 for splitting the light from the optics 11 into two beams: a first beam provided to a reference arm 14 and a second beam provided to a sample arm 16 .
- optics 11 may include, but are not limited to, various lenses or fiber optics components suitable for the purpose.
- Beam splitter 12 may comprise a 2 ⁇ 2 single-mode fiber coupler or any fiber coupler suitable for the purpose.
- Reference arm 14 may be configured to provide a reference light to a detection arm 30 (discussed more fully below), from the light provided by light source 10 , for producing a spectral interferogram in combination with backscattered light from sample 18 .
- Reference arm 14 may include optics 20 and a mirror 22 for reflecting light from light source 10 for providing the reference light.
- Optics 20 may include, but are not limited to, various lenses suitable for the purpose.
- Mirror 22 may be stationary or may be modulated. Modulation may be equivalent to frequency modulation of the detected signal at detection arm 30 . It has been observed that spectral interference signals (interferograms) may be modulated by a constant Doppler frequency by a modulated mirror 22 in the reference arm 14 . The spectral interference signal may then be recovered by de-modulating the modulated signal at the modulation frequency. De-modulation may be achieved using any suitable method including, for example, a digital or optical de-modulation method. Modulation and de-modulation of spectral interference signals may advantageously improve the signal-to-noise ratio, resulting in an improved image quality for structural, flow, and angiographic imaging.
- Sample arm 16 may be configured to provide light from light source 10 to a sample 18 by way of optics 24 , a scanner 26 , and optics 28 .
- Optics 24 may be used to couple the light from beam splitter 12 to scanner 26 .
- Optics 24 may include various optical lens, for example an optical collimator.
- Scanner 26 may include a pair of x-y galvanometer scanners for scanning sample 28 in an x-y direction.
- Optics 28 may comprise the appropriate optics for delivering the light from the scanner 26 onto sample 18 .
- scanner 26 may also receive backscattered light from sample 18 .
- the lateral imaging resolution may be approximately 16 micrometers ( ⁇ m) determined by an objective lens that focuses light onto sample 18 , with a light power on sample 18 being approximately 1 milliwatt (mW).
- detection arm 30 comprises a spectrometer 34 including one or more of various optics 36 including, but not limited to, one or more collimators, one or more diffracting/transmission gratings, and one or more lenses (not illustrated).
- optics 36 may include a 30-millimeter (mm) focal length collimator, a 1200 lines/mm diffracting grating, and an achromatic focusing lens with a 150 mm focal length.
- spectrometer 34 may have a designed spectral resolution of, for example, 0.055 nm, resulting in an optical range of approximately 6.4 mm in aft, where the positive frequency space is 3.2 mm and the negative frequency space is 3.2 mm.
- spectrometer 34 may have a designed spectral resolution of, for example, 0.055 nm, resulting in an optical range of approximately 6.4 mm in aft, where the positive frequency space is 3.2 mm and the negative frequency space is 3.2 mm.
- Such parameters are exemplary and may be modified in a variety of ways in accordance with the embodiments of the present invention.
- spectrometer 34 may include a detector such as a linear detector 38 configured to detect a spectral interference signal.
- Linear detector 38 may include one or more of a line-scan camera and an area scan camera.
- An exemplary suitable linear detector 38 may be a charge-coupled device (CCD).
- UHS-OMAG apparatus 100 may include a diffusion amplifier that may comprise one or more single element detectors rather than spectrometer 34 .
- a diffusion amplifier may comprise one or more single element detectors rather than spectrometer 34 .
- one or more dual-balanced photo-diode detectors may be used.
- reference arm 14 , sample arm 16 , and detection arm 30 may include polarization controllers (not illustrated). Polarization controllers may be configured to fine-tune the polarization states of light in UHS-OMAG apparatus 100 . Although an UHS-OMAG apparatus within the scope of the present invention may include more or less polarization controllers, inclusion of polarization controllers in reference arm 14 , sample arm 16 , and detection arm 30 , respectively, may advantageously maximize the spectral interference fringe contrast at linear detector 38 (or another suitable detector).
- UHS-OMAG apparatus 100 may include one or more user interfaces 40 for one or more purposes including controlling linear detector 38 and scanner 26 , computing data using algorithms, displaying images, input of data, and/or output of data, etc.
- UHS-OMAG apparatus 100 may be configured to build a 3-D data volume set by scanning sample 18 with a sample light in x, y, and ⁇ (z) directions to obtain a 3-D spectral interferogram data set.
- the scanner 26 may include an x-scanner and a y-scanner.
- the x-scanner may perform at least one fast scan along a fast scan axis
- the y-scanner may perform at least one slow scan along a slow scan axis.
- the fast scan axis may be orthogonal to the slow scan axis (i.e., the fast scan axis and slow scan axis may define an x-y plane).
- the fast scan may also be referred to herein as a B-scan, and the fast scan axis may also be referred to as the x-axis, the lateral axis, and/or the B-scan axis.
- the slow scan may also be referred to herein as a C-scan
- the slow scan axis may also be referred to as the y-axis, the elevational axis, and/or the C-scan axis.
- Each fast scan may be performed over a fast scan time interval, and each slow scan may be performed over a slow scan time interval, where the slow scan time interval is at least twice as long as the fast scan time interval.
- the scanner may perform the one or more fast scans contemporaneously with the one or more slow scans. In such embodiments, a plurality of fast scans may be performed during one slow scan.
- each B-scan there may be a number, N, of A-scans.
- An A-scan may be performed in the z-axis, orthogonal to both the x-axis and the y-axis.
- Each A-scan may include a number, K, of pixels, i.e., data points, that provide imaging depth information in the z-axis.
- a C-scan (slow scan) may include a number, M, of B-scans.
- the numbers N, K, and M may be at least two. Accordingly, during a composite scan, a three-dimensional (3D) data set may be produced.
- FIG. 2 shows an example of a 3D data set 200 in accordance with various embodiments.
- Data set 200 includes one C-scan 202 .
- C-scan 202 includes a number, M, of B-scans 204 .
- Each B-scan 204 includes a number, N, of A-scans 206 .
- Each A-scan 206 includes a number, K, of pixels 208 , i.e., data points 208 .
- an imaging algorithm may be applied to the 3D data set to produce at least one image.
- the imaging algorithm may be applied on the slow scan axis (i.e., y-axis).
- the imaging algorithm may separate a moving component from a structural component of the sample.
- the image may be a full range structural image and/or a separated structural/flow image.
- the image may be of blood flow, such as blood flow in a capillary and/or in a retina of an eye.
- the embodiments herein provide a higher sensitivity and a lower data acquisition time.
- the number, N, of A-scans in a B-scan may be decreased, while the number, M, of B-scans in a C-scan may be increased compared with the prior OMAG method.
- the imaging algorithm is applied on the slow axis (i.e., C-scan axis, elevational axis), rather than on the fast axis (i.e., B-scan axis, lateral axis) as in the prior OMAG method.
- This provides a higher sensitivity to the moving component of the image (i.e., can image slower speeds), while providing a lower data acquisition time.
- a high data acquisition time may be unsuitable for in vivo use, since involuntary movement of the subject is unavoidable.
- additional imaging algorithms are provided herein that differ from the prior OMAG method, although imaging algorithms used in the prior OMAG method may also be used in some embodiments.
- a red blood cell moves along the probe beam direction at a speed, v, of less than or equal to 200 ⁇ m/s (i.e., v ⁇ 200 ⁇ m/s)
- This time spacing translates into a scanning speed of about 643 A-scans/sec. Accordingly, the total imaging time to acquire a 3D capillary flow image of a tissue volume would be prohibitively long.
- the problem is exacerbated further under conditions where the probe beam is substantially perpendicular to the blow flow, such as when imaging the posterior segment of human eye, which makes the effective blood flow that is probed by the OMAG system very slow because the Doppler angle approaches 90 degrees.
- the algorithm is applied to the slow axis.
- the detectable flow signal is determined by the time spacing, ⁇ t B , between the adjacent B scans.
- the B frame rate may be 300 Hz (see example one below), so ⁇ t B may be about 3.3 milliseconds.
- the detectable flow would be about 140 ⁇ m/s, while the system scanning speed (i.e., A-scan rate) is not limited. Accordingly, the total imaging time to acquire a 3D capillary flow image of a tissue volume would be dramatically reduced.
- the fastest flow that can be detected by the UHS-OMAG system may determined by the system imaging speed.
- the imaging speed may have an A-scan rate of 47,000 per second (see example one below).
- the maximum detectable velocity may be about 30 mm/s.
- the slowest flow that can be detected by the UHS-OMAG system may be determined by the UHS-OMAG system phase noise floor.
- the system signal to noise ratio may be about 85 dB.
- the measured slowest detectable flow velocity may be about 4.0 ⁇ m/s.
- the x-scanner may be driven by a fast scan signal and the y-scanner may be driven by a slow scan signal. That is, the scanner may direct the probe beam along the fast scan axis and slow scan axis according to a property, such as a voltage, of the fast scan signal and the slow scan signal, respectively.
- the fast scan signal and slow scan signal may be any suitable waveforms, such as triangular, square, stepped, and/or sinusoidal.
- FIG. 3 shows an embodiment where the fast scan signal 302 and slow scan signal 304 are each a triangular waveform (i.e., continuous scan).
- the fast scan signal 302 has a higher frequency than the slow scan signal 304 .
- the slow scan is performed to obtain a C-scan
- a plurality of fast scans are performed to produce a plurality of B-scans.
- only one slow scan may be performed in the composite scan, as depicted in FIG. 3 .
- a plurality of slow scans may be performed and combined into the composite scans.
- FIG. 4 shows an alternative embodiment, in which fast scan signal 402 is a triangular waveform, and slow scan signal 404 is a stepped waveform.
- the slow scan signal may have a number, Q, of steps 406 to finish one C-scan.
- a number, P, of B-scans may be acquired.
- P may be at east two.
- any suitable imaging algorithm may be applied to the 3D data set to produce the image of the sample, such as a blood flow image.
- the imaging algorithm may be applied in the slow scan axis (i.e., the y-direction).
- the imaging algorithm may be designed to extract information on a moving component of the sample. Applying the imaging algorithm in the slow scan axis rather than the fast scan axis provides higher sensitivity, thereby allowing imaging of relatively slow movement, such as capillary blood flow within dermis which has typical flow in the range of about 100 to 900 ⁇ m/s in a resting condition, and even slower in diseased states.
- the imaging algorithm may include a differentiation operation in the y-direction, followed by an absolute operation.
- for j 1, 2, . . . , M).
- for j 1, 2, . . . , M).
- the imaging algorithm may include a high pass filter that may be applied to the 3D data set in the y-direction (slow scan axis) in order to separate the optical signals scattered by the moving particles from the optical signal scattered by the microstructures, i.e., the static particles.
- the high pass filter may use any suitable type of high pass filtering function.
- the prior OMAG method may be applied to the 3D data set in the y-direction (slow scan axis).
- averaging adjacent B-scans resulting from the imaging algorithm may improve the quality of final flow images.
- the number of B-scans used to average can be any number larger than, or equal to, 2.
- the imaging algorithm may be applied independently at each step. Then, the image at each step (i.e., I′(x i ,y j ,z k )), may be averaged to obtain a final B-scan image of the blood flow at that step. In some embodiments, the averaging may be done according to the following equation (Eq. (1)):
- the final 3D image may include a number, Q, of B-scan images in the y-direction.
- the algorithms described above may be applied to the 3D data set in the spectral/frequency domain (i.e., the data are in wavelength (or wavenumber) format). Further, in various embodiments, a logarithm operation may be applied to the 3D data set.
- the algorithms described above may be applied to the 3D data set in the time/distance domain (i.e., the data are in time (or distance) format).
- the time/distance domain signal and the spectral/frequency domain signal are a Fourier transform pair.
- an imaging tool for such an application must be able to resolve the capillary blood flows within dermis, which are normally very slow (in the range of about 100 to 900 ⁇ m/s at the resting condition, and even slower at diseased states).
- such tools must be able to provide depth information with an imaging resolution at a scale of capillary blood vessel (about 10 ⁇ m).
- the imaging tool must have a relatively low data acquisition time to allow in vivo use, since involuntary movement of the subject is unavoidable.
- the system setup for UHS-OMAG used in the experiment is similar to that described in R. K Wang, and L. An, “Doppler optical micro-angiography for volumetric imaging of vascular perfusion in vivo,” Opt. Express 17, 8926-8940 (2009) (hereinafter “Article 1”), which is hereby incorporated by reference in its entirety.
- the system used a superluminescent diode as the light source, which has a central wavelength of 1310 nm and a bandwidth of 65 nm that provided an axial resolution in air of about 12 ⁇ m.
- a 50 mm focal length objective lens was used to achieve a lateral resolution of about 16 ⁇ m.
- the output light from the interferometer was routed to a home-built spectrometer, which had a designed spectral resolution of about 0.141 nm that provided a detectable depth range of about 3 mm on each side of the zero delay line.
- the line scan rate of the camera was 47,000 per second (47 kHz). With this imaging speed, the signal to noise ratio was measured at about 85 dB with a light power on the sample of about 3 mW.
- the system applied a scanning protocol, similar to FIG. 3 (as discussed above), that was designed to achieve ultrahigh sensitive imaging to the blood flow.
- a scanning protocol similar to FIG. 3 (as discussed above), that was designed to achieve ultrahigh sensitive imaging to the blood flow.
- 128 A-lines were acquired with a spacing of about 15 ⁇ m between adjacent lines, thus covering a length of about 2 mm on the tissue.
- the imaging rate was 300 frames per second (fps). Note that with a 47 kHz line scan rate, the theoretical imaging rate should be 367 fps.
- the reduced imaging rate at 300 fps was due to data transfer limitations during the handshake between the camera and the computer.
- the self cross-correlation between the light backscattered from different positions within the sample is not considered in Eq. (2), because the light backscattered from the sample is weak compared to the light reflected from the reference mirror. Additionally, the DC signals are not considered, because they do not contribute to useful OMAG signals.
- Eq. (2) can still be used to represent the spectral interferogram signal captured by the system, except that the time variable, t, now corresponds to the B-scan numbers in one C-scan.
- the imaging rate is 300 fps, so ⁇ t B is about 3.3 ms. Note that the imaging speed is at 47,000 A scans per second, in contrast to 213 A scan per second required for prior OMAG method.
- the imaging algorithm first takes a differential operation on the captured B-scan spectral interferograrns along the C-scan direction, as in the following equation (Eq. (3)):
- i represents the index of the B-scans in the C-scan direction.
- the differential operation suppresses the optical scattering signals from the static elements within scanned tissue volume.
- high pass filtering may be used.
- FFT fast Fourier transform
- S intensity signal to noise ratio
- the minimum detectable flow velocity would be about 4.0 ⁇ m/s.
- the system described here would not provide a continuous trajectory for this blood cell in the 3D OMAG flow image, i.e., the trajectory would be seen as a broken line.
- the phase compensation method described in Article 1 is applied to the raw interference signal before applying the OMAG algorithm.
- the system was tested on the skin located on the backside of a hand of a male volunteer.
- the traditional OMAG and phase resolved optical Doppler tomography (PRODT) cross-sectional flow images were also obtained.
- the system captured 2000 A-scans over 2 mm at an imaging speed of 31,000 A-scans per second in order to fulfill the oversampling requirement for these previous methods.
- FIGS. 5A-F The results are shown in FIGS. 5A-F .
- the images in the top row ( FIGS. 5A-C ) are from the UHS-OMAG system while those in the bottom row ( FIGS. 5D-F ) are from the conventional prior OMAG system.
- FIGS. 5A and 5D are the Fourier domain optical coherence tomography (FDOCT) structural images obtained from the captured interferograms for the UHS-OMAG system and prior OMAG system, respectively. Although they are similar, they are not exactly the same due to the small subject movement when switching the system among the different approaches. However, it is sufficient to provide a fair comparison of their ability to extract slow flow information.
- FIG. 5B shows the image from FIG. 5A after it is processed by the UHS-OMAG system.
- FIG. 5B shows the image from FIG. 5A after it is processed by the UHS-OMAG system.
- FIG. 5C shows the PRODT image of the UHS-OMAG system, which is based on the phase difference between adjacent B-scans.
- FIG. 5E shows the image from FIG. 50 after it is processed by the prior OMAG system.
- FIG. 5F shows the conventional PRODT image of the prior OMAG system, which is based on the phase difference between adjacent A-scans in one B-scan. In FIGS. 50 and 5F , the phase differences are calculated only when the structural signal is 15 dB above the noise floor.
- FIG. 5B shows blood flows within the papillary dermis (indicated by white arrows), where only capillary blood vessels are present, as well as the blood flows within reticular dermis (indicated by red arrows), where both the capillary and larger blood vessels are present.
- the conventional OMAG requires oversampling in the fast scanning direction (i.e., B-scan), it is not sensitive to the slow blood flow within the capillaries, which are normally below 100 ⁇ m/s. Accordingly, the conventional PRODT approach totally failed in imaging any of blood vessels, as seen in FIG. 5F . It should be noted that there is a global noise ‘flow’ background in the UHS-OMAG flow image, e.g., in FIG. 5B , which may be caused by some ‘non-moving’ scatters, such as global motion, etc. In this case, de-noising filters may be used to further enhance the UHS-OMAG flow imaging quality.
- FIGS. 6A-B The imaging results are shown in FIGS. 6A-B .
- FIG. 6A is the OMAG/OCT microstructural image of the flow phantom
- FIG. 6B is the corresponding UHS-OMAG flow image. From this result, it is clear that UHS-OMAG is able to image the particle movements due to Brownian motion with almost no signals detected in the background region.
- FIG. 7A shows a plot of the calculated velocities across the center of the capillary tube at the position marked as the horizontal line in FIG. 7A , where the dashed box indicates the position of the capillary lumen.
- the velocity values of particle movements ranged from approximately ⁇ 50 to 100 ⁇ m/s at this cross-line position.
- the conventional PRODT image was also obtained of the same phantom.
- the system imaging rate was set at 31,000 A-scans per second.
- the A-line density across the B-scan of about 2.5 mm was set at 4000, equating to a spacing between adjacent A-scans of about 0.625 ⁇ m.
- FIGS. 7C and 7D respectively, showing that PRODT totally failed to image the Brownian motion of the particles under the current experimental setup. Note that the standard deviation of velocity values shown in FIG. 7D was about 180 ⁇ m/s, thus it is not surprising that PRODT is not able to achieve satisfactory imaging performance.
- FIG. 8A shows a schematic drawing of the blood vessel system of the human skin, in which an interconnected network of vessels is characterized by regular structures on all levels.
- the human skin is composed of the cutis and the subcutis (hypodermis [HD]).
- the cutis is further divided into the epidermis (EP) and the dermis (DR).
- the dermis and subcutis are pervaded with a complex system of blood vessels, while the epidermis is free of vessels.
- a superficial network comprises the interface between papillary (PD) and reticular dermis (RD), while a lower network is located on the border between dermis and subcutis.
- Vertical vessels connect both networks and thus make it complete.
- arteries are shown in red and veins in blue.
- 3D blood flow images were acquired over the palm of a healthy volunteer, as shown in FIG. 8B , where the black box indicates the scanning area (about 2 ⁇ 2 mm 2 ).
- FIG. 8C The 3D OMAG imaging result of blood vessel networks is shown in FIG. 8C , shown together with the 3D micro-structural image.
- FIG. 8D to shows a cross-sectional view, where the blood flows within blood vessel systems within the skin are clearly delineated. Because the UHS-OMAG sensitivity is as low as about 4 ⁇ m/s, even the dynamics of sweat glands are imaged.
- FIGS. 9A-D The projection views at different land-mark depths are shown in FIGS. 9A-D .
- FIG. 9A gives the projection view at the depths from 400 to 450 ⁇ m, which corresponds to the papillary dermis where the capillary vessels are dense (e.g., indicated by the arrows).
- FIG. 9B shows the depths from 450 to 650 ⁇ m, where the vessels are almost vertical that connect vessel networks between papillary dermis and reticular dermis, seen as the bright spots in FIG. 9B .
- the blood vessel network in reticular dermis (650 to 780 ⁇ m) and hypodermis (780 to 1100 ⁇ m) are shown in FIGS. 9C and 9D , respectively.
- the vessel diameter is smaller in the reticular dermis than in the hypodermis.
- FIG. 10 shows the setup of the UHS-OMAG system 1000 used to obtain the images.
- the UHS-OMAG system 1000 is similar to that described in L. An, and R. K. Wang, Optics Express, 16, 11438-11452 (2008) (hereinafter “Article 2”).
- the UHS-OMAG system 1000 used includes a light source 1002 to produce a probe beam with a wavelength centered at 842 nm and a bandwidth of 46 nm that provided an axial resolution in air of about 8 ⁇ m.
- a sample arm 1004 the light was delivered onto a human eye 1006 by a collimator 1008 , an objective lens 1010 , and an ocular lens 1012 .
- the sample arm further includes an x-scanner 1014 and a y-scanner 1016 for scanning the human eye 1006 along the x-axis and y-axis, respectively.
- a 20 mm water cell 1020 was used to compensate for the dispersion caused by the eye.
- a reference beam from the light source was reflected off a reference mirror 1022 to provide a reference light.
- An interferogram between the reference light and the light backscattered from the sample was sent to a home-built ultrafast spectrometer 1024 via an optical circulator 1026 .
- the spectrometer included a collimator 1027 , a transmission grating 1028 (with 1200 lines/mm), a camera lens 1030 with a focal length of 100 mm, and a 1024 element line scan CMOS detector 1032 capable of a 140 kHz line scan rate.
- the spectral resolution of the designed spectrometer 1024 was about 0.055 nm, which provided an imaging depth of about 3 mm in air.
- the system sensitivity was about 90 dB measured with about 900 ⁇ W power of light incident at the object and an exposure time of 6.9 ⁇ s with the time interval, ⁇ t A , between A-scans equal to about 7.4 ⁇ s.
- the system 1000 further includes polarization controllers 1034 , 1036 , and 1038 to control polarization of the light beam.
- the UHS-OMAG system 1000 employed a scanning protocol, similar to FIG. 4 (as discussed above), designed to achieve ultrahigh sensitive imaging of the flow.
- the x-scanner 1012 was driven by a 400 Hz sawtooth waveform, meaning that the imaging rate is 400 frames/sec (fps).
- the duty cycle for acquiring each B-scan i.e., x-direction scan
- the y-scanner was driven by a step function, and the entire C-scan includes 150 steps, with a spacing between adjacent steps of about 20 ⁇ m. In each step, 8 repeated B-scans were acquired. Accordingly, it required 3 seconds to acquire one 3D data set, covering an area of about 3 ⁇ 3 mm on the retina.
- the calculated OMAG signals at each step were collapsed into one B-scan through ensemble-averaging (e.g., Eq. (1)), resulting in 150 B-scans to form the final C-scan image, i.e. 3D OMAG blood flow distribution.
- FIGS. 11A-D show the in vivo imaging results produced by one volume dataset captured at the macular region towards the optic nerve head.
- FIG. 11A shows a typical cross-sectional image (B-scan) within the OMAG structural volume, which is identical to the conventional OCT image where the typical morphological features within retina and choroid are visualized.
- FIG. 11B gives the corresponding blood flow image obtained from the imaging algorithm, where the capillary flows within the cross-section of retina are abundant (indicated by arrows), as well as the blood flow signals in the choroid.
- the blood flows in the retina can be separated from those in the choroid.
- the segmentation algorithm described in Article 2 was first used to identify the retinal pigment epithelium (RPE) layer. Then, the flow signals from the retinal vessels were identified as the OMAG flow signals 50 ⁇ m above the RPE layer (to exclude the signals from photoreceptor inner and outer segments), while those below the RPE layer were identified as the choroidal vessels.
- the segmentation resulted in two volumetric flow images, one for retina and another for choroid, which are annotated in FIG. 11B by the labels of R and C, respectively.
- FIGS. 11C and 11 D the blood flow distribution maps shown in FIGS. 11C (retina) and 11 D (choroid), respectively.
- FIG. 11C there is seen a ring of blood vessels in the macular area around an avascular zone about 800 ⁇ m in diameter, denoting the fovea. This observation is in excellent agreement with the standard retinal pathology.
- the retina consists of three layers of capillary networks: the radial peripapillary capillaries (RPCs, R 1 ) and an inner (R 2 ) and an outer layer of capillaries (R 3 ).
- the RPCs are the most superficial layer of capillaries lying in the inner part of the nerve fiber layer.
- the inner capillaries lie in the ganglion cell layers under and parallel to the RPCs.
- the outer capillary network runs from the inner plexiform to the outer plexiform through the inner nuclear layers.
- the choroidal vessels in the macular region are not specialized like those in the retina.
- the arteries pierce the sclera around the optic nerve and fan out to form three vascular layers in the choroid: inner capillary bed (near RPE layer, C 1 ), medial arterioles and venules (C 2 ) and outer arteries and veins (C 3 ).
- inner capillary bed near RPE layer, C 1
- medial arterioles and venules C 2
- outer arteries and veins C 3
- FIGS. 12A-F FIGS. 12A-C for retina and FIGS. 12D-F for choroid. The results are correlated well with the descriptions found in the literature.
- the OMAG system running at 400 fps requires about 3 seconds to acquire one 3D blood flow image representing about 3 ⁇ 3 mm 2 area on the retina. At this speed, the effects of subject movement on the final results are clearly noticeable (the horizontal lines in FIGS. 11A-D and 12 A-F are the motion artifacts).
- the phase compensation algorithm developed in Article 2 may be used to minimize the motion artifacts before the imaging algorithm is applied. In doing so, however, it would inevitably increase the computational load required to obtain the meaningful 3D blood flow images because of the complexity of the compensation algorithms, and/or 2) the imaging speed may be further increased to minimize the motion artifacts.
- the current system speed is limited by the CMOS camera used in the spectrometer with a maximum line rate of 140 kHz.
- a report has shown that a 4096 element CMOS camera is capable of a line rate more than 240 kHz. Therefore, it would be expected that if this CMOS camera is used in the system, the motion artifacts may be mitigated.
- the UHS-OMAG system may be capable of imaging detailed ocular perfusion distributions within retina and choroid. Applying the imaging algorithm on the slow scanning axis, the system is sensitive to the ocular capillary flows. In addition, due to the depth-resolved nature, the system may be able to provide detailed micro-circulation within different land-marked depths, the results of which are in an excellent agreement with those described in the literature. The demonstrated superior imaging results show promise for future clinical applications for UHS-OMAG in ophthalmology.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Hematology (AREA)
- Physiology (AREA)
- Cardiology (AREA)
- Ophthalmology & Optometry (AREA)
- Dermatology (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
Embodiments herein provide an ultrahigh sensitive optical microangiography (OMAG) system that provides high sensitivity to slow flow information, such as that found in blood flow in capillaries, while also providing a relatively low data acquisition time. The system performs a plurality of fast scans (i.e., B-scans) on a fast scan axis, where each fast scan includes a plurality of A-scans. At the same time, the system performs a slow scan (i.e., C-scan), on a slow scan axis, where the slow scan includes the plurality of fast scans. A detector receives the spectral interference signal from the sample to produce a three dimensional (3D) data set. An imaging algorithm is then applied to the 3D data set in the slow scan axis to produce at least one image of the sample. In some embodiments, the imaging algorithm may separate flow information from structural information of the sample.
Description
- The present application is a continuation of U.S. patent application Ser. No. 13/577,857 filed Aug. 8, 2012 entitled “Method and Apparatus for Ultrahigh Sensitive Optical Microangiography,” which is a U.S. National Phase Application of PCT/US2011/024069 entitled “Method and Apparatus for Ultrahigh Sensitive Optical Microangiography,” which claims priority to U.S. Provisional Patent Application No. 61/302,409, filed Feb. 8, 2010, entitled “Method and Apparatus for Ultrahigh Sensitive Optical Microangiography,” the entire disclosures of which are hereby incorporated by reference in their entirety.
- The present application is related to U.S. Provisional Patent Application No. 61/175,229, filed May 4, 2009, entitled “Method and Apparatus for Quantitative Imaging of Blood Perfusion in Living Tissue” and Publication No. WO2008/039660, filed Sep. 18, 2007, entitled “In Vivo Structural and Flow Imaging,” the entire disclosures of which are hereby incorporated by reference in their entirety.
- This invention was made with Government support under Grant/Contract No. R01HL093140, R01EB009682 and R01DC010201 awarded by the US National Institute of Health. The Government has certain rights in the invention.
- Embodiments herein relate to the field of imaging, and, more specifically, to a method and apparatus for ultrahigh sensitive optical microangiography.
- The assessment of blood flow in living tissue provides important information for diagnostics, treatment and/or management of pathological conditions. For example, the assessment of cutaneous (skin) microcirculations may provide important information for pathological conditions in dermatology, such as skin cancer, port wine stain treatment, diabetes, and plastic surgery. Similarly, assessment of the ocular perfusion within the retina and choroid of the human eye is important in the diagnosis, treatment, and management of a number of pathological conditions in ophthalmology, such as age-related macular degeneration, diabetic retinopathy, and glaucoma. Accordingly, clinical and technical tools that can noninvasively image three dimensional (3D) micro-blood vessel networks in vivo are in demand.
- Several techniques have been developed to meet this need. However, current techniques suffer from various shortcomings which make them unsuitable for in vivo imaging in humans, such as low sensitivity to blood flow, insufficient resolution to provide useful depth information, and/or a long data acquisition time.
- Embodiments will be readily understood by the following detailed description in conjunction with the accompanying drawings. Embodiments are illustrated by way of example and not by way of limitation in the figures of the accompanying drawings.
-
FIG. 1 illustrates a functional block diagram of an imaging apparatus in accordance with various embodiments of the present invention -
FIG. 2 illustrates an example of a 3D data set, in accordance with various embodiments; -
FIG. 3 illustrates an example of drive signals for an x-scanner and y-scanner, in accordance with various embodiments; -
FIG. 4 illustrates another example of drive signals for an x-scanner and y-scanner, in accordance with various embodiments; -
FIGS. 5A-C show images of capillaries within human dermis, taken with a ultrahigh sensitive optical microangiography (UHS-OMAG) system according to various embodiments: (A) a Fourier domain optical coherence tomography (FDOCT) structural image; (B) a UHS-OMAG flow image; and (C) a phase resolved optical Doppler tomography (PRODT) cross-sectional flow image from the UHS-OMAG image ofFIG. 5B ; -
FIGS. 5D-F show images of capillaries within human dermis, taken with a prior optical microangiography (OMAG) system, for comparison with the images inFIGS. 5A-C : (D) a Fourier domain optical coherence tomography (FDOCT) structural image; (E) an OMAG flow image; and (F) a phase resolved optical Doppler tomography (PRODT) cross-sectional flow image from the OMAG image ofFIG. 5E ; -
FIG. 6A is a B-scan structural image of a flow phantom, and Hg. 6B is a corresponding UHS-OMAG flow image, in accordance with various embodiments. -
FIG. 7A is a velocity image obtained from a flow phantom assessed by UHS-OMAG, in accordance with various embodiments; -
FIG. 7B is a plot of the velocity data across the capillary tube at the position shown as the horizontal line inFIG. 7A , in accordance with various embodiments; -
FIG. 7C andFIG. 7D are the corresponding results toFIGS. 7A and 7B , respectively, obtained by PRODT imaging of the same phantom, in accordance with various embodiments; -
FIG. 8A is a schematic diagram of the blood vessel system of the human skin; -
FIG. 8B is a photograph of a palm showing the scanning area; -
FIG. 8C is a three dimensional (3D) rendered OMAG image of blood vessels together with the 3D structures, in accordance with various embodiments; -
FIG. 8D is a cross sectional view of the imaged blood vessels, in accordance with various embodiments; -
FIGS. 9A-D show UHS-OMAG detailed projection views of microcirculation network at different depths of skin obtained from: (A) 400-450 μm (closely representing papillary dermis); (B) 450-650 μm; (C) 650-780 μm (closely representing reticular dermis); and (D) 780-1100 μm (part of hypodermis), respectively, in accordance with various embodiments; -
FIG. 10 shows a schematic diagram of a UHS-OMAG system, in accordance with various embodiments; -
FIGS. 11A-B show in vivo UHS-OMAG imaging, in accordance with various embodiments, of the posterior segment of eye near the macular region towards the optic nerve head: (A) an OMAG B-scan of microstructures showing morphological features, and (B) the corresponding OMAG blood flow image; -
FIGS. 11C-D show projection maps of blood flow distribution within: (C) retina and (D) choroid, obtained from one 3D scan of an area of ˜3×3 mm2 near the macular region, in accordance with various embodiments; and -
FIGS. 12A-F show depth-resolved images, in accordance with various embodiments, of patent blood vessels within the retina (FIGS. 12A-C ) and choroid (FIGS. 12D-F ) at the land-marked depths annotated inFIG. 11B : (A) R1—beyond 425 μm above RPE; (B) R2—between 300 and 425 μm above RPE; (C) R3—between 50 and 300 μm above RPE; (D) C1—between 0 to 70 μm below RPE; (E) C2—between 70 to 200 μm below RPE; and (F) C3—beyond 200 μm below RPE. - In the following detailed description, reference is made to the accompanying drawings which form a part hereof, and in which are shown by way of illustration embodiments that may be practiced. It is to be understood that other embodiments may be utilized and structural or logical changes may be made without departing from the scope. Therefore, the following detailed description is not to be taken in a limiting sense, and the scope of embodiments is defined by the appended claims and their equivalents.
- Various operations may be described as multiple discrete operations in turn, in a manner that may be helpful in understanding embodiments; however, the order of description should not be construed to imply that these operations are order dependent.
- The description may use perspective-based descriptions such as up/down, back/front, and top/bottom. Such descriptions are merely used to facilitate the discussion and are not intended to restrict the application of disclosed embodiments.
- The terms “coupled” and “connected,” along with their derivatives, may be used. It should be understood that these terms are not intended as synonyms for each other. Rather, in particular embodiments, “connected” may be used to indicate that two or more elements are in direct physical or electrical contact with each other. “Coupled” may mean that two or more elements are in direct physical or electrical contact. However, “coupled” may also mean that two or more elements are not in direct contact with each other, but yet still cooperate or interact with each other.
- For the purposes of the description, a phrase in the form “NB” or in the form “A and/or B” means (A), (B), or (A and B). For the purposes of the description, a phrase in the form “at least one of A, B, and C” means (A), (B), (C), (A and B), (A and C), (B and C), or (A, B and C). For the purposes of the description, a phrase in the form “(A)B” means (B) or (AB) that is, A is an optional element.
- The description may use the terms “embodiment” or “embodiments,” which may each refer to one or more of the same or different embodiments. Furthermore, the terms “comprising,” “including,” “having,” and the like, as used with respect to embodiments, are synonymous, and are generally intended as “open” terms (e.g., the term “including” should be interpreted as “including but not limited to,” the term “having” should be interpreted as “having at least,” the term “includes” should be interpreted as “includes but is not limited to,” etc.).
- With respect to the use of any plural and/or singular terms herein, those having skill in the art can translate from the plural to the singular and/or from the singular to the plural as is appropriate to the context and/or application. The various singular/plural permutations may be expressly set forth herein for sake of clarity.
- In various embodiments, methods, apparatuses, and systems for ultrahigh sensitive optical microangiography (UHS-OMAG) are provided. In exemplary embodiments, a computing device may be endowed with one or more components of the disclosed apparatuses and/or systems and may be employed to perform one or more methods as disclosed herein.
- Embodiments herein provide an UHS-OMAG system that delivers high sensitivity with a relatively low data acquisition time. OMAG is an imaging modality that is a variation on optical coherence tomography (OCT). The imaging is based on the optical signals scattered by the moving particles. The light backscattered from a moving particle may carry a beating frequency that may be used to distinguish scattering signals by the moving elements from those by the static elements. Accordingly, OMAG can be used to image the flow of particles, such as blood flow.
- Various embodiments of the UHS-OMAG system may include a light source, a sample arm, a reference arm and a detection arm. Illustrated in
FIG. 1 is an exemplary embodiment of an UHS-OMAG apparatus 100 suitable for ultrahigh sensitive 2-D and 3-D flow imaging. The illustrated UHS-OMAG apparatus 100 may include some features known in the art, features which may not be explained in great length herein except where helpful in the understanding of embodiments of the present invention. - As illustrated, UHS-
OMAG apparatus 100 may include alight source 10.Light source 10 may comprise any light source suitable for the purpose including, but not limited to, a broadband light source or a tunable laser source. A suitable broadbandlight source 10 may include a superluminescent diode. In one embodiment,light source 10 comprises a superluminescent diode with a central wavelength of 1310 nanometers (nm) and a full-width-at-half-maximum bandwidth of 65 nm. In various embodiments,light source 10 may be a light source having one or more longer/shorter wavelengths, which may allow for deeper imaging. In various other embodiments,light source 10 may comprise a tunable laser source such as, for example, a swept laser source. - UHS-
OMAG apparatus 100 may includeoptics 11 to couple the light from thelight source 10 into the system. Theapparatus 100 may include abeam splitter 12 for splitting the light from theoptics 11 into two beams: a first beam provided to areference arm 14 and a second beam provided to asample arm 16. In various embodiments,optics 11 may include, but are not limited to, various lenses or fiber optics components suitable for the purpose.Beam splitter 12 may comprise a 2×2 single-mode fiber coupler or any fiber coupler suitable for the purpose. -
Reference arm 14 may be configured to provide a reference light to a detection arm 30 (discussed more fully below), from the light provided bylight source 10, for producing a spectral interferogram in combination with backscattered light fromsample 18.Reference arm 14 may includeoptics 20 and amirror 22 for reflecting light fromlight source 10 for providing the reference light.Optics 20 may include, but are not limited to, various lenses suitable for the purpose. -
Mirror 22 may be stationary or may be modulated. Modulation may be equivalent to frequency modulation of the detected signal atdetection arm 30. It has been observed that spectral interference signals (interferograms) may be modulated by a constant Doppler frequency by a modulatedmirror 22 in thereference arm 14. The spectral interference signal may then be recovered by de-modulating the modulated signal at the modulation frequency. De-modulation may be achieved using any suitable method including, for example, a digital or optical de-modulation method. Modulation and de-modulation of spectral interference signals may advantageously improve the signal-to-noise ratio, resulting in an improved image quality for structural, flow, and angiographic imaging. -
Sample arm 16 may be configured to provide light fromlight source 10 to asample 18 by way ofoptics 24, ascanner 26, andoptics 28.Optics 24 may be used to couple the light frombeam splitter 12 toscanner 26.Optics 24 may include various optical lens, for example an optical collimator.Scanner 26 may include a pair of x-y galvanometer scanners for scanningsample 28 in an x-y direction.Optics 28 may comprise the appropriate optics for delivering the light from thescanner 26 ontosample 18. In various embodiments,scanner 26 may also receive backscattered light fromsample 18. Although the characteristics of the light provided to sample 18 may depend on the particular application, in some embodiments, the lateral imaging resolution may be approximately 16 micrometers (μm) determined by an objective lens that focuses light ontosample 18, with a light power onsample 18 being approximately 1 milliwatt (mW). - The light returning from
reference arm 14 and the light returning from sample arm 16 (i.e., the spectral signal) may be recombined and coupled into thebeam splitter 12 for introduction todetection arm 30. As illustrated,detection arm 30 comprises aspectrometer 34 including one or more ofvarious optics 36 including, but not limited to, one or more collimators, one or more diffracting/transmission gratings, and one or more lenses (not illustrated). In exemplary embodiments,optics 36 may include a 30-millimeter (mm) focal length collimator, a 1200 lines/mm diffracting grating, and an achromatic focusing lens with a 150 mm focal length. In various embodiments,spectrometer 34 may have a designed spectral resolution of, for example, 0.055 nm, resulting in an optical range of approximately 6.4 mm in aft, where the positive frequency space is 3.2 mm and the negative frequency space is 3.2 mm. Such parameters are exemplary and may be modified in a variety of ways in accordance with the embodiments of the present invention. - In embodiments employing a broadband light source,
spectrometer 34 may include a detector such as alinear detector 38 configured to detect a spectral interference signal.Linear detector 38 may include one or more of a line-scan camera and an area scan camera. An exemplary suitablelinear detector 38 may be a charge-coupled device (CCD). - In embodiments wherein
light source 10 comprises a tunable laser rather than a broadband light source, however, UHS-OMAG apparatus 100 may include a diffusion amplifier that may comprise one or more single element detectors rather thanspectrometer 34. For example, one or more dual-balanced photo-diode detectors may be used. - In various embodiments,
reference arm 14,sample arm 16, anddetection arm 30 may include polarization controllers (not illustrated). Polarization controllers may be configured to fine-tune the polarization states of light in UHS-OMAG apparatus 100. Although an UHS-OMAG apparatus within the scope of the present invention may include more or less polarization controllers, inclusion of polarization controllers inreference arm 14,sample arm 16, anddetection arm 30, respectively, may advantageously maximize the spectral interference fringe contrast at linear detector 38 (or another suitable detector). - In various embodiments, UHS-
OMAG apparatus 100 may include one or more user interfaces 40 for one or more purposes including controllinglinear detector 38 andscanner 26, computing data using algorithms, displaying images, input of data, and/or output of data, etc. - As noted above, UHS-
OMAG apparatus 100 may be configured to build a 3-D data volume set by scanningsample 18 with a sample light in x, y, and λ (z) directions to obtain a 3-D spectral interferogram data set. - In various embodiments, the
scanner 26 may include an x-scanner and a y-scanner. During the composite scan, the x-scanner may perform at least one fast scan along a fast scan axis, and the y-scanner may perform at least one slow scan along a slow scan axis. The fast scan axis may be orthogonal to the slow scan axis (i.e., the fast scan axis and slow scan axis may define an x-y plane). The fast scan may also be referred to herein as a B-scan, and the fast scan axis may also be referred to as the x-axis, the lateral axis, and/or the B-scan axis. Similarly, the slow scan may also be referred to herein as a C-scan, and the slow scan axis may also be referred to as the y-axis, the elevational axis, and/or the C-scan axis. Each fast scan may be performed over a fast scan time interval, and each slow scan may be performed over a slow scan time interval, where the slow scan time interval is at least twice as long as the fast scan time interval. In some embodiments, the scanner may perform the one or more fast scans contemporaneously with the one or more slow scans. In such embodiments, a plurality of fast scans may be performed during one slow scan. - In each B-scan (fast scan), there may be a number, N, of A-scans. An A-scan may be performed in the z-axis, orthogonal to both the x-axis and the y-axis. Each A-scan may include a number, K, of pixels, i.e., data points, that provide imaging depth information in the z-axis. Similarly, a C-scan (slow scan) may include a number, M, of B-scans. In various embodiments, the numbers N, K, and M may be at least two. Accordingly, during a composite scan, a three-dimensional (3D) data set may be produced. The 3D data set may be represented as the complex function I(xi,yj,zk), where i=1, 2, . . . N, j=1, 2, . . . M, and k=1, 2, . . . , K. In various embodiments, the magnitude portion of the 3D data set may be represented by the scalar function A(xi,yj,zk), where I(xi,yj,zk)=A(xi,yj,zk)exp(iφ).
-
FIG. 2 shows an example of a3D data set 200 in accordance with various embodiments.Data set 200 includes one C-scan 202. C-scan 202 includes a number, M, of B-scans 204. Each B-scan 204 includes a number, N, ofA-scans 206. Each A-scan 206 includes a number, K, ofpixels 208, i.e., data points 208. - In various embodiments, an imaging algorithm may be applied to the 3D data set to produce at least one image. The imaging algorithm may be applied on the slow scan axis (i.e., y-axis). In some embodiments, the imaging algorithm may separate a moving component from a structural component of the sample. The image may be a full range structural image and/or a separated structural/flow image. In some embodiments, the image may be of blood flow, such as blood flow in a capillary and/or in a retina of an eye.
- In this disclosure, reference is made to a prior OMAG method developed by the present applicants and described in Publication No. WO2008/039660, filed Sep. 18, 2007, entitled “In Vivo Structural and Flow Imaging,” the entire disclosure of which is hereby incorporated by reference. As compared to the prior OMAG method, the embodiments herein provide a higher sensitivity and a lower data acquisition time. In some embodiments, the number, N, of A-scans in a B-scan may be decreased, while the number, M, of B-scans in a C-scan may be increased compared with the prior OMAG method. Additionally, the imaging algorithm is applied on the slow axis (i.e., C-scan axis, elevational axis), rather than on the fast axis (i.e., B-scan axis, lateral axis) as in the prior OMAG method. This provides a higher sensitivity to the moving component of the image (i.e., can image slower speeds), while providing a lower data acquisition time. A high data acquisition time may be unsuitable for in vivo use, since involuntary movement of the subject is unavoidable. Furthermore, additional imaging algorithms are provided herein that differ from the prior OMAG method, although imaging algorithms used in the prior OMAG method may also be used in some embodiments.
- In the prior OMAG method, high pass filtering is normally applied in the B scan frames (obtained in the fast scanning axis) to isolate the optical scattering signals between the static and moving scatters. The detectable flow velocity, v, is thus determined by the time spacing, ΔtA, between the adjacent A scans, i.e., v=λ/2nΔtA, where λ is the central wavelength of the light source, and n is the refractive index of the sample. If a red blood cell moves along the probe beam direction at a speed, v, of less than or equal to 200 μm/s (i.e., v≦200 μm/s), then it would require a time spacing, ΔtA, of at least about 1.5 milliseconds (i.e., ΔtA≧˜1.5 milliseconds) for the system to sample the moving blood cell (assuming λ=840 nm and n=1.35). This time spacing translates into a scanning speed of about 643 A-scans/sec. Accordingly, the total imaging time to acquire a 3D capillary flow image of a tissue volume would be prohibitively long. The problem is exacerbated further under conditions where the probe beam is substantially perpendicular to the blow flow, such as when imaging the posterior segment of human eye, which makes the effective blood flow that is probed by the OMAG system very slow because the Doppler angle approaches 90 degrees.
- In this disclosure, however, the algorithm is applied to the slow axis. Accordingly, the detectable flow signal is determined by the time spacing, ΔtB, between the adjacent B scans. In one exemplary embodiment, the B frame rate may be 300 Hz (see example one below), so ΔtB may be about 3.3 milliseconds. In this case, the detectable flow would be about 140 μm/s, while the system scanning speed (i.e., A-scan rate) is not limited. Accordingly, the total imaging time to acquire a 3D capillary flow image of a tissue volume would be dramatically reduced.
- In various embodiments, the fastest flow that can be detected by the UHS-OMAG system may determined by the system imaging speed. In one exemplary embodiment, the imaging speed may have an A-scan rate of 47,000 per second (see example one below). Under this circumstance, the maximum detectable velocity may be about 30 mm/s.
- In various embodiments, the slowest flow that can be detected by the UHS-OMAG system may be determined by the UHS-OMAG system phase noise floor. In one exemplary embodiment, with the imaging speed of 47,000 per second (see example one below), the system signal to noise ratio may be about 85 dB. The measured slowest detectable flow velocity may be about 4.0 μm/s.
- In various embodiments, the x-scanner may be driven by a fast scan signal and the y-scanner may be driven by a slow scan signal. That is, the scanner may direct the probe beam along the fast scan axis and slow scan axis according to a property, such as a voltage, of the fast scan signal and the slow scan signal, respectively. The fast scan signal and slow scan signal may be any suitable waveforms, such as triangular, square, stepped, and/or sinusoidal.
- Any suitable scanning protocol may be used to combine the fast scan signal and slow scan signal to produce a suitable 3D data set. For example,
FIG. 3 shows an embodiment where thefast scan signal 302 andslow scan signal 304 are each a triangular waveform (i.e., continuous scan). Thefast scan signal 302 has a higher frequency than theslow scan signal 304. As the slow scan is performed to obtain a C-scan, a plurality of fast scans are performed to produce a plurality of B-scans. In some embodiments, only one slow scan may be performed in the composite scan, as depicted inFIG. 3 . In other embodiments, a plurality of slow scans may be performed and combined into the composite scans. -
FIG. 4 shows an alternative embodiment, in which fast scansignal 402 is a triangular waveform, andslow scan signal 404 is a stepped waveform. The slow scan signal may have a number, Q, ofsteps 406 to finish one C-scan. At each step, a number, P, of B-scans may be acquired. Accordingly, the number, M, of B-scans in a C-scan may be P multiplied by Q (M=P×Q). In some embodiments, P may be at east two. - In various embodiments, any suitable imaging algorithm may be applied to the 3D data set to produce the image of the sample, such as a blood flow image. The imaging algorithm may be applied in the slow scan axis (i.e., the y-direction). The imaging algorithm may be designed to extract information on a moving component of the sample. Applying the imaging algorithm in the slow scan axis rather than the fast scan axis provides higher sensitivity, thereby allowing imaging of relatively slow movement, such as capillary blood flow within dermis which has typical flow in the range of about 100 to 900 μm/s in a resting condition, and even slower in diseased states.
- In some embodiments, the imaging algorithm may include a differentiation operation in the y-direction, followed by an absolute operation. In some embodiments, the differentiation and absolute operations may be performed on the complex function I′(xi, yj, zk) (i.e., I′(xi,yj,zk)=|I(xi,yj,zk)−I(xi,yj-1,zk)| for j=1, 2, . . . , M). In other embodiments, the differentiation and absolute operations may be performed on the magnitude portion A(xi, yj, zk) of the complex function I(xi, yj, zk) (i.e., I′(xi,yj,zk)=|A(xi,yj,zk)−A(xi,yj-1,zk)| for j=1, 2, . . . , M).
- In other embodiments, the imaging algorithm may include a high pass filter that may be applied to the 3D data set in the y-direction (slow scan axis) in order to separate the optical signals scattered by the moving particles from the optical signal scattered by the microstructures, i.e., the static particles. The high pass filter may use any suitable type of high pass filtering function.
- Furthermore, in some embodiments of the imaging algorithm, the prior OMAG method, as referenced above, may be applied to the 3D data set in the y-direction (slow scan axis).
- In embodiments where M is large compared to N, i.e., B-scans are densely taken in the y-direction, averaging adjacent B-scans resulting from the imaging algorithm may improve the quality of final flow images. The number of B-scans used to average can be any number larger than, or equal to, 2.
- In embodiments where the y-scanner is driven by a stepped function (as depicted in
FIG. 4 ), the imaging algorithm may be applied independently at each step. Then, the image at each step (i.e., I′(xi,yj,zk)), may be averaged to obtain a final B-scan image of the blood flow at that step. In some embodiments, the averaging may be done according to the following equation (Eq. (1)): -
I(x i ,y j ,z k)=[Σj=1 p I′(x i ,y j ,z k)]/P, where t=1, 2, . . . , Q (1) - Accordingly, the final 3D image may include a number, Q, of B-scan images in the y-direction.
- In various embodiments, the algorithms described above may be applied to the 3D data set in the spectral/frequency domain (i.e., the data are in wavelength (or wavenumber) format). Further, in various embodiments, a logarithm operation may be applied to the 3D data set.
- In various other embodiments, the algorithms described above may be applied to the 3D data set in the time/distance domain (i.e., the data are in time (or distance) format). In frequency/spectral domain optical coherence tomography, it has been observed that the time/distance domain signal and the spectral/frequency domain signal are a Fourier transform pair.
- Below are a description and results of an experiment conducted using a UHS-OMAG system, in accordance with various embodiments, as applied to imaging capillary blood flows within dermis. Ideally, an imaging tool for such an application must be able to resolve the capillary blood flows within dermis, which are normally very slow (in the range of about 100 to 900 μm/s at the resting condition, and even slower at diseased states). In addition, such tools must be able to provide depth information with an imaging resolution at a scale of capillary blood vessel (about 10 μm). Furthermore, the imaging tool must have a relatively low data acquisition time to allow in vivo use, since involuntary movement of the subject is unavoidable.
- The system setup for UHS-OMAG used in the experiment is similar to that described in R. K Wang, and L. An, “Doppler optical micro-angiography for volumetric imaging of vascular perfusion in vivo,” Opt. Express 17, 8926-8940 (2009) (hereinafter “
Article 1”), which is hereby incorporated by reference in its entirety. Here, the main parameters are briefly described. The system used a superluminescent diode as the light source, which has a central wavelength of 1310 nm and a bandwidth of 65 nm that provided an axial resolution in air of about 12 μm. In the sample arm, a 50 mm focal length objective lens was used to achieve a lateral resolution of about 16 μm. The output light from the interferometer was routed to a home-built spectrometer, which had a designed spectral resolution of about 0.141 nm that provided a detectable depth range of about 3 mm on each side of the zero delay line. The line scan rate of the camera was 47,000 per second (47 kHz). With this imaging speed, the signal to noise ratio was measured at about 85 dB with a light power on the sample of about 3 mW. - The system applied a scanning protocol, similar to
FIG. 3 (as discussed above), that was designed to achieve ultrahigh sensitive imaging to the blood flow. First, for each B-scan (i.e. x-direction scan), 128 A-lines were acquired with a spacing of about 15 μm between adjacent lines, thus covering a length of about 2 mm on the tissue. The imaging rate was 300 frames per second (fps). Note that with a 47 kHz line scan rate, the theoretical imaging rate should be 367 fps. The reduced imaging rate at 300 fps was due to data transfer limitations during the handshake between the camera and the computer. - Secondly, in the y-direction (i.e, C-scan direction), 1500 B-scans were captured over 2.0 mm on the tissue, which gave a spacing of about 1.3 μm between adjacent B-scans, equating to an oversampling factor of about 12 in the C-scan direction. The whole 3D data set was captured within 5 seconds.
- The imaging algorithm was applied on the slow axis (C-scan direction) rather than the fast axis (B-scan direction), As discussed in
Article 1, the interference signal of one B-scan captured by the CCD camera can be expressed as the following equation (Eq. (2)): -
I(t,k)=2S(k)EB [∫−∞ ∞ α(z,t)cos(2kn(t)z)dz+α(z1 )cos [2kn(t)(z1 −vt)]] (2) - where k is the wavenumber; t is the timing when an A-line was captured; ER is the light reflected from the reference mirror; S(k) is the spectral density of the light source used; n is the refractive index of tissue; z is the depth coordinate; a(z, t) is the amplitude of the back scattered light; v is the velocity of moving blood cells in a blood vessel, which is located at depth z1. The self cross-correlation between the light backscattered from different positions within the sample is not considered in Eq. (2), because the light backscattered from the sample is weak compared to the light reflected from the reference mirror. Additionally, the DC signals are not considered, because they do not contribute to useful OMAG signals.
- Prior OMAG systems used high pass filtering in the fast scanning axis, i.e. B-scan direction, to isolate the optical scattering signals between the static and moving scatters. Thus, the detect-able flow velocity, v, is determined by the time spacing, ΔtA, between the adjacent A-scans, i.e., v=λ/2nΔtA, wherein λ is the central wavelength of the light source, and n is the refractive index of the sample. If the flow velocity in a capillary is less than or equal to 100 μm/s, then it would require ΔtA to be greater than or equal to 4.7 ms for the system to have a chance to sample the blood cells flowing in the capillary. This time spacing translates into a scanning speed of about 213 A-scans per second. Therefore, the total data acquisition time to acquire a 3D capillary flow image of a tissue volume would be prohibitively long, and not ideal for in vivo imaging of capillary blood flows.
- In order to image the slow blood flow within capillary vessels while keeping the data acquisition time low, the imaging algorithm is performed in the C-scan direction (slow scan axis). In this case, Eq. (2) can still be used to represent the spectral interferogram signal captured by the system, except that the time variable, t, now corresponds to the B-scan numbers in one C-scan. With this modification, the requirement of oversampling in the B-scan direction, as was present in the prior OMAG system, is relaxed, making it possible to have a much faster B-scan imaging rate, provided that the line scan camera in the spectrometer is limited or fixed. The detectable flow velocity is determined by the time spacing, ΔtB, between adjacent B-scans, i.e., v=λ/2nΔtB. In the system setup described here, the imaging rate is 300 fps, so ΔtB is about 3.3 ms. Note that the imaging speed is at 47,000 A scans per second, in contrast to 213 A scan per second required for prior OMAG method.
- In the data processing, the imaging algorithm first takes a differential operation on the captured B-scan spectral interferograrns along the C-scan direction, as in the following equation (Eq. (3)):
-
I flow(t i ,k)=I(t i ,k),−I(t i-1 ,k), i=1, 2, 3 . . . 1500 (3) - where i represents the index of the B-scans in the C-scan direction. The differential operation suppresses the optical scattering signals from the static elements within scanned tissue volume. Alternatively, high pass filtering may be used. Then, a fast Fourier transform (FFT) is applied upon every wavenumber k (t is now constant) of Eq. (2) to obtain a depth-resolved OMAG flow image with ultrahigh sensitivity to the flow.
- The minimum detectable blood flow is determined by the system phase noise floor, which can be expressed by the intensity signal to noise ratio, S, of the OMAG/OCT system by σΔφ 2=1/S. Thus, with the system signal to noise ratio at 85 dB, the minimum detectable flow velocity would be about 4.0 μm/s. However, if a blood cell moves at 4 μm/s, the system described here would not provide a continuous trajectory for this blood cell in the 3D OMAG flow image, i.e., the trajectory would be seen as a broken line.
- Because the system is very sensitive to movement, the bulk motion of the sample may seriously degrade the final image result if the imaging algorithm is directly applied. To solve this problem, the phase compensation method described in
Article 1 is applied to the raw interference signal before applying the OMAG algorithm. - To test the performance of the UHS-OMAG system described above for imaging blood flow, the system was tested on the skin located on the backside of a hand of a male volunteer. For comparison, the traditional OMAG and phase resolved optical Doppler tomography (PRODT) cross-sectional flow images were also obtained. For these methods, the system captured 2000 A-scans over 2 mm at an imaging speed of 31,000 A-scans per second in order to fulfill the oversampling requirement for these previous methods.
- The results are shown in
FIGS. 5A-F . The images in the top row (FIGS. 5A-C ) are from the UHS-OMAG system while those in the bottom row (FIGS. 5D-F ) are from the conventional prior OMAG system.FIGS. 5A and 5D are the Fourier domain optical coherence tomography (FDOCT) structural images obtained from the captured interferograms for the UHS-OMAG system and prior OMAG system, respectively. Although they are similar, they are not exactly the same due to the small subject movement when switching the system among the different approaches. However, it is sufficient to provide a fair comparison of their ability to extract slow flow information.FIG. 5B shows the image fromFIG. 5A after it is processed by the UHS-OMAG system.FIG. 5C shows the PRODT image of the UHS-OMAG system, which is based on the phase difference between adjacent B-scans. Similarly,FIG. 5E shows the image fromFIG. 50 after it is processed by the prior OMAG system.FIG. 5F shows the conventional PRODT image of the prior OMAG system, which is based on the phase difference between adjacent A-scans in one B-scan. InFIGS. 50 and 5F , the phase differences are calculated only when the structural signal is 15 dB above the noise floor. - It is apparent that the UHS-OMAG approach outperforms the other methods.
FIG. 5B shows blood flows within the papillary dermis (indicated by white arrows), where only capillary blood vessels are present, as well as the blood flows within reticular dermis (indicated by red arrows), where both the capillary and larger blood vessels are present. By calculating the phase differences between adjacent B-scans, the blood flow velocities within the capillaries are within the reach of the UHS-OMAG (e.g., see the white arrows inFIG. 50 ). - Because the conventional OMAG requires oversampling in the fast scanning direction (i.e., B-scan), it is not sensitive to the slow blood flow within the capillaries, which are normally below 100 μm/s. Accordingly, the conventional PRODT approach totally failed in imaging any of blood vessels, as seen in
FIG. 5F . It should be noted that there is a global noise ‘flow’ background in the UHS-OMAG flow image, e.g., inFIG. 5B , which may be caused by some ‘non-moving’ scatters, such as global motion, etc. In this case, de-noising filters may be used to further enhance the UHS-OMAG flow imaging quality. - In order to test whether the flow sensitivity of UHS-OMAG approaches the system phase-noise floor (in this case about 4 μm/s as stared above), a second set of images was obtained using a highly scattering flow phantom as the imaging target. The phantom was made of gelatin mixed with about 1% milk to simulate the background optical heterogeneity of the tissue. In making this background tissue, precaution was taken so that the mixed gel was well solidified to minimize the possible Brownian motion of particles in the background. A capillary tube with an inner diameter of about 400 μm was submerged in this background tissue and about 2% TiO2 particle solution was flowing in it that was controlled by a precision syringe pump. Although such setup can precisely control the flow velocity in the capillary tube, a flow speed as low as about 4 μm/s is difficult to provide. This is especially true considering that if the flow is stopped, the Brownian motion of particles is unavoidable in the capillary tube. With this experimental condition, the motion speed of particles due to Brownian motion would be randomly distributed within a range of several tens of microns per second. Due to these reasons, the experiments tested the UHS-OMAG system's ability to measure the Brownian motion of particles.
- In the experiments, the capillary rube was made almost perpendicular to the incident sample beam to avoid free fall of the scattering particles within the tube. The imaging results are shown in
FIGS. 6A-B .FIG. 6A is the OMAG/OCT microstructural image of the flow phantom, whileFIG. 6B is the corresponding UHS-OMAG flow image. From this result, it is clear that UHS-OMAG is able to image the particle movements due to Brownian motion with almost no signals detected in the background region. - To examine in more detail, the phase-resolved technique was applied to the adjacent B-scans of the UHS-OMAG flow images to provide the velocity image of the flow phantom above. The result is shown in
FIG. 7A , where it can be seen that the velocity values in the background region are low while those within the capillary tube are contrasted out primarily due to the Brownian motion of the particles.FIG. 7B shows a plot of the calculated velocities across the center of the capillary tube at the position marked as the horizontal line inFIG. 7A , where the dashed box indicates the position of the capillary lumen. The velocity values of particle movements ranged from approximately −50 to 100 μm/s at this cross-line position. The standard deviation of the values outside the dashed-box region was evaluated to be about 4.5 μm/s, which is close to the theoretical value of about 4 μm/s. This experiment concluded that the UHS-OMAG system is sensitive to the flow as low as about 4 μm/s for the system setup used in this study. - The conventional PRODT image was also obtained of the same phantom. In doing so, the system imaging rate was set at 31,000 A-scans per second. Additionally, the A-line density across the B-scan of about 2.5 mm was set at 4000, equating to a spacing between adjacent A-scans of about 0.625 μm. The corresponding results are given in
FIGS. 7C and 7D , respectively, showing that PRODT totally failed to image the Brownian motion of the particles under the current experimental setup. Note that the standard deviation of velocity values shown inFIG. 7D was about 180 μm/s, thus it is not surprising that PRODT is not able to achieve satisfactory imaging performance. - After evaluation of the UHS-OMAG system's imaging performance, a subsequent experiment shows its capability to image the capillary blood flows within dermis in three dimensions.
FIG. 8A shows a schematic drawing of the blood vessel system of the human skin, in which an interconnected network of vessels is characterized by regular structures on all levels. The human skin is composed of the cutis and the subcutis (hypodermis [HD]). The cutis is further divided into the epidermis (EP) and the dermis (DR). The dermis and subcutis are pervaded with a complex system of blood vessels, while the epidermis is free of vessels. A superficial network comprises the interface between papillary (PD) and reticular dermis (RD), while a lower network is located on the border between dermis and subcutis. Vertical vessels connect both networks and thus make it complete. In the diagram, arteries are shown in red and veins in blue. To show whether the ultrahigh sensitive OMAG is able to image the blood flow within the patent blood vessels as described above, 3D blood flow images were acquired over the palm of a healthy volunteer, as shown inFIG. 8B , where the black box indicates the scanning area (about 2×2 mm2). - The 3D OMAG imaging result of blood vessel networks is shown in
FIG. 8C , shown together with the 3D micro-structural image.FIG. 8D to shows a cross-sectional view, where the blood flows within blood vessel systems within the skin are clearly delineated. Because the UHS-OMAG sensitivity is as low as about 4 μm/s, even the dynamics of sweat glands are imaged. - The projection views at different land-mark depths are shown in
FIGS. 9A-D .FIG. 9A gives the projection view at the depths from 400 to 450 μm, which corresponds to the papillary dermis where the capillary vessels are dense (e.g., indicated by the arrows).FIG. 9B shows the depths from 450 to 650 μm, where the vessels are almost vertical that connect vessel networks between papillary dermis and reticular dermis, seen as the bright spots inFIG. 9B . The blood vessel network in reticular dermis (650 to 780 μm) and hypodermis (780 to 1100 μm) are shown inFIGS. 9C and 9D , respectively. As shown, the vessel diameter is smaller in the reticular dermis than in the hypodermis. These observations from ultrahigh sensitive OMAG are almost identical to that described in the literature, demonstrating the power of the ultrahigh sensitive OMAG in the investigations of pathological conditions in dermatology. - These experiments demonstrated an ultrahigh sensitive OMAG system to image the volumetric microcirculation within the human skin. It was achieved by applying the OMAG algorithm along the slow scan axis (i.e., the C-scan direction), as opposed to the fast axis (i.e., the B scan direction). Comparing with the conventional OMAG flow image, the UHS-OMAG method delivers much better performance to extract slow flow information. Detailed 3D microvascular images obtained from the human skin by the UHS-OMAG system are comparable to those described in the standard textbook. Therefore, the ultrahigh sensitive OMAG may have great value in future clinical investigations of pathological conditions in human skin.
- In another example, the UHS-OMAG system was used to obtain depth-resolved images of the capillary networks in the retina and choroid of the human eye.
FIG. 10 shows the setup of the UHS-OMAG system 1000 used to obtain the images. The UHS-OMAG system 1000 is similar to that described in L. An, and R. K. Wang, Optics Express, 16, 11438-11452 (2008) (hereinafter “Article 2”). The UHS-OMAG system 1000 used includes alight source 1002 to produce a probe beam with a wavelength centered at 842 nm and a bandwidth of 46 nm that provided an axial resolution in air of about 8 μm. In asample arm 1004, the light was delivered onto ahuman eye 1006 by acollimator 1008, anobjective lens 1010, and anocular lens 1012. The sample arm further includes an x-scanner 1014 and a y-scanner 1016 for scanning thehuman eye 1006 along the x-axis and y-axis, respectively. - In a
reference arm 1018, a 20mm water cell 1020 was used to compensate for the dispersion caused by the eye. A reference beam from the light source was reflected off areference mirror 1022 to provide a reference light. An interferogram between the reference light and the light backscattered from the sample was sent to a home-builtultrafast spectrometer 1024 via anoptical circulator 1026. The spectrometer included acollimator 1027, a transmission grating 1028 (with 1200 lines/mm), acamera lens 1030 with a focal length of 100 mm, and a 1024 element linescan CMOS detector 1032 capable of a 140 kHz line scan rate. The spectral resolution of the designedspectrometer 1024 was about 0.055 nm, which provided an imaging depth of about 3 mm in air. The system sensitivity was about 90 dB measured with about 900 μW power of light incident at the object and an exposure time of 6.9 μs with the time interval, ΔtA, between A-scans equal to about 7.4 μs. Thesystem 1000 further includespolarization controllers - The UHS-
OMAG system 1000 employed a scanning protocol, similar toFIG. 4 (as discussed above), designed to achieve ultrahigh sensitive imaging of the flow. Firstly, the x-scanner 1012 was driven by a 400 Hz sawtooth waveform, meaning that the imaging rate is 400 frames/sec (fps). The duty cycle for acquiring each B-scan (i.e., x-direction scan) was about 75%, in which we acquired 256 A-lines with a spacing of about 12 μm between adjacent lines that covered a size of about 3 mm on the retina. Secondly, the y-scanner was driven by a step function, and the entire C-scan includes 150 steps, with a spacing between adjacent steps of about 20 μm. In each step, 8 repeated B-scans were acquired. Accordingly, it required 3 seconds to acquire one 3D data set, covering an area of about 3×3 mm on the retina. - An imaging algorithm was then applied on the 3D data set along the slow scan axis, i.e., C-scan direction. In this case, the detectable flow velocity is determined by the time spacing, ΔtB, between adjacent B-scans. Because the imaging speed was 400 fps, ΔtB=2.5 ms, which would be sufficient to image the slow flows in capillaries (see above). Finally, the calculated OMAG signals at each step were collapsed into one B-scan through ensemble-averaging (e.g., Eq. (1)), resulting in 150 B-scans to form the final C-scan image, i.e. 3D OMAG blood flow distribution.
- To demonstrate the performance of the US-OMAG system, the system was used to obtain images on healthy volunteers. To reduce the eye and head movement, the volunteer was asked to steer at a fixed position during the experiment.
FIGS. 11A-D show the in vivo imaging results produced by one volume dataset captured at the macular region towards the optic nerve head.FIG. 11A shows a typical cross-sectional image (B-scan) within the OMAG structural volume, which is identical to the conventional OCT image where the typical morphological features within retina and choroid are visualized.FIG. 11B gives the corresponding blood flow image obtained from the imaging algorithm, where the capillary flows within the cross-section of retina are abundant (indicated by arrows), as well as the blood flow signals in the choroid. - Because of the depth-resolved nature, the blood flows in the retina can be separated from those in the choroid. To do this, the segmentation algorithm described in
Article 2 was first used to identify the retinal pigment epithelium (RPE) layer. Then, the flow signals from the retinal vessels were identified as the OMAG flow signals 50 μm above the RPE layer (to exclude the signals from photoreceptor inner and outer segments), while those below the RPE layer were identified as the choroidal vessels. The segmentation resulted in two volumetric flow images, one for retina and another for choroid, which are annotated inFIG. 11B by the labels of R and C, respectively. Finally, the maximum amplitude projection (MAP) was performed on each segmented retinal and choroidal volumes, giving the blood flow distribution maps shown inFIGS. 11C (retina) and 11D (choroid), respectively. InFIG. 11C , there is seen a ring of blood vessels in the macular area around an avascular zone about 800 μm in diameter, denoting the fovea. This observation is in excellent agreement with the standard retinal pathology. - According to the literature, the retina consists of three layers of capillary networks: the radial peripapillary capillaries (RPCs, R1) and an inner (R2) and an outer layer of capillaries (R3). The RPCs are the most superficial layer of capillaries lying in the inner part of the nerve fiber layer. The inner capillaries lie in the ganglion cell layers under and parallel to the RPCs. The outer capillary network runs from the inner plexiform to the outer plexiform through the inner nuclear layers. However, the choroidal vessels in the macular region are not specialized like those in the retina. The arteries pierce the sclera around the optic nerve and fan out to form three vascular layers in the choroid: inner capillary bed (near RPE layer, C1), medial arterioles and venules (C2) and outer arteries and veins (C3). With these descriptions as the reference, the OMAG blood flow signals obtained from the retinal and choroidal layers were further separated, with respect to the RPE layer. The segmentation is pictorially illustrated in
FIG. 11B , with the labels of R1 (about 425 μm above RPE), R2 (between 300 and 425 μm above RPE) and R3 (between 50 and 300 μm above RPE) for retina, and C1 (from 0 to 70 μm below RPE), C2 (from 70 to 200 μm below RPE) and C3 (beyond 200 μm below RPE) for choroid, respectively. After segmentation, the blood flow MAPs within each land-marked depth are shown inFIGS. 12A-F (FIGS. 12A-C for retina andFIGS. 12D-F for choroid). The results are correlated well with the descriptions found in the literature. - The OMAG system running at 400 fps requires about 3 seconds to acquire one 3D blood flow image representing about 3×3 mm2 area on the retina. At this speed, the effects of subject movement on the final results are clearly noticeable (the horizontal lines in
FIGS. 11A-D and 12A-F are the motion artifacts). There may be several solutions to amend this problem, such as 1) the phase compensation algorithm developed inArticle 2 may be used to minimize the motion artifacts before the imaging algorithm is applied. In doing so, however, it would inevitably increase the computational load required to obtain the meaningful 3D blood flow images because of the complexity of the compensation algorithms, and/or 2) the imaging speed may be further increased to minimize the motion artifacts. The current system speed is limited by the CMOS camera used in the spectrometer with a maximum line rate of 140 kHz. However, a report has shown that a 4096 element CMOS camera is capable of a line rate more than 240 kHz. Therefore, it would be expected that if this CMOS camera is used in the system, the motion artifacts may be mitigated. - As shown above, the UHS-OMAG system may be capable of imaging detailed ocular perfusion distributions within retina and choroid. Applying the imaging algorithm on the slow scanning axis, the system is sensitive to the ocular capillary flows. In addition, due to the depth-resolved nature, the system may be able to provide detailed micro-circulation within different land-marked depths, the results of which are in an excellent agreement with those described in the literature. The demonstrated superior imaging results show promise for future clinical applications for UHS-OMAG in ophthalmology.
- Although certain embodiments have been illustrated and described herein, it will be appreciated by those of ordinary skill in the art that a wide variety of alternate and/or equivalent embodiments or implementations calculated to achieve the same purposes may be substituted for the embodiments shown and described without departing from the scope. Those with skill in the art will readily appreciate that embodiments may be implemented in a very wide variety of ways. This application is intended to cover any adaptations or variations of the embodiments discussed herein. Therefore, it is manifestly intended that embodiments be limited only by the claims and the equivalents thereof.
Claims (18)
1-20. (canceled)
21. A method of imaging, comprising:
receiving a three dimensional (3D) data set generated by one or more spectral interference signals obtained by scanning a sample with a probe beam from a light source, wherein the one or more spectral interference signals were generated by:
performing a plurality of fast scans (B-scans) on a fast scan axis, each fast scan comprising a plurality of A-scans, and
performing one or more slow scans (C-scans), contemporaneously with the fast scans, on a slow scan axis, the slow scan axis being orthogonal to the fast scan axis, and each slow scan comprising a plurality of fast scans; and
applying an imaging algorithm to the 3D data set along the slow scan axis to produce at least one image of the sample.
22. The method of claim 21 , wherein the fast scans are driven by a first triangular waveform and the slow scans are driven by a second triangular waveform, the first triangular waveform having a higher frequency than the second triangular waveform.
23. The method of claim 21 , wherein the slow scans are driven by a stepped waveform having a plurality of steps, and a plurality of fast scans are performed at each step of the stepped waveform.
24. The method of claim 21 , wherein the imaging algorithm is configured to operate along the slow scan axis to separate a moving component of the sample from a structural component of the sample.
25. The method of claim 24 , wherein the imaging algorithm comprises high pass filtering along the slow scan axis in the 3D data set.
26. The method of claim 21 , wherein the image of the sample comprises an image of a blood vessel network.
27. The method of claim 21 further comprises applying a phase-resolved technique to the adjacent B-scans of the image of the sample to provide blood flow velocity information.
28. The method of claim 21 wherein the method is applied in vivo.
29. A system for in vivo imaging, comprising:
an optical micro-angiography (OMAG) apparatus; and
one or more processors coupled to the OMAG apparatus and adapted to cause the OMAG apparatus to:
scan a sample with a probe beam from a light source;
perform a plurality of fast scans (B-scans) on a fast scan axis, each fast scan comprising a plurality of A-scans;
perform one or more slow scans (C-scans), contemporaneously with the fast scans, on a slow scan axis, the slow scan axis being orthogonal to the fast scan axis, and each slow scan comprising a plurality of fast scans;
detect one or more spectral interference signals from the sample during the scanning to generate a three dimensional (3D) data set; and
apply an imaging algorithm to the 3D data set along the slow scan axis to produce an image of the sample.
30. The system of claim 29 , wherein the OMAG apparatus includes an x-scanner to perform the plurality of fast scans and a y-scanner to perform the one or more slow scans, and the x-scanner is driven by a fast scan signal, and the y-scanner is driven by a slow scan signal.
31. The system of claim 30 , wherein the fast scan signal comprises a first triangular waveform and slow scan signal comprises a second triangular waveform, the first triangular waveform having a higher frequency than the second triangular waveform.
32. The system of claim 30 , wherein the slow scan signal comprises a stepped waveform having a plurality of steps, and the fast scan signal is configured to cause the x-scanner to perform a plurality of fast scans during each step of the stepped waveform.
33. The system of claim 29 , wherein the imaging algorithm is configured to operate along the slow scan axis to separate a moving component of the sample from a structural component of the sample.
34. The system of claim 33 , wherein the imaging algorithm comprises high pass filtering along the slow scan axis in the 3D data set.
35. The system of claim 29 , wherein the image of the sample comprises an image of a blood vessel network.
36. The system of claim 29 , further comprising applying the phase-resolved technique to the adjacent B-scans of the image of the sample to provide the blood flow velocity information.
37. The system of claim 29 , wherein the OMAG apparatus is a Fourier domain optical coherence tomography (FD-OCT) apparatus.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14/658,558 US20150211838A1 (en) | 2009-05-04 | 2015-03-16 | Method and apparatus for ultrahigh sensitive optical microangiography |
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17522909P | 2009-05-04 | 2009-05-04 | |
US30240910P | 2010-02-08 | 2010-02-08 | |
PCT/US2011/024069 WO2011097631A2 (en) | 2010-02-08 | 2011-02-08 | Method and apparatus for ultrahigh sensitive optical microangiography |
US201213577857A | 2012-08-08 | 2012-08-08 | |
US14/658,558 US20150211838A1 (en) | 2009-05-04 | 2015-03-16 | Method and apparatus for ultrahigh sensitive optical microangiography |
Related Parent Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2011/024069 Continuation WO2011097631A2 (en) | 2009-05-04 | 2011-02-08 | Method and apparatus for ultrahigh sensitive optical microangiography |
US13/577,857 Continuation US9013555B2 (en) | 2009-05-04 | 2011-02-08 | Method and apparatus for ultrahigh sensitive optical microangiography |
Publications (1)
Publication Number | Publication Date |
---|---|
US20150211838A1 true US20150211838A1 (en) | 2015-07-30 |
Family
ID=44356115
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/577,857 Active US9013555B2 (en) | 2009-05-04 | 2011-02-08 | Method and apparatus for ultrahigh sensitive optical microangiography |
US14/658,558 Abandoned US20150211838A1 (en) | 2009-05-04 | 2015-03-16 | Method and apparatus for ultrahigh sensitive optical microangiography |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US13/577,857 Active US9013555B2 (en) | 2009-05-04 | 2011-02-08 | Method and apparatus for ultrahigh sensitive optical microangiography |
Country Status (6)
Country | Link |
---|---|
US (2) | US9013555B2 (en) |
EP (1) | EP2538833A4 (en) |
JP (1) | JP2013518695A (en) |
CN (1) | CN103002794B (en) |
TW (1) | TW201140226A (en) |
WO (1) | WO2011097631A2 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11497397B2 (en) | 2018-03-27 | 2022-11-15 | Topcon Corporation | Ophthalmologic apparatus |
EP4218544A1 (en) * | 2021-11-12 | 2023-08-02 | Topcon Corporation | Apparatus and method for imaging an eye |
Families Citing this family (59)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2010129494A2 (en) | 2009-05-04 | 2010-11-11 | Oregon Health & Science University | Method and apparatus for quantitative imaging of blood perfusion in living tissue |
US9033510B2 (en) | 2011-03-30 | 2015-05-19 | Carl Zeiss Meditec, Inc. | Systems and methods for efficiently obtaining measurements of the human eye using tracking |
US8433393B2 (en) | 2011-07-07 | 2013-04-30 | Carl Zeiss Meditec, Inc. | Inter-frame complex OCT data analysis techniques |
JP6061554B2 (en) | 2012-01-20 | 2017-01-18 | キヤノン株式会社 | Image processing apparatus and image processing method |
JP6146951B2 (en) | 2012-01-20 | 2017-06-14 | キヤノン株式会社 | Image processing apparatus, image processing method, photographing apparatus, and photographing method |
ES2628937T3 (en) * | 2012-02-03 | 2017-08-04 | Oregon Health & Science University | Obtaining optical flow images in vivo |
US9357916B2 (en) | 2012-05-10 | 2016-06-07 | Carl Zeiss Meditec, Inc. | Analysis and visualization of OCT angiography data |
JP6316298B2 (en) * | 2012-09-10 | 2018-04-25 | オレゴン ヘルス アンド サイエンス ユニバーシティ | Quantification of local circulation by OCT angiography |
JP2014113207A (en) * | 2012-12-06 | 2014-06-26 | Tomey Corporation | Tomographic apparatus and tomographic image processing method |
WO2014160116A1 (en) * | 2013-03-13 | 2014-10-02 | University Of Washington Through Its Center For Commercialization | Methods and systems for imaging tissue motion using optical coherence tomography |
WO2014168930A1 (en) * | 2013-04-09 | 2014-10-16 | University Of Washington Through Its Center For Commercialization | Methods and systems for determining hemodynamic properties of a tissue |
CN103211588B (en) * | 2013-05-06 | 2014-10-22 | 王毅 | Blood flow velocity imaging method based on transverse velocity model |
JP6278295B2 (en) * | 2013-06-13 | 2018-02-14 | 国立大学法人 筑波大学 | Optical coherence tomography device for selectively visualizing and analyzing choroidal vascular network and its image processing program |
US9778021B2 (en) | 2013-08-29 | 2017-10-03 | Carl Zeiss Meditec, Inc. | Evaluation of optical coherence tomographic data prior to segmentation |
US10307056B2 (en) * | 2013-12-05 | 2019-06-04 | Bioptigen, Inc. | Systems and methods for quantitative doppler optical coherence tomography |
US9844320B2 (en) * | 2014-01-29 | 2017-12-19 | University Of Rochester | System and method for observing an object in a blood vessel |
WO2015120055A1 (en) * | 2014-02-04 | 2015-08-13 | University Of Southern California | Optical coherence tomography (oct) system with phase-sensitive b-scan registration |
EP3102091B1 (en) | 2014-02-04 | 2019-01-16 | University of Southern California | Optical coherence tomography (oct) system with improved motion contrast |
US9778018B2 (en) | 2014-02-14 | 2017-10-03 | Carl Zeiss Meditec, Inc. | Swept source interferometric imaging systems and methods |
US10398302B2 (en) | 2014-05-02 | 2019-09-03 | Carl Zeiss Meditec, Inc. | Enhanced vessel characterization in optical coherence tomograogphy angiography |
CN105310677B (en) * | 2014-08-01 | 2018-01-12 | 深圳市斯尔顿科技有限公司 | A kind of device and method for measuring blood flow |
WO2015172322A1 (en) * | 2014-05-14 | 2015-11-19 | 深圳市斯尔顿科技有限公司 | Blood flow measuring device and method |
WO2015172581A1 (en) * | 2014-05-14 | 2015-11-19 | 深圳市斯尔顿科技有限公司 | Device and method for measuring blood flow |
US10405793B2 (en) | 2014-07-01 | 2019-09-10 | University Of Washington | Systems and methods for in vivo visualization of lymphatic vessels with optical coherence tomography |
US20160007854A1 (en) * | 2014-07-09 | 2016-01-14 | Physical Sciences, Inc. | Apparatus and Method for Assessment of Interstitial Tissue |
US9759544B2 (en) | 2014-08-08 | 2017-09-12 | Carl Zeiss Meditec, Inc. | Methods of reducing motion artifacts for optical coherence tomography angiography |
JP6798095B2 (en) * | 2014-12-02 | 2020-12-09 | 株式会社ニデック | Optical coherence tomography equipment and control programs used for it |
CN104545872B (en) * | 2015-01-12 | 2017-04-19 | 南京理工大学 | Method and device for reconstructing three-dimensional micro blood flow distribution on basis of linearly dependent coefficients |
US9700206B2 (en) | 2015-02-05 | 2017-07-11 | Carl Zeiss Meditec, Inc. | Acquistion and analysis techniques for improved outcomes in optical coherence tomography angiography |
US10368734B2 (en) | 2015-02-19 | 2019-08-06 | Carl Zeiss Meditec, Inc. | Methods and systems for combined morphological and angiographic analyses of retinal features |
US11284792B2 (en) | 2015-03-06 | 2022-03-29 | University Of Washington | Methods and systems for enhancing microangiography image quality |
US9984459B2 (en) | 2015-04-15 | 2018-05-29 | Kabushiki Kaisha Topcon | OCT angiography calculation with optimized signal processing |
CN104825148B (en) * | 2015-05-27 | 2017-11-14 | 深圳市斯尔顿科技有限公司 | A kind of device and method for measuring blood vessel flow |
JP6843125B2 (en) | 2015-09-24 | 2021-03-17 | カール ツァイス メディテック インコーポレイテッドCarl Zeiss Meditec Inc. | High-sensitivity flow visualization method |
US10402965B1 (en) | 2015-11-12 | 2019-09-03 | Carl Zeiss Meditec, Inc. | Systems and methods for reducing artifacts in OCT angiography images |
CN107049242B (en) * | 2017-01-24 | 2023-08-29 | 青岛大学 | Scanning type human body microvascular ultrastructural three-dimensional imaging system |
US10441164B1 (en) | 2017-03-09 | 2019-10-15 | Carl Zeiss Meditec, Inc. | Correction of decorrelation tail artifacts in a whole OCT-A volume |
CN107862724B (en) * | 2017-12-01 | 2021-10-01 | 中国医学科学院生物医学工程研究所 | Improved microvascular blood flow imaging method |
US11751762B2 (en) * | 2018-12-19 | 2023-09-12 | Topcon Corporation | Method and apparatus for low coherence interferometry |
CN109691977B (en) * | 2018-12-29 | 2023-11-28 | 佛山科学技术学院 | Non-confocal adaptive optical imaging system |
US20220084210A1 (en) | 2019-02-08 | 2022-03-17 | Carl Zeiss Meditec, Inc. | Segmentation and classification of geographic atrophy patterns in patients with age related macular degeneration in widefield autofluorescence images |
CN113396440A (en) | 2019-02-14 | 2021-09-14 | 卡尔蔡司医疗技术公司 | System for OCT image conversion and ophthalmic image denoising and neural network thereof |
EP3941333A1 (en) | 2019-03-20 | 2022-01-26 | Carl Zeiss Meditec, Inc. | A patient tuned ophthalmic imaging system with single exposure multi-type imaging, improved focusing, and improved angiography image sequence display |
JP7202569B2 (en) * | 2019-05-15 | 2023-01-12 | 花王株式会社 | Imaging and analysis method of skin capillaries |
US20220257111A1 (en) | 2019-08-01 | 2022-08-18 | Carl Zeiss Meditec, Inc. | Ophthalmic imaging with k-mirror scanning, efficient interferometry, and pupil alignment through spatial frequency analysis |
JP2021037253A (en) * | 2019-08-27 | 2021-03-11 | 株式会社トプコン | Laminate, inspection device and model eye |
EP4025114A1 (en) | 2019-09-06 | 2022-07-13 | Carl Zeiss Meditec, Inc. | Machine learning methods for creating structure-derived visual field priors |
WO2021122562A1 (en) | 2019-12-18 | 2021-06-24 | Carl Zeiss Meditec Ag | Personalized patient interface for ophthalmic devices |
CN115135228A (en) | 2020-02-21 | 2022-09-30 | 卡尔蔡司医疗技术公司 | OCT suspension ligament imaging |
US20230091487A1 (en) | 2020-03-30 | 2023-03-23 | Carl Zeiss Meditec, Inc. | Correction of flow projection artifacts in octa volumes using neural networks |
JP7488090B2 (en) | 2020-04-08 | 2024-05-21 | 株式会社トプコン | Model eye and ophthalmic device |
EP4143781A1 (en) | 2020-04-29 | 2023-03-08 | Carl Zeiss Meditec, Inc. | Oct en face pathology segmentation using channel-coded slabs |
WO2021219773A1 (en) | 2020-04-29 | 2021-11-04 | Carl Zeiss Meditec, Inc. | Real-time ir fundus image tracking in the presence of artifacts using a reference landmark |
EP4143782A1 (en) | 2020-04-30 | 2023-03-08 | Carl Zeiss Meditec, Inc. | Bruch's membrane segmentation in oct volume |
JP2023551214A (en) | 2020-11-30 | 2023-12-07 | カール ツァイス メディテック インコーポレイテッド | Quality map for optical coherence tomography angiography |
US20240095876A1 (en) | 2020-12-04 | 2024-03-21 | Carl Zeiss Meditec, Inc. | Using multiple sub-volumes, thicknesses, and curvatures for oct/octa data registration and retinal landmark detection |
US20240081638A1 (en) | 2021-02-01 | 2024-03-14 | Carl Zeiss Meditec, Inc. | Micro-bench oct design |
WO2022180227A1 (en) | 2021-02-26 | 2022-09-01 | Carl Zeiss Meditec, Inc. | Semi-supervised fundus image quality assessment method using ir tracking |
WO2023126340A1 (en) | 2021-12-27 | 2023-07-06 | Carl Zeiss Meditec, Inc. | Fluid tracking in wet amd patients using thickness change analysis |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4751929A (en) * | 1985-01-23 | 1988-06-21 | Matsushita Electric Industrial Co., Ltd. | Ultrasonic bloodstream diagnostic apparatus with dual displays of velocity profiles and average flow velocity |
Family Cites Families (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4131023A (en) * | 1976-03-04 | 1978-12-26 | Rca Corporation | Pulse-echo ultrasonic-imaging display system |
US4131022A (en) * | 1976-03-04 | 1978-12-26 | Rca Corporation | Pulse-echo ultrasonic-imaging display system |
US6549801B1 (en) | 1998-06-11 | 2003-04-15 | The Regents Of The University Of California | Phase-resolved optical coherence tomography and optical doppler tomography for imaging fluid flow in tissue with fast scanning speed and high velocity sensitivity |
GB2363306B (en) * | 2000-05-05 | 2002-11-13 | Acoustical Tech Sg Pte Ltd | Acoustic microscope |
US7113818B2 (en) | 2002-04-08 | 2006-09-26 | Oti Ophthalmic Technologies Inc. | Apparatus for high resolution imaging of moving organs |
US7016048B2 (en) | 2002-04-09 | 2006-03-21 | The Regents Of The University Of California | Phase-resolved functional optical coherence tomography: simultaneous imaging of the stokes vectors, structure, blood flow velocity, standard deviation and birefringence in biological samples |
JP2003310689A (en) | 2002-04-18 | 2003-11-05 | Marutaka Co Ltd | Portable massager |
WO2004043245A1 (en) | 2002-11-07 | 2004-05-27 | Pawel Woszczyk | A method of fast imaging of objects by means of spectral optical coherence tomography |
JP2006516739A (en) | 2003-01-24 | 2006-07-06 | ザ・ジェネラル・ホスピタル・コーポレイション | System and method for identifying tissue using a low coherence interferometer |
US7359062B2 (en) | 2003-12-09 | 2008-04-15 | The Regents Of The University Of California | High speed spectral domain functional optical coherence tomography and optical doppler tomography for in vivo blood flow dynamics and tissue structure |
US7145661B2 (en) | 2003-12-31 | 2006-12-05 | Carl Zeiss Meditec, Inc. | Efficient optical coherence tomography (OCT) system and method for rapid imaging in three dimensions |
US7442918B2 (en) * | 2004-05-14 | 2008-10-28 | Microvision, Inc. | MEMS device having simplified drive |
JP2008542758A (en) * | 2005-05-31 | 2008-11-27 | ザ ジェネラル ホスピタル コーポレイション | System, method and apparatus capable of using spectrally encoded heterodyne interferometry for imaging |
US7847949B2 (en) * | 2005-09-29 | 2010-12-07 | The General Hospital Corporation | Method and apparatus for optical imaging via spectral encoding |
WO2008002839A2 (en) | 2006-06-26 | 2008-01-03 | California Institute Of Technology | Dynamic motion contrast and transverse flow estimation using optical coherence tomography |
US8180134B2 (en) | 2006-09-26 | 2012-05-15 | Oregon Health & Science University | In vivo structural and flow imaging |
US20090225324A1 (en) * | 2008-01-17 | 2009-09-10 | The General Hospital Corporation | Apparatus for providing endoscopic high-speed optical coherence tomography |
JP4627782B2 (en) * | 2008-03-05 | 2011-02-09 | 株式会社日立ハイテクノロジーズ | Edge detection method and charged particle beam apparatus |
WO2010129494A2 (en) | 2009-05-04 | 2010-11-11 | Oregon Health & Science University | Method and apparatus for quantitative imaging of blood perfusion in living tissue |
-
2011
- 2011-02-08 TW TW100104073A patent/TW201140226A/en unknown
- 2011-02-08 JP JP2012552927A patent/JP2013518695A/en not_active Withdrawn
- 2011-02-08 EP EP11740531.6A patent/EP2538833A4/en not_active Ceased
- 2011-02-08 WO PCT/US2011/024069 patent/WO2011097631A2/en active Application Filing
- 2011-02-08 CN CN201180018071.9A patent/CN103002794B/en not_active Expired - Fee Related
- 2011-02-08 US US13/577,857 patent/US9013555B2/en active Active
-
2015
- 2015-03-16 US US14/658,558 patent/US20150211838A1/en not_active Abandoned
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4751929A (en) * | 1985-01-23 | 1988-06-21 | Matsushita Electric Industrial Co., Ltd. | Ultrasonic bloodstream diagnostic apparatus with dual displays of velocity profiles and average flow velocity |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11497397B2 (en) | 2018-03-27 | 2022-11-15 | Topcon Corporation | Ophthalmologic apparatus |
EP4218544A1 (en) * | 2021-11-12 | 2023-08-02 | Topcon Corporation | Apparatus and method for imaging an eye |
Also Published As
Publication number | Publication date |
---|---|
US20120307014A1 (en) | 2012-12-06 |
EP2538833A2 (en) | 2013-01-02 |
TW201140226A (en) | 2011-11-16 |
EP2538833A4 (en) | 2014-03-12 |
US9013555B2 (en) | 2015-04-21 |
JP2013518695A (en) | 2013-05-23 |
CN103002794A (en) | 2013-03-27 |
CN103002794B (en) | 2016-08-17 |
WO2011097631A3 (en) | 2011-11-24 |
WO2011097631A2 (en) | 2011-08-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9013555B2 (en) | Method and apparatus for ultrahigh sensitive optical microangiography | |
Chen et al. | Optical coherence tomography based angiography | |
Grulkowski et al. | Scanning protocols dedicated to smart velocity ranging in spectral OCT | |
Zhang et al. | Wide-field imaging of retinal vasculature using optical coherence tomography-based microangiography provided by motion tracking | |
Braaf et al. | Angiography of the retina and the choroid with phase-resolved OCT using interval-optimized backstitched B-scans | |
Makita et al. | Comprehensive in vivo micro-vascular imaging of the human eye by dual-beam-scan Doppler optical coherence angiography | |
Wojtkowski | High-speed optical coherence tomography: basics and applications | |
An et al. | In vivo volumetric imaging of vascular perfusion within human retina and choroids with optical micro-angiography | |
JP6507615B2 (en) | Optical coherence tomography apparatus and program | |
Puyo et al. | In vivo laser Doppler holography of the human retina | |
Blatter et al. | Dove prism based rotating dual beam bidirectional Doppler OCT | |
Poddar et al. | Challenges and advantages in wide-field optical coherence tomography angiography imaging of the human retinal and choroidal vasculature at 1.7-MHz A-scan rate | |
US20150230708A1 (en) | Methods and systems for determining volumetric properties of a tissue | |
US20160066798A1 (en) | Methods and Systems for Determining Hemodynamic Properties of a Tissue | |
Ruminski et al. | OCT angiography by absolute intensity difference applied to normal and diseased human retinas | |
Poddar et al. | In vivo imaging of human vasculature in the chorioretinal complex using phase-variance contrast method with phase-stabilized 1-μ m swept-source optical coherence tomography | |
Mujat et al. | Visualizing the vasculature of the entire human eye posterior hemisphere without a contrast agent | |
Mujat et al. | Swept-source parallel OCT | |
Kumar et al. | Complementary capabilities of photoacoustic imaging to existing optical ocular imaging techniques | |
Haindl et al. | Absolute velocity profile measurement by 3-beam Doppler optical coherence tomography utilizing a MEMS scanning mirror | |
John et al. | Recent developments in optical coherence tomography angiography imaging for the diagnosis and assessment of diabetic retinopathy | |
Motaghiannezam et al. | In vivo human retinal and choroidal vasculature visualization using differential phase contrast swept source optical coherence tomography at 1060 nm | |
EP3858228A1 (en) | Ophthalmological imaging device, control method therefor, program, and recording medium | |
Li et al. | Optical Coherence Tomography Structural and Functional Imaging | |
Wang et al. | Optical microangiography: theory and application |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |