WO2017091370A1 - Tag reading using targeted spatial spectral detection - Google Patents

Tag reading using targeted spatial spectral detection Download PDF

Info

Publication number
WO2017091370A1
WO2017091370A1 PCT/US2016/061680 US2016061680W WO2017091370A1 WO 2017091370 A1 WO2017091370 A1 WO 2017091370A1 US 2016061680 W US2016061680 W US 2016061680W WO 2017091370 A1 WO2017091370 A1 WO 2017091370A1
Authority
WO
WIPO (PCT)
Prior art keywords
interest
spectral
region
data
determining
Prior art date
Application number
PCT/US2016/061680
Other languages
French (fr)
Inventor
Timothy Learmonth
Mark Hsu
Denis Ivanov
Original Assignee
Trutag Technologies, Inc.
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Trutag Technologies, Inc. filed Critical Trutag Technologies, Inc.
Priority to CN201680062190.7A priority Critical patent/CN108139201A/en
Priority to EP16869064.2A priority patent/EP3380806A4/en
Priority to JP2018521231A priority patent/JP2019502097A/en
Publication of WO2017091370A1 publication Critical patent/WO2017091370A1/en

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2803Investigating the spectrum using photoelectric array detector
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/12Generating the spectrum; Monochromators
    • G01J3/26Generating the spectrum; Monochromators using multiple reflection, e.g. Fabry-Perot interferometer, variable interference filters
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2823Imaging spectrometer
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2803Investigating the spectrum using photoelectric array detector
    • G01J2003/28132D-array
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J3/2823Imaging spectrometer
    • G01J2003/2826Multispectral imaging, e.g. filter imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J2003/283Investigating the spectrum computer-interfaced
    • G01J2003/2833Investigating the spectrum computer-interfaced and memorised spectra collection
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J2003/283Investigating the spectrum computer-interfaced
    • G01J2003/284Spectral construction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J2003/2859Peak detecting in spectrum
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J2003/2866Markers; Calibrating of scan
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • G01J2003/2866Markers; Calibrating of scan
    • G01J2003/2879Calibrating scan, e.g. Fabry Perot interferometer
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/10Terrestrial scenes
    • G06V20/194Terrestrial scenes using hyperspectral data, i.e. more or other wavelengths than RGB

Definitions

  • hyperspectral imaging is used for object detection, such as diseased crops, military targets, or geological formations. In those cases once the objects are detected, the job is done. So the number of frames in the hyperspectral cube is set at just enough to ensure that the detection job can be done. However for determining spectra, location and decoding must be achieved. So, just detection is not sufficient. Location and decoding require a much larger amount of data than detection, creating a problem where the storage and processing capabilities of the reader are challenged.
  • Figure 1 is a diagram illustrating an embodiment of an optical setup for spectral response detection.
  • Figure 2 is a diagram illustrating an embodiment of a Fabry -Perot interferometer.
  • Figure 3 is a diagram illustrating an embodiment of a spectral data cube.
  • Figure 4 is a diagram illustrating an embodiment of interferometer transmissivity.
  • Figure 5 is a diagram illustrating an embodiment of interferometer transmission peak wavelength locations.
  • Figure 6A is a diagram illustrating an embodiment of vectors in a parameter space.
  • Figure 6B is a diagram illustrating an embodiment of an angular measurement as a function of interferometer gap.
  • Figure 7 is a flow diagram illustrating an embodiment of a process for spectral reading using targeted spatial spectral detection.
  • Figure 8 is a flow diagram illustrating an embodiment of a process for determining a region of interest using a sample set of data.
  • Figure 9 is a flow diagram illustrating an embodiment of a process for determining a spectral peak for a region of interest.
  • the invention can be implemented in numerous ways, including as a process; an apparatus; a system; a composition of matter; a computer program product embodied on a computer readable storage medium; and/or a processor, such as a processor configured to execute instructions stored on and/or provided by a memory coupled to the processor.
  • these implementations, or any other form that the invention may take, may be referred to as techniques.
  • the order of the steps of disclosed processes may be altered within the scope of the invention.
  • a component such as a processor or a memory described as being configured to perform a task may be implemented as a general component that is temporarily configured to perform the task at a given time or a specific component that is manufactured to perform the task.
  • the term 'processor' refers to one or more devices, circuits, and/or processing cores configured to process data, such as computer program instructions.
  • a system for determining a spectrum comprises an interface and a processor.
  • the interface is configured to receive a sample set of intensity data for an array of spatial locations and a set of spectral configurations, where each spectral configuration detects a different wavelength or combination of wavelengths of light.
  • the processor is configured to determine a region of interest using the sample set of data and determine a spectral peak for the region of interest.
  • both finding an object and decoding a reflected spectrum is required.
  • Memory, time, and/or processing power can be saved by splitting the operations of detection and extraction/decoding.
  • memory can be saved by processing frames as they appear without saving them, so that the full set of frames required for spectral detection doesn't have to be kept in memory.
  • objects e.g., tags
  • objects can be found by looking for local maxima in the standard deviation of received intensity values along the spectral axis.
  • objects e.g., tags
  • filters can be found by looking for peaks by looking for maxima after filtering input data (e.g., running a convolutional filter - for example, a matched filter with suppression in the surround).
  • filtering and looking for a maxima in standard deviation are done in sequence to identify a region of interest.
  • filtering and looking for peaks in two ranges of frequencies are done in sequence to identify a region of interest.
  • each frame is a digitized image of a sample collected at a specific wavelength or set of wavelengths, and different frames are collected using different wavelengths or set of wavelengths.
  • the sum of the frames and the sum of the square of the frames are needed. So when a frame appears, each intensity value in each location of the frame and its square are each added to an accumulator. Each frame is then overwritten by the next frame. Once all the frames are collected, a simple operation is performed on the two accumulated values (e.g., the sum and the sum of the squares) for each location in the frame to get the standard deviation of the value in that location of the frame.
  • a new frame e.g., a map
  • a new frame e.g., a map
  • the hypercube is captured or to a subset of a hypercube.
  • spectral detection is achieved using filtering along the spectral axis, assuming that the filter kernel is shorter than the full spectral axis.
  • filtering With filtering, a subset of frames equal to the kernel length is kept in memory at one time, rather than just the one for the standard deviation. This is what we're calling sequential processing.
  • the kernel length used for this step we find the minimum number of frames sufficient to detect an object reliably by isolating a spectral feature in the filtering approach. For example, in the event that the filter kernel is tuned to couple strongly to a spectral feature with spectral width equal to 1/10 the length of the full spectral axis, the kernel length might be 1/10 the full spectral axis length.
  • the targeted spectral features might be a set of strong peaks or valleys in the spectrum that serve to increase the standard deviation along the spectral axis.
  • the frames in this step could be stored in memory and reused in the last step below. In that case processing power is saved by not processing the entire hyperspectral cube, and time is saved, by computing the regions of interest while capture of the hyperspectral cube is still ongoing. At the end of the sequential processing, a 2D array has been determined over the spatial axes that has a higher value where candidate objects are present, and a lower value where they are not.
  • regions of interest are chosen that are likely to contain an object. For example, regions of interest may be chosen by searching the 2D array for maxima (points where all neighbors are less than the point itself). Maxima may then be sorted in descending order of value, creating a list of points, with the first in the list most likely to correspond to a tag, the second in the list the second most likely, etc. Regions of interest can then be determined by including a fixed range of points around each point in the list, with the size of the range determined from the expected size of an object. The regions of interest are then the focus for analysis.
  • a full scan is run to extract spectral information with enough frames to extract the object information. But during the data collection, only the hyperspectral information in the regions of interest is saved. This greatly reduces the amount of data saved in memory. For example, in the event that the objects comprise 1% of the reader field of view, a data reduction on the order of 99% can be achieved in this manner.
  • the data collected is not a computed spectrum, so further processing is necessary to obtain a spectrum from the data collected.
  • the reader is based on a Fabry -Perot Interferometer in a low-finesse Fourier Transform mode, for example, then conversion from captured data to a spectrum requires a Fourier transform. If only the regions of interest are processed into spectra, the computational effort required to decode a spectrum of an object is further reduced. Once a spectrum from an object such as a tag is obtained, decoding can proceed.
  • a tag comprises a reflector with selective reflection.
  • the reflector comprises one or more of the following: a rugate tag, a Bragg reflector, or any other appropriate reflector.
  • tags comprise one of the following materials: silicon, silicon dioxide, silicon nitride, doped silicon, or any other appropriate material.
  • the unique optical signature of each tag can be read using an absolute or a relative spectral measurement device, apparatus, or system.
  • tags comprise the surface of a silicon wafer that is etched to have a spectral code encoded by the etching. A thin layer from the surface of the etched wafer is removed and divided into small tags, and the resultant tags contain a complex porous nanostructure that is programmed during electrochemical synthesis to display a unique reflectivity spectrum.
  • the tags are then oxidized by a high-temperature bake step to turn the crystalline, nanoporous silicon tags into amorphous, nanoporous silica.
  • This bake step stabilizes the nanoporous structure against further oxidation (thus stabilizing the spectral signature) and provides for the tags to be characterized as a GRAS excipient.
  • spectral configuration detects a different wavelength or combination of wavelengths of light by spectral filtering reflected light from an object before detecting the light in a detector.
  • the spectral filtering comprises a Fabry- Perot filter, an interference filter, or any other appropriate spectral filter.
  • the spectral filter is tunable.
  • the spectral filter is a set of fixed filters that are in front of separate detectors or are mechanically swapped or optically switched in front of a fixed detector.
  • the Fabry-Perot filter is mechanically tunable, electro-optically tunable (e.g., a change in index of refraction of the medium in between the partially reflecting mirrors), acousto-optically tunable, or any other appropriate filter.
  • Figure 1 is a diagram illustrating an embodiment of an optical setup for spectral response detection.
  • Figure 1 comprises measurement area 100 comprising one or more regions for producing optical spectra (e.g., region 102). Measurement area 100 is illuminated by light 104.
  • Reflected light 106 comprises light reflected by region 102. Reflected light 106 is focused by lens 108, passes through filter 110, and is captured by optical detector 112.
  • Filter 110 comprises a filter for transmitting light at some frequencies and not at other frequencies.
  • Optical detector 112 comprises an optical detector for detecting optical intensities.
  • Spectral response detector control system 114 comprises a spectral response detector control system for interacting with optical detector 112.
  • Spectral response detector control system 114 provides control information to optical detector 112 (e.g., control information indicating to capture data) and receives optical intensity data from optical detector 112.
  • regions for producing optical spectra comprise optical tags (e.g., optical tags engineered to reflect light with a recognizable spectrum).
  • light 104 comprises broadband light, narrowband light, filtered light, light from a light-emitting diode, laser light, or any other appropriate light.
  • light 104 is incident to measurement area 100 from a single point, at a single angle, from multiple angles, or incident in any other appropriate manner.
  • filter 110 comprises a tunable optical filter.
  • filter 110 comprises a Fabry -Perot interferometer.
  • optical detector 112 comprises an array of optical detector pixels for detecting an array of optical intensities.
  • the array of optical detector pixels comprises an x axis and a y axis.
  • each pixel of the array of optical detector pixels comprises a set of optical detectors, each optical detector comprising a color filter (e.g., each pixel comprises three detectors, a first detector comprising a red color filter, a second detector comprising a green color filter, and a third detector comprising a blue color filter).
  • spectral response detector control system 114 interacts with filter 110 (e.g., to indicate tunable filter properties).
  • spectral response detector control system 114 interacts with filter 110 to indicate a Fabry -Perot interferometer gap size.
  • spectral response detector control system 114 is configured to receive a calibration data.
  • the calibration data comprises a set of intensity data for an array of spatial locations and a range of spectral configurations for a monochromatic source (e.g., a set of gap sizes for a Fabry-Perot interferometer).
  • FIG. 2 is a diagram illustrating an embodiment of a Fabry-Perot interferometer.
  • Fabry-Perot interferometer 200 comprises filter 110 of Figure 1.
  • Fabry-Perot interferometer 200 comprises mirror 202 and mirror 204.
  • Each of mirror 202 and mirror 204 comprises a partially reflective mirror (e.g., some light is able to pass through and some light is reflected).
  • Each of mirror 202 and mirror 204 comprises a mirror silvering on one side (e.g., mirror 202 comprises a mirror silvering on its right side as shown in Figure 2 and mirror 204 comprises a mirror silvering on its left side as shown in Figure 2).
  • Piezo element 206 and piezo element 208 comprise piezo elements for changing size. Piezo element 206 and piezo element 208 change size according to an applied voltage. When piezo element 206 and piezo element 208 change size, the gap between the mirror silvering of mirror 202 and the mirror silvering of mirror 208 changes. Changing the gap between the mirror silvering of mirror 202 and the mirror silvering of mirror 208 causes the optical properties (e.g., the light transmission and reflectance properties) of Fabry-Perot interferometer 200 to change. In this way, Fabry-Perot interferometer 200 comprises a tunable optical filter.
  • mirror silvering comprises a partially reflective metal layer
  • mirror 202 and mirror 204 are held together by 1, 2, 3, 4, 5, or any other appropriate number of piezo elements.
  • the spatial positioning of the piezo elements separating mirror 202 and mirror 204 enable adjustment of the angle between mirror 202 and mirror 204.
  • the mirrors are separated with a fixed distance and a medium interior to the mirrors changes its index of refraction (e.g., electro- or acousto- optically) and instead of a set of gap sizes there is a set of different index-changed path lengths.
  • a medium interior to the mirrors changes its index of refraction (e.g., electro- or acousto- optically) and instead of a set of gap sizes there is a set of different index-changed path lengths.
  • FIG. 3 is a diagram illustrating an embodiment of a spectral data cube.
  • spectral data cube 300 of Figure 3 comprises a cube representing data measured by an optical detector (e.g., optical detector 112 of Figure 1).
  • Spectral data cube 300 comprises a cube of data recorded by an optical detector.
  • the x axis of spectral data cube 300 corresponds to the x axis of the optical detector and the y axis of spectral data cube 300 corresponds to the y axis of the optical detector.
  • the z axis of spectral data cube 300 corresponds to an interferometer gap size (e.g., the mirror gap size of a Fabry -Perot interferometer).
  • the data intensity shown at a given data location indicates the intensity of the light received by the optical detector.
  • the interferometer gap size comprises the gap size of Fabry -
  • spectral data cube 300 comprises a spectral data hypercube with one or more of the following properties: each data location comprising three light intensity measurements, each measurement taken by a separate optical detector at the given location, each optical detector associated with a different color filter, or any other appropriate properties. In some embodiments, multiple peaks associated with different interferometer gap sizes are measured at a given detector location.
  • Figure 4 is a diagram illustrating an embodiment of interferometer transmissivity.
  • curve 400 of Figure 4 illustrates the light transmissivity of a Fabry -Perot interferometer (e.g., Fabry -Perot interferometer 200 of Figure 2) vs. light wavelength.
  • curve 400 of Figure 4 illustrates the light transmissivity of a Fabry -Perot interferometer vs. light wavelength for a given gap size. Multiple light transmission peaks are seen for the given gap size. Three light transmission peaks are seen in curve 400.
  • Fabry -Perot interferometer e.g., Fabry -Perot interferometer 200 of Figure 2
  • curve 400 of Figure 4 illustrates the light transmissivity of a Fabry -Perot interferometer vs. light wavelength for a given gap size. Multiple light transmission peaks are seen for the given gap size. Three light transmission peaks are seen in curve 400.
  • FIG. 5 is a diagram illustrating an embodiment of interferometer transmission peak wavelength locations.
  • set of curves 500 illustrates a set of transmission peak locations for a Fabry-Perot interferometer (e.g., Fabry-Perot interferometer 200 of Figure 2).
  • set of curves 500 illustrates the location of each of a set of interferometer transmission peak locations as the interferometer gap of the Fabry-Perot interferometer is changed.
  • set of curves 500 comprises 4 curves. For a given gap size (e.g., gap size gi), there are a set of transmission peaks at different wavelengths (e.g., as shown by curve 400 of Figure
  • Figure 6A is a diagram illustrating an embodiment of vectors in a parameter space.
  • the vectors of Figure 6 A are two dimensional representations of intensity vectors (e.g., vectors comprising a first color intensity in a first direction as measured by a first color detector and a second color intensity in a second direction as measured by second color detector).
  • a vector is calculated and displayed in a three dimensional space (e.g., the three dimensions correspond to three color detectors).
  • vectors ⁇ and ⁇ 2 comprise calibration vectors for a Fabry-Perot interferometer (e.g., Fabry-Perot
  • a set of calibration vectors is associated with each gap size of the Fabry-Perot interferometer.
  • Each calibration vector of the set of calibration vectors is associated with a transmission peak associated with the gap size.
  • Determining a calibration vector comprises illuminating the Fabry-Perot interferometer with monochromatic light of the transmission peak wavelength and recording the intensity measured by each optical detector of a set of optical detectors forming a pixel (e.g., each optical detector associated with a different color filter).
  • Each optical detector is associated with an axis of the parameter space.
  • Measured intensity vector 600 comprises a vector associated with an optical measurement.
  • measured intensity vector 600 comprises a representation of data at one point within data cube 300 of Figure 3.
  • Angles ⁇ and ⁇ 2 comprise the angles formed between the calibration vectors and measured intensity vector 600.
  • intensity vectors comprise intensity vectors in a two- dimensional parameter space (e.g., representing two color intensities), intensity vectors in a three- dimensional parameter space (e.g., representing three color intensities), intensity vectors in a five- dimensional parameter space (e.g., representing five color intensities), or intensity vectors in any other appropriate parameter space.
  • Figure 6B is a diagram illustrating an embodiment of an angular measurement as a function of interferometer gap. In the example shown, ⁇ and ⁇ 2 are plotted as an interferometer gap is varied. ⁇ 2 is seen to dip at a gap size of interest, while ⁇ is seen to rise, indicating the measured intensity vector swings toward calibration vector ⁇ 2 . This swing indicates that the measured light is passing through the Fabry -Perot interferometer at the transmission peak associated with calibration vector ⁇ 2 and not the transmission peak associated with calibration vector ⁇ . In some
  • ⁇ and ⁇ 2 comprise ⁇ and ⁇ 2 as shown in Figure 6A.
  • angles are measured in three dimensions towards calibration vectors.
  • the angles shift as a function of gap towards one of the calibration vectors as the gap is incremented or decremented, it is determined that the peak corresponds to the peak associated with the calibration vector.
  • FIG. 7 is a flow diagram illustrating an embodiment of a process for spectral reading using targeted spatial spectral detection.
  • a sample set of intensity data is received for an array of spatial locations and a range of spectral configurations.
  • data in the array of spatial locations is transformed for rotation, translation, and/or key stoning.
  • a region of interest is determined using the sample set of data.
  • a spectral peak is determined for the region of interest.
  • the process of Figure 7 is executed by spectral response detector control system 114.
  • a set of regions of interest is determined in 702.
  • a spectral peak is determined in 704 for all regions of interest in the set of regions of interest.
  • more than one spectral peak is determined for a region of interest.
  • determining one or more spectral peaks comprises decoding a signature.
  • the signature comprises a tag signature.
  • spectral peak is identified within the spatial array of data associated with one spectral configuration.
  • the data of the spatial array is corrected for rotation, translation, and/or keystoning (e.g., trapezoidal scaling).
  • spatial array data is processed in an area larger than a region of interest to make sure that a peak that is associated with one spectral configuration can be matched to the peak associated with another spectral configuration. For example, a first spatial array of data is taken associated with one spectral configuration and then at a later time a spatial array of data is taken associated with another spectral configuration. Between the two times, it is possible that the data taking device is rotated, translated, and/or tilted making the data in the two corresponding spatial arrays of data related to each other with a rotation, translation, and/or key stoning translation.
  • Figure 8 is a flow diagram illustrating an embodiment of a process for determining a region of interest using a sample set of data.
  • the flow diagram of Figure 8 illustrates an
  • a data frame associated with a spectral configuration is received.
  • data frame is transformed for rotation, translation, and/or key stoning. For example, to get correspondence between spatial arrays of data associated with different spectral configurations, transformation is performed if necessary between spatial arrays to make sure that the data in the arrays are spatially associated with each other. In some embodiments, fiducials in the data frames are used to determine the transformations.
  • the data frame is added to a value accumulator. In the event the data frame comprises the first data frame, the value accumulator is set to the value of the data frame.
  • the data frame values are squared (e.g., each value of the data frame values is squared).
  • the squared data frame values are added to a squared value accumulator. In the event the data frame comprises the first data frame, the squared value accumulator is set to the value of the squared data frame.
  • the standard deviation is determined. Determining the standard deviation for a set of data frames comprises determining the standard deviation for each point of the data frame, across the set of data frames.
  • the standard deviation for the set of data frames can be determined by first determining the data mean by dividing the values of the value accumulator by the total number of data frames received.
  • a mean of squared data can additionally be determined by dividing the values of the squared value accumulator by the total number of data frames received.
  • the standard deviation for the set of data frames can then be determined by determining the square root of the difference between the mean of squared data and the square of the data mean.
  • one or more regions of interest are identified by determining local maxima in the standard deviation data (e.g., local maxima across the data frame).
  • a data frame comprises a two-dimensional array of data pixels.
  • each data pixel comprises three colored pixels, each colored pixel associated with a color filter of a different color.
  • one or more signal processing techniques other than standard deviation computation are used for processing the spectral axis of the data cube (e.g., thresholding, filtering, matched filtering, etc.).
  • Figure 9 is a flow diagram illustrating an embodiment of a process for determining a spectral peak for a region of interest.
  • the process of Figure 9 implements 704 of Figure 7.
  • an indication to capture further intensity data for one or more regions of interest is provided.
  • the further intensity data for the one or more regions of interest is received.
  • the intensity data is transformed into a parameter space (e.g., the parameter space shown in Figure 6A).
  • the difference between the intensity data in the parameter space and calibration data in the parameter space is determined.
  • the spectral peak is determined (e.g., by determining a calibration data vector in the parameter space that the intensity data in the parameter space exhibits a movement towards).

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Spectrometry And Color Measurement (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)
  • Solid State Image Pick-Up Elements (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

A system for determining a spectrum includes an interface and a processor. The interface is configured to receive a sample set of intensity data for an array of spatial locations and a set of spectral configurations. The processor is configured to determine a region of interest using the sample set of intensity data and determine a spectral peak for the region of interest.

Description

TAG READING USING TARGETED SPATIAL SPECTRAL DETECTION
CROSS REFERENCE TO OTHER APPLICATIONS
[0001] This application claims priority to U.S. Provisional Patent Application No.
62/259,244 entitled TAG READING USING TARGETED SPATIAL SPECTRAL DETECTION filed November 24, 2015 which is incorporated herein by reference for all purposes.
BACKGROUND OF THE INVENTION
[0002] Normally hyperspectral imaging is used for object detection, such as diseased crops, military targets, or geological formations. In those cases once the objects are detected, the job is done. So the number of frames in the hyperspectral cube is set at just enough to ensure that the detection job can be done. However for determining spectra, location and decoding must be achieved. So, just detection is not sufficient. Location and decoding require a much larger amount of data than detection, creating a problem where the storage and processing capabilities of the reader are challenged.
BRIEF DESCRIPTION OF THE DRAWINGS
[0003] Various embodiments of the invention are disclosed in the following detailed description and the accompanying drawings.
[0004] Figure 1 is a diagram illustrating an embodiment of an optical setup for spectral response detection.
[0005] Figure 2 is a diagram illustrating an embodiment of a Fabry -Perot interferometer.
[0006] Figure 3 is a diagram illustrating an embodiment of a spectral data cube.
[0007] Figure 4 is a diagram illustrating an embodiment of interferometer transmissivity.
[0008] Figure 5 is a diagram illustrating an embodiment of interferometer transmission peak wavelength locations.
[0009] Figure 6A is a diagram illustrating an embodiment of vectors in a parameter space. [0010] Figure 6B is a diagram illustrating an embodiment of an angular measurement as a function of interferometer gap.
[0011] Figure 7 is a flow diagram illustrating an embodiment of a process for spectral reading using targeted spatial spectral detection.
[0012] Figure 8 is a flow diagram illustrating an embodiment of a process for determining a region of interest using a sample set of data.
[0013] Figure 9 is a flow diagram illustrating an embodiment of a process for determining a spectral peak for a region of interest.
DETAILED DESCRIPTION
[0014] The invention can be implemented in numerous ways, including as a process; an apparatus; a system; a composition of matter; a computer program product embodied on a computer readable storage medium; and/or a processor, such as a processor configured to execute instructions stored on and/or provided by a memory coupled to the processor. In this specification, these implementations, or any other form that the invention may take, may be referred to as techniques. In general, the order of the steps of disclosed processes may be altered within the scope of the invention. Unless stated otherwise, a component such as a processor or a memory described as being configured to perform a task may be implemented as a general component that is temporarily configured to perform the task at a given time or a specific component that is manufactured to perform the task. As used herein, the term 'processor' refers to one or more devices, circuits, and/or processing cores configured to process data, such as computer program instructions.
[0015] A detailed description of one or more embodiments of the invention is provided below along with accompanying figures that illustrate the principles of the invention. The invention is described in connection with such embodiments, but the invention is not limited to any embodiment. The scope of the invention is limited only by the claims and the invention
encompasses numerous alternatives, modifications and equivalents. Numerous specific details are set forth in the following description in order to provide a thorough understanding of the invention. These details are provided for the purpose of example and the invention may be practiced according to the claims without some or all of these specific details. For the purpose of clarity, technical material that is known in the technical fields related to the invention has not been described in detail so that the invention is not unnecessarily obscured. [0016] A system for determining a spectrum comprises an interface and a processor. The interface is configured to receive a sample set of intensity data for an array of spatial locations and a set of spectral configurations, where each spectral configuration detects a different wavelength or combination of wavelengths of light. The processor is configured to determine a region of interest using the sample set of data and determine a spectral peak for the region of interest.
[0017] In some embodiments, for spectral decoding, both finding an object and decoding a reflected spectrum is required. Memory, time, and/or processing power can be saved by splitting the operations of detection and extraction/decoding.
[0018] In some embodiments, for determining a spectrum, memory can be saved by processing frames as they appear without saving them, so that the full set of frames required for spectral detection doesn't have to be kept in memory. For example, objects (e.g., tags) can be found by looking for local maxima in the standard deviation of received intensity values along the spectral axis. In another example, objects (e.g., tags) can be found by looking for peaks by looking for maxima after filtering input data (e.g., running a convolutional filter - for example, a matched filter with suppression in the surround). In some embodiments, filtering and looking for a maxima in standard deviation are done in sequence to identify a region of interest. In some embodiments, filtering and looking for peaks in two ranges of frequencies are done in sequence to identify a region of interest. In this case each frame is a digitized image of a sample collected at a specific wavelength or set of wavelengths, and different frames are collected using different wavelengths or set of wavelengths. To calculate the standard deviation, the sum of the frames and the sum of the square of the frames are needed. So when a frame appears, each intensity value in each location of the frame and its square are each added to an accumulator. Each frame is then overwritten by the next frame. Once all the frames are collected, a simple operation is performed on the two accumulated values (e.g., the sum and the sum of the squares) for each location in the frame to get the standard deviation of the value in that location of the frame. In this manner, a new frame (e.g., a map) is created that has the same spatial dimensions as the original frames but contains the standard deviation through the frames at each point. Then local maxima in the map correspond to locations of candidate tags on the full hyperspectral cube. In various embodiments, this method is applied either to a full hypercube as the hypercube is captured or to a subset of a hypercube.
[0019] In some embodiments, spectral detection is achieved using filtering along the spectral axis, assuming that the filter kernel is shorter than the full spectral axis. With filtering, a subset of frames equal to the kernel length is kept in memory at one time, rather than just the one for the standard deviation. This is what we're calling sequential processing. To determine the kernel length used for this step, we find the minimum number of frames sufficient to detect an object reliably by isolating a spectral feature in the filtering approach. For example, in the event that the filter kernel is tuned to couple strongly to a spectral feature with spectral width equal to 1/10 the length of the full spectral axis, the kernel length might be 1/10 the full spectral axis length. Or it could be done by capturing a subset of many spectral features in the standard deviation approach. For example, the targeted spectral features might be a set of strong peaks or valleys in the spectrum that serve to increase the standard deviation along the spectral axis. Also in the event that memory isn't a constraint, the frames in this step could be stored in memory and reused in the last step below. In that case processing power is saved by not processing the entire hyperspectral cube, and time is saved, by computing the regions of interest while capture of the hyperspectral cube is still ongoing. At the end of the sequential processing, a 2D array has been determined over the spatial axes that has a higher value where candidate objects are present, and a lower value where they are not. From this array, regions of interest are chosen that are likely to contain an object. For example, regions of interest may be chosen by searching the 2D array for maxima (points where all neighbors are less than the point itself). Maxima may then be sorted in descending order of value, creating a list of points, with the first in the list most likely to correspond to a tag, the second in the list the second most likely, etc. Regions of interest can then be determined by including a fixed range of points around each point in the list, with the size of the range determined from the expected size of an object. The regions of interest are then the focus for analysis.
[0020] In some embodiments, at this point a full scan is run to extract spectral information with enough frames to extract the object information. But during the data collection, only the hyperspectral information in the regions of interest is saved. This greatly reduces the amount of data saved in memory. For example, in the event that the objects comprise 1% of the reader field of view, a data reduction on the order of 99% can be achieved in this manner.
[0021] In some embodiments, the data collected is not a computed spectrum, so further processing is necessary to obtain a spectrum from the data collected. In the event that the reader is based on a Fabry -Perot Interferometer in a low-finesse Fourier Transform mode, for example, then conversion from captured data to a spectrum requires a Fourier transform. If only the regions of interest are processed into spectra, the computational effort required to decode a spectrum of an object is further reduced. Once a spectrum from an object such as a tag is obtained, decoding can proceed. In some embodiments, a tag comprises a reflector with selective reflection. In various embodiments, the reflector comprises one or more of the following: a rugate tag, a Bragg reflector, or any other appropriate reflector. In various embodiments, tags comprise one of the following materials: silicon, silicon dioxide, silicon nitride, doped silicon, or any other appropriate material. In some embodiments, the unique optical signature of each tag can be read using an absolute or a relative spectral measurement device, apparatus, or system. In some embodiments, tags comprise the surface of a silicon wafer that is etched to have a spectral code encoded by the etching. A thin layer from the surface of the etched wafer is removed and divided into small tags, and the resultant tags contain a complex porous nanostructure that is programmed during electrochemical synthesis to display a unique reflectivity spectrum. The tags are then oxidized by a high-temperature bake step to turn the crystalline, nanoporous silicon tags into amorphous, nanoporous silica. This bake step stabilizes the nanoporous structure against further oxidation (thus stabilizing the spectral signature) and provides for the tags to be characterized as a GRAS excipient.
[0022] In some embodiments, spectral configuration detects a different wavelength or combination of wavelengths of light by spectral filtering reflected light from an object before detecting the light in a detector. In various embodiments, the spectral filtering comprises a Fabry- Perot filter, an interference filter, or any other appropriate spectral filter. In some embodiments, the spectral filter is tunable. In some embodiments, the spectral filter is a set of fixed filters that are in front of separate detectors or are mechanically swapped or optically switched in front of a fixed detector. In various embodiments, the Fabry-Perot filter is mechanically tunable, electro-optically tunable (e.g., a change in index of refraction of the medium in between the partially reflecting mirrors), acousto-optically tunable, or any other appropriate filter.
[0023] Figure 1 is a diagram illustrating an embodiment of an optical setup for spectral response detection. Figure 1 comprises measurement area 100 comprising one or more regions for producing optical spectra (e.g., region 102). Measurement area 100 is illuminated by light 104. Reflected light 106 comprises light reflected by region 102. Reflected light 106 is focused by lens 108, passes through filter 110, and is captured by optical detector 112. Filter 110 comprises a filter for transmitting light at some frequencies and not at other frequencies. Optical detector 112 comprises an optical detector for detecting optical intensities. Spectral response detector control system 114 comprises a spectral response detector control system for interacting with optical detector 112. Spectral response detector control system 114 provides control information to optical detector 112 (e.g., control information indicating to capture data) and receives optical intensity data from optical detector 112.
[0024] In some embodiments, regions for producing optical spectra (e.g., region 102) comprise optical tags (e.g., optical tags engineered to reflect light with a recognizable spectrum). In various embodiments, light 104 comprises broadband light, narrowband light, filtered light, light from a light-emitting diode, laser light, or any other appropriate light. In various embodiments, light 104 is incident to measurement area 100 from a single point, at a single angle, from multiple angles, or incident in any other appropriate manner. In some embodiments, filter 110 comprises a tunable optical filter. In some embodiments, filter 110 comprises a Fabry -Perot interferometer. In some embodiments, optical detector 112 comprises an array of optical detector pixels for detecting an array of optical intensities. In some embodiments, the array of optical detector pixels comprises an x axis and a y axis. In some embodiments, each pixel of the array of optical detector pixels comprises a set of optical detectors, each optical detector comprising a color filter (e.g., each pixel comprises three detectors, a first detector comprising a red color filter, a second detector comprising a green color filter, and a third detector comprising a blue color filter). In some embodiments, spectral response detector control system 114 interacts with filter 110 (e.g., to indicate tunable filter properties). In some embodiments, spectral response detector control system 114 interacts with filter 110 to indicate a Fabry -Perot interferometer gap size. In some
embodiments, spectral response detector control system 114 is configured to receive a calibration data. In some embodiments, the calibration data comprises a set of intensity data for an array of spatial locations and a range of spectral configurations for a monochromatic source (e.g., a set of gap sizes for a Fabry-Perot interferometer).
[0025] Figure 2 is a diagram illustrating an embodiment of a Fabry-Perot interferometer. In some embodiments, Fabry-Perot interferometer 200 comprises filter 110 of Figure 1. Fabry-Perot interferometer 200 comprises mirror 202 and mirror 204. Each of mirror 202 and mirror 204 comprises a partially reflective mirror (e.g., some light is able to pass through and some light is reflected). Each of mirror 202 and mirror 204 comprises a mirror silvering on one side (e.g., mirror 202 comprises a mirror silvering on its right side as shown in Figure 2 and mirror 204 comprises a mirror silvering on its left side as shown in Figure 2). Mirror 202 and mirror 204 are held together by piezo element 206 and piezo element 208. Piezo element 206 and piezo element 208 comprise piezo elements for changing size. Piezo element 206 and piezo element 208 change size according to an applied voltage. When piezo element 206 and piezo element 208 change size, the gap between the mirror silvering of mirror 202 and the mirror silvering of mirror 208 changes. Changing the gap between the mirror silvering of mirror 202 and the mirror silvering of mirror 208 causes the optical properties (e.g., the light transmission and reflectance properties) of Fabry-Perot interferometer 200 to change. In this way, Fabry-Perot interferometer 200 comprises a tunable optical filter.
[0026] In some embodiments mirror silvering comprises a partially reflective metal layer
(e.g., a silver layer, an aluminum layer, a titanium layer, etc.). In various embodiments, mirror 202 and mirror 204 are held together by 1, 2, 3, 4, 5, or any other appropriate number of piezo elements. In some embodiments, the spatial positioning of the piezo elements separating mirror 202 and mirror 204 enable adjustment of the angle between mirror 202 and mirror 204.
[0027] In some embodiments, the mirrors are separated with a fixed distance and a medium interior to the mirrors changes its index of refraction (e.g., electro- or acousto- optically) and instead of a set of gap sizes there is a set of different index-changed path lengths.
[0028] Figure 3 is a diagram illustrating an embodiment of a spectral data cube. In some embodiments, spectral data cube 300 of Figure 3 comprises a cube representing data measured by an optical detector (e.g., optical detector 112 of Figure 1). Spectral data cube 300 comprises a cube of data recorded by an optical detector. The x axis of spectral data cube 300 corresponds to the x axis of the optical detector and the y axis of spectral data cube 300 corresponds to the y axis of the optical detector. The z axis of spectral data cube 300 corresponds to an interferometer gap size (e.g., the mirror gap size of a Fabry -Perot interferometer). The data intensity shown at a given data location (e.g., corresponding to a given x and y location on the optical detector and a given interferometer gap size) indicates the intensity of the light received by the optical detector.
[0029] In some embodiments, the interferometer gap size comprises the gap size of Fabry -
Perot interferometer 200 of Figure 2. In various embodiments, spectral data cube 300 comprises a spectral data hypercube with one or more of the following properties: each data location comprising three light intensity measurements, each measurement taken by a separate optical detector at the given location, each optical detector associated with a different color filter, or any other appropriate properties. In some embodiments, multiple peaks associated with different interferometer gap sizes are measured at a given detector location.
[0030] Figure 4 is a diagram illustrating an embodiment of interferometer transmissivity. In some embodiments, curve 400 of Figure 4 illustrates the light transmissivity of a Fabry -Perot interferometer (e.g., Fabry -Perot interferometer 200 of Figure 2) vs. light wavelength. In the example shown, curve 400 of Figure 4 illustrates the light transmissivity of a Fabry -Perot interferometer vs. light wavelength for a given gap size. Multiple light transmission peaks are seen for the given gap size. Three light transmission peaks are seen in curve 400. In various
embodiments, two light transmission peaks are seen as a function of gap size, three light transmission peaks are seen as a function of gap size, four light transmission peaks are seen as a function of gap size, or any other appropriate number of light transmission peaks are seen as a function of gap size. [0031] Figure 5 is a diagram illustrating an embodiment of interferometer transmission peak wavelength locations. In some embodiments, set of curves 500 illustrates a set of transmission peak locations for a Fabry-Perot interferometer (e.g., Fabry-Perot interferometer 200 of Figure 2). In the example shown, set of curves 500 illustrates the location of each of a set of interferometer transmission peak locations as the interferometer gap of the Fabry-Perot interferometer is changed. In the example shown, set of curves 500 comprises 4 curves. For a given gap size (e.g., gap size gi), there are a set of transmission peaks at different wavelengths (e.g., as shown by curve 400 of Figure
4).
[0032] Figure 6A is a diagram illustrating an embodiment of vectors in a parameter space.
In some embodiments, the vectors of Figure 6 A are two dimensional representations of intensity vectors (e.g., vectors comprising a first color intensity in a first direction as measured by a first color detector and a second color intensity in a second direction as measured by second color detector). In some embodiments, a vector is calculated and displayed in a three dimensional space (e.g., the three dimensions correspond to three color detectors). In the example shown, vectors λι and λ2 comprise calibration vectors for a Fabry-Perot interferometer (e.g., Fabry-Perot
interferometer 200 of Figure 2). A set of calibration vectors is associated with each gap size of the Fabry-Perot interferometer. Each calibration vector of the set of calibration vectors is associated with a transmission peak associated with the gap size. Determining a calibration vector comprises illuminating the Fabry-Perot interferometer with monochromatic light of the transmission peak wavelength and recording the intensity measured by each optical detector of a set of optical detectors forming a pixel (e.g., each optical detector associated with a different color filter). Each optical detector is associated with an axis of the parameter space. Measured intensity vector 600 comprises a vector associated with an optical measurement. In some embodiments, measured intensity vector 600 comprises a representation of data at one point within data cube 300 of Figure 3. Angles θι and θ2 comprise the angles formed between the calibration vectors and measured intensity vector 600.
[0033] In various embodiments, intensity vectors comprise intensity vectors in a two- dimensional parameter space (e.g., representing two color intensities), intensity vectors in a three- dimensional parameter space (e.g., representing three color intensities), intensity vectors in a five- dimensional parameter space (e.g., representing five color intensities), or intensity vectors in any other appropriate parameter space. [0034] Figure 6B is a diagram illustrating an embodiment of an angular measurement as a function of interferometer gap. In the example shown, θι and θ2 are plotted as an interferometer gap is varied. θ2 is seen to dip at a gap size of interest, while θι is seen to rise, indicating the measured intensity vector swings toward calibration vector λ2. This swing indicates that the measured light is passing through the Fabry -Perot interferometer at the transmission peak associated with calibration vector λ2 and not the transmission peak associated with calibration vector λι. In some
embodiments, θι and θ2 comprise θι and θ2 as shown in Figure 6A.
[0035] In some embodiments, in three dimensions the angles are measured in three dimensions towards calibration vectors. In some embodiments, in the event that the angles shift as a function of gap towards one of the calibration vectors as the gap is incremented or decremented, it is determined that the peak corresponds to the peak associated with the calibration vector.
[0036] Figure 7 is a flow diagram illustrating an embodiment of a process for spectral reading using targeted spatial spectral detection. In 700, a sample set of intensity data is received for an array of spatial locations and a range of spectral configurations. In 701, data in the array of spatial locations is transformed for rotation, translation, and/or key stoning. In 702, a region of interest is determined using the sample set of data. In 704, a spectral peak is determined for the region of interest.
[0037] In some embodiments, the process of Figure 7 is executed by spectral response detector control system 114. In some embodiments, a set of regions of interest is determined in 702. In some embodiments, a spectral peak is determined in 704 for all regions of interest in the set of regions of interest. In some embodiments, more than one spectral peak is determined for a region of interest. In some embodiments, determining one or more spectral peaks comprises decoding a signature. In some embodiments, the signature comprises a tag signature.
[0038] In some embodiments, spectral peak is identified within the spatial array of data associated with one spectral configuration. To correspond this peak with a peak identified within the spatial array of data associated with another spectral configuration, in some embodiments, the data of the spatial array is corrected for rotation, translation, and/or keystoning (e.g., trapezoidal scaling). In some embodiments, spatial array data is processed in an area larger than a region of interest to make sure that a peak that is associated with one spectral configuration can be matched to the peak associated with another spectral configuration. For example, a first spatial array of data is taken associated with one spectral configuration and then at a later time a spatial array of data is taken associated with another spectral configuration. Between the two times, it is possible that the data taking device is rotated, translated, and/or tilted making the data in the two corresponding spatial arrays of data related to each other with a rotation, translation, and/or key stoning translation.
[0039] Figure 8 is a flow diagram illustrating an embodiment of a process for determining a region of interest using a sample set of data. The flow diagram of Figure 8 illustrates an
embodiment of a process for determining a region of interest using a standard deviation calculation. In 800, a data frame associated with a spectral configuration is received. In 801, data frame is transformed for rotation, translation, and/or key stoning. For example, to get correspondence between spatial arrays of data associated with different spectral configurations, transformation is performed if necessary between spatial arrays to make sure that the data in the arrays are spatially associated with each other. In some embodiments, fiducials in the data frames are used to determine the transformations. In 802, the data frame is added to a value accumulator. In the event the data frame comprises the first data frame, the value accumulator is set to the value of the data frame. In 804, the data frame values are squared (e.g., each value of the data frame values is squared). In 806, the squared data frame values are added to a squared value accumulator. In the event the data frame comprises the first data frame, the squared value accumulator is set to the value of the squared data frame. In 808, it is determined whether there are more data frames. In the event there are more data frames, control passes to 800. In the event there are not more data frames, control passes to 810. In 810, the standard deviation is determined. Determining the standard deviation for a set of data frames comprises determining the standard deviation for each point of the data frame, across the set of data frames. The standard deviation for the set of data frames can be determined by first determining the data mean by dividing the values of the value accumulator by the total number of data frames received. A mean of squared data can additionally be determined by dividing the values of the squared value accumulator by the total number of data frames received. The standard deviation for the set of data frames can then be determined by determining the square root of the difference between the mean of squared data and the square of the data mean. In 812, one or more regions of interest are identified by determining local maxima in the standard deviation data (e.g., local maxima across the data frame).
[0040] In some embodiments, the process of Figure 8 implements 702 of Figure 7. In some embodiments, a data frame comprises a two-dimensional array of data pixels. In some
embodiments, each data pixel comprises three colored pixels, each colored pixel associated with a color filter of a different color. In some embodiments, one or more signal processing techniques other than standard deviation computation are used for processing the spectral axis of the data cube (e.g., thresholding, filtering, matched filtering, etc.).
[0041] Figure 9 is a flow diagram illustrating an embodiment of a process for determining a spectral peak for a region of interest. In some embodiments, the process of Figure 9 implements 704 of Figure 7. In 900, an indication to capture further intensity data for one or more regions of interest is provided. In 902, the further intensity data for the one or more regions of interest is received. In 904, the intensity data is transformed into a parameter space (e.g., the parameter space shown in Figure 6A). In 906, the difference between the intensity data in the parameter space and calibration data in the parameter space is determined. In 908, the spectral peak is determined (e.g., by determining a calibration data vector in the parameter space that the intensity data in the parameter space exhibits a movement towards).
[0042] Although the foregoing embodiments have been described in some detail for purposes of clarity of understanding, the invention is not limited to the details provided. There are many alternative ways of implementing the invention. The disclosed embodiments are illustrative and not restrictive.

Claims

1. A system for determining a spectrum, comprising:
an interface configured to receive a sample set of intensity data for an array of spatial locations and a set of spectral configurations; and
a processor configured to:
determine a region of interest using the sample set of intensity data; and determine a spectral peak for the region of interest.
2. The system of claim 1, wherein the sample set of intensity data for the array of spatial locations and the set of spectral configurations comprises an imaging measurement at each location in the array of spatial locations for each spectral configuration in the set of spectral configurations.
3. The system of claim 1, wherein the sample set of intensity data comprises three color intensity data.
4. The system of claim 1, wherein the array of spatial locations comprises a two dimensional array of spatial locations.
5. The system of claim 1, wherein the set of spectral configurations comprises a set of Fabry -Perot interferometer gap sizes.
6. The system of claim 1, wherein the set of spectral configurations comprises a sparse set of gap sizes across a broad gap range.
7. The system of claim 1, wherein the set of spectral configurations comprises a cluster of gap sizes within a predetermined narrow gap range.
8. The system of claim 7, wherein the predetermined narrow gap range is based at least in part on a known reference peak.
9. The system of claim 1, wherein the region of interest is determined using a filter.
10. The system of claim 1, wherein the region of interest is determined using a threshold.
11. The system of claim 1, wherein the region of interest is determined using a calculation of a standard deviation.
12. The system of claim 1, wherein the region of interest comprises one of a plurality of regions of interest.
13. The system of claim 1, wherein the interface is further configured to receive a calibration data.
14. The system of claim 13, wherein the calibration data comprises a calibration set of intensity data for a calibration array of spatial locations and a calibration set of spectral
configurations for a monochromatic source.
15. The system of claim 1, wherein determining the spectral peak for the region of interest comprises providing an indication to capture further intensity data for the region of interest.
16. The system of claim 1, wherein determining the spectral peak for the region of interest comprises receiving further intensity data for the region of interest.
17. The system of claim 1, wherein determining the spectral peak for the region of interest comprises:
transforming the intensity data into a parameter space;
determining a difference between the intensity data in the parameter space and a calibration data in the parameter space; and
determining the spectral peak.
18. The system of claim 1, wherein the region of interest comprises one of a set of regions of interest.
19. The system of claim 18, wherein the processor is further configured to determine a set of spectral peaks for corresponding to the set of regions of interest.
20. The system of claim 1, wherein the processor is further configured to determine more than one spectral peak for the region of interest.
21. The system of claim 1, wherein the processor is further configured to decode a signature.
22. The system of claim 21, wherein the signature comprises a tag signature.
23. A method for determining a spectrum, comprising:
receiving a sample set of intensity data for an array of spatial locations and a set of spectral configurations;
determining, using a processor, a region of interest using the sample set of data; and determining a spectral peak for the region of interest.
24. A computer program product for determining a spectrum, the computer program product being embodied in a non-transitory computer readable storage medium and comprising computer instructions for:
receiving a sample set of intensity data for an array of spatial locations and a set of spectral configurations;
determining a region of interest using the sample set of data; and
determining a spectral peak for the region of interest.
PCT/US2016/061680 2015-11-24 2016-11-11 Tag reading using targeted spatial spectral detection WO2017091370A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201680062190.7A CN108139201A (en) 2015-11-24 2016-11-11 Tag reading using targeted spatial spectral detection
EP16869064.2A EP3380806A4 (en) 2015-11-24 2016-11-11 Tag reading using targeted spatial spectral detection
JP2018521231A JP2019502097A (en) 2015-11-24 2016-11-11 Tag reading using target spatial spectral detection

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US201562259244P 2015-11-24 2015-11-24
US62/259,244 2015-11-24
US15/348,873 2016-11-10
US15/348,873 US10024717B2 (en) 2015-11-24 2016-11-10 Tag reading using targeted spatial spectral detection

Publications (1)

Publication Number Publication Date
WO2017091370A1 true WO2017091370A1 (en) 2017-06-01

Family

ID=58721655

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2016/061680 WO2017091370A1 (en) 2015-11-24 2016-11-11 Tag reading using targeted spatial spectral detection

Country Status (5)

Country Link
US (2) US10024717B2 (en)
EP (1) EP3380806A4 (en)
JP (1) JP2019502097A (en)
CN (1) CN108139201A (en)
WO (1) WO2017091370A1 (en)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10853758B2 (en) * 2017-04-25 2020-12-01 Tag-it Tech, Inc. Micro-taggant based agricultural product tracking system for licensed agricultural products and industries
WO2019127043A1 (en) * 2017-12-26 2019-07-04 深圳达闼科技控股有限公司 Terminal device control method and terminal device
DE102018103869B3 (en) * 2018-02-21 2019-05-09 Physik Instrumente (Pi) Gmbh & Co. Kg Measuring element for an optical measuring device
KR102663185B1 (en) * 2018-08-07 2024-05-03 삼성전자주식회사 Optical emission spectroscopy system and method of calibrating the same, and method of fabricating semiconductor device
KR20210158856A (en) 2019-05-23 2021-12-31 도쿄엘렉트론가부시키가이샤 Optical Diagnostics of Semiconductor Processes Using Hyperspectral Imaging
JP2022053088A (en) * 2020-09-24 2022-04-05 セイコーエプソン株式会社 Optical filter, spectroscopic module, and spectrometric method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080317311A1 (en) * 2004-09-11 2008-12-25 Koninklijke Philips Electronics, N.V. Coherent Scatter Imaging
WO2010019515A2 (en) * 2008-08-10 2010-02-18 Board Of Regents, The University Of Texas System Digital light processing hyperspectral imaging apparatus
US20120206813A1 (en) * 2009-10-25 2012-08-16 Elbit System Electro-Optics Elop Ltd. Tunable spectral filter comprising fabry-perot interferometer
US20140346235A1 (en) * 2009-12-19 2014-11-27 Trutag Technologies, Inc. Labeling and authenticating using a microtag

Family Cites Families (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005083396A1 (en) * 2004-03-01 2005-09-09 Art Advanced Research Technologies Inc. Continuous wave optical imaging assuming a scatter-law
WO2006086889A1 (en) * 2005-02-18 2006-08-24 American Dye Source Inc. Method for encoding materials with a luminescent tag and apparatus for reading same
US20120134582A1 (en) * 2007-01-16 2012-05-31 Chemimage Corporation System and Method for Multimodal Detection of Unknown Substances Including Explosives
US7649634B2 (en) * 2007-10-30 2010-01-19 Mountain View Optical Consultant Corp. Methods and systems for white light interferometry and characterization of films
SE532553C2 (en) * 2008-01-24 2010-02-23 Mikael Lindstrand Method and apparatus for obtaining high dynamic, spectral, spatial and angular resolved radius information
DE112009004707T5 (en) * 2009-04-22 2012-09-13 Hewlett-Packard Development Co., L.P. Spatially varying spectral response calibration data
JP6136357B2 (en) * 2013-02-25 2017-05-31 セイコーエプソン株式会社 Spectrometer, communication system and color management system
JP6003762B2 (en) * 2013-03-27 2016-10-05 セイコーエプソン株式会社 Spectral measurement method, spectroscopic instrument, and conversion matrix generation method
US9495753B2 (en) * 2013-05-30 2016-11-15 Canon Kabushiki Kaisha Spectral image data processing apparatus and two-dimensional spectral apparatus
JP6390090B2 (en) * 2013-11-19 2018-09-19 セイコーエプソン株式会社 Optical filter device, optical module, and electronic apparatus
CN103714546B (en) * 2013-12-27 2016-08-03 北京航空航天大学 A kind of data processing method of imaging spectrometer

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080317311A1 (en) * 2004-09-11 2008-12-25 Koninklijke Philips Electronics, N.V. Coherent Scatter Imaging
WO2010019515A2 (en) * 2008-08-10 2010-02-18 Board Of Regents, The University Of Texas System Digital light processing hyperspectral imaging apparatus
US20120206813A1 (en) * 2009-10-25 2012-08-16 Elbit System Electro-Optics Elop Ltd. Tunable spectral filter comprising fabry-perot interferometer
US20140346235A1 (en) * 2009-12-19 2014-11-27 Trutag Technologies, Inc. Labeling and authenticating using a microtag

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP3380806A4 *

Also Published As

Publication number Publication date
US20180292261A1 (en) 2018-10-11
EP3380806A4 (en) 2019-07-10
JP2019502097A (en) 2019-01-24
US20170146402A1 (en) 2017-05-25
EP3380806A1 (en) 2018-10-03
CN108139201A (en) 2018-06-08
US10024717B2 (en) 2018-07-17

Similar Documents

Publication Publication Date Title
US10024717B2 (en) Tag reading using targeted spatial spectral detection
JP7245835B2 (en) Light field image processing method for depth acquisition
US7336353B2 (en) Coding and modulation for hyperspectral imaging
KR102139858B1 (en) Hyperspectral Imaging Reconstruction Method Using Prism and System Therefor
US20160138975A1 (en) Imaging apparatus comprising coding element and spectroscopic system comprising the imaging apparatus
US8023724B2 (en) Apparatus and method of information extraction from electromagnetic energy based upon multi-characteristic spatial geometry processing
JP6945195B2 (en) Optical filters, photodetectors, and photodetectors
CN109490223B (en) Target detection and identification system and method based on programmable hyperspectral imaging
AU2015309700A1 (en) Imaging method and apparatus
JP7257644B2 (en) Photodetector, photodetector system, and filter array
WO2021085014A1 (en) Filter array and optical detection system
US20170180614A1 (en) Iris imaging
JP2020529602A (en) Coded aperture spectrum image analyzer
EP3707551A1 (en) Imaging method and apparatus using circularly polarized light
Bartlett et al. Anomaly detection of man-made objects using spectropolarimetric imagery
CN116739958B (en) Dual-spectrum polarization super-resolution fusion detection method and system
CN115993329A (en) Handheld multispectral imager
Wong et al. A novel snapshot polarimetric imager
CN114076637A (en) Hyperspectral acquisition method and system, electronic equipment and coding wide-spectrum imaging device
US20240012262A1 (en) Imaging systems and methods for dual depth and polarization sensing
JP7122636B2 (en) Filter array and photodetection system
KR20230113110A (en) Apparatus and Method for Acquiring Hyperspectral Images based on Snapshots using Lens-Less Camera
CN112163627A (en) Method, device and system for generating fusion image of target object
Roper Hyperspectral Image Acquisition and Calibration with Application to Skin Detection Systems

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16869064

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2018521231

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE