JPH0484280A - High precision peak recognition method - Google Patents
High precision peak recognition methodInfo
- Publication number
- JPH0484280A JPH0484280A JP2200279A JP20027990A JPH0484280A JP H0484280 A JPH0484280 A JP H0484280A JP 2200279 A JP2200279 A JP 2200279A JP 20027990 A JP20027990 A JP 20027990A JP H0484280 A JPH0484280 A JP H0484280A
- Authority
- JP
- Japan
- Prior art keywords
- peak
- value
- differential
- point
- inflection point
- 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
Landscapes
- Analysing Materials By The Use Of Radiation (AREA)
- Measurement Of Radiation (AREA)
- Image Analysis (AREA)
Abstract
(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。(57) [Summary] This bulletin contains application data before electronic filing, so abstract data is not recorded.
Description
本発明は、EPMA等によって得られるスペクトルデー
タからスペクトルピークを高精度(高信頼性)で認定す
る方法に関する。The present invention relates to a method for identifying spectral peaks with high accuracy (high reliability) from spectral data obtained by EPMA or the like.
EPMA等で得られるスペクトルには、測定試料からの
特性X線から得られる検出信号以外に、ノイズや周辺か
ら入射する他のX線からの検出信号が含まれており、ま
た、特性X線からの信号でも、ランダムに入射するため
に、統計的なゆらぎがあって、充分な積分時間をかけな
ければ、滑らかな曲線のスペクトルデータにはならない
。通常は、ピークが判定できる程度の積分時間で測定を
行うために、得られるスペクトルは、ギザギザ状となっ
ている。そこでスペクトルデータは、先ずスムージング
処理がなされギザギザのない理想的な曲線状とした後に
、バックグランドから突出した部分をスペクトルピーク
と認定していた。In addition to the detection signal obtained from the characteristic X-rays from the measurement sample, the spectrum obtained by EPMA etc. includes detection signals from noise and other X-rays incident from the surroundings, and also contains detection signals from the characteristic X-rays. Even in the case of a signal, there are statistical fluctuations due to random incidence, and unless a sufficient integration time is taken, spectral data with a smooth curve will not be obtained. Usually, measurement is performed with an integration time long enough to determine the peak, so the resulting spectrum has a jagged shape. Therefore, the spectral data was first subjected to smoothing processing to make it into an ideal curved shape without jagged edges, and then the portion that protruded from the background was recognized as the spectral peak.
上述した従来例は、単にスムージングされたスベクトル
データで突出している部分を、ピークと認定するだけな
ので、偶然存在したノイズのピークを眞のピークと誤認
したり、逆に、ピークと認定するための選別レベル以下
の眞のピークを見逃すと云うことがあって、信頼性が充
分でなかった本発明は、再度チエツクすることにより、
検出ピークの信頼性をより高めることを目的とする。In the conventional example described above, a prominent part of smoothed svector data is simply recognized as a peak, so a noise peak that happens to be present may be mistaken as a real peak, or conversely, it may be recognized as a peak. The present invention, which was not reliable enough because it sometimes missed true peaks below the screening level, could be improved by checking again.
The purpose is to further increase the reliability of detected peaks.
X線分光分析等で得てスムージングされたスペクトルデ
ータを微分して、微分値が十から−に変わる地点く微分
変曲点)を検索し、その微分変曲点をピーク頂点として
取り合えず認定し、同ピーク頂点における検出信号強度
値を(PK)とし、同(PK)から半値幅を推定し、上
記ピーク頂点から両側に上記半値幅のn倍離れた地点の
検出信号強度値を、バックグランド値[(BG−)、(
BG+)]とし、上記ピーク信号強度値(PK)が、
(F’K) −[kX E G−) 士、(BG
+)J/2±+ (B(、−) + (BG+) )
/2] >0を満足させる時、(PK)をピーク値とし
て認定するようにした。Differentiate the smoothed spectrum data obtained by X-ray spectroscopy, etc., search for the point where the differential value changes from 10 to - (differential inflection point), and identify that differential inflection point as the peak apex. , the detected signal strength value at the peak apex is (PK), the half-width is estimated from the same (PK), and the detected signal strength values at points n times the half-width on both sides from the peak apex are calculated as background. Value [(BG-), (
BG+)], and the above peak signal strength value (PK) is (F'K) - [kX E G-), (BG
+)J/2±+ (B(,-) + (BG+) )
/2] When >0 is satisfied, (PK) is recognized as the peak value.
スペクトルデータをスムージング処理し、スムージング
処理したスペクトルデータを微分して、微分変曲点を検
索し、検索された微分変曲点に存在するピークが、検出
ピークであるがどうかを判定する。この判定は、スムー
ジングされたスペクトルのベースラインより、一定値(
選別レベル)以上突出していればピークと認定すると云
う方法でピークを検出するのではなく、微分信号の正負
反転によってピークを検出するので、微小なピークでも
見逃さない。他方、これではピークの過剰検出の恐れも
ある。そこで、本発明では、更に、次のピーク検定を行
い、2重検出でピークを確認しているのである。
成る物質に電子線を照射し、発生する放射線を検出器で
検出する放射線測定では、検出信号には統計変動誤差が
存在する。そして、この変動誤差は、直値に対して、ガ
ウス分布をなす、従って、放射線の測定値に対して、誤
差標準偏差σは、σ=F10The spectral data is smoothed, the smoothed spectral data is differentiated, a differential inflection point is searched, and it is determined whether a peak present at the searched differential inflection point is a detected peak. This judgment is made by comparing the baseline of the smoothed spectrum with a constant value (
Rather than detecting peaks by recognizing them as peaks if they stand out above (selection level), peaks are detected by inverting the sign of the differential signal, so even minute peaks are not overlooked. On the other hand, this may lead to excessive detection of peaks. Therefore, in the present invention, the following peak verification is further performed to confirm the peak by double detection. In radiation measurement, in which a material is irradiated with an electron beam and the generated radiation is detected by a detector, there is a statistical fluctuation error in the detected signal. This fluctuation error has a Gaussian distribution with respect to the direct value. Therefore, the error standard deviation σ for the radiation measurement value is σ = F10
【「玉−
となる。このことより、測定値から1.5σの範囲内を
とると約9割(87%)の確からしさとなることから、
測定値の1.5σ範囲内がバックグランドの1.5σ範
囲内と交わらなくなった時に、測定値がバックグランド
と有為差があると判定されるので、以下の判定式で判定
行うことができる。
[<PK)−1,5f−でゴ「x−T] −E (
(BG−) +(BG十ン)/2+1.5x、/7r丁
丁=T下で]丁マ「1−Σ17二Σコ 〉O
また、微小ピークの場合、ピーク強度尋バックグランド
強度となるので、
J(PK) 岬 ン G−十 BG+))
/2であり、上式は、
(PK) −[3X −十 BG+
)’7/2+ + (BG−)+ (BG+))/2]
>0となり、実験等、経験的にこの式が最も判定精度
が高い。
上述したように、一応ピークと認定された点の検出信号
PKが、その両側の適当に離れたバックグランド部分と
見られる位置の検出信号から内挿されたピーク位置での
バックグランドに、そのバックグランドの平方根のに倍
を加えた値より突出しているか否かを判定することで、
当初認定したピークが眞のピークか否がを再度チエツク
しているのである。
【実施例】
第1図〜第4図に本発明のデータ処理方法の一実施例を
示す、第1図は電子線マイクロアナライザ(EPMA)
等で得られた未処理のスペクトルデータであるにの未処
理のスペクトルデータをスム−ジグ処理し、第2図に示
すようなスペクトルデータに変換する。このスペクトル
データを微分すると、第3図に示す微分曲線(増減関数
曲線)となる、ピーク頂上は増加から減少に変化する地
点であるから、微分曲線において十から−に変わる地点
く微分変曲点)にあると考えられる。
従って、微分曲線において、微分変曲点を検索し、その
微分変曲点における、検出信号強度値(PK)を調べて
記憶させると共に、第4図に示すように、ピーク点にお
ける検出信号強度(PK)に対応する半値幅を推定し、
ピーク頂点から尚早値幅の3倍離れた地点の検出信号強
度値を調べ、同強度値をバックグランド値[(BG−)
、(BG十)]とする、上記ピーク信号強度値(PK)
を(PK)−[3X (BG −+ BG+>7/
2+ + (BG =) + CBG +>
1 /2 コ 〉0に代入して、不等式が成立した
場合、(PK)をピーク値として認定する。
バックグランドの測定位置について説明する。
通常のスペクトルピーク波形も、ガウス分布をなす、し
かし、強度の点においては、直値から1゜5σ離れた地
点でも、ピーク強度の32%位の強度を有しており、3
σ離れた地点では、ピーク強度の1%位の強度となるの
で、EPMAによるX線強度測定では、装置安定性、試
料の平滑性、清浄性等の問題を考慮して、通常、測定誤
差は1%以下を目標としている。このため、バックグン
ド位置は、間に他のピークが存在しない時、3σ離れた
地点を設定することで、ピークによる影響を殆ど除くこ
とができる。なお、半値幅は約1.2σである。[It becomes ``Tama-.'' From this, it becomes about 90% (87%) certainty if it is within the range of 1.5σ from the measured value.
When the measured value within the 1.5σ range no longer intersects with the background within the 1.5σ range, it is determined that the measured value is significantly different from the background, so the determination can be made using the following determination formula. . [<PK) -1,5f- and Go ``x-T] -E (
(BG-) + (BG 10) / 2 + 1.5 , J (PK) Misaki G-10 BG+))
/2, and the above formula is (PK) −[3X −10 BG+
)'7/2+ + (BG-)+ (BG+))/2]
>0, and based on experience such as experiments, this formula has the highest determination accuracy. As mentioned above, the detection signal PK at a point that has been recognized as a peak is added to the background at the peak position, which is interpolated from the detection signals at positions on both sides of the peak that are considered to be background parts. By determining whether it is more prominent than the square root of the ground plus times,
We are checking again to see if the initially recognized peak is the true peak. [Example] Figures 1 to 4 show an example of the data processing method of the present invention. Figure 1 shows an electron beam microanalyzer (EPMA).
The unprocessed spectral data obtained in the above steps is subjected to smoothing processing and converted into spectral data as shown in FIG. Differentiating this spectral data results in a differential curve (increase/decrease function curve) as shown in Figure 3. Since the top of the peak is the point where the change changes from increase to decrease, the differential inflection point is the point where the differential curve changes from 10 to -. ). Therefore, in the differential curve, the differential inflection point is searched, the detected signal strength value (PK) at the differential inflection point is checked and stored, and the detected signal strength (PK) at the peak point is determined as shown in FIG. Estimate the half-width corresponding to PK),
Check the detected signal strength value at a point three times the premature value width from the peak apex, and use the same strength value as the background value [(BG-)
, (BG0)], the above peak signal strength value (PK)
(PK) − [3X (BG −+ BG+>7/
2+ + (BG =) + CBG +>
If the inequality is satisfied by substituting 1/2 ko>0, (PK) is recognized as the peak value. The background measurement position will be explained. A normal spectral peak waveform also has a Gaussian distribution, but in terms of intensity, even at a point 1°5σ away from the direct value, it has an intensity of about 32% of the peak intensity, 3
At a point σ away, the intensity is about 1% of the peak intensity, so when measuring X-ray intensity using EPMA, measurement errors are usually reduced by taking into account issues such as device stability, sample smoothness, and cleanliness. The target is 1% or less. Therefore, by setting the background position at a point 3σ apart when there is no other peak in between, the influence of the peak can be almost eliminated. Note that the half width is approximately 1.2σ.
本発明によれば、検出ピークの認定がより高精度になっ
た。According to the present invention, detection peak recognition has become more accurate.
第1図は本発明の一実施例の未処理のスペクトルデータ
、第2図は上記実施例のスムージング処理したスペクト
ルデータ、第3図は上記実施例の微分曲線(増減関数曲
線)、第4図は上記実施例のデータ処理説明図である。Figure 1 shows the unprocessed spectrum data of one embodiment of the present invention, Figure 2 shows the smoothed spectrum data of the above embodiment, Figure 3 shows the differential curve (increase/decrease function curve) of the above embodiment, and Figure 4 FIG. 2 is an explanatory diagram of data processing in the above embodiment.
Claims (1)
ータを微分して、微分値が+から−に変わる地点(微分
変曲点)を検索し、その微分変曲点をピーク頂点として
取り合えず認定し、同ピーク頂点における検出信号強度
値を(PK)とし、同(PK)から半値幅を推定し、上
記ピーク頂点から両側に上記半値幅のn倍離れた地点の
検出信号強度値を、バックグランド値[(BG−)、(
BG+)]とし、上記ピーク信号強度値(PK)が、 ▲数式、化学式、表等があります▼ を満足させる時、(PK)をピーク値として認定するよ
うにしたことを特徴とする高精度ピーク認定法。[Claims] Smoothed spectrum data obtained by X-ray spectroscopy, etc. is differentiated, a point where the differential value changes from + to - (differential inflection point) is searched, and the differential inflection point is peaked. The detection signal intensity value at the peak apex is set as (PK), the half-width is estimated from the peak apex, and the detected signal at a point n times the half-width on both sides from the peak apex is determined. The intensity value is divided into the background value [(BG-), (
BG+)], and when the above peak signal strength value (PK) satisfies ▲There are mathematical formulas, chemical formulas, tables, etc.▼, (PK) is recognized as the peak value. Certification law.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2200279A JPH07104956B2 (en) | 1990-07-26 | 1990-07-26 | High-accuracy peak certification method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2200279A JPH07104956B2 (en) | 1990-07-26 | 1990-07-26 | High-accuracy peak certification method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH0484280A true JPH0484280A (en) | 1992-03-17 |
| JPH07104956B2 JPH07104956B2 (en) | 1995-11-13 |
Family
ID=16421675
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2200279A Expired - Fee Related JPH07104956B2 (en) | 1990-07-26 | 1990-07-26 | High-accuracy peak certification method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH07104956B2 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010517015A (en) * | 2007-01-17 | 2010-05-20 | イノベイティブ アメリカン テクノロジー, インコーポレイテッド | Advanced pattern recognition system for spectral analysis. |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62103553A (en) * | 1985-10-31 | 1987-05-14 | Shimadzu Corp | Peak detection method in X-ray spectrometer |
| JPH01199276A (en) * | 1988-02-03 | 1989-08-10 | Nec Corp | Check device |
-
1990
- 1990-07-26 JP JP2200279A patent/JPH07104956B2/en not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS62103553A (en) * | 1985-10-31 | 1987-05-14 | Shimadzu Corp | Peak detection method in X-ray spectrometer |
| JPH01199276A (en) * | 1988-02-03 | 1989-08-10 | Nec Corp | Check device |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010517015A (en) * | 2007-01-17 | 2010-05-20 | イノベイティブ アメリカン テクノロジー, インコーポレイテッド | Advanced pattern recognition system for spectral analysis. |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH07104956B2 (en) | 1995-11-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7943915B2 (en) | Method of calibrating a wavelength-modulation spectroscopy apparatus | |
| KR19990071667A (en) | Surface Characterization Device and Method | |
| JP2841258B2 (en) | X-ray fluorescence qualitative analysis method | |
| Moctezuma et al. | Performance analysis of acoustic emission hit detection methods using time features | |
| US4920273A (en) | Z-axis measurement system | |
| US6222199B1 (en) | Ultrathin layer measurement having a controlled ambient of light path | |
| Hilden et al. | Optical quality assurance of GEM foils | |
| JPH0484280A (en) | High precision peak recognition method | |
| Jaffe et al. | A Westerbork survey of rich clusters of galaxies. II-The luminosity function of bright cluster galaxies at 1415 MHz | |
| JP3174186B2 (en) | Gas identification device | |
| CN117419847A (en) | Electromagnetic ultrasonic pre-tightening force testing system for bolt tightening tools | |
| US10408756B2 (en) | Analysis method and analysis apparatus | |
| JPH11304732A (en) | Identification method of analysis element by surface analysis equipment | |
| JPH07117457B2 (en) | Spectral measurement data processing method | |
| JPS63222241A (en) | Data processing method | |
| JPH09184761A (en) | Method for calculating coincidence of two spectra | |
| JP2000121333A (en) | Appearance inspection device and appearance inspection method | |
| JP2743975B2 (en) | Ground liquefaction determination method | |
| CN108827101A (en) | Part circumference runout error automatic detection device and method | |
| Snowden et al. | Radial Velocities of Late-Type Field Subgiant Stars | |
| Cooper et al. | Intensity correction factors for neutron powder diffraction measurements | |
| JPS59135349A (en) | System for measuring transmittance of infrared rays | |
| JP6167192B2 (en) | Lubricant spreading analysis device, method, program, recording medium | |
| JPH06213766A (en) | Method for measuring scattering light and distribution of scattering angle | |
| CN116295840A (en) | Method and system for detecting distortion of interferogram |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071113 Year of fee payment: 12 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081113 Year of fee payment: 13 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091113 Year of fee payment: 14 |
|
| LAPS | Cancellation because of no payment of annual fees |