WO2023178078A1 - Oct speckle velocimetry - Google Patents
Oct speckle velocimetry Download PDFInfo
- Publication number
- WO2023178078A1 WO2023178078A1 PCT/US2023/064299 US2023064299W WO2023178078A1 WO 2023178078 A1 WO2023178078 A1 WO 2023178078A1 US 2023064299 W US2023064299 W US 2023064299W WO 2023178078 A1 WO2023178078 A1 WO 2023178078A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- flow
- oct
- information
- structural
- images
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Ceased
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/102—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for optical coherence tomography [OCT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B3/00—Apparatus for testing the eyes; Instruments for examining the eyes
- A61B3/10—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions
- A61B3/12—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes
- A61B3/1225—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation
- A61B3/1233—Objective types, i.e. instruments for examining the eyes independent of the patients' perceptions or reactions for looking at the eye fundus, e.g. ophthalmoscopes using coherent radiation for measuring blood flow, e.g. at the retina
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/26—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring the direct influence of the streaming fluid on the properties of a detecting optical wave
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N20/00—Machine learning
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10016—Video; Image sequence
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10101—Optical tomography; Optical coherence tomography [OCT]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30041—Eye; Retina; Ophthalmic
-
- G—PHYSICS
- G06—COMPUTING OR CALCULATING; COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Definitions
- OCT optical coherence tomography
- OCT is a non-invasive imaging technique, often used in ophthalmology.
- OCT relies on principles of interferometry to image and collect information about an object (such as the eye of a subject). Particularly, light from a source is split into a sample arm where it is reflected by the object being imaged, and reference arm where it is reflected by a reference object such as a mirror. The reflected lights are then combined in a detection arm in a manner that produces an interference pattern that is detected by spectrometer, photodiode(s) or the like. The detected interference signal is processed to reconstruct the object and generate structural OCT images.
- structural OCT images and volumes are generated by combining numerous depth profiles (A-lines, e.g. along a Z-depth direction at an X-Y location) into a single cross-sectional image (B-scan, e.g., as an X-Z or Y-Z plane), and combining numerous B-scans into a volume.
- depth profiles are generated by scanning along the X and Y directions.
- En- face images in the X-Y plane may be generated by flattening a volume in all or a portion of the Z- depth direction
- C-scan images may be generated by extracting slices of a volume at a given depth.
- OCT images are acquired in an en-face manner in the X-Y plane with the Z-depth acquired sequentially.
- Cross-sectional images in the X-Z or Y-Z planes may be generated from the acquired volume.
- OCT imaging may be used to determine blood flow properties, such as velocity.
- Doppler OCT measures the Doppler shifts caused when blood cells scatter the OCT light beam.
- Doppler OCT requires extracting phase information from raw spectral data and determining Doppler phase shifts. These processes can be complicated and difficult to efficiently perform.
- a method comprises: capturing optical coherence tomography (OCT) data from an object; generating a first plurality of structural OCT images from a first location of the object based on the captured OCT data; extracting flow information from individual ones of the first plurality of structural OCT images; and generating a first time-series flow profile of the first location of the object, wherein the first flow profile is a relationship between the extracted flow information and a timing of the captured OCT data from which the corresponding individual one of the first plurality of structural OCT images was generated.
- OCT optical coherence tomography
- extracting flow information comprises applying a high frequency filter to frequency information of the individual one of the first plurality of structural OCT images, thereby producing high frequency information, and the flow information corresponds to the high frequency information; extracting flow information further comprises: applying a low frequency filter to the frequency information of the individual one of the first plurality of structural OCT images, thereby producing low frequency information, wherein the flow information is a ratio of the high frequency information to the low frequency information; extracting flow information comprises: applying two-dimensional Fourier transform to the individual one of the first plurality of structural OCT images, thereby producing the frequency information; the extracted flow information is a speckle density of the individual one of the first plurality of structural OCT images; extracting flow information comprises: applying a cooccurrence matrix to the first plurality of structural OCT images, and determining a correlation among the first plurality of structural OCT images based on the co-occurrence matrix; extracting flow information comprises: inputting the individual one of the first plurality of structural OCT images to a machine learning system
- Figure 1 illustrates an example optical coherence tomography (OCT) system.
- Figure 2 illustrates an example method according to the present disclosure.
- Figure 3 illustrates an example method for extracting blood flow information based on frequency.
- Figure 4A illustrates a first example method for extracting blood flow information based on speckle density.
- Figure 4B illustrates a second example method for extracting blood flow information based on speckle density.
- Figure 5 illustrates an example method for extracting blood flow information based on a co-occurrence matrix.
- Figure 6A comparatively illustrates a flow profile determined according to the present disclosure with a corresponding structural OCT signal from which the flow profile was determined.
- Figure 6B comparatively illustrates a blood flow map and structural OCT image for the frames identified in Fig. 6A.
- Figure 7 illustrates an example method for determining flow velocity.
- Figure 8A illustrates an en-face image identifying scan locations for determining a flow velocity.
- Figure 8B illustrates flow profiles of the locations indicated in Fig. 8A, and the determination of flow velocity therefrom.
- Figure 9 illustrates an example method of glaucoma testing.
- the present disclosure relates to determining blood flow information from optical coherence tomography (OCT) images without the need for complex processing. More particularly, the present disclosure relates to systems and methods for determining blood flow information from structural OCT images without a Doppler or like phase- related analysis. Still more particularly, the present disclosure determines blood flow information based on speckle information in the structural OCT images.
- OCT optical coherence tomography
- FIG. 1 An example OCT system 100 such as that of the present disclosure is illustrated in Fig. 1.
- the system 100 includes a light source 101.
- the light generated by the light source 101 is split by, for example, a beam splitter (as part of interferometer optics 108), and sent to a reference arm 104 and a sample arm 106.
- the light in the sample arm 106 is backscattered or otherwise reflected off an object, such as the retina of an eye 112.
- the light in the reference arm 104 is backscattered or otherwise reflected, by a mirror 110 or like object.
- Light from the sample arm 106 and the reference arm 104 is recombined at the optics 108 and a corresponding interference signal is detected by a detector 102.
- the “speckle” of an OCT image is understood to mean the characteristic interference pattern produced by the backscattered OCT light from a turbid media such as tissue microstructure.
- the detector 102 can be a spectrometer, photo detector, or any other light detecting device.
- the detector 102 outputs an electrical signal corresponding to the interference signal to one or more processors 114.
- the processor(s) 114 may process the electrical signal into OCT signal data, generate corresponding structural images, and/or further analyze the data and images. Particularly, the processor(s) 114 may implement any or all aspects of the present disclosure.
- the processor(s) 114 may also be associated with an input/output interface (not shown) including a display for outputting processed images, or information related to the analysis of those images.
- the input/output interface may also include hardware such as buttons, keys, or other controls for receiving user inputs to the system.
- the processor 114 may also be used to control the light source and imaging process.
- FIG. 2 An example method of the present disclosure is illustrated generally in Fig. 2.
- the method begins by generating a time-series of OCT structural images at a location of the object being imaged.
- the OCT system 100 of Fig. 1 images the object according to a scan pattern that produces OCT data from the same location at a plurality of times.
- the OCT system 100 may repeatedly scan a single cross-sectional plane location of the object so as to generate a plurality of B-scans of the cross-sectional plane location at different times.
- the OCT system 100 may repeatedly scan an area of the object so as to generate a plurality of volumes of the object, from which B-scans at different times can be derived.
- the scanning protocol of the OCT system 100 may collect data for a long enough period of time to capture one or more cardiac pulses.
- the OCT data may be collected for a time on the order of a few seconds.
- OCT data is collected for between two and three seconds.
- the rate and density of the OCT data capture may be oversampled to improve the amount of speckle information in each resulting structural OCT image.
- the OCT system 100 may operate at a frequency great enough to obtain at least 50 repeated structural OCT B-scans per second at a common cross-sectional plane location.
- each structural OCT B-scan is about 1 mm wide with between 500 to 1024 A-lines of information.
- the object may be tracked to ensure movement is limited during the capture period. If movement is too great, thereby causing too much noise or an inability to register the resulting OCT B-scans or volumes, the scanning protocol may be reset.
- blood flow information is extracted from the individually generated structural OCT images.
- blood flow information may be extracted from an entire OCT B-scan, or only from a region of interest (ROT).
- ROI region of interest
- the ROT may be determined manually or automatically.
- the ROT may be manually determined by a clinician selecting a region of an OCT B-scan corresponding to an area of blood flow or vasculature. Measurements of the vasculature, such as its size, may also be measured manually by a clinician from the structural OCT images. Such a region and measurements may instead be automatically determined by a variety of methods.
- the region of vasculature may be identified by segmentation techniques, machine learning techniques, thresholding techniques (blood flow and vasculature having a greater intensity and variation), and the like. Further, additional ROIs may be identified for other vessels.
- each of the repeated B-scans at the common cross-sectional plane location may be registered to each other. Therefore, the ROT may only be identified in a B-scan image and then extrapolated to the other repeated B-scan images at that location. In the event that images are acquired in an en-face manner in the X-Y plane, those X-Y plane images may also be registered to each other. Thus, the ROT would only need to be identified in a single X-Y plane image.
- Blood flow information may be extracted from individual OCT B-scans in a variety of ways.
- the techniques for extracting blood flow information from structural OCT B- scans are based on the recognition that blood flow produces greater speckle variation than static structural tissue.
- FIG. 3 A first example method for extracting blood flow information based on frequency analysis is illustrated in Fig. 3.
- an ROT is identified within an area of an OCT structural image corresponding to blood flow and vasculature (as indicated by the higher intensity signal and higher density speckle pattern).
- the ROT is then transformed (e.g., according to a 2D Fourier transform) to produce 2D frequency information of the ROT, which may then be frequency filtered.
- high and low frequency filters are applied to frequency information of the transformed ROT.
- the cutoff frequency of the filters may be determined based on a size of the ROT (or 2D Fourier transform region) in relation to the size of blood vessel speckles under analysis.
- the resulting signal is then statistically combined into a single value to produce a high frequency signal FH, and a low frequency signal FL.
- the statistical combination is an average.
- the high frequency signal FH can be understood to represent blood flow
- the low frequency signal FL represents static tissue. Therefore, in some embodiments, the high frequency signal FH alone may be determined and further processed, and correspond to the flow information. However, in some embodiments, normalizing the high frequency signal FH (for example, to the low frequency signal FL) may help account for natural variations in the imaged object and the OCT system 100, thereby improving the quality of the extracted flow information. In other words, the flow information may be determined as the ratio between the high and low frequency signals FH, FL.
- flow information may be extracted based on a speckle density.
- a threshold may be applied to the pixels within the ROI.
- a binary map of the ROI may then be generated based on an applied threshold. For example, all pixels greater than the threshold may be set to a value of 1, while all pixels less than threshold may be set to a value of 0.
- a speckle density may then be determined as a ratio of the number of pixels having a value of 1 to the total number of pixels of the ROI, or to the number of pixels having a value of 0. The flow information may thus be considered equal to or a function of speckle density.
- FIG. 4B A second example technique for determining speckle density is illustrated in Fig. 4B.
- an edge fdter is first applied to a structural OCT B-scan, which identifies edges of the speckles.
- An averaging filter is then applied to the speckle edges (the output of the edge filter) and to the structural OCT B-scan. Applying the averaging filter produces an edge strength and a mean intensity of the structural OCT B-scan, respectively.
- the averaging filter preferably has size that is close to the size of the blood vessel in the structural OCT B-scan from which the blood flow information is being determined.
- the speckle density is proportional to the density of visible edges
- the speckle density can be determined as a ratio between the edge strength and the mean intensity (e.g., edge strength divided by mean intensity).
- the flow information may be considered the speckle density (e.g., the ratio) or a function of the speckle density.
- speckle density could be determined by applying a threshold intensity to the speckle edges.
- the speckle density may be considered the number of speckle edges exceeding a threshold intensity.
- normalizing e.g., to the mean intensity as discussed above
- flow information may be determined by applying a cooccurrence matrix to the ROI.
- the co-occurrence matrix creates a dependence matrix by determining how often a pixel with pixel intensity value i occurs in the adjacent frame to a pixel with the value j.
- Each element (i, j) in the dependence matrix specifies the number of times that the pixel with a value i occurred in the adj acent frame to a pixel with value j.
- the flow information is then determined by applying the correlation of a pixel to its neighbor frames on the cooccurrence matrix.
- An example application of such a co-occurrence matrix application is illustrated in Fig. 5.
- the size of a moving window applied on temporally adjacent structural OCT B-scan images is first determined. According to one example, the moving window size is four.
- a co-occurrence matrix is then applied to the structural OCT B-scan image data inside this analysis window.
- the flow information is then determined by applying a correlation property on the cooccurrence matrix.
- flow may be determined by a machine learning system.
- a machine learning system may be trained to output one or more values representing flow information based on an input structural ROI, OCT B-scan, or OCT volume.
- the machine learning system may be trained in a supervised manner based on input structural OCT image information and a corresponding ground truth flow information.
- the ground truth flow information may be determined, for example, according to one of the above-described techniques.
- the ground truth flow information may additionally or alternatively be determined by other techniques, such as Doppler OCT or even non-OCT analysis techniques. Accordingly, the machine learning system is trained to recognize a relationship between a structural speckle signal and corresponding flow information.
- flow may be determined by applying known ultrasound and laser speckle techniques, such as speckle decorrelation, speckle contrast, and speckle auto-correlation. Each of these techniques could be applied on individual structural OCT B-scans and/or on temporally adjacent structural OCT B-scan images.
- flow may be determined by applying known amplitude-based OCT angiography techniques such as speckle variance, amplitude-decorrelation, or split-spectrum amplitude-decorrelation. As above, each of these techniques could be applied on individual structural OCT B-scans and/or on temporally adjacent structural OCT B-scan images.
- a flow profile may be determined as a time-series plot of the flow information.
- the ratio between the high and low frequency signals FH, FL, the speckle density, or like flow-related information values are plotted according to a time series in which the relative time for each flow-related information value corresponds to a time at which the corresponding repeated OCT image was captured.
- FIG. 6A An example flow profile is illustrated in Fig. 6A along with a corresponding structural OCT signal.
- the structural OCT signal represents an average pixel intensity within the ROI form which the flow information of the flow profile was extracted.
- the flow profile is similar to a traditional pulse waveform. Indeed, the dicrotic notch is even visible roughly between image frames 15 and 20.
- the structural OCT signal is relatively flat and bears no relationship to a pulse wave.
- Fig. 6A further identifies four image frames (9, 12, 29, and 36) which are illustrated in Fig. 6B. More particularly, Fig. 6B illustrates the original structural OCT image for each frame along with a blood flow map for each frame.
- the blood flow map is a map of extracted flow information for the entire area of the structural OCT image, while the time-series signals shown in Fig. 6A only represent the ROI noted in the image frames in Fig. 6B.
- a flow profile may be generated by a machine learning system.
- a machine learning system may be trained to output a flow profile based on an input time-series of structural OCT images.
- the first step of an example method for determining flow velocity involves simultaneously capturing OCT data from at least two different locations having a known distance therebetween.
- the OCT data capture occurs alternately between the different locations such that the resulting B-scan images from each location are temporally interleaved.
- Fig. 8A illustrates an en-face image in which two locations A, B are identified along a blood vessel therein. The locations A, B are separated by a distance d.
- OCT data may be captured according to circular scans having different radii. The difference in radius between each circular scan thus represents a known distance. These circular scans could be temporally interleaved in the same manner discussed above.
- flow profiles for the data from each location are generated. These flow profiles may be generated from any of the above-discussed methods.
- Fig. 8B illustrates flow profiles corresponding to the locations A, B of Fig. 8A.
- a time At between these two flow profiles can then be determined.
- the time t is determined between local maxima of the flow profiles; however, the time may instead be determined between other common locations of the flow profiles (e g., local minima).
- the local maxima may be determined, for example, by identifying the frame at which a maximum value in each cardiac cycle occurs.
- a derivative of the flow profile may be determined and then the frames at which the derivative equals (or crosses) zero correspond to local maxima and minima.
- the time At may then be determined based on the sampling rate at which OCT data was captured. In other words, if OCT B-scans are captured at 50 frames per second, a difference of 50 frames between local maxima of the different flow profiles represents a 1 second difference.
- a direction of the flow may be further determined according to whether a phase shift of the OCT signal between the two locations is positive for negative. This velocity may be determined multiple times between the same two locations and/or between multiple locations.
- the plurality of determined velocities may then be averaged (or combined according to another statistical determination) in order to identify a representative velocity.
- determined velocities can be compared between a plurality of patients to identify abnormalities, or between the same patient at multiple capture times to determine a change in the patient’s condition.
- any part of the above description may be incorporated with other analysis techniques and/or processing.
- the flow profdes may be averaged over several cardiac cycles or compared over a period of time (e.g., weeks, months, years) to monitor disease progression.
- a period of time e.g., weeks, months, years
- the determined distance of blood flow may be used to distinguish arteries from veins in the structural images.
- veins and arteries could further be distinguished by comparison of determined velocity.
- the above OCT data acquisition could be performed concurrently with other physiological measurements, such as pulse oximetry, electrocardiography, and the like. Further analysis of the flow profiles herein may be based on additional physiological information captured by the concurrent measurements.
- any vascular biomarkers of glaucomatous damage such as impaired vascular flow and autoregulation may be used to provide diagnostic information. Accordingly, it may be possible to detect glaucoma, or quantify its severity, by measuring changes in blood flow in response to an external stimulus. An example of such a method is illustrated in Fig. 9.
- an external stimulus is applied.
- a stimulus may be an applied pressure (e.g., an air puff) or other stimulus that produces a change in intraocular pressure.
- a flow profile may then be determined for the patient according to any of the above techniques. The resulting flow profile may be compared to known healthy patients or a previously determined baseline flow profile for that patient. Changes in the flow profile relative to the baseline may be indicative of a healthy eye, as flow in the eye (and thus the vessels of the eye) are affected by the applied pressure. In contrast, no (or relatively minimal) change in flow profile relative to the baseline may be indicative of glaucoma since the blood flow and vessel are not responsive to the applied pressure.
- the above analysis may also be applied vice versa depending on the applied external stimulus.
- a flow profile may be determined according to the present disclosure before, during, or after the application of any external stimulus. Response of the flow profile to the stimulus may then be analyzed for diagnostic or like purposes.
- the extracted flow information may be further analyzed and/or displayed.
- an analysis can include any known techniques for analyzing pulse waveforms Such techniques may include those for determining heart rate, rise time, flow skew, blood pressure, cardiac function, vessel stiffness, and the like.
- display of the extracted flow information may include the flow profile itself, the blood flow maps, and/or the structural OCT images.
- the extracted flow information may be mapped to pixel intensity and/or color in the blood flow maps.
- the blood flow maps may be shown as videos rather than static images.
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Ophthalmology & Optometry (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Physics & Mathematics (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Multimedia (AREA)
- Aviation & Aerospace Engineering (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Hematology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Eye Examination Apparatus (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Analysis (AREA)
Abstract
Description
Claims
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202380027971.2A CN119013693A (en) | 2022-03-14 | 2023-03-14 | OCT speckle velocimetry |
| EP23771572.7A EP4494094A4 (en) | 2022-03-14 | 2023-03-14 | OCT-SPECKLE VELOCITY MEASUREMENT |
| JP2024554731A JP2025509539A (en) | 2022-03-14 | 2023-03-14 | OCT speckle flow measurement |
| US18/846,538 US20250185908A1 (en) | 2022-03-14 | 2023-03-14 | Oct speckle velocimetry |
| DE23771572.7T DE23771572T1 (en) | 2022-03-14 | 2023-03-14 | OCT SPECKLE VELOCITY MEASUREMENT |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US202263269298P | 2022-03-14 | 2022-03-14 | |
| US63/269,298 | 2022-03-14 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| WO2023178078A1 true WO2023178078A1 (en) | 2023-09-21 |
Family
ID=88024366
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| PCT/US2023/064299 Ceased WO2023178078A1 (en) | 2022-03-14 | 2023-03-14 | Oct speckle velocimetry |
Country Status (6)
| Country | Link |
|---|---|
| US (1) | US20250185908A1 (en) |
| EP (1) | EP4494094A4 (en) |
| JP (1) | JP2025509539A (en) |
| CN (1) | CN119013693A (en) |
| DE (1) | DE23771572T1 (en) |
| WO (1) | WO2023178078A1 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117752298A (en) * | 2023-12-29 | 2024-03-26 | 山东大学 | Capillary blood flow velocity measurement method and system based on OCTA |
| WO2025231347A1 (en) * | 2024-05-03 | 2025-11-06 | Topcon Corporation | Optical coherence tomography speckle velocimetry |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US12543155B2 (en) * | 2022-05-31 | 2026-02-03 | Mediatek Inc. | Cooperation with digital frequency-translation repeater for uplink transmission and reception-UE behavior |
| US12542695B2 (en) * | 2022-05-31 | 2026-02-03 | Mediatek Inc. | Cooperation with digital frequency-translation repeater for uplink transmission and reception-base station behavior |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6173197B1 (en) * | 1996-11-09 | 2001-01-09 | Moor Instruments Limited | Apparatus for measuring microvascular blood flow |
| US20050254008A1 (en) * | 2002-06-14 | 2005-11-17 | Ferguson R D | Monitoring blood flow in the retina using a line-scanning laser ophthalmoscope |
| US9092691B1 (en) * | 2014-07-18 | 2015-07-28 | Median Technologies | System for computing quantitative biomarkers of texture features in tomographic images |
| US20170238798A1 (en) * | 2016-02-03 | 2017-08-24 | Nidek Co., Ltd. | Optical coherence tomography device |
| US20180315194A1 (en) * | 2017-04-28 | 2018-11-01 | Massachusetts Institute Of Technology | Systems and methods for generating and displaying oct angiography data using variable interscan time analysis |
| US20180317851A1 (en) * | 2017-05-08 | 2018-11-08 | Oregon Health & Science University | Bulk motion subtraction in optical coherence tomography angiography |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150371401A1 (en) * | 2013-03-13 | 2015-12-24 | University Of Washington | Methods and Systems for Imaging Tissue Motion Using Optical Coherence Tomography |
| JP6520099B2 (en) * | 2014-06-30 | 2019-05-29 | 株式会社ニデック | Optical coherence tomography apparatus and data processing program |
| EP3103386B1 (en) * | 2015-06-08 | 2019-04-10 | Tomey Corporation | Blood flow velocity measuring device, program, and method |
| DE102017001677B3 (en) * | 2016-12-30 | 2018-03-29 | I-Dynamometer GmbH | Method for determining the occurrence of vascular collapse in a blood vessel in or on the eye |
| JP7479030B2 (en) * | 2020-06-11 | 2024-05-08 | 国立大学法人旭川医科大学 | BLOOD FLOW ANALYSIS APPARATUS, BLOOD FLOW ANALYSIS METHOD, AND PROGRAM |
-
2023
- 2023-03-14 EP EP23771572.7A patent/EP4494094A4/en active Pending
- 2023-03-14 JP JP2024554731A patent/JP2025509539A/en active Pending
- 2023-03-14 CN CN202380027971.2A patent/CN119013693A/en active Pending
- 2023-03-14 DE DE23771572.7T patent/DE23771572T1/en active Pending
- 2023-03-14 WO PCT/US2023/064299 patent/WO2023178078A1/en not_active Ceased
- 2023-03-14 US US18/846,538 patent/US20250185908A1/en active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6173197B1 (en) * | 1996-11-09 | 2001-01-09 | Moor Instruments Limited | Apparatus for measuring microvascular blood flow |
| US20050254008A1 (en) * | 2002-06-14 | 2005-11-17 | Ferguson R D | Monitoring blood flow in the retina using a line-scanning laser ophthalmoscope |
| US9092691B1 (en) * | 2014-07-18 | 2015-07-28 | Median Technologies | System for computing quantitative biomarkers of texture features in tomographic images |
| US20170238798A1 (en) * | 2016-02-03 | 2017-08-24 | Nidek Co., Ltd. | Optical coherence tomography device |
| US20180315194A1 (en) * | 2017-04-28 | 2018-11-01 | Massachusetts Institute Of Technology | Systems and methods for generating and displaying oct angiography data using variable interscan time analysis |
| US20180317851A1 (en) * | 2017-05-08 | 2018-11-08 | Oregon Health & Science University | Bulk motion subtraction in optical coherence tomography angiography |
Non-Patent Citations (2)
| Title |
|---|
| BRAAF BOY, DONNER SABINE, URIBE-PATARROYO NÉSTOR, BOUMA BRETT E., VAKOC BENJAMIN J.: "A Neural Network Approach to Quantify Blood Flow from Retinal OCT Intensity Time-Series Measurements", SCIENTIFIC REPORTS, vol. 10, no. 1, XP093093326, DOI: 10.1038/s41598-020-66158-8 * |
| See also references of EP4494094A4 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117752298A (en) * | 2023-12-29 | 2024-03-26 | 山东大学 | Capillary blood flow velocity measurement method and system based on OCTA |
| CN117752298B (en) * | 2023-12-29 | 2024-05-24 | 山东大学 | Capillary blood flow velocity measurement method and system based on OCTA |
| WO2025231347A1 (en) * | 2024-05-03 | 2025-11-06 | Topcon Corporation | Optical coherence tomography speckle velocimetry |
Also Published As
| Publication number | Publication date |
|---|---|
| US20250185908A1 (en) | 2025-06-12 |
| EP4494094A4 (en) | 2026-03-18 |
| JP2025509539A (en) | 2025-04-11 |
| CN119013693A (en) | 2024-11-22 |
| DE23771572T1 (en) | 2025-03-27 |
| EP4494094A1 (en) | 2025-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20250185908A1 (en) | Oct speckle velocimetry | |
| CN104411233B (en) | In vivo optical flow imaging | |
| CN106943124B (en) | Quantification of local circulation with optical coherence tomography angiography | |
| CN107862724B (en) | An Improved Microvascular Blood Flow Imaging Method | |
| US20110141259A1 (en) | Image processing apparatus, method for image processing,image pickup system, and computer-readable storage medium | |
| CN106166058B (en) | One kind being applied to optical coherence tomography blood vessel imaging method and OCT system | |
| CN107452029A (en) | A kind of optics microvascular blood flow imaging method | |
| Belghith et al. | A hierarchical framework for estimating neuroretinal rim area using 3D spectral domain optical coherence tomography (SD-OCT) optic nerve head (ONH) images of healthy and glaucoma eyes | |
| JP7607820B2 (en) | Ophthalmology information processing device | |
| JP7191166B2 (en) | Ophthalmic photographing device and ophthalmic information processing device | |
| US12150711B2 (en) | Blood flow analysis apparatus, blood flow analysis method, and recording medium | |
| Greenspan et al. | Doppler echocardiography flow-velocity image analysis for patients with atrial fibrillation | |
| CN115844353B (en) | Imaging system and method for contactless detection of pulse wave spatial and temporal distribution and characteristics | |
| WO2025231347A1 (en) | Optical coherence tomography speckle velocimetry | |
| JP2019042263A (en) | Ophthalmologic system, ophthalmologic imaging apparatus, and ophthalmologic information processing apparatus | |
| JP2019042264A (en) | Ophthalmologic imaging apparatus and ophthalmologic information processing apparatus | |
| JP2021037088A (en) | Optical interference tomographic imaging device, operation method of optical interference tomographic imaging device, and program | |
| US20250194920A1 (en) | Vista de-noising | |
| JP7624619B2 (en) | BLOOD FLOW ANALYSIS APPARATUS, OPHTHALMOLOGIC APPARATUS, OPERATION METHOD OF BLOOD FLOW ANALYSIS APPARATUS, AND PROGRAM | |
| WO2025183024A1 (en) | Ophthalmic device, method for controlling same, program, and recording medium | |
| Mujat et al. | Retinal nerve fiber layer thickness map and blood flow pulsation measured with SDOCT |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| 121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 23771572 Country of ref document: EP Kind code of ref document: A1 |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 18846538 Country of ref document: US |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 2024554731 Country of ref document: JP |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 202380027971.2 Country of ref document: CN |
|
| WWE | Wipo information: entry into national phase |
Ref document number: 2023771572 Country of ref document: EP |
|
| NENP | Non-entry into the national phase |
Ref country code: DE |
|
| ENP | Entry into the national phase |
Ref document number: 2023771572 Country of ref document: EP Effective date: 20241014 |
|
| WWP | Wipo information: published in national office |
Ref document number: 18846538 Country of ref document: US |