US20210000352A1 - Fluorescence based flow imaging and measurements - Google Patents
Fluorescence based flow imaging and measurements Download PDFInfo
- Publication number
- US20210000352A1 US20210000352A1 US17/031,375 US202017031375A US2021000352A1 US 20210000352 A1 US20210000352 A1 US 20210000352A1 US 202017031375 A US202017031375 A US 202017031375A US 2021000352 A1 US2021000352 A1 US 2021000352A1
- Authority
- US
- United States
- Prior art keywords
- pixels
- frame
- time
- light
- calculate
- 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
- 238000003384 imaging method Methods 0.000 title description 8
- 238000005259 measurement Methods 0.000 title description 6
- 239000003550 marker Substances 0.000 claims abstract description 55
- 239000013598 vector Substances 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 24
- 210000001124 body fluid Anatomy 0.000 claims abstract description 20
- 239000012530 fluid Substances 0.000 claims abstract description 18
- 238000012544 monitoring process Methods 0.000 claims abstract description 6
- 238000012384 transportation and delivery Methods 0.000 claims description 29
- 238000009499 grossing Methods 0.000 claims description 16
- 238000012935 Averaging Methods 0.000 claims description 6
- 230000000630 rising effect Effects 0.000 claims description 4
- 239000000975 dye Substances 0.000 description 8
- 230000004044 response Effects 0.000 description 7
- 230000017531 blood circulation Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000005286 illumination Methods 0.000 description 5
- 239000008280 blood Substances 0.000 description 3
- 210000004369 blood Anatomy 0.000 description 3
- 230000000875 corresponding effect Effects 0.000 description 3
- MOFVSTNWEDAEEK-UHFFFAOYSA-M indocyanine green Chemical compound [Na+].[O-]S(=O)(=O)CCCCN1C2=CC=C3C=CC=CC3=C2C(C)(C)C1=CC=CC=CC=CC1=[N+](CCCCS([O-])(=O)=O)C2=CC=C(C=CC=C3)C3=C2C1(C)C MOFVSTNWEDAEEK-UHFFFAOYSA-M 0.000 description 3
- 229960004657 indocyanine green Drugs 0.000 description 3
- 230000004298 light response Effects 0.000 description 3
- 210000000056 organ Anatomy 0.000 description 3
- 210000004204 blood vessel Anatomy 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- GNBHRKFJIUUOQI-UHFFFAOYSA-N fluorescein Chemical compound O1C(=O)C2=CC=CC=C2C21C1=CC=C(O)C=C1OC1=CC(O)=CC=C21 GNBHRKFJIUUOQI-UHFFFAOYSA-N 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000002347 injection Methods 0.000 description 2
- 239000007924 injection Substances 0.000 description 2
- 238000001990 intravenous administration Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 238000005316 response function Methods 0.000 description 2
- 210000001321 subclavian vein Anatomy 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 238000001994 activation Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000004347 surface barrier Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- 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/0275—Measuring blood flow using tracers, e.g. dye dilution
-
- 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/0071—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by measuring fluorescence emission
-
- 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7246—Details of waveform analysis using correlation, e.g. template matching or determination of similarity
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61M—DEVICES FOR INTRODUCING MEDIA INTO, OR ONTO, THE BODY; DEVICES FOR TRANSDUCING BODY MEDIA OR FOR TAKING MEDIA FROM THE BODY; DEVICES FOR PRODUCING OR ENDING SLEEP OR STUPOR
- A61M5/00—Devices for bringing media into the body in a subcutaneous, intra-vascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, arm-rests
- A61M5/007—Devices for bringing media into the body in a subcutaneous, intra-vascular or intramuscular way; Accessories therefor, e.g. filling or cleaning devices, arm-rests for contrast media
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
- G01F1/708—Measuring the time taken to traverse a fixed distance
- G01F1/7086—Measuring the time taken to traverse a fixed distance using optical detecting arrangements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/001—Image restoration
- G06T5/002—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration by the use of local operators
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
- G06T7/248—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments involving reference images or patches
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/254—Analysis of motion involving subtraction of images
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/39—Markers, e.g. radio-opaque or breast lesions markers
- A61B2090/3933—Liquid markers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B90/00—Instruments, implements or accessories specially adapted for surgery or diagnosis and not covered by any of the groups A61B1/00 - A61B50/00, e.g. for luxation treatment or for protecting wound edges
- A61B90/39—Markers, e.g. radio-opaque or breast lesions markers
- A61B2090/3937—Visible markers
- A61B2090/3941—Photoluminescent markers
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
- G01F1/708—Measuring the time taken to traverse a fixed distance
- G01F1/712—Measuring the time taken to traverse a fixed distance using auto-correlation or cross-correlation detection means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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; CALCULATING OR 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/10024—Color image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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/10064—Fluorescence image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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/20212—Image combination
- G06T2207/20224—Image subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR 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/30204—Marker
Definitions
- the present disclosure relates to fluorescence based imaging and measurements. More particularly, the present disclosure relates to determining flow characteristics such as velocity in bodily vessels such as blood vessels.
- Fluorescent markers have been used for basic imaging of bodily structures, but improvements are needed in determining flow characteristics in such bodily fluids as blood.
- a method of fluorescence based tracking of a light-emitting marker in a bodily fluid stream includes: providing a light-emitting marker into a bodily fluid stream; monitoring, with a sensor, a region of interest traversed by the bodily fluid stream; recording data generated by the sensor; determining time characteristics of the recorded data; and calculating flow characteristics based on the time characteristics.
- the senor includes a camera
- the recorded data comprises motion video data
- the method further includes: dividing the motion video data into kernels; identifying which of the kernels receive some portion of the light-emitting marker using an intensity threshold; computing, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; performing smoothing on each intensity signal data set; calculating a lag time between the intensity signal data sets of neighboring identified kernels using cross-correlation; using a spatial resolution and the lag time, calculating velocity vectors; summing the velocity vectors of neighboring kernels to create a resultant velocity vector; and generating a velocity map from the resultant velocity vectors for all kernels.
- performing smoothing on each intensity signal data set includes time window averaging.
- performing smoothing on each intensity signal data set includes using a filter.
- performing smoothing on each intensity signal data set includes using a Gaussian filter.
- the method further includes: dividing the motion video data into kernels; identifying which of the kernels receive some portion of the light-emitting marker using an intensity threshold; computing, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; performing smoothing on each intensity signal data set; for each particular identified kernel, finding segments in which a slope of the intensity signal data set rises for a minimum consecutive number of frames or falls for a minimum consecutive number of frames, which segments occur when a leading edge or falling edge of a portion of the light-emitting marker passes through the identified kernel; searching the intensity signal data sets of neighboring identified kernels for a rising or falling segment of similar length; calculating a lag time between segments in the particular identified kernel and segments in the neighboring identified kernels; using a spatial resolution and the lag time, calculating velocity vectors; summing the velocity vectors of neighboring kernels to create a resultant velocity vector; and generating a velocity map from the resultant velocity vectors for all kernels.
- performing smoothing on each intensity signal data set includes time window averaging. In at least one example, performing smoothing on each intensity signal data set includes using a filter. In at least one example, performing smoothing on each intensity signal data set includes using a Gaussian filter.
- the method further includes: calculating a difference frame by subtracting a frame of the motion video data from a consecutive frame of the motion video data; applying a threshold the difference frame to eliminate pixels therein below a specified intensity value; calculating a pixel size of a remaining blob in the difference frame in a direction of blood flow; calculating a size of the remaining blob using the pixel size and a spatial resolution; and calculating a velocity by using a distance traveled by the remaining and a time between frames.
- the method further includes: dividing the motion video data into frames each including pixels; identifying which of the pixels receive some portion of the light-emitting marker using an intensity threshold; creating a logical frame in which a respective indicator for each pixel can be set as true or false; setting the indicators of the identified pixels as true; setting the indicators of all other pixels as false; calculating a difference frame by subtracting a first logical frame from a second logical frame such that the difference frame includes pixels that reached the specified threshold after a time of the first logical frame; finding length in pixels of the remaining blob in the difference frame in a direction of blood flow; converting the length in pixels of the difference frame to physical distance using the spatial resolution; and calculating velocity by dividing the physical distance by a time between frames.
- a system for fluorescence based tracking of a light-emitting marker in a bodily fluid stream includes: a delivery apparatus configured to provide a light-emitting marker into a bodily fluid stream; a sensor configured to monitor a region of interest traversed by the bodily fluid stream; and a computing device configured to: record data generated by the sensor; determine time characteristics of the recorded data; and calculate flow characteristics based on the time characteristics.
- the senor includes a camera
- the recorded data includes motion video data
- the computing device is further configured to: divide the motion video data into kernels; identify which of the kernels receive some portion of the light-emitting marker using an intensity threshold; compute, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; perform smoothing on each intensity signal data set; calculate a lag time between the intensity signal data sets of neighboring identified kernels using cross-correlation; using a spatial resolution and the lag time, calculate velocity vectors; sum the velocity vectors of neighboring kernels to create a resultant velocity vector; and generate a velocity map from the resultant velocity vectors for all kernels.
- the computing device performs smoothing on each intensity signal data set by time window averaging. In at least one example, the computing device performs smoothing on each intensity signal data set by using a Gaussian filter.
- the computing device is further configured to: divide the motion video data into kernels; identify which of the kernels receive some portion of the light-emitting marker using an intensity threshold; compute, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; perform smoothing on each intensity signal data set; for each particular identified kernel, find segments in which a slope of the intensity signal data set rises for a minimum consecutive number of frames or falls for a minimum consecutive number of frames, which segments occur when a leading edge or falling edge of a portion of the light-emitting marker passes through the identified kernel; search the intensity signal data sets of neighboring identified kernels for a rising or falling segment of similar length; calculate a lag time between segments in the particular identified kernel and segments in the neighboring identified kernels; use a spatial resolution and the lag time to calculate velocity vectors; sum the velocity vectors of neighboring kernels to create a resultant velocity vector; and generate a velocity map from the resultant velocity vectors for all kernels.
- the computing device is further configured to: calculate a difference frame by subtracting a frame of the motion video data from a consecutive frame of the motion video data; apply a threshold the difference frame to eliminate pixels therein below a specified intensity value; calculate a pixel size of a remaining blob in the difference frame in a direction of blood flow; calculate a size of the remaining blob using the pixel size and a spatial resolution; and calculate a velocity by using a distance traveled by the remaining blob and a time between frames.
- the computing device is further configured to: divide the motion video data into frames each including pixels; identify which of the pixels receive some portion of the light-emitting marker using an intensity threshold; create a logical frame in which a respective indicator for each pixel can be set as true or false; set the indicators of the identified pixels as true; set the indicators of all other pixels as false; calculate a difference frame by subtracting a first logical frame from a second logical frame such that the difference frame includes pixels that reached the specified threshold after a time of the first logical frame; find length in pixels of the remaining blob in the difference frame in a direction of blood flow; convert the length in pixels of the difference frame to physical distance using the spatial resolution; and calculate velocity by dividing the physical distance by a time between frames.
- FIG. 1 shows a fluorescent marker delivery time plot in a fluid stream and a corresponding response time plot of light intensity measured downstream in a fixed field of view according to at least one embodiment.
- FIG. 2 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using a cross correlation technique.
- FIG. 3 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using leading and falling edge considerations.
- FIG. 4 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using subtraction.
- FIG. 5 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using thresholding.
- FIG. 6 shows a system, according to at least one embodiment, by which at least the methods of FIGS. 2-5 are implemented.
- Fluorescence based tracking includes the providing of a marker such as a glowing dye into a fluid stream, such as a bloodstream, and making measurements and generating imagery based on the arrival, movement, and departure of the marker downstream as detected by sensor(s) to characterize the flow of the fluid stream and vessels or structures within which the flow travels.
- the marker is provided into a fluid stream for example by direct injection or via a port as discrete bolus deliveries separated over time.
- a bolus refers to the administration of a discrete amount of a fluid substance, in this case the marker into a bodily fluid stream such as blood, in order to provide a concentration of the substance to gain a response.
- a bolus can be delivered by active pumping or by passive gravity based delivery such as via an intravenous drip line.
- a central line delivery arrangement is used, in which a port is placed in fluid communication with the subclavian vein and bolus deliveries are injected into the port.
- the dye briefly fluoresces when excited by an illumination source that emits a particular range of wavelengths.
- the dye is illuminated over the Region of Interest (ROI) where imaging of the fluorescence is also performed.
- ROI Region of Interest
- a field of view monitoring is established by placement of a sensor, such as a high speed camera, at a region of interest.
- the field of view can be established and held generally fixed as the marker enters and traverses the field of view of a high-speed camera sensitive to the light emitted by the marker. Time characteristics of the light output of the marker traversing the field of view can be deducted from the light output intensity as recorded by the camera.
- a field of view may be established for example at the heart or other organ where flow diagnostics are wanted.
- the visual response in the field of view indicates presence of the marker, with the intensity of the light response being correlated with the time evolving concentration of the marker in the stream as the marker diffuses and travels with the host fluid.
- the light intensity in the field of view may typically have both rise and fall characteristics.
- the rise characteristics correspond to the arrival and increasing concentration of the marker in the field of view.
- the fall characteristics correspond to the departure or diffusion of the marker and/or the reduction of its light output.
- rise time may be faster generally than fall time such that response time curves typically show steeper climbs than falls.
- FIG. 1 shows a fluorescent marker delivery time plot 102 in a fluid stream and a corresponding response time plot 112 of light intensity measured downstream in a fixed field of view according to at least one embodiment.
- the plots 102 and 112 are shown on a common time axis 100 .
- the delivery time plot 102 records several bolus deliveries 104 as step functions separated in time.
- the response plot 112 records the corresponding response functions 114 of light intensity in a field of view downstream from the marker delivery location into a fluid stream.
- Each response function 114 is shown as a wave or envelope having a rise side 116 representing the arrival of the marker in the field of view, and a fall side 118 indicating the departure or diffusion of the marker and/or the reduction of its light output.
- a marker delivery system in at least one embodiment includes a controller and a delivery pump in wired or wireless communication.
- a data acquisition system including the camera and a computing device for recording, analyzing and visualizing the data and/or field of view are also in communication with the controller in at least one embodiment.
- the time of initiation, delivered volume, and duration of each bolus delivery can be controlled.
- the time interval between consecutive bolus deliveries is also controlled.
- multiple parameters for bolus delivery can be adjusted to ultimately query and determine varied flow characteristics within a region of interest subjected to field of view monitoring. Shorter time gaps can be used for slower moving fluids and longer time gaps can be used for faster moving fluids within the region of interest. Injection times can be varied to address different anatomical depths and tissue surface barriers.
- the marker from a bolus delivery enters an established field of view
- the light response of a bolus is captured by a high speed camera.
- the time domain dynamics of the light response is analyzed to arrive at velocity vectors representing movement of the host fluid in the field of view.
- FIGS. 2-5 Several embodiments of generating velocity vectors using the data from fluorescence based tracking are described in the following with reference to FIGS. 2-5 .
- a video including multiple frames (images) taken of a region of interest is analyzed to track pixels in or across the field of view.
- a method 200 of generating a velocity vector map using a cross correlation technique is represented as a flow chart in FIG. 2 .
- step 202 divide the video into spatiotemporal cubes, which are termed “kernels” in these descriptions.
- step 204 which is optional, isolate kernels in which fluorescence appears at any point during the length of the video using an intensity threshold.
- step 206 compute a signal for each kernel of mean intensity vs. time.
- the signal may be 1D, 2D, or 3D.
- Vessel(s) may cross under other vessel(s) and at different angles, and multiple perspectives from different cameras can be used. The potential use of coded apertures is within the scope of these descriptions.
- step 208 which is optional, perform smoothing on the intensity signal for each kernel using time window averaging, Gaussian filter, etc.
- step 210 for each kernel, calculate the lag time between the intensity signal of the kernel and the time intensity signals of its neighboring kernels using cross-correlation.
- step 212 using the known spatial resolution of the image, convert lag time in each direction into velocity.
- step 214 sum the resulting velocity vectors to each neighboring kernel to create one resultant velocity vector.
- step 216 generate a velocity map from the resultant velocity vectors for all kernels.
- a method 300 of generating a velocity vector map using leading and falling edge considerations is represented as a flow chart in FIG. 3 .
- step 302 divide video into kernels.
- step 304 which is optional, isolate kernels in which fluorescence appears at any point during the length of the video using an intensity threshold.
- step 306 compute a signal for each kernel of mean intensity vs. time. As described above with reference to step 206 of FIG. 2 , the signal may be 1D, 2D, or 3D.
- step 308 which is optional, perform smoothing on the intensity signal for each kernel using time window averaging, Gaussian filter, etc.
- step 310 for each kernel, find segments in which the slope of the intensity signal rises for a minimum consecutive number of frames or falls for a minimum consecutive number of frames. These segments occur when the leading or falling edges of the ICG (fluorescein or other glow dye) bolus pass through the kernel.
- step 312 search the intensity signals of neighboring kernels for a rising or falling segment of similar length.
- step 314 calculate the temporal offset (lag time) between the segment in the original kernel and the segment in the neighboring kernels.
- step 316 using the known spatial resolution of the image, convert lag time in each direction into velocity.
- step 318 sum the resulting velocity vectors to each neighboring kernel to create one resultant velocity vector.
- step 320 generate a velocity map from the resultant velocity vectors for all kernels.
- a method 400 of generating a velocity vector map using subtraction is represented as a flow chart in FIG. 4 .
- step 402 subtract two consecutive frames, resulting in a difference frame.
- step 404 threshold difference frame to eliminate pixels below a specified intensity value.
- step 406 calculate pixel size of a remaining blob in the difference frame in the direction of blood flow.
- step 408 convert size of blob in pixels to physical distance using spatial resolution.
- step 410 calculate velocity by dividing distance traveled by time between frames.
- a method 500 of generating a velocity vector map using thresholding is represented as a flow chart in FIG. 5 .
- step 502 for each frame in video, isolate pixels with intensities above specified threshold. Create logical frame and set all pixels at or above threshold to true. Set all other pixels to false.
- step 504 for two consecutive frames, subtract the first logical frame from the second logical frame. The resulting frame contains the pixels that reached the specified threshold between frames.
- step 506 find the length in pixels of the remaining blob in the difference frame in the direction of blood flow.
- step 508 convert the pixel length of the difference frame to physical distance using the spatial resolution.
- step 510 calculate velocity by dividing physical distance by the time between frames.
- FIG. 6 shows a system 600 , according to at least one embodiment, by which at least the methods of FIGS. 2-5 are implemented.
- the system 600 includes a computing device 602 , a delivery apparatus 604 configured to provide a light-emitting marker 606 into a bodily fluid stream, and a sensor 610 configured to monitor a region of interest traversed by the bodily fluid stream.
- the computing device 602 records data generated by the sensor 610 , determines time characteristics of the recorded data; and calculates flow characteristics based on the time characteristics.
- the computing device 602 is illustrated as a laptop or other personal computer. Other computing devices including local and remote servers are within the scope of these descriptions.
- the delivery apparatus 604 provides a light-emitting marker 606 .
- the delivery apparatus 604 is in communication with and/or under the control of the computing device 602 .
- the delivery apparatus 604 may include a powered pump or a gravity based arrangement.
- the light-emitting marker 606 may be delivered to the bodily fluid stream 620 via a catheter, an intravenous drip line, a central line delivery, or other needle or port device.
- the delivery apparatus 604 delivers the light-emitting marker 606 in discrete bolus deliveries separated over time.
- the light-emitting marker 606 may include Indocyanine green (ICG), Fluorescein or other glow dye. Two or more dyes, each having a different respective color, may be used. For example, each bolus of two or more may include a different dye and thus the presence and response of each can be determined separately by color distinction.
- the bodily fluid stream 620 in FIG. 6 is represented as flowing along a blood vessel or other biological conduit 622 .
- Several downstream organs or tissue areas 624 , 626 and 628 are represented.
- a region of interest can be selected for observation downstream of the bodily fluid stream 620 carrying with it the light-emitting marker 606 .
- the sensor 610 monitors for fluorescence or other indication of the presence of the light-emitting marker 606 in the selected field of view.
- the sensor 610 can be a high-speed and high-resolution camera for example.
- a first exemplary field of view 630 the sensor 610 observes an area where the bodily fluid stream 620 is divided into several downstream flows.
- the sensor 610 observes an area downstream of the division to isolate monitoring to a single branch of downstream flow.
- the sensor 610 observes a particular organ or tissue area 628 .
- a user such as a physician can deliver a light-emitting marker 606 at any selected location and then observe the time evolving arrival and dispersion or other activity of the marker downstream of the selected location at any selected field of view.
- a central line delivery arrangement is used, in which a port is placed in fluid communication with the subclavian vein and bolus deliveries of the light-emitting marker 606 are injected into the port.
- the delivery apparatus 604 and sensor 610 are shown as connected to the computing device 602 by respective cables 612 and 614 , however wireless connections may be used as well.
- the light-emitting marker 606 briefly fluoresces when excited by an illumination source 640 that emits a particular range of wavelengths upon the region of interest within the field of view of the sensor 610 .
- the illumination source 640 is also shown as connected to the computing device 602 by a cable 642 , however a wireless connection may be used as well.
- the computing device correlates activations of the delivery apparatus 604 , the illumination source 640 , and the sensor 610 , and collects data from the sensor 610 as the light-emitting marker 606 in the field of view responds to the excitation from the illumination source 640 .
- the computing device 602 is configured to record data generated by the sensor 610 ; determine time characteristics of the recorded data; and calculate flow characteristics based on the time characteristics. Further embodiments and examples of fluorescence based imaging and data analysis conducted by the computing device 602 are described above with reference to FIGS. 2-5 .
Abstract
Fluorescence based tracking of a light-emitting marker in a bodily fluid stream is conducted by: providing a light-emitting marker into a fluid stream; establishing field of view monitoring by placement of a sensor, such as a high speed camera, at a region of interest; recording image data of light emitted by the marker at the region of interest; determining time characteristics of the light output of the marker traversing the field of view; and calculating flow characteristics based on the time characteristics. Furthermore generating a velocity vector map may be conducted using a cross correlation technique, leading and falling edge considerations, subtraction, and/or thresholding.
Description
- This application is a continuation U.S. Non-Provisional patent application Ser. No. 15/863,338, entitled “FLUORESCENCE BASED FLOW IMAGING AND MEASUREMENTS,” filed on Jan. 5, 2018, which is a continuation of International Patent Application Serial No. PCT/US2016/041045, entitled “FLUORESCENCE BASED FLOW IMAGING AND MEASUREMENTS,” filed on Jul. 6, 2016, which claims the benefit of priority of U.S. Provisional Patent Application No. 62/189,126, entitled “FLUORESCENCE BASED FLOW IMAGING AND MEASUREMENTS,” filed on Jul. 6, 2015. Each above-reference application is incorporated herein in its entirety by this reference.
- The present disclosure relates to fluorescence based imaging and measurements. More particularly, the present disclosure relates to determining flow characteristics such as velocity in bodily vessels such as blood vessels.
- Fluorescent markers have been used for basic imaging of bodily structures, but improvements are needed in determining flow characteristics in such bodily fluids as blood.
- This summary is provided to introduce in a simplified form concepts that are further described in the following detailed descriptions. This summary is not intended to identify key features or essential features of the claimed subject matter, nor is it to be construed as limiting the scope of the claimed subject matter.
- According to at least one embodiment, a method of fluorescence based tracking of a light-emitting marker in a bodily fluid stream includes: providing a light-emitting marker into a bodily fluid stream; monitoring, with a sensor, a region of interest traversed by the bodily fluid stream; recording data generated by the sensor; determining time characteristics of the recorded data; and calculating flow characteristics based on the time characteristics.
- In at least one example, the sensor includes a camera, and the recorded data comprises motion video data.
- In at least one example, the method further includes: dividing the motion video data into kernels; identifying which of the kernels receive some portion of the light-emitting marker using an intensity threshold; computing, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; performing smoothing on each intensity signal data set; calculating a lag time between the intensity signal data sets of neighboring identified kernels using cross-correlation; using a spatial resolution and the lag time, calculating velocity vectors; summing the velocity vectors of neighboring kernels to create a resultant velocity vector; and generating a velocity map from the resultant velocity vectors for all kernels.
- In at least one example, performing smoothing on each intensity signal data set includes time window averaging.
- In at least one example, performing smoothing on each intensity signal data set includes using a filter.
- In at least one example, wherein performing smoothing on each intensity signal data set includes using a Gaussian filter.
- In at least one example, the method further includes: dividing the motion video data into kernels; identifying which of the kernels receive some portion of the light-emitting marker using an intensity threshold; computing, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; performing smoothing on each intensity signal data set; for each particular identified kernel, finding segments in which a slope of the intensity signal data set rises for a minimum consecutive number of frames or falls for a minimum consecutive number of frames, which segments occur when a leading edge or falling edge of a portion of the light-emitting marker passes through the identified kernel; searching the intensity signal data sets of neighboring identified kernels for a rising or falling segment of similar length; calculating a lag time between segments in the particular identified kernel and segments in the neighboring identified kernels; using a spatial resolution and the lag time, calculating velocity vectors; summing the velocity vectors of neighboring kernels to create a resultant velocity vector; and generating a velocity map from the resultant velocity vectors for all kernels.
- In at least one example, performing smoothing on each intensity signal data set includes time window averaging. In at least one example, performing smoothing on each intensity signal data set includes using a filter. In at least one example, performing smoothing on each intensity signal data set includes using a Gaussian filter.
- In at least one example, the method further includes: calculating a difference frame by subtracting a frame of the motion video data from a consecutive frame of the motion video data; applying a threshold the difference frame to eliminate pixels therein below a specified intensity value; calculating a pixel size of a remaining blob in the difference frame in a direction of blood flow; calculating a size of the remaining blob using the pixel size and a spatial resolution; and calculating a velocity by using a distance traveled by the remaining and a time between frames.
- In at least one example, the method further includes: dividing the motion video data into frames each including pixels; identifying which of the pixels receive some portion of the light-emitting marker using an intensity threshold; creating a logical frame in which a respective indicator for each pixel can be set as true or false; setting the indicators of the identified pixels as true; setting the indicators of all other pixels as false; calculating a difference frame by subtracting a first logical frame from a second logical frame such that the difference frame includes pixels that reached the specified threshold after a time of the first logical frame; finding length in pixels of the remaining blob in the difference frame in a direction of blood flow; converting the length in pixels of the difference frame to physical distance using the spatial resolution; and calculating velocity by dividing the physical distance by a time between frames.
- According to at least one embodiment, a system for fluorescence based tracking of a light-emitting marker in a bodily fluid stream includes: a delivery apparatus configured to provide a light-emitting marker into a bodily fluid stream; a sensor configured to monitor a region of interest traversed by the bodily fluid stream; and a computing device configured to: record data generated by the sensor; determine time characteristics of the recorded data; and calculate flow characteristics based on the time characteristics.
- In at least one example, the sensor includes a camera, and the recorded data includes motion video data.
- In at least one example, the computing device is further configured to: divide the motion video data into kernels; identify which of the kernels receive some portion of the light-emitting marker using an intensity threshold; compute, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; perform smoothing on each intensity signal data set; calculate a lag time between the intensity signal data sets of neighboring identified kernels using cross-correlation; using a spatial resolution and the lag time, calculate velocity vectors; sum the velocity vectors of neighboring kernels to create a resultant velocity vector; and generate a velocity map from the resultant velocity vectors for all kernels.
- In at least one example, the computing device performs smoothing on each intensity signal data set by time window averaging. In at least one example, the computing device performs smoothing on each intensity signal data set by using a Gaussian filter.
- In at least one example, the computing device is further configured to: divide the motion video data into kernels; identify which of the kernels receive some portion of the light-emitting marker using an intensity threshold; compute, for each identified kernel, an intensity signal data set including information of mean light intensity versus time; perform smoothing on each intensity signal data set; for each particular identified kernel, find segments in which a slope of the intensity signal data set rises for a minimum consecutive number of frames or falls for a minimum consecutive number of frames, which segments occur when a leading edge or falling edge of a portion of the light-emitting marker passes through the identified kernel; search the intensity signal data sets of neighboring identified kernels for a rising or falling segment of similar length; calculate a lag time between segments in the particular identified kernel and segments in the neighboring identified kernels; use a spatial resolution and the lag time to calculate velocity vectors; sum the velocity vectors of neighboring kernels to create a resultant velocity vector; and generate a velocity map from the resultant velocity vectors for all kernels.
- In at least one example, the computing device is further configured to: calculate a difference frame by subtracting a frame of the motion video data from a consecutive frame of the motion video data; apply a threshold the difference frame to eliminate pixels therein below a specified intensity value; calculate a pixel size of a remaining blob in the difference frame in a direction of blood flow; calculate a size of the remaining blob using the pixel size and a spatial resolution; and calculate a velocity by using a distance traveled by the remaining blob and a time between frames.
- In at least one example, wherein the computing device is further configured to: divide the motion video data into frames each including pixels; identify which of the pixels receive some portion of the light-emitting marker using an intensity threshold; create a logical frame in which a respective indicator for each pixel can be set as true or false; set the indicators of the identified pixels as true; set the indicators of all other pixels as false; calculate a difference frame by subtracting a first logical frame from a second logical frame such that the difference frame includes pixels that reached the specified threshold after a time of the first logical frame; find length in pixels of the remaining blob in the difference frame in a direction of blood flow; convert the length in pixels of the difference frame to physical distance using the spatial resolution; and calculate velocity by dividing the physical distance by a time between frames.
- The previous summary and the following detailed descriptions are to be read in view of the drawings, which illustrate particular exemplary embodiments and features as briefly described below. The summary and detailed descriptions, however, are not limited to only those embodiments and features explicitly illustrated.
-
FIG. 1 shows a fluorescent marker delivery time plot in a fluid stream and a corresponding response time plot of light intensity measured downstream in a fixed field of view according to at least one embodiment. -
FIG. 2 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using a cross correlation technique. -
FIG. 3 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using leading and falling edge considerations. -
FIG. 4 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using subtraction. -
FIG. 5 shows a flowchart representing a method, according to at least one embodiment, of generating a velocity vector map using thresholding. -
FIG. 6 shows a system, according to at least one embodiment, by which at least the methods ofFIGS. 2-5 are implemented. - These descriptions are presented with sufficient details to provide an understanding of one or more particular embodiments of broader inventive subject matters. These descriptions expound upon and exemplify particular features of those particular embodiments without limiting the inventive subject matters to the explicitly described embodiments and features. Considerations in view of these descriptions will likely give rise to additional and similar embodiments and features without departing from the scope of the inventive subject matters. Although the term “step” may be expressly used or implied relating to features of processes or methods, no implication is made of any particular order or sequence among such expressed or implied steps unless an order or sequence is explicitly stated.
- Any dimensions expressed or implied in the drawings and these descriptions are provided for exemplary purposes. Thus, not all embodiments within the scope of the drawings and these descriptions are made according to such exemplary dimensions. The drawings are not made necessarily to scale. Thus, not all embodiments within the scope of the drawings and these descriptions are made according to the apparent scale of the drawings with regard to relative dimensions in the drawings. However, for each drawing, at least one embodiment is made according to the apparent relative scale of the drawing.
- Fluorescence based tracking according to several embodiments described herein includes the providing of a marker such as a glowing dye into a fluid stream, such as a bloodstream, and making measurements and generating imagery based on the arrival, movement, and departure of the marker downstream as detected by sensor(s) to characterize the flow of the fluid stream and vessels or structures within which the flow travels. The marker is provided into a fluid stream for example by direct injection or via a port as discrete bolus deliveries separated over time. A bolus refers to the administration of a discrete amount of a fluid substance, in this case the marker into a bodily fluid stream such as blood, in order to provide a concentration of the substance to gain a response. A bolus can be delivered by active pumping or by passive gravity based delivery such as via an intravenous drip line. In at least one embodiment, a central line delivery arrangement is used, in which a port is placed in fluid communication with the subclavian vein and bolus deliveries are injected into the port. The dye briefly fluoresces when excited by an illumination source that emits a particular range of wavelengths. The dye is illuminated over the Region of Interest (ROI) where imaging of the fluorescence is also performed.
- A field of view monitoring is established by placement of a sensor, such as a high speed camera, at a region of interest. The field of view can be established and held generally fixed as the marker enters and traverses the field of view of a high-speed camera sensitive to the light emitted by the marker. Time characteristics of the light output of the marker traversing the field of view can be deducted from the light output intensity as recorded by the camera. A field of view may be established for example at the heart or other organ where flow diagnostics are wanted.
- The visual response in the field of view indicates presence of the marker, with the intensity of the light response being correlated with the time evolving concentration of the marker in the stream as the marker diffuses and travels with the host fluid. The light intensity in the field of view may typically have both rise and fall characteristics. The rise characteristics correspond to the arrival and increasing concentration of the marker in the field of view. The fall characteristics correspond to the departure or diffusion of the marker and/or the reduction of its light output. In the case of a dye marker in a blood stream as injected by a bolus, rise time may be faster generally than fall time such that response time curves typically show steeper climbs than falls.
-
FIG. 1 shows a fluorescent markerdelivery time plot 102 in a fluid stream and a correspondingresponse time plot 112 of light intensity measured downstream in a fixed field of view according to at least one embodiment. Theplots common time axis 100. Thedelivery time plot 102 recordsseveral bolus deliveries 104 as step functions separated in time. Theresponse plot 112 records the corresponding response functions 114 of light intensity in a field of view downstream from the marker delivery location into a fluid stream. Eachresponse function 114 is shown as a wave or envelope having arise side 116 representing the arrival of the marker in the field of view, and afall side 118 indicating the departure or diffusion of the marker and/or the reduction of its light output. In order to correlate bolus deliveries and data acquisition, a marker delivery system in at least one embodiment includes a controller and a delivery pump in wired or wireless communication. A data acquisition system including the camera and a computing device for recording, analyzing and visualizing the data and/or field of view are also in communication with the controller in at least one embodiment. - The time of initiation, delivered volume, and duration of each bolus delivery can be controlled. The time interval between consecutive bolus deliveries is also controlled. Thus, multiple parameters for bolus delivery can be adjusted to ultimately query and determine varied flow characteristics within a region of interest subjected to field of view monitoring. Shorter time gaps can be used for slower moving fluids and longer time gaps can be used for faster moving fluids within the region of interest. Injection times can be varied to address different anatomical depths and tissue surface barriers.
- In at least one embodiment, as the marker from a bolus delivery enters an established field of view, the light response of a bolus is captured by a high speed camera. The time domain dynamics of the light response is analyzed to arrive at velocity vectors representing movement of the host fluid in the field of view. Several embodiments of generating velocity vectors using the data from fluorescence based tracking are described in the following with reference to
FIGS. 2-5 . In each, a video including multiple frames (images) taken of a region of interest is analyzed to track pixels in or across the field of view. - A
method 200 of generating a velocity vector map using a cross correlation technique, according to at least one embodiment, is represented as a flow chart inFIG. 2 . Instep 202, divide the video into spatiotemporal cubes, which are termed “kernels” in these descriptions. Instep 204, which is optional, isolate kernels in which fluorescence appears at any point during the length of the video using an intensity threshold. Instep 206, compute a signal for each kernel of mean intensity vs. time. The signal may be 1D, 2D, or 3D. Vessel(s) may cross under other vessel(s) and at different angles, and multiple perspectives from different cameras can be used. The potential use of coded apertures is within the scope of these descriptions. Instep 208, which is optional, perform smoothing on the intensity signal for each kernel using time window averaging, Gaussian filter, etc. Instep 210, for each kernel, calculate the lag time between the intensity signal of the kernel and the time intensity signals of its neighboring kernels using cross-correlation. Instep 212, using the known spatial resolution of the image, convert lag time in each direction into velocity. Instep 214, sum the resulting velocity vectors to each neighboring kernel to create one resultant velocity vector. Instep 216, generate a velocity map from the resultant velocity vectors for all kernels. - A
method 300 of generating a velocity vector map using leading and falling edge considerations, according to at least one embodiment, is represented as a flow chart inFIG. 3 . Instep 302, divide video into kernels. Instep 304, which is optional, isolate kernels in which fluorescence appears at any point during the length of the video using an intensity threshold. Instep 306, compute a signal for each kernel of mean intensity vs. time. As described above with reference to step 206 ofFIG. 2 , the signal may be 1D, 2D, or 3D. Instep 308, which is optional, perform smoothing on the intensity signal for each kernel using time window averaging, Gaussian filter, etc. Instep 310, for each kernel, find segments in which the slope of the intensity signal rises for a minimum consecutive number of frames or falls for a minimum consecutive number of frames. These segments occur when the leading or falling edges of the ICG (fluorescein or other glow dye) bolus pass through the kernel. Instep 312, search the intensity signals of neighboring kernels for a rising or falling segment of similar length. Instep 314, calculate the temporal offset (lag time) between the segment in the original kernel and the segment in the neighboring kernels. Instep 316, using the known spatial resolution of the image, convert lag time in each direction into velocity. Instep 318, sum the resulting velocity vectors to each neighboring kernel to create one resultant velocity vector. Instep 320, generate a velocity map from the resultant velocity vectors for all kernels. - A
method 400 of generating a velocity vector map using subtraction, according to at least one embodiment, is represented as a flow chart inFIG. 4 . Instep 402, subtract two consecutive frames, resulting in a difference frame. Instep 404, threshold difference frame to eliminate pixels below a specified intensity value. Instep 406, calculate pixel size of a remaining blob in the difference frame in the direction of blood flow. Instep 408, convert size of blob in pixels to physical distance using spatial resolution. Instep 410, calculate velocity by dividing distance traveled by time between frames. - A
method 500 of generating a velocity vector map using thresholding, according to at least one embodiment, is represented as a flow chart inFIG. 5 . Instep 502, for each frame in video, isolate pixels with intensities above specified threshold. Create logical frame and set all pixels at or above threshold to true. Set all other pixels to false. Instep 504, for two consecutive frames, subtract the first logical frame from the second logical frame. The resulting frame contains the pixels that reached the specified threshold between frames. Instep 506, find the length in pixels of the remaining blob in the difference frame in the direction of blood flow. Instep 508, convert the pixel length of the difference frame to physical distance using the spatial resolution. Instep 510, calculate velocity by dividing physical distance by the time between frames. -
FIG. 6 shows asystem 600, according to at least one embodiment, by which at least the methods ofFIGS. 2-5 are implemented. Thesystem 600 includes acomputing device 602, adelivery apparatus 604 configured to provide a light-emittingmarker 606 into a bodily fluid stream, and asensor 610 configured to monitor a region of interest traversed by the bodily fluid stream. Thecomputing device 602 records data generated by thesensor 610, determines time characteristics of the recorded data; and calculates flow characteristics based on the time characteristics. - The
computing device 602 is illustrated as a laptop or other personal computer. Other computing devices including local and remote servers are within the scope of these descriptions. Thedelivery apparatus 604 provides a light-emittingmarker 606. Thedelivery apparatus 604 is in communication with and/or under the control of thecomputing device 602. Thedelivery apparatus 604 may include a powered pump or a gravity based arrangement. The light-emittingmarker 606 may be delivered to thebodily fluid stream 620 via a catheter, an intravenous drip line, a central line delivery, or other needle or port device. Thedelivery apparatus 604 delivers the light-emittingmarker 606 in discrete bolus deliveries separated over time. The light-emittingmarker 606 may include Indocyanine green (ICG), Fluorescein or other glow dye. Two or more dyes, each having a different respective color, may be used. For example, each bolus of two or more may include a different dye and thus the presence and response of each can be determined separately by color distinction. - The
bodily fluid stream 620 inFIG. 6 is represented as flowing along a blood vessel or otherbiological conduit 622. Several downstream organs ortissue areas sensor 610, a region of interest can be selected for observation downstream of thebodily fluid stream 620 carrying with it the light-emittingmarker 606. Thesensor 610 monitors for fluorescence or other indication of the presence of the light-emittingmarker 606 in the selected field of view. Thesensor 610 can be a high-speed and high-resolution camera for example. - Several fields of view are represented. In a first exemplary field of
view 630, thesensor 610 observes an area where thebodily fluid stream 620 is divided into several downstream flows. In a second exemplary field ofview 632, thesensor 610 observes an area downstream of the division to isolate monitoring to a single branch of downstream flow. In a third exemplary field ofview 634, thesensor 610 observes a particular organ ortissue area 628. These examples represent that a user such as a physician can deliver a light-emittingmarker 606 at any selected location and then observe the time evolving arrival and dispersion or other activity of the marker downstream of the selected location at any selected field of view. In at least one embodiment, a central line delivery arrangement is used, in which a port is placed in fluid communication with the subclavian vein and bolus deliveries of the light-emittingmarker 606 are injected into the port. - The
delivery apparatus 604 andsensor 610 are shown as connected to thecomputing device 602 byrespective cables marker 606 briefly fluoresces when excited by anillumination source 640 that emits a particular range of wavelengths upon the region of interest within the field of view of thesensor 610. Theillumination source 640 is also shown as connected to thecomputing device 602 by acable 642, however a wireless connection may be used as well. The computing device correlates activations of thedelivery apparatus 604, theillumination source 640, and thesensor 610, and collects data from thesensor 610 as the light-emittingmarker 606 in the field of view responds to the excitation from theillumination source 640. - In various embodiments, the
computing device 602 is configured to record data generated by thesensor 610; determine time characteristics of the recorded data; and calculate flow characteristics based on the time characteristics. Further embodiments and examples of fluorescence based imaging and data analysis conducted by thecomputing device 602 are described above with reference toFIGS. 2-5 . - Particular embodiments and features have been described with reference to the drawings. It is to be understood that these descriptions are not limited to any single embodiment or any particular set of features, and that similar embodiments and features may arise or modifications and additions may be made without departing from the scope of these descriptions and the spirit of the appended claims.
Claims (10)
1. A system for fluorescence based tracking of a light-emitting marker in a bodily fluid stream, the system comprising:
a delivery apparatus configured to provide a light-emitting marker into the bodily fluid stream;
a camera configured to monitor a region of interest traversed by the bodily fluid stream; and
a computing device configured to:
record motion video data generated by the camera;
determine time characteristics of the recorded data; and
calculate flow characteristics based on the time characteristics.
2. The system according to claim 1 , wherein the computing device is further configured to:
divide the motion video data into kernels;
identify which of the kernels receive some portion of the light-emitting marker using an intensity threshold;
compute, for each identified kernel, an intensity signal data set comprising information of mean light intensity versus time;
perform smoothing on each intensity signal data set; and
calculate a lag time between the intensity signal data sets of neighboring identified kernels using cross-correlation.
3. The system according to claim 1 , wherein the computing device is further configured to:
using a spatial resolution and the lag time, calculate velocity vectors;
sum the velocity vectors of neighboring kernels to create a resultant velocity vector; and
generate a velocity map from the resultant velocity vectors for all kernels.
4. The system according to claim 3 , wherein the computing device performs smoothing on each intensity signal data set by time window averaging or by using a filter.
5. The system according to claim 3 , wherein the computing device is further configured to:
for each particular identified kernel, find segments in which a slope of the intensity signal data set rises for a minimum consecutive number of frames or falls for a minimum consecutive number of frames, which segments occur when a leading edge or falling edge of a portion of the light-emitting marker passes through the identified kernel;
search the intensity signal data sets of neighboring identified kernels for a rising or falling segment of similar length; and
calculate a lag time between segments in the particular identified kernel and segments in the neighboring identified kernels.
6. The system according to claim 3 , wherein the computing device is further configured to:
calculate a difference frame by subtracting a frame of the motion video data from a consecutive frame of the motion video data;
apply a threshold to the difference frame to eliminate pixels therein below a specified intensity value;
calculate a pixel size of a remaining blob in the difference frame in a direction of bodily fluid flow;
calculate a size of the remaining blob using the pixel size and a spatial resolution; and
calculate a velocity by using a distance traveled by the remaining blob and a time between frames.
7. The system according to claim 3 , wherein the computing device is further configured to:
create a logical frame in which a respective indicator for each pixel can be set as true or false;
set the indicators of the identified pixels as true;
set the indicators of all other pixels as false;
calculate a difference frame by subtracting a first logical frame from a second logical frame such that the difference frame comprises pixels that reached the specified threshold after a time of the first logical frame;
find length in pixels of the remaining blob in the difference frame in a direction of bodily fluid flow;
convert the length in pixels of the difference frame to physical distance using the spatial resolution; and
calculate velocity by dividing the physical distance by a time between frames.
8. A method of fluorescence based tracking of a light-emitting marker in a fluid stream, the method comprising:
monitoring, with a camera, a region of interest traversed by fluid stream into which a light-emitting marker has been introduced;
recording motion video data generated by the sensor;
dividing the motion video data into frames each comprising pixels;
identifying which of the pixels receive some portion of the light-emitting marker using an intensity threshold;
calculating a difference frame by subtracting a frame of the motion video data from a consecutive frame of the motion video data; and
applying a threshold to the difference frame to eliminate pixels therein below a specified intensity value.
9. The method according to claim 8 , further comprising:
calculating a pixel size of a remaining blob in the difference frame in a direction of fluid flow;
calculating a size of the remaining blob using the pixel size and a spatial resolution; and
calculating a velocity by using a distance traveled by the remaining blob and a time between frames.
10. The method according to claim 9 , further comprising:
creating a logical frame in which a respective indicator for each pixel can be set as true or false;
setting the indicators of the identified pixels as true;
setting the indicators of all other pixels as false;
calculating a difference frame by subtracting a first logical frame from a second logical frame such that the difference frame comprises pixels that reached the specified threshold after a time of the first logical frame;
finding length in pixels of a remaining blob in the difference frame in a direction of fluid flow;
converting the length in pixels of the difference frame to physical distance using the spatial resolution; and
calculating velocity by dividing the physical distance by a time between frames.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US17/031,375 US20210000352A1 (en) | 2015-07-06 | 2020-09-24 | Fluorescence based flow imaging and measurements |
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201562189126P | 2015-07-06 | 2015-07-06 | |
PCT/US2016/041045 WO2017007791A1 (en) | 2015-07-06 | 2016-07-06 | Fluorescence based flow imaging and measurements |
US15/863,338 US10813563B2 (en) | 2015-07-06 | 2018-01-05 | Fluorescence based flow imaging and measurements |
US17/031,375 US20210000352A1 (en) | 2015-07-06 | 2020-09-24 | Fluorescence based flow imaging and measurements |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/863,338 Continuation US10813563B2 (en) | 2015-07-06 | 2018-01-05 | Fluorescence based flow imaging and measurements |
Publications (1)
Publication Number | Publication Date |
---|---|
US20210000352A1 true US20210000352A1 (en) | 2021-01-07 |
Family
ID=57686151
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/863,338 Active 2037-07-12 US10813563B2 (en) | 2015-07-06 | 2018-01-05 | Fluorescence based flow imaging and measurements |
US17/031,375 Abandoned US20210000352A1 (en) | 2015-07-06 | 2020-09-24 | Fluorescence based flow imaging and measurements |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/863,338 Active 2037-07-12 US10813563B2 (en) | 2015-07-06 | 2018-01-05 | Fluorescence based flow imaging and measurements |
Country Status (5)
Country | Link |
---|---|
US (2) | US10813563B2 (en) |
EP (1) | EP3319515B1 (en) |
CA (1) | CA3028084A1 (en) |
ES (1) | ES2798759T3 (en) |
WO (1) | WO2017007791A1 (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2648029C2 (en) * | 2016-08-10 | 2018-03-21 | Самсунг Электроникс Ко., Лтд. | Device and method of blood pressure measurement |
BR112022011870A2 (en) * | 2019-12-19 | 2022-09-06 | Perfusion Tech Aps | METHOD IMPLEMENTED BY COMPUTER, COMPUTER PROGRAM, AND SYSTEM TO IDENTIFY BLOOD VESSELS |
CN112232223A (en) * | 2020-10-19 | 2021-01-15 | 北京百度网讯科技有限公司 | Multi-modal image processing method, device, equipment and storage medium |
Family Cites Families (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050182434A1 (en) * | 2000-08-11 | 2005-08-18 | National Research Council Of Canada | Method and apparatus for performing intra-operative angiography |
DE10120980B4 (en) * | 2001-05-01 | 2009-12-03 | Pulsion Medical Systems Ag | A method, apparatus and computer program for determining blood flow in a tissue or organ region |
AU2002216861A1 (en) * | 2001-12-19 | 2003-06-30 | Institut De Cardiologie De Montreal | Non-invasive detection of endothelial dysfunction by blood flow measurement in opposed limbs |
US7697972B2 (en) * | 2002-11-19 | 2010-04-13 | Medtronic Navigation, Inc. | Navigation system for cardiac therapies |
EP1713400B1 (en) * | 2004-01-16 | 2012-12-26 | Boston Scientific Scimed, Inc. | Apparatus for medical imaging |
US7491502B2 (en) * | 2004-12-17 | 2009-02-17 | The General Hospital Corporation | In vivo flow cytometry system and method |
US20080019921A1 (en) * | 2006-06-30 | 2008-01-24 | Invitrogen Corporation | Uniform fluorescent microsphere with hydrophobic surfaces |
US8155730B2 (en) * | 2006-10-24 | 2012-04-10 | The Research Foundation Of State University Of New York | Composition, method, system, and kit for optical electrophysiology |
CA2682275C (en) * | 2007-03-28 | 2017-05-09 | Bionanomatrix, Inc. | Methods of macromolecular analysis using nanochannel arrays |
US8285366B2 (en) * | 2007-10-04 | 2012-10-09 | The Invention Science Fund I, Llc | Vasculature and lymphatic system imaging and ablation associated with a local bypass |
US9144383B2 (en) * | 2007-12-13 | 2015-09-29 | The Board Of Trustees Of The University Of Arkansas | Device and method for in vivo noninvasive magnetic manipulation of circulating objects in bioflows |
US8406860B2 (en) * | 2008-01-25 | 2013-03-26 | Novadaq Technologies Inc. | Method for evaluating blush in myocardial tissue |
US9173554B2 (en) * | 2008-03-18 | 2015-11-03 | Novadaq Technologies, Inc. | Imaging system for combined full-color reflectance and near-infrared imaging |
US8574859B2 (en) * | 2008-03-27 | 2013-11-05 | The General Hospital Corporation | In vivo flow cytometry based on cellular autofluorescence |
RU2506896C2 (en) * | 2008-04-14 | 2014-02-20 | Новадак Текнолоджиз Инк. | Determination of location and analysis of perforated flaps for plastic and rehabilitative surgery |
US10219742B2 (en) * | 2008-04-14 | 2019-03-05 | Novadaq Technologies ULC | Locating and analyzing perforator flaps for plastic and reconstructive surgery |
DE102008040803A1 (en) * | 2008-07-28 | 2010-02-04 | Carl Zeiss Surgical Gmbh | Method for the quantitative representation of the blood flow |
JP5395396B2 (en) * | 2008-10-15 | 2014-01-22 | 株式会社東芝 | Ultrasonic diagnostic apparatus, medical image processing apparatus, and medical image processing program |
DE102009010446B4 (en) * | 2009-02-26 | 2015-07-09 | Carl Zeiss Meditec Ag | Method and apparatus for quantifying blood flow in a tissue area |
US8615311B2 (en) * | 2010-03-15 | 2013-12-24 | Rutgers, The State University Of New Jersey | Microelectorode array, methods for preparing the same and uses thereof |
JP2014524121A (en) * | 2011-07-14 | 2014-09-18 | ザ・ジョージ・ワシントン・ユニバーシティ | Plume collimation for laser ablation and electrospray ionization mass spectrometry |
US9523682B2 (en) * | 2011-11-16 | 2016-12-20 | Becton, Dickinson And Company | Methods and systems for detecting an analyte in a sample |
WO2013116472A1 (en) * | 2012-01-31 | 2013-08-08 | Massachusetts Institute Of Technology | Imaging probe including nanoparticle |
US10180398B2 (en) * | 2012-02-29 | 2019-01-15 | National Technology & Engineering Solutions Of Sandia, Llc | Trajectory-based triggering system for hyperspectral imaging flow cytometer |
US9501620B2 (en) | 2012-03-20 | 2016-11-22 | Oxford University Innovation Limited | Quantification of blood volume flow rates from dynamic angiography data |
WO2014081419A2 (en) * | 2012-11-20 | 2014-05-30 | The Board Of Trustees Of The Leland Stanford Junior University | High resolution imaging using near-infrared-ii fluorescence |
US8977331B2 (en) * | 2012-12-13 | 2015-03-10 | General Electric Company | Systems and methods for nerve imaging |
WO2015166503A1 (en) * | 2014-05-01 | 2015-11-05 | Yeda Research And Development Co. Ltd. | Multimodal transcranial brain optical imaging |
-
2016
- 2016-07-06 CA CA3028084A patent/CA3028084A1/en not_active Abandoned
- 2016-07-06 EP EP16821893.1A patent/EP3319515B1/en active Active
- 2016-07-06 WO PCT/US2016/041045 patent/WO2017007791A1/en active Application Filing
- 2016-07-06 ES ES16821893T patent/ES2798759T3/en active Active
-
2018
- 2018-01-05 US US15/863,338 patent/US10813563B2/en active Active
-
2020
- 2020-09-24 US US17/031,375 patent/US20210000352A1/en not_active Abandoned
Also Published As
Publication number | Publication date |
---|---|
CA3028084A1 (en) | 2017-01-12 |
WO2017007791A1 (en) | 2017-01-12 |
ES2798759T3 (en) | 2020-12-14 |
EP3319515A4 (en) | 2019-03-13 |
EP3319515B1 (en) | 2020-03-18 |
US20180125378A1 (en) | 2018-05-10 |
US10813563B2 (en) | 2020-10-27 |
EP3319515A1 (en) | 2018-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20210000352A1 (en) | Fluorescence based flow imaging and measurements | |
US11564583B2 (en) | Method for evaluating blush in myocardial tissue | |
CA2750760C (en) | Method for evaluating blush in myocardial tissue | |
US20230355121A1 (en) | Facilitating assessment of blood flow and tissue perfusion using fluorescence-mediated photoplethysmography | |
CN107361791B (en) | Rapid super-resolution blood flow imaging method | |
EP3513381A1 (en) | Systems and methods for generating 3d images based on fluorescent illumination | |
EP3488417B1 (en) | A system and method for automated characterization of solid tumors using medical imaging | |
KR20150128405A (en) | Method, apparatus and system for providing medical image | |
WO2021191909A1 (en) | Functional measurements of vessels using a temporal feature | |
CN110636798A (en) | Method and apparatus for physiological function parameter determination | |
Guo et al. | Frame rate effects and their compensation on super-resolution microvessel imaging using ultrasound localization microscopy | |
US10803601B2 (en) | Rapid assessment and visual reporting of local particle velocity | |
Lo et al. | Pseudo three-dimensional vision-based nail-fold morphological and hemodynamic analysis | |
Kurochkin et al. | Toward label-free imaging of brain vasculature: frame-by-frame spatial adaptive filtration and adaptive PIV approaches | |
You et al. | A novel line detection method in space-time images for microvascular blood flow analysis in sublingual microcirculatory videos | |
US20220211258A1 (en) | Contrast enhancement for medical imaging | |
Tichauer et al. | A dual-reporter fluorescent imaging approach can be used to estimate sentinel lymph node tumor burden | |
Bárdossy et al. | Determination of venous blood flow velocity using digital videomicroscopy (A short methodical communication) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: SCINOVIA CORP., NORTH CAROLINA Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:SUND, JAMES BRADLEY, SR.;COHEN, DAVID S.;REEL/FRAME:053876/0263 Effective date: 20160907 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: APPLICATION DISPATCHED FROM PREEXAM, NOT YET DOCKETED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: DOCKETED NEW CASE - READY FOR EXAMINATION |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |