US20120323109A1 - Photon Measurement Method and Apparatus - Google Patents
Photon Measurement Method and Apparatus Download PDFInfo
- Publication number
- US20120323109A1 US20120323109A1 US13/598,324 US201213598324A US2012323109A1 US 20120323109 A1 US20120323109 A1 US 20120323109A1 US 201213598324 A US201213598324 A US 201213598324A US 2012323109 A1 US2012323109 A1 US 2012323109A1
- Authority
- US
- United States
- Prior art keywords
- signal
- pseudorandom binary
- detector
- optical
- coupled
- 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
- 238000000691 measurement method Methods 0.000 title description 2
- 238000000034 method Methods 0.000 claims abstract description 19
- 230000003287 optical effect Effects 0.000 claims description 76
- 238000005259 measurement Methods 0.000 claims description 24
- 238000005070 sampling Methods 0.000 claims description 19
- 230000002123 temporal effect Effects 0.000 claims description 18
- 239000000463 material Substances 0.000 claims description 16
- 238000012545 processing Methods 0.000 claims description 5
- 238000010521 absorption reaction Methods 0.000 claims description 3
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 claims description 2
- 239000010931 gold Substances 0.000 claims description 2
- 229910052737 gold Inorganic materials 0.000 claims description 2
- 238000009826 distribution Methods 0.000 abstract description 9
- 239000000523 sample Substances 0.000 description 57
- 230000004044 response Effects 0.000 description 24
- 230000002596 correlated effect Effects 0.000 description 10
- 238000001514 detection method Methods 0.000 description 10
- 238000005286 illumination Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000000835 fiber Substances 0.000 description 6
- 238000012546 transfer Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 239000013307 optical fiber Substances 0.000 description 5
- 210000000481 breast Anatomy 0.000 description 4
- 230000003111 delayed effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000001934 delay Effects 0.000 description 3
- 238000002596 diffuse optical imaging Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 239000000872 buffer Substances 0.000 description 2
- 230000000875 corresponding effect Effects 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 239000002872 contrast media Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000009543 diffuse optical tomography Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000009607 mammography Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000006833 reintegration Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/62—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
- G01N21/63—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
- G01N21/64—Fluorescence; Phosphorescence
- G01N21/645—Specially adapted constructive features of fluorimeters
- G01N21/6456—Spatial resolved fluorescence measurements; Imaging
-
- 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
-
- 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/0082—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes
- A61B5/0084—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for introduction into the body, e.g. by catheters
- A61B5/0086—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence adapted for particular medical purposes for introduction into the body, e.g. by catheters using infrared radiation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4795—Scattering, i.e. diffuse reflection spatially resolved investigating of object in scattering medium
-
- 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
Definitions
- the field of the present invention pertains generally to systems and methods for detecting and measuring photons, including, more specifically, the measurement of the time-of-flight of photons traveling through a scattering media such as tissue.
- Diffuse optical imaging techniques are known in medical and biological applications. Overviews of diffuse optical imaging techniques can be found in “Recent Advances in Diffusion Optical Imaging” by Gibson, et al, Phys. Med. Biology, vol. 50 (2005), R1-R43 and in “Near-infrared Diffuse Optical Tomography,” by Hielscher, et al, Disease Markers, Vol. 18 (2002), 313-337. Briefly, diffuse optical imaging involves the use of near-infrared light incident upon a sample of interest. An example in the medical and biological field is optical mammography where near infrared light is used to illuminate breast tissue. A detector is placed on the opposite side of the breast from the incident light some distance away and collects scattered light from the breast tissue.
- the scattered light of interest that is detected may be directly scattered incident light or scattered fluorescence light caused by the excitation of an injected fluorescing material that fluoresces when exposed to the incident light.
- Measurements of the photon flight-time distributions are typically carried out using either the time-domain or the frequency-domain technique.
- the time-domain technique the sample is excited with a pulse of light from a pulsed laser and the scattered light is measured using a detector with single-photon sensitivity.
- the detector measures the time delay between the excitation pulse and the first detected photon.
- the flight-time distribution is determined by using many repeated pulses and building up a histogram of the measured time delays.
- the pulsed laser sources and single-photon detectors are relatively expensive. Because detection is typically done at the single-photon level, it can require a significant amount of time to build-up enough data to approximate the flight-time distribution.
- One disadvantage of the frequency-domain approach is that it is not a direct measurement of the photon flight time.
- U.S. Pat. No. 5,565,982 discloses a time-resolved spectroscopy system using digital processing techniques and two low power, continuous wave light sources.
- the disclosed system requires two light transmitters of different wavelengths modulated with separate codes for interrogating a sample of interest. Properties of the sample are inferred by differential comparison of the return signals from each of the two light sources. It is undesirable to have two distinct light sources due to the cost and complexity involved. Furthermore, the noise level associated with a measurement made with two separate light sources will be higher than with a single source even if the codes used to drive the two sources are orthogonal.
- a system and method capable of addressing these disadvantages while providing acceptable temporal information for whatever application the photon measurement is being used is needed.
- a preferred system uses information from a single low-power continuous-wave light source detected by a low-cost detector.
- the measurement system is implemented with digital electronics.
- One embodiment of the system and methods disclosed comprises an information from a single continuous-wave light source modulated with a digital waveform for interrogating a sample, a photo-sensitive detector for measuring the scattered light from the sample, and electronics for sampling the detector output and performing a correlation of the output signal with the modulation waveform.
- Other embodiments include electronics and software for determining the parameters of the flight-time distribution from the measured correlation.
- FIG. 1 is a functional block diagram of the major components of a preferred photon measurement system of the present invention.
- FIG. 2 is a diagram of preferred Analog-to-Digital converters and their interface to the signal detector.
- FIG. 3 is a functional block diagram of a preferred signal generator.
- FIG. 4 depicts an implementation of a preferred Linear Feedback Shift Register.
- FIG. 5 is a functional block diagram of a preferred signal detector.
- FIG. 6 is a functional block diagram of a preferred frame accumulator.
- FIG. 7 is a functional block diagram of a preferred frame correlator.
- FIG. 8 is a mechanical view of a preferred embodiment.
- FIG. 9 is a functional representation of a preferred embodiment using a needle-based optical probe coupled to a source and a detector with a single optical fiber.
- FIG. 10 depicts a preferred needle-based probe.
- FIG. 1 A functional block diagram of a preferred photon measurement system 100 is depicted in FIG. 1 .
- the photon measurement system can be used to measure the interaction of photons with a sample 5 .
- the sample 5 may be human breast tissue or fat tissue but it could just as well be any semitransparent material.
- the photon measurement system 100 preferably includes Temporal Response Analysis Engine 11 .
- the Temporal Response Analysis Engine 11 generates a digital modulation signal for driving an illumination light source that is used to interrogate the sample.
- the Temporal Response Engine 11 also provides a means for processing a detected optical signal from the sample 5 to extract information about the sample 5 .
- a digital modulation signal 16 is generated in the signal generator 1 and transmitted to the transmit signal conditioner 2 .
- the digital modulation signal 16 is the digital representation of a chosen code sequence.
- the code sequence is preferably chosen from the known pseudorandom binary sequences (PRBS), Gold codes, Golay codes, Kasami codes, Walsh codes, or other codes that possess the preferred desirable property of large auto-correlation values and low cross-correlation values.
- PRBS pseudorandom binary sequences
- the digital modulation signal 16 may represent a single code pattern or multiple repeats of the same pattern. A single complete set of code patterns is designated a modulation frame or code pattern frame.
- the digital modulation signal 16 is preferably transmitted to the signal detector 10 as an electronic reference signal 17 .
- the transmit signal conditioner 2 formats the digital modulation signal 16 as necessary to drive an optical illumination source 3 .
- the optical illumination source is a 785 nm continuous wave diode laser made by Hitachi Corp. Formatting of the digital modulation signal 16 in the preferred embodiment involves converting the digital modulation signal 16 to an analog voltage waveform that is coupled through a 50-ohm bias-T to the DC drive current of the optical illumination source 3 .
- the optical illumination source 3 may be a different laser diode, a light-emitting diode, or a light source used together with an external optical modulator.
- the optical illumination source 3 generates the modulated optical wave 20 which is preferably transmitted to the sample 5 by light delivery optics 4 .
- the preferred light delivery optics 4 is a 3 mm diameter fiber bundle located between the optical illumination source 3 and the sample 5 to deliver the modulated optical wave 20 from the optical illumination source 3 to the sample 5 .
- the light delivery optics 4 comprises other arrangements of optical fibers, lenses, mirrors or other optical delivery components.
- the fluorescent material is preferably an exogenous contrast agent injected into the sample 5 or alternatively it is preferably some constituent component of a material that exhibits endogenous fluorescence.
- the detection optics 6 are situated so that a portion of the modulated optical waves 21 are detected by the detection optics 6 .
- the detection optics 6 include an optical filter for separating the fluorescing scattered optical waves 21 from the modulated optical waves 20 .
- the optical filter preferably transmits the higher wavelength fluorescence and blocks the lower wavelength illumination light. In applications where the scattered optical waves 21 of interest are not fluorescing, an optical filter is not required.
- the detection optics 6 preferably include a second 3 mm diameter fiber bundle located between the optical filter and the optical detector 7 .
- the optical detector 7 converts the scattered optical waves 21 to an electronic signal.
- the optical detector 7 is preferably a photomultiplier tube, model R7400U-20 from Hamamatsu Corp.
- the optical detector 7 may be a PIN photodiode, an avalanche photodiode, a charge-couple device, or other suitable photosensitive element.
- the optical detector 7 preferably converts detected scattered optical waves 21 into an electronic signal which is communicated to the detected signal conditioner 8 .
- the detected signal conditioner 8 preferably formats the signal so it may be converted to discrete samples by an Analog to Digital (A/D) converter 9 .
- the A/D converter 9 outputs a detected response signal 19 .
- the detected response signal 19 is communicated to a signal detector 10 , where it is preferably correlated with the electronic reference signal 17 to extract a sample transfer characteristic.
- Information about the temporal properties of the photons is preferably calculated from the sample transfer characteristic. This information preferably includes such properties as direct measurements of photon time-of-flight and the fluorescence lifetime. The estimate of photon times-of-flight is then preferably used to estimate characteristics of the tissue such as the absorption coefficient, scattering coefficient, or location of fluorescing material.
- the photon measurement system 100 includes an optical reference generator 22 .
- the optical reference generator 22 preferably includes an optical splitter 12 A or 12 B that routes a portion of the modulated optical wave 20 to a secondary optical detector 13 .
- the position of the optical splitter 12 A or 12 B can be either before or after the light delivery optics.
- the output of the secondary optical detector 13 is preferably routed to a secondary signal conditioner 14 whose output is communicated to a secondary A/D converter 15 .
- the secondary A/D converter 15 preferably outputs a source reference signal 18 which can be correlated with the detected response 19 to extract the sample transfer characteristic.
- Using the source reference signal 18 as opposed to the electronic reference signal 17 allows the filtering of the temporal properties of the signal conditioner 2 and the modulated optical source 3 from the measured transfer characteristic.
- FIG. 2 The preferred hardware implementation of the A/D converter module and its interfaces to the signal detector 10 are shown in FIG. 2 .
- An array of N A/D converters 90 preferably receives the analog signal 95 in parallel from the signal conditioner 8 or 14 .
- the output samples 18 or 19 from the A/D converters 90 are preferably communicated to the First-In-First-Out buffers (FIFOs) 91 where they are buffered for distribution to the internal components of the signal detector 10 .
- the A/D converters 90 are eight MAX 108 integrated circuits made by Maxim operating at 250 Msample/sec and outputting two data samples at a time in parallel at 125 MHz.
- the FIFOs 91 are preferably implemented within a Xilinx 4 FPGA.
- the acquisition synchronizer 92 preferably controls signal acquisition and digital data distribution through the conversion clock (CClk) signals 96 .
- the acquisition synchronizer 92 is preferably synchronized with an externally provided synchronization clock (SClk) 40 which is also preferably used to synchronize the signal generator 1 .
- SClk synchronization clock
- the signals CClk[ 1 . . . N] are preferably generated within the acquisition synchronizer 92 and preferably have the same frequency as SClk 40 but are offset in phase from SClk 40 in N fixed increments of (360/N)°, with the phase of CClk[ 1 ] set to the fixed offset of Z°.
- the internal clock generation capabilities of the Xilinx FPGA are used to implement the acquisition synchronizer 92 directly.
- the A/D converters 90 preferably perform their conversions in sync with the conversion clocks 96 such that they generate samples at N discrete sample times spread evenly throughout the fundamental sample interval defined by the period of SClk 40 .
- the effective sample rate for the array of converters is preferably N times the rate defined by SClk 40 . This process of using multiple A/D converters sampling out of phase to increase the effective sample rate is what we call parallel over-sampling. In the preferred photon measuring system, parallel over-sampling results in an effective sample rate of 2 Gsamples/sec.
- the offset value Z allows the entire sample set to be offset by some phase from the synchronization clock 40 .
- the acquisition synchronizer 92 preferably is configured such that the value of Z can be varied synchronously with the modulation frame, or with a block of frames called a frame block. This allows Z to follow a sequence of K values smaller than (360/N)° such that on successive modulation frames/frame blocks the effective sampling phases (relative to the synchronization clock) take on K values intermediate to those created by the N conversion clocks in any given frame.
- the input signal at any given A/D converter 90 will be sampled at K discrete phases over K blocks.
- the detected response 19 is preferably assumed to be stationary with respect to the start of the code pattern block over that time interval.
- the preferred K discrete sampling phases correspond to K discrete sample times and the effective temporal resolution of the sampling process is preferably increased by a factor of K. This process is referred to as temporal over-sampling.
- the value of Z is always zero and temporal over-sampling is achieved by adjusting the phase of the modulation as described below rather than by adjusting the phase of the A/D converter sampling.
- the FIFOs latch input data to the A/D converters 90 synchronously with the corresponding conversion clock 96 .
- the FIFO 91 output data is preferably provided to the internal components of the signal detector 10 synchronously with the synchronization clock 40 such that all further processing is synchronized with the synchronization clock 40 .
- the preferred implementation of the Temporal Response Analysis Engine 11 are shown in FIGS. 3 through 7 ; the preferred signal generator 1 is shown in FIGS. 3 and 4 , while the preferred signal detector 10 is shown in FIGS. 5 , 6 , and 7 .
- the Temporal Response Analysis Engine 11 is implemented as logic blocks within a Xilinx 4 FPGA.
- the functional blocks of the preferred signal generator 1 are shown in FIG. 3 .
- the top 41 and bottom 42 signal paths are two preferred variants for generating different code patterns for the modulation signal 16 .
- a Linear Feedback Shift Register (LFSR) 30 is preferably used to create a PRBS code.
- the specific code pattern is preferably determined by the number of state bits within the LFSR 30 and the gain code 36 input to the LFSR 30 .
- the gain code 36 is stored in a gain memory 31 , which is preferably configured to allow the code pattern 16 to be changed during operation either by selecting one of several gain codes from a read-only memory or by setting a new gain code into a writable memory.
- the gain code 36 may be hard-wired into the LFSR 30 , or a code-specific state-machine designed to generate a desired code through a series of state transformations may be used in place of the LFSR 30 .
- the entire code pattern is preferably stored as a bit sequence in a pattern memory 32 .
- the sequence in which pattern bits are presented is preferably determined by an address sequencer 33 which preferably provides the cell addresses 37 for the memory.
- the address sequencer 33 is preferably configured to allow changing the code pattern 16 during operation either by selecting one of several patterns stored in a read-only memory or by inputting a new pattern into a writable memory.
- the modulation signal 16 for both the LFSR 30 or pattern memory implementation is preferably buffered by an output buffer 35 to make the signals 16 more robust when driving external components.
- Timing for presentation of the code pattern bits is preferably controlled by a generation synchronizer 34 which preferably generates the master clock (MClk) 38 for the LFSR 30 and the address sequencer 33 .
- the master clock 38 is preferably synchronized to a system synchronization clock (SClk) 40 which preferably controls both code pattern generation and response signal acquisition.
- SClk system synchronization clock
- MClk 38 preferably operates at the same frequency as SClk 40 but is preferably offset in phase by an amount specified by the phase input 39 , which is preferably an externally programmable parameter.
- phase offset allows the relative phase between the modulation signal 16 and the detected response 19 to be adjusted. If the phase is adjusted by some increment, (360/K)°, at the end of each code pattern block or set of blocks the detected response resulting from the modulation signal will preferably be sampled at K discrete phases over K blocks.
- the detected response 19 is assumed to be stationary with respect to the start of the code pattern block over that time interval so that the K discrete sampling phases correspond to K discrete sample times and the effective temporal resolution of the sampling process is increased by a factor of K.
- This temporal over-sampling is functionally equivalent to the technique described for temporal over-sampling in the A/D converter embodiment.
- the external phase specification may represent the phase increment rather than the absolute phase, and the generation synchronizer 34 may increment the phase internally.
- the LFSR 30 is preferably a state-machine comprising M standard LFSR cells 48 which hold and transform the state.
- the LFSR cells 48 are preferably linked in a numbered sequence, and the output from the LFSR 30 is the current state of cell number zero.
- Each cell preferably comprises a state latch 45 which holds a single bit of state information, a gain element 46 to control the feedback gain for the cell based on the externally provided gain code 36 , and an accumulator 47 .
- the accumulator 47 preferably adds the feedback from the cell to the cumulative feedback from all previous cells.
- the state for a cell is updated to match the previous state from the next higher cell in the chain; the state of the last cell in the chain is updated with the accumulated feedback from all the previous cells.
- the accumulator 47 for the last cell in the chain may be omitted if desired.
- the pattern generated by the LFSR 30 is preferably determined by the number of cells in the chain and by the gain code.
- the gain code is provided from an external source to allow the code pattern to be modified.
- the gain code may be a fixed value.
- the implementation of the gain elements and accumulators for each cell may be optimized for the specific gain code for that cell rather than implemented in the generalized fashion shown.
- the clock for the LFSR 30 and for all its internal latches is preferably the signal generator master clock 38 .
- the preferred functional blocks for the signal detector 10 are shown in FIG. 5 .
- the detected response 19 and either the electronic reference signal 17 or the source reference signal 18 are received at two frame accumulators 50 and 51 , where the samples for each discrete sample time are accumulated with samples from identical sample times from different modulation frames to form the aggregated detected response 58 and the aggregated reference signal 59 .
- the effective data rate at which samples are preferably processed in following blocks is reduced by a factor equal to the number of frames aggregated into each sample point.
- the frame accumulators 50 and 51 are preferably replicated N times to handle the N channels of the A/D converter independently.
- the internal details of the frame accumulators 50 and 51 for the detected response and the reference signal may differ, depending on the digital format of the two signals.
- the reference signal used for analysis is the electronic reference signal 17 rather than the source reference signal 18 its value for each sample time is known a priori to be identical for every frame and to take on only two possible binary values, 0 or 1.
- the frame accumulator 51 for the reference signal 17 need only store one bit per sample time, equal to the value of the modulation signal for that sample time.
- the N acquisition/accumulation channels are preferably re-interleaved into a single data stream.
- two multiplexers 52 and 53 perform this reintegration at the output of the frame accumulators 50 and 51 . In other embodiments this re-integration may take place at any other desired point in the signal processing chain.
- the aggregated detected response 58 and the aggregated reference signal 59 are routed to the frame correlator 55 where the two signals 58 and 59 are preferably combined by a cross-correlation algorithm to produce the correlated signal 61 which preferably comprises a single value for each complete aggregated frame of samples.
- the correlated signal 61 represents the degree to which the aggregated response signal 58 contains components matching the aggregated reference signal 59 .
- the correlated signal 61 represents the degree to which the aggregated response signal 58 contains components of the delayed version of the aggregated reference signal 60 .
- the sample transfer characteristic 57 comprises a set of correlated signals calculated for a range of J such delay times.
- Phase delay blocks 54 generate the delayed versions of the aggregated reference signal 60 .
- the J phase delay blocks 54 are illustrated as discrete components operating in parallel and each providing the complete delay required for one correlated signal. In one preferred embodiment they comprise a cascade of J phase delay blocks each providing the time increment between one correlated signal and the next.
- the phase delays for the correlated signals are preferably discrete and correspond to integral multiples of the synchronization clock 40 period.
- the phase delay blocks 54 are preferably implemented as shift registers or FIFOs of the appropriate depth and clocked by the synchronization clock 40 .
- the time delay may be implemented using other methods.
- each phase delay is processed by a corresponding frame correlator 55 .
- a single frame correlator 55 may be used to calculate the correlated signal 61 for multiple phase delays by presenting the detected response data to its input multiple times, using a different phase delayed version of the reference signal 60 for each iteration. In this case fewer frame correlators 55 are required.
- the details of the preferred frame accumulator 50 or 51 are shown in FIG. 6 .
- Samples from the signal 17 , 18 , or 19 are preferably accumulated in the adder 70 by summing them with values taken from the memory 71 ; the resulting aggregated signal 58 or 59 is routed to the output of the accumulator and stored back into the memory at the same location from which the original data was taken.
- Each discrete sample time for the channel is represented by a single addressed cell within the memory.
- the size of the memory is preferably determined by two parameters, K and R, which preferably encode the sampling scheme.
- K represents the number of discrete phases at which samples are preferably taken in various frames during temporal over-sampling.
- R is the ratio of the number of samples in a modulation frame to the number of sampling channels provided in the A/D converter 90 for parallel over-sampling and signifies the number of samples that must be accommodated by each channel within a single frame.
- a preferred sample enable gate 72 is provided to restart the accumulation process at the beginning of each set of frames by clearing the cells in the memory.
- the address sequencer 73 selects the cell of the memory to be addressed for each sample point.
- the frame accumulators 50 or 51 preferably run synchronously with the synchronization clock 40 (although out of phase), so only a single address sequencer is required to address all the frame accumulators.
- the details of the preferred frame correlator 55 are shown in FIG. 7 .
- the ideal method for correlating the signals is to take the integral of the detected response 19 weighted by the reference signal 17 or 18 . Because the preferred embodiment is a sampled system the integration is approximated by summation over all the samples within a frame set using the adder 81 to generate the correlation signal 61 .
- the weighting of the aggregated detected response 58 by the aggregated reference signal 59 is preferably performed by a multiplier 80 . Other embodiments may employ other weighting and integration schemes, including scaling and integration in the analog domain directly on the detected signals.
- a sample enable gate 82 is preferably provided to restart the accumulation process at the beginning of each set of frames by clearing the correlator.
- FIG. 8 depicts a mechanical view of a preferred embodiment of the photon measurement system 100 .
- An electronics unit 85 preferably includes the modulated optical source 3 , optical detector 7 , temporal response analysis engine 11 , and associated electronics.
- the modulated optical wave 20 is delivered to the sample 5 using the delivery fiber bundle 86 .
- the scattered optical wave 21 is delivered from the sample to the electronic unit 85 using the detection fiber bundle 87 .
- the fiber bundles are flexible and easily repositioned with respect to the sample.
- FIG. 9 depicts a functional representation of another preferred embodiment in which a single optical fiber serves as part of both the light delivery optics and the detection optics.
- the light 91 output by the modulated optical source 90 is reflected by an optical filter 97 toward coupling optics 92 .
- the coupling optics focuses the light into an optical fiber 93 which delivers the light to an optical probe 94 with a removable needle 95 .
- the removable needle includes a section of optical fiber allowing the light to be transmitted from the probe to the end of the needle.
- the needle can be inserted into a tissue sample 96 so that greater depth can be interrogated than in the case of a surface-based probe. A portion of light scattered by the tissue is coupled back into the fiber and delivered back to the coupling optics 92 .
- a portion of this return light 98 is transmitted by the optical filter to an optical detector 99 .
- the modulated optical source 90 and the optical detector 99 are connected to the temporal response analysis engine and associated electronics as described previously.
- a mechanical view of the needle-based probe is given in FIG. 10 .
- the needle assembly 95 connects to the probe 94 with a standard SMA connector 101 .
- other geometric arrangements of light delivery optics, detection optics, and sample are possible.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Engineering & Computer Science (AREA)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Physics & Mathematics (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
A system and method for measuring photons utilizing a low-power light source modulated with a code sequence to interrogate a sample of interest. Preferably a portion of the scattered light from the sample is detected by a photo-detector. A correlation of the photo-detector signal and the code sequence produces an estimate of the distribution of flight times for photons traveling from the source to the detector.
Description
- This application is a continuation application claiming priority from the co-pending U.S. non-provisional patent application Ser. No. 11/381,450, entitled “Photon Measurement Method and Apparatus,” with filing date May 3, 2006, which is hereby incorporated by reference in its entirety.
- The field of the present invention pertains generally to systems and methods for detecting and measuring photons, including, more specifically, the measurement of the time-of-flight of photons traveling through a scattering media such as tissue.
- Diffuse optical imaging techniques are known in medical and biological applications. Overviews of diffuse optical imaging techniques can be found in “Recent Advances in Diffusion Optical Imaging” by Gibson, et al, Phys. Med. Biology, vol. 50 (2005), R1-R43 and in “Near-infrared Diffuse Optical Tomography,” by Hielscher, et al, Disease Markers, Vol. 18 (2002), 313-337. Briefly, diffuse optical imaging involves the use of near-infrared light incident upon a sample of interest. An example in the medical and biological field is optical mammography where near infrared light is used to illuminate breast tissue. A detector is placed on the opposite side of the breast from the incident light some distance away and collects scattered light from the breast tissue. The scattered light of interest that is detected may be directly scattered incident light or scattered fluorescence light caused by the excitation of an injected fluorescing material that fluoresces when exposed to the incident light. By measuring the amplitude of the light of interest at the detector and the distribution of photon arrival times at the detector for various source and detector positions, a reconstruction of the underlying tissue optical properties can be made. An overview of image reconstruction techniques can be found in the citations given in the aforementioned review articles.
- Measurements of the photon flight-time distributions are typically carried out using either the time-domain or the frequency-domain technique. In the time-domain technique, the sample is excited with a pulse of light from a pulsed laser and the scattered light is measured using a detector with single-photon sensitivity. The detector measures the time delay between the excitation pulse and the first detected photon. The flight-time distribution is determined by using many repeated pulses and building up a histogram of the measured time delays. Unfortunately, the pulsed laser sources and single-photon detectors are relatively expensive. Because detection is typically done at the single-photon level, it can require a significant amount of time to build-up enough data to approximate the flight-time distribution. One disadvantage of the frequency-domain approach is that it is not a direct measurement of the photon flight time. Rather, it provides an estimate of the mean flight time based on the phase shift between a detected signal and the excitation signal. In some cases, more accurate image reconstructions can be obtained using more complete measurements of the flight-time distributions. This data is not readily obtained with frequency-domain instrumentation. A further disadvantage of the frequency-domain approach is the need for accurate high-frequency analog electronics. An overview of both the time-domain and frequency-domain techniques can be found in the above-referenced article by Hielscher, et al.
- U.S. Pat. No. 5,565,982 discloses a time-resolved spectroscopy system using digital processing techniques and two low power, continuous wave light sources. The disclosed system requires two light transmitters of different wavelengths modulated with separate codes for interrogating a sample of interest. Properties of the sample are inferred by differential comparison of the return signals from each of the two light sources. It is undesirable to have two distinct light sources due to the cost and complexity involved. Furthermore, the noise level associated with a measurement made with two separate light sources will be higher than with a single source even if the codes used to drive the two sources are orthogonal.
- A system and method capable of addressing these disadvantages while providing acceptable temporal information for whatever application the photon measurement is being used is needed.
- The inventions presented herein provide for direct measurements of photon flight-time using any single light source modulated with a known digital pattern. A preferred system uses information from a single low-power continuous-wave light source detected by a low-cost detector. Preferably the measurement system is implemented with digital electronics. One embodiment of the system and methods disclosed comprises an information from a single continuous-wave light source modulated with a digital waveform for interrogating a sample, a photo-sensitive detector for measuring the scattered light from the sample, and electronics for sampling the detector output and performing a correlation of the output signal with the modulation waveform. Other embodiments include electronics and software for determining the parameters of the flight-time distribution from the measured correlation.
-
FIG. 1 is a functional block diagram of the major components of a preferred photon measurement system of the present invention. -
FIG. 2 is a diagram of preferred Analog-to-Digital converters and their interface to the signal detector. -
FIG. 3 is a functional block diagram of a preferred signal generator. -
FIG. 4 depicts an implementation of a preferred Linear Feedback Shift Register. -
FIG. 5 is a functional block diagram of a preferred signal detector. -
FIG. 6 is a functional block diagram of a preferred frame accumulator. -
FIG. 7 is a functional block diagram of a preferred frame correlator. -
FIG. 8 is a mechanical view of a preferred embodiment. -
FIG. 9 is a functional representation of a preferred embodiment using a needle-based optical probe coupled to a source and a detector with a single optical fiber. -
FIG. 10 depicts a preferred needle-based probe. - A functional block diagram of a preferred
photon measurement system 100 is depicted inFIG. 1 . The photon measurement system can be used to measure the interaction of photons with a sample 5. In certain applications, the sample 5 may be human breast tissue or fat tissue but it could just as well be any semitransparent material. Thephoton measurement system 100 preferably includes TemporalResponse Analysis Engine 11. The TemporalResponse Analysis Engine 11 generates a digital modulation signal for driving an illumination light source that is used to interrogate the sample. TheTemporal Response Engine 11 also provides a means for processing a detected optical signal from the sample 5 to extract information about the sample 5. Preferably adigital modulation signal 16 is generated in thesignal generator 1 and transmitted to thetransmit signal conditioner 2. Thedigital modulation signal 16 is the digital representation of a chosen code sequence. The code sequence is preferably chosen from the known pseudorandom binary sequences (PRBS), Gold codes, Golay codes, Kasami codes, Walsh codes, or other codes that possess the preferred desirable property of large auto-correlation values and low cross-correlation values. Thedigital modulation signal 16 may represent a single code pattern or multiple repeats of the same pattern. A single complete set of code patterns is designated a modulation frame or code pattern frame. Thedigital modulation signal 16 is preferably transmitted to thesignal detector 10 as anelectronic reference signal 17. The transmitsignal conditioner 2 formats thedigital modulation signal 16 as necessary to drive an optical illumination source 3. In the preferredphoton measurement system 100, the optical illumination source is a 785 nm continuous wave diode laser made by Hitachi Corp. Formatting of thedigital modulation signal 16 in the preferred embodiment involves converting thedigital modulation signal 16 to an analog voltage waveform that is coupled through a 50-ohm bias-T to the DC drive current of the optical illumination source 3. In other embodiments, the optical illumination source 3 may be a different laser diode, a light-emitting diode, or a light source used together with an external optical modulator. The optical illumination source 3 generates the modulatedoptical wave 20 which is preferably transmitted to the sample 5 bylight delivery optics 4. The preferredlight delivery optics 4 is a 3 mm diameter fiber bundle located between the optical illumination source 3 and the sample 5 to deliver the modulatedoptical wave 20 from the optical illumination source 3 to the sample 5. In other embodiments thelight delivery optics 4 comprises other arrangements of optical fibers, lenses, mirrors or other optical delivery components. When the modulatedoptical wave 20 illuminates the sample 5, scatteredoptical waves 21 are generated. In the preferred photon measurement system, the sample 5 is treated with a fluorescent material that will fluoresce when it is illuminated by the modulated optical waves 20. In the preferred system the scatteredoptical waves 21 are fluorescence generated from a fluorescent material within the sample 5. The fluorescent material is preferably an exogenous contrast agent injected into the sample 5 or alternatively it is preferably some constituent component of a material that exhibits endogenous fluorescence. The detection optics 6 are situated so that a portion of the modulatedoptical waves 21 are detected by the detection optics 6. In the preferredphoton measurement system 100, the detection optics 6 include an optical filter for separating the fluorescing scatteredoptical waves 21 from the modulated optical waves 20. The optical filter preferably transmits the higher wavelength fluorescence and blocks the lower wavelength illumination light. In applications where the scatteredoptical waves 21 of interest are not fluorescing, an optical filter is not required. - In the preferred
photon measurement system 100, the detection optics 6 preferably include a second 3 mm diameter fiber bundle located between the optical filter and theoptical detector 7. Theoptical detector 7 converts the scatteredoptical waves 21 to an electronic signal. In the preferredphoton measurement system 100, theoptical detector 7 is preferably a photomultiplier tube, model R7400U-20 from Hamamatsu Corp. In other embodiments, theoptical detector 7 may be a PIN photodiode, an avalanche photodiode, a charge-couple device, or other suitable photosensitive element. As previously stated, theoptical detector 7 preferably converts detected scatteredoptical waves 21 into an electronic signal which is communicated to the detected signal conditioner 8. The detected signal conditioner 8 preferably formats the signal so it may be converted to discrete samples by an Analog to Digital (A/D)converter 9. The A/D converter 9 outputs a detectedresponse signal 19. The detectedresponse signal 19 is communicated to asignal detector 10, where it is preferably correlated with theelectronic reference signal 17 to extract a sample transfer characteristic. - Information about the temporal properties of the photons is preferably calculated from the sample transfer characteristic. This information preferably includes such properties as direct measurements of photon time-of-flight and the fluorescence lifetime. The estimate of photon times-of-flight is then preferably used to estimate characteristics of the tissue such as the absorption coefficient, scattering coefficient, or location of fluorescing material.
- Another embodiment of the
photon measurement system 100 includes anoptical reference generator 22. Theoptical reference generator 22 preferably includes an 12A or 12B that routes a portion of the modulatedoptical splitter optical wave 20 to a secondaryoptical detector 13. The position of the 12A or 12B can be either before or after the light delivery optics. The output of the secondaryoptical splitter optical detector 13 is preferably routed to asecondary signal conditioner 14 whose output is communicated to a secondary A/D converter 15. The secondary A/D converter 15 preferably outputs asource reference signal 18 which can be correlated with the detectedresponse 19 to extract the sample transfer characteristic. Using thesource reference signal 18 as opposed to theelectronic reference signal 17 allows the filtering of the temporal properties of thesignal conditioner 2 and the modulated optical source 3 from the measured transfer characteristic. - The preferred hardware implementation of the A/D converter module and its interfaces to the
signal detector 10 are shown inFIG. 2 . An array of N A/D converters 90 preferably receives theanalog signal 95 in parallel from thesignal conditioner 8 or 14. The 18 or 19 from the A/output samples D converters 90 are preferably communicated to the First-In-First-Out buffers (FIFOs) 91 where they are buffered for distribution to the internal components of thesignal detector 10. In the preferred photon measurement system the A/D converters 90 are eight MAX 108 integrated circuits made by Maxim operating at 250 Msample/sec and outputting two data samples at a time in parallel at 125 MHz. The FIFOs 91 are preferably implemented within aXilinx 4 FPGA. Theacquisition synchronizer 92 preferably controls signal acquisition and digital data distribution through the conversion clock (CClk) signals 96. - The
acquisition synchronizer 92 is preferably synchronized with an externally provided synchronization clock (SClk) 40 which is also preferably used to synchronize thesignal generator 1. The signals CClk[1 . . . N] are preferably generated within theacquisition synchronizer 92 and preferably have the same frequency asSClk 40 but are offset in phase fromSClk 40 in N fixed increments of (360/N)°, with the phase of CClk[1] set to the fixed offset of Z°. In the preferred system the internal clock generation capabilities of the Xilinx FPGA are used to implement theacquisition synchronizer 92 directly. The A/D converters 90 preferably perform their conversions in sync with the conversion clocks 96 such that they generate samples at N discrete sample times spread evenly throughout the fundamental sample interval defined by the period ofSClk 40. The effective sample rate for the array of converters is preferably N times the rate defined bySClk 40. This process of using multiple A/D converters sampling out of phase to increase the effective sample rate is what we call parallel over-sampling. In the preferred photon measuring system, parallel over-sampling results in an effective sample rate of 2 Gsamples/sec. The offset value Z allows the entire sample set to be offset by some phase from thesynchronization clock 40. Theacquisition synchronizer 92 preferably is configured such that the value of Z can be varied synchronously with the modulation frame, or with a block of frames called a frame block. This allows Z to follow a sequence of K values smaller than (360/N)° such that on successive modulation frames/frame blocks the effective sampling phases (relative to the synchronization clock) take on K values intermediate to those created by the N conversion clocks in any given frame. In this case preferably the input signal at any given A/D converter 90 will be sampled at K discrete phases over K blocks. The detectedresponse 19 is preferably assumed to be stationary with respect to the start of the code pattern block over that time interval. The preferred K discrete sampling phases correspond to K discrete sample times and the effective temporal resolution of the sampling process is preferably increased by a factor of K. This process is referred to as temporal over-sampling. - In the preferred photon measuring system the value of Z is always zero and temporal over-sampling is achieved by adjusting the phase of the modulation as described below rather than by adjusting the phase of the A/D converter sampling. Preferably the FIFOs latch input data to the A/
D converters 90 synchronously with thecorresponding conversion clock 96. TheFIFO 91 output data is preferably provided to the internal components of thesignal detector 10 synchronously with thesynchronization clock 40 such that all further processing is synchronized with thesynchronization clock 40. - The preferred implementation of the Temporal
Response Analysis Engine 11 are shown inFIGS. 3 through 7 ; thepreferred signal generator 1 is shown inFIGS. 3 and 4 , while thepreferred signal detector 10 is shown inFIGS. 5 , 6, and 7. In the preferred system the TemporalResponse Analysis Engine 11 is implemented as logic blocks within aXilinx 4 FPGA. - The functional blocks of the
preferred signal generator 1 are shown inFIG. 3 . The top 41 and bottom 42 signal paths are two preferred variants for generating different code patterns for themodulation signal 16. In the top path 41 a Linear Feedback Shift Register (LFSR) 30 is preferably used to create a PRBS code. The specific code pattern is preferably determined by the number of state bits within theLFSR 30 and thegain code 36 input to theLFSR 30. In one preferred implementation thegain code 36 is stored in again memory 31, which is preferably configured to allow thecode pattern 16 to be changed during operation either by selecting one of several gain codes from a read-only memory or by setting a new gain code into a writable memory. In other embodiments thegain code 36 may be hard-wired into theLFSR 30, or a code-specific state-machine designed to generate a desired code through a series of state transformations may be used in place of theLFSR 30. In the bottom path 42 the entire code pattern is preferably stored as a bit sequence in a pattern memory 32. The sequence in which pattern bits are presented is preferably determined by an address sequencer 33 which preferably provides the cell addresses 37 for the memory. The address sequencer 33 is preferably configured to allow changing thecode pattern 16 during operation either by selecting one of several patterns stored in a read-only memory or by inputting a new pattern into a writable memory. - The
modulation signal 16 for both theLFSR 30 or pattern memory implementation is preferably buffered by anoutput buffer 35 to make thesignals 16 more robust when driving external components. Timing for presentation of the code pattern bits is preferably controlled by a generation synchronizer 34 which preferably generates the master clock (MClk) 38 for theLFSR 30 and the address sequencer 33. Themaster clock 38 is preferably synchronized to a system synchronization clock (SClk) 40 which preferably controls both code pattern generation and response signal acquisition.MClk 38 preferably operates at the same frequency asSClk 40 but is preferably offset in phase by an amount specified by thephase input 39, which is preferably an externally programmable parameter. This phase offset allows the relative phase between themodulation signal 16 and the detectedresponse 19 to be adjusted. If the phase is adjusted by some increment, (360/K)°, at the end of each code pattern block or set of blocks the detected response resulting from the modulation signal will preferably be sampled at K discrete phases over K blocks. In this embodiment of the photon measuring system as with the preferred embodiment, the detectedresponse 19 is assumed to be stationary with respect to the start of the code pattern block over that time interval so that the K discrete sampling phases correspond to K discrete sample times and the effective temporal resolution of the sampling process is increased by a factor of K. - This temporal over-sampling is functionally equivalent to the technique described for temporal over-sampling in the A/D converter embodiment. In other embodiments the external phase specification may represent the phase increment rather than the absolute phase, and the generation synchronizer 34 may increment the phase internally.
- The preferred implementation of the
LFSR 30 is shown inFIG. 4 . TheLFSR 30 is preferably a state-machine comprising Mstandard LFSR cells 48 which hold and transform the state. TheLFSR cells 48 are preferably linked in a numbered sequence, and the output from theLFSR 30 is the current state of cell number zero. Each cell preferably comprises astate latch 45 which holds a single bit of state information, again element 46 to control the feedback gain for the cell based on the externally providedgain code 36, and anaccumulator 47. Theaccumulator 47 preferably adds the feedback from the cell to the cumulative feedback from all previous cells. At each clock increment the state for a cell is updated to match the previous state from the next higher cell in the chain; the state of the last cell in the chain is updated with the accumulated feedback from all the previous cells. Theaccumulator 47 for the last cell in the chain may be omitted if desired. The pattern generated by theLFSR 30 is preferably determined by the number of cells in the chain and by the gain code. In a preferred embodiment the gain code is provided from an external source to allow the code pattern to be modified. In other embodiments the gain code may be a fixed value. In embodiments in which the gain code is fixed, the implementation of the gain elements and accumulators for each cell may be optimized for the specific gain code for that cell rather than implemented in the generalized fashion shown. The clock for theLFSR 30 and for all its internal latches is preferably the signalgenerator master clock 38. - The preferred functional blocks for the
signal detector 10 are shown inFIG. 5 . The detectedresponse 19 and either theelectronic reference signal 17 or thesource reference signal 18 are received at two 50 and 51, where the samples for each discrete sample time are accumulated with samples from identical sample times from different modulation frames to form the aggregated detectedframe accumulators response 58 and the aggregatedreference signal 59. As a result of this aggregation, the effective data rate at which samples are preferably processed in following blocks is reduced by a factor equal to the number of frames aggregated into each sample point. The frame accumulators 50 and 51 are preferably replicated N times to handle the N channels of the A/D converter independently. The internal details of the 50 and 51 for the detected response and the reference signal may differ, depending on the digital format of the two signals. For example, if the reference signal used for analysis is theframe accumulators electronic reference signal 17 rather than thesource reference signal 18 its value for each sample time is known a priori to be identical for every frame and to take on only two possible binary values, 0 or 1. In that case preferably theframe accumulator 51 for thereference signal 17 need only store one bit per sample time, equal to the value of the modulation signal for that sample time. At some point between the output of the frame accumulators and final output of the sample transfer characteristic 57 the N acquisition/accumulation channels are preferably re-interleaved into a single data stream. In one preferred embodiment two 52 and 53 perform this reintegration at the output of themultiplexers 50 and 51. In other embodiments this re-integration may take place at any other desired point in the signal processing chain. With or without re-integration the aggregated detectedframe accumulators response 58 and the aggregatedreference signal 59 are routed to theframe correlator 55 where the two 58 and 59 are preferably combined by a cross-correlation algorithm to produce the correlatedsignals signal 61 which preferably comprises a single value for each complete aggregated frame of samples. The correlatedsignal 61 represents the degree to which the aggregatedresponse signal 58 contains components matching the aggregatedreference signal 59. If the aggregatedreference signal 59 is delayed by a time t before presentation to thecorrelator 55 then the correlatedsignal 61 represents the degree to which the aggregatedresponse signal 58 contains components of the delayed version of the aggregatedreference signal 60. The sample transfer characteristic 57 comprises a set of correlated signals calculated for a range of J such delay times. Phase delay blocks 54 generate the delayed versions of the aggregatedreference signal 60. For simplicity the J phase delay blocks 54 are illustrated as discrete components operating in parallel and each providing the complete delay required for one correlated signal. In one preferred embodiment they comprise a cascade of J phase delay blocks each providing the time increment between one correlated signal and the next. The phase delays for the correlated signals are preferably discrete and correspond to integral multiples of thesynchronization clock 40 period. The phase delay blocks 54 are preferably implemented as shift registers or FIFOs of the appropriate depth and clocked by thesynchronization clock 40. In other embodiments the time delay may be implemented using other methods. In one preferred embodiment each phase delay is processed by acorresponding frame correlator 55. In other embodiments asingle frame correlator 55 may be used to calculate the correlatedsignal 61 for multiple phase delays by presenting the detected response data to its input multiple times, using a different phase delayed version of thereference signal 60 for each iteration. In this casefewer frame correlators 55 are required. - The details of the
50 or 51 are shown inpreferred frame accumulator FIG. 6 . Samples from the 17, 18, or 19 are preferably accumulated in the adder 70 by summing them with values taken from thesignal memory 71; the resulting aggregated 58 or 59 is routed to the output of the accumulator and stored back into the memory at the same location from which the original data was taken. Each discrete sample time for the channel is represented by a single addressed cell within the memory. The size of the memory is preferably determined by two parameters, K and R, which preferably encode the sampling scheme. K represents the number of discrete phases at which samples are preferably taken in various frames during temporal over-sampling. R is the ratio of the number of samples in a modulation frame to the number of sampling channels provided in the A/signal D converter 90 for parallel over-sampling and signifies the number of samples that must be accommodated by each channel within a single frame. A preferred sample enablegate 72 is provided to restart the accumulation process at the beginning of each set of frames by clearing the cells in the memory. Theaddress sequencer 73 selects the cell of the memory to be addressed for each sample point. The frame accumulators 50 or 51 preferably run synchronously with the synchronization clock 40 (although out of phase), so only a single address sequencer is required to address all the frame accumulators. - The details of the preferred
frame correlator 55 are shown inFIG. 7 . The ideal method for correlating the signals is to take the integral of the detectedresponse 19 weighted by the 17 or 18. Because the preferred embodiment is a sampled system the integration is approximated by summation over all the samples within a frame set using thereference signal adder 81 to generate thecorrelation signal 61. The weighting of the aggregated detectedresponse 58 by the aggregatedreference signal 59 is preferably performed by amultiplier 80. Other embodiments may employ other weighting and integration schemes, including scaling and integration in the analog domain directly on the detected signals. A sample enablegate 82 is preferably provided to restart the accumulation process at the beginning of each set of frames by clearing the correlator. - The geometric relationship between the
light delivery optics 4, the sample 5, and the detection optics 8 of the preferredphoton measurement system 100 depicted inFIG. 1 is schematic and not intended to limit the possible actual physical geometry of the system. For example, the delivery optics and the detection optics may be placed on the same side of the sample 5, on opposite sides of the sample 5, or at arbitrary positions with respect to the sample 5 so long as the scatteredoptical wave 21 is detectable.FIG. 8 depicts a mechanical view of a preferred embodiment of thephoton measurement system 100. Anelectronics unit 85 preferably includes the modulated optical source 3,optical detector 7, temporalresponse analysis engine 11, and associated electronics. Preferably the modulatedoptical wave 20 is delivered to the sample 5 using thedelivery fiber bundle 86. The scatteredoptical wave 21 is delivered from the sample to theelectronic unit 85 using thedetection fiber bundle 87. The fiber bundles are flexible and easily repositioned with respect to the sample. -
FIG. 9 depicts a functional representation of another preferred embodiment in which a single optical fiber serves as part of both the light delivery optics and the detection optics. The light 91 output by the modulatedoptical source 90 is reflected by anoptical filter 97 towardcoupling optics 92. The coupling optics focuses the light into anoptical fiber 93 which delivers the light to anoptical probe 94 with aremovable needle 95. The removable needle includes a section of optical fiber allowing the light to be transmitted from the probe to the end of the needle. The needle can be inserted into atissue sample 96 so that greater depth can be interrogated than in the case of a surface-based probe. A portion of light scattered by the tissue is coupled back into the fiber and delivered back to thecoupling optics 92. A portion of thisreturn light 98 is transmitted by the optical filter to anoptical detector 99. The modulatedoptical source 90 and theoptical detector 99 are connected to the temporal response analysis engine and associated electronics as described previously. A mechanical view of the needle-based probe is given inFIG. 10 . Theneedle assembly 95 connects to theprobe 94 with astandard SMA connector 101. In addition to the embodiments described here, other geometric arrangements of light delivery optics, detection optics, and sample are possible.
Claims (20)
1. A cost-effective and accurate photon measurement system for human tissue comprising:
a memory device configured to store a plurality of values corresponding to a plurality of pseudorandom binary code patterns wherein said plurality of values is located by a plurality of cell addresses;
a code selector configured to select a pseudorandom binary code pattern from said plurality of pseudorandom binary code patterns during operation of said photon measurement system;
an address sequencer coupled to said memory device and said code selector configured to provide a cell address from said plurality of cell addresses corresponding to said selected pseudorandom binary code pattern;
a pseudorandom binary code signal generator coupled to the code selector configured to generate a pseudorandom binary code output in accordance with said selected pseudorandom binary code pattern;
a continuous wave light source coupled to the pseudorandom binary code signal generator configured to generate a pseudorandom binary modulated optical wave modulated in accordance with said pseudorandom binary code output from the pseudorandom binary code signal generator;
delivery optics coupled to the continuous wave light source to direct the pseudorandom binary modulated optical wave at a human tissue translucent material;
an optical splitter coupled to the delivery optics between the delivery optics and the translucent material; the optical splitter coupled to an optical detector configured to direct the pseudorandom binary modulated optical wave to the optical detector; the optical detector configured to generate a pseudorandom binary optical reference signal; the optical detector coupled to a signal detector such that when the pseudorandom binary optical reference signal is generated it is communicated to the signal detector;
a corresponding detector configured to detect scattered optical waves from the translucent material; the corresponding detector having as an output a signal indicative of a characteristic of the detected scattered optical wave; the corresponding detector output signal communicated to the signal detector;
the signal detector configured to generate and provide an output indicative of a characteristic of the translucent material, based in part on the pseudorandom binary optical reference signal;
a correlator coupled to the signal detector configured to generate a correlation of said signal indicative of said characteristic of the detected scattered optical wave with said pseudorandom binary optical reference signal; and
a processor coupled to said correlator configured to calculate photon time-of-flight information of said scattered optical wave based on said correlation.
2. The system of claim 1 further comprising:
an acquisition synchronizer coupled to the corresponding detector to control acquisition of said signal indicative of said characteristic of the detected scattered optical wave; and
a synchronization clock coupled to said acquisition synchronizer to control phase of said acquisition synchronizer and wherein the optical splitter is positioned between the continuous wave light source and the delivery optics.
3. The system of claim 1 wherein a first analog to digital converter is positioned between the optical detector and the signal detector to convert the optical reference signal into a digital signal;
a second analog to digital converter is positioned between the corresponding detector and the signal detector to convert the corresponding detector output signal into a digital detector signal; and
the signal detector configured to use digital processing techniques to generate the output indicative of the characteristic of the translucent material, based in part on the digitized optical reference signal.
4. The system of claim 1 wherein the signal detector is configured to employ temporal over-sampling to generate and provide the output indicative of the characteristic of the translucent material.
5. The system of claim 1 wherein the pseudorandom binary code pattern represents a single code pattern.
6. The system of claim 1 wherein the pseudorandom binary code pattern represents multiple repeats of a same pattern.
7. The system of claim 1 wherein the pseudorandom binary code pattern is a Gold code pattern.
8. The system of claim 1 wherein the pseudorandom binary code pattern is a Kasami code pattern.
9. The system of claim 1 wherein the pseudorandom binary code pattern is a Walsh code pattern.
10. The system of claim 1 wherein said characteristic of the translucent material is an absorption coefficient.
11. The system of claim 1 wherein said characteristic of the translucent material is a scattering coefficient.
12. The system of claim 1 wherein said characteristic of the translucent material is an absorption coefficient and a scattering coefficient.
13. The system of claim 1 wherein said pseudorandom binary code signal generator comprises a linear feedback shift register.
14. The system of claim 13 wherein said pseudorandom binary code signal generator comprises a gain code input.
15. The system of claim 14 wherein said pseudorandom binary code signal generator comprises a gain code memory.
16. The system of claim 1 further comprising:
a linear feedback shift register coupled to said memory device.
17. The system of claim 16 further comprising:
a gain code input coupled to said linear feedback shift register.
18. The system of claim 17 further comprising:
a gain code memory coupled to said gain code input.
19. The system of claim 1 further comprising:
a state machine coupled to said memory device.
20. The system of claim 1 wherein said pseudorandom binary code signal generator comprises a state machine.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US13/598,324 US20120323109A1 (en) | 2006-05-03 | 2012-08-29 | Photon Measurement Method and Apparatus |
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US11/381,450 US8269964B2 (en) | 2006-05-03 | 2006-05-03 | Photon measurement method and apparatus |
| US13/598,324 US20120323109A1 (en) | 2006-05-03 | 2012-08-29 | Photon Measurement Method and Apparatus |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US11/381,450 Continuation US8269964B2 (en) | 2006-05-03 | 2006-05-03 | Photon measurement method and apparatus |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| US20120323109A1 true US20120323109A1 (en) | 2012-12-20 |
Family
ID=38662029
Family Applications (2)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US11/381,450 Expired - Fee Related US8269964B2 (en) | 2006-05-03 | 2006-05-03 | Photon measurement method and apparatus |
| US13/598,324 Abandoned US20120323109A1 (en) | 2006-05-03 | 2012-08-29 | Photon Measurement Method and Apparatus |
Family Applications Before (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| US11/381,450 Expired - Fee Related US8269964B2 (en) | 2006-05-03 | 2006-05-03 | Photon measurement method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| US (2) | US8269964B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107798378A (en) * | 2017-09-22 | 2018-03-13 | 中国科学院西安光学精密机械研究所 | Optimal code rapid generation and decoding method based on photon counting device |
| US20210374458A1 (en) * | 2020-05-26 | 2021-12-02 | Hi Llc | Systems and Methods for Data Representation in an Optical Measurement System |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7846391B2 (en) * | 2006-05-22 | 2010-12-07 | Lumencor, Inc. | Bioanalytical instrumentation using a light source subsystem |
| US10908084B2 (en) * | 2008-10-14 | 2021-02-02 | Timothy M. Ragan | Devices and methods for direct-sampling analog time-resolved detection |
| US8766161B2 (en) * | 2009-12-02 | 2014-07-01 | Nucript LLC | System for controling and calibrating single photon detection devices |
| US8527212B2 (en) | 2011-02-14 | 2013-09-03 | Honeywell Asca Inc. | Increased absorption-measurement accuracy through windowing of photon-transit times to account for scattering in continuous webs and powders |
| DE102011055330B4 (en) * | 2011-11-14 | 2024-10-31 | Leica Microsystems Cms Gmbh | Method for measuring the lifetime of an excited state in a sample |
| US9002112B2 (en) * | 2013-08-27 | 2015-04-07 | Trimble Navigation Limited | Video alignment system |
| CN110226943B (en) * | 2019-07-05 | 2023-08-15 | 上海联影医疗科技股份有限公司 | Parameter Calculation Method, Device and Computer Equipment for Photon Arrival at Detector |
| KR102747828B1 (en) * | 2020-03-02 | 2024-12-31 | 삼성전자주식회사 | Electronic device for measuring bio-signal using light and operation mehtod thereof |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5000568A (en) * | 1986-11-26 | 1991-03-19 | Hewlett-Packard Company | Spread spectrum optical time domain reflectometer |
| US5447159A (en) * | 1993-02-03 | 1995-09-05 | Massachusetts Institute Of Technology | Optical imaging for specimens having dispersive properties |
| US5565982A (en) * | 1994-05-31 | 1996-10-15 | Recon Exploration | Apparatus and method for time resolved spectroscopy |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5460182A (en) * | 1992-09-14 | 1995-10-24 | Sextant Medical Corporation | Tissue penetrating apparatus and methods |
-
2006
- 2006-05-03 US US11/381,450 patent/US8269964B2/en not_active Expired - Fee Related
-
2012
- 2012-08-29 US US13/598,324 patent/US20120323109A1/en not_active Abandoned
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5000568A (en) * | 1986-11-26 | 1991-03-19 | Hewlett-Packard Company | Spread spectrum optical time domain reflectometer |
| US5447159A (en) * | 1993-02-03 | 1995-09-05 | Massachusetts Institute Of Technology | Optical imaging for specimens having dispersive properties |
| US5565982A (en) * | 1994-05-31 | 1996-10-15 | Recon Exploration | Apparatus and method for time resolved spectroscopy |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107798378A (en) * | 2017-09-22 | 2018-03-13 | 中国科学院西安光学精密机械研究所 | Optimal code rapid generation and decoding method based on photon counting device |
| US20210374458A1 (en) * | 2020-05-26 | 2021-12-02 | Hi Llc | Systems and Methods for Data Representation in an Optical Measurement System |
| US11941857B2 (en) * | 2020-05-26 | 2024-03-26 | Hi Llc | Systems and methods for data representation in an optical measurement system |
Also Published As
| Publication number | Publication date |
|---|---|
| US20070260145A1 (en) | 2007-11-08 |
| US8269964B2 (en) | 2012-09-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20120323109A1 (en) | Photon Measurement Method and Apparatus | |
| US8433384B2 (en) | Method and apparatus for cerebral oximetry | |
| US20130338458A1 (en) | Method and apparatus for total hemoglobin measurement | |
| JP5275830B2 (en) | Optical ultrasonic tomographic imaging apparatus and optical ultrasonic tomographic imaging method | |
| EP1303756B1 (en) | Apparatus and method for probing light absorbing agents in biological tissues | |
| US7847950B2 (en) | Method and a system for generating three- or two-dimensional images | |
| US9408538B2 (en) | Method and apparatus for pressure sore detection | |
| EP2340429B1 (en) | Devices and methods for direct sampling analog time-resolved detection | |
| JPH1043166A (en) | Optimum diagnostic point search apparatus for blood-free diagnosis of blood components and blood-free diagnostic apparatus using the same | |
| EP0442295A2 (en) | Frequency domain fluorometry using coherent sampling | |
| EP0834277A1 (en) | Measuring absorption of a scattering medium | |
| EP0801727A1 (en) | Apparatus and method for time resolved spectroscopy | |
| WO2003050518A2 (en) | System and method for time correlated multi-photon counting measurements | |
| JP2006226698A (en) | Fluorescence detection device using intensity-modulated laser light | |
| US11788960B2 (en) | Device and method for detecting time-resolved optical data | |
| US20070259451A1 (en) | Fluorescence measurement method and apparatus | |
| US8755866B2 (en) | Method and apparatus for lymph node mapping | |
| US20060149150A1 (en) | Apparatus and method for acquiring time-resolved measurements utilizing direct digitization of the temporal point spread function of the detected light | |
| US8649010B2 (en) | Integral transformed optical measurement method and apparatus | |
| EP2021777A2 (en) | Fluorescence measurement and optical imaging method and apparatus | |
| JP2006266905A (en) | Chlorophyll analyzer and chlorophyll analysis method | |
| JP4167991B2 (en) | Fluorescence detection device using coded modulation laser | |
| JP5654509B2 (en) | Fluorescence detection apparatus and fluorescence detection method | |
| RU2370202C1 (en) | Automated device for diagnosing oncopathologies in gynecology | |
| CN114705661B (en) | A hyperspectral fluorescence lifetime imaging system and method based on streak camera |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| AS | Assignment |
Owner name: COVIDIEN LP, MASSACHUSETTS Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNOR:NELLCOR PURITAN BENNETT LLC;REEL/FRAME:029432/0260 Effective date: 20120929 |
|
| STCB | Information on status: application discontinuation |
Free format text: ABANDONED -- FAILURE TO RESPOND TO AN OFFICE ACTION |