JPH07134703A - 測定カーブのベースライン決定方法 - Google Patents

測定カーブのベースライン決定方法

Info

Publication number
JPH07134703A
JPH07134703A JP30579193A JP30579193A JPH07134703A JP H07134703 A JPH07134703 A JP H07134703A JP 30579193 A JP30579193 A JP 30579193A JP 30579193 A JP30579193 A JP 30579193A JP H07134703 A JPH07134703 A JP H07134703A
Authority
JP
Japan
Prior art keywords
baseline
points
curve
measurer
measurement curve
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
Application number
JP30579193A
Other languages
English (en)
Inventor
Kazuya Shinya
和也 新屋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP30579193A priority Critical patent/JPH07134703A/ja
Publication of JPH07134703A publication Critical patent/JPH07134703A/ja
Pending legal-status Critical Current

Links

Landscapes

  • Complex Calculations (AREA)

Abstract

(57)【要約】 【目的】 測定者の負担を少なくし、しかも再現性の良
い正しいピーク高さの決定を可能とする。 【構成】 測定カーブの仮ベースライン上の点を複数個
指定し、隣接する2個の指定点から成る組の各組に対し
てg(x)={((yi−yi-1)/(xi−xi-1))・(x−xi)+
yi}・f(x);f(x)={1/(1+exp(−a・(x−xi-
1))}・{1/(1+exp(a・(x−xi))}で定義される単位カ
ーブを生成し、全ての組の単位カーブを加算したものを
実ベースラインとする。

Description

【発明の詳細な説明】
【0001】
【産業上の利用分野】本発明は、デンシトメータ、液体
クロマトグラフィ、ガスクロマトグラフィ等の検出器の
出力カーブから、ベースラインの変動の影響を除いたピ
ーク高さを決定する方法に関する。
【0002】
【従来の技術】クロマトグラフィの検出器の調整を予め
正しく行なっておけば、本来は出力カーブのピークでな
い部分(ベースライン)はゼロライン上に出力されるは
ずであるが、キャリヤの分離等の種々の要因により、ベ
ースラインがゼロラインを外れて変動する。各成分に対
応するピークの高さはベースラインから測定しなければ
ならないため、クロマトグラフの解析を行なう際は、ベ
ースラインを正しく決定することが必要になる。
【0003】従来、このようなベースラインの決定には
次の2つの方法が用いられていた。一つは、CRT等の
ディスプレイ上に測定データのカーブをそのまま表示
し、測定者に、マウスによって手動でベースラインを描
いて貰うという方法である。もう一つは、測定データの
カーブの接線の傾きが所定の角度以上になる点をピーク
の開始点とし、所定の角度以下になる点をピークの終了
点として、それらピーク開始点・終了点で挟まれるピー
ク以外の部分の測定データのカーブをベースラインと定
義する、という方法である。
【0004】
【発明が解決しようとする課題】上記の第1の方法で
は、測定者がベースラインの全ての点を決定しなければ
ならないため、滑らかなラインを生成するために測定者
はマウス等の入力機器を慎重に操作する必要があり、時
間がかかると共に測定者に大きな負担となっていた。ま
た、再現性が低いという問題もある。第2の方法では、
一旦パラメータ(今の場合、ピーク開始点の接線角、終
了点の接線角)を決めてしまえばコンピュータ等によっ
て自動的に行なうことができ、再現性も良いという特長
はあるが、そのパラメータをどのような値に設定するか
により結果が変わってくるため、測定者にとって今度は
パラメータ値決定が難しい問題となる。また、ピーク開
始点・終了点を確実に判断するためには或る程度大きい
角度を設定する必要があるが、この場合、ピーク高さが
低めに出るという欠点がある。
【0005】本発明はこのような課題を解決するために
成されたものであり、その目的とするところは、測定者
にとって負担が少なく、しかも、再現性の良い、正しい
ピーク高さを測定するためのベースラインの決定方法を
提供することにある。
【0006】
【課題を解決するための手段】上記課題を解決するため
に成された本発明に係る測定カーブのベースライン決定
方法は、 a)測定カーブの仮ベースライン上の点を複数個指定し、 b)隣接する2個の指定点から成る組の各組に対して、両
点の間では両点を結ぶ線分に近く、両点の外側ではほぼ
0となる無限階微分可能な関数を生成し、 c)全ての組の上記関数を加算したものを実ベースライン
とすることを特徴としている。
【0007】
【作用】測定カーブの仮ベースライン上の点とは、測定
カーブより一見してベースライン上に在ると考えられる
点のことであり、これらの点は測定者が指定してもよい
し、コンピュータにより自動的に指定させてもよい。測
定者が指定する場合は、測定カーブ中の明らかにベース
ラインであると見られる部分に含まれる点を指定する
他、或るピーク(ベースラインから立ち上がった部分)
の途中でベースラインが屈曲していると考えられる場合
には、そのピークの中(x軸の区間内)に在ると考えら
れる仮想的なベースライン上の点を指定してもよい。こ
の場合、従来のマウス等によりベースラインを全て描く
場合と比較すると、少数の点の指定は測定者にとって遙
かに負担が少ない。
【0008】コンピュータにより指定する場合は、例え
ば次のような方法が考えられる。最初に、測定カーブか
ら接線の傾きが所定の角度以下の部分のみを抜き出し、
その部分に含まれる点の中で、等間隔の点を指定する。
丁度等間隔の点が存在しない場合は、欠けたままとする
か、最も近い点を取る。この場合、指定された点が丁度
ピークの裾(ベースラインからやや高くなる)に当たる
という確率は低いため、上記従来の方法に比べると、ピ
ーク高さが減少するということが少なく、また、指定点
の位置(高さ)は、パラメータ(今の場合、上記角度)
の多少の変動にはそれ程影響されないため、パラメータ
として比較的ラフな値を決めておくことができ、測定者
の負担が少ない。
【0009】隣接する2点の間では両点を結ぶ直線に近
く、その外側ではほぼ0となる関数の例は、下記実施例
において挙げる。各隣接2点間においてこのような関数
を生成し、それらを加算することにより、全区間におい
て無限階微分可能な、滑らかなベースラインが生成され
る。
【0010】
【実施例】本発明の一実施例を図1〜図4により説明す
る。図1は、ベースラインが変動している測定データを
プロットしたカーブの例である。本実施例では、まず、
この測定カーブをCRT等のディスプレイ上に表示し、
測定者に仮ベースライン(図1破線。但し、ディスプレ
イ上には現われない。)上の点を複数個指定してもら
う。
【0011】測定者がマウス、デジタイザ等の入力装置
で点(例えば図1の×印で示した点)を指定すると、コ
ンピュータはその点の座標(…,(xi-1,yi-1),(xi,y
i),(xi+1,yi+1),…)を入力装置から読み取る。そし
て、これら複数の指定点より、隣接する2点毎の組
({(xi-1,yi-1),(xi,yi)},{(xi,yi),(xi+1,
yi+1)}等)を構成し、各組(2点)について以下の単
位カーブ生成を行なう。
【0012】本実施例では、単位カーブの基本となるカ
ーブとして次の関数f(x)を用いる。 f(x)={1/(1+exp(−a・(x−x1))}・{1/
(1+exp(a・(x−x2))} ここで、x1<x2である。関数f(x)は無限階微分可
能な、滑らかな関数であり、そのカーブは図3(a)に
示すように、x1<<x<<x2の部分でほぼ1、x<<
x1及びx>>x2の部分でほぼ0となる台形形状を有す
るカーブである。x=x1の近傍及びx=x2の近傍では
両値(0と1)の間を遷移するが、この遷移領域の広さ
はパラメータaに依存し、aの値が小さい場合は図3
(a)に示すように遷移領域Tr1が狭くなり、aの値が
大きい場合は図3(b)に示すように遷移領域Tr2が広
くなる。従って、比較的小さいaの値を有する上記関数
f(x)を用いることにより、区間(x1,x2)内での
み1、その外では0となる基本カーブを生成することが
できる。
【0013】この基本カーブのf(x)=1の部分(台
の部分)を用いて、関数g(x)を次のように定義する
ことにより、図3(c)に示すように、平面上の任意の
2点{(x1,y1),(x2,y2)}の間でのみ両点を結
ぶ線分として存在し、その他の部分では0となる単位カ
ーブを生成することができる。 g(x)={((y2−y1)/(x2−x1))・(x−x1)+y
1}・f(x)
【0014】従って、複数の点(…,(xi-1,yi-1),(x
i,yi),(xi+1,yi+1),…)の各隣接2点間の区間(…,
[(xi-1,yi-1),(xi,yi)),[(xi,yi),(xi+1,yi
+1)),…)でそれぞれこのような単位カーブを生成し、
これらを加算することにより、図4に示すように、これ
らの点を通る滑らかな(無限階微分可能な)カーブを生
成することができる。
【0015】図1の測定カーブにおいて測定者により複
数の点が指定されたとき、このように各隣接点間で単位
カーブを生成し、それらを加算することにより、ベース
ラインを表わす滑らかなカーブ(図1点線)を決定する
ことができる。こうして決定されたベースラインの各点
のy座標値を測定データの各点の値(y座標値)から減
算することにより、図2に示すように、各ピークの真の
高さを測定することができるようになる。
【0016】なお、上記実施例では基本カーブとして上
記関数f(x)を用いたが、これは一例にしか過ぎず、
無限階微分可能であり、指定区間内でのみ1、その他の
区間で0となる関数はその他にも各種存在する。例え
ば、次の関数f3(x)も上記関数f(x)と全く同様
に使用することができる。 f3(x)={1/(1+exp(−a・(x−x1)3)}・{1
/(1+exp(a・(x−x2)3)}
【0017】
【発明の効果】本発明では、測定者は直観的に決定する
ことが困難なパラメータを指定する必要がなく、少数の
点を指定するのみでよい。それ以降は自動的に滑らかな
ベースラインが決定されるため、測定者の労力が大幅に
軽減される。また、指定点はベースラインが屈曲するよ
うな特徴点のみでよく、このような特徴点の場合には指
定点の位置の変動が少ないため、再現性も比較的良好で
ある。更に、ベースラインは滑らかなカーブとなり、ピ
ークに入り込むことがないため、正しいピーク高さを得
ることができる。
【図面の簡単な説明】
【図1】 本発明の一実施例で処理する測定カーブの一
例を示すグラフ。
【図2】 実施例の測定カーブからベースラインを引い
たカーブのグラフ。
【図3】 実施例で用いた基本カーブのグラフ(a)、
(b)及び単位カーブのグラフ(c)。
【図4】 複数の単位カーブを加算して生成された滑ら
かなカーブのグラフ。

Claims (1)

    【特許請求の範囲】
  1. 【請求項1】 a)測定カーブの仮ベースライン上の点を
    複数個指定し、 b)隣接する2個の指定点から成る組の各組に対して、両
    点の間では両点を結ぶ線分に近く、両点の外側ではほぼ
    0となる無限階微分可能な関数を生成し、 c)全ての組の上記関数を加算したものを実ベースライン
    とすることを特徴とする測定カーブのベースライン決定
    方法。
JP30579193A 1993-11-10 1993-11-10 測定カーブのベースライン決定方法 Pending JPH07134703A (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP30579193A JPH07134703A (ja) 1993-11-10 1993-11-10 測定カーブのベースライン決定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP30579193A JPH07134703A (ja) 1993-11-10 1993-11-10 測定カーブのベースライン決定方法

Publications (1)

Publication Number Publication Date
JPH07134703A true JPH07134703A (ja) 1995-05-23

Family

ID=17949405

Family Applications (1)

Application Number Title Priority Date Filing Date
JP30579193A Pending JPH07134703A (ja) 1993-11-10 1993-11-10 測定カーブのベースライン決定方法

Country Status (1)

Country Link
JP (1) JPH07134703A (ja)

Similar Documents

Publication Publication Date Title
JP3857328B2 (ja) ディスプレイ編集システム
JPH11149562A (ja) 曲線セグメントのコーナーの特性付け
Beardah et al. MATLAB routines for kernel density estimation and the graphical representation of archaeological data
JP2003529114A (ja) 光学スキャナとそのためのソフトウエア
US6675107B2 (en) Data processing apparatus for chromatograph
JP2873883B2 (ja) 線画像の線幅検出方法
JPH07134703A (ja) 測定カーブのベースライン決定方法
JP4570207B2 (ja) 3次元モデル解析装置および記録媒体
JPH0621982B2 (ja) 3次元スペクトル表示装置
JPH07160762A (ja) 見積り支援装置
US5442735A (en) Method and apparatus for interactively manipulating models
US20030083820A1 (en) Extended recursive f-k migration
JP6555414B2 (ja) データ処理方法及び装置
US7702158B2 (en) Image processing apparatus, image processing method and record medium for the image processing apparatus
JPH05272993A (ja) 測定データのベースライン決定方法
JP2005025605A (ja) 物理データフィッティング係数生成システムおよび方法
JPH06103274A (ja) 原価見積りシステム
Kristofferson et al. An automated method for defining microtubule length distributions
JPH11272888A (ja) Cad/cae装置
JPH0793522A (ja) 画像情報処理装置
JPH07239711A (ja) 過渡応答波形特徴量算出方法
JP3980165B2 (ja) 波形データ作成装置
JPH06266848A (ja) 測定した特性曲線の接線を求める方法
JP2817845B2 (ja) 図形処理装置及び方法
JPH01311218A (ja) 分析データ処理装置