JP5623061B2 - 検査装置及び検査方法 - Google Patents
検査装置及び検査方法 Download PDFInfo
- Publication number
- JP5623061B2 JP5623061B2 JP2009260108A JP2009260108A JP5623061B2 JP 5623061 B2 JP5623061 B2 JP 5623061B2 JP 2009260108 A JP2009260108 A JP 2009260108A JP 2009260108 A JP2009260108 A JP 2009260108A JP 5623061 B2 JP5623061 B2 JP 5623061B2
- Authority
- JP
- Japan
- Prior art keywords
- expansion coefficient
- wavelet
- terahertz wave
- terahertz
- time waveform
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims description 41
- 238000007689 inspection Methods 0.000 title claims description 29
- 230000003595 spectral effect Effects 0.000 claims description 27
- 238000000605 extraction Methods 0.000 claims description 19
- 230000009466 transformation Effects 0.000 claims description 18
- 238000001514 detection method Methods 0.000 claims description 16
- 238000006243 chemical reaction Methods 0.000 claims description 15
- 230000001131 transforming effect Effects 0.000 claims description 13
- 230000001678 irradiating effect Effects 0.000 claims description 5
- 238000000862 absorption spectrum Methods 0.000 description 42
- 238000001228 spectrum Methods 0.000 description 34
- 238000003384 imaging method Methods 0.000 description 7
- 239000004038 photonic crystal Substances 0.000 description 7
- 238000010586 diagram Methods 0.000 description 5
- 239000000463 material Substances 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 229940088597 hormone Drugs 0.000 description 3
- 239000005556 hormone Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- RJKFOVLPORLFTN-LEKSSAKUSA-N Progesterone Chemical compound C1CC2=CC(=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H](C(=O)C)[C@@]1(C)CC2 RJKFOVLPORLFTN-LEKSSAKUSA-N 0.000 description 2
- CCGSUNCLSOWKJO-UHFFFAOYSA-N cimetidine Chemical compound N#CNC(=N/C)\NCCSCC1=NC=N[C]1C CCGSUNCLSOWKJO-UHFFFAOYSA-N 0.000 description 2
- 229960001380 cimetidine Drugs 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- VYFYYTLLBUKUHU-UHFFFAOYSA-N dopamine Chemical compound NCCC1=CC=C(O)C(O)=C1 VYFYYTLLBUKUHU-UHFFFAOYSA-N 0.000 description 2
- 239000003814 drug Substances 0.000 description 2
- HYYBABOKPJLUIN-UHFFFAOYSA-N mefenamic acid Chemical compound CC1=CC=CC(NC=2C(=CC=CC=2)C(O)=O)=C1C HYYBABOKPJLUIN-UHFFFAOYSA-N 0.000 description 2
- 229960003464 mefenamic acid Drugs 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 235000000346 sugar Nutrition 0.000 description 2
- 150000008163 sugars Chemical class 0.000 description 2
- 238000001328 terahertz time-domain spectroscopy Methods 0.000 description 2
- OWEGMIWEEQEYGQ-UHFFFAOYSA-N 100676-05-9 Natural products OC1C(O)C(O)C(CO)OC1OCC1C(O)C(O)C(O)C(OC2C(OC(O)C(O)C2O)CO)O1 OWEGMIWEEQEYGQ-UHFFFAOYSA-N 0.000 description 1
- KLSJWNVTNUYHDU-UHFFFAOYSA-N Amitrole Chemical compound NC1=NC=NN1 KLSJWNVTNUYHDU-UHFFFAOYSA-N 0.000 description 1
- 238000012935 Averaging Methods 0.000 description 1
- RKWGIWYCVPQPMF-UHFFFAOYSA-N Chloropropamide Chemical compound CCCNC(=O)NS(=O)(=O)C1=CC=C(Cl)C=C1 RKWGIWYCVPQPMF-UHFFFAOYSA-N 0.000 description 1
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 1
- GUBGYTABKSRVRQ-PICCSMPSSA-N Maltose Natural products O[C@@H]1[C@@H](O)[C@H](O)[C@@H](CO)O[C@@H]1O[C@@H]1[C@@H](CO)OC(O)[C@H](O)[C@H]1O GUBGYTABKSRVRQ-PICCSMPSSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- OIPILFWXSMYKGL-UHFFFAOYSA-N acetylcholine Chemical compound CC(=O)OCC[N+](C)(C)C OIPILFWXSMYKGL-UHFFFAOYSA-N 0.000 description 1
- 229960004373 acetylcholine Drugs 0.000 description 1
- RWCCWEUUXYIKHB-UHFFFAOYSA-N benzophenone Chemical compound C=1C=CC=CC=1C(=O)C1=CC=CC=C1 RWCCWEUUXYIKHB-UHFFFAOYSA-N 0.000 description 1
- 239000012965 benzophenone Substances 0.000 description 1
- WQZGKKKJIJFFOK-VFUOTHLCSA-N beta-D-glucose Chemical compound OC[C@H]1O[C@@H](O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-VFUOTHLCSA-N 0.000 description 1
- GUBGYTABKSRVRQ-QUYVBRFLSA-N beta-maltose Chemical compound OC[C@H]1O[C@H](O[C@H]2[C@H](O)[C@@H](O)[C@H](O)O[C@@H]2CO)[C@H](O)[C@@H](O)[C@@H]1O GUBGYTABKSRVRQ-QUYVBRFLSA-N 0.000 description 1
- 229960001761 chlorpropamide Drugs 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000875 corresponding effect Effects 0.000 description 1
- 229960003638 dopamine Drugs 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 239000008103 glucose Substances 0.000 description 1
- 229940048662 kwai Drugs 0.000 description 1
- MJIHNNLFOKEZEW-UHFFFAOYSA-N lansoprazole Chemical compound CC1=C(OCC(F)(F)F)C=CN=C1CS(=O)C1=NC2=CC=CC=C2N1 MJIHNNLFOKEZEW-UHFFFAOYSA-N 0.000 description 1
- 229960003174 lansoprazole Drugs 0.000 description 1
- 239000000186 progesterone Substances 0.000 description 1
- 229960003387 progesterone Drugs 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3581—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light using far infrared light; using Terahertz radiation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
- G06F2218/04—Denoising
- G06F2218/06—Denoising by applying a scale-space analysis, e.g. using wavelet analysis
-
- 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/16—Classification; Matching by matching signal segments
- G06F2218/18—Classification; Matching by matching signal segments by plotting the signal segments against each other, e.g. analysing scattergrams
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Toxicology (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
本発明は、テラヘルツ波を用いて対象物を検査するための検査装置及び検査方法に関する。なお、本明細書では、３０ＧＨｚ以上３０ＴＨｚ以下のうち少なくとも一部を含む周波数帯域を有する電磁波のことをテラヘルツ波と呼ぶ。
従来、サンプル（対象物）を透過または反射したテラヘルツ波の時間波形からサンプルの吸収率と屈折率の周波数依存性を調べるために、フーリエ解析が用いられていた。
しかしながら、テラヘルツ波の時間波形に含まれる広い周波数帯域において、吸収率と屈折率の周波数依存性は材料毎に特有である。このとき、比較するデータ量が多くなってしまうという問題があった。
この問題を解決するために、テラヘルツ波の時間波形に含まれる分光に関する情報を、少ないデータ量に圧縮し、該圧縮したデータから材料の構成成分を特定する方法が、特許文献１に開示されている。特許文献１には、この圧縮方法として、ウェーブレット解析を用いることが開示されている。
しかし、特許文献１には、具体的にどのようにウェーブレット解析を行うかに関しては開示されていない。特に、ウェーブレット展開係数についての開示もない。
本発明の目的は、比較的少ないデータ量の展開係数を用いて、対象物を検査（例えば、対象物の同定やイメージング）することである。
本発明の一側面としての検査装置は、対象物を検査するための検査装置であって、前記対象物にテラヘルツ波を照射する照射手段と、前記対象物から得られたテラヘルツ波を検出する検出手段と、前記検出手段で検出したテラヘルツ波を用いて得た、該テラヘルツ波の時間波形をウェーブレット変換する変換手段と、前記ウェーブレット変換における第１の展開係数から、前記対象物が持つ特徴的なスペクトル部分に関連する展開係数として該対象物の参照物体を用いて予め取得され記憶され且つ該第１の展開係数に含まれる第２の展開係数を選択する選択手段と、前記第２の展開係数における前記参照物体を用いて予め取得され記憶された第１の値と、前記第２の展開係数における前記検出手段で検出したテラヘルツ波の時間波形をウェーブレット変換することで取得した第２の値とを比較するための比較手段と、前記比較手段で得た結果から、前記対象物が前記参照物体であるかどうかを判定する判定手段と、を有することを特徴とする。
本発明の装置及び方法（検査装置及び検査方法）によれば、比較的少ないデータ量のウェーブレット展開係数を用いて対象物（サンプル）を検査すること（例えば、対象物の同定やイメージング）ができる。比較的少ないデータ量を扱うことにより、検査スピードを上げることができる。
（検査装置）
本実施形態に係る検査装置（対象物を検査するための検査装置）について、図１（ａ）を用いて説明する。
本実施形態に係る検査装置（対象物を検査するための検査装置）について、図１（ａ）を用いて説明する。
本実施形態に係る検査装置は、対象物１２にテラヘルツ波を照射する照射手段１１（例えば、光伝導素子などにより構成）を備える。
また、テラヘルツ波が照射された前記対象物１２からのテラヘルツ波を検出する検出手段１３（例えば、光伝導素子などにより構成）を備える。このとき、前記検出手段１３は、対象物１２から得られたテラヘルツ波を検出する。即ち、前記対象物１２で反射、散乱したテラヘルツ波、あるいは対象物１２を透過したテラヘルツ波を検出する。
また、前記検出手段１３で検出された信号に基づいて得られたテラヘルツ時間波形（例えば、図２（ａ）の波形のこと。）をウェーブレット変換する変換手段１４を備える。このとき、前記変換手段１４は、前記検出手段１３で検出したテラヘルツ波を用いて得た、該テラヘルツ波の時間波形（例えば、図２（ａ）の波形のこと。）をウェーブレット変換する。
さらに、前記ウェーブレット変換における第１の展開係数（例えば、図３（ａ）の横軸のこと。）から、予め記憶され（例えば、半導体メモリーに記憶される。）且つ該第１の展開係数に含まれる第２の展開係数（例えば、図３（ｂ）で求めた３０９番目のウェーブレット展開係数のこと。これは、吸収スペクトルなどの特徴的なスペクトル部分を持つ参照物体を用いて得た、該特徴的なスペクトル部分に関連する展開係数である。）を選択する選択手段１５を備える。
そして、前記第２の展開係数における第１の値（例えば、図３（ｂ）で求めた３０９番目のウェーブレット展開係数の場合は、３．１７のこと。）と、前記選択手段１５で選択された第２の展開係数における第２の値（前記対象物１２における３０９番目のウェーブレット展開係数における値のこと。）とを比較するための比較手段１６を備える。
これにより、典型的には、比較的少ないデータ量のウェーブレット展開係数を用いて対象物１２を検査すること（例えば、対象物の同定やイメージング）ができる。比較的少ないデータ量を扱うことにより、検査スピードを上げることができる。なお、上記検査とは、参照物体（例えば、医薬品）と、該参照物体と同一の種類である別の物体（上記対象物１２）と、を比較することを言う。これにより、対象物の同定を行うことができる。
このとき、前記比較手段１６で得た結果から、前記対象物１２が前記参照物体であるかどうかを判定する判定手段１７を備えることが好ましい。
（抽出手段）
この構成において、典型的には、前記特徴的なスペクトル部分に関連するウェーブレット展開係数を抽出するための抽出手段を有することが好ましい。この抽出手段により、予め記憶する前記特徴的なスペクトル部分を持つスペクトルのある物の該特徴的なスペクトル部分に関連するウェーブレット展開係数の係数位置と値を抽出する。
この構成において、典型的には、前記特徴的なスペクトル部分に関連するウェーブレット展開係数を抽出するための抽出手段を有することが好ましい。この抽出手段により、予め記憶する前記特徴的なスペクトル部分を持つスペクトルのある物の該特徴的なスペクトル部分に関連するウェーブレット展開係数の係数位置と値を抽出する。
また、この場合、抽出手段を用いて、前記対象物についても特徴的なスペクトル部分に関連するウェーブレット展開係数を抽出し、この抽出したウェーブレット展開係数の値と、前記予め記憶しておいたウェーブレット展開係数の値を比較してもよい。
前記抽出手段は、例えば、前記変換手段によるウェーブレット変換により分解されたウェーブレット展開係数の或る基準で選択された一部をゼロに置き換えてウェーブレット逆変換を行なう。そして、該ウェーブレット逆変換で得られたテラヘルツ時間波形に対してフーリエ変換を行ない、前記特徴的なスペクトル部分に変化を与えるウェーブレット展開係数を見つけ出す。
また、特徴的なスペクトルとの相関（影響度合い）の大きいウェーブレット展開係数を抽出しても良い。これは、特徴的なスペクトルとの相関の大きいウェーブレット展開係数の有無により、ウェーブレット逆変換により取得した時間波形が特徴的なスペクトルを有するか否かが決まる、という性質を利用する手法である。
前記比較や抽出が好適に行われる様に、前記変換手段は、例えば、テラヘルツ波形と相関の高いマザーウェーブレットを用い、かつ前記特徴的なスペクトル部分がウェーブレット展開係数の一部に関連付けられる様にテラヘルツ時間波形をウェーブレット変換する。
前記装置を用いて、対象物が特徴的なスペクトル部分を持つスペクトルのある物であるかどうかを判定し、この判定に基づいて対象物をスクリーニングするスクリーニングシステムを構築することができる。
（検査方法）
また、別の本実施形態に係る検査方法（対象物を検査するための検査方法）は、次のａ）からｅ）の各工程を含む。
ａ）対象物１２にテラヘルツ波を照射する照射工程１１。
ｂ）テラヘルツ波が照射された前記対象物１２からのテラヘルツ波を検出する検出工程１３。あるいは、前記対象物１２で反射、散乱あるいは透過したテラヘルツ波を検出する検出工程１３。
ｃ）前記検出工程１３で検出された信号に基づいて得られたテラヘルツ時間波形をウェーブレット変換する変換工程１４。あるいは、前記検出工程１３で検出したテラヘルツ波を用いて得た、該テラヘルツ波の時間波形をウェーブレット変換する変換工程１４。
ｄ）前記変換工程１４で得た第１の展開係数から、予め記憶した第２の展開係数（特徴的なスペクトル部分を持つ参照物体を用いて得た、該特徴的なスペクトル部分に関連する展開係数。）を選択する選択工程１５。
ｅ）前記選択工程１５で選択したウェーブレット展開係数の値と、前記変換工程１４で得られたウェーブレット展開係数の値のうちの前記抽出したウェーブレット展開係数と同じ係数位置のものの値を比較する比較工程１６。あるいは、前記予め記憶した第２の展開係数における第１の値と、前記選択工程１５で選択された第２の展開係数における第２の値とを比較するための比較工程１６。
また、別の本実施形態に係る検査方法（対象物を検査するための検査方法）は、次のａ）からｅ）の各工程を含む。
ａ）対象物１２にテラヘルツ波を照射する照射工程１１。
ｂ）テラヘルツ波が照射された前記対象物１２からのテラヘルツ波を検出する検出工程１３。あるいは、前記対象物１２で反射、散乱あるいは透過したテラヘルツ波を検出する検出工程１３。
ｃ）前記検出工程１３で検出された信号に基づいて得られたテラヘルツ時間波形をウェーブレット変換する変換工程１４。あるいは、前記検出工程１３で検出したテラヘルツ波を用いて得た、該テラヘルツ波の時間波形をウェーブレット変換する変換工程１４。
ｄ）前記変換工程１４で得た第１の展開係数から、予め記憶した第２の展開係数（特徴的なスペクトル部分を持つ参照物体を用いて得た、該特徴的なスペクトル部分に関連する展開係数。）を選択する選択工程１５。
ｅ）前記選択工程１５で選択したウェーブレット展開係数の値と、前記変換工程１４で得られたウェーブレット展開係数の値のうちの前記抽出したウェーブレット展開係数と同じ係数位置のものの値を比較する比較工程１６。あるいは、前記予め記憶した第２の展開係数における第１の値と、前記選択工程１５で選択された第２の展開係数における第２の値とを比較するための比較工程１６。
前記比較工程１６で得た結果から、前記対象物１２が前記参照物体であるかどうかを判定する判定工程を含むことが好ましい。
（抽出工程）
また、前記変換工程でのウェーブレット変換により分解されたウェーブレット展開係数の一部をゼロに置き換えてウェーブレット逆変換を行なう抽出工程を含むことが好ましい。このとき、該ウェーブレット逆変換で得られたテラヘルツ時間波形に対してフーリエ変換を行ない、前記特徴的なスペクトル部分に変化を与えるウェーブレット展開係数を見つけ出す。
また、前記変換工程でのウェーブレット変換により分解されたウェーブレット展開係数の一部をゼロに置き換えてウェーブレット逆変換を行なう抽出工程を含むことが好ましい。このとき、該ウェーブレット逆変換で得られたテラヘルツ時間波形に対してフーリエ変換を行ない、前記特徴的なスペクトル部分に変化を与えるウェーブレット展開係数を見つけ出す。
前記抽出工程では、例えば、しきい値以上或いはしきい値以下のウェーブレット展開係数の値をゼロに置き換えて、ウェーブレット逆変換などの前記処理を実行し、しきい値を変化させる毎に該処理を実行する。このことで前記特徴的なスペクトル部分に変化を与えるウェーブレット展開係数を見つけ出す。
また、特徴的なスペクトルとの相関（影響度合い）の大きいウェーブレット展開係数を抽出しても良い。これは、特徴的なスペクトルとの相関の大きいウェーブレット展開係数の有無により、ウェーブレット逆変換により取得した時間波形が特徴的なスペクトルを有するか否かが決まる、という性質を利用する手法である。
また、抽出工程では、指定する係数位置の範囲のウェーブレット展開係数の値をゼロに置き換えて前記処理を実行し、該範囲を変化させる毎に該処理を実行することで前記特徴的なスペクトル部分に変化を与えるウェーブレット展開係数を見つけ出すこともできる。
前記方法は、更に、次の様な保存工程と復元工程を含むことができる。保存工程では、前記比較工程で前記対象物が前記特徴的なスペクトル部分を持つスペクトルのある物であると判定したとき、信号成分を殆ど含まずノイズを表現しているウェーブレット展開係数の値をゼロもしくは小さい値に置き換える。そして、情報圧縮されたウェーブレット展開係数を保存する。復元工程では、前記保存工程で保存した、情報圧縮されたウェーブレット展開係数をウェーブレット逆変換し、前記対象物のテラヘルツ時間波形を復元する。この復元工程で復元した前記対象物のテラヘルツ時間波形に基づき、吸収率や屈折率などを求めて、対象物のイメージングを実行するイメージング工程を実行することもできる。この復元したテラヘルツ時間波形は、ＳＮ比が良くない高速に取得したデータを用いても、ＳＮ比が向上したものとできる。よって、高速に、良好なイメージングを実行することが可能となる。
（抽出方法）
また、別の本発明の実施形態について説明する。
また、別の本発明の実施形態について説明する。
特徴的なスペクトル部分（あるいは指紋スペクトル）に関連するウェーブレット展開係数を適当な方法で抽出する。この適当な方法は、ウェーブレット展開係数の変化が特徴的なスペクトル部分の変化をもたらすか否かを見て、もたらすウェーブレット展開係数を抽出すると言う方法を採る。すなわち、候補としてのウェーブレット展開係数をゼロにして、ウェーブレット逆変換し、これで得られたテラヘルツ時間波形をフーリエ変換する。そして、このフーリエ変換して得たスペクトルにおいて、特徴的なスペクトル部分が変化しているか否かを判定し、変化をもたらしたウェーブレット展開係数を目標のウェーブレット展開係数として抽出する。
また、別の本実施形態に係る抽出方法（特徴的なスペクトル部分を持つ参照物体から、該特徴的なスペクトル部分に関連する展開係数を抽出するための抽出方法。）は、次の様な工程を有する。
この方法は、取得工程と、変換工程と、抽出工程とを有する。前記取得工程では、特徴的なスペクトル部分を持つスペクトルのある対象物からのテラヘルツ波の時間波形を取得する。前記変換工程では、前記取得工程で得られたテラヘルツ時間波形をウェーブレット変換する（第１の変換工程）。前記抽出工程では、特徴的なスペクトル部分を持つスペクトルのある物の該特徴的なスペクトル部分に関連するウェーブレット展開係数を抽出する。また、抽出工程では、前記変換工程でのウェーブレット変換により分解されたウェーブレット展開係数の一部をゼロに置き換えてウェーブレット逆変換を行なう（逆変換工程）。そして、該ウェーブレット逆変換で得られたテラヘルツ時間波形に対してフーリエ変換を行なう（第２の変換工程）。そして、前記展開係数の一部を変えながら前記逆変換工程と前記第２の変換工程とを繰り返す。これにより、前記特徴的なスペクトル部分に変化を与えるウェーブレット展開係数を見つけ出す。
次に、特徴的なスペクトル部分を持つスペクトルの例である、特徴的な吸収スペクトルに関連するウェーブレット展開係数を用いた検査装置及び方法の具体的な実施例について、説明する。なお、特徴的な吸収スペクトルは、テラヘルツ領域では指紋スペクトルと呼ばれることもある。指紋スペクトルを得る場合、典型的には、３０ＧＨｚ以上３０ＴＨｚ以下の範囲内の０．１ＴＨｚ以上１０ＴＨｚ以下の周波数の成分を含むテラヘルツ波を用いる。
以下の実施例では、高速にデータ取得したテラヘルツ時間波形をウェーブレット変換によりウェーブレット展開係数に分解し、ウェーブレット展開係数のうち、特徴的な吸収スペクトルに関連するウェーブレット展開係数を抽出する方法について説明する。そして、この特徴的な吸収スペクトルに関連するウェーブレット展開係数を用いて、サンプルが特徴的な吸収スペクトルのある物質であるかどうかを判定する検査装置及び方法について説明する。
後述する本発明の実施例は、テラヘルツ領域に特徴的な吸収スペクトルを持つサンプルの検査に係り、ここではそのサンプルの一つとしてフォトニック結晶を例に説明する。使用したフォトニック結晶は、テラヘルツ領域に２つの特徴的な吸収スペクトルが存在するものである。
以下ではフォトニック結晶を用いた例について説明するが、本発明は、一般に、指紋スペクトルを持つ材料について適用することができる。例えば、医薬品（シメチジン、メフェナム酸、クロルプロパミド、ランソプラゾール）、ホルモン分子（プロゲステロン、ドーパミン、アセチルコリン、エステリオール）、環境ホルモン（アミトロール、ベンゾフェノン）、糖類（マルトース、グルコース）等にも適用できる。例として、シメチジンとメフェナム酸の指紋スペクトルを図４（ａ）、（ｂ）にそれぞれ示す。これらの例では、０．１ＴＨｚ付近から２．５ＴＨｚ付近の周波数の領域に特徴的な吸収スペクトルを持つことが分かる。
（実施例１）
図１（ｂ）は、実施例１の構成の概略図である。
図１（ｂ）は、実施例１の構成の概略図である。
テラヘルツ波照射手段１は、発生したテラヘルツ波をサンプル２に照射する。また、テラヘルツ波検出手段３は、サンプルを透過したテラヘルツ波を検出する。そして、パソコン４は、テラヘルツ波検出手段３で検出した信号を取り込む。
テラヘルツ波照射手段１としては、光伝導素子を用いてテラヘルツ波を発生させる手段を使用することができる。また、テラヘルツ波検出手段３についても、光伝導素子を用いてテラヘルツ波を検出する手段を使用することができる。ここで使用した光伝導素子では、周波数が０．１ＴＨｚ以上２．５ＴＨｚ以下の範囲のテラヘルツ波を発生、検出することができる。
光伝導素子を用いたテラヘルツ波の発生・検出方法は、ＴＨｚ−ＴＤＳ（テラヘルツ時間領域分光法）として知られており、テラヘルツ波検出手段３で得られる信号はテラヘルツ波の時間波形である。テラヘルツ波の時間波形は、ステージからなる遅延系（図１（ｂ）に明示せず）を走査することによって得られるが、ＳＮ比を向上させるため、ステージをゆっくり走査して、ステージ各点での信号を積算することが行なわれる。そのため、テラヘルツ波の時間波形を取得するのには時間を要する。
或いは、まず、ステージを高速に走査する。このとき、データを取得する各点において、平均化するために取得するデータ点数の比較的少ない（ＳＮ比の良くない）テラヘルツ波の時間波形を高速に取得することになる。次に、このデータ点数の比較的少ないテラヘルツ波の時間波形を繰り返し取得する。そして、これらのテラヘルツ波の時間波形を足し合わせて平均化し、ノイズを除去する（ＳＮ比を向上させる）ことも行なわれている。この場合、テラヘルツ波の時間波形を繰り返し測定するので、平均化して最終的にＳＮ比の良いテラヘルツ波の時間波形にするのには時間を要する。
従来は、測定したサンプルに特徴的な吸収スペクトルがあるかどうかを判定するためには、この様にして時間をかけて測定したＳＮ比の良いテラヘルツ波の時間波形をフーリエ変換し、スペクトル分布を比較することが行われる。しかし、この方法は、周波数範囲が広いため、比較されるべきデータ量が多くなり応用上実際的でない。この様に従来の方法は、フーリエ変換に使用するためのテラヘルツ波の時間波形を取得するのに時間がかかる。また、測定したサンプルに特徴的な吸収スペクトルがあるかどうかを判定するために使用する比較されるべきデータ量が多くなってしまう。
そこで、本発明では高速に取得したテラヘルツ波の時間波形でも、測定したサンプルに特徴的な吸収スペクトルがあるかどうかを判定することができる様にする。すなわち、例えば、遅延系のステージを高速に走査して、データ点数の比較的少ない（ＳＮ比の良くない）テラヘルツ波の時間波形を短時間に取得し、そのデータをパソコン４に取り込み、ウェーブレット変換する手段５（変換手段）でウェーブレット展開係数に分解する。テラヘルツ時間波形と相関の高いマザーウェーブレットを使用することで、信号成分であるテラヘルツ波とノイズを比較的容易に分離することができる。「相関の高い」とは、できるだけ少ない数のウェーブレット展開係数で表せる様にテラヘルツ時間波形と似た形状を持つことを意味する。ここでは、テラヘルツ波の時間波形のデータの点数ＬｓはＬｓ＝４０９６（＝２１２）点のものを使用し、マザーウェーブレットはＣｏｉｆｌｅｔ４を用いた。また、ウェーブレット展開係数の数を決めるレベルは、この場合の最大数であるｌｏｇ２Ｌｓ＝１２を使用した。
次に、特徴的な吸収スペクトルに関連するウェーブレット展開係数を抽出する手段６（抽出手段）における抽出方法について述べる。ウェーブレット変換により得られたウェーブレット展開係数をそのまま（各ウェーブレット展開係数の値を何も変えずに）ウェーブレット逆変換すると、ウェーブレット変換する前のテラヘルツ波の時間波形が得られる。しかし、ここでは次のような操作を行なう。まず、ゼロに近い値（比較的小さい値）にしきい値を設定する。そして、例えば、設定したしきい値と同じか、それよりも小さいウェーブレットの展開係数を用いてウェーブレット逆変換する。すなわち、設定したしきい値よりも大きい値をゼロに置き換えて、値を持つ（ゼロに置き換えられなかった）ウェーブレットの展開係数を用いてウェーブレット逆変換する。そして、ウェーブレット逆変換により得られた波形をフーリエ変換する。
同様に、設定したしきい値の値を少し大きくして同じ処理を行なう。しきい値の値を少しずつ大きくして、フーリエ変換によるスペクトルの変化を見ていくと、しきい値の値が変わったときに特徴的な吸収スペクトルに大きな変化、例えば、特徴的な吸収スペクトルが消える、あるいは現れる、といったことが起きる場合がある。そのときの、前後のしきい値に対応するウェーブレット展開係数の差を算出すると、特徴的な吸収スペクトルに大きな変化を与えているウェーブレット展開係数を見つけ出すことができる。この様にして見つけた、特徴的な吸収スペクトルに大きな変化を与えているウェーブレット展開係数の位置と値をメモリーに保存しておく。
図２（ａ）は、サンプルのないときのテラヘルツ波の時間波形である参照波形（破線）と、テラヘルツ領域に特徴的な吸収スペクトルを持つフォトニック結晶をサンプルにしたときのテラヘルツ波の時間波形（実線）を示したものである。
図２（ｂ）は、これらの波形をそれぞれフーリエ変換したときのスペクトルを示したものである。１ＴＨｚ付近と１．６ＴＨｚ付近に落ち込みのあるスペクトルが、フォトニック結晶をサンプルとしたときのテラヘルツ波の時間波形をフーリエ変換したもの（実線）で、吸収スペクトルのないスペクトルが、参照波形をフーリエ変換したもの（破線）である。
図３（ａ）は、フォトニック結晶をサンプルにしたときのテラヘルツ波の時間波形をウェーブレット変換したときのウェーブレット展開係数の一部を図示したものである。前述した様に、マザーウェーブレットはＣｏｉｆｌｅｔ４を使用し、レベルはｌｏｇ２ Ｌｓ＝１２（信号のデータの点数Ｌｓ＝４０９６）を使用した。図３（ａ）における縦軸の５と−５付近にある四角で囲った領域は、絶対値が同じ値のしきい値を少しずつ変化させていったとき、その変化前後のウェーブレット展開係数の差を模式的に図示したものである。もし、絶対値が同じ値のしきい値を変化させた前後のスペクトルに差があれば、そのスペクトルの変化は、この四角で囲った領域に含まれるウェーブレット展開係数によるものである。そのため、しきい値をいろいろと変化させてこの四角で囲った領域を変化させ、その中に含まれるウェーブレット展開係数を抽出することで、ウェーブレット展開係数とスペクトルの変化の関係を調べることができる。しきい値の変化の幅が大きければこの領域に含まれるウェーブレット展開係数は複数になることもある。しきい値の変化の幅を小さく取ることによって、四角で囲った領域を小さくすることができ、特徴的な吸収スペクトルに関連するウェーブレット展開係数の数を絞り込むことができる。特徴的な吸収スペクトルに関連するウェーブレット展開係数の数は、データ処理速度の向上の観点からは少なければ少ないほど都合がよい。
前記フォトニック結晶について、特徴的な吸収スペクトルに関連するウェーブレット展開係数をこの様にして調べた結果、ウェーブレット展開係数３０９番目の値が、特徴的な吸収スペクトルに関連するウェーブレット展開係数であることが分かった。特徴的な吸収スペクトルがあるときは展開係数３０９番目の値は３．１７であるが、この値をゼロにすると特徴的な吸収スペクトルの吸収スペクトルは劇的に消滅する。図３（ｂ）の実線は、展開係数３０９番目の値が３．１７のときで、特徴的な吸収スペクトルがあるときのスペクトルを表している。
一方、図３（ｂ）の破線は展開係数３０９番目の値をゼロにしたときのスペクトルであり、スペクトルが劇的に消滅している様子が分かる。従って、この様にまず、特徴的な吸収スペクトルのあるサンプルについて、予め特徴的な吸収スペクトルに関連するウェーブレット展開係数の位置と値を調べておく。そして、未知のサンプルについて、同じマザーウェーブレットを使用して同じウェーブレット変換を行ない、予め調べておいた係数位置の特徴的な吸収スペクトルに関連するウェーブレット展開係数の値を取り出す。その値が特徴的な吸収スペクトルのある物の値とほぼ同じ値かどうかを比較手段７で調べることによって、サンプルの検査を行なうことができる。
検査によって、未知のサンプルが特徴的な吸収スペクトルのあるサンプルであると判断した場合にのみ、ノイズを表現しているウェーブレット展開係数の値をゼロもしくは小さい値に置き換えて、ウェーブレット展開係数をパソコンのメモリーに保存すればよい。こうすることで、特徴的な吸収スペクトルのあるサンプルかどうかの検査と、ノイズ除去、情報圧縮が同時に実現できる。この保存した、情報圧縮されたウェーブレット展開係数をウェーブレット逆変換すれば、ノイズの除去された特徴的な吸収スペクトルを持つテラヘルツ波の時間波形を復元することができる。これに基づき、ＳＮの良い対象物のイメージングなどを実現することができる。
なお、本実施例ではサンプル２を透過したテラヘルツ波をテラヘルツ波検出手段３で検出する構成としたが、サンプル２で反射または散乱したテラヘルツ波をテラヘルツ波検出手段３で検出するような構成としてもよい。また、マザーウェーブレットとしてＣｏｉｆｌｅｔ４を使用したが、テラヘルツ時間波形と相関の高いマザーウェーブレットが存在すれば他のマザーウェーブレットを使用してもよい。さらに、ウェーブレットの展開レベルも、上記例ではｌｏｇ２Ｌｓ＝１２（Ｌｓは信号のデータの点数）を使用したが、特徴的な吸収スペクトルに関連するウェーブレット展開係数を見つけ出すことができれば、レベルは１２より小さい値でもよい。
先に述べた医薬品、ホルモン、糖類などの指紋スペクトルを持つ材料に本ウェーブレット変換を用いた検査装置及び方法を適用することもできる。
（実施例２）
実施例２は、特徴的な吸収スペクトルに関連するウェーブレット展開係数の位置と値を、実施例１のようなしきい値を使わずに求めるものである。実施例１ではしきい値を使用して特徴的な吸収スペクトルに関連するウェーブレット展開係数を求めたが、これは、テラヘルツ波の時間波形とマザーウェーブレットとの相関度すなわちウェーブレット展開係数の大きさに着目した方法といえる。実施例２では、ウェーブレット展開係数における周波数（すなわち係数の位置）に着目して、特徴的な吸収スペクトルに関連するウェーブレット展開係数の位置と値を見つけていく。
実施例２は、特徴的な吸収スペクトルに関連するウェーブレット展開係数の位置と値を、実施例１のようなしきい値を使わずに求めるものである。実施例１ではしきい値を使用して特徴的な吸収スペクトルに関連するウェーブレット展開係数を求めたが、これは、テラヘルツ波の時間波形とマザーウェーブレットとの相関度すなわちウェーブレット展開係数の大きさに着目した方法といえる。実施例２では、ウェーブレット展開係数における周波数（すなわち係数の位置）に着目して、特徴的な吸収スペクトルに関連するウェーブレット展開係数の位置と値を見つけていく。
図３（ｃ）において、ウェーブレット展開係数は位置の値が大きくなればなるほど波の周波数は高くなる。そのため図３（ｃ）において、５００付近の位置にあるウェーブレット展開係数と１００付近の位置にあるウェーブレット展開係数では、５００付近にあるウェーブレット展開係数の方が、波に含まれる周波数は高い。そこで、ウェーブレット展開係数の値ではなくウェーブレット展開係数の位置を指定して、その位置のウェーブレット展開係数の値をゼロにしてウェーブレット逆変換を行ない、その後、フーリエ変換で特徴的な吸収スペクトルの変化があるかどうかを調べる。このとき、ゼロにするウェーブレット展開係数の位置は１つでもよいし、例えば５５０から６００の位置にある全てのウェーブレット展開係数の値をゼロにする様に、位置の範囲で指定してもよい。或いは、指定する範囲も一ヵ所だけでなく複数の範囲を指定してもよい。いずれにしても、特徴的な吸収スペクトルに関連するウェーブレット展開係数をこのような手順で絞り込んでいき、特徴的な吸収スペクトルに関連するウェーブレット展開係数を見つけ出す。
実施例２では特徴的な吸収スペクトルに関連するウェーブレット展開係数の位置と値を見つけ出す方法が実施例１とは異なるが、その他については実施例１と同じである。
（その他の実施例）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。なお、記憶媒体は、コンピュータに実行させるためのプログラムを格納できるものであれば何でも良い。
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。なお、記憶媒体は、コンピュータに実行させるためのプログラムを格納できるものであれば何でも良い。
１１ 照射手段
１２ 対象物
１３ 検出手段
１４ 変換手段
１５ 選択手段
１６ 比較手段
１７ 判定手段
１２ 対象物
１３ 検出手段
１４ 変換手段
１５ 選択手段
１６ 比較手段
１７ 判定手段
Claims (4)
- 対象物を検査するための検査装置であって、
前記対象物にテラヘルツ波を照射する照射手段と、
前記対象物から得られたテラヘルツ波を検出する検出手段と、
前記検出手段で検出したテラヘルツ波を用いて得た、該テラヘルツ波の時間波形をウェーブレット変換する変換手段と、
前記ウェーブレット変換における第１の展開係数から、前記対象物が持つ特徴的なスペクトル部分に関連する展開係数として該対象物の参照物体を用いて予め取得され記憶され且つ該第１の展開係数に含まれる第２の展開係数を選択する選択手段と、
前記第２の展開係数における前記参照物体を用いて予め取得され記憶された第１の値と、前記第２の展開係数における前記検出手段で検出したテラヘルツ波の時間波形をウェーブレット変換することで取得した第２の値とを比較するための比較手段と、
前記比較手段で得た結果から、前記対象物が前記参照物体であるかどうかを判定する判定手段と、
を有することを特徴とする検査装置。 - 対象物を検査するための検査方法であって、
前記対象物にテラヘルツ波を照射する照射工程と、
前記対象物から得られたテラヘルツ波を検出する検出工程と、
前記検出工程で検出したテラヘルツ波を用いて得た、該テラヘルツ波の時間波形をウェーブレット変換する変換工程と、
前記ウェーブレット変換における第１の展開係数から、前記対象物が持つ特徴的なスペクトル部分に関連する展開係数として該対象物の参照物体を用いて予め取得され記憶され且つ該第１の展開係数に含まれる第２の展開係数を選択する選択工程と、
前記第２の展開係数における前記参照物体を用いて予め取得され記憶された第１の値と、前記第２の展開係数における前記検出工程で検出したテラヘルツ波の時間波形をウェーブレット変換することで取得した第２の値とを比較するための比較工程と、
前記比較工程で得た結果から、前記対象物が前記参照物体であるかどうかを判定する判定工程と、
を含むことを特徴とする検査方法。 - 請求項２に記載の検査方法をコンピュータに実行させることを特徴とするプログラム。
- 特徴的なスペクトル部分を持つ参照物体から、該特徴的なスペクトル部分に関連する展開係数を抽出するための抽出方法であって、
前記参照物体で反射、散乱または透過したテラヘルツ波の時間波形を取得する取得工程と、
前記取得工程で得たテラヘルツ時間波形をウェーブレット変換する第１の変換工程と、前記ウェーブレット変換における展開係数のうち、該展開係数の一部をウェーブレット逆変換する逆変換工程と、
前記逆変換工程で得たテラヘルツ時間波形に対してフーリエ変換を行なう第２の変換工程と、を含み、
前記逆変換工程では、前記展開係数の一部を指定するときにしきい値を設定し、該設定されたしきい値以上またはしきい値以下の展開係数の値をゼロに置き換え、
前記展開係数の一部を置き換えながら前記逆変換工程と前記第２の変換工程とを繰り返すことを特徴とする抽出方法。
Priority Applications (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009260108A JP5623061B2 (ja) | 2008-12-12 | 2009-11-13 | 検査装置及び検査方法 |
| US12/626,251 US8258477B2 (en) | 2008-12-12 | 2009-11-25 | Inspection apparatus and inspection method |
| EP20090178633 EP2196791A3 (en) | 2008-12-12 | 2009-12-10 | Terahertz inspection apparatus and method |
| CN2009102585146A CN101769864B (zh) | 2008-12-12 | 2009-12-11 | 检查装置和检查方法 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2008316415 | 2008-12-12 | ||
| JP2008316415 | 2008-12-12 | ||
| JP2009260108A JP5623061B2 (ja) | 2008-12-12 | 2009-11-13 | 検査装置及び検査方法 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2010160136A JP2010160136A (ja) | 2010-07-22 |
| JP2010160136A5 JP2010160136A5 (ja) | 2012-12-27 |
| JP5623061B2 true JP5623061B2 (ja) | 2014-11-12 |
Family
ID=41528704
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2009260108A Expired - Fee Related JP5623061B2 (ja) | 2008-12-12 | 2009-11-13 | 検査装置及び検査方法 |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US8258477B2 (ja) |
| EP (1) | EP2196791A3 (ja) |
| JP (1) | JP5623061B2 (ja) |
| CN (1) | CN101769864B (ja) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5552321B2 (ja) * | 2010-01-08 | 2014-07-16 | キヤノン株式会社 | 電磁波の測定装置及び方法 |
| JP5596982B2 (ja) * | 2010-01-08 | 2014-10-01 | キヤノン株式会社 | 電磁波の測定装置及び方法 |
| JP5743453B2 (ja) * | 2010-05-18 | 2015-07-01 | キヤノン株式会社 | テラヘルツ波の測定装置及び測定方法 |
| BG66598B1 (bg) * | 2011-02-25 | 2017-08-31 | Габрово, Технически Университет | Метод за откриване и определяне на вещества или смеси и/или определяне на техни характеристики |
| US8913688B1 (en) | 2012-04-30 | 2014-12-16 | Xilinx, Inc. | Pre-distortion for a phase interpolator with nonlinearity |
| JP6195052B2 (ja) * | 2013-03-08 | 2017-09-13 | セイコーエプソン株式会社 | 標本検査装置 |
| JP2014209094A (ja) * | 2013-03-29 | 2014-11-06 | キヤノン株式会社 | テラヘルツ波を用いて試料の情報を取得する情報取得装置および情報取得方法 |
| US20200296265A1 (en) * | 2019-03-14 | 2020-09-17 | Canon Kabushiki Kaisha | Processing system |
| CN111950351B (zh) * | 2020-06-29 | 2024-03-08 | 北京农业智能装备技术研究中心 | 基于太赫兹和可见光的农机劳损早期诊断巡检系统 |
| JP2023047720A (ja) * | 2021-09-27 | 2023-04-06 | キヤノン株式会社 | 検査システム、検査システムの制御方法、プログラムおよび記憶媒体 |
| CN115082339B (zh) * | 2022-06-22 | 2024-10-18 | 西安电子科技大学 | 基于连通域自适应划分的太赫兹图像去噪方法 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5939721A (en) * | 1996-11-06 | 1999-08-17 | Lucent Technologies Inc. | Systems and methods for processing and analyzing terahertz waveforms |
| US6078047A (en) * | 1997-03-14 | 2000-06-20 | Lucent Technologies Inc. | Method and apparatus for terahertz tomographic imaging |
| JP3561151B2 (ja) * | 1998-07-27 | 2004-09-02 | 株式会社東芝 | 異常診断装置および異常診断方法 |
| JP4283391B2 (ja) * | 1999-09-29 | 2009-06-24 | 高砂熱学工業株式会社 | ベアリングの異常検出方法 |
| JP2006275910A (ja) * | 2005-03-30 | 2006-10-12 | Canon Inc | 位置センシング装置及び位置センシング方法 |
| CN101210873A (zh) * | 2006-12-31 | 2008-07-02 | 清华大学 | 一种利用太赫兹时域光谱快速检测植物油纯度的方法及设备 |
| CN101251492B (zh) * | 2008-01-02 | 2011-06-22 | 阮双琛 | 一种连续波太赫兹实时成像装置及其方法 |
| US20110184654A1 (en) * | 2008-09-17 | 2011-07-28 | Opticul Diagnostics Ltd. | Means and Methods for Detecting Bacteria in an Aerosol Sample |
-
2009
- 2009-11-13 JP JP2009260108A patent/JP5623061B2/ja not_active Expired - Fee Related
- 2009-11-25 US US12/626,251 patent/US8258477B2/en not_active Expired - Fee Related
- 2009-12-10 EP EP20090178633 patent/EP2196791A3/en not_active Withdrawn
- 2009-12-11 CN CN2009102585146A patent/CN101769864B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| US20100148071A1 (en) | 2010-06-17 |
| CN101769864A (zh) | 2010-07-07 |
| CN101769864B (zh) | 2012-08-08 |
| JP2010160136A (ja) | 2010-07-22 |
| US8258477B2 (en) | 2012-09-04 |
| EP2196791A2 (en) | 2010-06-16 |
| EP2196791A3 (en) | 2011-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5623061B2 (ja) | 検査装置及び検査方法 | |
| CN103822879B (zh) | 一种基于高光谱成像技术的猕猴桃膨大果无损检测方法 | |
| CN103837492B (zh) | 一种基于近红外光谱技术的猕猴桃膨大果无损检测方法 | |
| Reddy et al. | Accurate histopathology from low signal-to-noise ratio spectroscopic imaging data | |
| JP2010160136A5 (ja) | ||
| CN108444934A (zh) | 一种复合材料脱粘缺陷的自动分割及量化方法 | |
| WO2017134057A1 (fr) | Procédé et dispositif de reconstruction d'un signal utile à partir d'un signal acquis bruité | |
| JP5596982B2 (ja) | 電磁波の測定装置及び方法 | |
| CN108827904B (zh) | 基于太赫兹光谱的物质识别方法、装置、设备及存储介质 | |
| CN116399850A (zh) | 一种用于光信号处理的光谱探测识别系统及其探测方法 | |
| CN116148210A (zh) | 水果品质检测方法及装置 | |
| CN105203495A (zh) | 一种基于希尔伯特-黄变换的光谱信号去噪方法 | |
| Kažys et al. | Ultrasonic detection of defects in strongly attenuating structures using the Hilbert–Huang transform | |
| CN114002203A (zh) | 一种基于拉曼光谱分析木质成分含量的方法及装置 | |
| CN116385296B (zh) | 基于自适应软阈值收缩的太赫兹成像方法 | |
| JP5552321B2 (ja) | 電磁波の測定装置及び方法 | |
| FR3105527A1 (fr) | PROCEDE DE surveillance d’UNe structure par traitement d’imageS | |
| Roth et al. | Signal processing approaches for Terahertz data obtained from inspection of the shuttle external tank thermal protection system foam | |
| Melchiorre et al. | Compact Convolutional Transformer Fourier analysis for GPR tunnels assessment | |
| CN118329833B (zh) | 水果品质连续化检测方法 | |
| CN119274698A (zh) | 基于高光谱成像的果蔬糖度检测模型构建及糖度检测方法 | |
| Burrascano et al. | Near-infrared nondestrutive evaluation of food based on multi-frequency Lock-in detection | |
| Zhang et al. | Nondestructive determination of the soluble solid content for tomato using hyperspectral diffuse transmittance imaging | |
| Zhou et al. | Pork grade evaluation using hyperspectral imaging techniques | |
| CN119757267A (zh) | 一种绝缘油故障特征气体原位无损检测方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20100630 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20121113 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20121113 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20130830 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140107 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140310 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140603 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140623 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20140826 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140924 |
|
| LAPS | Cancellation because of no payment of annual fees |