US7728986B2 - Systems and methods for digital detection of a tomographic signal - Google Patents
Systems and methods for digital detection of a tomographic signal Download PDFInfo
- Publication number
- US7728986B2 US7728986B2 US12/257,825 US25782508A US7728986B2 US 7728986 B2 US7728986 B2 US 7728986B2 US 25782508 A US25782508 A US 25782508A US 7728986 B2 US7728986 B2 US 7728986B2
- Authority
- US
- United States
- Prior art keywords
- signal
- analog
- digital
- phase
- frequency
- 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.)
- Active, expires
Links
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/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
Definitions
- This invention relates in general to optical tomography, and in particular to systems, methods and apparatuses that perform digital detection for use in optical tomography.
- Optical tomography or optical tomographic imaging, is often performed by measuring the amplitude attenuation of light that has been passed through tissue. Propagation of injected photons is determined by the spatially-varying absorption and scattering characteristics of the tissue being probed. In biological tissues, scattering interactions are often the principal mechanisms affecting the light trajectory. As a result of this highly scattering nature, these photons do not navigate in a straight procession but rather diffuse throughout the medium. The photon flux exiting the tissue at any single point is the net effect of the incident light source, and discrete absorption and scattering interactions throughout their pathlength.
- the transmitted intensity measured along the target surface maintains the same frequency with respect to the source, however, it will exhibit an amplitude attenuation and an induced phase shift. This amplitude attenuation and phase shift provide spatial information regarding the absorption and scattering distribution inside the tissue.
- the light is either illuminated at a constant amplitude or modulated by a low frequency sine wave (up to a few kHz), and the decay in amplitude relative to the incident source is measured.
- a synchronous or homodyne detection technique is often employed to extract the zero-frequency amplitude information.
- This method requires generating a reference signal whose frequency is equal to and phase locked with the input waveform. Multiplying the input waveform by its reference signal produces an output waveform that is a composite of two independent contributions; one component located at zero-frequency and the other component straddling twice the modulation frequency. This resulting mixed signal is then sent through a low-pass filter to eliminate the higher frequency component, leaving only the remaining DC constituent whose amplitude is directly proportional to the amplitude of the detected optical signal.
- each wavelength must be modulated at distinct frequencies and/or phase in order to isolate the individual signals and their respective amplitudes.
- analog phase sensitive detection has many problems associated with it that adversely affect their performance and restrict subsequent applications.
- Some primary deficiencies include, for example, signal drift, output offsets, gain error, limited dynamic reserve, and harmonic rejection.
- external parameters such as temperature or age contribute to analog noise and, consequently, measurement uncertainty.
- analog processing is notably sensitive to component tolerances thereby limiting functional utility.
- coupling can occur, causing fluctuations along the analog lines.
- a direct consequence of these undesirable attributes is that the instrument noise floor is elevated, causing a reduction in the detection sensitivity, diminished dynamic range for the overall system, and a slowing of the data acquisition.
- Analog detection systems suffer from other deficiencies and limitations, as well.
- the invention provides systems, methods and apparatuses that perform digital detection for use in optical tomography.
- Methods and systems are provided in which digital lock-in detection is performed using an algorithm that employs a phase-independent quadrature technique; however, the invention also contemplates other phase-sensitive detection or lock-in techniques.
- a unique manipulation of an ordinary averaging filter optimized for sources discrimination and the consequent sampling constraints is presented as a novel filtering scheme for the lock-in detection; however, the invention also contemplates other filtering schemes.
- Systems and apparatuses are provided which include a digital signal processor which performs digital lock-in detection, data handling, and system management, and an integrated complex programmable logic device for timing control.
- Apparatuses are provided which include an instrument having an integrated digital signal processor for performing digital processing and detection for use in optical tomography.
- the invention provides a digital detection method for use in optical tomography, including obtaining an analog signal to contain information.
- the method further includes performing digital processing and detection on the digital signal to facilitate obtaining information from the digital signal.
- the invention provides a digital detection method for use in optical tomography that employs embedded signal processing by means of an integrated digital signal processor.
- the invention provides a digital data-acquisition method for use in optical tomography.
- the invention provides a lock-in detection method for use in optical tomography, including obtaining an analog signal modulated to contain information.
- the method further includes performing digital lock-in detection on the digital signal to facilitate obtaining the information from the digital signal.
- the invention provides a lock-in detection method for use in optical tomography, including obtaining an analog signal modulated to contain optical tomographic information.
- the method further includes performing digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- the invention provides a lock-in detection method for use in optical tomography, including obtaining an analog signal modulated to contain optical tomographic information.
- the method further includes amplifying the analog signal.
- the method further includes filtering the amplified analog signal using a low-pass filter.
- the method further includes converting the filtered analog signal into a digital signal.
- the method further includes performing digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- the invention provides a lock-in detection method for use in optical tomography, including obtaining an analog signal modulated to contain optical tomographic information.
- the method further includes amplifying the analog signal.
- the method further includes filtering the amplified analog signal using a low-pass filter.
- the method further includes converting the filtered analog signal into a digital signal.
- the method further includes performing digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- Performing digital lock-in detection includes using a lock-in algorithm employing a phase-independent quadrature technique and sampling constraints that optimize an averaging filter for source discrimination.
- the invention provides a lock-in detection system for use in optical tomography, including a low-pass anti-aliasing filter, an analog to digital converter, and a digital signal processor.
- An analog signal modulated to contain optical tomographic information is input to the low pass filter to obtain a filtered output analog signal.
- the filtered output analog signal is input to the analog to digital converter to be converted into an output digital signal.
- the output digital signal is input to the digital signal processor to perform digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- the invention provides a lock-in detection system for use in optical tomography, including a trans-impedance amplifier, a programmable gain amplifier, a low-pass anti-aliasing filter, an analog to digital converter, and a digital signal processor.
- An analog signal modulated to contain optical tomographic information is input to the trans-impedance amplifier for amplification to obtain a first output amplified analog signal.
- the first output amplified analog signal is input to the programmable gain amplifier for amplification to obtain a second output amplified analog signal.
- the second output amplified analog signal is input to the low pass filter to obtain a filtered output analog signal.
- the filtered output analog signal is input to the analog to digital converter to be converted into an output digital signal.
- the output digital signal is input to the digital signal processor to perform digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- the invention provides an apparatus for use in optical tomography, including a low-pass anti-aliasing filter, an analog to digital converter; and a digital signal processor.
- An analog signal modulated to contain optical tomographic information is input to the low pass anti-aliasing filter to obtain a filtered output analog signal.
- the filtered output analog signal is input to the analog to digital converter to be converted into an output digital signal.
- the output digital signal is input to the digital signal processor to perform digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- the invention provides an apparatus for use in optical tomography, including a trans-impedance amplifier, a programmable gain amplifier, a low-pass anti-aliasing filter, an analog to digital converter; and a digital signal processor.
- An analog signal modulated to contain optical tomographic information is input to the trans-impedance amplifier for amplification to obtain a first output amplified analog signal.
- the first output amplified analog signal is input to the programmable gain amplifier for amplification to obtain an second output amplified analog signal.
- the second output amplified analog signal is input to the low pass filter to obtain a filtered output analog signal.
- the filtered output analog signal is input to the analog to digital converter to be converted into an output digital signal.
- the output digital signal is input to the digital signal processor to perform digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- the invention provides an apparatus for use in optical tomography, including a low-pass anti-aliasing filter for receiving an analog signal modulated to contain optical tomographic information, and for filtering the analog signal to generate a filtered analog signal.
- the apparatus further includes an analog to digital converter, communicatively connected to the low-pass anti-aliasing filter, for receiving the filtered analog signal, and for converting the filtered analog signal into a digital signal.
- the apparatus further includes a digital signal processor, communicatively connected to the analog to digital converter, for receiving the digital signal, and for performing digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- the invention provides an apparatus for use in optical tomography, including a trans-impedance amplifier for receiving an analog signal modulated to contain optical tomographic information, and for amplifying the analog signal to generate a first amplified analog signal.
- the apparatus further includes a programmable gain amplifier, communicatively connected to the trans-impedance amplifier, for receiving the first amplified analog signal, and for amplifying the first amplified analog signal to generate a second amplified analog signal.
- the apparatus further includes a low-pass anti-aliasing filter, communicatively connected to the programmable gain amplifier, for receiving the second amplified analog signal, and for filtering the second amplified analog signal to generate a filtered analog signal.
- the apparatus further includes an analog to digital converter, communicatively connected to the low-pass anti-aliasing filter, for receiving the filtered analog signal, and for converting the filtered analog signal into a digital signal.
- the apparatus further includes a digital signal processor, communicatively connected to the analog to digital converter, for receiving the digital signal, and for performing digital lock-in detection on the digital signal to facilitate obtaining the optical tomographic information from the digital signal.
- FIG. 1 is a block diagram depicting a prior art analog detection system used in optical tomography
- FIG. 2 is a block diagram of an optical tomography system that uses digital detection, according to one embodiment of the invention
- FIG. 3 is a block diagram of an instrument layout for an optical tomography system that uses digital detection, according to one embodiment of the invention
- FIG. 4 is a graph depicting filter frequency response of an anti-aliasing filter, according to one embodiment of the invention.
- FIG. 5 is a flow diagram depicting a digital signal processor-based Imaging routine, according to one embodiment of the invention.
- FIGS. 6A and 6B are a table depicting a timing scheme for an imaging routine, according to one embodiment of the invention.
- FIG. 7 is a flow diagram depicting a lock in algorithm, according to one embodiment of the invention.
- FIG. 8 is a set of graphs depicting lock-in filter wavelength discrimination, according to one embodiment of the invention.
- FIG. 9 is a block diagram depicting hardware components of a digital CW instrument, according to one embodiment of the invention.
- FIG. 10 is a set of graphs depicting frequency response of filter, original, and modulated signals, according to one embodiment of the invention.
- FIG. 11 is a set of graphs that depict plots of calculated signals along with the measured signals, according to one embodiment of the invention.
- FIG. 12 is a block diagram depicting system layout and signal communication, according to one embodiment of the invention.
- FIG. 13 is a graph 1300 depicting frequency response of a 20 point averaging filter, with the horizontal frequency axis running from zero to Nyquist frequency;
- FIG. 14 is a flow diagram depicting a method 1400 according to one embodiment of the invention.
- FIG. 15 is a block diagram depicting an optical switch, according to one embodiment of the invention.
- Some embodiments of the invention provide, for example, digital processing and acquisition methods which are better than existing analog detection systems.
- FIG. 1 (prior art) and FIG. 2 are referenced in comparing the conventional analog detection scheme with a digital detection scheme according to embodiments of the invention.
- FIG. 1 is a block diagram depicting a prior art analog detection system 100 used in optical tomography. Even though the analog layout shown in FIG. 1 is typical for many measurement devices, a focus will be on a direct comparison between systems according to embodiments of the invention and instrument described in Schmitz 2002.
- the instrument described in Schmitz 2002 is essentially a first generation optical tomographic imager, and reference to it will aid in understanding some embodiments of the present invention, or aspects thereof.
- a transducer is generally required to convert modified parameters into a detectable measurement.
- the modified parameter is photon flux.
- a light source may be directed at a specific tissue or limb of interest and the attenuated light exiting the tissue is detected at multiple positions along the surface. These measurements are then fed through an image reconstruction algorithm whose results impart insight into anatomic structure or physiologic conditions. Imaging with multiple wavelengths allows one to perform spectral analysis and quantify changes in hemoglobin parameters such as oxy and deoxy-hemoglobin, or other biological chromophores.
- the light sources utilized with the instrument described in Schmitz 2002 can be, for example, laser diodes emitting light at wavelengths in the 700 nm-850 nm range.
- the laser diodes are amplitude-modulated at specific frequencies. Modulating the light source permits the use of lock-in (homodyne) detection which is a technique used to extricate signals of specific frequencies buried in noise of potentially greater magnitude.
- lock-in homogene
- a photodiode is used to convert the light energy into a current that can be measured.
- a gain stage is required to bring the incoming signal to a detectable level, i.e. above system noise floor. Additionally, the amplification serves to increase the dynamic range of the instrument and maximize the sensitivity. The appropriate gain for a given channel is dependent on the target geometry, tissue density, and the detector's proximity to the light source and must be updated as the source position changes.
- a variable trans-impedance stage followed by a programmable amplifier offer a multitude of possible signal gains. Once the signal is adequately amplified, the DC value or amplitude can be extracted from it.
- the incoming signal After the incoming signal has been amplified appropriately, it goes through an analog lock-in filter which consists of several stages. Initially, the signal is mixed with a ‘reference’ signal of identical frequency to the one being extracted (one reference for each wavelength). The affect of this mixing procedure is that the resulting signal now has two components. The first component is sitting at DC. The second component is a sinusoid at twice the modulation/reference frequency. In the second stage, the mixed signal goes through a low-pass filter to remove the higher frequency component. All that remains is the DC component whose magnitude is directly proportional to the magnitude of the input signal and the phase difference between the input signal and reference.
- the Instrument described in Schmitz 2002 employs a phase-shifter to minimize this phase difference and maximize the output. This procedure is executed for each wavelength.
- another lock-in method which includes a quadrature-based detection technique can be used that would eliminate the need for this phase-shifter. To fully appreciate these techniques, a mathematical synopsis is provided herein.
- the magnitude information in some embodiments, it must be read into a computer and stored in a file. This is accomplished via a commercially available data-acquisition board. These boards have numerous analog input channels. For example, the instrument described in Schmitz 2002 uses one with 64, (32 channels ⁇ 2 wavelengths).
- the DAQ board utilizes on-board multiplexers to sequence the sampling events for all inputs. Before it is digitized by an analog-to-digital converter (ADC), it must pass through an anti-aliasing filter to ensure that the signal conforms to the Nyquist constraints. The digitized samples are temporarily stored in a buffer and streamlined into the memory space of the host computer. This process is controlled by the DAQ software, LabVIEW, available from National Instruments Corp.
- FIG. 2 is a block diagram of an optical tomography system 200 (a system used in or for optical tomography, or aspects thereof) that uses digital detection, according to one embodiment of the invention.
- the modulated signal is to be digitized immediately following programmable gain segment 202 .
- an anti-aliasing filter must be included, and is depicted as included in segment 204 .
- ADC 204 samples the analog channels two at a time. Data samples get clocked out of the ADC and into the DSP 208 via a serial port connection. Once in the DSP 208 , it goes through a digital lock-in filter which also consists of a mixing process and low-pass filter. Upon completion of all computations, the DSP 208 sends out the data to the host PC 210 for storage, display, and image reconstruction.
- Some embodiments of systems according to the invention include components or design features that are replaced or re-designed with different components, relative to components or design features of the instrument and system described in Schmitz 2002, or relative to components or features which may be used in a typical analog system 100 as depicted in FIG. 1 (Prior Art).
- This noise consists of two components: one is the intrinsic noise generated by the electronics themselves and the other is extrinsic noise, radiated from some outside source.
- Another change is that the settling time of each channel after a change in source positions is decreased thereby increasing the temporal resolution. This settling time is the limiting factor determining overall image frame rate.
- the input channels of the instrument described in Schmitz 2002 settle in about 10 ms due to the source switch and the hardware lock-in amplifier.
- An instrument according to some embodiments of the invention has one analog detection channel per circuit board using through-hole components. By using surface-mount components, four analog channels per board are to be used, decreasing the number of total boards.
- Another change is that the input dynamic range of the system is increased. It can be difficult to saturate the detector channels when the optical signal is large.
- FIG. 3 is a block diagram of an instrument layout for an optical tomography system 300 that uses digital detection, according to one embodiment of the invention.
- ADC analog-to-digital converter
- Each ADC has an internal 4-channel MUX, so that four detector channels can be included on each detector board.
- Each detector channel will consist of a trans-impedance amplifier (TIA), a programmable-gain amplifier (PGA), a low-pass filter (LPF), and the ADC 302 .
- TIA trans-impedance amplifier
- PGA programmable-gain amplifier
- LPF low-pass filter
- the detector channel 304 includes TIA 306 , PGA 308 , LPF 310 , and the ADC 302 .
- FIG. 4 is a graph 400 depicting filter frequency response of a low pass anti-aliasing filter according to one embodiment of the invention.
- analog signals were passed through a backplane and eventually to an ADC card on a host PC.
- ADC Analog Detection Boards
- Another system change is to move the lock-in amplification (LIA) function from hardware to software.
- LIA lock-in amplification
- One benefit is a decrease in the settling time of the detector electronics.
- Preliminary LIA code shows settling times in the 1-2 ms range as opposed to 7 ms in the hardware LIA of the instrument described in Schmitz 2002.
- the rest of the electronics in a digital detector board take from 2-4 ms to settle (determined by the high-pass filter at the TIA output—this can probably settle simultaneously with the source switch).
- the overall settling time of the system's detector electronics is between 3 ms and 6 ms (easily adjustable within this range by minor differences in the design).
- Another benefit to moving the LIA to software is a decrease in the complexity of the ADB.
- significant performance benefits can be realized by implementing the LIA in the digital domain. This is a result of eliminating the adverse analog effects of the analog electronics, analog signal processing, and filter chips—and being able to customize the filter for user-specific applications.
- a third gain setting for the TIA is included. This 1-k ⁇ gain setting will be included with the previous 10-k ⁇ and 10-M ⁇ gain settings.
- TIA TIA according to some embodiments of the invention, such as the TIA 306 depicted in FIG. 3 , as referenced above.
- Rf trans-impedance feedback resistor
- TIA Rf there are three values of TIA Rf to choose from (1 k ⁇ , 10 k ⁇ and 10 M ⁇ ).
- the PGA 308 of system 300 only has three gain settings (1, 10 and 100). This allows for overall gains ranging from 1 kA/V to 1 G ⁇ A/V spaced by a factor of 10. Using the 1 k ⁇ resistor results in a 10 times larger input dynamic range at the lowest gain setting.
- One motivation for this additional gain flexibility is the realization that when performing measurements on small geometries (1 cm-3 cm) as is common with small animal imaging, the TIA stage will readily saturate for small optode separations and intensity filters must be introduced.
- the PGA gain setting of 1000 has been omitted in the new system because the noise levels at this gain setting become too high to achieve 1% linearity in many cases.
- Another change or improvement, in systems according to embodiments of the invention, is to extend the bandwidth of the TIA 306 at its highest gain setting (10 M ⁇ ).
- the BW was limited to 16-kHz by a 1-pF feedback capacitor in parallel with the 10 M ⁇ resistor. This capacitance cannot be reduced much because parasitic board capacitances can be between 0.5 pF and 3 pF anyways.
- reed relays K 1 and K 2
- R 1 , R 2 and R 3 the three feedback resistors
- K 1 and K 2 the reed relays are controlled by an 8-1 analog MUX and three gain bits.
- the corresponding reed relay coil is grounded through the MUX. This allows the relay coil to conduct, closing the switch.
- K 1 or K 2 is activated, R 1 or R 2 is placed in parallel with R 3 .
- R 1 and R 2 are both at least 1000 times lower than R 3 , it follows that R 1 or R 2 sets the gain of the TIA when K 1 or K 2 is activated.
- only one reed relay switch can be closed at a time (because only one channel of the 8-1 MUX can be selected at a time).
- the reed relays consume more current than any other circuit element, it is significant that only one can be activated at a given time.
- R 3 sets the TIA gain.
- the TIA circuit consists of an op-amp, two reed relays, a MUX, and various passive components. The following includes an explanation of component selections.
- a photodiode such as that used and described in Schmitz 2002 is utilized.
- the op-amp (IC 1 ) utilized has a very low input bias current. This is important when trying to accurately measure photodiode currents in the tens of picoamps range.
- the op-amp also has low voltage noise and low current noise. With a 10 M ⁇ feedback resistor, this current noise generates a voltage noise of 7 uV. This is less than the resistive thermal noise of the 10 M ⁇ resistor. It is believed that the current noise increases with temperature so it is important to keep the temperature not much higher than 25° C.
- PGA gains are to be had of 1, 10 and 100, instead of the PGA gains of the instrument described in Schmitz 2002 of 1, 10, 100 and 1000.
- an instrumentation amplifier gain programmed by a resistor
- 8-1 MUX 8-1 MUX
- various passive components The integrated PGA as used in the instrument described in Schmitz 2002 is not used, in some embodiments of the invention, for at least a few reasons, as follows. First, the integrated PGAs tend to have fairly high supply currents. Second, many PGAs have relatively low input impedances. Low input impedance is not good for the design because the input to the PGA is a 4.7 k ⁇ resistor and voltage division would result. Also, the PGA as included in the system described in Schmitz 2002 is relatively non-linear (0.015%, while 16-bit accuracy is about 0.0015%).
- the 8-1 MUX connects a different resistor between pins 1 and 8 of the instrumentation amplifier based on two control bits.
- the particular in-amp used for the PGA 308 is available from Texas Instruments, Inc., although other in-amps may be utilized in other embodiments of the invention.
- This amplifier has low quiescent current, low input bias current, low input offset voltage, low non-linearity, dual supply operation and a high gain-bandwidth product.
- the maximum quiescent current of 0.75 mA is almost a tenfold decrease in current consumption compared to the PGA 308 .
- Analog Devices Model ADG608B MUX available from Analog Devices, Inc., can be used as the analog MUX because it operates on ⁇ 5-V supplies and has an acceptable on-resistance (measured to be about 16 ⁇ ). Operating on ⁇ 5-V supplies is important because the input can be positive or negative.
- Resistors are chosen for PGA gains of 1, 10 and 100 according to Equation (1):
- the 8-1 MUX governed by two control bits—connects a different resistor R G across the appropriate pins of the instrumentation amplifier.
- the following chart (Table 1) tabulates the total effective signal gain based on the values of all three control bits.
- the appropriate gain values are determined prior to making a measurement and depend on the target geometry, tissue density, and the detector's proximity to the light source. These values for all detector channels must be updated for each source position.
- the corresponding bits are stored in local memory banks located on each detector card and made available for immediate access.
- Including the on-board memory cache into the modular structure of the detector cards has numerous advantages. First, it eliminates the need for a designated gainbit-routing board. Because of the substantial fan-out generated by such a routing board, it necessitates a sizable allotment of space and produces many parallel lines running throughout the instrument. Secondly, this modular format makes it very simple to add additional detector channels without requiring a system restructure. Furthermore, it is the most compact and efficient way of distributing the gainbit values.
- the following description relates to aspects and specifications of an anti-aliasing LPF according to embodiments of the invention, such as the LPF 310 depicted in FIG. 3 .
- the system 300 is designed to provide a signal passband from 3 kHz to 9 kHz (with 1% flatness). Since the TIA 306 sets the attenuation at 9 kHz to 1%, the attenuation at 9 kHz caused by the LPF stage is to be minimized. In some embodiments, the system 300 is designed to provide a 0.1% LPF attenuation at 9-kHz and as steep of a roll-off as possible after 9 kHz while preserving flatness in the passband. This roll-off attenuates any noise above 9 kHz, aiding the LPF function implemented in the LIA. It also provides anti-aliasing before sampling of the signal with the ADC.
- the system 300 is designed to provide at least this much attenuation at any frequencies that could alias to 9 kHz or below.
- the LPF 310 is designed to have a good time response to a step input (for when the source position and gain settings switch) and to settle to 16-bit accuracy in less than 1 ms with a full-scale step input. Also, the LPF 310 IC is designed to provide low-power and small size.
- the LPF 310 is implemented with an 8 th order Butterworth filter.
- the Butterworth filter topology tends to have the flattest passband and relatively steep roll-off. They also tend to have a fairly good time response to step inputs.
- An 8 th order filter is used because it is practical with the currently available filter ICs. Higher orders would require more filter ICs in series (increasing size and power consumption), while lower orders do not provide the necessary attenuation at frequencies just above 9 kHz.
- a cutoff frequency of 12.5 kHz is used to meet the specification of 0.1% LPF attenuation at 9 kHz.
- the system 300 utilizes two 4th order filter ICs, cascaded to implement the 8 th order Butterworth LPF.
- These filter ICs have low noise and distortion (specified for 16-bit systems). They are continuous-time active filters so they do not suffer from clock noise as switched-capacitor filters do.
- these 4 th order filter Ics are used instead of a single 8 th order IC because these have a low-power mode with a quiescent current of 2.2 mA each. This is considerably less current consumption than the lowest single-IC 8 th order continuous-time filters.
- the system 300 incorporates a 2.5V offset into the input buffer of the ADC 302 .
- Signals exiting the LPF 310 are centered about zero volts whose amplitude has an upper limit of 5V p-p .
- the ADC though, operates between 0V and +5V so we must add a 2.5V offset to form a pseudo-bipolar signal. This preserves the 5V p-p now centered around 2.5 volts.
- the driver must be able to settle for a full-scale step of the capacitor array at a 16-bit level (0.0015%).
- the noise generated by the driver must be kept to a minimum in order to preserve the SNR performance of the ADC.
- the following description relates to aspects and specifications of a DSP according to embodiments of the invention, such as the DSP 208 depicted in FIG. 2 , as referenced above.
- the system 300 is designed to provide the highest possible resolution ADC while retaining the ability to sample at ⁇ 50 ksps for each of four multiplexed channels.
- the ADC 302 would preferably have at least a four-channel input multiplexer built-in. It should also have a serial output to the DSP 208 . This will result in less digital data lines than a parallel data output. Low current consumption would also be helpful, but is not critical since there is only one ADC 302 per four detector channels.
- a 4-channel, 16-bit, 1 Msps ADC is used. This is the fastest 4-channel, 16-bit ADC available at present. The effective number of bits of this ADC is about 14. This is sufficient because the noise floor of the analog electronics does not permit better than 14-bit resolution except for two gain settings (noise is dominated in most gain settings by the TIA feedback resistor).
- the 1 Msps sampling rate allows sampling of each channel at a sufficiently high rate (sampling is being done at ⁇ 75-80 ksps). It also has a relatively low current consumption.
- the ADC's sample two channels simultaneously. Initially it samples the first two channels, whose samples are read into the DSP 208 over the serial port line. The MUX is then switched and it samples and clocks out the remaining two channels.
- the following includes a description of aspects and functions of a DSP and a complex programmable logic device (CPLD), according to some embodiments of the invention. It is through these devices where the multitude of tasks takes place for each source position and each imaging frame. These two devices require meticulously-orchestrated timed operations and must be in continuous communication with one another.
- CPLD complex programmable logic device
- the Digital Signal Processor plays a substantial role in data acquisition timing and overall system control, providing or helping to provide a number of advantages to an instrument into which it is integrated, which can include computational power, flexibility, and speed.
- the DSP upon power up, goes through a boot-up sequence which loads the programming code into internal memory.
- This sequence is a set of instruction used to initialize and direct the loading of the DSP.
- This boot-sequence is stored in an on-board EEPROM and the DSP is programmed to download and read its contents when power is applied.
- the DSP goes through a series of system initialization processes so that registers are appropriately prepared for running experiments. These processes include:
- the DSP is continuously polling the digital I/O fines connected to the PC waiting for an instruction.
- Gain-bit values are received from the PC, sent over the data bus through a relevant transfer protocol, and stored in DSP onboard memory.
- DSP sends instruction word to CPLD and tells it to begin going through gain-upload sequence.
- the gainbits for all channels are sent out from the DSP one source at a time and are routed to the detector boards by the CPLD.
- the DSP increments through the number of sources and then stops when all bits are uploaded.
- the gainbits are sent out over the serial port transmit line using a DMA protocol which has an independent processor so the DSP's core is not occupied with I/O responsibilities.
- CPLD controls sampling process and DSP waits for the data to come in over the serial port. Data coming into the DSP is time-division-multiplexed. The DMA channel is setup such that it collates the samples for each channel so that when sampling is complete, data is organized in the DSP's memory block.
- the DSP mixes the reference signals and the data samples. It filters the resulting signal with an averaging filter (more detail on that below) to discard the higher frequency components.
- Steps b-f are repeated for all sources until the imaging frame is complete. After data from last source is received, the DSP check the digital I/O status line to see if it should image another frame.
- control lines indicate termination of data acquisition, the DSP completes the frame it is processing and then stops.
- the CPLD does the same.
- the CPLD acts as a multi-level state machine and can execute a variety of standard logic operations. Once enabled, the CPLD's main state machine enters standby mode. There it waits for an instruction from the DSP telling it to enter either the Gainbit or Imaging state-machines.
- the CPLD after entering a Gainbit state machine, performs the following:
- the CPLD after entering a Gainbit state machine, performs the following:
- FIG. 5 is a flow diagram depicting a DSP-based Imaging routine 500 , according to one embodiment of the invention.
- step 502 it is inquired whether an imaging request has been issued.
- routine 500 proceeds to steps 504 , 506 , and 508 . These steps include, respectively, activate imaging protocol in the CPLD, initialize timer for settling time of electronics, and setup SPORT (Serial Port) receive DMA for the detection data.
- steps 504 , 506 , and 508 include, respectively, activate imaging protocol in the CPLD, initialize timer for settling time of electronics, and setup SPORT (Serial Port) receive DMA for the detection data.
- step 506 the routine 500 proceeds to step 510 , at which it is queried whether a timer interrupt has been issued. If not, then the routine 500 remains in query step 510 .
- routine 500 proceeds to steps 512 and 514 , at which include, respectively, disable timer and collate incoming time-division multiplex samples.
- step 514 the routine 500 proceeds to step 516 , at which it is queried whether a SPORT interrupt has been issued. If not, then the routine remains at step 514 .
- step 520 which includes, mix optical signals with references for each wavelength.
- step 520 the routine 500 proceeds to step 522 , which includes utilization of the LPF, and also returns to step 518 .
- step 518 it is inquired whether the source is the last source. If yes, then the routine 500 returns to step 502 . If no, then the routine 500 returns to step 506 .
- step 522 the low pass filter is utilized. Following step 522 , the routine 500 proceeds to step 524 , which includes sending & storing measurements in external memory.
- step 526 includes providing additional control for streamlining data to a computer or computerized device, such as a personal computer (PC), or, in some embodiments, over one or more networks to one or more computing devices.
- a computer or computerized device such as a personal computer (PC)
- PC personal computer
- the PC could be any computing device or devices, and could be reachable through or include one or more networks.
- FIGS. 6A and 6B are a table depicting a timing scheme 600 for an imaging routine, according to one embodiment of the invention.
- the vertical axis 602 indicates time passing as the axis 602 proceeds downwardly, while each task 606 is indicated by an entry in the table 602 in a vertical position indicating the time or time period at which the task is performed or accomplished.
- FIG. 7 is a flow diagram depicting one type of lock in algorithm 700 , according to one embodiment of the invention.
- the algorithm 700 extracts from the measured signal M[k] the amplitude A ⁇ i of the component due to wavelength ⁇ i.
- the algorithm stores K values for cosine (Inphase) and sine (Quadrature) wave at frequencies f ⁇ i .
- M[k] is modulated by the stored signals forming quadrature components and then passed through an averaging filter in accordance with the following equations:
- FIG. 8 is a set of graphs 802 , 804 , 806 depicting lock-in filter wavelength discrimination, according to one embodiment of the invention.
- the lock in filter can be, but is not limited to, an ordinary averaging filter.
- the averaging or mean filter is presented and described in detail such that its unique exploitation, according to some embodiments of the invention, is demonstrated.
- other digital filters including but not limited to, cascaded RC filters, Bessel filters, critically damped filters etc., are utilized in some embodiments of the invention.
- Digital low-pass filter characteristics include, but are not limited to, Bessel and critically damped (CD) filters. Both are characterized by a good time-domain response, i.e. fast settling and little overshoot and ringing.
- a multi-order CD filter is used in some embodiments of the invention because it has superior step response characteristics (i.e. fast settling without overshoot) compared to other filter types.
- the following provides a brief mathematical synopsis of the lock-in algorithm 700 , as referenced above, and is followed by a unique exploitation of a simple finite impulse response (FIR) low pass filter, according to one embodiment of the invention.
- FIR finite impulse response
- phase sensitive detection One stipulation imposed on by phase sensitive detection is that the source intensity must be modulated with a sinusoidal waveform. In our instance, this is accomplished by modulating the bias current feeding the laser diode.
- V mj (t) A mj ⁇ sin(2 ⁇ f m t+ ⁇ mj ), where A mj is the amplitude of the optical signal at detector j, f m is the source (modulation) frequency, and ⁇ mj is the phase shift.
- V mj ⁇ S ref A mj ⁇ A r ⁇ ( 1 / 2 ) ⁇ ⁇ cos ⁇ ( [ 2 ⁇ ⁇ ⁇ ⁇ f m - 2 ⁇ ⁇ ⁇ ⁇ f r ] ⁇ ( n / f s ) + ⁇ mj - ⁇ ref ) - cos ⁇ ( [ 2 ⁇ ⁇ ⁇ ⁇ f m + 2 ⁇ ⁇ ⁇ ⁇ f r ] ⁇ ( n / f s ) + ⁇ mj + ⁇ ref ) ⁇ ( 5 )
- the two components provide unique frequency contributions to the composite waveform; one is sitting at zero frequency (DC) whose magnitude is proportional to the phase difference between the optical and reference signals, and another component that is positioned at twice the modulation frequency on the Fourier spectrum.
- DC zero frequency
- the amplitude of the output signal is contingent on the amplitudes of the detected signal, the references signal, and the phase difference between them. If a fixed phase relationship is maintained and doesn't vary over time, then the final DC signal is directly proportional to the input signal. In an effort to maximize the DC value, one would normally be required to incorporate various phase shifting schemes. This problem is compounded in our application because each channel retains a distinctive phase, a consequence of the sampling process. Therefore, employing independent phase shifters would be imposed on the system for each detector channel and each wavelength. In some embodiments, to avoid the complexity of having to integrate many individual phase shifters, a solution is provided that makes use of an additional mixing stage with a quadrature reference signal.
- This method yields a value that is independent of the phase differences and is directly related to the amplitude of channel j.
- the lock in detection can be, but is not limited to, a phase-insensitive quadrature method.
- Digital detection or lock-in techniques employing phase-diversity or phase-dependent detection, with or without utilizing a phase shifter, are also provided in some embodiments of the invention.
- the only quantity we wish to preserve is the one straddling zero frequency. As has been shown before, this value is directly proportional to the amplitude of the optical input signal.
- the other frequency components should be eliminated via the low-pass filter. Elimination of the other frequency components can be achieved via a low-pass filter. Effective extraction of a single frequency band while suppressing all others, has direct implications on the bandwidth of the filter chosen and of the spacing between two modulation frequencies to ensure that NONE of the unwanted fractions falls into the passband of the filter.
- the lock-in filter can be an ordinary averaging filter.
- the following description, which references FIGS. 9-11 relates generally to a digital lock-in amplification scheme that is optimal for lock-in detection and optical imaging, including continuous wave (CW) imaging, according to some embodiments of the invention.
- the digital lock-in amplification scheme is for simultaneously measuring amplitude and phase of multiple modulated optical signals used in continuous wave optical imaging.
- a digital lock-in algorithm according to some embodiments of the invention is described, including being described under certain sampling and modulation conditions, and results are shown for sample data.
- the detection scheme, after careful choice of sampling and modulation frequencies, can be excellent or ideal for isolating signals modulated by different frequencies. Its filter also provides good or the best non-DC rejection vs. time response characteristics.
- FIG. 9 is a block diagram depicting hardware components 900 of a digital CW instrument, according to one embodiment of the invention.
- Photo detector 902 outputs a current signal, I(t), proportional to the incident light signal, which is then converted to a voltage, V(t), by a trans impedance I-V converter 904 and then amplified by a voltage amplifier 906 with gain G.
- An antialiasing filter 908 , h(t) makes sure there are no frequencies above the Nyquist frequency.
- ADC 910 samples the signal, F(t), at a frequency f s controlled by a CPLD controller 912 .
- the digital signal F[n] is then sent though the CPLD 912 on the way to DSP 914 where the lock-in algorithm, according to some embodiments of the invention, is performed and outputs the amplitude and phase for every modulation frequency.
- the lock-in algorithm every time an imaging system wants to measure the amplitude and phase of a signal, the lock-in algorithm actually takes N s samples at sampling frequency f s .
- the magnitude response of the averaging filter g[n] has a stopband attenuation of 20*log(N s ) and can be approximated as:
- (15) which clearly has zeroes at f f m .
- FIG. 10 is a set 1000 of graphs 1002 , 1004 depicting frequency response of filter, original, and modulated signals, according to one embodiment of the invention.
- FIG. 10 shows that when a set of modulated signals are modulated by a reference, the only signal that does not fall into a zero bin is when the two signals have the same frequency. When frequencies other than specified are used, inaccurate measures of amplitude and phase result.
- Red and green LED's were modulated in a very noisy environment by 4 kHz and 6 kHz wave respectively and the shown onto a single photo detector, such as the photodetector 902 depicted in FIG. 9 .
- a single photo detector such as the photodetector 902 depicted in FIG. 9 .
- For convenience analog signals where sampled with an oscilloscope after the filter stage in FIG. 9 . Data was saved to a disk and the algorithm was performed using MATLAB. The amplitudes of the signals were kept at a constant value or turned off completely. Measurements were taken with just the red LED on, just the green LED on, and both of them on. The calculated amplitude and phase were used to reconstruct the measured signals in order to verify that the results were correct.
- FIG. 11 depicts a set 1100 of graphs 1102 , 1104 , 1106 that depict plots of these calculated signals, as described above, along with the measured signals.
- the graphs 1102 , 1104 , 1106 depict measured signals, calculated amplitudes, and reconstructed signals. Calculated amplitudes are shown by each plot.
- circuit construction and signal conditioning to accommodate desired bandwidth, sensitivity, and temporal response, with regard to embedded digital signal processing techniques used in optical tomography.
- the architecture introduced herein provides a versatile, compact, high precision optical-measurement system by applying integrated digital data acquisition and digital signal processing techniques.
- digital-based systems offer distinct functional benefits over their analog counterpart including noise reduction, maintaining signal integrity, and enhancing overall performance.
- a detection device unlike existing digital systems, implements the data-acquisition and filtering via a signal processor embedded in-system conserving costs, reducing power and size, increasing speed, and generating real-time measurements.
- the circuitry was fashioned for swift data transfers and a settling time reduction of both the analog front-end electronics, and the digital filter compared to the former system. This affords faster source-switching speeds, yielding high acquisition rates which facilitate the study of functional imaging.
- FIG. 12 is a block diagram 1200 depicting system layout and signal communication, according to one embodiment of the invention.
- DSP 1202 is the component by which the signals are collected, processed and filtered, and finally routed to a host PC 1208 .
- the system uses Model ADSP-21161 with Super Harvard Architecture Computer (SHARC), available from Analog Devices, Inc., whose key features include 600 MFLOPS peak, 32-bit and 40-bit floating point arithmetic and user-configurable 1 MBits on-chip SRAM memory.
- SHARC Super Harvard Architecture Computer
- the system employs parallel detection to assist in achieving a high temporal resolution.
- the design accepts up to 32 independent detectors simultaneously whose distribution comprises 8 Printed Circuit Boards (PCB), each accommodating 4 channels.
- PCB Printed Circuit Boards
- a single 4-channel high-speed 16-bit Analog to-Digital Converter (ADC) 1204 is mounted on each board 1206 .
- a timing board 1210 consisting of the Complex Programmable Logic Device (CPLD), clock, and some support circuitry is interfaced with the DSP 1202 and detector boards 1206 .
- the USB protocol is used to transfer data between the DSP 1202 and the host PC 1208 .
- the USB chip with its support circuitry and interface glue logic is mounted on a separate PCB.
- the photodiode and trans-impedance amplifier 1210 of the system of FIG. 12 are similar to those of the system described in Schmitz 2002. Silicone-based photodiodes provide the sensitivity, frequency response, and linearity requisite for the application. These attributes combined with their low cost and low dark current make them optimal for the design.
- the trans-impedance amplifier (TIA) 1210 was carefully selected to keep the existing and inherent noise minimal.
- the feedback network of the TIA 1210 and a subsequent PGA 1212 stage offer a multitude of signal gain stages to provide a large global dynamic range for each channel and to maximize the resolution of the 16-bit ADC.
- a gain of 1, 10, or 10,000 kV/A is available from the TIA 1210 feedback while an additional gain of 1, 10, or 100 V/V is offered by the PGA yielding a dynamic range exceeding 120 dB and an absolute maximum signal gain of 10 9 V/A (180 dB).
- a 16-bit ADC offers up to 0.0015% resolution full-scale, so the system can theoretically detect 0.76% fluctuations for signals as small as 10 pA, the approximate noise floor.
- the analog signal before it must be digitized, it must be sent through a low-pass anti-aliasing filter 1214 so that the frequency spectrum conforms to the Nyquist limitation. Since the sampling rate is ⁇ 75 k samples/sec and the modulation frequencies of the lasers are from 1-9 kHz, the filter's impulse response can be enhanced at the expense of the frequency roll-off.
- the filter chosen is an 8th order Butterworth, which afford a balance between the time versus frequency response tradeoff. With a cutoff frequency of 12.5 kHz, one achieves 0.1% attenuation at the upper limit of modulation (9 kHz). By sampling at 75 ksps, noise above 67 kHz can alias into the passband.
- the filter attenuation at this frequency is ⁇ 117 dB which is adequate noise suppression.
- An impulse at the input of the LPF yields a filter settling time of 450 ⁇ sec, not a limiting factor in the general detection temporal response.
- the ADC has a 4-channel multiplexed input and a serial output which directly interfaces with the DSP 1202 . Its maximum sampling rate of 1 Msps assures that every channel is sampled at the desired 75 ksps.
- the detection scheme assumes an amplitude-modulated light source, which allows one to execute the lock-in operation.
- the offset errors, nonlinearities, and drifts are greatly reduced or virtually eliminated.
- the digitized data signals from all detectors are time-division-multiplexed onto a double serial line.
- the DSP 1202 receives the data from the ADC 1204 , it collates the signals and performs the homodyne detection algorithm. This technique mixes the optical signal with a DSP-generated reference signal at an identical frequency and sends the output through a low-pass filter.
- the result is a DC signal whose magnitude is proportional to magnitudes of the optical power and reference voltage, and their phase difference.
- the final output is independent of phase difference and is only contingent on the amplitude of the inputs. In some embodiments, this is an attractive option for some systems as described herein, for example, since all channels do not share a common phase, an effect of the sampling sequence.
- the homodyne algorithm can be easily extended to accommodate multiple wavelengths, and the low-pass-filter 1214 is specifically designed with this consideration.
- a low-pass filter 1214 used for the lock-in stage is an ordinary averaging filter; however, in other embodiments, other filters may be used.
- N is the number of samples acquired
- k is any positive integer
- f Nyquist f s /2
- f m is the modulation frequency.
- f s and N s are chosen based on based on various resulting tradeoffs.
- T acq is the data acquisition period for a single channel
- T x is the period of parameter x.
- FIG. 13 is a graph 1300 depicting frequency response of a 20 point averaging filter, with the horizontal frequency axis running from zero to Nyquist frequency.
- the DSP With respect to collecting data, once the DSP has completed its calculation, the data must be streamlined into the host computer.
- the DSP has many responsibilities and processes accessing its data busses and internal memory, immediately after performing the lock-in computations the DSP deposits the data values in an external memory bank. There, it can be accessed by a timed transfer protocol and directed to the host PC without causing a bottleneck in the digital processor.
- data latches and control circuitry are used to coordinate the gain-bits being downloaded from the host PC into the signal processor.
- bi-directional driver and transceiver devices are implemented to all data transmissions. Streamlining measurement results or gain-bit values can be accomplished via a microcontroller-based USB protocol.
- USB offer many distinct benefits, making it the interface-bus of choice. It makes the instrument portable, eliminates the need for a data-acquisition board, presents a plug-and-play design, and can support transfer speeds up to 480 Mb/s. This electrically simple user experience requires some inner complexity and adds a whole protocol layer to the system interface.
- LabVIEW is the PC control software, available from National Instruments Corp., which integrates embedded C-code to communicate with the USB's microcontroller. The timing sequences required for data transfers between the host PC and embedded DSP are synchronized by a designated controller or programmable logic device.
- lock-in detection is performed through software algorithms, enhancing signal integrity by eliminating the analog shortcomings, and making it configurable to the user's specific application.
- the resulting instrument eliminates many adverse analog affects and offer greater functionality than related systems.
- the invention provides systems, methods, and apparatuses associated with an illuminating source for a digital optical imager, such as a digital optical imager usable in optical tomography.
- these methods, systems, and apparatuses can be used, for example, with systems and methods described previously herein, including, for example, the system 200 depicted in FIG. 2 , as described in detail above.
- FIG. 14 is a flow diagram depicting a method 1400 according to one embodiment of the invention, for generating multi-wavelength illuminating sources at multiple locales, as required for attaining spectroscopic tomography.
- Step 1402 includes establishing modulation frequency for each wavelength.
- Step 1404 includes driving lasers.
- Step 1404 includes performing multiple wavelength coupling.
- Step 1406 includes optical time-division multiplexing.
- phase sensitive detection scheme since a phase sensitive detection scheme is employed to extract the amplitude of the optical signal, it is required to amplitude-modulate the light source that probes the tissue and provide a reference frequency equal to and phase-locked with the modulation frequency.
- each requires their own distinct modulation frequency to encode the optical intensity and reference signals to decode it.
- Direct digital synthesis is a method of producing an analog sine waveform by generating a time-varying signal in digital form and then performing a digital-to-analog conversion. Because operations within a DDS device are primarily digital, it can offer fast switching between output frequencies, fine frequency resolution, and operation over a broad spectrum of frequencies. With advances in design and process technology, modern day DDS devices are very compact and consume little power. It therefore serves as a stable and accurate waveform generator to produce the frequency stimulus required for the application. Each DDS device outputs one frequency and therefore supplies only a single illuminating wavelength.
- the design uses a synthesizer that provides two sine waves, in-phase and quadrature, which are 90 degrees phase shifted from each other. This provides the option of sampling both components to use as the reference signals in quadrature lock-in detection, as described in detail above.
- the first choice is to digitally generate these reference signals internally in the DSP thereby utilizing the DDS synthesizers only to modulate the lasers. In some embodiments, however, other options may be utilized, such as digitizing the analog reference waveforms instead.
- the DDS contains dual integrated 12-bit digital-to-analog (DAC) converters for high resolution waveforms and demonstrates excellent dynamic performance.
- a low pass filter is applied to the produced waveform in an effort to lower the harmonic content even further and make efficient use of the higher precision 16-bit analog-to-digital converters (ADC) used for sampling the detector channels.
- ADC analog-to-digital converters
- the function of the digital synthesizer is determined and controlled through a set of register values together with specified commands.
- a designated microcontroller is assigned to interfacing directly with the DDS chips. Programming the relevant instructions required for generating specific waveforms as well as timing organization will be executed via this controller and/or the digital signal processor. Selecting the modulation frequencies can be preprogrammed into the controller or modified spontaneously through computer software or autonomous devices, for example.
- Some embodiments provide an instrument that uses laser diodes in the near-infrared region as the radiative source for optical imaging.
- the lasers are regulated by select OEM laser diode controllers that combine a low noise, low drift current source with a precise thermoelectric cooler (TEC).
- TEC thermoelectric cooler Operating in constant current mode creates a laser output whose power is proportional to the driving current. Therefore, once the respective frequencies are generated, one can use them to modulate the driving current of the light source, which in effect, modulates the laser output.
- light transmission data must be collected on multiple wavelengths.
- the minimum number of wavelengths for performing spectral analysis is two. By combining these wavelengths together and making simultaneous measurements, one can shorten the data acquisition period and increase the temporal response of the instrument.
- the all that's left is to spatially encode multiple illumination locales around the target being probed.
- a method frequently employed to separate source locations is through time-division multiplexing of the combined light source. This can be realized by a 1 ⁇ N optical switch where N is the number of different source positions.
- a single input is routed to a distinct output which is controlled by a specific address. When the switch receives an appropriate address, the input is aligned with the respective output.
- the settling time of the switch is on the order of 2-5 msec (a significant feature for acquiring dynamic images).
- a problem with commercially available switches is that the inter-channel crosstalk is only specified to about ⁇ 50 to ⁇ 60 dB while the dynamic range of the instrument is effectively around 90 dB of optical power (when comparably adjusted). This means that light can leak from one channel to another contaminating the measured quantity.
- one uses a configuration of a 1 ⁇ N switch in series with a 2 ⁇ 2 optical switch that suppresses crosstalk by over ⁇ 100 dB.
- FIG. 15 depicts one embodiment of such a configuration.
- FIG. 15 is a block diagram 1500 depicting a 1 ⁇ 16 optical switch 1502 in series with 8-2 ⁇ 2 switches 1504 to provide a composite 1 ⁇ 16 multiplexed source architecture with an effective crosstalk elimination of over 100 dB.
- the output of the optical switch then interfaces directly with the tissue or measurement probe.
- This entire apparatus can function independently or be controlled by a personal computer, for example. Offering variable parameter control will allow for spontaneous user command.
- an LCD display may be integrated into this setup to provide the user with some system monitoring and feedback.
Abstract
Description
TABLE 1 | |||||
G2 | G1 | G0 | TIA gain (Ω) | PGA gain (V/V) | Overall Gain (Ω) |
0 | 0 | 0 | |
1 | |
0 | 0 | 1 | |
1 | |
0 | 1 | 0 | |
10 | |
0 | 1 | 1 | |
100 | 1 | |
1 | 0 | 0 | 10 |
1 | 10 | |
1 | 0 | 1 | 10 |
1 | 10 | |
1 | 1 | 0 | 10 |
10 | 100 | |
1 | 1 | 1 | 10 M | 100 | 1 | G |
f λi = m i ( f s /K), where: (2)
Filter Settling time T F =K/f s, where: (3)
A λi=2√{square root over (I λi 2 +Q λi 2)}. (4b)
I mj [n]=V mj ·S ref =A mj A r(½){cos(φmj−φref)−cos([4πf m(n/f s)]+φmj+φref)} (6)
X mj [n]=A mj A r(½){cos(φmj−φref)}=A mj A r(½){cos(θj)} (7)
Q mj [n]=A mj sin(2πf m n/f s+Φmj)·A r cos(2πf r n/f s+Φref)
Y mj [n]=A mj A r(½){sin(φmj−φref)}=A mj A r(½){sin(θj)} (8)
V f1 +V f2 =V 1 sin(f 1 t)+V 2 sin(f 2 t)
f1−f1=0 Hz=DC component (I)
f1+f1=3 kHz+3 kHz=6 kHz (II)
f2−f1=5 kHz−3 kHz=2 kHz (III)
f2+f1=5 kHz+3 kHz=8 kHz (IV)
T acq =N s /f s (11)
Since the sampling is controlled by some digital controller with a clock frequency fclk. Only certain values of fs are allowed:
f s =f clk /N c (12)
where Nc is the integer number of clock cycles it takes the digital controller to run through a sampling period.
T acq −1 s f m =k·f s /N s (13)
g[n]=1/N s (14)
to both remove noise and isolate signals. The magnitude response of the averaging filter g[n] has a stopband attenuation of 20*log(Ns) and can be approximated as:
|G(f)|=|sin c(N s πf/2f s)| (15)
which clearly has zeroes at f=fm.
F m(t)=A m cos(2πf m t+φ m) (16)
then the digital signal is:
Fm[n]=A m cos [2ππm n/fs+φ m]. (17)
The first step of the algorithm is to modulate the signal by the stored sequences for every modulation frequency to produce 2 new sequences:
I m [n]=C m [n]·F[n]=A m/2·cos(φm)+A m/2·cos [4πf m n/f s φm] (17)
and:
Q m [n]=S m [n]·F[n]=A m/2·cos(φm)+A m/2·sin [4πf m n/f s +φm]. (19)
These sequences are then filtered by averaging filter. This is the same thing as taking the mean of the sequence to result in the following values:
A m=2√{square root over ((X m 2 +Y m 2))}. (22)
The phase φm is found by taking the quadrant dependent inverse tan function
G(f i +f j)=G(f i −f j)=0 (i≠j). (25)
f zero =k[f Nyquist/(N s/2)]=k[f s/(N s)],f m =f zero (26)
T m =T zero=[(N s /k)T s ], T acq =N s T ss so T m =T acq /k or T acq /T m =k (27)
Claims (27)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/257,825 US7728986B2 (en) | 2004-04-13 | 2008-10-24 | Systems and methods for digital detection of a tomographic signal |
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US56198904P | 2004-04-13 | 2004-04-13 | |
US11/106,348 US7463362B2 (en) | 2004-04-13 | 2005-04-13 | Digital signal processor-based detection system, method, and apparatus for optical tomography |
US12/257,825 US7728986B2 (en) | 2004-04-13 | 2008-10-24 | Systems and methods for digital detection of a tomographic signal |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/106,348 Continuation US7463362B2 (en) | 2004-04-13 | 2005-04-13 | Digital signal processor-based detection system, method, and apparatus for optical tomography |
Publications (2)
Publication Number | Publication Date |
---|---|
US20090103102A1 US20090103102A1 (en) | 2009-04-23 |
US7728986B2 true US7728986B2 (en) | 2010-06-01 |
Family
ID=36060458
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/106,348 Active 2026-05-24 US7463362B2 (en) | 2004-04-13 | 2005-04-13 | Digital signal processor-based detection system, method, and apparatus for optical tomography |
US12/257,825 Active 2025-05-19 US7728986B2 (en) | 2004-04-13 | 2008-10-24 | Systems and methods for digital detection of a tomographic signal |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US11/106,348 Active 2026-05-24 US7463362B2 (en) | 2004-04-13 | 2005-04-13 | Digital signal processor-based detection system, method, and apparatus for optical tomography |
Country Status (2)
Country | Link |
---|---|
US (2) | US7463362B2 (en) |
WO (1) | WO2006031258A2 (en) |
Cited By (69)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090046295A1 (en) * | 2007-07-12 | 2009-02-19 | Volcano Corporation | Apparatus and methods for uniform sample clocking |
US20100292569A1 (en) * | 2008-03-14 | 2010-11-18 | The Trustees Of Columbia University In The City Of New York | Systems and methods for dynamic imaging of tissue using digital optical tomography |
US20110255863A1 (en) * | 2010-04-14 | 2011-10-20 | Nanguang Chen | Integral transformed optical measurement method and apparatus |
US8395781B2 (en) | 2007-07-12 | 2013-03-12 | Volcano Corporation | Automatic calibration systems and methods of use |
US20130080097A1 (en) * | 2011-09-26 | 2013-03-28 | Belkin International Inc. | Systems and methods to isolate lower amplitude signals for analysis in the presence of large amplitude transients |
US20130338496A1 (en) * | 2010-12-13 | 2013-12-19 | The Trustees Of Columbia University In The City New York | Medical imaging devices, methods, and systems |
US9286673B2 (en) | 2012-10-05 | 2016-03-15 | Volcano Corporation | Systems for correcting distortions in a medical image and methods of use thereof |
US9292918B2 (en) | 2012-10-05 | 2016-03-22 | Volcano Corporation | Methods and systems for transforming luminal images |
US9301687B2 (en) | 2013-03-13 | 2016-04-05 | Volcano Corporation | System and method for OCT depth calibration |
US9307926B2 (en) | 2012-10-05 | 2016-04-12 | Volcano Corporation | Automatic stent detection |
US9324141B2 (en) | 2012-10-05 | 2016-04-26 | Volcano Corporation | Removal of A-scan streaking artifact |
US9360630B2 (en) | 2011-08-31 | 2016-06-07 | Volcano Corporation | Optical-electrical rotary joint and methods of use |
US9367965B2 (en) | 2012-10-05 | 2016-06-14 | Volcano Corporation | Systems and methods for generating images of tissue |
US9383263B2 (en) | 2012-12-21 | 2016-07-05 | Volcano Corporation | Systems and methods for narrowing a wavelength emission of light |
US9478940B2 (en) | 2012-10-05 | 2016-10-25 | Volcano Corporation | Systems and methods for amplifying light |
US9486143B2 (en) | 2012-12-21 | 2016-11-08 | Volcano Corporation | Intravascular forward imaging device |
US9596993B2 (en) | 2007-07-12 | 2017-03-21 | Volcano Corporation | Automatic calibration systems and methods of use |
US9612105B2 (en) | 2012-12-21 | 2017-04-04 | Volcano Corporation | Polarization sensitive optical coherence tomography system |
US9622706B2 (en) | 2007-07-12 | 2017-04-18 | Volcano Corporation | Catheter for in vivo imaging |
US9709379B2 (en) | 2012-12-20 | 2017-07-18 | Volcano Corporation | Optical coherence tomography system that is reconfigurable between different imaging modes |
US9730613B2 (en) | 2012-12-20 | 2017-08-15 | Volcano Corporation | Locating intravascular images |
US9766277B2 (en) | 2009-09-25 | 2017-09-19 | Belkin International, Inc. | Self-calibrating contactless power consumption sensing |
US9770172B2 (en) | 2013-03-07 | 2017-09-26 | Volcano Corporation | Multimodal segmentation in intravascular images |
US9858668B2 (en) | 2012-10-05 | 2018-01-02 | Volcano Corporation | Guidewire artifact removal in images |
US9857449B2 (en) | 2010-07-02 | 2018-01-02 | Belkin International, Inc. | System and method for monitoring electrical power usage in an electrical power infrastructure of a building |
US9867530B2 (en) | 2006-08-14 | 2018-01-16 | Volcano Corporation | Telescopic side port catheter device with imaging system and method for accessing side branch occlusions |
US10058284B2 (en) | 2012-12-21 | 2018-08-28 | Volcano Corporation | Simultaneous imaging, monitoring, and therapy |
US10070827B2 (en) | 2012-10-05 | 2018-09-11 | Volcano Corporation | Automatic image playback |
US10166003B2 (en) | 2012-12-21 | 2019-01-01 | Volcano Corporation | Ultrasound imaging with variable line density |
US10191220B2 (en) | 2012-12-21 | 2019-01-29 | Volcano Corporation | Power-efficient optical circuit |
US10219780B2 (en) | 2007-07-12 | 2019-03-05 | Volcano Corporation | OCT-IVUS catheter for concurrent luminal imaging |
US10219887B2 (en) | 2013-03-14 | 2019-03-05 | Volcano Corporation | Filters with echogenic characteristics |
US10226597B2 (en) | 2013-03-07 | 2019-03-12 | Volcano Corporation | Guidewire with centering mechanism |
US10238367B2 (en) | 2012-12-13 | 2019-03-26 | Volcano Corporation | Devices, systems, and methods for targeted cannulation |
US10247765B2 (en) | 2007-09-18 | 2019-04-02 | Georgia Tech Research Corporation | Detecting actuation of electrical devices using electrical noise over a power line |
US10292677B2 (en) | 2013-03-14 | 2019-05-21 | Volcano Corporation | Endoluminal filter having enhanced echogenic properties |
US10332228B2 (en) | 2012-12-21 | 2019-06-25 | Volcano Corporation | System and method for graphical processing of medical data |
US10413317B2 (en) | 2012-12-21 | 2019-09-17 | Volcano Corporation | System and method for catheter steering and operation |
US10420530B2 (en) | 2012-12-21 | 2019-09-24 | Volcano Corporation | System and method for multipath processing of image signals |
US10426590B2 (en) | 2013-03-14 | 2019-10-01 | Volcano Corporation | Filters with echogenic characteristics |
US10459012B2 (en) | 2010-07-02 | 2019-10-29 | Belkin International, Inc. | System for monitoring electrical power usage of a structure and method of same |
US10506181B2 (en) | 2018-03-31 | 2019-12-10 | Open Water Internet Inc. | Device for optical imaging |
US10568586B2 (en) | 2012-10-05 | 2020-02-25 | Volcano Corporation | Systems for indicating parameters in an imaging data set and methods of use |
US10595820B2 (en) | 2012-12-20 | 2020-03-24 | Philips Image Guided Therapy Corporation | Smooth transition catheters |
US10638939B2 (en) | 2013-03-12 | 2020-05-05 | Philips Image Guided Therapy Corporation | Systems and methods for diagnosing coronary microvascular disease |
US10724082B2 (en) | 2012-10-22 | 2020-07-28 | Bio-Rad Laboratories, Inc. | Methods for analyzing DNA |
US10758207B2 (en) | 2013-03-13 | 2020-09-01 | Philips Image Guided Therapy Corporation | Systems and methods for producing an image from a rotational intravascular ultrasound device |
US10778911B2 (en) | 2018-03-31 | 2020-09-15 | Open Water Internet Inc. | Optical transformation device for imaging |
US10778912B2 (en) | 2018-03-31 | 2020-09-15 | Open Water Internet Inc. | System and device for optical transformation |
US10772574B2 (en) | 2016-09-13 | 2020-09-15 | Open Water Internet Inc. | Imaging with infrared imaging signals |
US10874370B2 (en) | 2019-01-28 | 2020-12-29 | Open Water Internet Inc. | Pulse measurement in optical imaging |
US10942022B2 (en) | 2012-12-20 | 2021-03-09 | Philips Image Guided Therapy Corporation | Manual calibration of imaging system |
US10939826B2 (en) | 2012-12-20 | 2021-03-09 | Philips Image Guided Therapy Corporation | Aspirating and removing biological material |
US10955406B2 (en) | 2019-02-05 | 2021-03-23 | Open Water Internet Inc. | Diffuse optical imaging with multiple beams |
US10962929B2 (en) | 2018-09-14 | 2021-03-30 | Open Water Internet Inc. | Interference optics for optical imaging device |
US10966612B2 (en) | 2018-06-14 | 2021-04-06 | Open Water Internet Inc. | Expanding beam optical element |
US10993694B2 (en) | 2012-12-21 | 2021-05-04 | Philips Image Guided Therapy Corporation | Rotational ultrasound imaging catheter with extended catheter body telescope |
US11026591B2 (en) | 2013-03-13 | 2021-06-08 | Philips Image Guided Therapy Corporation | Intravascular pressure sensor calibration |
US11040140B2 (en) | 2010-12-31 | 2021-06-22 | Philips Image Guided Therapy Corporation | Deep vein thrombosis therapeutic methods |
US11141063B2 (en) | 2010-12-23 | 2021-10-12 | Philips Image Guided Therapy Corporation | Integrated system architectures and methods of use |
US11154313B2 (en) | 2013-03-12 | 2021-10-26 | The Volcano Corporation | Vibrating guidewire torquer and methods of use |
US11259706B2 (en) | 2020-05-19 | 2022-03-01 | Open Water Internet Inc. | Dual wavelength imaging and out of sample optical imaging |
US11272845B2 (en) | 2012-10-05 | 2022-03-15 | Philips Image Guided Therapy Corporation | System and method for instant and automatic border detection |
US11320370B2 (en) | 2019-06-26 | 2022-05-03 | Open Water Internet Inc. | Apparatus for directing optical and acoustic signals |
US11406498B2 (en) | 2012-12-20 | 2022-08-09 | Philips Image Guided Therapy Corporation | Implant delivery system and implants |
US11559208B2 (en) | 2020-05-19 | 2023-01-24 | Open Water Internet Inc. | Imaging with scattering layer |
US11581696B2 (en) | 2019-08-14 | 2023-02-14 | Open Water Internet Inc. | Multi-channel laser |
US11622686B2 (en) | 2019-11-22 | 2023-04-11 | Open Water Internet, Inc. | Optical imaging with unshifted reference beam |
US11819318B2 (en) | 2020-04-27 | 2023-11-21 | Open Water Internet Inc. | Optical imaging from light coherence |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006031258A2 (en) * | 2004-04-13 | 2006-03-23 | The Trustees Of Columbia University In The City Of New York | Digital signal processor-based detection system, method, and apparatus for optical tomography |
US7616311B2 (en) * | 2005-05-02 | 2009-11-10 | Jmar Llc | Systems and methods for a multiple angle light scattering (MALS) instrument having two-dimensional detector array |
US10884712B1 (en) | 2005-12-30 | 2021-01-05 | The Mathworks, Inc. | Component-based framework for generating device driver model elements |
JP5138701B2 (en) * | 2006-12-22 | 2013-02-06 | ワシントン・ユニバーシティ | High performance imaging system for diffuse optical tomography and associated methods of use |
US9784561B2 (en) * | 2007-01-19 | 2017-10-10 | Thorlabs, Inc. | Optical coherence tomography imaging system and method |
US8705047B2 (en) | 2007-01-19 | 2014-04-22 | Thorlabs, Inc. | Optical coherence tomography imaging system and method |
US7936462B2 (en) * | 2007-01-19 | 2011-05-03 | Thorlabs, Inc. | Optical coherence tomography imaging system and method |
US9347765B2 (en) * | 2007-10-05 | 2016-05-24 | Volcano Corporation | Real time SD-OCT with distributed acquisition and processing |
US7709794B2 (en) * | 2008-02-05 | 2010-05-04 | Kla-Tencor Corporation | Defect detection using time delay lock-in thermography (LIT) and dark field LIT |
US7973688B2 (en) * | 2008-04-07 | 2011-07-05 | Nec Laboratories America, Inc. | Photonic filter bank for high-speed analog-to-digital conversion |
TWI412940B (en) * | 2009-10-06 | 2013-10-21 | Univ Nat Chiao Tung | Image reconstruction method, device and computer program for diffuse optical tomography |
JP5538822B2 (en) * | 2009-11-06 | 2014-07-02 | キヤノン株式会社 | Ultrasonic detection device and ultrasonic diagnostic device |
JP6373005B2 (en) | 2010-11-11 | 2018-08-15 | ザ トラスティーズ オブ コロンビア ユニバーシティ イン ザ シティ オブ ニューヨーク | Dynamic optical tomography imaging apparatus, method and system |
US8988684B1 (en) * | 2011-09-08 | 2015-03-24 | Lawrence Livermore National Security, Llc | System and method for measuring fluorescence of a sample |
USD763938S1 (en) | 2014-04-02 | 2016-08-16 | Cephalogics, LLC | Optical sensor array |
USD763939S1 (en) | 2014-04-02 | 2016-08-16 | Cephalogics, LLC | Optical sensor array liner with optical sensor array pad |
US9869542B2 (en) * | 2014-04-21 | 2018-01-16 | Axsun Technologies, Inc. | System and method for resampling optical coherence tomography signals in segments |
ES2701334T3 (en) * | 2016-08-01 | 2019-02-21 | Knowledge Dev For Pof Sl | LED controller for high speed optical communications based on linear modulations |
US10627349B2 (en) * | 2018-02-21 | 2020-04-21 | Olive Healthcare Inc. | Signal processing device of analyzing bio-signal and bio-signal analyzing apparatus using the same |
DE102018206434A1 (en) * | 2018-04-25 | 2019-10-31 | Aktiebolaget Skf | Signal processing method and apparatus |
CN111505701B (en) * | 2020-05-27 | 2023-03-03 | 中国电子科技集团公司第三十八研究所 | Compact accelerator beam phase and position measuring system and method |
US20230401765A1 (en) * | 2020-11-03 | 2023-12-14 | Institut National De La Recherche Scientifique | A method and a system for compressed ultrafast tomographic imaging |
CN114631832B (en) * | 2022-01-25 | 2023-11-17 | 中国科学院自动化研究所 | Muscle signal acquisition device |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4807146A (en) * | 1986-02-19 | 1989-02-21 | Louisiana State University | Digital lock-in amplifier |
US5746211A (en) * | 1993-11-30 | 1998-05-05 | Bell Communications Research, Inc. | Absorption imaging system and method using direct reconstruction of scattered radiation |
US6011401A (en) * | 1997-05-20 | 2000-01-04 | Raytheon Company | Reference-free direct digital lock-in method and apparatus |
US6549284B1 (en) * | 1999-09-17 | 2003-04-15 | The General Hospital Corporation | Calibration methods and systems for diffuse optical tomography and spectroscopy |
US6795195B1 (en) * | 1999-09-14 | 2004-09-21 | Research Foundation Of State University Of New York | System and method for tomographic imaging of dynamic properties of a scattering medium |
US6816102B2 (en) * | 2002-06-20 | 2004-11-09 | Dakota Technologies, Inc. | System for digitizing transient signals |
US20050243322A1 (en) * | 2004-04-13 | 2005-11-03 | Lasker Joseph M | Digital signal processor-based detection system, method, and apparatus for optical tomography |
US20080154126A1 (en) * | 2006-12-22 | 2008-06-26 | Washington University | High Performance Imaging System for Diffuse Optical Tomography and Associated Method of Use |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63318811A (en) * | 1987-06-22 | 1988-12-27 | Toshiba Corp | Digital filter device |
US5137355A (en) * | 1988-06-08 | 1992-08-11 | The Research Foundation Of State University Of New York | Method of imaging a random medium |
US5010455A (en) * | 1990-05-07 | 1991-04-23 | General Motors Corporation | Headlamp assembly |
US5220589A (en) * | 1991-07-18 | 1993-06-15 | General Electric Company | Correction circuit for a floating-point amplifier |
US5528365A (en) * | 1994-03-01 | 1996-06-18 | The Trustees Of The University Of Pennsylvania | Methods and apparatus for imaging with diffuse light |
EP1611411A2 (en) * | 2003-03-26 | 2006-01-04 | Southwest Sciences Incorporated | Method and apparatus for imaging internal structures of transparent and translucent materials |
-
2005
- 2005-04-13 WO PCT/US2005/012692 patent/WO2006031258A2/en active Application Filing
- 2005-04-13 US US11/106,348 patent/US7463362B2/en active Active
-
2008
- 2008-10-24 US US12/257,825 patent/US7728986B2/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4807146A (en) * | 1986-02-19 | 1989-02-21 | Louisiana State University | Digital lock-in amplifier |
US5746211A (en) * | 1993-11-30 | 1998-05-05 | Bell Communications Research, Inc. | Absorption imaging system and method using direct reconstruction of scattered radiation |
US6011401A (en) * | 1997-05-20 | 2000-01-04 | Raytheon Company | Reference-free direct digital lock-in method and apparatus |
US6795195B1 (en) * | 1999-09-14 | 2004-09-21 | Research Foundation Of State University Of New York | System and method for tomographic imaging of dynamic properties of a scattering medium |
US6549284B1 (en) * | 1999-09-17 | 2003-04-15 | The General Hospital Corporation | Calibration methods and systems for diffuse optical tomography and spectroscopy |
US6816102B2 (en) * | 2002-06-20 | 2004-11-09 | Dakota Technologies, Inc. | System for digitizing transient signals |
US20050243322A1 (en) * | 2004-04-13 | 2005-11-03 | Lasker Joseph M | Digital signal processor-based detection system, method, and apparatus for optical tomography |
US7463362B2 (en) * | 2004-04-13 | 2008-12-09 | The Trustees Of Columbia University In The City Of New York | Digital signal processor-based detection system, method, and apparatus for optical tomography |
US20090103102A1 (en) * | 2004-04-13 | 2009-04-23 | The Trustees Of Columbia University In The City Of New York | Systems and methods for digital detection of a tomographic signal |
US20080154126A1 (en) * | 2006-12-22 | 2008-06-26 | Washington University | High Performance Imaging System for Diffuse Optical Tomography and Associated Method of Use |
Cited By (92)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9867530B2 (en) | 2006-08-14 | 2018-01-16 | Volcano Corporation | Telescopic side port catheter device with imaging system and method for accessing side branch occlusions |
US20090046295A1 (en) * | 2007-07-12 | 2009-02-19 | Volcano Corporation | Apparatus and methods for uniform sample clocking |
US9622706B2 (en) | 2007-07-12 | 2017-04-18 | Volcano Corporation | Catheter for in vivo imaging |
US8049900B2 (en) * | 2007-07-12 | 2011-11-01 | Volcano Corporation | Apparatus and methods for uniform frequency sample clocking |
US8395781B2 (en) | 2007-07-12 | 2013-03-12 | Volcano Corporation | Automatic calibration systems and methods of use |
US11350906B2 (en) | 2007-07-12 | 2022-06-07 | Philips Image Guided Therapy Corporation | OCT-IVUS catheter for concurrent luminal imaging |
US8593641B2 (en) | 2007-07-12 | 2013-11-26 | Volcano Corporation | Apparatus and methods for uniform frequency sample clocking |
US10219780B2 (en) | 2007-07-12 | 2019-03-05 | Volcano Corporation | OCT-IVUS catheter for concurrent luminal imaging |
US9596993B2 (en) | 2007-07-12 | 2017-03-21 | Volcano Corporation | Automatic calibration systems and methods of use |
US10247765B2 (en) | 2007-09-18 | 2019-04-02 | Georgia Tech Research Corporation | Detecting actuation of electrical devices using electrical noise over a power line |
US11119141B2 (en) | 2007-09-18 | 2021-09-14 | Georgia Tech Research Corporation | Detecting actuation of electrical devices using electrical noise over a power line |
US20150223697A1 (en) * | 2008-03-14 | 2015-08-13 | The Trustees Of Columbia University In The City Of New York | Systems and methods for dynamic imaging of tissue using digital optical tomography |
US9037216B2 (en) * | 2008-03-14 | 2015-05-19 | The Trustees Of Columbia University In The City Of New York | Systems and methods for dynamic imaging of tissue using digital optical tomography |
US20100292569A1 (en) * | 2008-03-14 | 2010-11-18 | The Trustees Of Columbia University In The City Of New York | Systems and methods for dynamic imaging of tissue using digital optical tomography |
US10371728B2 (en) | 2009-09-25 | 2019-08-06 | Belkin International, Inc. | Self-calibrating contactless power consumption sensing |
US9766277B2 (en) | 2009-09-25 | 2017-09-19 | Belkin International, Inc. | Self-calibrating contactless power consumption sensing |
US8649010B2 (en) * | 2010-04-14 | 2014-02-11 | Nanguang Chen | Integral transformed optical measurement method and apparatus |
US20110255863A1 (en) * | 2010-04-14 | 2011-10-20 | Nanguang Chen | Integral transformed optical measurement method and apparatus |
US10459012B2 (en) | 2010-07-02 | 2019-10-29 | Belkin International, Inc. | System for monitoring electrical power usage of a structure and method of same |
US10345423B2 (en) | 2010-07-02 | 2019-07-09 | Belkin International Inc. | System and method for monitoring electrical power usage in an electrical power infrastructure of a building |
US9857449B2 (en) | 2010-07-02 | 2018-01-02 | Belkin International, Inc. | System and method for monitoring electrical power usage in an electrical power infrastructure of a building |
US20130338496A1 (en) * | 2010-12-13 | 2013-12-19 | The Trustees Of Columbia University In The City New York | Medical imaging devices, methods, and systems |
US9486142B2 (en) * | 2010-12-13 | 2016-11-08 | The Trustees Of Columbia University In The City Of New York | Medical imaging devices, methods, and systems |
US11141063B2 (en) | 2010-12-23 | 2021-10-12 | Philips Image Guided Therapy Corporation | Integrated system architectures and methods of use |
US11040140B2 (en) | 2010-12-31 | 2021-06-22 | Philips Image Guided Therapy Corporation | Deep vein thrombosis therapeutic methods |
US9360630B2 (en) | 2011-08-31 | 2016-06-07 | Volcano Corporation | Optical-electrical rotary joint and methods of use |
US9130393B2 (en) * | 2011-09-26 | 2015-09-08 | Belkin International, Inc. | Systems and methods to isolate lower amplitude signals for analysis in the presence of large amplitude transients |
US20130080097A1 (en) * | 2011-09-26 | 2013-03-28 | Belkin International Inc. | Systems and methods to isolate lower amplitude signals for analysis in the presence of large amplitude transients |
US11890117B2 (en) | 2012-10-05 | 2024-02-06 | Philips Image Guided Therapy Corporation | Systems for indicating parameters in an imaging data set and methods of use |
US9307926B2 (en) | 2012-10-05 | 2016-04-12 | Volcano Corporation | Automatic stent detection |
US9478940B2 (en) | 2012-10-05 | 2016-10-25 | Volcano Corporation | Systems and methods for amplifying light |
US11864870B2 (en) | 2012-10-05 | 2024-01-09 | Philips Image Guided Therapy Corporation | System and method for instant and automatic border detection |
US9367965B2 (en) | 2012-10-05 | 2016-06-14 | Volcano Corporation | Systems and methods for generating images of tissue |
US9858668B2 (en) | 2012-10-05 | 2018-01-02 | Volcano Corporation | Guidewire artifact removal in images |
US10070827B2 (en) | 2012-10-05 | 2018-09-11 | Volcano Corporation | Automatic image playback |
US9292918B2 (en) | 2012-10-05 | 2016-03-22 | Volcano Corporation | Methods and systems for transforming luminal images |
US11272845B2 (en) | 2012-10-05 | 2022-03-15 | Philips Image Guided Therapy Corporation | System and method for instant and automatic border detection |
US11510632B2 (en) | 2012-10-05 | 2022-11-29 | Philips Image Guided Therapy Corporation | Systems for indicating parameters in an imaging data set and methods of use |
US9286673B2 (en) | 2012-10-05 | 2016-03-15 | Volcano Corporation | Systems for correcting distortions in a medical image and methods of use thereof |
US10568586B2 (en) | 2012-10-05 | 2020-02-25 | Volcano Corporation | Systems for indicating parameters in an imaging data set and methods of use |
US9324141B2 (en) | 2012-10-05 | 2016-04-26 | Volcano Corporation | Removal of A-scan streaking artifact |
US10724082B2 (en) | 2012-10-22 | 2020-07-28 | Bio-Rad Laboratories, Inc. | Methods for analyzing DNA |
US10238367B2 (en) | 2012-12-13 | 2019-03-26 | Volcano Corporation | Devices, systems, and methods for targeted cannulation |
US9709379B2 (en) | 2012-12-20 | 2017-07-18 | Volcano Corporation | Optical coherence tomography system that is reconfigurable between different imaging modes |
US10595820B2 (en) | 2012-12-20 | 2020-03-24 | Philips Image Guided Therapy Corporation | Smooth transition catheters |
US11141131B2 (en) | 2012-12-20 | 2021-10-12 | Philips Image Guided Therapy Corporation | Smooth transition catheters |
US11892289B2 (en) | 2012-12-20 | 2024-02-06 | Philips Image Guided Therapy Corporation | Manual calibration of imaging system |
US9730613B2 (en) | 2012-12-20 | 2017-08-15 | Volcano Corporation | Locating intravascular images |
US10939826B2 (en) | 2012-12-20 | 2021-03-09 | Philips Image Guided Therapy Corporation | Aspirating and removing biological material |
US10942022B2 (en) | 2012-12-20 | 2021-03-09 | Philips Image Guided Therapy Corporation | Manual calibration of imaging system |
US11406498B2 (en) | 2012-12-20 | 2022-08-09 | Philips Image Guided Therapy Corporation | Implant delivery system and implants |
US11253225B2 (en) | 2012-12-21 | 2022-02-22 | Philips Image Guided Therapy Corporation | System and method for multipath processing of image signals |
US10166003B2 (en) | 2012-12-21 | 2019-01-01 | Volcano Corporation | Ultrasound imaging with variable line density |
US9486143B2 (en) | 2012-12-21 | 2016-11-08 | Volcano Corporation | Intravascular forward imaging device |
US10332228B2 (en) | 2012-12-21 | 2019-06-25 | Volcano Corporation | System and method for graphical processing of medical data |
US9383263B2 (en) | 2012-12-21 | 2016-07-05 | Volcano Corporation | Systems and methods for narrowing a wavelength emission of light |
US11786213B2 (en) | 2012-12-21 | 2023-10-17 | Philips Image Guided Therapy Corporation | System and method for multipath processing of image signals |
US10058284B2 (en) | 2012-12-21 | 2018-08-28 | Volcano Corporation | Simultaneous imaging, monitoring, and therapy |
US10993694B2 (en) | 2012-12-21 | 2021-05-04 | Philips Image Guided Therapy Corporation | Rotational ultrasound imaging catheter with extended catheter body telescope |
US10191220B2 (en) | 2012-12-21 | 2019-01-29 | Volcano Corporation | Power-efficient optical circuit |
US9612105B2 (en) | 2012-12-21 | 2017-04-04 | Volcano Corporation | Polarization sensitive optical coherence tomography system |
US10420530B2 (en) | 2012-12-21 | 2019-09-24 | Volcano Corporation | System and method for multipath processing of image signals |
US10413317B2 (en) | 2012-12-21 | 2019-09-17 | Volcano Corporation | System and method for catheter steering and operation |
US10226597B2 (en) | 2013-03-07 | 2019-03-12 | Volcano Corporation | Guidewire with centering mechanism |
US9770172B2 (en) | 2013-03-07 | 2017-09-26 | Volcano Corporation | Multimodal segmentation in intravascular images |
US11154313B2 (en) | 2013-03-12 | 2021-10-26 | The Volcano Corporation | Vibrating guidewire torquer and methods of use |
US10638939B2 (en) | 2013-03-12 | 2020-05-05 | Philips Image Guided Therapy Corporation | Systems and methods for diagnosing coronary microvascular disease |
US10758207B2 (en) | 2013-03-13 | 2020-09-01 | Philips Image Guided Therapy Corporation | Systems and methods for producing an image from a rotational intravascular ultrasound device |
US9301687B2 (en) | 2013-03-13 | 2016-04-05 | Volcano Corporation | System and method for OCT depth calibration |
US11026591B2 (en) | 2013-03-13 | 2021-06-08 | Philips Image Guided Therapy Corporation | Intravascular pressure sensor calibration |
US10292677B2 (en) | 2013-03-14 | 2019-05-21 | Volcano Corporation | Endoluminal filter having enhanced echogenic properties |
US10426590B2 (en) | 2013-03-14 | 2019-10-01 | Volcano Corporation | Filters with echogenic characteristics |
US10219887B2 (en) | 2013-03-14 | 2019-03-05 | Volcano Corporation | Filters with echogenic characteristics |
US10772574B2 (en) | 2016-09-13 | 2020-09-15 | Open Water Internet Inc. | Imaging with infrared imaging signals |
US11547370B2 (en) | 2016-09-13 | 2023-01-10 | Open Water Internet Inc. | Method of infrared imaging |
US10506181B2 (en) | 2018-03-31 | 2019-12-10 | Open Water Internet Inc. | Device for optical imaging |
US11252343B2 (en) | 2018-03-31 | 2022-02-15 | Open Water Internet Inc. | Optical imaging through display |
US10778912B2 (en) | 2018-03-31 | 2020-09-15 | Open Water Internet Inc. | System and device for optical transformation |
US10778911B2 (en) | 2018-03-31 | 2020-09-15 | Open Water Internet Inc. | Optical transformation device for imaging |
US10880497B2 (en) | 2018-03-31 | 2020-12-29 | Open Water Internet Inc. | Optical imaging with ultrasonic signal |
US10966612B2 (en) | 2018-06-14 | 2021-04-06 | Open Water Internet Inc. | Expanding beam optical element |
US11320783B2 (en) | 2018-09-14 | 2022-05-03 | Open Water Internet Inc. | Flexible tip optical imaging |
US10962929B2 (en) | 2018-09-14 | 2021-03-30 | Open Water Internet Inc. | Interference optics for optical imaging device |
US10874370B2 (en) | 2019-01-28 | 2020-12-29 | Open Water Internet Inc. | Pulse measurement in optical imaging |
US10955406B2 (en) | 2019-02-05 | 2021-03-23 | Open Water Internet Inc. | Diffuse optical imaging with multiple beams |
US11320370B2 (en) | 2019-06-26 | 2022-05-03 | Open Water Internet Inc. | Apparatus for directing optical and acoustic signals |
US11846586B2 (en) | 2019-06-26 | 2023-12-19 | Open Water Internet, Inc. | Common axis for optical and acoustic signals |
US11581696B2 (en) | 2019-08-14 | 2023-02-14 | Open Water Internet Inc. | Multi-channel laser |
US11622686B2 (en) | 2019-11-22 | 2023-04-11 | Open Water Internet, Inc. | Optical imaging with unshifted reference beam |
US11819318B2 (en) | 2020-04-27 | 2023-11-21 | Open Water Internet Inc. | Optical imaging from light coherence |
US11559208B2 (en) | 2020-05-19 | 2023-01-24 | Open Water Internet Inc. | Imaging with scattering layer |
US11259706B2 (en) | 2020-05-19 | 2022-03-01 | Open Water Internet Inc. | Dual wavelength imaging and out of sample optical imaging |
Also Published As
Publication number | Publication date |
---|---|
WO2006031258A2 (en) | 2006-03-23 |
US7463362B2 (en) | 2008-12-09 |
US20090103102A1 (en) | 2009-04-23 |
US20050243322A1 (en) | 2005-11-03 |
WO2006031258A3 (en) | 2007-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US7728986B2 (en) | Systems and methods for digital detection of a tomographic signal | |
US9037216B2 (en) | Systems and methods for dynamic imaging of tissue using digital optical tomography | |
Khan et al. | FPGA-based voltage and current dual drive system for high frame rate electrical impedance tomography | |
US5774213A (en) | Techniques for measuring difference of an optical property at two wavelengths by modulating two sources to have opposite-phase components at a common frequency | |
Siegel et al. | Design and evaluation of a continuous-wave diffuse optical tomography system | |
Hartov et al. | A multichannel continuously selectable multifrequency electrical impedance spectroscopy measurement system | |
US7400919B2 (en) | Oximeter ambient light cancellation | |
US7215985B2 (en) | Oximeter cross-talk reduction | |
EP0689398A1 (en) | Pathlength corrected oximeter and the like | |
Lasker et al. | Digital-signal-processor-based dynamic imaging system for optical tomography | |
IL260889A (en) | Heterodyne-mimicking adapter | |
CN110768661B (en) | Phase-locked amplifier based on neural network | |
JP4112882B2 (en) | Biological light measurement device | |
Gervasoni et al. | FPGA-based lock-in amplifier with sub-ppm resolution working up to 6 MHz | |
Teren et al. | Direct comparison of analogue and digital FGPA-based approaches to synchronous detection | |
Pollastrone et al. | Fully-digital low-frequency lock-in amplifier for photoluminescence measurements | |
Lasker | A digital-signal-processor-based optical tomographic system for dynamic imaging of joint diseases | |
Masciotti et al. | Digital lock-in algorithm for biomedical spectroscopy and imaging instruments with multiple modulated sources | |
Gubello et al. | Modular 76-channel instrument for Broadband Raman Spectroscopy | |
Lasker et al. | Digital detection techniques for dynamic optical tomographic imaging | |
Lasker et al. | Dynamic optical tomographic imager with optimized digital lock-in filtering | |
SU1708298A1 (en) | Device for determining the bioobject hydration degree | |
Li et al. | Design of a digital optical tomography system for dynamic breast imaging | |
WO2015051016A1 (en) | System, method and computer-accessible medium for utilizing discrete fourier-transform for frequency near-infrared spectroscopy | |
KR101820798B1 (en) | Method and apparatus of Homodyne-based diffuse optical spectroscopic imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
AS | Assignment |
Owner name: NATIONAL INSTITUTES OF HEALTH (NIH), U.S. DEPT. OF Free format text: CONFIRMATORY LICENSE;ASSIGNOR:COLUMBIA UNIV NEW YORK MORNINGSIDE;REEL/FRAME:025121/0741 Effective date: 20091217 |
|
FPAY | Fee payment |
Year of fee payment: 4 |
|
SULP | Surcharge for late payment | ||
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.) |
|
FEPP | Fee payment procedure |
Free format text: 7.5 YR SURCHARGE - LATE PMT W/IN 6 MO, SMALL ENTITY (ORIGINAL EVENT CODE: M2555) |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2552) Year of fee payment: 8 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
FEPP | Fee payment procedure |
Free format text: 11.5 YR SURCHARGE- LATE PMT W/IN 6 MO, SMALL ENTITY (ORIGINAL EVENT CODE: M2556); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 12TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2553); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 12 |