US20190073561A2 - Polarimeter with multiple independent tunable channels and method for material and object classification and recognition - Google Patents
Polarimeter with multiple independent tunable channels and method for material and object classification and recognition Download PDFInfo
- Publication number
- US20190073561A2 US20190073561A2 US15/542,064 US201615542064A US2019073561A2 US 20190073561 A2 US20190073561 A2 US 20190073561A2 US 201615542064 A US201615542064 A US 201615542064A US 2019073561 A2 US2019073561 A2 US 2019073561A2
- Authority
- US
- United States
- Prior art keywords
- polarization
- independent
- modulators
- polarimeter
- modulator
- 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.)
- Granted
Links
- 239000000463 material Substances 0.000 title claims abstract description 42
- 238000000034 method Methods 0.000 title claims abstract description 35
- 229940125730 polarisation modulator Drugs 0.000 claims abstract description 92
- 230000010287 polarization Effects 0.000 claims abstract description 51
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 44
- 238000007635 classification algorithm Methods 0.000 claims abstract description 16
- 239000002131 composite material Substances 0.000 claims abstract description 14
- 230000005670 electromagnetic radiation Effects 0.000 claims description 75
- 239000013598 vector Substances 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 12
- 230000001360 synchronised effect Effects 0.000 claims description 11
- 230000003595 spectral effect Effects 0.000 claims description 4
- 230000002708 enhancing effect Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 abstract description 23
- 238000005259 measurement Methods 0.000 abstract description 18
- 230000000712 assembly Effects 0.000 abstract 1
- 238000000429 assembly Methods 0.000 abstract 1
- 230000001419 dependent effect Effects 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 11
- 238000013461 design Methods 0.000 description 9
- 238000003384 imaging method Methods 0.000 description 8
- 239000004033 plastic Substances 0.000 description 8
- 229920003023 plastic Polymers 0.000 description 8
- 238000012545 processing Methods 0.000 description 7
- 239000000523 sample Substances 0.000 description 7
- 238000012549 training Methods 0.000 description 6
- 230000001965 increasing effect Effects 0.000 description 5
- 230000003466 anti-cipated effect Effects 0.000 description 4
- 238000012512 characterization method Methods 0.000 description 4
- 238000010801 machine learning Methods 0.000 description 4
- 239000002184 metal Substances 0.000 description 4
- 229910052751 metal Inorganic materials 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 239000013078 crystal Substances 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000005286 illumination Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- 150000002739 metals Chemical class 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013146 percutaneous coronary intervention Methods 0.000 description 2
- 230000005855 radiation Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012706 support-vector machine Methods 0.000 description 2
- 241000195493 Cryptophyta Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 239000005427 atmospheric aerosol Substances 0.000 description 1
- -1 but not limited to Substances 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000000356 contaminant Substances 0.000 description 1
- 239000013068 control sample Substances 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000000572 ellipsometry Methods 0.000 description 1
- 239000002657 fibrous material Substances 0.000 description 1
- 231100001261 hazardous Toxicity 0.000 description 1
- 239000000976 ink Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000000386 microscopy Methods 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
- 239000013077 target material Substances 0.000 description 1
- 239000004753 textile Substances 0.000 description 1
- 239000010409 thin film Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/60—Extraction of image or video features relating to illumination properties, e.g. using a reflectance or lighting model
-
- G06K9/6232—
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J4/00—Measuring polarisation of light
- G01J4/04—Polarimeters using electric detection means
-
- 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/21—Polarisation-affecting properties
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B26/00—Optical devices or arrangements for the control of light using movable or deformable optical elements
- G02B26/08—Optical devices or arrangements for the control of light using movable or deformable optical elements for controlling the direction of light
- G02B26/10—Scanning systems
- G02B26/105—Scanning systems with one or more pivoting mirrors or galvano-mirrors
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
- G02B27/28—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 for polarising
- G02B27/283—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00 for polarising used for beam splitting or combining
-
- G—PHYSICS
- G02—OPTICS
- G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
- G02B5/00—Optical elements other than lenses
- G02B5/30—Polarising elements
- G02B5/3083—Birefringent or phase retarding elements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/10—Image acquisition
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/70—Arrangements for image or video recognition or understanding using pattern recognition or machine learning
- G06V10/77—Processing image or video features in feature spaces; using data integration or data reduction, e.g. principal component analysis [PCA] or independent component analysis [ICA] or self-organising maps [SOM]; Blind source separation
- G06V10/7715—Feature extraction, e.g. by transforming the feature space, e.g. multi-dimensional scaling [MDS]; Mappings, e.g. subspace methods
Definitions
- Embodiments of the present invention are related to the field of remote sensing or testing by polarized light examination, which comprises art from the field of polarimeters and art from the field of automated or aided discrimination, classification, and recognition of objects, materials, and material states through application of learning algorithms to polarized-light measurements.
- the object, material, or material state of interest is often called the target and this field of remote-sensing is often termed automated target recognition (ATR.)
- polarimeters can be categorized according to their radiation source as either a) active, which employ a laser or other controlled electromagnetic-radiation (EMR) source or b) passive, which employ solar or thermal electromagnetic radiation, or c) hybrid, which employ both active and passive sources. Polarimeters can be further subdivided according to the measurements they make, their polarization modulators, and their applications. Stokes polarimeters measure the Stokes parameters (or Stokes vector), which completely characterize the polarization properties of the EMR. If the source is unpolarized, as are most passive sources, then the Stokes parameters reveal all obtainable polarization properties of the illuminated object, hence Stokes polarimeters can be applied for passive remote sensing.
- EMR electromagnetic-radiation
- Mueller-matrix polarimeters employ a controlled polarized source, often a laser, to measure the Mueller matrix, which characterizes the complete linear response of a material or object to an electromagnetic wave.
- U.S. Pat. No. 7,218,398 teaches application of the data obtained with a Mueller-matrix polarimeter to an inverse problem for material characterization.
- Ellipsometers are a subclass of Mueller-matrix polarimeters that solve the inverse problem under the assumption of no scattering or birefringence, which results in a sparse Mueller matrix, although generalized ellipsometry, as taught for instance in U.S. Pat. No. 7,218,398, relaxes certain of these assumptions.
- Most polarimeters designed for characterization are laboratory instruments that also comprise a sample stage to accurately control the position and orientation of the material or object sample. Polarimeters can also be categorized according to their polarization modulators.
- PEM photoelastic modulators
- LCVR liquid-crystal variable retarders
- microgrid arrays have also been employed.
- Polarimeters for remote sensing must be compact, lightweight, and typically fast, and their sensing schemes must work for materials and objects in generally arbitrary locations and orientations within the field-of-view. The most common approach is to limit the number of measurements or channels to 2, achievable in real-time using a polarizing prism or beamsplitter in the receiver, and to form the difference of the 2 channels.
- This class of polarimeter for remote sensing is taught for instance in U.S. Pat. No. 8,116,000 for passive sources and in U.S. Pat. No. 7,333,897 for active sources.
- most polarimeters of this type measure 2 orthogonal channels.
- the channels are often aligned parallel to and orthogonal to (or crossed with) the polarization direction of the illuminating EMR.
- Certain active polarimeters for instance those taught in U.S. Pat Nos. 5,956,147 and 6,060,710, employ photoelastic modulators (PEM), which require complex calibration and demodulation routines but can measure the Mueller matrix sufficiently fast for remote sensing.
- PEM photoelastic modulators
- detect, discriminate, classify, recognize, and identify have formal meanings, although the terms are casually interchanged in many publications.
- the term “detect” means to record EMR reflected or transmitted by the target;
- the term “discriminate” means to distinguish the target from the background;
- the term “classify” means to distinguish the target or targets among various classes including the background and other non-targets, which are often referred to as clutter;
- the terms “recognize” and “identify” imply even more specificity. Distinguishing materials (eg, metal, plastic, or wood) is generally considered classification, while distinguishing an object as a known combination of materials and possibly shape may be considered recognition.
- Electromagnetic sensors that can rapidly and non-destructively discriminate, classify, and recognize materials, material states, and objects composed of those materials and states have many useful applications depending on the material sensitivities of the measured characteristics of the electromagnetic wave. Polarization at optical and near-optical frequencies is sensitive to broad ranges of materials and material states including metals, plastics, fibrous materials including fiber composites, crystals, chiral molecules, and biological tissues, and even broader ranges if patterned, stressed, or damaged materials are considered. Polarimeters have been demonstrated for the discrimination and classification of such materials and states, usually by increasing the contrast between a target material or material state and other materials or material states in the observation of reflected or transmitted EMR.
- polarimeters increase contrast by differencing irradiance channels defined by different settings of their polarization modulators. Ideally the irradiances of the target are different in the two channels, while the irradiances of non-targets are similar in the two channels, such that the difference channel provides increased target contrast.
- the polarimeter irradiance channels can be in the form of digital images or in the form of detector signals, multiplexed in time and/or space. For high-speed performance, polarimeters are usually limited to two simultaneous channels, as in the case of polarization-difference imaging (PDI), as taught for instance in U.S. Pat. Nos. 4,881,818, 5,929,443, and 8,116,000.
- PDI polarization-difference imaging
- Active polarimeters comprise a source of controlled, polarized EMR, often a laser, that is typically formed into a beam that illuminates the sample or scene of interest. Compared to passive polarimeters, which utilize solar illumination or thermal emission, active polarimeters provide access to additional signature components and perform independently of weather and time-of-day. Active polarimeters designed for characterization, as taught for optical and infrared frequencies in U.S. Pat. Nos. 4,306,809, 5,247,176, 5,956,147, and 7,218,398, typically measure the complete polarization signature, also known as the Mueller matrix, for illumination at given angles and electromagnetic frequencies.
- the complete polarization signature also known as the Mueller matrix
- S 0 is the irradiance
- S 1 is the degree of linear polarization in the horizontal-vertical coordinate system
- S 2 is the degree of linear polarization at ⁇ 45° to the horizontal-vertical coordinates
- S 3 is the degree of circular polarization
- [. . . ] T denotes a matrix transpose.
- the Stokes parameters define the polarization state of the EMR.
- the Mueller matrix is the linear transformation of the Stokes parameters upon interaction of the electromagnetic wave with a material object, written mathematically as
- the Stokes parameters are often referred to as the Stokes vector, which is correct only in the computing sense of a one-dimensional array of numbers. Since irradiance cannot be negative, both Stokes parameters and Mueller matrices lack additive inverses, so are not vectors in the algebraic sense.
- Active polarimeters designed for remote sensing typically do not measure, the complete Mueller matrix, but only parts of the signature that show contrast between the target of interest and the anticipated background and clutter, allowing SWaP to be reduced and speed increased due to fewer components and fewer measurements.
- most atmospheric aerosols and clouds exhibit a diagonal Mueller matrix with a single independent parameter that can be obtained from measurement of M 11 , which can be accomplished with an active polarimeter by forming the difference of two orthogonal channels.
- active polarization signatures of many relevant materials are actually multi-dimensional, and discrimination, classification, and recognition performance can be improved by applying learning algorithms to measured multi-dimensional signatures.
- multi-dimensional measurements are usually needed to achieve classification performance that is relatively invariant to the target orientation or pose, which is usually a critical requirement for remote sensing.
- Fielding an active polarimeter for high-speed, high-performance classification or ATR therefore requires a new class of polarimeter that can measure specific multi-dimensional signatures while maintaining the low SWaP and high speed achievable by 2-channel polarimeters.
- One embodiment of the present invention provides for a registered-channel multiplexer (RCM) having a scanning element positioned to direct a beam from a controlled source of electromagnetic radiation (EMR) serially among a set of independent modulators.
- RCM registered-channel multiplexer
- EMR electromagnetic radiation
- FIG. 3 One version of the RCM is illustrated in FIG. 3 .
- Each member of the set of independent modulators defines a channel that is configured to modulate a state of the beam among a set of independent states to produce a modulated beam.
- An assembly of reflectors is configured to redirect the modulated beam to a co-registered path that is pointed at a field-of-view.
- the scanning element can be a scanner mirror attached to a first mechanical scanner, for instance a galvanometer scanner, that comprises a rotation axis, and the assembly of reflectors can be an assembly of mirrors.
- the scanner mirror can be a composite mirror comprising a flat mirror parallel to the rotation axis of the first mechanical scanner and a second wedge mirror mounted at an angle to the flat mirror, for example, at an angle of about 45 degrees.
- the assembly of mirrors may comprise one or more flat scanner mirrors attached to a second set of mechanical scanners.
- the set of independent modulators can be selected from the group consisting of lenses, polarization modulators, spectral filters, diffractive array generators, holograms, amplitude-phase masks and any combination thereof.
- One or more members of the set of independent modulators can be tunable, and one or more different members can be a beam block as illustrated in FIG. 4 .
- the present invention provides for a sensor for classifying or recognizing a target or targets within a field-of-view (FOV) using multiple independent channels.
- the sensor comprises a source of controlled electromagnetic radiation (EMR) positioned in a path with a detector and the FOV positioned there between.
- the sensor includes a registered-channel multiplexer (RCM) comprising a scanning element positioned to direct a beam from the source of controlled EMR serially among a first set of independent modulators, wherein each member of the first set of independent modulators corresponds to a channel and the scanner and modulators are configured to modulate a state of the beam among a set of independent states to produce a modulated beam.
- EMR electromagnetic radiation
- RCM registered-channel multiplexer
- An assembly of reflectors is configured to redirect the modulated beam to a co-registered path that is pointed at the FOV.
- An EMR collector is positioned to direct EMR reflected from or transmitted by objects within the FOV to a second modulator independent of the first set of independent modulators.
- the detector is positioned to receive the EMR from the second modulator, wherein the detector produces a set of signals that are synchronized with the set of channels formed by the first set of modulators and the second modulator.
- the first set of independent modulators and the second independent modulator can be selected from the group consisting of lenses, polarization modulators, spectral filters, diffractive array generators, holograms, amplitude-phase masks, and any combination thereof.
- the sensor includes a processor that is connected with a memory, wherein the processor is configured to execute a classification algorithm stored in the memory by applying a user-adjustable threshold to the set of detector signals to assign a class label (which may also be stored in memory) to the set of detector signals and wherein the classification algorithm specifies the function of the threshold in terms of the set of channels formed by the first set of modulators and the second modulator.
- the class label can be either binary, assuming the values 0 or 1 only, or multi-valued.
- the first set of independent modulators and the second independent modulator can be tunable, and one or more different members of the first set of independent modulators can be a beam block.
- the EMR directed by the collector to the second modulator may include EMR from a passive source illuminating or emanating from objects within the FOV, in which case the EMR directed by the collector to the second modulator while the scanning element is directed at the beam block can produce a detector signal corresponding to a passive channel as a member of the set of detector signals.
- the set of detector signals is a set of registered digital images and a class label is assigned to each pixel to form a digital classification image.
- the class labels can be assigned at a rate of 20 Hz or greater and for imaging the class labels can be assigned and the classification image electronically displayed at a rate of 20 frames-per-second or greater.
- the FOV can be scanned or swept over an area or a volume by a gimbal or a moving platform on which the sensor is mounted. Further, the set of detector signals can be synchronized with the set of channels formed by the first set of modulators and the second modulator by a lock-in amplifier using a timing signal provided by the scanning element.
- a polarimeter for classifying or recognizing a target or targets within a field-of-view (FOV) using multiple independent polarization channels.
- the polarimeter comprises a source of controlled electromagnetic radiation (EMR) positioned in a path with a detector and the FOV positioned there between.
- the polarimeter includes a registered-channel multiplexer (RCM) comprising a scanning element positioned to direct a beam from the source of controlled EMR serially among a first set of independent polarization modulators.
- EMR electromagnetic radiation
- RCM registered-channel multiplexer
- the first set of independent polarization modulators wherein each member of the set corresponds to a polarization channel, is configured to modulate the polarization state of the beam among a set of independent polarization states to produce a modulated beam.
- An assembly of reflectors is configured to redirect the modulated beam to a co-registered path that is pointed at the FOV.
- An EMR collector is positioned to direct EMR reflected from or transmitted by objects within the FOV to a second polarization modulator independent of the first set of polarization modulators.
- the detector is positioned to receive the EMR from the second polarization modulator, wherein the detector produces a set of signals that are synchronized with the set of channels formed by the first set of polarization modulators and the second polarization modulator.
- the polarimeter includes a processor connected with a memory, wherein the processor is configured to execute a classification algorithm stored in the memory by applying a user-adjustable threshold to the set of detector signals to assign a class label (which may also be stored in memory) to the set of detector signals, wherein the classification algorithm specifies the function of the threshold in terms of the set of channels formed by the first set of polarization modulators and the second polarization modulator.
- the class label can be either binary, assuming the values 0 or 1 only, or multi-valued.
- the first set of independent polarization modulators and the second independent polarization modulator can be tunable, and one or more member of the first set of polarization modulators can be a beam block.
- the EMR directed by the collector to the second modulator may include EMR from a passive source illuminating or emanating from the field-of-view, in which case the EMR directed by the collector to the second modulator while the scanning element is directed at the beam block can produce a detector signal corresponding to a passive channel as a member of the set of detector signals.
- the set of detector signals is a set of registered digital images and a class label is assigned to each pixel to form a digital polarization classification image (PCI).
- the class labels can be assigned at a rate of 20 Hz or greater and PCIs can be electronically displayed at a rate of 20 frames-per-second or greater.
- the FOV can be scanned or swept over an area or a volume by a gimbal or a moving platform on which the polarimeter is mounted.
- the set of detector signals can be synchronized with the set of channels formed by the first set of polarization modulators and the second polarization modulator by a lock-in amplifier using a timing signal provided by the scanning element.
- the channels are defined by a subspace-projection algorithm that applies a gradient operator on a vector space of reduced Mueller matrices and maximizes a resulting contrast parameter
- the user-adjustable threshold is a function of the channels derived by the subspace-projection algorithm.
- the contrast parameter and the function of the user-adjustable threshold are defined by an external machine-learning algorithm.
- Yet another embodiment of the present invention provides a method for classifying or recognizing a target or targets within a field-of-view (FOV) by using a polarimeter with multiple independent channels.
- the method includes positioning a source of controlled electromagnetic radiation (EMR) in a path with a detector and the FOV positioned there between.
- EMR controlled electromagnetic radiation
- the method employs a registered-channel multiplexer (RCM), which uses a scanning element to serially direct a beam from the controlled source of EMR among a first set of independent polarization modulators, producing a modulated beam, wherein each member of the set of independent polarization modulators corresponds to a channel and the set is configured to modulate the polarization state of the beam among a set of independent polarization states, and an assembly of reflectors to redirect the modulated beam to a co-registered path pointed at the FOV.
- RCM registered-channel multiplexer
- the method further includes collecting a portion of the EMR reflected from or transmitted by objects within the FOV with an EMR collector that directs the EMR to a second polarization modulator independent of the first set of polarization modulators, and receiving the EMR from the second polarization modulator at the detector, wherein the detector produces a set of detector signals that are synchronized with the set of channels formed by the first set of polarization modulators and the second polarization modulator.
- the method applies a classification algorithm and a user-adjustable threshold to the set of detector signals to assign a class label to the set of detector signals, wherein the classification algorithm specifies the function of the threshold in terms of the set of channels formed by the first set of polarization modulators and the second polarization modulator and the class label can be either binary, assuming the values 0 or 1 only, or multi-valued.
- the first set of independent polarization modulators and the second independent polarization modulator can be tunable, and one or more members of the first set of polarization modulators can be a beam block.
- the EMR directed by the collector to the second modulator may include EMR from a passive source illuminating or emanating from objects within the FOV, in which case the EMR directed by the collector to the second modulator while the scanning element is directed at the beam block can produce a detector signal corresponding to a passive channel as a member of the set of detector signals.
- the set of detector signals is a set of registered digital images and a class label is assigned to each pixel to form a digital polarization classification image (PCI).
- the class labels can be assigned at a rate of 20 Hz or greater and PCIs can be electronically displayed at a rate of 20 frames-per-second or greater.
- the method may also include scanning or sweeping the FOV through an area or a volume using a gimbal or a moving platform on which the polarimeter is mounted.
- the method may also include synchronizing the set of detector signals with the set of channels formed by the first set of polarization modulators and the second polarization modulator using a lock-in amplifier that reads a timing signal from the scanning element.
- One embodiment includes defining the polarization channels using a subspace-projection algorithm that applies a gradient operator on a vector space of reduced Mueller matrices and maximizes a resulting contrast parameter, and the user-adjustable threshold can be a function of the channels derived by the subspace-projection algorithm.
- Another embodiment includes applying an external machine-learning algorithm to define the contrast parameter and the function of the user-adjustable threshold.
- Another embodiment comprises a hybrid polarimeter and method of use that can achieve data rates up to video-rate in imaging mode with sufficiently small SWaP for field implementations.
- the new subspace-projection algorithm specifies the set of polarization-modulator settings or channels that maximize contrast between the target or state of interest and the anticipated backgrounds and other objects or states that might be observed.
- Each polarimeter channel corresponds to an irradiance measurement or a detector signal, and measurements in multiple channels are generally combined to form a set of detector signals or a projection that is used to test whether the object or state corresponds to the target or not.
- the projection can be either a one-dimensional or multi-dimensional irradiance signal depending on the type of classification algorithm used.
- a simple irradiance threshold serves as the classifier, such that an irradiance above or below the threshold, or between two thresholds, corresponds to a positive target classification.
- every pixel is an independent sensor, and positive classifications can be represented by white pixels with all other pixels black, providing binary polarization-classification images (PCI) that are advantageous for fast interpretation and practical data transmission and storage.
- PCI binary polarization-classification images
- the subspace-projection algorithm derives the polarimeter channels from the available polarization modulators and the Mueller matrices of the anticipated objects, materials, and/or material states, with the Mueller matrices known a priori from a model and/or training characterization measurements.
- the subspace-projection algorithm provides the maximally-biased one-dimensional classifier if only a single target. Mueller matrix is used as input, and provides more generalized classifiers if non-target Mueller matrices are incorporated into the input (see Detailed Description.)
- the method can be applied to an existing polarimeter, to derive the optimal set of realizable channels, or at the polarimeter design stage to derive the minimum set of channels that can achieve the required classification or ATR performance.
- Another embodiment of the present invention provides a method of enhancing the contrast of a target or targets recorded by a polarimeter.
- This method includes obtaining modeled or measured Mueller matrices of the target or targets and other materials or objects to be recorded by the polarimeter and applying the subspace-projection algorithm to the Mueller matrices and to the set of polarization modulators of the polarimeter to specify the settings of the polarization modulators, wherein the subspace-projection algorithm applies a gradient operator on a vector space of reduced Mueller matrices and maximizes a resulting contrast parameter.
- Another embodiment includes applying an external machine-learning algorithm to define the contrast parameter.
- the current invention includes new hardware designs for a polarization multiplexer to achieve data rates up to video-rate in imaging mode. While passive polarimeters contain polarization modulators in the receiver only, active polarimeters comprise modulators in both the transmitter and receiver.
- One embodiment of the current invention employs a conventional polarization modulator, a polarizing prism, in the receiver, and a new temporally-multiplexed polarization modulator (or polarization multiplexer) in the transmitter, which is a form of the more general registered-channel multiplexer (RCM.)
- the RCM can assume various geometries, all of which utilize a scanner to switch or scan the beam among multiple independent and tunable states at high speed, possibly including a null state that corresponds to a passive channel.
- Achieving video-rate multiplexing with a mechanical scanner requires a rotor of lower inertia than the hollow-core motors used in rotating-retarder polarimeters like those taught in U.S. Pat. Nos. 4,306,809, 5,247,176, and 7,218,398.
- One embodiment of the invention employs a low-inertia galvanometer scanner that switches a laser beam among a set of tunable polarization modulators.
- Each modulator in the transmitter, in combination with a static polarizing prism in the receiver constitutes 2 orthogonal polarimeter channels, while overall the polarimeter may have 8 or more channels.
- Each polarization modulator in the transmitter can be independently tuned to realize multiple arbitrary channels.
- Classification of various materials on natural backgrounds has been demonstrated in field tests, at approximately 20 m horizontal range, employing the polarimeter main embodiment of the invention.
- Applications of the main embodiment include, but are not limited to, classification of hazardous objects for route clearance, classification of materials in mining, manufacturing, waste, and recycling streams, and real-time non-destructive testing of vehicle and structural materials including, but not limited to, fiber composites and non-cubic metals.
- the invention can be customized for manufacturing applications including, but not limited to, monitoring of thin-film, textile, polymer, chiral pharmaceutical, composite material, and non-cubic metal fabrication and processing, including mechanical and laser machining.
- Smaller embodiments of the invention are applicable to recognition of security inks for asset tagging and product authentication, quantification of micro-plastics, biomedical-tissue assays, and general polarization microscopy, while larger embodiments that utilize high-power lasers are applicable for long-range material and object classification and recognition from resolved or unresolved imagery, for instance surveillance and reconnaissance, damage assessment, search and rescue, mapping and, assessment of contaminants and debris, particularly plastics, and space-object classification.
- FIG. 1 is a taxonomy of polarimeter types based on electromagnetic-radiation (EMR) source, measurements made, polarization modulators employed, and applications.
- EMR electromagnetic-radiation
- FIG. 2 is an illustration of active monostatic, passive bistatic, and hybrid active-passive sensor geometries.
- FIG. 3 is an illustration of the polarization multiplexer of an embodiment of the present invention that provides 8 active tunable polarization channels or 4 independent active tunable polarization-difference channels.
- FIG. 4 is an illustration of the polarization multiplexer of an embodiment of the present invention that provides 6 active tunable polarization channels and 2 passive tunable polarization channels.
- FIG. 5 is an illustration of the composite wedge scanner mirror employed in an embodiment of the present invention.
- FIG. 6A-B is a comparison of a conventional image and the corresponding polarization-classification image (PCI) of plastic bags on a beach obtained in field tests of an embodiment of the invention at 15 m horizontal range.
- PCI polarization-classification image
- FIG. 7 is a flowchart of an embodiment of the invention designed to map the distribution of a selected material at a test site.
- FIG. 8A-D is a set of polarization-classification images (PCI) of plastic bags on a beach obtained using different classifier thresholds.
- PCI polarization-classification images
- the invention can be described in detail, with reference to the figures, by first describing the subspace-projection algorithm used to derive polarimeter channels and corresponding classifiers from the available polarization modulators and a priori Mueller matrices, then describing the specific hardware by which the new polarimeter design can be reduced to practice in a hybrid polarimeter that can achieve data rates up to video-rate in imaging mode with sufficiently small size, weight, and power (SWaP) to allow field implementations.
- SWaP size, weight, and power
- the subspace-projection algorithm derives what is known in the technical literature as a partial Mueller-matrix polarimeter (pMMP.)
- pMMP partial Mueller-matrix polarimeter
- the newly disclosed subspace-projection algorithm has advantages over previously published pMMP design methods in that it avoids use of matrix pseudoinverses, which are prone to errors, and imposes no limitations on the objects, materials, or material states that can be observed by the polarimeter.
- the subspace-projection algorithm is based on application of vector algebra on the vector space of reduced Mueller matrices formed by reducing the measured irradiance as
- I is the irradiance actually measured by the polarimeter and I U is the irradiance that a conventional non-polarizing camera would measure with unpolarized illumination.
- the reduced irradiance I′ can be positive or negative. It can be shown that the reduced irradiance is expressible as a vector product of the reduced polarimeter instrument matrix and the reduced Mueller matrix of the illuminated object
- the parameters ⁇ are retarder-waveplate angular orientations, but ⁇ can represent any combination of polarization-modulator settings.
- the algorithm is also general enough to derive channels and classifiers for passive polarimeters.
- the polarimeter instrument matrix In order to maximize target contrast in a polarimeter measurement, the polarimeter instrument matrix must be chosen to maximize its projection onto the subspace of the reduced target matrix ⁇ tilde over (M) ⁇ ′ T .
- This subspace-projection can be formalized by defining a differentiable target matrix and applying the gradient operator on the space of reduced Mueller matrices, defined as
- the representation of the target in the polarimeter measurement is then maximized by maximizing the contrast parameter, which is a vector dot product, over the set of modulator parameters ⁇ , that is
- FIG. 2 Field sensors that utilize different electromagnetic-radiation (EMR) sources are illustrated in FIG. 2 .
- the sensor field-of-view (FOV) may be scanned or swept over an area or volume, for instance if the sensor is mounted to a vehicle, or objects may move through a fixed FOV, for instance on a conveyor belt.
- Each recorded frame may or may not contain a target, although for the purposes of illustration a target is assumed in FIG. 2 .
- the target 1 is illuminated by the probe beam 2 , typically a laser beam, that is projected by the sensor transmitter 3 .
- the target may also be illuminated by a passive source 4 , which is typically the sun but can also be thermal radiation from the object itself if infrared frequencies are recorded.
- the reflected probe EMR 5 is recorded by the sensor receiver 6
- the receiver 7 records only reflected passive EMR 8
- a hybrid active-passive sensor illuminates the target with the probe beam 9 projected by the transmitter 10 and records both the reflected probe EMR 11 and reflected passive EMR 12 in the receiver 13 .
- the polarization multiplexer is enclosed in the sensor transmitter, embodiments of which are illustrated at 3 and 10 in FIG. 2 .
- the polarization multiplexer is a form of registered-channel multiplexer (RCM). Birefringent crystal waveplates are assumed for the polarization modulators, although other polarization modulators can be used.
- Conventional polarimeters employ waveplates in rotating hollow-core motors, which can be very accurate but not fast enough to achieve polarization video.
- an embodiment of the present invention employs low-inertia scanners, for instance galvanometer scanners, to switch among multiple arms that comprise independent tunable polarization modulators.
- a laser 14 emits a collimated beam of light that strikes a mirror 15 that is scanned by a galvanometer scanner 16 serially among multiple independent arms 17 a - d , each of which, in combination with the static polarizing prism in the receiver (not shown), defines a pair of orthogonal channels.
- the scanner directs the beam through the first arm 17 a , which comprises the first polarization modulator 18 , which is set to a channel specified by the subspace-projection algorithm.
- the modulator 18 may be a rotatable crystal waveplate or a pair of waveplates, but can be other types of polarization modulator as well.
- the beam is reflected by a fold-mirror assembly 19 back to the scanner mirror 15 and then to the field mirror 20 , which reflects the beam toward the field-of-view.
- the beam is generally expanded or focused by projector lenses located at beam position 21 .
- the scanner dwells on the first arm 17 a for a predetermined time interval, then switches to the second arm 17 b for a second time interval, generally of the same duration as the first time interval, then switches through the remaining arms 17 c - d before rewinding to the first arm 17 a and repeating the scan.
- the subsystem illustrated in FIG. 3 is generally called a polarization-state generator (PSG) and may have more or fewer than 4 arms.
- PSG polarization-state generator
- the set of irradiances or images collected for the set of scanner positions forms the projection vector used for classification.
- the dwell time on each arm is determined by the number of arms, the EMR power collected from the FOV, the detector sensitivity, and the sensor speed requirement, and is generally limited by the scanner switching speed.
- the scanner has a characteristic switching time.
- a typical sum of dwell and switching times with 4 arms and a high-performance galvanometer scanner is 12.5 ms, which implies 50 ms for collection of each image or projection vector. With sufficiently fast data-processing this enables PCI video at 20 frames-per-second. Faster scanners would enable higher data rates.
- FIG. 4 illustrates a PSG with a null channel recorded while the scanner dwells on a beam block 22 .
- the null channel records a pair of passive channels that are combined with the active channels obtained from the remaining arms to realize a hybrid active-passive polarimeter.
- the dwell time on the null channel may be the same, longer, or shorter than the dwell times on the active channels if the reflected power from the passive source is less than or greater than the reflected laser power.
- the scanner timing signal can provide a reference frequency for lock-in detection.
- FIG. 5 illustrates the composite scanner mirror of an embodiment of the present invention.
- the main parts of the composite scanner mirror are the mounting boss 23 , the flat mirror 24 , and the wedge mirror 25 .
- the mirrors 24 and 25 must have high optical quality, at least in their central areas, their surfaces must be flush with and centered on the rotation axis, and the entire assembly must have low inertia in order to achieve video-rate data.
- the composite mirror used for field tests was constructed with 3D-printed plastic parts and has a moment of inertia of 6.4 g/cm 2 .
- Alternative embodiments of a galvanometer-based multiplexer can be used, for instance embodiments that do not require the composite wedge mirror.
- the disclosed polarization multiplexer is a form of registered-channel multiplexer (RCM), which generally switches in time among different channels that produce spatially- and angularly-registered beams at the output.
- RCM registered-channel multiplexer
- Different forms of RCM can be realized by substituting other modulators for the polarization modulators of the illustrated embodiment of the present invention.
- the classification performance and data-processing speed depend on the classification algorithm, a great variety of which may be employed.
- the classification algorithm, or simply classifier can be one-dimensional or multi-dimensional, corresponding to the number of measured channels.
- Two-channel polarimeters typically use one-dimensional classifiers based on the difference between the measured channels, hence the term polarization-difference imaging (PDI.)
- PDI polarization-difference imaging
- a one-dimensional classifier makes a positive target classification when the irradiance projection measured in a specified channel is above or below a user-defined threshold, or between two thresholds, and a negative classification otherwise.
- Classifiers can be defined in terms of the channels derived by the subspace-projection algorithm of the present invention.
- Classifiers are described in terms of the bias-variance tradeoff, which balances performance on training datasets with the ability to generalize to unforeseen data.
- the smallest possible training dataset comprises the target only and results in a high-bias classifier that can perform well on scenes containing the target only but may suffer frequent false-alarms when clutter objects are observed.
- a one-dimensional classifier based on a channel derived from Eq. 7 of the subspace-projection algorithm is a high-bias classifier.
- Variance is the ability of the classifier to perform well on scenes that contain target variations and/or background or clutter not represented in the training dataset. Increasing the classifier dimensionally can increase the variance as the target signature can be projected into a higher dimensional space where its projection is more unique.
- the classifier dimension can be increased by combining two or more channels derived from Eq. 7 of the subspace-projection algorithm.
- the classifier bias can be lowered by using a larger training dataset that includes signatures of the background and/or anticipated clutter objects.
- the background signature can be incorporated by maximizing the generalized contrast parameter R,
- the classifier In image format the classifier is applied to each pixel independently. Positive classifications can be represented by white pixels with all other pixels black, resulting in binary polarization classification images (PCI) as illustrated in FIG. 6 .
- Binary images are advantageous for fast interpretation and practical data transmission and storage.
- the target In the conventional unpolarized image of FIG. 6 the target, the plastic bags 26 , are nearly indiscernible from the algae, twigs, and other clutter objects 27 , especially if scanning over a large area at video rate, while the PCI reveals the target clearly.
- the PCI shown in FIG. 6 was obtained in field tests with the main embodiment of the invention in full daylight at 15 m horizontal range.
- FIG. 7 is a flowchart of an embodiment of the invention designed to map the distribution of a selected material at a test site using a classifier derived by the subspace-projection pMMP-design algorithm.
- the components inside the dashed box 28 are on-board the field sensor and the flowchart corresponds to a single image or video frame.
- the sensor field-of-view may be scanned or swept over the test area by a motion system not shown.
- the subspace-projection pMMP-design algorithm defines the classifier, which is stored in computer memory 29 , and the settings of the polarization modulators, which can be changed by the user via the tuner 30 .
- the PSG scanner 31 switches the beam among the independent modulators in the RCM and provides a timing signal to trigger the camera 32 to record corresponding images in the channels 33 defined by the subspace-projection algorithm.
- These images and the corresponding channel specifications constitute the test data 34 that is input to a real-time processor such as a general-purpose graphics-processing unit (GPGPU) 35 .
- GPGPU general-purpose graphics-processing unit
- a parallel-processing architecture such as GPGPU is advantageous since each pixel is processed as an independent sensor, although other real-time architectures may be employed.
- the user can adjust the classifier threshold via the control 36 . From these inputs the processor produces the class labels 37 , which are 0 or 1 at each pixel for binary PCI.
- the frame is stored on the site map with the corresponding global-positioning system (GPS) coordinates provided by an external GPS module 38 .
- GPS global-positioning system
- each frame can be assigned a class label, 0 or 1, to produce a density map without images.
- the GPS module and stored site map may also be on-board the sensor.
- FIG. 8A-D is a set of binary PCIS of the scene illustrated in FIG. 6A obtained by the main embodiment of the invention as the classifier threshold is adjusted.
- Both the classifier training and the polarimeter testing and ordinary use are executed by a general or specific-purpose computer or distributed system programmed with computer software implementing the steps described above, which computer software may be in any appropriate computer language, including C++, Python, FORTRAN, BASIC, Java, assembly language, distributed programming languages, etc.
- the apparatus may also include a plurality of such computers/distributed systems (e.g., connected over the Internet and/or one or more intranets) in a variety of hardware implementations.
- data processing can be performed by an appropriately programmed microprocessor, computing cloud, field-programmable gate array (FPGA), general-purpose graphics-processing unit (GPGPU), or the like, in conjunction with appropriate memory, network, and bus elements.
- All computer software disclosed herein may be embodied on any computer-readable medium (including combinations of mediums), including without limitation CD-ROMs, DVD-ROMs, hard drives (local or network storage device), USB keys, other removable drives, ROM, and firmware.
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Optics & Photonics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Software Systems (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Databases & Information Systems (AREA)
- Computing Systems (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Artificial Intelligence (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Mechanical Optical Scanning Systems (AREA)
- Mechanical Light Control Or Optical Switches (AREA)
Abstract
Description
- This invention was made with Government support under contract W909MY-12-C-0023 awarded by the United States Army. The United States Government has certain rights in the invention.
- Embodiments of the present invention are related to the field of remote sensing or testing by polarized light examination, which comprises art from the field of polarimeters and art from the field of automated or aided discrimination, classification, and recognition of objects, materials, and material states through application of learning algorithms to polarized-light measurements. The object, material, or material state of interest is often called the target and this field of remote-sensing is often termed automated target recognition (ATR.)
- Referring to
FIG. 1 , polarimeters can be categorized according to their radiation source as either a) active, which employ a laser or other controlled electromagnetic-radiation (EMR) source or b) passive, which employ solar or thermal electromagnetic radiation, or c) hybrid, which employ both active and passive sources. Polarimeters can be further subdivided according to the measurements they make, their polarization modulators, and their applications. Stokes polarimeters measure the Stokes parameters (or Stokes vector), which completely characterize the polarization properties of the EMR. If the source is unpolarized, as are most passive sources, then the Stokes parameters reveal all obtainable polarization properties of the illuminated object, hence Stokes polarimeters can be applied for passive remote sensing. Mueller-matrix polarimeters, as taught for instance in U.S. Pat. Nos. 4,306,809, 5,247,176, and 5,956,147, employ a controlled polarized source, often a laser, to measure the Mueller matrix, which characterizes the complete linear response of a material or object to an electromagnetic wave. U.S. Pat. No. 7,218,398 teaches application of the data obtained with a Mueller-matrix polarimeter to an inverse problem for material characterization. Ellipsometers are a subclass of Mueller-matrix polarimeters that solve the inverse problem under the assumption of no scattering or birefringence, which results in a sparse Mueller matrix, although generalized ellipsometry, as taught for instance in U.S. Pat. No. 7,218,398, relaxes certain of these assumptions. Most polarimeters designed for characterization are laboratory instruments that also comprise a sample stage to accurately control the position and orientation of the material or object sample. Polarimeters can also be categorized according to their polarization modulators. The majority of polarimeters rely on some combination of polarizing prisms or beamsplitters and rotating retarder waveplates, although alternatives including photoelastic modulators (PEM), liquid-crystal variable retarders (LCVR), and microgrid arrays have also been employed. - Polarimeters for remote sensing must be compact, lightweight, and typically fast, and their sensing schemes must work for materials and objects in generally arbitrary locations and orientations within the field-of-view. The most common approach is to limit the number of measurements or channels to 2, achievable in real-time using a polarizing prism or beamsplitter in the receiver, and to form the difference of the 2 channels. This class of polarimeter for remote sensing is taught for instance in U.S. Pat. No. 8,116,000 for passive sources and in U.S. Pat. No. 7,333,897 for active sources. Because most polarizing prisms and beamsplitters separate EMR into orthogonal polarization states (eg, vertical and horizontal linear states), most polarimeters of this type measure 2 orthogonal channels. For active 2-channel polarimeters the channels are often aligned parallel to and orthogonal to (or crossed with) the polarization direction of the illuminating EMR. Certain active polarimeters, for instance those taught in U.S. Pat Nos. 5,956,147 and 6,060,710, employ photoelastic modulators (PEM), which require complex calibration and demodulation routines but can measure the Mueller matrix sufficiently fast for remote sensing.
- In remote sensing the terms detect, discriminate, classify, recognize, and identify have formal meanings, although the terms are casually interchanged in many publications. As used herein the term “detect” means to record EMR reflected or transmitted by the target; the term “discriminate” means to distinguish the target from the background; the term “classify” means to distinguish the target or targets among various classes including the background and other non-targets, which are often referred to as clutter; the terms “recognize” and “identify” imply even more specificity. Distinguishing materials (eg, metal, plastic, or wood) is generally considered classification, while distinguishing an object as a known combination of materials and possibly shape may be considered recognition. The variability of measured data, and especially field data, requires the application of learning algorithms to anticipate the distributions of measured signatures. Learning algorithms also help exclude clutter and reduce false-alarms. Algorithms can be adapted from the field of machine-learning. Learning algorithms of varying sophistication have been applied to classify materials based on polarimeter measurements, as taught for instance in U.S. Pat. Nos. 6,060,710 and 7,333,897 and by Jones, et al., Proc. SPIE 6240, 62400A (2006) and by Hoover and Tyo, Applied Optics 46, 8364-8373 (2007).
- Electromagnetic sensors that can rapidly and non-destructively discriminate, classify, and recognize materials, material states, and objects composed of those materials and states have many useful applications depending on the material sensitivities of the measured characteristics of the electromagnetic wave. Polarization at optical and near-optical frequencies is sensitive to broad ranges of materials and material states including metals, plastics, fibrous materials including fiber composites, crystals, chiral molecules, and biological tissues, and even broader ranges if patterned, stressed, or damaged materials are considered. Polarimeters have been demonstrated for the discrimination and classification of such materials and states, usually by increasing the contrast between a target material or material state and other materials or material states in the observation of reflected or transmitted EMR. Most polarimeters increase contrast by differencing irradiance channels defined by different settings of their polarization modulators. Ideally the irradiances of the target are different in the two channels, while the irradiances of non-targets are similar in the two channels, such that the difference channel provides increased target contrast. The polarimeter irradiance channels can be in the form of digital images or in the form of detector signals, multiplexed in time and/or space. For high-speed performance, polarimeters are usually limited to two simultaneous channels, as in the case of polarization-difference imaging (PDI), as taught for instance in U.S. Pat. Nos. 4,881,818, 5,929,443, and 8,116,000.
- Active polarimeters comprise a source of controlled, polarized EMR, often a laser, that is typically formed into a beam that illuminates the sample or scene of interest. Compared to passive polarimeters, which utilize solar illumination or thermal emission, active polarimeters provide access to additional signature components and perform independently of weather and time-of-day. Active polarimeters designed for characterization, as taught for optical and infrared frequencies in U.S. Pat. Nos. 4,306,809, 5,247,176, 5,956,147, and 7,218,398, typically measure the complete polarization signature, also known as the Mueller matrix, for illumination at given angles and electromagnetic frequencies. These are typically large laboratory instruments with stages to control sample orientation and position and polarization modulators with high size, weight, and power (SWaP) requirements and data rates well below video rate. For instance, U.S. Pat. No. 7,218,398 specifies the maximum data rate of this type of polarimeter on the order of 1 Hz.
- Mathematically, the Stokes parameters are written as
-
S=[S 0 S 1 S 2 S 3]T, (1) - where S0 is the irradiance, S1 is the degree of linear polarization in the horizontal-vertical coordinate system, S2 is the degree of linear polarization at ±45° to the horizontal-vertical coordinates, S3 is the degree of circular polarization, and [. . . ]T denotes a matrix transpose. The Stokes parameters define the polarization state of the EMR. The Mueller matrix is the linear transformation of the Stokes parameters upon interaction of the electromagnetic wave with a material object, written mathematically as
-
- For mathematical analysis and design the Mueller matrix is often written in one-dimensional form as
-
{tilde over (M)}=[M 00 M 01 M 02 M 03 M 10 . . . M 31 M 32 M 33]. (3) - The Stokes parameters are often referred to as the Stokes vector, which is correct only in the computing sense of a one-dimensional array of numbers. Since irradiance cannot be negative, both Stokes parameters and Mueller matrices lack additive inverses, so are not vectors in the algebraic sense.
- Active polarimeters designed for remote sensing typically do not measure, the complete Mueller matrix, but only parts of the signature that show contrast between the target of interest and the anticipated background and clutter, allowing SWaP to be reduced and speed increased due to fewer components and fewer measurements. For example, as taught in U.S. Pat. No. 7,580,127, most atmospheric aerosols and clouds exhibit a diagonal Mueller matrix with a single independent parameter that can be obtained from measurement of M11, which can be accomplished with an active polarimeter by forming the difference of two orthogonal channels. Nearly all active polarimeters designed for remote sensing are based on measurement of one or two polarization channels, with this limitation imposed by requirements for low SWaP and/or high speed Vannier, et al., Applied Optics 55, 2881-2891 (2016) demonstrates the current public-domain state-of-the-art in polarimeters for optical remote sensing, employing LCVR modulators and an optimization routine to empirically derive and implement a single-channel active polarimeter that provides enhanced-contrast imaging of a specific class of materials. However, as shown in Hoover and Tyo, Applied Optics 46, 8364-8373 (2007), active polarization signatures of many relevant materials are actually multi-dimensional, and discrimination, classification, and recognition performance can be improved by applying learning algorithms to measured multi-dimensional signatures. In particular, multi-dimensional measurements are usually needed to achieve classification performance that is relatively invariant to the target orientation or pose, which is usually a critical requirement for remote sensing. Fielding an active polarimeter for high-speed, high-performance classification or ATR therefore requires a new class of polarimeter that can measure specific multi-dimensional signatures while maintaining the low SWaP and high speed achievable by 2-channel polarimeters.
- Most active field sensors, including polarimeters, are configured for monostatic or near-monostatic geometries, due to the practical advantage of placing the source-transmitter and the receiver-detector on the same platform. Passive sensors, on the other hand, often operate at large bistatic angles, for instance as illustrated in
FIG. 2 . While the dimensionally of passive polarization signatures is generally too low to achieve pose invariance or match the classification performance of active polarimeters, for certain materials and objects the passive signature contains unique information at large bistatic angles that can augment the active signature and improve performance. There is therefore a need for field polarimeters to combine active and passive signature measurements in the same hybrid polarimeter in order to utilize more information and improve remote-sensing performance. - One embodiment of the present invention provides for a registered-channel multiplexer (RCM) having a scanning element positioned to direct a beam from a controlled source of electromagnetic radiation (EMR) serially among a set of independent modulators. One version of the RCM is illustrated in
FIG. 3 . Each member of the set of independent modulators defines a channel that is configured to modulate a state of the beam among a set of independent states to produce a modulated beam. An assembly of reflectors is configured to redirect the modulated beam to a co-registered path that is pointed at a field-of-view. The scanning element can be a scanner mirror attached to a first mechanical scanner, for instance a galvanometer scanner, that comprises a rotation axis, and the assembly of reflectors can be an assembly of mirrors. Further still the scanner mirror can be a composite mirror comprising a flat mirror parallel to the rotation axis of the first mechanical scanner and a second wedge mirror mounted at an angle to the flat mirror, for example, at an angle of about 45 degrees. The assembly of mirrors may comprise one or more flat scanner mirrors attached to a second set of mechanical scanners. The set of independent modulators can be selected from the group consisting of lenses, polarization modulators, spectral filters, diffractive array generators, holograms, amplitude-phase masks and any combination thereof. One or more members of the set of independent modulators can be tunable, and one or more different members can be a beam block as illustrated inFIG. 4 . - Another embodiment of the present invention provides for a sensor for classifying or recognizing a target or targets within a field-of-view (FOV) using multiple independent channels. The sensor comprises a source of controlled electromagnetic radiation (EMR) positioned in a path with a detector and the FOV positioned there between. The sensor includes a registered-channel multiplexer (RCM) comprising a scanning element positioned to direct a beam from the source of controlled EMR serially among a first set of independent modulators, wherein each member of the first set of independent modulators corresponds to a channel and the scanner and modulators are configured to modulate a state of the beam among a set of independent states to produce a modulated beam. An assembly of reflectors is configured to redirect the modulated beam to a co-registered path that is pointed at the FOV. An EMR collector is positioned to direct EMR reflected from or transmitted by objects within the FOV to a second modulator independent of the first set of independent modulators. The detector is positioned to receive the EMR from the second modulator, wherein the detector produces a set of signals that are synchronized with the set of channels formed by the first set of modulators and the second modulator. The first set of independent modulators and the second independent modulator can be selected from the group consisting of lenses, polarization modulators, spectral filters, diffractive array generators, holograms, amplitude-phase masks, and any combination thereof. The sensor includes a processor that is connected with a memory, wherein the processor is configured to execute a classification algorithm stored in the memory by applying a user-adjustable threshold to the set of detector signals to assign a class label (which may also be stored in memory) to the set of detector signals and wherein the classification algorithm specifies the function of the threshold in terms of the set of channels formed by the first set of modulators and the second modulator. The class label can be either binary, assuming the
values 0 or 1 only, or multi-valued. The first set of independent modulators and the second independent modulator can be tunable, and one or more different members of the first set of independent modulators can be a beam block. The EMR directed by the collector to the second modulator may include EMR from a passive source illuminating or emanating from objects within the FOV, in which case the EMR directed by the collector to the second modulator while the scanning element is directed at the beam block can produce a detector signal corresponding to a passive channel as a member of the set of detector signals. In one embodiment the set of detector signals is a set of registered digital images and a class label is assigned to each pixel to form a digital classification image. The class labels can be assigned at a rate of 20 Hz or greater and for imaging the class labels can be assigned and the classification image electronically displayed at a rate of 20 frames-per-second or greater. The FOV can be scanned or swept over an area or a volume by a gimbal or a moving platform on which the sensor is mounted. Further, the set of detector signals can be synchronized with the set of channels formed by the first set of modulators and the second modulator by a lock-in amplifier using a timing signal provided by the scanning element. - Yet another embodiment of the present invention provides for a polarimeter for classifying or recognizing a target or targets within a field-of-view (FOV) using multiple independent polarization channels. The polarimeter comprises a source of controlled electromagnetic radiation (EMR) positioned in a path with a detector and the FOV positioned there between. The polarimeter includes a registered-channel multiplexer (RCM) comprising a scanning element positioned to direct a beam from the source of controlled EMR serially among a first set of independent polarization modulators. The first set of independent polarization modulators, wherein each member of the set corresponds to a polarization channel, is configured to modulate the polarization state of the beam among a set of independent polarization states to produce a modulated beam. An assembly of reflectors is configured to redirect the modulated beam to a co-registered path that is pointed at the FOV. An EMR collector is positioned to direct EMR reflected from or transmitted by objects within the FOV to a second polarization modulator independent of the first set of polarization modulators. The detector is positioned to receive the EMR from the second polarization modulator, wherein the detector produces a set of signals that are synchronized with the set of channels formed by the first set of polarization modulators and the second polarization modulator. The polarimeter includes a processor connected with a memory, wherein the processor is configured to execute a classification algorithm stored in the memory by applying a user-adjustable threshold to the set of detector signals to assign a class label (which may also be stored in memory) to the set of detector signals, wherein the classification algorithm specifies the function of the threshold in terms of the set of channels formed by the first set of polarization modulators and the second polarization modulator. The class label can be either binary, assuming the
values 0 or 1 only, or multi-valued. The first set of independent polarization modulators and the second independent polarization modulator can be tunable, and one or more member of the first set of polarization modulators can be a beam block. The EMR directed by the collector to the second modulator may include EMR from a passive source illuminating or emanating from the field-of-view, in which case the EMR directed by the collector to the second modulator while the scanning element is directed at the beam block can produce a detector signal corresponding to a passive channel as a member of the set of detector signals. In one embodiment the set of detector signals is a set of registered digital images and a class label is assigned to each pixel to form a digital polarization classification image (PCI). The class labels can be assigned at a rate of 20 Hz or greater and PCIs can be electronically displayed at a rate of 20 frames-per-second or greater. The FOV can be scanned or swept over an area or a volume by a gimbal or a moving platform on which the polarimeter is mounted. Further, the set of detector signals can be synchronized with the set of channels formed by the first set of polarization modulators and the second polarization modulator by a lock-in amplifier using a timing signal provided by the scanning element. In one embodiment the channels are defined by a subspace-projection algorithm that applies a gradient operator on a vector space of reduced Mueller matrices and maximizes a resulting contrast parameter, and the user-adjustable threshold is a function of the channels derived by the subspace-projection algorithm. In another embodiment the contrast parameter and the function of the user-adjustable threshold are defined by an external machine-learning algorithm. - Yet another embodiment of the present invention provides a method for classifying or recognizing a target or targets within a field-of-view (FOV) by using a polarimeter with multiple independent channels. The method includes positioning a source of controlled electromagnetic radiation (EMR) in a path with a detector and the FOV positioned there between. The method employs a registered-channel multiplexer (RCM), which uses a scanning element to serially direct a beam from the controlled source of EMR among a first set of independent polarization modulators, producing a modulated beam, wherein each member of the set of independent polarization modulators corresponds to a channel and the set is configured to modulate the polarization state of the beam among a set of independent polarization states, and an assembly of reflectors to redirect the modulated beam to a co-registered path pointed at the FOV. The method further includes collecting a portion of the EMR reflected from or transmitted by objects within the FOV with an EMR collector that directs the EMR to a second polarization modulator independent of the first set of polarization modulators, and receiving the EMR from the second polarization modulator at the detector, wherein the detector produces a set of detector signals that are synchronized with the set of channels formed by the first set of polarization modulators and the second polarization modulator. The method applies a classification algorithm and a user-adjustable threshold to the set of detector signals to assign a class label to the set of detector signals, wherein the classification algorithm specifies the function of the threshold in terms of the set of channels formed by the first set of polarization modulators and the second polarization modulator and the class label can be either binary, assuming the
values 0 or 1 only, or multi-valued. The first set of independent polarization modulators and the second independent polarization modulator can be tunable, and one or more members of the first set of polarization modulators can be a beam block. The EMR directed by the collector to the second modulator may include EMR from a passive source illuminating or emanating from objects within the FOV, in which case the EMR directed by the collector to the second modulator while the scanning element is directed at the beam block can produce a detector signal corresponding to a passive channel as a member of the set of detector signals. In one embodiment the set of detector signals is a set of registered digital images and a class label is assigned to each pixel to form a digital polarization classification image (PCI). The class labels can be assigned at a rate of 20 Hz or greater and PCIs can be electronically displayed at a rate of 20 frames-per-second or greater. The method may also include scanning or sweeping the FOV through an area or a volume using a gimbal or a moving platform on which the polarimeter is mounted. The method may also include synchronizing the set of detector signals with the set of channels formed by the first set of polarization modulators and the second polarization modulator using a lock-in amplifier that reads a timing signal from the scanning element. One embodiment includes defining the polarization channels using a subspace-projection algorithm that applies a gradient operator on a vector space of reduced Mueller matrices and maximizes a resulting contrast parameter, and the user-adjustable threshold can be a function of the channels derived by the subspace-projection algorithm. Another embodiment includes applying an external machine-learning algorithm to define the contrast parameter and the function of the user-adjustable threshold. - Another embodiment comprises a hybrid polarimeter and method of use that can achieve data rates up to video-rate in imaging mode with sufficiently small SWaP for field implementations.
- The new subspace-projection algorithm specifies the set of polarization-modulator settings or channels that maximize contrast between the target or state of interest and the anticipated backgrounds and other objects or states that might be observed. Each polarimeter channel corresponds to an irradiance measurement or a detector signal, and measurements in multiple channels are generally combined to form a set of detector signals or a projection that is used to test whether the object or state corresponds to the target or not. The projection can be either a one-dimensional or multi-dimensional irradiance signal depending on the type of classification algorithm used. For a one-dimensional projection a simple irradiance threshold serves as the classifier, such that an irradiance above or below the threshold, or between two thresholds, corresponds to a positive target classification. In imaging mode every pixel is an independent sensor, and positive classifications can be represented by white pixels with all other pixels black, providing binary polarization-classification images (PCI) that are advantageous for fast interpretation and practical data transmission and storage.
- The subspace-projection algorithm derives the polarimeter channels from the available polarization modulators and the Mueller matrices of the anticipated objects, materials, and/or material states, with the Mueller matrices known a priori from a model and/or training characterization measurements. The subspace-projection algorithm provides the maximally-biased one-dimensional classifier if only a single target. Mueller matrix is used as input, and provides more generalized classifiers if non-target Mueller matrices are incorporated into the input (see Detailed Description.) The method can be applied to an existing polarimeter, to derive the optimal set of realizable channels, or at the polarimeter design stage to derive the minimum set of channels that can achieve the required classification or ATR performance.
- Another embodiment of the present invention provides a method of enhancing the contrast of a target or targets recorded by a polarimeter. This method includes obtaining modeled or measured Mueller matrices of the target or targets and other materials or objects to be recorded by the polarimeter and applying the subspace-projection algorithm to the Mueller matrices and to the set of polarization modulators of the polarimeter to specify the settings of the polarization modulators, wherein the subspace-projection algorithm applies a gradient operator on a vector space of reduced Mueller matrices and maximizes a resulting contrast parameter. Another embodiment includes applying an external machine-learning algorithm to define the contrast parameter.
- The current invention includes new hardware designs for a polarization multiplexer to achieve data rates up to video-rate in imaging mode. While passive polarimeters contain polarization modulators in the receiver only, active polarimeters comprise modulators in both the transmitter and receiver. One embodiment of the current invention employs a conventional polarization modulator, a polarizing prism, in the receiver, and a new temporally-multiplexed polarization modulator (or polarization multiplexer) in the transmitter, which is a form of the more general registered-channel multiplexer (RCM.) The RCM can assume various geometries, all of which utilize a scanner to switch or scan the beam among multiple independent and tunable states at high speed, possibly including a null state that corresponds to a passive channel. Achieving video-rate multiplexing with a mechanical scanner requires a rotor of lower inertia than the hollow-core motors used in rotating-retarder polarimeters like those taught in U.S. Pat. Nos. 4,306,809, 5,247,176, and 7,218,398. One embodiment of the invention employs a low-inertia galvanometer scanner that switches a laser beam among a set of tunable polarization modulators. Each modulator in the transmitter, in combination with a static polarizing prism in the receiver, constitutes 2 orthogonal polarimeter channels, while overall the polarimeter may have 8 or more channels. Each polarization modulator in the transmitter can be independently tuned to realize multiple arbitrary channels.
- Classification of various materials on natural backgrounds has been demonstrated in field tests, at approximately 20 m horizontal range, employing the polarimeter main embodiment of the invention. Applications of the main embodiment include, but are not limited to, classification of hazardous objects for route clearance, classification of materials in mining, manufacturing, waste, and recycling streams, and real-time non-destructive testing of vehicle and structural materials including, but not limited to, fiber composites and non-cubic metals. By changing the laser frequency and/or the projection and collection optics, the invention can be customized for manufacturing applications including, but not limited to, monitoring of thin-film, textile, polymer, chiral pharmaceutical, composite material, and non-cubic metal fabrication and processing, including mechanical and laser machining. Smaller embodiments of the invention are applicable to recognition of security inks for asset tagging and product authentication, quantification of micro-plastics, biomedical-tissue assays, and general polarization microscopy, while larger embodiments that utilize high-power lasers are applicable for long-range material and object classification and recognition from resolved or unresolved imagery, for instance surveillance and reconnaissance, damage assessment, search and rescue, mapping and, assessment of contaminants and debris, particularly plastics, and space-object classification.
- The accompanying drawings, which are incorporated into and form a part of the specification, illustrate one or more embodiments of the present invention and, together with the description, serve to explain the principles of the invention. The drawings are only for the purpose of illustrating one or more embodiments of the invention and are not to be construed as limiting the invention.
-
FIG. 1 is a taxonomy of polarimeter types based on electromagnetic-radiation (EMR) source, measurements made, polarization modulators employed, and applications. -
FIG. 2 is an illustration of active monostatic, passive bistatic, and hybrid active-passive sensor geometries. -
FIG. 3 is an illustration of the polarization multiplexer of an embodiment of the present invention that provides 8 active tunable polarization channels or 4 independent active tunable polarization-difference channels. -
FIG. 4 is an illustration of the polarization multiplexer of an embodiment of the present invention that provides 6 active tunable polarization channels and 2 passive tunable polarization channels. -
FIG. 5 is an illustration of the composite wedge scanner mirror employed in an embodiment of the present invention. -
FIG. 6A-B is a comparison of a conventional image and the corresponding polarization-classification image (PCI) of plastic bags on a beach obtained in field tests of an embodiment of the invention at 15 m horizontal range. -
FIG. 7 is a flowchart of an embodiment of the invention designed to map the distribution of a selected material at a test site. -
FIG. 8A-D is a set of polarization-classification images (PCI) of plastic bags on a beach obtained using different classifier thresholds. - The invention can be described in detail, with reference to the figures, by first describing the subspace-projection algorithm used to derive polarimeter channels and corresponding classifiers from the available polarization modulators and a priori Mueller matrices, then describing the specific hardware by which the new polarimeter design can be reduced to practice in a hybrid polarimeter that can achieve data rates up to video-rate in imaging mode with sufficiently small size, weight, and power (SWaP) to allow field implementations.
- The subspace-projection algorithm derives what is known in the technical literature as a partial Mueller-matrix polarimeter (pMMP.) The newly disclosed subspace-projection algorithm has advantages over previously published pMMP design methods in that it avoids use of matrix pseudoinverses, which are prone to errors, and imposes no limitations on the objects, materials, or material states that can be observed by the polarimeter. The subspace-projection algorithm is based on application of vector algebra on the vector space of reduced Mueller matrices formed by reducing the measured irradiance as
-
I′=I−I U, (4) - where I is the irradiance actually measured by the polarimeter and IU is the irradiance that a conventional non-polarizing camera would measure with unpolarized illumination. The reduced irradiance I′ can be positive or negative. It can be shown that the reduced irradiance is expressible as a vector product of the reduced polarimeter instrument matrix and the reduced Mueller matrix of the illuminated object
-
{tilde over (M)}′=[M 01 M 02 M 03 M 10 . . . M 31 M 32 M 33]. (5) - Construction of reduced Mueller matrices is necessary because conventional Mueller matrices do not form an algebraic vector space. The general pMMP design problem is then stated as follows: Given a target Mueller matrix MT and a polarimeter configuration defined by a parameter set Θ={θ1,θ2, . . . ,θQ}, determine the set Θ that maximizes representation of the target matrix in the set of M×N polarimeter irradiance measurements {Imdcn}, m=1, . . . , M and n=1, . . . , N. In one embodiment of the invention the parameters Θ are retarder-waveplate angular orientations, but Θ can represent any combination of polarization-modulator settings. The algorithm is also general enough to derive channels and classifiers for passive polarimeters. In order to maximize target contrast in a polarimeter measurement, the polarimeter instrument matrix must be chosen to maximize its projection onto the subspace of the reduced target matrix {tilde over (M)}′T. This subspace-projection can be formalized by defining a differentiable target matrix and applying the gradient operator on the space of reduced Mueller matrices, defined as
-
- to the reduced irradiances. The representation of the target in the polarimeter measurement is then maximized by maximizing the contrast parameter, which is a vector dot product, over the set of modulator parameters Θ, that is
-
- This new polarimeter system and method of use has been validated with field data collected using the main embodiment of the invention. This embodiment can be described in more detail with reference to the figures.
- Field sensors that utilize different electromagnetic-radiation (EMR) sources are illustrated in
FIG. 2 . The sensor field-of-view (FOV) may be scanned or swept over an area or volume, for instance if the sensor is mounted to a vehicle, or objects may move through a fixed FOV, for instance on a conveyor belt. Each recorded frame may or may not contain a target, although for the purposes of illustration a target is assumed inFIG. 2 . Thetarget 1 is illuminated by theprobe beam 2, typically a laser beam, that is projected by thesensor transmitter 3. The target may also be illuminated by apassive source 4, which is typically the sun but can also be thermal radiation from the object itself if infrared frequencies are recorded. For active sensors the reflectedprobe EMR 5 is recorded by thesensor receiver 6, while for passive sensors thereceiver 7 records only reflectedpassive EMR 8. A hybrid active-passive sensor illuminates the target with theprobe beam 9 projected by thetransmitter 10 and records both the reflected probe EMR 11 and reflected passive EMR 12 in thereceiver 13. - New hardware designs as discussed below allow the channels specified by the subspace-projection algorithm to be implemented and corresponding data to be obtained at video rate. The polarization multiplexer is enclosed in the sensor transmitter, embodiments of which are illustrated at 3 and 10 in
FIG. 2 . The polarization multiplexer is a form of registered-channel multiplexer (RCM). Birefringent crystal waveplates are assumed for the polarization modulators, although other polarization modulators can be used. Conventional polarimeters employ waveplates in rotating hollow-core motors, which can be very accurate but not fast enough to achieve polarization video. As illustrated inFIG. 3 , an embodiment of the present invention employs low-inertia scanners, for instance galvanometer scanners, to switch among multiple arms that comprise independent tunable polarization modulators. - Referring to
FIG. 3 , in one embodiment alaser 14 emits a collimated beam of light that strikes amirror 15 that is scanned by agalvanometer scanner 16 serially among multiple independent arms 17 a-d, each of which, in combination with the static polarizing prism in the receiver (not shown), defines a pair of orthogonal channels. During the initial time interval the scanner directs the beam through thefirst arm 17 a, which comprises thefirst polarization modulator 18, which is set to a channel specified by the subspace-projection algorithm. Themodulator 18 may be a rotatable crystal waveplate or a pair of waveplates, but can be other types of polarization modulator as well. Following the polarization modulator the beam is reflected by a fold-mirror assembly 19 back to thescanner mirror 15 and then to thefield mirror 20, which reflects the beam toward the field-of-view. The beam is generally expanded or focused by projector lenses located atbeam position 21. The scanner dwells on thefirst arm 17 a for a predetermined time interval, then switches to thesecond arm 17 b for a second time interval, generally of the same duration as the first time interval, then switches through the remainingarms 17 c-d before rewinding to thefirst arm 17 a and repeating the scan. The subsystem illustrated inFIG. 3 is generally called a polarization-state generator (PSG) and may have more or fewer than 4 arms. - Generally the set of irradiances or images collected for the set of scanner positions, for example 4 positions for the embodiment illustrated in
FIG. 3 , forms the projection vector used for classification. The dwell time on each arm is determined by the number of arms, the EMR power collected from the FOV, the detector sensitivity, and the sensor speed requirement, and is generally limited by the scanner switching speed. In addition to dwell time the scanner has a characteristic switching time. A typical sum of dwell and switching times with 4 arms and a high-performance galvanometer scanner is 12.5 ms, which implies 50 ms for collection of each image or projection vector. With sufficiently fast data-processing this enables PCI video at 20 frames-per-second. Faster scanners would enable higher data rates. -
FIG. 4 illustrates a PSG with a null channel recorded while the scanner dwells on abeam block 22. The null channel records a pair of passive channels that are combined with the active channels obtained from the remaining arms to realize a hybrid active-passive polarimeter. The dwell time on the null channel may be the same, longer, or shorter than the dwell times on the active channels if the reflected power from the passive source is less than or greater than the reflected laser power. The scanner timing signal can provide a reference frequency for lock-in detection. -
FIG. 5 illustrates the composite scanner mirror of an embodiment of the present invention. The main parts of the composite scanner mirror are the mountingboss 23, theflat mirror 24, and thewedge mirror 25. Themirrors - The classification performance and data-processing speed depend on the classification algorithm, a great variety of which may be employed. The classification algorithm, or simply classifier, can be one-dimensional or multi-dimensional, corresponding to the number of measured channels. Two-channel polarimeters typically use one-dimensional classifiers based on the difference between the measured channels, hence the term polarization-difference imaging (PDI.) A one-dimensional classifier makes a positive target classification when the irradiance projection measured in a specified channel is above or below a user-defined threshold, or between two thresholds, and a negative classification otherwise. Classifiers can be defined in terms of the channels derived by the subspace-projection algorithm of the present invention.
- Classifiers are described in terms of the bias-variance tradeoff, which balances performance on training datasets with the ability to generalize to unforeseen data. The smallest possible training dataset comprises the target only and results in a high-bias classifier that can perform well on scenes containing the target only but may suffer frequent false-alarms when clutter objects are observed. A one-dimensional classifier based on a channel derived from Eq. 7 of the subspace-projection algorithm is a high-bias classifier. Variance is the ability of the classifier to perform well on scenes that contain target variations and/or background or clutter not represented in the training dataset. Increasing the classifier dimensionally can increase the variance as the target signature can be projected into a higher dimensional space where its projection is more unique. The classifier dimension can be increased by combining two or more channels derived from Eq. 7 of the subspace-projection algorithm. Alternatively or in combination, the classifier bias can be lowered by using a larger training dataset that includes signatures of the background and/or anticipated clutter objects. For example, the background signature can be incorporated by maximizing the generalized contrast parameter R,
-
- where the denominator quantifies representation of the background signature in the measured channel, rather than the simple target contrast of Eq. 7. Many other versions of generalized contrast parameter based on channels derived by the subspace-projection algorithm are possible incorporating the background and any number of clutter polarization signatures. It will also be obvious to those skilled in the art of learning algorithms that the subspace-projection algorithm of the current invention can be combined with established algorithms including but not limited to principal-components analysis (PCA) and support-vector machines (SVM) to further generalize the classifier.
- In image format the classifier is applied to each pixel independently. Positive classifications can be represented by white pixels with all other pixels black, resulting in binary polarization classification images (PCI) as illustrated in
FIG. 6 . Binary images are advantageous for fast interpretation and practical data transmission and storage. In the conventional unpolarized image ofFIG. 6 the target, theplastic bags 26, are nearly indiscernible from the algae, twigs, and other clutter objects 27, especially if scanning over a large area at video rate, while the PCI reveals the target clearly. The PCI shown inFIG. 6 was obtained in field tests with the main embodiment of the invention in full daylight at 15 m horizontal range. -
FIG. 7 is a flowchart of an embodiment of the invention designed to map the distribution of a selected material at a test site using a classifier derived by the subspace-projection pMMP-design algorithm. The components inside the dashedbox 28 are on-board the field sensor and the flowchart corresponds to a single image or video frame. The sensor field-of-view may be scanned or swept over the test area by a motion system not shown. The subspace-projection pMMP-design algorithm defines the classifier, which is stored incomputer memory 29, and the settings of the polarization modulators, which can be changed by the user via thetuner 30. ThePSG scanner 31 switches the beam among the independent modulators in the RCM and provides a timing signal to trigger thecamera 32 to record corresponding images in thechannels 33 defined by the subspace-projection algorithm. These images and the corresponding channel specifications constitute thetest data 34 that is input to a real-time processor such as a general-purpose graphics-processing unit (GPGPU) 35. A parallel-processing architecture such as GPGPU is advantageous since each pixel is processed as an independent sensor, although other real-time architectures may be employed. The user can adjust the classifier threshold via thecontrol 36. From these inputs the processor produces the class labels 37, which are 0 or 1 at each pixel for binary PCI. If the target is present in the frame, that is if the frame contains a sufficient number of white pixels, then the frame is stored on the site map with the corresponding global-positioning system (GPS) coordinates provided by anexternal GPS module 38. With sufficient memory all frames can be stored on the site map. Alternatively each frame can be assigned a class label, 0 or 1, to produce a density map without images. The GPS module and stored site map may also be on-board the sensor. -
FIG. 8A-D is a set of binary PCIS of the scene illustrated inFIG. 6A obtained by the main embodiment of the invention as the classifier threshold is adjusted. - Both the classifier training and the polarimeter testing and ordinary use are executed by a general or specific-purpose computer or distributed system programmed with computer software implementing the steps described above, which computer software may be in any appropriate computer language, including C++, Python, FORTRAN, BASIC, Java, assembly language, distributed programming languages, etc. The apparatus may also include a plurality of such computers/distributed systems (e.g., connected over the Internet and/or one or more intranets) in a variety of hardware implementations. For example, data processing can be performed by an appropriately programmed microprocessor, computing cloud, field-programmable gate array (FPGA), general-purpose graphics-processing unit (GPGPU), or the like, in conjunction with appropriate memory, network, and bus elements. All computer software disclosed herein may be embodied on any computer-readable medium (including combinations of mediums), including without limitation CD-ROMs, DVD-ROMs, hard drives (local or network storage device), USB keys, other removable drives, ROM, and firmware.
- Note that in the specification and the claims the words “a”, “an”, and “the” mean one or more unless otherwise specified, and the words “about” or “approximately” mean within twenty percent (20%) of the numerical value cited.
- Although the invention has been described in detail with particular reference to these embodiments, other embodiments can achieve the same results. Variations and modifications of the present invention will be obvious to those skilled in the art and it is intended to cover all such modifications and equivalents. The entire disclosures of all references, applications, patents, and publications cited above and/or in the attachments, and of the corresponding application(s), are hereby incorporated by reference.
Claims (62)
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2016/068411 WO2018118073A1 (en) | 2016-12-22 | 2016-12-22 | Polarimeter with multiple independent tunable channels and method for material and object classification and recognition |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2016/068411 A-371-Of-International WO2018118073A1 (en) | 2016-12-22 | 2016-12-22 | Polarimeter with multiple independent tunable channels and method for material and object classification and recognition |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US16/748,463 Continuation-In-Part US11867891B2 (en) | 2016-12-22 | 2020-01-21 | Polarimeter with multiple independent tunable channels and method for material orientation imaging |
Publications (3)
Publication Number | Publication Date |
---|---|
US20180181839A1 US20180181839A1 (en) | 2018-06-28 |
US20190073561A2 true US20190073561A2 (en) | 2019-03-07 |
US10540571B2 US10540571B2 (en) | 2020-01-21 |
Family
ID=62626931
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US15/542,064 Active US10540571B2 (en) | 2016-12-22 | 2016-12-22 | Polarimeter with multiple independent tunable channels and method for material and object classification and recognition |
Country Status (5)
Country | Link |
---|---|
US (1) | US10540571B2 (en) |
EP (1) | EP3371631B1 (en) |
JP (1) | JP6748648B2 (en) |
AU (2) | AU2016388716A1 (en) |
WO (1) | WO2018118073A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021150275A1 (en) * | 2020-01-21 | 2021-07-29 | Advanced Optical Technologies, Inc. | Polarimeter with multiple independent tunable channels and method for material orientation imaging |
US11867891B2 (en) | 2016-12-22 | 2024-01-09 | Advanced Optical Technologies, Inc. | Polarimeter with multiple independent tunable channels and method for material orientation imaging |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2015115041A (en) * | 2013-12-16 | 2015-06-22 | ソニー株式会社 | Image processor, and image processing method |
CN108801930B (en) * | 2018-05-30 | 2020-09-08 | 华中科技大学 | High-time-resolution muller matrix ellipsometry measurement device and method |
JP7435464B2 (en) * | 2018-11-21 | 2024-02-21 | ソニーグループ株式会社 | Workpiece discrimination system, workpiece discrimination device, and workpiece discrimination method |
CN111670372A (en) * | 2019-01-09 | 2020-09-15 | 深圳市大疆创新科技有限公司 | Scanning module, distance measuring device and mobile platform |
WO2020202221A1 (en) * | 2019-04-04 | 2020-10-08 | Lovely Professional University | A method and system for detecting and modeling objects in space using lidar |
CN110163303B (en) * | 2019-06-03 | 2022-01-11 | 中国农业大学 | Grid-based remote sensing image parallel classification method and system |
US11499911B2 (en) * | 2019-11-15 | 2022-11-15 | Montana State University | Systems and methods for detecting thermodynamic phase of clouds with optical polarization |
WO2022015621A1 (en) * | 2020-07-13 | 2022-01-20 | Board Of Regents, The University Of Texas System | Devices, systems, and methods for chiral sensing |
CN114419452B (en) * | 2022-03-30 | 2022-07-05 | 中国人民解放军火箭军工程大学 | High-resolution dual-polarization SAR anti-corner reflector interference target identification method |
Family Cites Families (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4306809A (en) | 1979-03-26 | 1981-12-22 | The Board Of Regents Of The University Of Nebraska | Polarimeter |
US4688939A (en) * | 1985-12-27 | 1987-08-25 | At&T Technologies, Inc. | Method and apparatus for inspecting articles |
US4881818A (en) | 1986-04-22 | 1989-11-21 | The University Of New Mexico | Differential imaging device |
JP2969009B2 (en) * | 1991-02-22 | 1999-11-02 | 株式会社リコー | Axial mirror deflector |
US5247176A (en) | 1991-07-22 | 1993-09-21 | The United States Of America As Represented By The Secretary Of The Air Force | Infrared laser polarimeter |
US5929443A (en) | 1995-12-18 | 1999-07-27 | The Research Foundation City College Of New York | Imaging of objects based upon the polarization or depolarization of light |
EP0902885A4 (en) * | 1996-05-16 | 2006-09-27 | Affymetrix Inc | Systems and methods for detection of labeled materials |
US5956147A (en) | 1997-06-13 | 1999-09-21 | Lockheed Martin Energy Research Corporation | Two modulator generalized ellipsometer for complete mueller matrix measurement |
US6317170B1 (en) * | 1997-09-13 | 2001-11-13 | Samsung Electronics Co., Ltd. | Large screen compact image projection apparatus using a hybrid video laser color mixer |
WO1999038046A1 (en) | 1998-01-23 | 1999-07-29 | Burger Robert J | Lenslet array systems and methods |
US6060710A (en) | 1998-12-21 | 2000-05-09 | The United States Of America As Represented By The Secretary Of The Army | Infrared Mueller matrix detection and ranging system |
US6351324B1 (en) | 2000-03-09 | 2002-02-26 | Photera Technologies, Inc. | Laser imaging system with progressive multi-beam scan architecture |
JP2002151774A (en) * | 2000-09-01 | 2002-05-24 | Nidek Co Ltd | Laser equipment |
JP2002148554A (en) | 2000-11-03 | 2002-05-22 | Samsung Electronics Co Ltd | Optical scanner and laser video projector applying the same and driving method for the same |
US6927888B2 (en) | 2002-05-13 | 2005-08-09 | Juan Manuel Bueno Garcia | Method and apparatus for imaging using polarimetry and matrix based image reconstruction |
US7420675B2 (en) * | 2003-06-25 | 2008-09-02 | The University Of Akron | Multi-wavelength imaging system |
EP1660849A2 (en) | 2003-08-06 | 2006-05-31 | Russell Chipman | Advanced polarization imaging method, apparatus, and computer program product for retinal imaging, liquid crystal testing, active remote sensing, and other applications |
US7218398B2 (en) | 2003-11-26 | 2007-05-15 | Smith Matthew H | Method and apparatus for determining liquid crystal cell parameters from full Mueller matrix measurements |
JP2005181927A (en) | 2003-12-24 | 2005-07-07 | Canon Inc | Optical deflection device and image display device using same |
US7333897B2 (en) | 2005-12-07 | 2008-02-19 | Motorola, Inc. | Method and system for identifying material composition based upon polarization trajectories |
US8116000B2 (en) | 2006-08-18 | 2012-02-14 | Q5 Innovations Inc. | Tandem polarization beam splitters for real-time polarization difference sensing |
US20090116518A1 (en) | 2007-11-02 | 2009-05-07 | Pranalytica, Inc. | Multiplexing of optical beams using reversed laser scanning |
EP2225699A1 (en) | 2007-11-27 | 2010-09-08 | University of Washington | Eliminating illumination crosstalk while imaging using multiple imaging devices |
US8570650B2 (en) * | 2009-12-09 | 2013-10-29 | Applied Precision, Inc. | Method and system for fast three-dimensional structured-illumination-microscopy imaging |
DE102011013613A1 (en) * | 2010-10-01 | 2012-04-05 | Carl Zeiss Microimaging Gmbh | Microscope and microscopy method |
EP2682776B1 (en) * | 2012-07-07 | 2017-06-28 | Airbus DS Electronics and Border Security GmbH | Method for detecting pulsed laser radiation and imaging laser detector |
DE102012017922B4 (en) * | 2012-09-11 | 2024-03-14 | Carl Zeiss Microscopy Gmbh | Optical arrangement and light microscope |
US9983260B2 (en) * | 2012-10-12 | 2018-05-29 | The United States Of America As Represented By The Secretary Of The Air Force | Dual-phase interferometry for charge modulation mapping in ICS |
DE102013227105A1 (en) | 2013-09-03 | 2015-03-05 | Leica Microsystems Cms Gmbh | Microscope and acousto-optical beam combiner for a microscope |
-
2016
- 2016-12-22 JP JP2017537482A patent/JP6748648B2/en not_active Expired - Fee Related
- 2016-12-22 US US15/542,064 patent/US10540571B2/en active Active
- 2016-12-22 EP EP16920449.2A patent/EP3371631B1/en active Active
- 2016-12-22 WO PCT/US2016/068411 patent/WO2018118073A1/en active Application Filing
- 2016-12-22 AU AU2016388716A patent/AU2016388716A1/en not_active Abandoned
-
2018
- 2018-08-29 AU AU2018222946A patent/AU2018222946A1/en not_active Abandoned
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11867891B2 (en) | 2016-12-22 | 2024-01-09 | Advanced Optical Technologies, Inc. | Polarimeter with multiple independent tunable channels and method for material orientation imaging |
WO2021150275A1 (en) * | 2020-01-21 | 2021-07-29 | Advanced Optical Technologies, Inc. | Polarimeter with multiple independent tunable channels and method for material orientation imaging |
Also Published As
Publication number | Publication date |
---|---|
US10540571B2 (en) | 2020-01-21 |
EP3371631A4 (en) | 2019-08-14 |
JP2019509463A (en) | 2019-04-04 |
WO2018118073A1 (en) | 2018-06-28 |
US20180181839A1 (en) | 2018-06-28 |
EP3371631B1 (en) | 2021-08-18 |
AU2018222946A1 (en) | 2018-09-20 |
JP6748648B2 (en) | 2020-09-02 |
AU2016388716A1 (en) | 2018-07-12 |
EP3371631A1 (en) | 2018-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10540571B2 (en) | Polarimeter with multiple independent tunable channels and method for material and object classification and recognition | |
Tyo et al. | Review of passive imaging polarimetry for remote sensing applications | |
US6353673B1 (en) | Real-time opto-electronic image processor | |
EP2223083B1 (en) | Sensing using polarization diversity and wavelength dependent backscatter | |
US20190325586A1 (en) | Method and device of measuring infrared spectral characteristics of moving target | |
US11860034B2 (en) | Systems, devices, and methods for hyperspectral imaging | |
CN109490223B (en) | Target detection and identification system and method based on programmable hyperspectral imaging | |
Goudail et al. | When is polarimetric imaging preferable to intensity imaging for target detection? | |
CN107203743B (en) | Face depth tracking device and implementation method | |
WO2017020402A1 (en) | Mail detection device and method | |
CN113252168A (en) | Polarization spectrum imaging system based on four-phase modulation | |
CN107356333B (en) | A kind of infrared spectrum polarized imaging system and its parameter optimization configuration method | |
Martino et al. | Material recognition by feature classification using time-of-flight camera | |
US7023546B1 (en) | Real-time imaging spectropolarimeter based on an optical modulator | |
Hoover et al. | Polarization components analysis for invariant discrimination | |
Chen et al. | Camouflage target detection via hyperspectral imaging plus information divergence measurement | |
JP4252252B2 (en) | Imaging device | |
Eriksson et al. | Characterization and performance of a LWIR polarimetric imager | |
WO2020021306A1 (en) | Method for material discrimination and respective implementation system | |
JP6524794B2 (en) | Object identification apparatus, object identification system, object identification method and object identification program | |
Chun et al. | Polarimetric imaging system for automatic target detection and recognition | |
US20190369421A1 (en) | Analyzing microtextured regions of optically anisotropic materials | |
Xu et al. | Integrated zigzag Vander Lugt correlators incorporating an optimal trade-off synthetic discriminant filter for invariant pattern recognition | |
US20230400720A1 (en) | Analyzing microtextured regions of optically anisotropic materials | |
Arı et al. | Scene Text Recognition Using Permutated Autoregressive Sequence Model and YOLOv8 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
FEPP | Fee payment procedure |
Free format text: ENTITY STATUS SET TO SMALL (ORIGINAL EVENT CODE: SMAL); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
AS | Assignment |
Owner name: ADVANCED OPTICAL TECHNOLOGIES, INC., NEW MEXICO Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:HOOVER, BRIAN G.;REYES, PABLO A.;TALIAFERRO, DAVID E.;AND OTHERS;SIGNING DATES FROM 20170427 TO 20170508;REEL/FRAME:045413/0721 Owner name: ADVANCED OPTICAL TECHNOLOGIES, INC., NEW MEXICO Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:HOOVER, BRIAN G.;REYES, PABLO A.;TALIAFERRO, DAVID E.;AND OTHERS;SIGNING DATES FROM 20170427 TO 20170508;REEL/FRAME:045413/0806 |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NON FINAL ACTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE TO NON-FINAL OFFICE ACTION ENTERED AND FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: FINAL REJECTION MAILED |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: RESPONSE AFTER FINAL ACTION FORWARDED TO EXAMINER |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: AWAITING TC RESP., ISSUE FEE NOT PAID |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: NOTICE OF ALLOWANCE MAILED -- APPLICATION RECEIVED IN OFFICE OF PUBLICATIONS |
|
STPP | Information on status: patent application and granting procedure in general |
Free format text: PUBLICATIONS -- ISSUE FEE PAYMENT VERIFIED |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
CC | Certificate of correction | ||
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: SURCHARGE FOR LATE PAYMENT, SMALL ENTITY (ORIGINAL EVENT CODE: M2554); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YR, SMALL ENTITY (ORIGINAL EVENT CODE: M2551); ENTITY STATUS OF PATENT OWNER: SMALL ENTITY Year of fee payment: 4 |