JP2022529624A - サンプルを分類する方法およびシステム - Google Patents
サンプルを分類する方法およびシステム Download PDFInfo
- Publication number
- JP2022529624A JP2022529624A JP2021560551A JP2021560551A JP2022529624A JP 2022529624 A JP2022529624 A JP 2022529624A JP 2021560551 A JP2021560551 A JP 2021560551A JP 2021560551 A JP2021560551 A JP 2021560551A JP 2022529624 A JP2022529624 A JP 2022529624A
- Authority
- JP
- Japan
- Prior art keywords
- measurement
- sample
- data
- vector
- group
- 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.)
- Pending
Links
- 238000005259 measurement Methods 0.000 claims abstract description 242
- 239000013598 vector Substances 0.000 claims abstract description 161
- 238000000034 method Methods 0.000 claims abstract description 69
- 230000003595 spectral effect Effects 0.000 claims abstract description 60
- 238000012545 processing Methods 0.000 claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims abstract description 11
- 239000000523 sample Substances 0.000 claims description 161
- 238000001228 spectrum Methods 0.000 claims description 152
- 230000006870 function Effects 0.000 claims description 95
- 230000005855 radiation Effects 0.000 claims description 36
- 239000013074 reference sample Substances 0.000 claims description 33
- 239000002131 composite material Substances 0.000 claims description 26
- 238000012512 characterization method Methods 0.000 claims description 7
- 238000007405 data analysis Methods 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 5
- 239000010432 diamond Substances 0.000 claims description 4
- 230000005251 gamma ray Effects 0.000 claims description 4
- 239000010437 gem Substances 0.000 claims description 4
- 230000004044 response Effects 0.000 claims description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 3
- 239000011707 mineral Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 claims 2
- 239000000470 constituent Substances 0.000 claims 1
- 238000013479 data entry Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 9
- 238000004891 communication Methods 0.000 description 6
- 238000004876 x-ray fluorescence Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000013459 approach Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 239000000126 substance Substances 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000005461 Bremsstrahlung Effects 0.000 description 1
- 238000002441 X-ray diffraction Methods 0.000 description 1
- 238000002083 X-ray spectrum Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 239000002178 crystalline material Substances 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000026954 response to X-ray Effects 0.000 description 1
- 230000031264 response to gamma radiation Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/389—Precious stones; Pearls
-
- 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
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/44—Raman spectrometry; Scattering spectrometry ; Fluorescence spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/20—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by using diffraction of the radiation by the materials, e.g. for investigating crystal structure; by using scattering of the radiation by the materials, e.g. for investigating non-crystalline materials; by using reflection of the radiation by the materials
- G01N23/207—Diffractometry using detectors, e.g. using a probe in a central position and one or more displaceable detectors in circumferential positions
- G01N23/2076—Diffractometry using detectors, e.g. using a probe in a central position and one or more displaceable detectors in circumferential positions for spectrometry, i.e. using an analysing crystal, e.g. for measuring X-ray fluorescence spectrum of a sample with wavelength-dispersion, i.e. WDXFS
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/22—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material
- G01N23/223—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by measuring secondary emission from the material by irradiating the sample with X-rays or gamma-rays and by measuring X-ray fluorescence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/07—Investigating materials by wave or particle radiation secondary emission
- G01N2223/076—X-ray fluorescence
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/30—Accessories, mechanical or electrical features
- G01N2223/304—Accessories, mechanical or electrical features electric circuits, signal processing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/30—Accessories, mechanical or electrical features
- G01N2223/305—Accessories, mechanical or electrical features computer simulations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/616—Specific applications or type of materials earth materials
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/14—Classification; Matching by matching peak patterns
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medicinal Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Food Science & Technology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Dispersion Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
対象のサンプルをモデルベース分析しサンプル分類を管理するための方法およびシステムが提供される。所定のモデル化データが提供され、これは、スペクトル線の形状を有する所定の関数に基づくそれぞれのK個の測定スキームに対するK個のモデルを示すデータと、異なるサンプルが関連するM個の所定のグループのM個の特性ベクトルを示すデータと、M個のグループに対する重みの共通ベクトルを示すデータとを含む所定のモデル化データとを含む。データプロセッサがこのデータを利用し、前記所定のモデル化データを用いて対象のサンプルの測定スペクトルデータにモデルベース処理を適用し、前記M個の所定のグループの1つに対する前記対象の特定のサンプルの関係を示す分類データを生成するように動作する。
Description
本発明は、サンプルの測定値をモデル化およびモデルに基づいた分析、およびサンプルの分類の分野にある。
本発明者らは、様々な産業、特に、鉱物やダイヤモンドなどの貴石のような対象物の製造および流通を扱う産業において、興味のある対象物/サンプルを、共通または類似の特性を持つ対象物/サンプルの特定のグループに関連するものとして識別/分類する必要があるかもしれないことを見出した。これには、対象物の原産地および/または対象物の原産地の地理的な位置の1つまたは複数の構造パラメータが含まれ得る。
本発明者らはまた、同一のグループ(すなわち、予め定義されたグループ関連特性またはグループ固有の特性を有するグループ)に属する対象物が、その対象物のスペクトルによって、1つまたは複数の他のグループと区別される態様で分類できることを見出した。例えば、そのようなスペクトルデータは、X線またはガンマ線放射に対する対象物/サンプルの蛍光X線(XRF)応答を示すものでありうる。
したがって、本発明の技術は、測定されたスペクトルデータに基づいて対象のサンプルを関連するグループに分類する際に使用される、新規のモデルデータの作成を可能にする新規のモデル化技術を提供する。言い換えれば、本発明は、測定された1つまたは複数のスペクトルに基づいて、対象のサンプルを複数の予め定められたグループの1つに関連付けるための、特定の新規のモデルベースのアプローチを提供する。
したがって、本発明の1つの広範な態様によれば、対象のサンプルをモデルベース分析する方法であって、所定の異なる特性を有する2つ以上の異なるグループに属する参照サンプルに関連する特定の参照スペクトルデータの分析、およびモデル化されたデータの作成を含む方法が提供される。より詳細には、本方法は以下の工程を含む:
所定の異なる特性を有するM個のグループに関連する複数のN個の参照サンプルにおいて実行された数K個の測定スキームのスペクトル測定を示す参照データを提供する工程であって、該参照データは、複数の(NxK)個の測定参照スペクトルを含む未加工測定データを含み、かつ、前記参照サンプルのそれぞれが前記M個のグループのそれぞれ1つに対応することを示すデータを含む、工程;
前記複数の(NxK)個の測定基準スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定する工程であって、前記モデルは、スペクトル線の形状を有する所定の関数に基づいており、かつ、それぞれの測定スキームに関連している、工程;
前記K個のモデルのそれぞれに、それぞれの測定スキームに対応するN個の測定参照スペクトルのそれぞれをフィッティングし、かつ、前記参照サンプルのそれぞれについて、前記数K個の測定スキームについてサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルのそれぞれをそれぞれの成分のベクトルで表現する工程;
前記M個のグループのそれぞれ1つに対する前記サンプルの各々の対応を示すデータを利用し、各グループについて、そのグループに関連するサンプルの成分のベクトルを分析し、かつ、そのグループの特性ベクトルを示すデータを決定する工程;および
前記参照サンプルの成分のベクトルと前記グループの特性ベクトルとの間の距離関数に基づいて、前記参照サンプルの成分のベクトルのすべてをそれぞれのグループに関連付けるための複合尤度(combined likelihood)を最大にする距離関数の重みパラメータを決定し、それによって前記距離関数の前記重みパラメータの共通ベクトルを提供する工程;
それぞれのK個の測定スキームに対するK個のモデルを示すデータ、グループの特性ベクトルを示すデータ、およびM個のグループに対する重みの共通ベクトルを示すデータを含むモデル化データを保存する工程であって、それによって、前記モデル化データを用いて対象のサンプルの未加工測定スペクトルデータをモデルベース分析することにより、対象のサンプルを分類して前記M個のグループの1つに関連付けることを可能にする、工程。
所定の異なる特性を有するM個のグループに関連する複数のN個の参照サンプルにおいて実行された数K個の測定スキームのスペクトル測定を示す参照データを提供する工程であって、該参照データは、複数の(NxK)個の測定参照スペクトルを含む未加工測定データを含み、かつ、前記参照サンプルのそれぞれが前記M個のグループのそれぞれ1つに対応することを示すデータを含む、工程;
前記複数の(NxK)個の測定基準スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定する工程であって、前記モデルは、スペクトル線の形状を有する所定の関数に基づいており、かつ、それぞれの測定スキームに関連している、工程;
前記K個のモデルのそれぞれに、それぞれの測定スキームに対応するN個の測定参照スペクトルのそれぞれをフィッティングし、かつ、前記参照サンプルのそれぞれについて、前記数K個の測定スキームについてサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルのそれぞれをそれぞれの成分のベクトルで表現する工程;
前記M個のグループのそれぞれ1つに対する前記サンプルの各々の対応を示すデータを利用し、各グループについて、そのグループに関連するサンプルの成分のベクトルを分析し、かつ、そのグループの特性ベクトルを示すデータを決定する工程;および
前記参照サンプルの成分のベクトルと前記グループの特性ベクトルとの間の距離関数に基づいて、前記参照サンプルの成分のベクトルのすべてをそれぞれのグループに関連付けるための複合尤度(combined likelihood)を最大にする距離関数の重みパラメータを決定し、それによって前記距離関数の前記重みパラメータの共通ベクトルを提供する工程;
それぞれのK個の測定スキームに対するK個のモデルを示すデータ、グループの特性ベクトルを示すデータ、およびM個のグループに対する重みの共通ベクトルを示すデータを含むモデル化データを保存する工程であって、それによって、前記モデル化データを用いて対象のサンプルの未加工測定スペクトルデータをモデルベース分析することにより、対象のサンプルを分類して前記M個のグループの1つに関連付けることを可能にする、工程。
1つまたは複数の測定スキームを使用して実行された前記サンプルの未加工測定スペクトルデータからの対象のサンプル(いわゆる「未知のサンプル」)の分類は、以下のように行うことができる:
対象のサンプルの未加工測定スペクトルデータに基づいて、K個のデータピースが、K個の測定スキームの下で対象のサンプルのK個の測定スペクトルに対応して決定され、
モデルベース分析は、K個のデータピースに適用され、以下の工程を含む:
保存されたK個のモデルを使用して、対象のサンプルに対する前記K個の測定スペクトルのそれぞれを保存されたK個のモデルのそれぞれ1つにフィッティングし、かつ、前記K個の測定スペクトルのそれぞれに対するベストフィット条件に基づいて、前記K個の測定スキームのすべてのサンプルの複合ベクトル表現を作成する工程;
前記重みの共通ベクトルに前記距離関数を適用して、前記サンプルの複合ベクトル表現の前記グループの特性ベクトルのそれぞれまでの距離を決定し、かつ、決定された距離が最小となるグループに前記サンプルを関連付ける工程。
対象のサンプルの未加工測定スペクトルデータに基づいて、K個のデータピースが、K個の測定スキームの下で対象のサンプルのK個の測定スペクトルに対応して決定され、
モデルベース分析は、K個のデータピースに適用され、以下の工程を含む:
保存されたK個のモデルを使用して、対象のサンプルに対する前記K個の測定スペクトルのそれぞれを保存されたK個のモデルのそれぞれ1つにフィッティングし、かつ、前記K個の測定スペクトルのそれぞれに対するベストフィット条件に基づいて、前記K個の測定スキームのすべてのサンプルの複合ベクトル表現を作成する工程;
前記重みの共通ベクトルに前記距離関数を適用して、前記サンプルの複合ベクトル表現の前記グループの特性ベクトルのそれぞれまでの距離を決定し、かつ、決定された距離が最小となるグループに前記サンプルを関連付ける工程。
概して、本発明の技術は、適切なモデル化データを作成するとともに、未知のサンプルを適切に分類する。しかしながら、好ましくは、2つ以上の異なる測定スキームが使用される。測定スキーム/条件は、1つまたは複数のパラメータで互いに異なっていてもよい。いくつかの実施形態では、そのようなパラメータは、以下の1つまたは複数を含むことができる:一次放射線強度、一次放射線中の光子のエネルギー分布(これは、一次放射線を放出する管の電流および電圧によって設定されうる)、および/または放射線放出源におけるフィルタ。追加的に、または代替的に、異なる測定スキームでは、以下の1つまたは複数のバリエーションを使用することができ:一次放射線信号のコリメーション、照射スポットのサイズ、検出器での放射線応答信号のフィルタリング、放射線源の幾何学的構成、サンプルの表面および放射線源および/または放射線検出器の相対的な向きおよび配置(例えば、サンプルの表面、放射線源、および検出器の間の角度および距離)、これらは、スペクトルの測定に影響を与える可能性があり、異なる測定条件を作成するために変化させることができる。さらに、測定される一部またはすべてのスペクトルについて、サンプルを1つまたは複数の軸の周りに回転させ、様々なサンプルの方向付けの間にサンプルからの放出された放射線部分の計数を単一のスペクトルに収集してもよい。
測定スキームごとに作成されたモデルは、スペクトル線の形状の所定の関数、および、ピースワイズ線形またはピースワイズ多項式関数である特定のピースワイズ(またはハイブリッド)関数に基づいて、混合モデルとして構成されている。このようなスペクトル線の形状を持つ関数としては、ローレンツ関数、ガウス関数および/またはフォークト関数が挙げられる。
グループ特性ベクトルには、同一グループの参照サンプルを表す成分ベクトル内の成分の平均値が含まれている。この平均値および標準偏差に距離関数を対応させることで、成分ベクトル内の成分値の広がりの量が表現される。
K個のモデルを決定するための、複数の(NxK)の測定された参照スペクトルの処理は、以下のように行われてもよい:
i番目の測定スキーム(i=1、…、K)に対応するN個の参照サンプルの複数の測定された参照スペクトルの各i番目に対して、平均測定参照スペクトルが決定され;かつ
スペクトル線の形状を有する所定の関数に従って各i番目の平均測定基準スペクトルに対して所定の変換を適用し、i番目の測定スキームに対応するそれぞれのi番目のモデルを取得し、それによりK個の測定スキームについてK個のモデルを取得する。
i番目の測定スキーム(i=1、…、K)に対応するN個の参照サンプルの複数の測定された参照スペクトルの各i番目に対して、平均測定参照スペクトルが決定され;かつ
スペクトル線の形状を有する所定の関数に従って各i番目の平均測定基準スペクトルに対して所定の変換を適用し、i番目の測定スキームに対応するそれぞれのi番目のモデルを取得し、それによりK個の測定スキームについてK個のモデルを取得する。
本発明の別の広範な態様によれば、サンプル上の測定をモデル化するためのデータ分析システムが提供される。このシステムは、測定システム(例えば、X線またはガンマ線放射に対するサンプルの蛍光X線(XRF)応答を測定するためのもの)、および、測定された参照データに基づいて、対象のサンプルのさらなる分類を可能にするモデル化されたデータを決定するように構成および動作可能なコントロールを含む。より詳細には、測定システムは、数Kの測定スキームの下で、所定の異なる特性のM個のグループに関連する複数のN個の参照サンプルに対してスペクトル測定を行い、前記M個のグループに関連する複数の(NxK)個の測定参照スペクトルを含む測定参照データを生成するように構成および動作可能である。コントロールシステムは、以下を含む:
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュールであって、前記モデルは、スペクトル線の形状を有する所定の関数に基づいており、それぞれの測定スキームに関連する;
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個のモデルのそれぞれについて、それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々をモデルにフィッティングし;参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現する;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの共通ベクトルを提供するように構成および動作可能な重み付けモジュール;および
それぞれのK個の測定スキームに対するK個のモデルを示すデータであって、グループの特性ベクトルを示すデータ、および、M個のグループに対する重みの共通ベクトルを示すデータを含む、保存されるモデル化データを生成するように構成および動作可能な出力ユーティリティ。
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュールであって、前記モデルは、スペクトル線の形状を有する所定の関数に基づいており、それぞれの測定スキームに関連する;
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個のモデルのそれぞれについて、それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々をモデルにフィッティングし;参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現する;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの共通ベクトルを提供するように構成および動作可能な重み付けモジュール;および
それぞれのK個の測定スキームに対するK個のモデルを示すデータであって、グループの特性ベクトルを示すデータ、および、M個のグループに対する重みの共通ベクトルを示すデータを含む、保存されるモデル化データを生成するように構成および動作可能な出力ユーティリティ。
本発明は、さらに広範な態様において、以下を含むサンプル分類システムを提供する:
数Kの測定スキームの下でサンプルのスペクトル測定を行い、かつ、測定されたサンプルのそれぞれについて、K個の測定スキームにそれぞれ対応する測定スペクトルを示すK個の測定データピースを含む測定スペクトルデータを生成するように構成および動作可能な測定システム;
前記測定システムと通信し対象の測定スペクトルデータを受信するように構成および動作可能であり、かつ、スペクトル線の形状を有する所定の関数に基づくそれぞれのK個の測定スキームに対するK個のモデルを示すデータと、異なるサンプルが関連するM個の所定のグループのM個の特性ベクトルを示すデータと、M個のグループに対する重みの共通ベクトルを示すデータとを含む所定のモデル化データを格納するメモリと通信するように構成および動作可能なコントロールシステムであって、前記所定のモデル化データを用いて対象のサンプルの受信した測定スペクトルデータにモデルベース処理を適用し、前記M個の所定のグループの1つに対する前記対象の特定のサンプルの関係を示す分類データを生成するように構成および動作可能なデータプロセッサを備える、コントロールシステム。
数Kの測定スキームの下でサンプルのスペクトル測定を行い、かつ、測定されたサンプルのそれぞれについて、K個の測定スキームにそれぞれ対応する測定スペクトルを示すK個の測定データピースを含む測定スペクトルデータを生成するように構成および動作可能な測定システム;
前記測定システムと通信し対象の測定スペクトルデータを受信するように構成および動作可能であり、かつ、スペクトル線の形状を有する所定の関数に基づくそれぞれのK個の測定スキームに対するK個のモデルを示すデータと、異なるサンプルが関連するM個の所定のグループのM個の特性ベクトルを示すデータと、M個のグループに対する重みの共通ベクトルを示すデータとを含む所定のモデル化データを格納するメモリと通信するように構成および動作可能なコントロールシステムであって、前記所定のモデル化データを用いて対象のサンプルの受信した測定スペクトルデータにモデルベース処理を適用し、前記M個の所定のグループの1つに対する前記対象の特定のサンプルの関係を示す分類データを生成するように構成および動作可能なデータプロセッサを備える、コントロールシステム。
いくつかの実施形態では、コントロールシステムは以下を含む:
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個の測定スペクトルのそれぞれについて、測定スペクトルをそれぞれのモデルにフィッティングし、K個のベストフィット条件スペクトルを得る;および、前記K個のベストフィット条件スペクトルを用いて、前記K個の測定スキームすべてについて関心対象サンプルの複合ベクトル表現を作成する;
前記共通の重みベクトルを用いて所定の距離関数を利用し、前記M個のグループの前記M個の特性ベクトルのそれぞれまでの対象のサンプルの前記複合ベクトル表現の距離を決定し、決定された距離が最小となるグループに前記対象のサンプルを関連付けるように構成および動作可能な分類モジュール。
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個の測定スペクトルのそれぞれについて、測定スペクトルをそれぞれのモデルにフィッティングし、K個のベストフィット条件スペクトルを得る;および、前記K個のベストフィット条件スペクトルを用いて、前記K個の測定スキームすべてについて関心対象サンプルの複合ベクトル表現を作成する;
前記共通の重みベクトルを用いて所定の距離関数を利用し、前記M個のグループの前記M個の特性ベクトルのそれぞれまでの対象のサンプルの前記複合ベクトル表現の距離を決定し、決定された距離が最小となるグループに前記対象のサンプルを関連付けるように構成および動作可能な分類モジュール。
いくつかの実施形態では、コントロールシステムは、前記M個のグループに関連する複数のN個の参照サンプル上で実行された前記測定スキームの数Kのスペクトル参照測定に対応する測定スペクトルデータに基づいて、所定のモデル化データを決定するようにさらに構成および動作可能であり、ここで、スペクトル参照データは複数の(NxK)個の測定参照スペクトルを含み、前記M個のグループのそれぞれの1つに対する参照サンプルのそれぞれの対応関係を示すデータを含む。コントロールシステムは、以下を含む:
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュール;
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個のモデルのそれぞれについて、それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々をモデルにフィッティングし;参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現する;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;および
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする前記所定の距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの前記共通ベクトルを提供するように構成および動作可能な重み付けモジュール。
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュール;
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個のモデルのそれぞれについて、それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々をモデルにフィッティングし;参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現する;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;および
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする前記所定の距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの前記共通ベクトルを提供するように構成および動作可能な重み付けモジュール。
本発明のさらに広範な態様によれば、サンプル分類の管理に使用するコントロールシステムが提供される。コントロールシステムは、測定データプロバイダと通信し対象のサンプルの測定スペクトルデータを受信するように構成および動作可能であり、かつ、スペクトル線の形状を有する所定の関数に基づくそれぞれのK個の測定スキームに対するK個のモデルを示すデータと、異なるサンプルが関連するM個の所定のグループのM個の特性ベクトルを示すデータと、M個のグループに対する重みの共通ベクトルを示すデータとを含む所定のモデル化データを格納するメモリと通信するように構成および動作可能である。コントロールシステムは、前記所定のモデル化データを用いて対象のサンプルの受信した測定スペクトルデータにモデルベース処理を適用し、前記M個の所定のグループの1つに対する前記対象の特定のサンプルの関係を示す分類データを生成するように構成および動作可能なデータプロセッサを備える。
本明細書に開示されている主題をよりよく理解し、それが実際にどのように実施され得るかを例示するために、添付の図面を参照しながら非限定的な例としてのみ実施形態が説明される。
本発明は、サンプルの測定スペクトルに基づいて、類似/関連サンプルの特性グループに関する/関連するものとしてサンプルを分類するための新しいアプローチを提供する。上述したように、鉱物や貴石(特にダイヤモンド)のようないくつかのタイプのサンプル/対象物は、特定のグループとの関連性/関係性によって識別される必要がある。このグループは、サンプルの原産地および/またはサンプルの原産地の地理的位置の1つまたは複数の構造パラメータによって記述されうる。同じグループ(すなわち、予め定められたグループ関連のおよびグループ固有の特性を有するグループ)に関連するサンプルは、1つまたは複数の他のグループのサンプル/スペクトルと区別する態様で、そのスペクトルデータに基づいて分類することができる。
本発明は、サンプルの未加工測定スペクトルデータに基づいて、対象のサンプルを関連するグループに分類するために使用される、新規のモデル化データを作成するための新規の技術を提供するものである。
図1を参照すると、ブロック図によって、サンプルを分類するためにさらに使用されるモデル化データを作成するための本発明のデータ分析システム10が示されている。システム10は、測定データプロバイダ12とのデータ通信のために構成されたコントロールシステムである。コントロールシステム10は、典型的にはコンピュータシステムであり、測定データプロバイダの一部/またはそれと一体であってもよく、あるいは、例えばクラウドコンピューティング技術を用いて、任意の既知の適切な通信技術およびデータプロトコルを用いて、通信ネットワークを介して測定データプロバイダと通信してもよい。リモートエンティティ間のデータ通信ネットワークおよびプロトコルの構築および運用は、それ自体がよく知られており、本発明の一部を構成するものではないため、詳細に説明する必要はない。
測定データプロバイダは、この非限定的な実施例に示すように、測定システム14自体によって構成されてもよいし、任意の既知の適切な通信技術を用いて、測定システムとデータ通信する別個の記憶装置であってもよい。この特定の実施例に示すように、測定システム14は、放射線源14A、放射線検出器14B、コントローラ14C、並びにサンプル支持ユニット14Dを含む。
具体的に示されていないが、測定システムは、測定手順を管理するための様々な他のユニットおよびハードウェア/ソフトウェア・ユーティリティを含みうるが、これらは本発明の一部を構成するものではないため、以下の点を除いて詳細に説明する必要がないことを理解されたい:本発明の目的のために、モデル化データの作成に必要な測定データは、各サンプルについて、異なる測定条件/スキームで得られた所定の数K(K≧1)の測定スペクトルを含む。概して、モデル化データの作成には、単一の測定スキーム(K=1)を用いた測定が十分でありうる。しかしながら、スペクトル測定、さらに様々な形状および幾何形状の容積測定サンプルを扱う場合には、異なる測定スキームに対応する複数のスペクトルを提供することが好ましい。
いくつかの実施形態では、表面上および/または容積内に様々なマーキングを有しうる貴石、特にダイヤモンドの測定に適しているスペクトルデータは、X線またはガンマ線放射に対するサンプルの蛍光X線(XRF)反応を示しうる。したがって、放射線源14Aは、一次励起放射線によってサンプルを照射して、サンプルから二次蛍光X線(XRF)反応の放出を誘発するように構成されたX線またはガンマ線の放射線源であってもよく、放射線検出器14Bは、蛍光X線(XRF)を検出し、検出された放射線を示す測定スペクトルデータを生成するように構成されている。このような測定システムは、例えば、国際特許出願公開第2016/157185号、同第2017/175219号、同第2018/051353号に記載されており、これらは全て本出願の譲受人に譲渡されており、参照により本明細書に組み込まれる。
異なる測定スキームを設定するパラメータ/条件には、以下の1つまたは複数が含まれうる:一次放射線のパラメータ(例えば、強度、コリメーション、スポットサイズ、一次放射線中の光子のエネルギー分布など);検出される二次放射線のフィルタリングパラメータ/条件;並びに、例えばサンプルの支持ユニット14Dを1つまたは複数の軸の周りで回転させることによって達成される(様々なサンプルの向きからの計数が単一のスペクトルに収集されるように)、放射線源および/または検出器に対するサンプルの向き。したがって、支持ユニット14Dは、放射線源および/または検出器に対する測定平面の位置を調整するだけでなく、測定平面内の自身の位置を調整するための1つまたは複数のドライバと関連していてもよいことを理解すべきである。また、放射線源14Aは、放射線源の動作パラメータ(例えば、一次放射線を放出する管の電流および/または電圧;および/またはフィルタ)を調整/変化させるための1つまたは複数のドライバと関連していてもよい;並びに、検出器14Bは、放射線検出器の入力におけるフィルタを動作/変化させるためのフィルタリングアセンブリと関連していてもよい。さらに、システムの性能を改善/最適化するために、放射線源および検出器の幾何学的特性は可変/調整可能でもよい。このような幾何学的特性は、以下の1つまたは複数を含みうる:X線源からサンプルの所定の表面領域までの距離;この表面領域から検出器までの距離(検出面);照射チャネルの角度方向(X線源から伝搬する一次X線ビーム(一次ビーム伝搬軸)とサンプルの表面との間の角度);および収集/検出チャネルの角度方向(サンプルから検出器に向かってくる二次X線放射(二次ビーム軸)とサンプルの表面との間の角度)。
したがって、システムコントローラ14Cは、K個の測定スキームのそれぞれを画定し、したがって各サンプル上で測定セッションを動作するために、測定システムの要素の上記に例示されたパラメータ/条件のいずれかを変化/調整するように構成され、動作可能である。
モデル化データの作成時には、いわゆる「参照サンプル」に対してスペクトル測定が行われるため、図中では測定データを「参照データ」と称する。参照サンプルとは、特定のグループとの関連性が知られているサンプルである。
したがって、測定システム12は、N個の参照サンプルにスペクトル測定を適用するように動作し、各サンプルはK個の異なる測定スキームで測定される。これらのN個の参照サンプルは、M個のグループに関連するサンプルを含み、各g番目のグループ(g=1、...、M)は、所定の異なる(グループ固有の/関連する)特性を有する。したがって、概して、第1グループG1はn1サンプルを含み、第2グループG2はn2サンプルを含み、…、M番目のグループGMはnMサンプルを含み、ここでは以下である。
各データピースは、参照サンプルRSの分光感度を示す/対応するものである。したがって、図2Aにも示されているように、参照測定データには以下が含まれる:
グループG1について:
グループG2について:
...
グループGmについて:
グループG1について:
グループGmについて:
ここでは、指標は以下のようであることを理解されたい:(RS(i)
n)g、ここで、上付きのインデックスiはi番目の測定スキーム(i=1、…、K)に対応し,下付きのインデックスnおよびgはg番目のグループのn番目のサンプルに対応する。したがって、例えば(RS(2)
(3))4は、測定スキーム2に従って測定されたグループ4のサンプル3の参照スペクトルを称する。
分類される対象の未知のサンプルに対して同様の測定が行われるが、その場合、グループとのサンプルの関連性は未知であり決定されるものであることが理解されるべきであり、以下にさらに説明される。したがって、このような未知のサンプルの場合、測定されたスペクトルデータは、モデル化データの作成に使用された異なる測定スキームに対応するKスペクトルを含む。
上述したように、コントロールシステム10は、コンピュータシステムとして構成されており、これは、データ入力および出力ユーティリティ16、18;メモリ20;およびデータプロセッサ22のような主要な構造および機能パーツ/ユーティリティを含む。データプロセッサは、モデル作成モジュール22A、フィッティングモジュール22B、グループ特性化モジュール22Cおよび重み付けモジュール22Dを含む。受信される参照スペクトルデータは、典型的にはメモリ20に格納され、その後、プロセッサ22によって使用されてモデル化データを作成する。
モデル作成モジュール22Aは、(NxK)個の測定された基準スペクトルを処理して、K個の測定スキームのそれぞれについてモデルを決定する、すなわち、サンプルのスペクトル応答を記述するK個のモデルを決定するように予めプログラムされている。このモデルは、スペクトル線の形状を有し、それぞれの測定スキームに関連する所定の関数に基づいている。そのようなスペクトル線の形状の所定の関数は、例えば、ローレンツ関数、ガウス関数またはフォークト関数を含むことができ、そのパラメータは、線の位置、最大の高さおよび幅(または半値幅)を含む。以下にさらに説明されるように、モデルは、このようなスペクトル線の形状の所定の関数、および、特定のピースワイズ線形関数を含んでいてもよい。モデル作成処理については、図3を参照してさらに以下でより具体的に説明する。
フィッティングモジュール22Bは、反復的なフィッティング手順において、測定された参照スペクトルのそれぞれを、対応する測定スキームのモデルと比較するように構成されている。フィッティングの際には、ベストフィット条件を介してモデルパラメータが最適化され、各参照スペクトルに対してそのベクトル表現が決定される。言い換えれば、参照サンプルのそれぞれは、それぞれの成分ベクトルによって表される。このようなベクトル成分表現は、すべてのK個の測定スキームに対して結合されたものであり、これがサンプルの表現であることが理解されるべきであり、以下でさらに具体的に説明される。
グループ特性評価モジュール22Cは、グループの特性ベクトルを決定するように動作する。この目的のために、モジュールは、M個のグループのそれぞれの1つに対する参照サンプルのそれぞれの対応関係を示すデータに基づいて、サンプルの成分ベクトルを分析する。
重み付けモジュール22Dは、参照サンプルのすべての成分ベクトルをそれぞれのグループに関連付けるための複合尤度の最大値に対応するベクトル成分の重みパラメータを決定するように構成される。これにより、重みの共通ベクトルが決定される(すべてのグループに共通)。
そのようにして決定されたデータは、以下を含むモデル化データを形成する:(i)K個の測定スキームのそれぞれに対するK個のモデルを示すデータ、(ii)グループの特性ベクトルを示すデータ、および(iii)すべてのグループに対して重みの共通ベクトルを示すデータ。
次に図3を参照し、これは、上述のようにして得ることができかつM個のグループに関連するN個の参照サンプルについてK個の測定スキームで得られた参照スペクトルを含む、測定された参照データからモデル化データを生成/作成するための本発明の方法のフロー図100を例示する。概して、同じサンプルから、同じ測定スキームを用いて複数のスペクトルが得られることに留意すべきである。
したがって、参照測定データが提供され(工程102)、測定システムまたは別個の記憶装置(すなわち、測定データプロバイダ)のいずれかでアクセスすることができる。必要に応じて、測定スペクトルの前処理を行ってもよい。これは、各スペクトルにおいて、モデル化および/または分類が行われる関心領域を決定すること、および/または、スペクトルからバックグラウンドノイズおよび/またはアーチファクト信号を識別および除去することを目的とすることができる。選択された関心領域は、概して、スペクトルが測定される測定条件によって影響を受ける可能性がある。ノイズおよびアーチファクト信号には、例えば、結晶性の材料で作製されたサンプルの場合、サンプルの結晶構造に起因するX線回折ピークが含まれうる。さらに、XRFスペクトルの場合、これらのアーチファクト信号には、放射線源、検出器、またはサンプルの周辺に見られる物質(サンプル自体にはない)に由来するピーク、並びに、パイルアップピークおよびバックグラウンド計数または他のプロセスに由来する信号が含まれうる。ノイズ、および/またはアーチファクト信号を除去するためにスペクトルを処理する目的で、任意の既知の適切な技術を使用することができ、例えば、本出願の譲受人に譲渡され参照により本明細書に組み込まれている上述の国際公開第2016/157185号に記載されている方法を使用することができる。
したがって、モデル化データ作成のために処理される参照測定スペクトルは、前処理されたスペクトルであってもよく、また、サンプルに関連するスペクトルまたはサンプル内の以前に画定された関心領域のスペクトルであってもよい。このような前処理されたまたはされていない参照スペクトルデータが処理および分析されて、参照スペクトルの取得に使用されるK個の測定スキームに対応するK個のモデルを作成する(工程104)。この目的のために、各測定スキームについて、平均化されたスペクトルが得られる、すなわち、同じ測定スキームに対応する参照スペクトルについて、これらのスペクトルをすべて合計しサンプル数で割ることによって平均化が行われる。より詳細には、各i番目の測定スキーム(i=1、...、K)について以下である:
ここで、
は以下の測定スキームに対応する合計スペクトルである:
したがって、K個のそのような平均化されたスペクトルが決定される。各グループの平均化されたスペクトルはさらに、スペクトル線の形状(例えばガウス)を有する所定の基底関数BFおよびバックグラウンド関数AF(例えばピースワイズ線形関数またはピースワイズ多項式関数)に従って変換Tを平均化スペクトルに適用することによって、対応するモデル(いわゆる「混合モデル」)を作成するよう処理される。より詳細には、各i番目の測定スキームに対して以下である:
例えば、そのような変換の結果は、以下のようであり:
ここで、AF=B(x)はバックグラウンド関数であり、BF=P(x)は基底関数であり、これは通常は主関数の領域の間隔xにおいて異なるピークを持つ複数のサブ関数(例えばガウス)の形をしており、インデックスjは基底関数のj番目のサブ関数(特定のガウス/ピークを持つ)に対応している。
本発明の目的のために、スペクトル測定データを考慮する場合、モデルはピーク関数およびバックグラウンド関数を持つように選択される。ピーク関数は、対応する平均化スペクトルのピークを表し、これは一般的にサンプル内の物質および元素に関連するが、さらに、サンプル、サンプルの近傍(例えば、サンプルカップ内)、放射線源または検出器内の様々な他の現象およびプロセスに関連しうる。例えば、放射線源にプリセットされた異物に対応しうるアーチファクトピークである。
特定の非限定的な実施例では、測定されたスペクトルはX線スペクトルであり、アーチファクトピークには、コンプトンピーク、レイリーピーク、パイルアップピーク、制動放射、並びに他のプロセスに由来するピークが含まれうる。バックグラウンド関数は、対応する平均化スペクトルのバックグラウンドを表す。
したがって、特定のi番目の測定条件/スキームの下で測定された平均化スペクトルに対応するスペクトルモデルは、以下のような形式であってもよい:
ここで、B(x)はバックグラウンド関数であり、(入射光子の)エネルギーxについてカウント数またはカウント毎秒(CPS)に対するバックグラウンドの寄与を表す;Pj(x)はピーク関数であり、光子エネルギーxのカウント数またはCPSに対するピークの寄与を表す。
別の実施例では、ピーク関数はローレンツ関数である。ある実施例では、バックグラウンド関数B(x)は、区分けされた多項式関数によって定義されるスプラインである。ある実施例では、バックグラウンド関数は、指数多項式である。
次いで、このようにして決定されたK個のモデルを用いて、各参照スペクトルについて、対応する成分ベクトルを決定する(工程106)。これは、i番目の測定スキームに対応するg番目のグループのn番目のサンプルの各参照スペクトル(R(i)
n)g(i=1、…、K)を、最良のフィッティング条件が得られるまで、選択されたモデルパラメータ(例えば、hj、参照スペクトルのピークに主に対応するピークの高さ)の値を変化させながら、それぞれのi番目のモデルにフィッティングすることによって行われる。これにより、特定の測定スキームの特定のサンプルの参照スペクトルに対応するパラメータセットが得られる。次いで、特定のサンプルに対応するK個のパラメータセットをすべて組み合わせて、参照サンプルごとに単一のパラメータベクトルを作成する。これは、サンプルに適用されるすべての測定スキームについて、参照サンプルに関連する/参照サンプルを表すパラメータの「複合」ベクトルであることが理解されたい。
より詳細には、モデルスペクトルのピークのパラメータを測定スペクトルに調整することによってフィッティングを行う。そのために、ピーク関数のパラメータの1つまたは複数を選択し、測定された参照スペクトルとモデルとの間の一致が得られるように設定する。これは、(所定の測定条件の)モデルと、選択されたピーク関数のパラメータによって決定されかつこれらのパラメータの不確実性にも依存しうる測定されたスペクトルとの間の距離の尺度を最小化するように、選択されたパラメータを設定することによって行うことができる。
選択されたパラメータがピーク関数の高さである実施例では、モデルと測定されたスペクトル(どちらも同じ測定条件に対応する)の間の距離は次のように定義される:
ここで、yrはエネルギーrにおけるスペクトルの測定値である;Trは同じエネルギーにおけるモデル(変換関数)の対応する値である;そして、Δyrは測定値の不確実性(測定の種類による)である;Tr(モデル)の値はベストフィット条件によって最適化される。カウントまたはカウント毎秒で測定されたピーク高さについて、不確実性は√yjである。
ある実施例において、フィッティングは、例えば非線形最小化によって反復的に行われる。設定されるピーク関数Pj(モデルTに含まれる)の1つまたは複数のパラメータは、特定のサンプルから特定の測定スキームの下で取られたスペクトルに対応するパラメータのベクトルにおける成分jとして定義される。サンプルsに対応する成分ベクトルは、サンプsに対応しK個の異なる測定条件およびバックグラウンドを特徴づけるパラメータ下で取られたすべてのスペクトルからのすべてのパラメータ/成分を、単一の複合成分ベクトルに結合することによって得られる。
ある実施例において、モデルのピークを表すピーク関数はガウス関数であり、サンプルのスペクトルをモデルにフィッティングさせるために設定されるパラメータはガウスの高さhjである。したがって、n番目のサンプルに対応する成分のベクトルは、次のような形になる:
ここで、パラメータ/コンポーネントセットhp、hf、hqのそれぞれは、異なる測定条件で測定されたスペクトルに対応していてもよく、blはバックグラウンドパラメータである。
したがって、N個の測定された参照サンプルを表すN個の成分ベクトル
が得られる(工程106)。これは図2Bにも示されており、これは、参照サンプルをグループに関連付けて、サンプルのスペクトルデータを、サンプルを表す成分のベクトルに変換したことを示す。
このようにして得られたサンプルに関連する成分ベクトル、および、参照サンプルのグループへの関連性に関する既知のデータを用いて、各グループに対する特性ベクトルCVを決定する、すなわち、M個のグループに対するM個のそのような特性ベクトルCV1、CV2、…、CVMを決定する(工程108)。この目的のために、成分のベクトルを処理して、各サンプルがグループ(サンプルのクラスタ)に属する可能性を推定するための式を得る。この推定は以下のように行うことができる。
基準となる分類されたサンプル(g番目のグループに属する)に対応する成分のベクトルの各成分jについて、グループ平均(vj、g)およびグループ標準偏差(σj、g)を評価する。平均値および標準偏差は、成分ベクトル中の各成分の値の広がりを表す距離関数を画定する。上述したように、グループの特性ベクトルは、同一グループの参照サンプルを表す成分ベクトル中の成分の平均値を含んでいる。距離関数は、平均値および標準偏差に関連付けられており、次いでこれらの値を用いて、分類された各サンプルsが各グループに属する尤度Ls(g)について第1の値を計算する。これは、成分ごとに行うことができ、尤度は、成分ベクトルの各成分(サンプルsに関する)がg番目のグループに属する確率、Ps(j、g、wi)の積として定義される:
確率Ps(j、g、wi)は、平均値vj、gおよび標準偏差σj、gに依存し、非負の重みwjにも依存するが、初期値は1に設定されている。
次に、すべてのグループに共通する重み付けベクトルを決定する(工程110)。この目的のために、距離関数の重みパラメータwjは、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする条件に対応して決定される。これは、参照サンプルの成分のベクトルとグループの特性ベクトルとの間の距離関数に基づいて決定される。
より具体的には、重みwjの最適(最終)値は、分類されたサンプルがグループに正しく分類される確率Pcorrを最適化することで得られる。正しく分類される確率は、サンプルがグループに属する確率の、グループ内のすべてのサンプルに対する積のすべてのグループに対する積で表すことができる:
ここで、サンプルsがグループgに属する確率は、正規化された尤度として定義される:
言い換えれば、Pcorrの値を最大化するように重みの値を設定する。最適化処理は、任意の非線形最適化手法(例えば、Levenberg-Marquardt、BFGS、GRG、進化的手法)によって行うことができる。
上述のように、重みのベクトルwjは、グループのM個の特性ベクトルCV1、CV2、...、CVMと、K個の測定スキームに対応するK個のモデルとともに、未知/未分類の関心対象サンプルの分類に使用するモデル化データとして格納される。
これに関連して、次に図4が参照され、これは、未分類のサンプルを分類されたサンプルのグループに関連付けるための、本発明の例示的な方法のフロー図200を示す。
このために、K個の測定スキームに対応して、対象のサンプルの未加工測定スペクトルデータが提供される(工程202)。このような測定データは、上述のように、測定システム14を用いて得られてもよい。測定データは、測定システムから直接、または別の記憶装置から(概して、測定データプロバイダ12から)提供されてもよい。測定データは、K個の測定スキームの下で対象のサンプルのK個の測定スペクトル:MS(1)、MS(2)、...MS(K)にそれぞれ対応するK個のデータピースを含む。
測定データは、上述のモデル化データを用いて、モデルベースの解析/処理が行われる。より詳細には、K個の測定スペクトルから各i番目の測定スペクトルMS(i)を、記憶されたK個のモデルのうちそれぞれi番目のモデルに、ベストフィット条件が得られるまでフィッティングし、K個の測定スペクトルに対するこれらのベストフィット条件のパラメータを用いて、K個のすべての測定スキームに対するサンプルの複合ベクトル表現CVRを作成する-工程204。次に、この複合ベクトル表現CVRを、グループの特性ベクトルCV1、CV2、…CVMにフィッティングして、グループ関連の最大尤度を決定する-工程206。より詳細には、サンプルの複合ベクトル表現CVRについて、各グループに属する尤度Ls(g)が決定され(重みに対する最終値を用いて)、その尤度が最大となるグループがサンプルの関連/相関グループとして選択される(工程208)。この目的のために、重みの共通ベクトルを持つ上述の距離関数を使用して、サンプルの複合ベクトル表現の各グループの特性ベクトルへの距離を決定し、決定された距離が最小となるグループにサンプルを関連付ける。
モデル(モデルスペクトル)を使用することで、次元を下げることができることを理解されたい。実際、未加工データ(測定されたスペクトル)には、約2000個のスペクトルチャネルのカウントまたはカウント毎秒が含まれ、それぞれが(入射光子の)エネルギーバンドに対応している。このモデルでは、特定のピークに属するこれらのチャンネルをすべてグループ化することで、ピークの数を大幅に減らすことができる(それぞれのピークは、例えばガウス関数として上述される)。パラメータの数を大幅に減らすことで、計算能力、時間などのリソースを削減することができる。さらに、モデルベースのアプローチによって、ノイズを低減することができる。チャネル内のカウント数hのノイズは√hであるため、複数のチャンネルのカウントを取得すれば、信号対雑音比を向上させることができる。
本発明はまた、サンプル間の対応関係や相互関係に関する事前の知識を持たずに、サンプルをクラスタリングする、すなわちサンプルをグループやクラスタに分類するための新規技術を提供する。この技術では、「既知の」参照サンプルのグループ/クラスタへの関連付けを使用して準備されたモデル化データは存在しない。サンプルは、サンプルから放出される電磁信号の1つまたは複数のスペクトルを調べることによって分類される。これは、例えば、X線やガンマ線放射に対するサンプルの蛍光X線反応でもよい。
これに関連して、図5を参照し、これは、サンプルをクラスタリングするための本発明の方法のフロー図300を示す。サンプルの測定データは、各サンプルからの1つまたは複数のスペクトルを含んで提供され、ここで、上述したモデリングおよび分類技術と同様に、サンプルごとの測定データは、K個の異なる測定条件/スキームの下で測定されたK個のスペクトルを含む(工程302)。したがって、N個のサンプルに関する測定データは、(NxK)個のスペクトルを含む:
MS1 (1)、MS1 (K)、...、MSN (1)、MSN (K)。
MS1 (1)、MS1 (K)、...、MSN (1)、MSN (K)。
必要に応じて、上述の技術と同様に、測定スペクトルを処理して、各スペクトルに、クラスタリングが処理される関心領域を画定し、スペクトルからバックグラウンドノイズおよび/またはアーチファクト信号を識別および除去する。
測定データは、上述の手法と同様に、平均化スペクトルを決定するために処理される(工程304)。この目的のために、1つまたは複数の和スペクトルが決定され、これらはそれぞれ、同じ測定スキームで測定されたすべてのスペクトルのカウント数(検出器での光子のカウント数)またはカウント毎秒(CPS)の合計、対、サンプルから到達する入射光子の測定周波数(エネルギー)に対応する。
平均化スペクトルを使用して、図1~3を参照して上述した方法で、K個の測定スキームに対応するモデルを作成する(工程306)。各測定スペクトルは、対応するモデル(同じ測定スキームのスペクトルモデルである)にフィッティングされ、各サンプルについて、成分ベクトルが(上述した技術と同様に)決定される-工程308。
これらの成分のベクトルを使用して、サンプルをグループに反復的に分類する-工程310。分類は、クラスタリングアルゴリズムを用いて実施してもよい。ある実施例では、クラスタリングは、セントロイドベースのクラスタリングアルゴリズムによって実施されてもよい。より詳細には、セットされたサンプルは、グループに分割され、グループの数Mは、サンプルに関する何らかの事前の知識(例えば、サンプルは既知の数のソースに由来しうる)に基づいて、またはランダムに決定される。サンプルのグループへの割り当ては、ランダムに行われてもよい。サンプルの各クラスタのセントロイドは、クラスタ内のサンプルに関連する成分のベクトルの各成分の平均値を評価することによって決定される。その平均値のベクトルをクラスタのセントロイドと画定する。
特定の実施例では、K-means型アルゴリズムによってクラスタリングが行われ、クラスタリングは反復的に進行する。各反復において、パラメータの各ベクトルの各セントロイドへの距離が評価される。グループvjgのセントロイドからのベクトルvjの距離は、ユークリッド距離または正規化ユークリッド距離として画定され、例えば、各成分の距離は、その成分のグループ標準偏差で正規化される:
その後、成分のベクトルは、そのクラスタ(すなわちセントロイド)までの距離が最短であれば、異なるクラスタに再割り当てされてもよい。ベクトル間の距離は、ユークリッド距離として画定されてもよい。さらに、階層型クラスタリング、密度ベースのクラスタリングなどのような、他のクラスタリング方法を使用してもよい。
したがって、本発明は、サンプルの測定スペクトルデータをモデルベースで分析して、サンプルを関連/類似サンプルのグループに分類/関連付けるための新規技術、並びに、モデル化データ作成のための新規技術を提供する。本発明の技術は、サンプル/オブジェクトのクラスタリング/グループ化を扱う様々なアプリケーションに使用することができる。データ解析システムは、スペクトル測定システムと一体化していても、または別のコントロールシステムにあってもよく、また、データ解析プロセスは、いわゆる「オンライン」またはオフラインモードで実行されてもよい。
Claims (20)
- 対象のサンプルをモデルベース分析する方法であって、
所定の異なる特性を有するM個のグループに関連する複数のN個の参照サンプルにおいて実行された数K個の測定スキームのスペクトル測定を示す参照データを提供する工程であって、該参照データは、複数の(NxK)個の測定参照スペクトルを含む未加工測定データを含み、かつ、前記参照サンプルのそれぞれが前記M個のグループのそれぞれ1つに対応することを示すデータを含む、工程;
前記複数の(NxK)個の測定基準スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定する工程であって、前記モデルは、スペクトル線の形状を有する所定の関数に基づいており、かつ、それぞれの測定スキームに関連している、工程;
前記K個のモデルのそれぞれに、それぞれの測定スキームに対応するN個の測定参照スペクトルのそれぞれをフィッティングし、かつ、前記参照サンプルのそれぞれについて、前記数K個の測定スキームについて前記サンプルの参照スペクトルのベクトル表現を作成し、それによって前記各参照サンプルをそれぞれの成分のベクトルで表現する工程;
前記M個のグループのそれぞれ1つに対する前記サンプルの各々の対応を示す前記データを利用し、各グループについて、前記グループに関連するサンプルの成分のベクトルを分析し、かつ、前記グループの特性ベクトルを示すデータを決定する工程;および
前記参照サンプルの成分のベクトルと前記グループの特性ベクトルとの間の距離関数に基づいて、前記参照サンプルの成分のベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする距離関数の重みパラメータを決定し、それによって前記距離関数の前記重みパラメータの共通ベクトルを提供する工程;
それぞれのK個の測定スキームに対するK個のモデルを示すデータ、前記グループの特性ベクトルを示すデータ、および前記M個のグループに対する重みの共通ベクトルを示すデータを含むモデル化データを保存する工程であって、それによって、前記モデル化データを用いて対象のサンプルの未加工測定スペクトルデータをモデルベース分析することにより、対象のサンプルを分類して前記M個のグループの1つに関連付けることを可能にする、工程
を含む、方法。 - 前記対象のサンプルの分類を行う工程をさらに含み、該工程が、
前記対象のサンプルの未加工測定スペクトルデータに基づいて、前記K個の測定スキームの下で前記対象のサンプルのK個の測定スペクトルに対応するK個のデータピースをそれぞれ決定する工程;
前記モデルベース分析を前記K個のデータピースに適用する工程
を含み、該適用する工程が、
前記保存されたK個のモデルを使用して、対象のサンプルに対する前記K個の測定スペクトルのそれぞれを前記保存されたK個のモデルのそれぞれ1つにフィッティングし、かつ、前記K個の測定スペクトルのそれぞれに対するベストフィット条件に基づいて、前記K個の測定スキームのすべてのサンプルの複合ベクトル表現を作成する工程;
前記重みの共通ベクトルに前記距離関数を適用して、前記サンプルの複合ベクトル表現の前記グループの特性ベクトルのそれぞれまでの距離を決定し、かつ、前記決定された距離が最小となるグループに前記サンプルを関連付ける工程
を含むことを特徴とする、請求項1に記載の方法。 - 前記測定スキームの前記数Kが少なくとも2であることを特徴とする、請求項1または2に記載の方法。
- 前記モデルが、スペクトル線の形状の前記所定の関数および特定のピースワイズ多項式関数に基づく混合モデルとして構成されることを特徴とする、請求項1~3のいずれか一項に記載の方法。
- 前記距離関数が統計関数であることを特徴とする、請求項1~4のいずれか一項に記載の方法。
- 前記グループの特性ベクトルが、同一グループの前記参照サンプルを表す成分ベクトル内の成分の平均値を含むことを特徴とする、請求項1~5のいずれか一項に記載の方法。
- 前記距離関数が、前記ベクトルの成分の平均値および標準偏差に関連付けられ、それにより、前記成分ベクトル内の成分値の広がりの量が表現されることを特徴とする、請求項6に記載の方法。
- 前記K個のモデルを決定するための前記複数の(NxK)個の測定された参照スペクトルの処理が、
i番目の測定スキーム(i=1、…、K)に対応するN個の参照サンプルの複数の測定された参照スペクトルの各i番目に対して、平均測定参照スペクトルを決定する工程;および
スペクトル線の形状を有する前記所定の関数に従って各i番目の平均測定基準スペクトルに対して所定の変換を適用し、前記i番目の測定スキームに対応するそれぞれのi番目のモデルを取得し、それによりK個の測定スキームについてK個のモデルを取得する工程
を含むことを特徴とする、請求項1~7のいずれか一項に記載の方法。 - 前記所定の関数がガウス関数を含むことを特徴とする、請求項1~8のいずれか一項に記載の方法。
- 前記サンプルが、以下のタイプ:鉱物、貴石、ダイヤモンドの少なくとも1つであることを特徴とする、請求項1~9のいずれか一項に記載の方法。
- 前記M個のグループの所定の異なる特性が、以下:サンプルの原産地、およびサンプルの原産地の地理的な位置の1つまたは複数の構造パラメータ
の1つまたは複数を含むことを特徴とする、請求項10に記載の方法。 - 前記サンプルの測定スペクトルデータが、X線またはガンマ線放射に対する前記サンプルの蛍光X線(XRF)応答を示すことを特徴とする、請求項1~11のいずれか一項に記載の方法。
- サンプルにおける測定をモデル化するデータ分析システムであって、
数Kの測定スキームの下で、所定の異なる特性のM個のグループに関連する複数のN個の参照サンプルに対してスペクトル測定を行い、前記M個のグループに関連する複数の(NxK)個の測定参照スペクトルを含む測定参照データを生成するように構成および動作可能である測定システム;
前記測定参照データに基づいて、対象のサンプルのさらなる分類を可能にするモデル化データを決定するように構成および動作可能なコントロールシステム
を含み、前記コントロールシステムが、
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュールであって、前記モデルは、スペクトル線の形状を有する所定の関数に基づいており、それぞれの測定スキームに関連する;
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個のモデルのそれぞれについて、それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々をモデルにフィッティングし;参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現する;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの共通ベクトルを提供するように構成および動作可能な重み付けモジュール;および
それぞれのK個の測定スキームに対するK個のモデルを示すデータであって、グループの特性ベクトルを示すデータ、および、M個のグループに対する重みの共通ベクトルを示すデータを含む、保存されるモデル化データを生成するように構成および動作可能な出力ユーティリティ
を含むことを特徴とする、データ分析システム。 - サンプル分類システムであって、
数Kの測定スキームの下でサンプルのスペクトル測定を行い、かつ、測定されたサンプルのそれぞれについて、K個の測定スキームにそれぞれ対応する測定スペクトルを示すK個の測定データピースを含む測定スペクトルデータを生成するように構成および動作可能な測定システム;
前記測定システムと通信し対象の測定スペクトルデータを受信するように構成および動作可能であり、かつ、スペクトル線の形状を有する所定の関数に基づくそれぞれのK個の測定スキームに対するK個のモデルを示すデータと、異なるサンプルが関連するM個の所定のグループのM個の特性ベクトルを示すデータと、M個のグループに対する重みの共通ベクトルを示すデータとを含む所定のモデル化データを格納するメモリと通信するように構成および動作可能なコントロールシステムであって、前記所定のモデル化データを用いて対象のサンプルの受信した測定スペクトルデータにモデルベース処理を適用し、前記M個の所定のグループの1つに対する前記対象の特定のサンプルの関係を示す分類データを生成するように構成および動作可能なデータプロセッサを備える、コントロールシステム
を含むことを特徴とする、サンプル分類システム。 - 前記コントロールシステムが、
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個の測定スペクトルのそれぞれについて、測定スペクトルをそれぞれのモデルにフィッティングし、K個のベストフィット条件スペクトルを得る;および、前記K個のベストフィット条件スペクトルを用いて、前記K個の測定スキームすべてについて関心対象サンプルの複合ベクトル表現を作成する;
前記共通の重みベクトルを用いて所定の距離関数を利用し、前記M個のグループの前記M個の特性ベクトルのそれぞれまでの対象のサンプルの前記複合ベクトル表現の距離を決定し、決定された距離が最小となるグループに前記対象のサンプルを関連付けるように構成および動作可能な分類モジュール
を含むことを特徴とする、請求項14に記載のシステム。 - 前記コントロールシステムが、前記M個のグループに関連する複数のN個の参照サンプル上で実行された前記測定スキームの数Kのスペクトル参照測定に対応する測定スペクトルデータに基づいて、所定のモデル化データを決定するようにさらに構成および動作可能であり、ここで、前記スペクトル参照データは複数の(NxK)個の測定参照スペクトルを含み、前記M個のグループのそれぞれの1つに対する参照サンプルのそれぞれの対応関係を示すデータを含み、
前記コントロールシステムは、
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュール;
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個のモデルのそれぞれについて、それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々をモデルにフィッティングし;参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現する;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;および
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする前記所定の距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの前記共通ベクトルを提供するように構成および動作可能な重み付けモジュール
を含むことを特徴とする、請求項14または15に記載のシステム。 - サンプル分類の管理に使用するコントロールシステムであって、該コントロールシステムは、測定データプロバイダと通信し対象のサンプルの測定スペクトルデータを受信するように構成および動作可能であり、かつ、スペクトル線の形状を有する所定の関数に基づくそれぞれのK個の測定スキームに対するK個のモデルを示すデータと、異なるサンプルが関連するM個の所定のグループのM個の特性ベクトルを示すデータと、M個のグループに対する重みの共通ベクトルを示すデータとを含む所定のモデル化データを格納するメモリと通信するように構成および動作可能であり、前記コントロールシステムは、前記所定のモデル化データを用いて対象のサンプルの受信した測定スペクトルデータにモデルベース処理を適用し、前記M個の所定のグループの1つに対する前記対象の特定のサンプルの関係を示す分類データを生成するように構成および動作可能なデータプロセッサを備えることを特徴とする、コントロールシステム。
- 以下を実行するように構成および動作可能なフィッティングモジュール:前記K個の測定スペクトルのそれぞれについて、測定スペクトルをそれぞれのモデルにフィッティングし、K個のベストフィット条件スペクトルを得る;および、前記K個のベストフィット条件スペクトルを用いて、前記K個の測定スキームすべてについて関心対象サンプルの複合ベクトル表現を作成する;および
前記共通の重みベクトルを用いて所定の距離関数を利用し、前記M個のグループの前記M個の特性ベクトルのそれぞれまでの対象のサンプルの前記複合ベクトル表現の距離を決定し、決定された距離が最小となるグループに前記対象のサンプルを関連付けるように構成および動作可能な分類モジュール
を含むことを特徴とする、請求項17に記載のコントロールシステム。 - 前記M個のグループに関連する複数のN個の参照サンプル上で実行された前記測定スキームの数Kのスペクトル参照測定に対応する測定スペクトルデータに基づいて、所定のモデル化データを決定するようにさらに構成および動作可能であり、ここで、前記スペクトル参照データは複数の(NxK)個の測定参照スペクトルを含み、前記M個のグループのそれぞれの1つに対する参照サンプルのそれぞれの対応関係を示すデータを含む、コントロールシステムであって、
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュール;
以下を実行するように構成および動作可能なフィッティングモジュール:前記K個のモデルのそれぞれについて、それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々をモデルにフィッティングし;参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現する;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;および
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする前記所定の距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの前記共通ベクトルを提供するように構成および動作可能な重み付けモジュール
を含むことを特徴とする、請求項17または18に記載のコントロールシステム。 - 対象のサンプルをモデルベース分析するコントロールシステムであって、
所定の異なる特性を有するM個のグループに関連する複数のN個の参照サンプルに対して行われる数Kの測定スキームのスペクトル測定を示す参照データを受信するように構成および動作可能であるデータ入力ユーティリティであって、前記参照データは、複数の(NxK)個の測定された参照スペクトルを含む未加工測定データを含み、かつ、前記M個のグループのそれぞれ1つに対する前記各参照サンプルの対応関係を示すデータを含む、データ入力ユーティリティ;
前記複数の(NxK)個の測定された参照スペクトルを処理して、前記K個の測定スキームにそれぞれ対応するK個のモデルを決定するように構成および動作可能なモデル作成モジュールであって、前記モデルは、スペクトル線の形状を有する所定の関数に基づいており、それぞれの測定スキームに関連する;
前記それぞれの測定スキームに対応するN個の測定された参照スペクトルの各々を前記K個のモデルの各々にフィッティングするように構成および動作可能であり、かつ、参照サンプルのそれぞれについて、前記数Kの測定スキームに対するサンプルの参照スペクトルのベクトル表現を作成し、それによって参照サンプルの各々をそれぞれの成分ベクトルで表現するフィッティングモジュール;
前記M個のグループのそれぞれ1つに対する前記参照サンプルの各々の対応関係を示すデータを利用して、各グループについて、そのグループに関連するサンプルの成分ベクトルを分析し、そのグループの特性ベクトルを示すデータを決定するように構成および動作可能なグループ特性評価モジュール;
参照サンプルの成分ベクトルとグループの特性ベクトルとの間の距離関数に基づいて、参照サンプルの成分ベクトルのすべてをそれぞれのグループに関連付けるための複合尤度を最大にする距離関数の重みパラメータを決定し、それによって距離関数の前記重みパラメータの共通ベクトルを提供するように構成および動作可能な重み付けモジュール;
前記それぞれのK個の測定スキームに対するK個のモデルを示すデータと、前記グループの特性ベクトルを示すデータと、前記M個のグループに対する重みの共通ベクトルを示すデータとを含むモデル化データを格納するためのストレージユーティリティ;および
前記モデル化データを用いて対象のサンプルの未加工測定スペクトルデータをモデルベース分析することにより、対象のサンプルを分類し前記M個のグループの1つに関連付けるように構成および動作可能である分類モジュール
を含むことを特徴とする、コントロールシステム。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201962833956P | 2019-04-15 | 2019-04-15 | |
US62/833,956 | 2019-04-15 | ||
PCT/IL2020/050409 WO2020212969A1 (en) | 2019-04-15 | 2020-04-05 | Method and system for classification of samples |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2022529624A true JP2022529624A (ja) | 2022-06-23 |
JPWO2020212969A5 JPWO2020212969A5 (ja) | 2023-04-04 |
Family
ID=72837097
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021560551A Pending JP2022529624A (ja) | 2019-04-15 | 2020-04-05 | サンプルを分類する方法およびシステム |
Country Status (8)
Country | Link |
---|---|
US (1) | US20220317069A1 (ja) |
EP (1) | EP3956811A4 (ja) |
JP (1) | JP2022529624A (ja) |
KR (1) | KR20210151184A (ja) |
AU (1) | AU2020259413A1 (ja) |
CA (1) | CA3136899A1 (ja) |
IL (1) | IL287180A (ja) |
WO (1) | WO2020212969A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023236540A1 (zh) * | 2023-01-17 | 2023-12-14 | 山东大学 | 矿石成分分析设备及方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190103258A1 (en) * | 2016-03-07 | 2019-04-04 | Micromass Uk Limited | Spectrometric Analysis |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6266390B1 (en) * | 1998-09-21 | 2001-07-24 | Spectramet, Llc | High speed materials sorting using x-ray fluorescence |
US6140643A (en) * | 1999-03-09 | 2000-10-31 | Exxonmobil Upstream Research Company | Method for identification of unknown substances |
RU2008128455A (ru) * | 2005-12-12 | 2010-01-20 | Платформ Девелопмент Энд Инвестмент Лтд. (Il) | Оценка цвета алмазов |
WO2011159269A1 (en) * | 2010-06-17 | 2011-12-22 | Spectramet, Llc | Sorting pieces of material based on optical and x - ray photon emissions |
JP5986646B2 (ja) * | 2012-02-07 | 2016-09-06 | マテリアリティクス,エルエルシー | 試料を分析する方法およびシステム |
NL2009015C2 (en) * | 2012-04-10 | 2013-10-15 | Biosparq B V | Method for classification of a sample on the basis of spectral data, method for creating a database and method for using this database, and corresponding computer program, data storage medium and system. |
ES2929396T3 (es) | 2015-04-02 | 2022-11-29 | Soreq Nuclear Res Ct | Sistema y método para la lectura de marcación por fluorescencia de rayos X |
EP3472599B1 (en) | 2016-04-04 | 2022-06-01 | Soreq Nuclear Research Center | A method and a system for xrf marking and reading xrf marks of electronic systems |
JP7336382B2 (ja) | 2016-09-19 | 2023-08-31 | ソレク ニュークリア リサーチ センター | 蛍光x線システムおよび試料を識別する方法 |
DE102017217543B4 (de) * | 2017-10-02 | 2020-01-09 | Siemens Healthcare Gmbh | Verfahren und System zur Klassifikation von Materialien mittels maschinellen Lernens |
-
2020
- 2020-04-05 WO PCT/IL2020/050409 patent/WO2020212969A1/en unknown
- 2020-04-05 US US17/594,406 patent/US20220317069A1/en active Pending
- 2020-04-05 AU AU2020259413A patent/AU2020259413A1/en not_active Abandoned
- 2020-04-05 EP EP20791992.9A patent/EP3956811A4/en active Pending
- 2020-04-05 JP JP2021560551A patent/JP2022529624A/ja active Pending
- 2020-04-05 CA CA3136899A patent/CA3136899A1/en active Pending
- 2020-04-05 KR KR1020217036810A patent/KR20210151184A/ko unknown
-
2021
- 2021-10-11 IL IL287180A patent/IL287180A/en unknown
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190103258A1 (en) * | 2016-03-07 | 2019-04-04 | Micromass Uk Limited | Spectrometric Analysis |
Also Published As
Publication number | Publication date |
---|---|
KR20210151184A (ko) | 2021-12-13 |
EP3956811A4 (en) | 2023-01-11 |
IL287180A (en) | 2021-12-01 |
US20220317069A1 (en) | 2022-10-06 |
AU2020259413A1 (en) | 2021-11-04 |
WO2020212969A1 (en) | 2020-10-22 |
EP3956811A1 (en) | 2022-02-23 |
CA3136899A1 (en) | 2020-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104380088B (zh) | 对sem-eds数据集中的未知物的聚类分析 | |
US9128055B2 (en) | Information processing apparatus, information processing method, program, and method of correcting intensity of fluorescence spectrum | |
JP5870851B2 (ja) | 情報処理装置、情報処理方法、及びプログラム | |
KR102137184B1 (ko) | 자동 및 수동 결함 분류의 통합 | |
Cooper et al. | Measuring galaxy environments with deep redshift surveys | |
JP5264061B2 (ja) | 物質同定のための方法及び装置 | |
Tinker et al. | Halo histories versus Galaxy properties at z= 0–I. The quenching of star formation | |
TW201350836A (zh) | 用於自動缺陷分類之不明缺陷拒絕率之最佳化 | |
CN103959426B (zh) | 用于通过质谱术识别微生物的方法 | |
Walser et al. | On the parametrization of optical particle counter response including instrument-induced broadening of size spectra and a self-consistent evaluation of calibration measurements | |
CN109521725A (zh) | 检测异常数据的方法、装置和设备以及机器可读介质 | |
CN105593967A (zh) | 减小目标样品的厚度的方法 | |
JP2022529624A (ja) | サンプルを分類する方法およびシステム | |
CN108614296B (zh) | 观测系统重复性确定方法及装置 | |
EP2332073B1 (en) | Shape parameter for hematology instruments | |
CN106383135B (zh) | 基于全反射x射线荧光光谱的元素定量分析方法和系统 | |
US9476923B2 (en) | Method and device for identifying a material by the spectral analysis of electromagnetic radiation passing through said material | |
CN105844169A (zh) | 信息安全度量方法及装置 | |
CN112040405A (zh) | 一种基于核极限学习机与粒子滤波的室内定位方法 | |
US10677788B2 (en) | Reagent calibration system and method | |
CN108280429B (zh) | 基于特征峰匹配的核素识别方法、装置和计算机可读介质 | |
WO2022049854A1 (ja) | 全反射蛍光x線分析装置及び推定方法 | |
Mazumder et al. | A study on the clustering properties of radio-selected sources in the Lockman Hole region at 325 MHz | |
JP6949034B2 (ja) | 信号分析装置、信号分析方法、コンピュータプログラム、測定装置及び測定方法 | |
Mantz et al. | Erratum: Weighing the giants–V. Galaxy cluster scaling relations |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230327 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230327 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231122 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20231127 |